astronoby 0.6.0 → 0.8.0

data/lib/astronoby/data/constellations/sorted_declinations.dat

@@ -0,0 +1,202 @@
+-85.0
+-82.5
+-76.0
+-75.0
+-70.0
+-67.5
+-65.0
+-64.0
+-63.583333333333336
+-61.0
+-60.0
+-59.0
+-58.5
+-58.0
+-57.5
+-57.0
+-56.5
+-55.0
+-54.5
+-54.0
+-53.5
+-53.166666666666664
+-53.0
+-52.5
+-51.5
+-51.0
+-50.75
+-50.0
+-49.0
+-48.166666666666664
+-48.0
+-46.5
+-46.0
+-45.5
+-44.0
+-43.0
+-42.0
+-40.0
+-39.75
+-39.583333333333336
+-37.0
+-36.75
+-36.0
+-35.0
+-33.0
+-31.166666666666668
+-30.0
+-29.5
+-29.166666666666668
+-28.0
+-27.25
+-26.5
+-25.5
+-24.583333333333332
+-24.5
+-24.383333333333333
+-24.0
+-22.0
+-20.0
+-19.25
+-19.0
+-18.25
+-17.0
+-16.0
+-15.0
+-14.5
+-12.033333333333333
+-11.666666666666666
+-11.0
+-10.0
+-9.0
+-8.0
+-7.0
+-6.0
+-4.0
+-3.25
+-1.75
+0.0
+1.5
+1.75
+2.0
+2.75
+3.0
+4.0
+4.5
+5.5
+6.0
+6.25
+7.0
+7.5
+8.0
+8.5
+9.916666666666666
+10.0
+11.0
+11.833333333333334
+12.0
+12.5
+12.833333333333334
+13.5
+14.0
+14.333333333333334
+15.0
+15.5
+15.75
+16.0
+16.166666666666668
+17.5
+18.0
+18.5
+19.0
+19.166666666666668
+19.5
+20.0
+20.5
+21.0
+21.083333333333332
+21.25
+21.5
+22.0
+22.833333333333332
+23.5
+23.75
+25.0
+25.5
+26.0
+27.0
+27.25
+27.5
+28.0
+28.5
+29.0
+30.0
+30.666666666666668
+30.75
+31.333333333333332
+32.0
+32.083333333333336
+33.0
+33.5
+34.0
+34.5
+35.0
+35.5
+36.0
+36.5
+36.75
+39.75
+40.0
+42.0
+43.5
+43.75
+44.0
+44.5
+45.0
+46.0
+47.0
+47.5
+48.0
+48.5
+50.0
+50.5
+51.5
+52.5
+52.75
+53.0
+54.0
+54.833333333333336
+55.0
+55.5
+56.0
+56.25
+57.0
+57.5
+58.0
+58.5
+59.083333333333336
+59.5
+60.0
+60.916666666666664
+61.5
+62.0
+63.0
+64.0
+66.0
+66.5
+67.0
+68.0
+70.0
+73.5
+75.0
+77.0
+80.0
+82.0
+85.0
+86.0
+86.16666666666667
+86.5
+88.0
+89.0
+89.5
+90.0

data/lib/astronoby/data/constellations/sorted_right_ascensions.dat

@@ -0,0 +1,237 @@
+0.06666666666666667
+0.14166666666666666
+0.16666666666666666
+0.3333333333333333
+0.7166666666666667
+0.75
+0.85
+0.8666666666666667
+1.1166666666666667
+1.3333333333333333
+1.3666666666666667
+1.4083333333333334
+1.5833333333333333
+1.6666666666666667
+1.7
+1.8333333333333333
+1.9083333333333334
+1.9166666666666667
+2.0
+2.0416666666666665
+2.1666666666666665
+2.3333333333333335
+2.4166666666666665
+2.433333333333333
+2.5166666666666666
+2.566666666666667
+2.65
+2.6666666666666665
+2.716666666666667
+3.0
+3.1
+3.1666666666666665
+3.2
+3.283333333333333
+3.3333333333333335
+3.3666666666666667
+3.4166666666666665
+3.5
+3.5083333333333333
+3.5833333333333335
+3.75
+3.8333333333333335
+3.8666666666666667
+4.0
+4.083333333333333
+4.266666666666667
+4.333333333333333
+4.5
+4.583333333333333
+4.616666666666666
+4.666666666666667
+4.691666666666666
+4.7
+4.75
+4.833333333333333
+4.916666666666667
+4.966666666666667
+5.0
+5.083333333333333
+5.333333333333333
+5.5
+5.6
+5.7
+5.766666666666667
+5.833333333333333
+5.883333333333334
+6.0
+6.1
+6.116666666666666
+6.166666666666667
+6.216666666666667
+6.241666666666666
+6.308333333333334
+6.5
+6.533333333333333
+6.583333333333333
+6.8
+6.833333333333333
+6.933333333333334
+7.0
+7.016666666666667
+7.2
+7.366666666666666
+7.5
+7.666666666666667
+7.75
+7.808333333333334
+7.883333333333334
+7.925
+7.966666666666667
+8.0
+8.083333333333334
+8.166666666666666
+8.366666666666667
+8.416666666666666
+8.45
+8.583333333333334
+8.833333333333334
+9.033333333333333
+9.083333333333334
+9.166666666666666
+9.25
+9.366666666666667
+9.583333333333334
+9.75
+9.883333333333333
+10.166666666666666
+10.25
+10.5
+10.583333333333334
+10.666666666666666
+10.75
+10.783333333333333
+10.833333333333334
+11.0
+11.25
+11.333333333333334
+11.5
+11.516666666666667
+11.833333333333334
+11.866666666666667
+12.0
+12.083333333333334
+12.25
+12.333333333333334
+12.583333333333334
+12.833333333333334
+13.0
+13.25
+13.5
+13.583333333333334
+13.666666666666666
+13.958333333333334
+14.0
+14.033333333333333
+14.166666666666666
+14.25
+14.416666666666666
+14.5
+14.533333333333333
+14.666666666666666
+14.75
+14.916666666666666
+15.05
+15.083333333333334
+15.166666666666666
+15.183333333333334
+15.25
+15.333333333333334
+15.433333333333334
+15.666666666666666
+15.75
+15.916666666666666
+16.0
+16.033333333333335
+16.083333333333332
+16.166666666666668
+16.266666666666666
+16.333333333333332
+16.375
+16.420833333333334
+16.533333333333335
+16.583333333333332
+16.75
+16.833333333333332
+17.0
+17.166666666666668
+17.25
+17.5
+17.583333333333332
+17.6
+17.666666666666668
+17.833333333333332
+17.966666666666665
+18.0
+18.175
+18.233333333333334
+18.25
+18.366666666666667
+18.425
+18.583333333333332
+18.662222222222223
+18.866666666666667
+19.0
+19.083333333333332
+19.166666666666668
+19.25
+19.258333333333333
+19.358333333333334
+19.4
+19.416666666666668
+19.666666666666668
+19.766666666666666
+19.833333333333332
+20.0
+20.141666666666666
+20.166666666666668
+20.25
+20.3
+20.333333333333332
+20.416666666666668
+20.533333333333335
+20.536666666666665
+20.566666666666666
+20.6
+20.666666666666668
+20.833333333333332
+20.875
+20.916666666666668
+21.0
+21.05
+21.116666666666667
+21.25
+21.333333333333332
+21.416666666666668
+21.466666666666665
+21.666666666666668
+21.733333333333334
+21.866666666666667
+21.875
+21.908333333333335
+21.966666666666665
+22.0
+22.133333333333333
+22.316666666666666
+22.75
+22.816666666666666
+22.866666666666667
+23.0
+23.166666666666668
+23.333333333333332
+23.5
+23.583333333333332
+23.75
+23.833333333333332
+23.916666666666668
+24.0

data/lib/astronoby/deflection.rb

@@ -0,0 +1,187 @@
+# frozen_string_literal: true
+
+module Astronoby
+  class Deflection
+    # Solar gravitational constant (m^3/s^2)
+    SOLAR_GRAVITATION_CONSTANT = 1.32712440017987e+20
+
+    # @param instant [Astronoby::Instant] the instant of the observation
+    # @param target_astrometric_position [Astronoby::Vector<Distance>] the
+    #   astrometric position of the target
+    # @param ephem [Astronoby::Ephemeris] the ephemeris to use for the
+    #   computation
+    def initialize(instant:, target_astrometric_position:, ephem:)
+      @instant = instant
+      @target_astrometric_position = target_astrometric_position
+      @ephem = ephem
+    end
+
+    # @return [Astronoby::Vector<Distance>] corrected position of the target
+    def corrected_position
+      Astronoby::Vector[
+        *(target_position + deflection_vector)
+          .map { |au| Astronoby::Distance.from_au(au) }
+      ]
+    end
+
+    private
+
+    def tt
+      @tt ||= @instant.tt
+    end
+
+    def earth_geometric
+      Earth.compute_geometric(ephem: @ephem, instant: @instant)
+    end
+
+    def sun_geometric
+      Sun.compute_geometric(ephem: @ephem, instant: @instant)
+    end
+
+    def observer_position
+      @observer_position ||= earth_geometric.position.map(&:au)
+    end
+
+    def target_position
+      @target_position ||= @target_astrometric_position.map(&:au)
+    end
+
+    def sun_position
+      @sun_position ||= sun_geometric.position.map(&:au)
+    end
+
+    def light_travel_time
+      @light_travel_time ||=
+        target_position.magnitude / Astronoby::Velocity.light_speed.aupd
+    end
+
+    def deflector_to_observer_position
+      sun_position - observer_position
+    end
+
+    def closest_deflector_to_target_position
+      observer_position + target_position - deflector_position_at_closest
+    end
+
+    def closest_deflector_to_observer_position
+      @closest_deflector_to_observer_position ||=
+        observer_position - deflector_position_at_closest
+    end
+
+    # Compute light-time difference for the point where light ray is closest to
+    # deflector
+    def closest_approach_time_diff
+      @closest_approach_time_diff ||=
+        Util::Maths.dot_product(
+          target_position / target_position.magnitude,
+          deflector_to_observer_position
+        ) / Astronoby::Velocity.light_speed.aupd
+    end
+
+    # Determine time when incoming photons were closest to the deflecting body
+    def time_at_closest_approach
+      @time_at_closest_approach ||= if closest_approach_time_diff > 0.0
+        tt - closest_approach_time_diff
+      elsif light_travel_time < closest_approach_time_diff
+        tt - light_travel_time
+      else
+        tt
+      end
+    end
+
+    # Get position of deflecting body at the time of closest approach
+    def deflector_position_at_closest
+      @deflector_position_at_closest ||=
+        Sun.compute_geometric(
+          ephem: @ephem,
+          instant: Instant.from_terrestrial_time(time_at_closest_approach)
+        ).position.map(&:au)
+    end
+
+    def observer_to_target_distance
+      @observer_to_target_distance ||= target_position.magnitude
+    end
+
+    def closest_deflector_to_target_distance
+      @closest_deflector_to_target_distance ||=
+        closest_deflector_to_target_position.magnitude
+    end
+
+    def closest_deflector_to_observer_distance
+      @closest_deflector_to_observer_distance ||=
+        closest_deflector_to_observer_position.magnitude
+    end
+
+    def observer_to_target_unit
+      @observer_to_target_unit ||=
+        target_position / observer_to_target_distance
+    end
+
+    def closest_deflector_to_target_unit
+      @closest_deflector_to_target_unit ||=
+        closest_deflector_to_target_position / closest_deflector_to_target_distance
+    end
+
+    def closest_deflector_to_observer_unit
+      @closest_deflector_to_observer_unit ||=
+        closest_deflector_to_observer_position / closest_deflector_to_observer_distance
+    end
+
+    def cos_angle_target
+      Util::Maths.dot_product(
+        observer_to_target_unit,
+        closest_deflector_to_target_unit
+      )
+    end
+
+    def cos_angle_deflector
+      Util::Maths.dot_product(
+        closest_deflector_to_target_unit,
+        closest_deflector_to_observer_unit
+      )
+    end
+
+    def cos_angle_observer
+      @cos_angle_observer ||= Util::Maths.dot_product(
+        closest_deflector_to_observer_unit,
+        observer_to_target_unit
+      )
+    end
+
+    # Calculate the relativistic deflection coefficient
+    # This implements Einstein's light deflection formula: α = 4GM/(c²r)
+    # where:
+    # - G is the gravitational constant
+    # - M is the mass of the deflecting body
+    # - c is the speed of light
+    # - r is the impact parameter (closest approach distance)
+    def relativistic_factor
+      2.0 * SOLAR_GRAVITATION_CONSTANT / (
+        Astronoby::Velocity.light_speed.mps *
+          Astronoby::Velocity.light_speed.mps *
+          closest_deflector_to_observer_distance *
+          Constants::ASTRONOMICAL_UNIT_IN_METERS
+      )
+    end
+
+    def geometry_factor
+      1.0 + cos_angle_deflector
+    end
+
+    def deflection_vector
+      return 0 if colinear?
+
+      relativistic_factor *
+        (
+          cos_angle_target * closest_deflector_to_observer_unit -
+            cos_angle_observer * closest_deflector_to_target_unit
+        ) / geometry_factor * observer_to_target_distance
+    end
+
+    # If deflector is nearly in line with target, make no correction
+    # (avoids numerical instability in nearly-collinear cases)
+    def colinear?
+      cos_angle_observer.abs > 0.99999999999
+    end
+  end
+end

data/lib/astronoby/distance.rb

@@ -25,6 +25,11 @@ module Astronoby
         from_meters(meters)
       end
       alias_method :from_au, :from_astronomical_units
+
+      def vector_from_meters(array)
+        Vector.elements(array.map { from_meters(_1) })
+      end
+      alias_method :vector_from_m, :vector_from_meters
     end
 
     attr_reader :meters
@@ -69,6 +74,10 @@ module Astronoby
       meters.zero?
     end
 
+    def abs2
+      meters**2
+    end
+
     def hash
       [meters, self.class].hash
     end

data/lib/astronoby/ephem.rb

@@ -0,0 +1,39 @@
+# frozen_string_literal: true
+
+require "ephem"
+
+module Astronoby
+  class Ephem
+    # Download an ephemeris file.
+    #
+    # @param name [String] Name of the ephemeris file, supported by the Ephem
+    #   gem
+    # @param target [String] Location where to store the file
+    # @return [Boolean] true if the download was successful, false otherwise
+    #
+    # @example Downloading de440t SPK from NASA JPL
+    #   Astronoby::Ephem.download(name: "de440t.bsp", target: "tmp/de440t.bsp")
+    def self.download(name:, target:)
+      ::Ephem::Download.call(name: name, target: target)
+    end
+
+    # Load an ephemeris file.
+    #
+    # @param target [String] Path of the ephemeris file
+    # @return [::Ephem::SPK] Ephemeris object from the Ephem gem
+    #
+    # @example Loading previously downloaded de440t SPK from NASA JPL
+    #   Astronoby::Ephem.load("tmp/de440t.bsp")
+    def self.load(target)
+      spk = ::Ephem::SPK.open(target)
+      unless ::Ephem::SPK::TYPES.include?(spk&.type)
+        raise(
+          EphemerisError,
+          "#{target} is not a valid type. Accepted: #{::Ephem::SPK::TYPES.join(", ")}"
+        )
+      end
+
+      spk
+    end
+  end
+end