astronoby 0.6.0 → 0.8.0

data/lib/astronoby/angle.rb

@@ -22,6 +22,10 @@ module Astronoby
         from_radians(radians)
       end
 
+      def from_degree_arcseconds(arcseconds)
+        from_dms(0, 0, arcseconds)
+      end
+
       def from_hours(hours)
         radians = hours * Constants::RADIAN_PER_HOUR
         from_radians(radians)
@@ -120,10 +124,10 @@ module Astronoby
     end
     alias_method :eql?, :==
 
-    def str(format)
+    def str(format, precision: 4)
       case format
-      when :dms then to_dms(degrees).format
-      when :hms then to_hms(hours).format
+      when :dms then to_dms.format(precision: precision)
+      when :hms then to_hms.format(precision: precision)
       else
         raise UnsupportedFormatError.new(
           "Expected a format between #{FORMATS.join(", ")}, got #{format}"
@@ -131,36 +135,36 @@ module Astronoby
       end
     end
 
-    def to_dms(deg)
-      sign = deg.negative? ? "-" : "+"
-      absolute_degrees = deg.abs
-      degrees = absolute_degrees.floor
+    def to_dms
+      sign = degrees.negative? ? "-" : "+"
+      absolute_degrees = degrees.abs
+      deg = absolute_degrees.floor
       decimal_minutes = Constants::MINUTES_PER_DEGREE *
-        (absolute_degrees - degrees)
+        (absolute_degrees - deg)
       absolute_decimal_minutes = (
-        Constants::MINUTES_PER_DEGREE * (absolute_degrees - degrees)
+        Constants::MINUTES_PER_DEGREE * (absolute_degrees - deg)
       ).abs
       minutes = decimal_minutes.floor
       seconds = Constants::SECONDS_PER_MINUTE * (
         absolute_decimal_minutes - absolute_decimal_minutes.floor
       )
 
-      Dms.new(sign, degrees, minutes, seconds.floor(4))
+      Dms.new(sign, deg, minutes, seconds)
     end
 
-    def to_hms(hrs)
-      absolute_hours = hrs.abs
-      hours = absolute_hours.floor
-      decimal_minutes = Constants::MINUTES_PER_HOUR * (absolute_hours - hours)
+    def to_hms
+      absolute_hours = hours.abs
+      hrs = absolute_hours.floor
+      decimal_minutes = Constants::MINUTES_PER_HOUR * (absolute_hours - hrs)
       absolute_decimal_minutes = (
-        Constants::MINUTES_PER_HOUR * (absolute_hours - hours)
+        Constants::MINUTES_PER_HOUR * (absolute_hours - hrs)
       ).abs
       minutes = decimal_minutes.floor
       seconds = Constants::SECONDS_PER_MINUTE * (
         absolute_decimal_minutes - absolute_decimal_minutes.floor
       )
 
-      Hms.new(hours, minutes, seconds.floor(4))
+      Hms.new(hrs, minutes, seconds)
     end
   end
 end

data/lib/astronoby/angles/dms.rb

@@ -11,8 +11,8 @@ module Astronoby
       @seconds = seconds
     end
 
-    def format
-      "#{sign}#{degrees}° #{minutes}′ #{seconds}″"
+    def format(precision: 4)
+      "#{sign}#{degrees}° #{minutes}′ #{seconds.floor(precision)}″"
     end
   end
 end

data/lib/astronoby/angles/hms.rb

@@ -10,8 +10,8 @@ module Astronoby
       @seconds = seconds
     end
 
-    def format
-      "#{hours}h #{minutes}m #{seconds}s"
+    def format(precision: 4)
+      "#{hours}h #{minutes}m #{seconds.floor(precision)}s"
     end
   end
 end

data/lib/astronoby/bodies/earth.rb

@@ -0,0 +1,62 @@
+# frozen_string_literal: true
+
+module Astronoby
+  class Earth < SolarSystemBody
+    def self.ephemeris_segments(ephem_source)
+      if ephem_source == ::Ephem::SPK::JPL_DE
+        [
+          [SOLAR_SYSTEM_BARYCENTER, EARTH_MOON_BARYCENTER],
+          [EARTH_MOON_BARYCENTER, EARTH]
+        ]
+      elsif ephem_source == ::Ephem::SPK::INPOP
+        [
+          [SOLAR_SYSTEM_BARYCENTER, EARTH]
+        ]
+      end
+    end
+
+    private
+
+    # Attributes that require Sun data like phase angle or magnitude are not
+    # applicable for Earth.
+    def requires_sun_data?
+      false
+    end
+
+    def compute_astrometric(_ephem)
+      Astrometric.new(
+        position: Vector[
+          Distance.zero,
+          Distance.zero,
+          Distance.zero
+        ],
+        velocity: Vector[
+          Velocity.zero,
+          Velocity.zero,
+          Velocity.zero
+        ],
+        instant: @instant,
+        center_identifier: EARTH,
+        target_body: self.class
+      )
+    end
+
+    def compute_mean_of_date(_ephem)
+      MeanOfDate.new(
+        position: Vector[
+          Distance.zero,
+          Distance.zero,
+          Distance.zero
+        ],
+        velocity: Vector[
+          Velocity.zero,
+          Velocity.zero,
+          Velocity.zero
+        ],
+        instant: @instant,
+        center_identifier: EARTH,
+        target_body: self.class
+      )
+    end
+  end
+end

data/lib/astronoby/bodies/jupiter.rb

@@ -0,0 +1,28 @@
+# frozen_string_literal: true
+
+module Astronoby
+  class Jupiter < SolarSystemBody
+    EQUATORIAL_RADIUS = Distance.from_meters(71_492_000)
+    ABSOLUTE_MAGNITUDE = -9.395
+
+    def self.ephemeris_segments(_ephem_source)
+      [[SOLAR_SYSTEM_BARYCENTER, JUPITER_BARYCENTER]]
+    end
+
+    def self.absolute_magnitude
+      ABSOLUTE_MAGNITUDE
+    end
+
+    private
+
+    # Source:
+    #  Title: Computing Apparent Planetary Magnitudes for The Astronomical
+    #    Almanac (2018)
+    #  Authors: Anthony Mallama and James L. Hilton
+    def magnitude_correction_term
+      phase_angle_degrees = phase_angle.degrees
+      -3.7 * 10**-4 * phase_angle_degrees +
+        6.16 * 10**-4 * phase_angle_degrees * phase_angle_degrees
+    end
+  end
+end

data/lib/astronoby/bodies/mars.rb

@@ -0,0 +1,28 @@
+# frozen_string_literal: true
+
+module Astronoby
+  class Mars < SolarSystemBody
+    EQUATORIAL_RADIUS = Distance.from_meters(3_396_200)
+    ABSOLUTE_MAGNITUDE = -1.601
+
+    def self.ephemeris_segments(_ephem_source)
+      [[SOLAR_SYSTEM_BARYCENTER, MARS_BARYCENTER]]
+    end
+
+    def self.absolute_magnitude
+      ABSOLUTE_MAGNITUDE
+    end
+
+    private
+
+    # Source:
+    #  Title: Computing Apparent Planetary Magnitudes for The Astronomical
+    #    Almanac (2018)
+    #  Authors: Anthony Mallama and James L. Hilton
+    def magnitude_correction_term
+      phase_angle_degrees = phase_angle.degrees
+      2.267 * 10**-2 * phase_angle_degrees -
+        1.302 * 10**-4 * phase_angle_degrees * phase_angle_degrees
+    end
+  end
+end

data/lib/astronoby/bodies/mercury.rb

@@ -0,0 +1,32 @@
+# frozen_string_literal: true
+
+module Astronoby
+  class Mercury < SolarSystemBody
+    EQUATORIAL_RADIUS = Distance.from_meters(2_439_700)
+    ABSOLUTE_MAGNITUDE = -0.613
+
+    def self.ephemeris_segments(_ephem_source)
+      [[SOLAR_SYSTEM_BARYCENTER, MERCURY_BARYCENTER]]
+    end
+
+    def self.absolute_magnitude
+      ABSOLUTE_MAGNITUDE
+    end
+
+    private
+
+    # Source:
+    #  Title: Computing Apparent Planetary Magnitudes for The Astronomical
+    #    Almanac (2018)
+    #  Authors: Anthony Mallama and James L. Hilton
+    def magnitude_correction_term
+      phase_angle_degrees = phase_angle.degrees
+      6.328 * 10**-2 * phase_angle_degrees -
+        1.6336 * 10**-3 * phase_angle_degrees * phase_angle_degrees +
+        3.3634 * 10**-5 * phase_angle_degrees**3 -
+        3.4265 * 10**-7 * phase_angle_degrees**4 +
+        1.6893 * 10**-9 * phase_angle_degrees**5 -
+        3.0334 * 10**-12 * phase_angle_degrees**6
+    end
+  end
+end

data/lib/astronoby/bodies/moon.rb

@@ -1,9 +1,24 @@
 # frozen_string_literal: true
 
 module Astronoby
-  class Moon
+  class Moon < SolarSystemBody
     SEMIDIAMETER_VARIATION = 0.7275
-    MEAN_GEOCENTRIC_DISTANCE = Astronoby::Distance.from_meters(385_000_560)
+    EQUATORIAL_RADIUS = Distance.from_meters(1_737_400)
+    ABSOLUTE_MAGNITUDE = 0.28
+
+    def self.ephemeris_segments(ephem_source)
+      if ephem_source == ::Ephem::SPK::JPL_DE
+        [
+          [SOLAR_SYSTEM_BARYCENTER, EARTH_MOON_BARYCENTER],
+          [EARTH_MOON_BARYCENTER, MOON]
+        ]
+      elsif ephem_source == ::Ephem::SPK::INPOP
+        [
+          [SOLAR_SYSTEM_BARYCENTER, EARTH],
+          [EARTH, MOON]
+        ]
+      end
+    end
 
     # Source:
     #  Title: Astronomical Algorithms
@@ -18,323 +33,61 @@ module Astronoby
       Events::MoonPhases.phases_for(year: year, month: month)
     end
 
-    def initialize(time:)
-      @time = time
+    def self.absolute_magnitude
+      ABSOLUTE_MAGNITUDE
     end
 
-    # @return [Astronoby::Coordinates::Ecliptic] Apparent ecliptic coordinates
-    #   of the Moon
-    def apparent_ecliptic_coordinates
-      @ecliptic_coordinates ||= begin
-        latitude = Astronoby::Angle.from_degrees(
-          (latitude_terms + additive_latitude_terms) *
-            EphemerideLunaireParisienne::DEGREES_UNIT
-        )
-
-        longitude = mean_longitude + Astronoby::Angle.from_degrees(
-          (longitude_terms + additive_longitude_terms) *
-            EphemerideLunaireParisienne::DEGREES_UNIT
-        ) + nutation
-
-        Coordinates::Ecliptic.new(
-          latitude: latitude,
-          longitude: longitude
-        )
-      end
+    # @return [Float] Phase fraction, from 0 to 1
+    def current_phase_fraction
+      mean_elongation.degrees / Constants::DEGREES_PER_CIRCLE
     end
 
+    private
+
     # Source:
     #  Title: Astronomical Algorithms
     #  Author: Jean Meeus
     #  Edition: 2nd edition
-    #  Chapter: 47 - Position of the Moon
-
-    # @return [Astronoby::Coordinates::Equatorial] Apparent geocentric
-    #   equatorial coordinates of the Moon
-    def apparent_equatorial_coordinates
-      @apparent_equatorial_coordinates ||=
-        apparent_ecliptic_coordinates
-          .to_apparent_equatorial(epoch: Epoch.from_time(@time))
-    end
-
-    def horizontal_coordinates(observer:)
-      apparent_topocentric_equatorial_coordinates =
-        Astronoby::GeocentricParallax.for_equatorial_coordinates(
-          observer: observer,
-          time: @time,
-          coordinates: apparent_equatorial_coordinates,
-          distance: distance
-        )
-
-      apparent_topocentric_equatorial_coordinates.to_horizontal(
-        observer: observer,
-        time: @time
-      )
-    end
-
-    # @return [Astronoby::Distance] Distance between the Earth and the Moon centers
-    def distance
-      @distance ||= Astronoby::Distance.from_meters(
-        (MEAN_GEOCENTRIC_DISTANCE.meters + distance_terms).round
-      )
-    end
-
-    # @return [Angle] Moon's mean longitude
-    def mean_longitude
-      @mean_longitude ||= Angle.from_degrees(
-        (
-          218.3164477 +
-          481267.88123421 * elapsed_centuries -
-          0.0015786 * elapsed_centuries**2 +
-          elapsed_centuries**3 / 538841 -
-          elapsed_centuries**4 / 65194000
-        ) % 360
-      )
-    end
-
-    # @return [Angle] Moon's mean elongation
+    #  Chapter: 48 - Illuminated Fraction of the Moon's Disk
     def mean_elongation
       @mean_elongation ||= Angle.from_degrees(
         (
           297.8501921 +
-          445267.1114034 * elapsed_centuries -
-          0.0018819 * elapsed_centuries**2 +
-          elapsed_centuries**3 / 545868 -
-          elapsed_centuries**4 / 113065000
-        ) % 360
-      )
-    end
-
-    # @return [Angle] Moon's mean anomaly
-    def mean_anomaly
-      @mean_anomaly ||= Angle.from_degrees(
-        (
-          134.9633964 +
-          477198.8675055 * elapsed_centuries +
-          0.0087414 * elapsed_centuries**2 +
-          elapsed_centuries**3 / 69699 -
-          elapsed_centuries**4 / 14712000
+            445267.1114034 * elapsed_centuries -
+            0.0018819 * elapsed_centuries**2 +
+            elapsed_centuries**3 / 545868 -
+            elapsed_centuries**4 / 113065000
         ) % 360
       )
     end
 
-    # @return [Angle] Moon's argument of latitude
-    def argument_of_latitude
-      @argument_of_latitude ||= Angle.from_degrees(
-        (
-          93.2720950 +
-          483202.0175233 * elapsed_centuries -
-          0.0036539 * elapsed_centuries**2 -
-          elapsed_centuries**3 / 3526000
-        ) % 360
-      )
-    end
-
-    # Source:
-    #  Title: Astronomical Algorithms
-    #  Author: Jean Meeus
-    #  Edition: 2nd edition
-    #  Chapter: 48 - Illuminated Fraction of the Moon's Disk
-
-    # @return [Angle] Moon's phase angle
-    def phase_angle
-      @phase_angle ||= begin
-        sun_apparent_equatorial_coordinates = sun
-          .apparent_ecliptic_coordinates
-          .to_apparent_equatorial(epoch: Epoch.from_time(@time))
-        moon_apparent_equatorial_coordinates = apparent_equatorial_coordinates
-        geocentric_elongation = Angle.acos(
-          sun_apparent_equatorial_coordinates.declination.sin *
-            moon_apparent_equatorial_coordinates.declination.sin +
-            sun_apparent_equatorial_coordinates.declination.cos *
-            moon_apparent_equatorial_coordinates.declination.cos *
-            (
-              sun_apparent_equatorial_coordinates.right_ascension -
-              moon_apparent_equatorial_coordinates.right_ascension
-            ).cos
-        )
-
-        term1 = sun.earth_distance.km * geocentric_elongation.sin
-        term2 = distance.km - sun.earth_distance.km * geocentric_elongation.cos
-        angle = Angle.atan(term1 / term2)
-        Astronoby::Util::Trigonometry
-          .adjustement_for_arctangent(term1, term2, angle)
-      end
-    end
-
-    # @return [Float] Moon's illuminated fraction
-    def illuminated_fraction
-      @illuminated_fraction ||= (1 + phase_angle.cos) / 2
-    end
-
-    # @return [Float] Phase fraction, from 0 to 1
-    def current_phase_fraction
-      mean_elongation.degrees / Constants::DEGREES_PER_CIRCLE
-    end
-
-    # @param observer [Astronoby::Observer] Observer of the event
-    # @return [Astronoby::Events::ObservationEvents] Moon's observation events
-    def observation_events(observer:)
-      today = @time.to_date
-      leap_seconds = Util::Time.terrestrial_universal_time_delta(today)
-      yesterday = today.prev_day
-      yesterday_midnight_terrestrial_time =
-        Time.utc(yesterday.year, yesterday.month, yesterday.day) - leap_seconds
-      today_midnight_terrestrial_time =
-        Time.utc(today.year, today.month, today.day) - leap_seconds
-      tomorrow = today.next_day
-      tomorrow_midnight_terrestrial_time =
-        Time.utc(tomorrow.year, tomorrow.month, tomorrow.day) - leap_seconds
-
-      coordinates_of_the_previous_day = self.class
-        .new(time: yesterday_midnight_terrestrial_time)
-        .apparent_equatorial_coordinates
-      coordinates_of_the_day = self.class
-        .new(time: today_midnight_terrestrial_time)
-        .apparent_equatorial_coordinates
-      coordinates_of_the_next_day = self.class
-        .new(time: tomorrow_midnight_terrestrial_time)
-        .apparent_equatorial_coordinates
-      additional_altitude = -Angle.from_degrees(
-        SEMIDIAMETER_VARIATION *
-          GeocentricParallax.angle(distance: distance).degrees
-      )
-
-      Events::ObservationEvents.new(
-        observer: observer,
-        date: today,
-        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,
-        additional_altitude: additional_altitude
-      )
+    def elapsed_centuries
+      (@instant.tt - JulianDate::DEFAULT_EPOCH) / Constants::DAYS_PER_JULIAN_CENTURY
     end
 
     private
 
-    def terrestrial_time
-      @terrestrial_time ||= begin
-        delta = Util::Time.terrestrial_universal_time_delta(@time)
-        @time + delta
+    # Source:
+    #  Title: Computing Apparent Planetary Magnitudes for The Astronomical
+    #    Almanac (2018)
+    #  Authors: Anthony Mallama and James L. Hilton
+    def magnitude_correction_term
+      phase_angle_degrees = phase_angle.degrees
+      if phase_angle_degrees <= 150 && current_phase_fraction <= 0.5
+        2.9994 * 10**-2 * phase_angle_degrees -
+          1.6057 * 10**-4 * phase_angle_degrees**2 +
+          3.1543 * 10**-6 * phase_angle_degrees**3 -
+          2.0667 * 10**-8 * phase_angle_degrees**4 +
+          6.2553 * 10**-11 * phase_angle_degrees**5
+      elsif phase_angle_degrees <= 150 && current_phase_fraction > 0.5
+        3.3234 * 10**-2 * phase_angle_degrees -
+          3.0725 * 10**-4 * phase_angle_degrees**2 +
+          6.1575 * 10**-6 * phase_angle_degrees**3 -
+          4.7723 * 10**-8 * phase_angle_degrees**4 +
+          1.4681 * 10**-10 * phase_angle_degrees**5
+      else
+        super
       end
     end
-
-    def julian_date
-      Epoch.from_time(terrestrial_time)
-    end
-
-    def elapsed_centuries
-      (julian_date - Epoch::DEFAULT_EPOCH) / Constants::DAYS_PER_JULIAN_CENTURY
-    end
-
-    def sun
-      @sun ||= Sun.new(time: @time)
-    end
-
-    def sun_mean_anomaly
-      @sun_mean_anomaly ||= sun.mean_anomaly
-    end
-
-    def a1
-      @a1 ||= Angle.from_degrees(
-        (119.75 + 131.849 * elapsed_centuries) % 360
-      )
-    end
-
-    def a2
-      @a2 ||= Angle.from_degrees(
-        (53.09 + 479264.290 * elapsed_centuries) % 360
-      )
-    end
-
-    def a3
-      @a3 ||= Angle.from_degrees(
-        (313.45 + 481266.484 * elapsed_centuries) % 360
-      )
-    end
-
-    def eccentricity_correction
-      @eccentricity_correction ||=
-        1 - 0.002516 * elapsed_centuries - 0.0000074 * elapsed_centuries**2
-    end
-
-    def latitude_terms
-      @latitude_terms ||=
-        Astronoby::EphemerideLunaireParisienne
-          .periodic_terms_for_moon_latitude
-          .inject(0) do |sum, terms|
-            value = terms[4] * Math.sin(
-              terms[0] * mean_elongation.radians +
-              terms[1] * sun_mean_anomaly.radians +
-              terms[2] * mean_anomaly.radians +
-              terms[3] * argument_of_latitude.radians
-            )
-
-            value *= eccentricity_correction if terms[1].abs == 1
-            value *= eccentricity_correction**2 if terms[1].abs == 2
-
-            sum + value
-          end
-    end
-
-    def additive_latitude_terms
-      @additive_latitude_terms ||=
-        -2235 * mean_longitude.sin +
-        382 * a3.sin +
-        175 * (a1 - argument_of_latitude).sin +
-        175 * (a1 + argument_of_latitude).sin +
-        127 * (mean_longitude - mean_anomaly).sin -
-        115 * (mean_longitude + mean_anomaly).sin
-    end
-
-    def longitude_terms
-      @longitude_terms ||=
-        Astronoby::EphemerideLunaireParisienne
-          .periodic_terms_for_moon_longitude_and_distance
-          .inject(0) do |sum, terms|
-            value = terms[4] * Math.sin(
-              terms[0] * mean_elongation.radians +
-              terms[1] * sun_mean_anomaly.radians +
-              terms[2] * mean_anomaly.radians +
-              terms[3] * argument_of_latitude.radians
-            )
-
-            value *= eccentricity_correction if terms[1].abs == 1
-            value *= eccentricity_correction**2 if terms[1].abs == 2
-
-            sum + value
-          end
-    end
-
-    def additive_longitude_terms
-      @additive_longitude_terms ||=
-        3958 * a1.sin +
-        1962 * (mean_longitude - argument_of_latitude).sin +
-        318 * a2.sin
-    end
-
-    def distance_terms
-      @distance_terms ||=
-        EphemerideLunaireParisienne
-          .periodic_terms_for_moon_longitude_and_distance
-          .inject(0) do |sum, terms|
-            value = terms[5] * Math.cos(
-              terms[0] * mean_elongation.radians +
-              terms[1] * sun_mean_anomaly.radians +
-              terms[2] * mean_anomaly.radians +
-              terms[3] * argument_of_latitude.radians
-            )
-
-            value *= eccentricity_correction if terms[1].abs == 1
-            value *= eccentricity_correction**2 if terms[1].abs == 2
-
-            sum + value
-          end
-    end
-
-    def nutation
-      @nutation ||=
-        Astronoby::Nutation.for_ecliptic_longitude(epoch: julian_date)
-    end
   end
 end

data/lib/astronoby/bodies/neptune.rb

@@ -0,0 +1,32 @@
+# frozen_string_literal: true
+
+module Astronoby
+  class Neptune < SolarSystemBody
+    EQUATORIAL_RADIUS = Distance.from_meters(24_764_000)
+    ABSOLUTE_MAGNITUDE = -7.0
+
+    def self.ephemeris_segments(_ephem_source)
+      [[SOLAR_SYSTEM_BARYCENTER, NEPTUNE_BARYCENTER]]
+    end
+
+    def self.absolute_magnitude
+      ABSOLUTE_MAGNITUDE
+    end
+
+    private
+
+    # Source:
+    #  Title: Computing Apparent Planetary Magnitudes for The Astronomical
+    #    Almanac (2018)
+    #  Authors: Anthony Mallama and James L. Hilton
+    def magnitude_correction_term
+      phase_angle_degrees = phase_angle.degrees
+      if phase_angle_degrees < 133 && @instant.tt > JulianDate::J2000
+        7.944 * 10**-3 * phase_angle_degrees +
+          9.617 * 10**-5 * phase_angle_degrees * phase_angle_degrees
+      else
+        super
+      end
+    end
+  end
+end