astronoby 0.7.0 → 0.9.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/distance.rb

@@ -26,6 +26,12 @@ module Astronoby
       end
       alias_method :from_au, :from_astronomical_units
 
+      def from_parsecs(parsecs)
+        meters = parsecs * Constants::PARSEC_IN_METERS
+        from_meters(meters)
+      end
+      alias_method :from_pc, :from_parsecs
+
       def vector_from_meters(array)
         Vector.elements(array.map { from_meters(_1) })
       end

data/lib/astronoby/equinox_solstice.rb

@@ -110,7 +110,7 @@ module Astronoby
     private
 
     def compute
-      t = (julian_day - Epoch::J2000) / Constants::DAYS_PER_JULIAN_CENTURY
+      t = (julian_day - JulianDate::J2000) / Constants::DAYS_PER_JULIAN_CENTURY
       w = Angle.from_degrees(35999.373 * t) - Angle.from_degrees(2.47)
       delta = 1 +
         0.0334 * w.cos +
@@ -123,7 +123,7 @@ module Astronoby
       delta_days = 0.00001 * s / delta
       epoch = julian_day + delta_days
 
-      Epoch.to_utc(epoch)
+      Instant.from_terrestrial_time(epoch).to_time
     end
 
     def julian_day

data/lib/astronoby/events/extremum_calculator.rb

@@ -0,0 +1,233 @@
+# frozen_string_literal: true
+
+module Astronoby
+  # Calculates extrema (minima and maxima) in the distance between two
+  # celestial bodies over a given time range using adaptive sampling and golden
+  # section search refinement.
+  class ExtremumCalculator
+    # Mathematical constants
+    PHI = (1 + Math.sqrt(5)) / 2
+    INVPHI = 1 / PHI
+    GOLDEN_SECTION_TOLERANCE = 1e-5
+
+    # Algorithm parameters
+    MIN_SAMPLES_PER_PERIOD = 20
+    DUPLICATE_THRESHOLD_DAYS = 0.5
+    BOUNDARY_BUFFER_DAYS = 0.01
+
+    # Orbital periods
+    ORBITAL_PERIODS = {
+      "Astronoby::Moon" => 27.504339,
+      "Astronoby::Mercury" => 87.969,
+      "Astronoby::Venus" => 224.701,
+      "Astronoby::Earth" => 365.256,
+      "Astronoby::Mars" => 686.98,
+      "Astronoby::Jupiter" => 4332.59,
+      "Astronoby::Saturn" => 10759.22,
+      "Astronoby::Uranus" => 30688.5,
+      "Astronoby::Neptune" => 60182.0
+    }.freeze
+
+    # @param ephem [::Ephem::SPK] Ephemeris data source
+    # @param body [Astronoby::SolarSystemBody] The celestial body to track
+    # @param primary_body [Astronoby::SolarSystemBody] The reference body
+    #   (e.g., Sun for planetary orbits)
+    # @param samples_per_period [Integer] Number of samples to take per orbital
+    #   period
+    def initialize(
+      body:,
+      primary_body:,
+      ephem:,
+      samples_per_period: 60
+    )
+      @ephem = ephem
+      @body = body
+      @primary_body = primary_body
+      @orbital_period = ORBITAL_PERIODS.fetch(body.name)
+      @samples_per_period = samples_per_period
+    end
+
+    # Finds all apoapsis events between two times
+    # @param start_time [Time] Start time
+    # @param end_time [Time] End time
+    # @return [Array<Astronoby::ExtremumEvent>] Array of apoapsis events
+    def apoapsis_events_between(start_time, end_time)
+      find_extrema(
+        Astronoby::Instant.from_time(start_time).tt,
+        Astronoby::Instant.from_time(end_time).tt,
+        type: :maximum
+      )
+    end
+
+    # Finds all periapsis events between two times
+    # @param start_time [Time] Start time
+    # @param end_time [Time] End time
+    # @return [Array<Astronoby::ExtremumEvent>] Array of periapsis events
+    def periapsis_events_between(start_time, end_time)
+      find_extrema(
+        Astronoby::Instant.from_time(start_time).tt,
+        Astronoby::Instant.from_time(end_time).tt,
+        type: :minimum
+      )
+    end
+
+    # Finds extrema (minima or maxima) in the distance between bodies within
+    #   a time range
+    # @param start_jd [Float] Start time in Julian Date (Terrestrial Time)
+    # @param end_jd [Float] End time in Julian Date (Terrestrial Time)
+    # @param type [Symbol] :maximum or :minimum
+    # @return [Array<Astronoby::ExtremumEvent>] Array of extrema events
+    def find_extrema(start_jd, end_jd, type: :maximum)
+      # 1: Find extrema candidates through adaptive sampling
+      candidates = find_extrema_candidates(start_jd, end_jd, type)
+
+      # 2: Refine each candidate using golden section search
+      refined_extrema = candidates
+        .map { |candidate| refine_extremum(candidate, type) }
+        .compact
+
+      # 3: Remove duplicates and boundary artifacts
+      refined_extrema = remove_duplicates(refined_extrema)
+      filter_boundary_artifacts(refined_extrema, start_jd, end_jd)
+    end
+
+    private
+
+    def distance_at(jd)
+      instant = Instant.from_terrestrial_time(jd)
+      body_geometric = @body.geometric(ephem: @ephem, instant: instant)
+      primary_geometric = @primary_body
+        .geometric(ephem: @ephem, instant: instant)
+
+      distance_vector = body_geometric.position - primary_geometric.position
+      distance_vector.magnitude
+    end
+
+    def find_extrema_candidates(start_jd, end_jd, type)
+      samples = collect_samples(start_jd, end_jd)
+      find_local_extrema_in_samples(samples, type)
+    end
+
+    def collect_samples(start_jd, end_jd)
+      duration = end_jd - start_jd
+      sample_count = calculate_sample_count(duration)
+      step = duration / sample_count
+
+      (0..sample_count).map do |i|
+        jd = start_jd + (i * step)
+        {jd: jd, value: distance_at(jd)}
+      end
+    end
+
+    def calculate_sample_count(duration)
+      # Adaptive sampling: scale with duration and orbital period
+      periods_in_range = duration / @orbital_period
+      base_samples = (periods_in_range * @samples_per_period).to_i
+
+      # Ensure minimum sample density for short ranges
+      [base_samples, MIN_SAMPLES_PER_PERIOD].max
+    end
+
+    def find_local_extrema_in_samples(samples, type)
+      candidates = []
+
+      # Check each interior point for local extrema
+      (1...samples.length - 1).each do |i|
+        if local_extremum?(samples, i, type)
+          candidates << {
+            start_jd: samples[i - 1][:jd],
+            end_jd: samples[i + 1][:jd],
+            center_jd: samples[i][:jd]
+          }
+        end
+      end
+
+      candidates
+    end
+
+    def local_extremum?(samples, index, type)
+      current_val = samples[index][:value].m
+      prev_val = samples[index - 1][:value].m
+      next_val = samples[index + 1][:value].m
+
+      if type == :maximum
+        current_val > prev_val && current_val > next_val
+      else
+        current_val < prev_val && current_val < next_val
+      end
+    end
+
+    def refine_extremum(candidate, type)
+      golden_section_search(candidate[:start_jd], candidate[:end_jd], type)
+    end
+
+    def golden_section_search(a, b, type)
+      return nil if b <= a
+
+      tol = GOLDEN_SECTION_TOLERANCE * (b - a).abs
+
+      # Initial points using golden ratio
+      x1 = a + (1 - INVPHI) * (b - a)
+      x2 = a + INVPHI * (b - a)
+
+      f1 = distance_at(x1).m
+      f2 = distance_at(x2).m
+
+      while (b - a).abs > tol
+        should_keep_left = (type == :maximum) ? (f1 > f2) : (f1 < f2)
+
+        if should_keep_left
+          b = x2
+          x2 = x1
+          f2 = f1
+          x1 = a + (1 - INVPHI) * (b - a)
+          f1 = distance_at(x1).m
+        else
+          a = x1
+          x1 = x2
+          f1 = f2
+          x2 = a + INVPHI * (b - a)
+          f2 = distance_at(x2).m
+        end
+      end
+
+      mid = (a + b) / 2
+      ExtremumEvent.new(
+        Instant.from_terrestrial_time(mid),
+        distance_at(mid)
+      )
+    end
+
+    def remove_duplicates(extrema)
+      return extrema if extrema.length <= 1
+
+      cleaned = [extrema.first]
+
+      extrema.each_with_index do |current, i|
+        next if i == 0
+
+        is_duplicate = cleaned.any? do |existing|
+          time_diff = (
+            current.instant.tt - existing.instant.tt
+          ).abs
+          time_diff < DUPLICATE_THRESHOLD_DAYS
+        end
+
+        cleaned << current unless is_duplicate
+      end
+
+      cleaned
+    end
+
+    def filter_boundary_artifacts(extrema, start_jd, end_jd)
+      extrema.reject do |extreme|
+        start_diff = (extreme.instant.tt - start_jd).abs
+        end_diff = (extreme.instant.tt - end_jd).abs
+
+        too_close_to_start = start_diff < BOUNDARY_BUFFER_DAYS
+        too_close_to_end = end_diff < BOUNDARY_BUFFER_DAYS
+        too_close_to_start || too_close_to_end
+      end
+    end
+  end
+end

data/lib/astronoby/events/extremum_event.rb

@@ -0,0 +1,15 @@
+# frozen_string_literal: true
+
+module Astronoby
+  # Represents an extremum event with its timing and value
+  class ExtremumEvent
+    attr_reader :instant, :value
+
+    # @param instant [Astronoby::Instant] When the event occurs
+    # @param value [Object] The extreme value
+    def initialize(instant, value)
+      @instant = instant
+      @value = value
+    end
+  end
+end