astronoby 0.0.1 → 0.2.0

data/lib/astronoby/refraction.rb

@@ -0,0 +1,73 @@
+# frozen_string_literal: true
+
+module Astronoby
+  class Refraction
+    LOW_ALTITUDE_BODY_ANGLE = Angle.as_degrees(15)
+    ZENITH = Angle.as_degrees(90)
+
+    def self.angle(coordinates:, observer:)
+      new(coordinates, observer).refraction_angle
+    end
+
+    def self.correct_horizontal_coordinates(coordinates:, observer:)
+      new(coordinates, observer).refract
+    end
+
+    def initialize(coordinates, observer)
+      @coordinates = coordinates
+      @observer = observer
+    end
+
+    # Source:
+    #  Title: Practical Astronomy with your Calculator or Spreadsheet
+    #  Authors: Peter Duffett-Smith and Jonathan Zwart
+    #  Edition: Cambridge University Press
+    #  Chapter: 37 - Refraction
+    def refract
+      Coordinates::Horizontal.new(
+        azimuth: @coordinates.azimuth,
+        altitude: @coordinates.altitude + refraction_angle,
+        latitude: @coordinates.latitude,
+        longitude: @coordinates.longitude
+      )
+    end
+
+    def refraction_angle
+      if @coordinates.altitude > LOW_ALTITUDE_BODY_ANGLE
+        high_altitude_angle
+      else
+        low_altitude_angle
+      end
+    end
+
+    private
+
+    def pressure
+      @_pressure ||= @observer.pressure
+    end
+
+    def temperature
+      @_temperature ||= @observer.temperature
+    end
+
+    def altitude_in_degrees
+      @_altitude_in_degrees ||= @coordinates.altitude.degrees
+    end
+
+    def high_altitude_angle
+      zenith_angle = ZENITH - @coordinates.altitude
+      Angle.as_degrees(0.00452 * pressure * zenith_angle.tan / temperature)
+    end
+
+    def low_altitude_angle
+      term1 = pressure * (
+        0.1594 + 0.0196 * altitude_in_degrees + 0.00002 * altitude_in_degrees**2
+      )
+      term2 = temperature * (
+        1 + 0.505 * altitude_in_degrees + 0.0845 * altitude_in_degrees**2
+      )
+
+      Angle.as_degrees(term1 / term2)
+    end
+  end
+end

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/true_obliquity.rb

@@ -0,0 +1,12 @@
+# frozen_string_literal: true
+
+module Astronoby
+  class TrueObliquity
+    def self.for_epoch(epoch)
+      mean_obliquity = MeanObliquity.for_epoch(epoch)
+      nutation = Nutation.for_obliquity_of_the_ecliptic(epoch: epoch)
+
+      mean_obliquity + nutation
+    end
+  end
+end

data/lib/astronoby/util/astrodynamics.rb

@@ -0,0 +1,60 @@
+# frozen_string_literal: true
+
+module Astronoby
+  module Util
+    module Astrodynamics
+      class << self
+        # Source:
+        #  Title: Celestial Calculations
+        #  Author: J. L. Lawrence
+        #  Edition: MIT Press
+        #  Chapter: 4 - Orbits and Coordinate Systems
+        def eccentric_anomaly_newton_raphson(
+          mean_anomaly,
+          orbital_eccentricity,
+          precision,
+          maximum_iteration_count,
+          current_iteration = 0,
+          solution_on_previous_interation = nil
+        )
+          previous_solution = solution_on_previous_interation&.radians
+
+          solution = if previous_solution.nil?
+            if orbital_eccentricity <= 0.75
+              mean_anomaly.radians
+            else
+              Math::PI
+            end
+          else
+            previous_solution - (
+              (
+                previous_solution -
+                orbital_eccentricity * Math.sin(previous_solution) -
+                mean_anomaly.radians
+              ) / (
+                1 - orbital_eccentricity * Math.cos(previous_solution)
+              )
+            )
+          end
+
+          if current_iteration >= maximum_iteration_count ||
+              (
+                previous_solution &&
+                (solution - previous_solution).abs <= precision
+              )
+            return Angle.as_radians(solution)
+          end
+
+          eccentric_anomaly_newton_raphson(
+            mean_anomaly,
+            orbital_eccentricity,
+            precision,
+            maximum_iteration_count,
+            current_iteration + 1,
+            Angle.as_radians(solution)
+          )
+        end
+      end
+    end
+  end
+end

data/lib/astronoby/util/trigonometry.rb

@@ -0,0 +1,26 @@
+# frozen_string_literal: true
+
+require "bigdecimal/math"
+
+module Astronoby
+  module Util
+    module Trigonometry
+      class << self
+        # Source:
+        #  Title: Celestial Calculations
+        #  Author: J. L. Lawrence
+        #  Edition: MIT Press
+        #  Chapter: 4 - Orbits and Coordinate Systems
+        def adjustement_for_arctangent(y, x, angle)
+          return angle if y.positive? && x.positive?
+
+          if y.negative? && x.positive?
+            return Angle.as_degrees(angle.degrees + 360)
+          end
+
+          Angle.as_degrees(angle.degrees + 180)
+        end
+      end
+    end
+  end
+end

data/lib/astronoby/version.rb

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

data/lib/astronoby.rb

@@ -1,11 +1,26 @@
 # frozen_string_literal: true
 
-require "astronoby/version"
 require "astronoby/angle"
-require "astronoby/angles/degree"
-require "astronoby/angles/radian"
-
-module Astronoby
-  class Error < StandardError; end
-  # Your code goes here...
-end
+require "astronoby/angles/dms"
+require "astronoby/angles/hms"
+require "astronoby/epoch"
+require "astronoby/body"
+require "astronoby/bodies/sun"
+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/trigonometry"
+require "astronoby/true_obliquity"
+require "astronoby/version"

metadata

@@ -1,57 +1,57 @@
 --- !ruby/object:Gem::Specification
 name: astronoby
 version: !ruby/object:Gem::Version
-  version: 0.0.1
+  version: 0.2.0
 platform: ruby
 authors:
 - Rémy Hannequin
-autorequire: 
+autorequire:
 bindir: exe
 cert_chain: []
-date: 2022-04-20 00:00:00.000000000 Z
+date: 2024-03-24 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
-  name: rake
+  name: matrix
   requirement: !ruby/object:Gem::Requirement
     requirements:
     - - "~>"
       - !ruby/object:Gem::Version
-        version: '13.0'
+        version: 0.4.2
   type: :runtime
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
     - - "~>"
       - !ruby/object:Gem::Version
-        version: '13.0'
+        version: 0.4.2
 - !ruby/object:Gem::Dependency
-  name: rspec
+  name: rake
   requirement: !ruby/object:Gem::Requirement
     requirements:
     - - "~>"
       - !ruby/object:Gem::Version
-        version: '3.0'
+        version: '13.0'
   type: :runtime
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
     - - "~>"
       - !ruby/object:Gem::Version
-        version: '3.0'
+        version: '13.0'
 - !ruby/object:Gem::Dependency
-  name: prettier
+  name: rspec
   requirement: !ruby/object:Gem::Requirement
     requirements:
     - - "~>"
       - !ruby/object:Gem::Version
-        version: '2.0'
-  type: :development
+        version: '3.0'
+  type: :runtime
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
     - - "~>"
       - !ruby/object:Gem::Version
-        version: '2.0'
+        version: '3.0'
 - !ruby/object:Gem::Dependency
   name: standard
   requirement: !ruby/object:Gem::Requirement
@@ -66,31 +66,51 @@ dependencies:
     - - "~>"
       - !ruby/object:Gem::Version
         version: '1.3'
-description: Ruby version of the calculations from Astronomical Algorithms by Jean
-  Meeus
+description: Ruby version of the calculations from various books like Celestial Calculations
+  by J. L. Lawrence, Practical Astronomy with your Calculator or Spreadsheet by Peter
+  Duffett-Smith and Jonathan Zwart, or Astronomical Algorithms by Jean Meeus
 email:
 - hello@rhannequ.in
 executables: []
 extensions: []
 extra_rdoc_files: []
 files:
-- ".prettierrc"
 - ".rspec"
-- ".standard.yml"
+- ".tool-versions"
 - CHANGELOG.md
 - CODE_OF_CONDUCT.md
+- CONTRIBUTING.md
 - Gemfile
 - Gemfile.lock
 - LICENSE.txt
 - README.md
 - Rakefile
-- astronoby.gemspec
+- UPGRADING.md
 - lib/astronoby.rb
+- lib/astronoby/aberration.rb
 - lib/astronoby/angle.rb
-- lib/astronoby/angles/degree.rb
-- lib/astronoby/angles/radian.rb
+- lib/astronoby/angles/dms.rb
+- lib/astronoby/angles/hms.rb
+- lib/astronoby/bodies/sun.rb
+- lib/astronoby/body.rb
+- lib/astronoby/coordinates/ecliptic.rb
+- 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/trigonometry.rb
 - lib/astronoby/version.rb
-- sig/astronoby.rbs
 homepage: https://github.com/rhannequin/astronoby
 licenses:
 - MIT
@@ -98,7 +118,7 @@ metadata:
   homepage_uri: https://github.com/rhannequin/astronoby
   source_code_uri: https://github.com/rhannequin/astronoby
   changelog_uri: https://github.com/rhannequin/astronoby/blob/main/CHANGELOG.md
-post_install_message: 
+post_install_message:
 rdoc_options: []
 require_paths:
 - lib
@@ -106,15 +126,15 @@ required_ruby_version: !ruby/object:Gem::Requirement
   requirements:
   - - ">="
     - !ruby/object:Gem::Version
-      version: 2.7.6
+      version: 3.2.0
 required_rubygems_version: !ruby/object:Gem::Requirement
   requirements:
   - - ">="
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubygems_version: 3.3.7
-signing_key: 
+rubygems_version: 3.5.3
+signing_key:
 specification_version: 4
 summary: Astronomical calculations
 test_files: []

data/.prettierrc

@@ -1,11 +0,0 @@
-{
-  "printWidth": 80,
-  "requirePragma": false,
-  "rubyArrayLiteral": true,
-  "rubyHashLabel": true,
-  "rubyModifier": false,
-  "rubySingleQuote": false,
-  "rubyToProc": true,
-  "tabWidth": 2,
-  "trailingComma": "all"
-}

data/.standard.yml

@@ -1,3 +0,0 @@
-# For available configuration options, see:
-#   https://github.com/testdouble/standard
-ruby_version: 2.6

data/astronoby.gemspec

@@ -1,40 +0,0 @@
-# frozen_string_literal: true
-
-require_relative "lib/astronoby/version"
-
-Gem::Specification.new do |spec|
-  spec.name = "astronoby"
-  spec.version = Astronoby::VERSION
-  spec.authors = ["Rémy Hannequin"]
-  spec.email = ["hello@rhannequ.in"]
-
-  spec.summary = "Astronomical calculations"
-  spec.description = "Ruby version of the calculations from Astronomical Algorithms by Jean Meeus"
-  spec.homepage = "https://github.com/rhannequin/astronoby"
-  spec.license = "MIT"
-  spec.required_ruby_version = ">= 2.7.6"
-
-  spec.metadata["homepage_uri"] = spec.homepage
-  spec.metadata["source_code_uri"] = spec.homepage
-  spec.metadata["changelog_uri"] = "#{spec.homepage}/blob/main/CHANGELOG.md"
-
-  # Specify which files should be added to the gem when it is released.
-  # The `git ls-files -z` loads the files in the RubyGem that have been added into git.
-  spec.files = Dir.chdir(__dir__) do
-    `git ls-files -z`.split("\x0").reject do |f|
-      (f == __FILE__) || f.match(%r{\A(?:(?:bin|test|spec|features)/|\.(?:git|travis|circleci)|appveyor)})
-    end
-  end
-  spec.bindir = "exe"
-  spec.executables = spec.files.grep(%r{\Aexe/}) { |f| File.basename(f) }
-  spec.require_paths = ["lib"]
-
-  spec.add_dependency "rake", "~> 13.0"
-  spec.add_dependency "rspec", "~> 3.0"
-
-  spec.add_development_dependency "prettier", "~> 2.0"
-  spec.add_development_dependency "standard", "~> 1.3"
-
-  # For more information and examples about making a new gem, check out our
-  # guide at: https://bundler.io/guides/creating_gem.html
-end

data/lib/astronoby/angles/degree.rb

@@ -1,17 +0,0 @@
-# frozen_string_literal: true
-
-module Astronoby
-  class Degree < Angle
-    def initialize(angle)
-      super(angle, unit: DEGREES)
-    end
-
-    def to_degrees
-      self
-    end
-
-    def to_radians
-      self.class.as_radians(@angle / 180 * PI)
-    end
-  end
-end

data/lib/astronoby/angles/radian.rb

@@ -1,17 +0,0 @@
-# frozen_string_literal: true
-
-module Astronoby
-  class Radian < Angle
-    def initialize(angle)
-      super(angle, unit: RADIANS)
-    end
-
-    def to_degrees
-      self.class.as_degrees(@angle * 180 / PI)
-    end
-
-    def to_radians
-      self
-    end
-  end
-end

data/sig/astronoby.rbs

@@ -1,4 +0,0 @@
-module Astronoby
-  VERSION: String
-  # See the writing guide of rbs: https://github.com/ruby/rbs#guides
-end