astronoby 0.4.0 → 0.5.0

data/lib/astronoby/events/rise_transit_set_iteration.rb

@@ -0,0 +1,215 @@
+# frozen_string_literal: true
+
+module Astronoby
+  class RiseTransitSetIteration
+    EARTH_SIDEREAL_ROTATION_RATE = 360.98564736629
+
+    # Source:
+    #  Title: Astronomical Algorithms
+    #  Author: Jean Meeus
+    #  Edition: 2nd edition
+    #  Chapter: 15 - Rising, Transit, and Setting
+
+    # @param observer [Astronoby::Observer] Observer
+    # @param date [Date] Date of the event
+    # @param coordinates_of_the_previous_day [Astronoby::Coordinates::Equatorial]
+    #  Coordinates of the body of the previous day
+    # @param coordinates_of_the_day [Astronoby::Coordinates::Equatorial]
+    #   Coordinates of the body of the day
+    # @param coordinates_of_the_next_day [Astronoby::Coordinates::Equatorial]
+    #   Coordinates of the body of the next day
+    # @param shift [Astronoby::Angle] Altitude shift
+    # @param initial_rising [Float] Initial rising
+    # @param initial_transit [Float] Initial transit
+    # @param initial_setting [Float] Initial setting
+    def initialize(
+      observer:,
+      date:,
+      coordinates_of_the_previous_day:,
+      coordinates_of_the_day:,
+      coordinates_of_the_next_day:,
+      shift:,
+      initial_rising:,
+      initial_transit:,
+      initial_setting:
+    )
+      @observer = observer
+      @date = date
+      @coordinates_of_the_previous_day = coordinates_of_the_previous_day
+      @coordinates_of_the_day = coordinates_of_the_day
+      @coordinates_of_the_next_day = coordinates_of_the_next_day
+      @shift = shift
+      @initial_rising = initial_rising
+      @initial_transit = initial_transit
+      @initial_setting = initial_setting
+    end
+
+    # @return [Array<Float>] Iteration results
+    def iterate
+      [
+        delta_m_rising,
+        delta_m_transit,
+        delta_m_setting
+      ]
+    end
+
+    private
+
+    def delta_m_rising
+      (local_horizontal_altitude_rising - @shift).degrees./(
+        Constants::DEGREES_PER_CIRCLE *
+          declination_rising.cos *
+          @observer.latitude.cos *
+          local_hour_angle_rising.sin
+      )
+    end
+
+    def delta_m_transit
+      -local_hour_angle_transit.degrees / Constants::DEGREES_PER_CIRCLE
+    end
+
+    def delta_m_setting
+      (local_horizontal_altitude_setting - @shift).degrees./(
+        Constants::DEGREES_PER_CIRCLE *
+          declination_setting.cos *
+          @observer.latitude.cos *
+          local_hour_angle_setting.sin
+      )
+    end
+
+    def observer_longitude
+      # Longitude must be treated positively westwards from the meridian of
+      # Greenwich, and negatively to the east
+      @observer_longitude ||= -@observer.longitude
+    end
+
+    def apparent_gst_at_midnight
+      @apparent_gst_at_midnight ||= Angle.from_hours(
+        GreenwichSiderealTime.from_utc(
+          Time.utc(@date.year, @date.month, @date.day)
+        ).time
+      )
+    end
+
+    def gst_rising
+      @gst_rising ||= Angle.from_degrees(
+        apparent_gst_at_midnight.degrees +
+          EARTH_SIDEREAL_ROTATION_RATE * @initial_rising
+      )
+    end
+
+    def gst_transit
+      @gst_transit ||= Angle.from_degrees(
+        apparent_gst_at_midnight.degrees +
+          EARTH_SIDEREAL_ROTATION_RATE * @initial_transit
+      )
+    end
+
+    def gst_setting
+      @gst_setting ||= Angle.from_degrees(
+        apparent_gst_at_midnight.degrees +
+          EARTH_SIDEREAL_ROTATION_RATE * @initial_setting
+      )
+    end
+
+    def leap_day_portion
+      @leap_day_portion ||= begin
+        leap_seconds = Util::Time.terrestrial_universal_time_delta(@date)
+        leap_seconds / Constants::SECONDS_PER_DAY
+      end
+    end
+
+    def local_hour_angle_rising
+      @local_hour_angle_rising ||=
+        gst_rising - observer_longitude - right_ascension_rising
+    end
+
+    def local_hour_angle_transit
+      gst_transit - observer_longitude - right_ascension_transit
+    end
+
+    def local_hour_angle_setting
+      @local_hour_angle_setting ||=
+        gst_setting - observer_longitude - right_ascension_setting
+    end
+
+    def local_horizontal_altitude_rising
+      Angle.asin(
+        @observer.latitude.sin * declination_rising.sin +
+          @observer.latitude.cos * declination_rising.cos * local_hour_angle_rising.cos
+      )
+    end
+
+    def local_horizontal_altitude_setting
+      Angle.asin(
+        @observer.latitude.sin * declination_setting.sin +
+          @observer.latitude.cos * declination_setting.cos * local_hour_angle_setting.cos
+      )
+    end
+
+    def right_ascension_rising
+      Angle.from_degrees(
+        Util::Maths.interpolate(
+          [
+            @coordinates_of_the_previous_day.right_ascension.degrees,
+            @coordinates_of_the_day.right_ascension.degrees,
+            @coordinates_of_the_next_day.right_ascension.degrees
+          ],
+          @initial_rising + leap_day_portion
+        )
+      )
+    end
+
+    def right_ascension_transit
+      Angle.from_degrees(
+        Util::Maths.interpolate(
+          [
+            @coordinates_of_the_previous_day.right_ascension.degrees,
+            @coordinates_of_the_day.right_ascension.degrees,
+            @coordinates_of_the_next_day.right_ascension.degrees
+          ],
+          @initial_transit + leap_day_portion
+        )
+      )
+    end
+
+    def right_ascension_setting
+      Angle.from_degrees(
+        Util::Maths.interpolate(
+          [
+            @coordinates_of_the_previous_day.right_ascension.degrees,
+            @coordinates_of_the_day.right_ascension.degrees,
+            @coordinates_of_the_next_day.right_ascension.degrees
+          ],
+          @initial_setting + leap_day_portion
+        )
+      )
+    end
+
+    def declination_rising
+      Angle.from_degrees(
+        Util::Maths.interpolate(
+          [
+            @coordinates_of_the_previous_day.declination.degrees,
+            @coordinates_of_the_day.declination.degrees,
+            @coordinates_of_the_next_day.declination.degrees
+          ],
+          @initial_rising + leap_day_portion
+        )
+      )
+    end
+
+    def declination_setting
+      Angle.from_degrees(
+        Util::Maths.interpolate(
+          [
+            @coordinates_of_the_previous_day.declination.degrees,
+            @coordinates_of_the_day.declination.degrees,
+            @coordinates_of_the_next_day.declination.degrees
+          ],
+          @initial_setting + leap_day_portion
+        )
+      )
+    end
+  end
+end

data/lib/astronoby/geocentric_parallax.rb

@@ -8,104 +8,82 @@ module Astronoby
     #  Edition: Cambridge University Press
     #  Chapter: 39 - Calculating correction for parallax
 
-    ASTRONOMICAL_UNIT_IN_METERS = 149_597_870_700
-    EARTH_FLATTENING_CORRECTION = 0.996647
-    EARTH_EQUATORIAL_RADIUS = 6378140.0
-
     # Equatorial horizontal parallax
-    # @param distance [Numeric] Distance of the body from the center of the
-    #   Earth, in meters
+    # @param distance [Astronoby::Distance] Distance of the body from the center
+    #   of the Earth
     # @return [Astronoby::Angle] Equatorial horizontal parallax angle
     def self.angle(distance:)
-      distance_in_earth_radius = distance / EARTH_EQUATORIAL_RADIUS
-      Angle.asin(1 / distance_in_earth_radius)
+      Angle.from_radians(Angle.from_dms(0, 0, 8.794).sin / distance.au)
     end
 
     # Correct equatorial coordinates with the equatorial horizontal parallax
-    # @param latitude [Astronoby::Angle] Observer's latitude
-    # @param longitude [Astronoby::Angle] Observer's longitude
-    # @param elevation [Numeric] Observer's elevation above sea level in meters
+    # @param observer [Astronoby::Observer] Observer
     # @param time [Time] Date-time of the observation
     # @param coordinates [Astronoby::Coordinates::Equatorial]
     #   Equatorial coordinates of the observed body
-    # @param distance [Numeric] Distance of the observed body from the center of
-    #   the Earth, in meters
+    # @param distance [Astronoby::Distance] Distance of the observed body from
+    #   the center of the Earth
     # @return [Astronoby::Coordinates::Equatorial] Apparent equatorial
     #   coordinates with equatorial horizontal parallax
     def self.for_equatorial_coordinates(
-      latitude:,
-      longitude:,
-      elevation:,
+      observer:,
       time:,
       coordinates:,
       distance:
     )
       new(
-        latitude,
-        longitude,
-        elevation,
+        observer,
         time,
         coordinates,
         distance
       ).apply
     end
 
-    # @param latitude [Astronoby::Angle] Observer's latitude
-    # @param longitude [Astronoby::Angle] Observer's longitude
-    # @param elevation [Numeric] Observer's elevation above sea level in meters
+    # @param observer [Astronoby::Observer] Observer
     # @param time [Time] Date-time of the observation
     # @param coordinates [Astronoby::Coordinates::Equatorial] Equatorial
     #   coordinates of the observed body
-    # @param distance [Numeric] Distance of the observed body from the center of
-    #   the Earth, in meters
+    # @param distance [Astronoby::Distance] Distance of the observed body from
+    #   the center of the Earth
     def initialize(
-      latitude,
-      longitude,
-      elevation,
+      observer,
       time,
       coordinates,
       distance
     )
-      @latitude = latitude
-      @longitude = longitude
-      @elevation = elevation
+      @observer = observer
       @time = time
       @coordinates = coordinates
       @distance = distance
     end
 
     def apply
-      term1 = Angle.atan(EARTH_FLATTENING_CORRECTION * @latitude.tan)
-      quantity1 = term1.cos + elevation_ratio * @latitude.cos
-      quantity2 = EARTH_FLATTENING_CORRECTION * term1.sin +
-        elevation_ratio * @latitude.sin
-
-      term2 = quantity1 * hour_angle.sin
-      term3 = distance_in_earth_radius * declination.cos -
-        quantity1 * hour_angle.cos
+      term1 = Angle.atan(Constants::EARTH_FLATTENING_CORRECTION * latitude.tan)
+      quantity1 = term1.cos + elevation_ratio * latitude.cos
+      quantity2 = Constants::EARTH_FLATTENING_CORRECTION * term1.sin +
+        elevation_ratio * latitude.sin
 
-      delta = Angle.atan(term2 / term3)
+      term1 = -quantity1 * equatorial_horizontal_parallax.sin * hour_angle.sin
+      term2 = declination.cos - quantity1 * equatorial_horizontal_parallax.sin * hour_angle.cos
+      delta_right_ascension = Angle.atan(term1 / term2)
 
-      apparent_hour_angle = hour_angle + delta
-      apparent_right_ascension = right_ascension - delta
-      apparent_declination = Angle.atan(
-        (
-          apparent_hour_angle.cos *
-            (distance_in_earth_radius * declination.sin - quantity2)
-        ) / (
-          distance_in_earth_radius * declination.cos * hour_angle.cos - quantity1
-        )
-      )
+      term1 = (declination.sin - quantity2 * equatorial_horizontal_parallax.sin) * delta_right_ascension.cos
+      term2 = declination.cos - quantity1 * equatorial_horizontal_parallax.sin * hour_angle.cos
+      new_declination = Angle.atan(term1 / term2)
 
       Coordinates::Equatorial.new(
-        right_ascension: apparent_right_ascension,
-        declination: apparent_declination,
+        right_ascension: delta_right_ascension + right_ascension,
+        declination: new_declination,
         epoch: @coordinates.epoch
       )
     end
 
     private
 
+    def latitude
+      @observer.latitude
+    end
+
     def right_ascension
       @coordinates.right_ascension
     end
@@ -115,16 +93,18 @@ module Astronoby
     end
 
     def hour_angle
-      @_hour_angle ||=
-        @coordinates.compute_hour_angle(time: @time, longitude: @longitude)
+      @_hour_angle ||= @coordinates.compute_hour_angle(
+        time: @time,
+        longitude: @observer.longitude
+      )
     end
 
     def elevation_ratio
-      @elevation / EARTH_EQUATORIAL_RADIUS
+      @observer.elevation.meters / Constants::EARTH_EQUATORIAL_RADIUS_IN_METERS.to_f
     end
 
-    def distance_in_earth_radius
-      @distance / EARTH_EQUATORIAL_RADIUS
+    def equatorial_horizontal_parallax
+      self.class.angle(distance: @distance)
     end
   end
 end

data/lib/astronoby/moon_phase.rb

@@ -0,0 +1,43 @@
+# frozen_string_literal: true
+
+module Astronoby
+  class MoonPhase
+    NEW_MOON = :new_moon
+    FIRST_QUARTER = :first_quarter
+    FULL_MOON = :full_moon
+    LAST_QUARTER = :last_quarter
+
+    attr_reader :time, :phase
+
+    # @param time [Time] Time of the Moon phase
+    # @return [Astronoby::MoonPhase] New Moon phase
+    def self.new_moon(time)
+      new(time: time, phase: NEW_MOON)
+    end
+
+    # @param time [Time] Time of the Moon phase
+    # @return [Astronoby::MoonPhase] First quarter Moon phase
+    def self.first_quarter(time)
+      new(time: time, phase: FIRST_QUARTER)
+    end
+
+    # @param time [Time] Time of the Moon phase
+    # @return [Astronoby::MoonPhase] Full Moon phase
+    def self.full_moon(time)
+      new(time: time, phase: FULL_MOON)
+    end
+
+    # @param time [Time] Time of the Moon phase
+    # @return [Astronoby::MoonPhase] Last quarter Moon phase
+    def self.last_quarter(time)
+      new(time: time, phase: LAST_QUARTER)
+    end
+
+    # @param time [Time] Time of the Moon phase
+    # @param phase [Symbol] Moon phase
+    def initialize(time:, phase:)
+      @time = time
+      @phase = phase
+    end
+  end
+end

data/lib/astronoby/observer.rb

@@ -2,7 +2,7 @@
 
 module Astronoby
   class Observer
-    DEFAULT_ELEVATION = 0
+    DEFAULT_ELEVATION = Distance.zero
     DEFAULT_TEMPERATURE = 283.15
     PRESSURE_AT_SEA_LEVEL = 1013.25
     PASCAL_PER_MILLIBAR = 0.01
@@ -14,8 +14,8 @@ module Astronoby
 
     # @param latitude [Angle] geographic latitude of the observer
     # @param longitude [Angle] geographic longitude of the observer
-    # @param elevation [Numeric] geographic elevation (or altitude) of the
-    #   observer above sea level in meters
+    # @param elevation [Astronoby::Distance] geographic elevation (or altitude)
+    #   of the observer above sea level
     # @param temperature [Numeric] temperature at the observer's location in
     #   kelvins
     # @param pressure [Numeric] atmospheric pressure at the observer's
@@ -67,7 +67,9 @@ module Astronoby
     # Barometric formula
     # https://en.wikipedia.org/wiki/Barometric_formula
     def pressure_ratio
-      term1 = EARTH_GRAVITATIONAL_ACCELERATION * MOLAR_MASS_OF_AIR * @elevation
+      term1 = EARTH_GRAVITATIONAL_ACCELERATION *
+        MOLAR_MASS_OF_AIR *
+        @elevation.meters
       term2 = UNIVERSAL_GAS_CONSTANT * @temperature
 
       Math.exp(-term1 / term2)

data/lib/astronoby/refraction.rb

@@ -5,17 +5,16 @@ module Astronoby
     LOW_ALTITUDE_BODY_ANGLE = Angle.from_degrees(15)
     ZENITH = Angle.from_degrees(90)
 
-    def self.angle(coordinates:, observer:)
-      new(coordinates, observer).refraction_angle
+    def self.angle(coordinates:)
+      new(coordinates).refraction_angle
     end
 
-    def self.correct_horizontal_coordinates(coordinates:, observer:)
-      new(coordinates, observer).refract
+    def self.correct_horizontal_coordinates(coordinates:)
+      new(coordinates).refract
     end
 
-    def initialize(coordinates, observer)
+    def initialize(coordinates)
       @coordinates = coordinates
-      @observer = observer
     end
 
     # Source:
@@ -27,8 +26,7 @@ module Astronoby
       Coordinates::Horizontal.new(
         azimuth: @coordinates.azimuth,
         altitude: @coordinates.altitude + refraction_angle,
-        latitude: @coordinates.latitude,
-        longitude: @coordinates.longitude
+        observer: @coordinates.observer
       )
     end
 
@@ -43,11 +41,11 @@ module Astronoby
     private
 
     def pressure
-      @_pressure ||= @observer.pressure
+      @_pressure ||= @coordinates.observer.pressure
     end
 
     def temperature
-      @_temperature ||= @observer.temperature
+      @_temperature ||= @coordinates.observer.temperature
     end
 
     def altitude_in_degrees

data/lib/astronoby/util/maths.rb

@@ -16,7 +16,7 @@ module Astronoby
         def interpolate(values, factor)
           unless factor.between?(0, 1)
             raise IncompatibleArgumentsError,
-              "Interpolation factor must be between 0 and 1"
+              "Interpolation factor must be between 0 and 1, got #{factor}"
           end
 
           if values.length == 3

data/lib/astronoby/version.rb

@@ -1,5 +1,5 @@
 # frozen_string_literal: true
 
 module Astronoby
-  VERSION = "0.4.0"
+  VERSION = "0.5.0"
 end

data/lib/astronoby.rb

@@ -4,7 +4,11 @@ require "astronoby/constants"
 require "astronoby/angle"
 require "astronoby/angles/dms"
 require "astronoby/angles/hms"
+require "astronoby/distance"
 require "astronoby/epoch"
+require "astronoby/astronomical_models/ephemeride_lunaire_parisienne"
+require "astronoby/astronomical_models/moon_phases_periodic_terms"
+require "astronoby/bodies/moon"
 require "astronoby/bodies/sun"
 require "astronoby/coordinates/ecliptic"
 require "astronoby/coordinates/equatorial"
@@ -12,10 +16,13 @@ require "astronoby/coordinates/horizontal"
 require "astronoby/aberration"
 require "astronoby/equinox_solstice"
 require "astronoby/errors"
+require "astronoby/events/moon_phases"
 require "astronoby/events/observation_events"
+require "astronoby/events/rise_transit_set_iteration"
 require "astronoby/events/twilight_events"
 require "astronoby/geocentric_parallax"
 require "astronoby/mean_obliquity"
+require "astronoby/moon_phase"
 require "astronoby/nutation"
 require "astronoby/observer"
 require "astronoby/precession"

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: astronoby
 version: !ruby/object:Gem::Version
-  version: 0.4.0
+  version: 0.5.0
 platform: ruby
 authors:
 - Rémy Hannequin
 autorequire:
 bindir: exe
 cert_chain: []
-date: 2024-04-29 00:00:00.000000000 Z
+date: 2024-06-11 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: matrix
@@ -91,18 +91,25 @@ files:
 - lib/astronoby/angle.rb
 - lib/astronoby/angles/dms.rb
 - lib/astronoby/angles/hms.rb
+- lib/astronoby/astronomical_models/ephemeride_lunaire_parisienne.rb
+- lib/astronoby/astronomical_models/moon_phases_periodic_terms.rb
+- lib/astronoby/bodies/moon.rb
 - lib/astronoby/bodies/sun.rb
 - lib/astronoby/constants.rb
 - lib/astronoby/coordinates/ecliptic.rb
 - lib/astronoby/coordinates/equatorial.rb
 - lib/astronoby/coordinates/horizontal.rb
+- lib/astronoby/distance.rb
 - lib/astronoby/epoch.rb
 - lib/astronoby/equinox_solstice.rb
 - lib/astronoby/errors.rb
+- lib/astronoby/events/moon_phases.rb
 - lib/astronoby/events/observation_events.rb
+- lib/astronoby/events/rise_transit_set_iteration.rb
 - lib/astronoby/events/twilight_events.rb
 - lib/astronoby/geocentric_parallax.rb
 - lib/astronoby/mean_obliquity.rb
+- lib/astronoby/moon_phase.rb
 - lib/astronoby/nutation.rb
 - lib/astronoby/observer.rb
 - lib/astronoby/precession.rb