astronoby 0.1.0 → 0.2.0

data/lib/astronoby/time/greenwich_sidereal_time.rb

@@ -0,0 +1,86 @@
+# frozen_string_literal: true
+
+module Astronoby
+  class GreenwichSiderealTime
+    JULIAN_CENTURIES_EXPONENTS = [
+      BigDecimal("6.697374558"),
+      BigDecimal("2400.051336"),
+      BigDecimal("0.000025862")
+    ].freeze
+
+    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 - Epoch::J2000) / Epoch::DAYS_PER_JULIAN_CENTURY
+      t0 = (
+        (JULIAN_CENTURIES_EXPONENTS[0] +
+          (JULIAN_CENTURIES_EXPONENTS[1] * t) +
+          (JULIAN_CENTURIES_EXPONENTS[2] * t * t)) % 24
+      ).abs
+
+      ut_in_hours = utc.hour +
+        utc.min / 60.0 +
+        (utc.sec + utc.subsec) / 3600.0
+
+      gmst = BigDecimal("1.002737909") * ut_in_hours + t0
+      gmst += 24 if gmst.negative?
+      gmst -= 24 if gmst > 24
+
+      new(date: date, time: gmst)
+    end
+
+    def initialize(date:, time:)
+      @date = date
+      @time = time
+    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 - Epoch::J2000) / Epoch::DAYS_PER_JULIAN_CENTURY
+
+      t0 = (
+        (JULIAN_CENTURIES_EXPONENTS[0] +
+          (JULIAN_CENTURIES_EXPONENTS[1] * t) +
+          (JULIAN_CENTURIES_EXPONENTS[2] * t * t)) % 24
+      ).abs
+
+      a = @time - t0
+      a += 24 if a.negative?
+      a -= 24 if a > 24
+
+      utc = BigDecimal("0.9972695663") * a
+
+      decimal_hour_to_time(date, utc)
+    end
+
+    def to_lst(longitude:)
+      LocalSiderealTime.from_gst(gst: self, longitude: longitude)
+    end
+
+    private
+
+    def decimal_hour_to_time(date, decimal)
+      absolute_hour = decimal.abs
+      hour = absolute_hour.floor
+      decimal_minute = 60 * (absolute_hour - hour)
+      absolute_decimal_minute = (60 * (absolute_hour - hour)).abs
+      minute = decimal_minute.floor
+      second = 60 * (absolute_decimal_minute - absolute_decimal_minute.floor)
+
+      ::Time.utc(date.year, date.month, date.day, hour, minute, second).round
+    end
+  end
+end

data/lib/astronoby/time/local_sidereal_time.rb

@@ -0,0 +1,41 @@
+# frozen_string_literal: true
+
+module Astronoby
+  class LocalSiderealTime
+    attr_reader :date, :time, :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 += 24 if time.negative?
+      time -= 24 if time > 24
+
+      new(date: date, time: time, longitude: longitude)
+    end
+
+    def initialize(date:, time:, longitude:)
+      @date = date
+      @time = time
+      @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 += 24 if time.negative?
+      time -= 24 if time > 24
+
+      GreenwichSiderealTime.new(date: date, time: time)
+    end
+  end
+end

data/lib/astronoby/version.rb

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

data/lib/astronoby.rb

@@ -10,13 +10,17 @@ require "astronoby/coordinates/ecliptic"
 require "astronoby/coordinates/equatorial"
 require "astronoby/coordinates/horizontal"
 require "astronoby/aberration"
+require "astronoby/equinox_solstice"
 require "astronoby/errors"
+require "astronoby/geocentric_parallax"
 require "astronoby/mean_obliquity"
 require "astronoby/nutation"
+require "astronoby/observer"
 require "astronoby/precession"
 require "astronoby/refraction"
+require "astronoby/time/greenwich_sidereal_time"
+require "astronoby/time/local_sidereal_time"
 require "astronoby/util/astrodynamics"
-require "astronoby/util/time"
 require "astronoby/util/trigonometry"
 require "astronoby/true_obliquity"
 require "astronoby/version"

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: astronoby
 version: !ruby/object:Gem::Version
-  version: 0.1.0
+  version: 0.2.0
 platform: ruby
 authors:
 - Rémy Hannequin
 autorequire:
 bindir: exe
 cert_chain: []
-date: 2024-02-27 00:00:00.000000000 Z
+date: 2024-03-24 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: matrix
@@ -85,6 +85,7 @@ files:
 - LICENSE.txt
 - README.md
 - Rakefile
+- UPGRADING.md
 - lib/astronoby.rb
 - lib/astronoby/aberration.rb
 - lib/astronoby/angle.rb
@@ -96,14 +97,18 @@ files:
 - lib/astronoby/coordinates/equatorial.rb
 - lib/astronoby/coordinates/horizontal.rb
 - lib/astronoby/epoch.rb
+- lib/astronoby/equinox_solstice.rb
 - lib/astronoby/errors.rb
+- lib/astronoby/geocentric_parallax.rb
 - lib/astronoby/mean_obliquity.rb
 - lib/astronoby/nutation.rb
+- lib/astronoby/observer.rb
 - lib/astronoby/precession.rb
 - lib/astronoby/refraction.rb
+- lib/astronoby/time/greenwich_sidereal_time.rb
+- lib/astronoby/time/local_sidereal_time.rb
 - lib/astronoby/true_obliquity.rb
 - lib/astronoby/util/astrodynamics.rb
-- lib/astronoby/util/time.rb
 - lib/astronoby/util/trigonometry.rb
 - lib/astronoby/version.rb
 homepage: https://github.com/rhannequin/astronoby

data/lib/astronoby/util/time.rb

@@ -1,93 +0,0 @@
-# frozen_string_literal: true
-
-module Astronoby
-  module Util
-    module Time
-      class << self
-        # 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 ut_to_gmst(universal_time)
-          julian_day = universal_time.to_date.ajd
-          t = (julian_day - Epoch::J2000) / Epoch::DAYS_PER_JULIAN_CENTURY
-          t0 = (
-            (BigDecimal("6.697374558") +
-            (BigDecimal("2400.051336") * t) +
-            (BigDecimal("0.000025862") * t * t)) % 24
-          ).abs
-
-          ut_in_hours = universal_time.hour +
-            universal_time.min / 60.0 +
-            (universal_time.sec + universal_time.subsec) / 3600.0
-
-          gmst = BigDecimal("1.002737909") * ut_in_hours + t0
-
-          # If gmst negative, add 24 hours to the date
-          # If gmst is greater than 24, subtract 24 hours from the date
-          gmst += 24 if gmst.negative?
-          gmst -= 24 if gmst > 24
-
-          gmst
-        end
-
-        def lst_to_ut(date:, longitude:, lst:)
-          gmst = lst_to_gmst(lst: lst, longitude: longitude)
-          julian_day = date.ajd
-          jd0 = ::DateTime.new(date.year, 1, 1, 0, 0, 0).ajd - 1
-          days_into_the_year = julian_day - jd0
-
-          t = (jd0 - Epoch::J1900) / Epoch::DAYS_PER_JULIAN_CENTURY
-          r = BigDecimal("6.6460656") +
-            BigDecimal("2400.051262") * t +
-            BigDecimal("0.00002581") * t * t
-          b = 24 - r + 24 * (date.year - 1900)
-
-          # If t0 negative, add 24 hours to the date
-          # If t0 is greater than 24, subtract 24 hours from the date
-          t0 = BigDecimal("0.0657098") * days_into_the_year - b
-          t0 += 24 if t0.negative?
-
-          a = gmst - t0
-          a += 24 if a.negative?
-
-          ut = BigDecimal("0.99727") * a
-
-          absolute_hour = ut.abs
-          hour = absolute_hour.floor
-          decimal_minute = 60 * (absolute_hour - hour)
-          absolute_decimal_minute = (60 * (absolute_hour - hour)).abs
-          minute = decimal_minute.floor
-          second = (
-            60 * (absolute_decimal_minute - absolute_decimal_minute.floor)
-          ).round
-
-          ::Time.utc(
-            date.year,
-            date.month,
-            date.day,
-            hour,
-            minute,
-            second
-          )
-        end
-
-        def local_sidereal_time(time:, longitude:)
-          ut_to_gmst(time.utc) + longitude.hours
-        end
-
-        def lst_to_gmst(lst:, longitude:)
-          gmst = lst - longitude.hours
-
-          # If gmst negative, add 24 hours to the date
-          # If gmst is greater than 24, subtract 24 hours from the date
-          gmst += 24 if gmst.negative?
-          gmst -= 24 if gmst > 24
-
-          gmst
-        end
-      end
-    end
-  end
-end