astronoby 0.8.0 → 0.9.0

data/lib/astronoby/time/greenwich_apparent_sidereal_time.rb

@@ -0,0 +1,22 @@
+# frozen_string_literal: true
+
+module Astronoby
+  class GreenwichApparentSiderealTime < GreenwichSiderealTime
+    def self.from_utc(utc)
+      gmst = GreenwichMeanSiderealTime.from_utc(utc)
+      instant = Instant.from_time(utc)
+      nutation = Nutation.new(instant: instant)
+      mean_obliquity = MeanObliquity.at(instant)
+
+      equation_of_equinoxes = nutation.nutation_in_longitude.hours *
+        mean_obliquity.cos
+      gast_time = normalize_time(gmst.time + equation_of_equinoxes)
+
+      new(date: gmst.date, time: gast_time)
+    end
+
+    def initialize(date:, time:)
+      super(date: date, time: time, type: APPARENT)
+    end
+  end
+end

data/lib/astronoby/time/greenwich_mean_sidereal_time.rb

@@ -0,0 +1,64 @@
+# frozen_string_literal: true
+
+module Astronoby
+  class GreenwichMeanSiderealTime < GreenwichSiderealTime
+    JULIAN_CENTURIES_EXPONENTS = [
+      6.697374558,
+      2400.051336,
+      0.000025862
+    ].freeze
+
+    SIDEREAL_MINUTE_IN_UT_MINUTE = 0.9972695663
+
+    # Source:
+    #  Title: Practical Astronomy with your Calculator or Spreadsheet
+    #  Authors: Peter Duffett-Smith and Jonathan Zwart
+    #  Edition: Cambridge University Press
+    #  Chapter: 12 - Conversion of UT to Greenwich sidereal time (GST)
+    def self.from_utc(utc)
+      date = utc.to_date
+      julian_day = utc.to_date.ajd
+      t = (julian_day - JulianDate::J2000) / Constants::DAYS_PER_JULIAN_CENTURY
+      t0 = (
+        (JULIAN_CENTURIES_EXPONENTS[0] +
+          (JULIAN_CENTURIES_EXPONENTS[1] * t) +
+          (JULIAN_CENTURIES_EXPONENTS[2] * t * t)) % Constants::HOURS_PER_DAY
+      ).abs
+
+      ut_in_hours = utc.hour +
+        utc.min / Constants::MINUTES_PER_HOUR +
+        (utc.sec + utc.subsec) / Constants::SECONDS_PER_HOUR
+
+      gmst = normalize_time(1.002737909 * ut_in_hours + t0)
+
+      new(date: date, time: gmst)
+    end
+
+    def initialize(date:, time:)
+      super(date: date, time: time, type: MEAN)
+    end
+
+    # Source:
+    #  Title: Practical Astronomy with your Calculator or Spreadsheet
+    #  Authors: Peter Duffett-Smith and Jonathan Zwart
+    #  Edition: Cambridge University Press
+    #  Chapter: 13 - Conversion of GST to UT
+    def to_utc
+      date = @date
+      julian_day = @date.ajd
+      t = (julian_day - JulianDate::J2000) / Constants::DAYS_PER_JULIAN_CENTURY
+
+      t0 = (
+        (JULIAN_CENTURIES_EXPONENTS[0] +
+          (JULIAN_CENTURIES_EXPONENTS[1] * t) +
+          (JULIAN_CENTURIES_EXPONENTS[2] * t * t)) % Constants::HOURS_PER_DAY
+      ).abs
+
+      a = normalize_time(@time - t0)
+
+      utc = SIDEREAL_MINUTE_IN_UT_MINUTE * a
+
+      Util::Time.decimal_hour_to_time(date, 0, utc)
+    end
+  end
+end

data/lib/astronoby/time/greenwich_sidereal_time.rb

@@ -1,71 +1,33 @@
 # frozen_string_literal: true
 
 module Astronoby
-  class GreenwichSiderealTime
-    JULIAN_CENTURIES_EXPONENTS = [
-      6.697374558,
-      2400.051336,
-      0.000025862
-    ].freeze
-
-    SIDEREAL_MINUTE_IN_UT_MINUTE = 0.9972695663
-
-    attr_reader :date, :time
-
-    # Source:
-    #  Title: Practical Astronomy with your Calculator or Spreadsheet
-    #  Authors: Peter Duffett-Smith and Jonathan Zwart
-    #  Edition: Cambridge University Press
-    #  Chapter: 12 - Conversion of UT to Greenwich sidereal time (GST)
-    def self.from_utc(utc)
-      date = utc.to_date
-      julian_day = utc.to_date.ajd
-      t = (julian_day - JulianDate::J2000) / Constants::DAYS_PER_JULIAN_CENTURY
-      t0 = (
-        (JULIAN_CENTURIES_EXPONENTS[0] +
-          (JULIAN_CENTURIES_EXPONENTS[1] * t) +
-          (JULIAN_CENTURIES_EXPONENTS[2] * t * t)) % Constants::HOURS_PER_DAY
-      ).abs
-
-      ut_in_hours = utc.hour +
-        utc.min / Constants::MINUTES_PER_HOUR +
-        (utc.sec + utc.subsec) / Constants::SECONDS_PER_HOUR
-
-      gmst = 1.002737909 * ut_in_hours + t0
-      gmst += Constants::HOURS_PER_DAY if gmst.negative?
-      gmst -= Constants::HOURS_PER_DAY if gmst > Constants::HOURS_PER_DAY
+  class GreenwichSiderealTime < SiderealTime
+    def self.from_utc(utc, type: MEAN)
+      validate_type!(type)
+      case type
+      when MEAN
+        GreenwichMeanSiderealTime.from_utc(utc)
+      when APPARENT
+        GreenwichApparentSiderealTime.from_utc(utc)
+      end
+    end
 
-      new(date: date, time: gmst)
+    def self.mean_from_utc(utc)
+      GreenwichMeanSiderealTime.from_utc(utc)
     end
 
-    def initialize(date:, time:)
-      @date = date
-      @time = time
+    def self.apparent_from_utc(utc)
+      GreenwichApparentSiderealTime.from_utc(utc)
     end
 
-    # Source:
-    #  Title: Practical Astronomy with your Calculator or Spreadsheet
-    #  Authors: Peter Duffett-Smith and Jonathan Zwart
-    #  Edition: Cambridge University Press
-    #  Chapter: 13 - Conversion of GST to UT
     def to_utc
-      date = @date
-      julian_day = @date.ajd
-      t = (julian_day - JulianDate::J2000) / Constants::DAYS_PER_JULIAN_CENTURY
-
-      t0 = (
-        (JULIAN_CENTURIES_EXPONENTS[0] +
-          (JULIAN_CENTURIES_EXPONENTS[1] * t) +
-          (JULIAN_CENTURIES_EXPONENTS[2] * t * t)) % Constants::HOURS_PER_DAY
-      ).abs
-
-      a = @time - t0
-      a += Constants::HOURS_PER_DAY if a.negative?
-      a -= Constants::HOURS_PER_DAY if a > Constants::HOURS_PER_DAY
-
-      utc = SIDEREAL_MINUTE_IN_UT_MINUTE * a
+      unless mean?
+        raise NotImplementedError,
+          "UTC conversion only supported for mean sidereal time"
+      end
 
-      Util::Time.decimal_hour_to_time(date, 0, utc)
+      gmst = GreenwichMeanSiderealTime.new(date: @date, time: @time)
+      gmst.to_utc
     end
 
     def to_lst(longitude:)

data/lib/astronoby/time/local_apparent_sidereal_time.rb

@@ -0,0 +1,42 @@
+# frozen_string_literal: true
+
+module Astronoby
+  class LocalApparentSiderealTime < LocalSiderealTime
+    # Source:
+    #  Title: Practical Astronomy with your Calculator or Spreadsheet
+    #  Authors: Peter Duffett-Smith and Jonathan Zwart
+    #  Edition: Cambridge University Press
+    #  Chapter: 14 - Local sidereal time (LST)
+    def self.from_gst(gst:, longitude:)
+      unless gst.apparent?
+        raise ArgumentError, "GST must be apparent sidereal time"
+      end
+
+      date = gst.date
+      time = normalize_time(gst.time + longitude.hours)
+
+      new(date: date, time: time, longitude: longitude)
+    end
+
+    def self.from_utc(utc, longitude:)
+      gast = GreenwichApparentSiderealTime.from_utc(utc)
+      from_gst(gst: gast, longitude: longitude)
+    end
+
+    def initialize(date:, time:, longitude:)
+      super(date: date, time: time, longitude: longitude, type: APPARENT)
+    end
+
+    # Source:
+    #  Title: Practical Astronomy with your Calculator or Spreadsheet
+    #  Authors: Peter Duffett-Smith and Jonathan Zwart
+    #  Edition: Cambridge University Press
+    #  Chapter: 15 - Converting LST to GST
+    def to_gst
+      date = @date
+      time = normalize_time(@time - @longitude.hours)
+
+      GreenwichApparentSiderealTime.new(date: date, time: time)
+    end
+  end
+end

data/lib/astronoby/time/local_mean_sidereal_time.rb

@@ -0,0 +1,42 @@
+# frozen_string_literal: true
+
+module Astronoby
+  class LocalMeanSiderealTime < LocalSiderealTime
+    # Source:
+    #  Title: Practical Astronomy with your Calculator or Spreadsheet
+    #  Authors: Peter Duffett-Smith and Jonathan Zwart
+    #  Edition: Cambridge University Press
+    #  Chapter: 14 - Local sidereal time (LST)
+    def self.from_gst(gst:, longitude:)
+      unless gst.mean?
+        raise ArgumentError, "GST must be mean sidereal time"
+      end
+
+      date = gst.date
+      time = normalize_time(gst.time + longitude.hours)
+
+      new(date: date, time: time, longitude: longitude)
+    end
+
+    def self.from_utc(utc, longitude:)
+      gmst = GreenwichMeanSiderealTime.from_utc(utc)
+      from_gst(gst: gmst, longitude: longitude)
+    end
+
+    def initialize(date:, time:, longitude:)
+      super(date: date, time: time, longitude: longitude, type: MEAN)
+    end
+
+    # Source:
+    #  Title: Practical Astronomy with your Calculator or Spreadsheet
+    #  Authors: Peter Duffett-Smith and Jonathan Zwart
+    #  Edition: Cambridge University Press
+    #  Chapter: 15 - Converting LST to GST
+    def to_gst
+      date = @date
+      time = normalize_time(@time - @longitude.hours)
+
+      GreenwichMeanSiderealTime.new(date: date, time: time)
+    end
+  end
+end

data/lib/astronoby/time/local_sidereal_time.rb

@@ -1,41 +1,50 @@
 # frozen_string_literal: true
 
 module Astronoby
-  class LocalSiderealTime
-    attr_reader :date, :time, :longitude
+  class LocalSiderealTime < SiderealTime
+    attr_reader :longitude
 
-    # Source:
-    #  Title: Practical Astronomy with your Calculator or Spreadsheet
-    #  Authors: Peter Duffett-Smith and Jonathan Zwart
-    #  Edition: Cambridge University Press
-    #  Chapter: 14 - Local sidereal time (LST)
     def self.from_gst(gst:, longitude:)
-      date = gst.date
-      time = gst.time + longitude.hours
-      time += Constants::HOURS_PER_DAY if time.negative?
-      time -= Constants::HOURS_PER_DAY if time > Constants::HOURS_PER_DAY
+      case gst.type
+      when MEAN
+        LocalMeanSiderealTime.from_gst(gst: gst, longitude: longitude)
+      when APPARENT
+        LocalApparentSiderealTime.from_gst(gst: gst, longitude: longitude)
+      end
+    end
 
-      new(date: date, time: time, longitude: longitude)
+    def self.from_utc(utc, longitude:, type: MEAN)
+      validate_type!(type)
+      case type
+      when MEAN
+        LocalMeanSiderealTime.from_utc(utc, longitude: longitude)
+      when APPARENT
+        LocalApparentSiderealTime.from_utc(utc, longitude: longitude)
+      end
     end
 
-    def initialize(date:, time:, longitude:)
-      @date = date
-      @time = time
+    def initialize(date:, time:, longitude:, type: MEAN)
+      super(date: date, time: time, type: type)
       @longitude = longitude
     end
 
-    # Source:
-    #  Title: Practical Astronomy with your Calculator or Spreadsheet
-    #  Authors: Peter Duffett-Smith and Jonathan Zwart
-    #  Edition: Cambridge University Press
-    #  Chapter: 15 - Converting LST to GST
     def to_gst
-      date = @date
-      time = @time - @longitude.hours
-      time += Constants::HOURS_PER_DAY if time.negative?
-      time -= Constants::HOURS_PER_DAY if time > Constants::HOURS_PER_DAY
-
-      GreenwichSiderealTime.new(date: date, time: time)
+      case @type
+      when MEAN
+        lst = LocalMeanSiderealTime.new(
+          date: @date,
+          time: @time,
+          longitude: @longitude
+        )
+        lst.to_gst
+      when APPARENT
+        last = LocalApparentSiderealTime.new(
+          date: @date,
+          time: @time,
+          longitude: @longitude
+        )
+        last.to_gst
+      end
     end
   end
 end

data/lib/astronoby/time/sidereal_time.rb

@@ -0,0 +1,42 @@
+# frozen_string_literal: true
+
+module Astronoby
+  class SiderealTime
+    TYPES = [
+      MEAN = :mean,
+      APPARENT = :apparent
+    ].freeze
+
+    attr_reader :date, :time, :type
+
+    def self.normalize_time(time)
+      time += Constants::HOURS_PER_DAY if time.negative?
+      time -= Constants::HOURS_PER_DAY if time > Constants::HOURS_PER_DAY
+      time
+    end
+
+    def self.validate_type!(type)
+      unless TYPES.include?(type)
+        raise ArgumentError, "Invalid type: #{type}. Must be one of #{TYPES}"
+      end
+    end
+
+    def initialize(date:, time:, type: MEAN)
+      @date = date
+      @time = time
+      @type = type
+    end
+
+    def mean?
+      @type == MEAN
+    end
+
+    def apparent?
+      @type == APPARENT
+    end
+
+    def normalize_time(time)
+      self.class.normalize_time(time)
+    end
+  end
+end

data/lib/astronoby/util/time.rb

@@ -4,7 +4,7 @@ module Astronoby
   module Util
     module Time
       OLD_LEAP_SECONDS = {
-        2378496.5 => 13.188148343557259,
+        2378496.5 => 13.188148343557259, # 1800-01-01T00:00:00+00:00
         2378677.5 => 12.977779959648615,
         2378861.5 => 12.789307379498496,
         2379042.5 => 12.627219619724201,
@@ -400,52 +400,70 @@ module Astronoby
         2450265.5 => 61.2837935622083,
         2450449.5 => 61.96410528309934,
         2450630.5 => 62.7581330776884,
-        2450814.5 => 63,
-        2450995.5 => 63,
-        2451179.5 => 64
+        2450814.5 => 62.966,
+        2450995.5 => 63.284,
+        2451179.5 => 63.467,
+        2451360.5 => 63.664 # 1999-07-01T00:00:00+00:00
       }.freeze
 
       RECENT_LEAP_SECONDS = {
-        2451360.5 => 64,
-        2451544.5 => 64,
-        2451726.5 => 64,
-        2451910.5 => 64,
-        2452091.5 => 64,
-        2452275.5 => 64,
-        2452456.5 => 64,
-        2452640.5 => 64,
-        2452821.5 => 64,
-        2453005.5 => 64,
-        2453187.5 => 64,
-        2453371.5 => 64,
-        2453552.5 => 64,
-        2453736.5 => 65,
-        2453917.5 => 65,
-        2454101.5 => 65,
-        2454282.5 => 65,
-        2454466.5 => 65,
-        2454648.5 => 65,
-        2454832.5 => 66,
-        2455013.5 => 66,
-        2455197.5 => 66,
-        2455378.5 => 66,
-        2455562.5 => 66,
-        2455743.5 => 66,
-        2455927.5 => 66,
-        2456109.5 => 67,
-        2456293.5 => 67,
-        2456474.5 => 67,
-        2456658.5 => 67,
-        2456839.5 => 67,
-        2457023.5 => 67,
-        2457204.5 => 68,
-        2457388.5 => 68,
-        2457570.5 => 68,
-        2457754.5 => 69
+        2451544.5 => 63.829, # 2000-01-01T00:00:00+00:00
+        2451726.5 => 63.980,
+        2451910.5 => 64.091,
+        2452091.5 => 64.212,
+        2452275.5 => 64.300,
+        2452456.5 => 64.413,
+        2452640.5 => 64.473,
+        2452821.5 => 64.551,
+        2453005.5 => 64.574,
+        2453187.5 => 64.653,
+        2453371.5 => 64.688,
+        2453552.5 => 64.799,
+        2453736.5 => 64.845,
+        2453917.5 => 64.989,
+        2454101.5 => 65.146,
+        2454282.5 => 65.341,
+        2454466.5 => 65.457,
+        2454648.5 => 65.629,
+        2454832.5 => 65.777,
+        2455013.5 => 65.951,
+        2455197.5 => 66.070,
+        2455378.5 => 66.241,
+        2455562.5 => 66.325,
+        2455743.5 => 66.475,
+        2455927.5 => 66.603,
+        2456109.5 => 66.771,
+        2456293.5 => 66.907,
+        2456474.5 => 67.127,
+        2456658.5 => 67.281,
+        2456839.5 => 67.486,
+        2457023.5 => 67.644,
+        2457204.5 => 67.861,
+        2457388.5 => 68.102,
+        2457570.5 => 68.396,
+        2457754.5 => 68.593,
+        2457935.5 => 68.824,
+        2458119.5 => 68.968,
+        2458300.5 => 69.113,
+        2458484.5 => 69.220,
+        2458665.5 => 69.358,
+        2458849.5 => 69.361,
+        2459031.5 => 69.424,
+        2459215.5 => 69.359,
+        2459396.5 => 69.351,
+        2459580.5 => 69.294,
+        2459761.5 => 69.253,
+        2459945.5 => 69.204,
+        2460126.5 => 69.220,
+        2460310.5 => 69.175,
+        2460492.5 => 69.188,
+        2460676.5 => 69.138,
+        2460857.5 => 69.139 # 2025-07-01T00:00:00+00:00
       }.freeze
 
       OLDEST_JD = 2378496.5
-      MOST_RECENT_JD = 2457754.5
+      FIRST_RECENT_JD = 2451544.5
+      MOST_RECENT_JD = 2460857.5
 
       # @param date [Date]
       # @param decimal [Numeric] Hour of the day, in decimal hours
@@ -502,10 +520,10 @@ module Astronoby
             "Expected a Numeric, Time, Date or DateTime object, got #{instant.class}"
         end
 
-        return RECENT_LEAP_SECONDS[MOST_RECENT_JD] if jd >= 2457754.5
+        return RECENT_LEAP_SECONDS[MOST_RECENT_JD] if jd >= MOST_RECENT_JD
         return 0 if jd < OLDEST_JD
 
-        if jd >= 2451360.5
+        if jd >= FIRST_RECENT_JD
           closest_jd = RECENT_LEAP_SECONDS.keys.bsearch { |key| key >= jd }
           leap_seconds = RECENT_LEAP_SECONDS[closest_jd]
         else

data/lib/astronoby/velocity.rb

@@ -40,6 +40,11 @@ module Astronoby
       end
       alias_method :vector_from_mps, :vector_from_meters_per_second
 
+      def vector_from_astronomical_units_per_day(array)
+        Vector.elements(array.map { from_aupd(_1) })
+      end
+      alias_method :vector_from_aupd, :vector_from_astronomical_units_per_day
+
       def light_speed
         from_meters_per_second(Constants::LIGHT_SPEED_M_PER_S)
       end

data/lib/astronoby/version.rb

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

data/lib/astronoby.rb

@@ -6,6 +6,7 @@ require "astronoby/constants"
 require "astronoby/angle"
 require "astronoby/angles/dms"
 require "astronoby/angles/hms"
+require "astronoby/angular_velocity"
 require "astronoby/constellation"
 require "astronoby/constellations/repository"
 require "astronoby/constellations/data"
@@ -33,9 +34,13 @@ require "astronoby/coordinates/equatorial"
 require "astronoby/coordinates/horizontal"
 require "astronoby/corrections/light_time_delay"
 require "astronoby/aberration"
+require "astronoby/bodies/deep_sky_object"
+require "astronoby/bodies/deep_sky_object_position"
 require "astronoby/deflection"
 require "astronoby/equinox_solstice"
 require "astronoby/errors"
+require "astronoby/events/extremum_event"
+require "astronoby/events/extremum_calculator"
 require "astronoby/events/moon_phases"
 require "astronoby/events/rise_transit_set_event"
 require "astronoby/events/rise_transit_set_events"
@@ -56,8 +61,14 @@ require "astronoby/reference_frames/mean_of_date"
 require "astronoby/reference_frames/apparent"
 require "astronoby/reference_frames/topocentric"
 require "astronoby/refraction"
+require "astronoby/stellar_propagation"
+require "astronoby/time/sidereal_time"
 require "astronoby/time/greenwich_sidereal_time"
+require "astronoby/time/greenwich_mean_sidereal_time"
+require "astronoby/time/greenwich_apparent_sidereal_time"
 require "astronoby/time/local_sidereal_time"
+require "astronoby/time/local_mean_sidereal_time"
+require "astronoby/time/local_apparent_sidereal_time"
 require "astronoby/util/maths"
 require "astronoby/util/time"
 require "astronoby/util/trigonometry"

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: astronoby
 version: !ruby/object:Gem::Version
-  version: 0.8.0
+  version: 0.9.0
 platform: ruby
 authors:
 - Rémy Hannequin
@@ -154,9 +154,9 @@ files:
 - UPGRADING.md
 - docs/README.md
 - docs/angles.md
-- docs/celestial_bodies.md
 - docs/configuration.md
 - docs/coordinates.md
+- docs/deep_sky_bodies.md
 - docs/ephem.md
 - docs/equinoxes_solstices_times.md
 - docs/glossary.md
@@ -165,13 +165,17 @@ files:
 - docs/observer.md
 - docs/reference_frames.md
 - docs/rise_transit_set_times.md
+- docs/solar_system_bodies.md
 - docs/twilight_times.md
 - lib/astronoby.rb
 - lib/astronoby/aberration.rb
 - lib/astronoby/angle.rb
 - lib/astronoby/angles/dms.rb
 - lib/astronoby/angles/hms.rb
+- lib/astronoby/angular_velocity.rb
 - lib/astronoby/astronomical_models/moon_phases_periodic_terms.rb
+- lib/astronoby/bodies/deep_sky_object.rb
+- lib/astronoby/bodies/deep_sky_object_position.rb
 - lib/astronoby/bodies/earth.rb
 - lib/astronoby/bodies/jupiter.rb
 - lib/astronoby/bodies/mars.rb
@@ -204,6 +208,8 @@ files:
 - lib/astronoby/ephem.rb
 - lib/astronoby/equinox_solstice.rb
 - lib/astronoby/errors.rb
+- lib/astronoby/events/extremum_calculator.rb
+- lib/astronoby/events/extremum_event.rb
 - lib/astronoby/events/moon_phases.rb
 - lib/astronoby/events/rise_transit_set_calculator.rb
 - lib/astronoby/events/rise_transit_set_event.rb
@@ -226,8 +232,14 @@ files:
 - lib/astronoby/reference_frames/mean_of_date.rb
 - lib/astronoby/reference_frames/topocentric.rb
 - lib/astronoby/refraction.rb
+- lib/astronoby/stellar_propagation.rb
+- lib/astronoby/time/greenwich_apparent_sidereal_time.rb
+- lib/astronoby/time/greenwich_mean_sidereal_time.rb
 - lib/astronoby/time/greenwich_sidereal_time.rb
+- lib/astronoby/time/local_apparent_sidereal_time.rb
+- lib/astronoby/time/local_mean_sidereal_time.rb
 - lib/astronoby/time/local_sidereal_time.rb
+- lib/astronoby/time/sidereal_time.rb
 - lib/astronoby/true_obliquity.rb
 - lib/astronoby/util/maths.rb
 - lib/astronoby/util/time.rb