astronoby 0.6.0 → 0.7.0

data/lib/astronoby/events/rise_transit_set_iteration.rb

@@ -1,218 +0,0 @@
-# frozen_string_literal: true
-
-module Astronoby
-  class RiseTransitSetIteration
-    EARTH_SIDEREAL_ROTATION_RATE = 360.98564736629
-
-    # Source:
-    #  Title: Astronomical Algorithms
-    #  Author: Jean Meeus
-    #  Edition: 2nd edition
-    #  Chapter: 15 - Rising, Transit, and Setting
-
-    # @param observer [Astronoby::Observer] Observer
-    # @param date [Date] Date of the event
-    # @param coordinates_of_the_previous_day [Astronoby::Coordinates::Equatorial]
-    #  Coordinates of the body of the previous day
-    # @param coordinates_of_the_day [Astronoby::Coordinates::Equatorial]
-    #   Coordinates of the body of the day
-    # @param coordinates_of_the_next_day [Astronoby::Coordinates::Equatorial]
-    #   Coordinates of the body of the next day
-    # @param shift [Astronoby::Angle] Altitude shift
-    # @param initial_rising [Float] Initial rising
-    # @param initial_transit [Float] Initial transit
-    # @param initial_setting [Float] Initial setting
-    def initialize(
-      observer:,
-      date:,
-      coordinates_of_the_previous_day:,
-      coordinates_of_the_day:,
-      coordinates_of_the_next_day:,
-      shift:,
-      initial_rising:,
-      initial_transit:,
-      initial_setting:
-    )
-      @observer = observer
-      @date = date
-      @coordinates_of_the_previous_day = coordinates_of_the_previous_day
-      @coordinates_of_the_day = coordinates_of_the_day
-      @coordinates_of_the_next_day = coordinates_of_the_next_day
-      @shift = shift
-      @initial_rising = initial_rising
-      @initial_transit = initial_transit
-      @initial_setting = initial_setting
-    end
-
-    # @return [Array<Float>] Iteration results
-    def iterate
-      [
-        delta_m_rising,
-        delta_m_transit,
-        delta_m_setting
-      ]
-    end
-
-    private
-
-    def delta_m_rising
-      (local_horizontal_altitude_rising - @shift).degrees./(
-        Constants::DEGREES_PER_CIRCLE *
-          declination_rising.cos *
-          @observer.latitude.cos *
-          local_hour_angle_rising.sin
-      )
-    end
-
-    def delta_m_transit
-      -local_hour_angle_transit.degrees / Constants::DEGREES_PER_CIRCLE
-    end
-
-    def delta_m_setting
-      (local_horizontal_altitude_setting - @shift).degrees./(
-        Constants::DEGREES_PER_CIRCLE *
-          declination_setting.cos *
-          @observer.latitude.cos *
-          local_hour_angle_setting.sin
-      )
-    end
-
-    def observer_longitude
-      # Longitude must be treated positively westwards from the meridian of
-      # Greenwich, and negatively to the east
-      @observer_longitude ||= -@observer.longitude
-    end
-
-    def apparent_gst_at_midnight
-      @apparent_gst_at_midnight ||= Angle.from_hours(
-        GreenwichSiderealTime.from_utc(
-          Time.utc(@date.year, @date.month, @date.day)
-        ).time
-      )
-    end
-
-    def gst_rising
-      @gst_rising ||= Angle.from_degrees(
-        apparent_gst_at_midnight.degrees +
-          EARTH_SIDEREAL_ROTATION_RATE * @initial_rising
-      )
-    end
-
-    def gst_transit
-      @gst_transit ||= Angle.from_degrees(
-        apparent_gst_at_midnight.degrees +
-          EARTH_SIDEREAL_ROTATION_RATE * @initial_transit
-      )
-    end
-
-    def gst_setting
-      @gst_setting ||= Angle.from_degrees(
-        apparent_gst_at_midnight.degrees +
-          EARTH_SIDEREAL_ROTATION_RATE * @initial_setting
-      )
-    end
-
-    def local_hour_angle_rising
-      @local_hour_angle_rising ||=
-        gst_rising - observer_longitude - right_ascension_rising
-    end
-
-    def local_hour_angle_transit
-      gst_transit - observer_longitude - right_ascension_transit
-    end
-
-    def local_hour_angle_setting
-      @local_hour_angle_setting ||=
-        gst_setting - observer_longitude - right_ascension_setting
-    end
-
-    def local_horizontal_altitude_rising
-      Angle.asin(
-        @observer.latitude.sin * declination_rising.sin +
-          @observer.latitude.cos * declination_rising.cos * local_hour_angle_rising.cos
-      )
-    end
-
-    def local_horizontal_altitude_setting
-      Angle.asin(
-        @observer.latitude.sin * declination_setting.sin +
-          @observer.latitude.cos * declination_setting.cos * local_hour_angle_setting.cos
-      )
-    end
-
-    def right_ascension_rising
-      Angle.from_degrees(
-        Util::Maths.interpolate(
-          Util::Maths.normalize_angles_for_interpolation(
-            [
-              @coordinates_of_the_previous_day.right_ascension.degrees,
-              @coordinates_of_the_day.right_ascension.degrees,
-              @coordinates_of_the_next_day.right_ascension.degrees
-            ]
-          ),
-          @initial_rising
-        )
-      )
-    end
-
-    def right_ascension_transit
-      Angle.from_degrees(
-        Util::Maths.interpolate(
-          Util::Maths.normalize_angles_for_interpolation(
-            [
-              @coordinates_of_the_previous_day.right_ascension.degrees,
-              @coordinates_of_the_day.right_ascension.degrees,
-              @coordinates_of_the_next_day.right_ascension.degrees
-            ]
-          ),
-          @initial_transit
-        )
-      )
-    end
-
-    def right_ascension_setting
-      Angle.from_degrees(
-        Util::Maths.interpolate(
-          Util::Maths.normalize_angles_for_interpolation(
-            [
-              @coordinates_of_the_previous_day.right_ascension.degrees,
-              @coordinates_of_the_day.right_ascension.degrees,
-              @coordinates_of_the_next_day.right_ascension.degrees
-            ]
-          ),
-          @initial_setting
-        )
-      )
-    end
-
-    def declination_rising
-      Angle.from_degrees(
-        Util::Maths.interpolate(
-          Util::Maths.normalize_angles_for_interpolation(
-            [
-              @coordinates_of_the_previous_day.declination.degrees,
-              @coordinates_of_the_day.declination.degrees,
-              @coordinates_of_the_next_day.declination.degrees
-            ]
-          ),
-          @initial_rising
-        )
-      )
-    end
-
-    def declination_setting
-      Angle.from_degrees(
-        Util::Maths.interpolate(
-          Util::Maths.normalize_angles_for_interpolation(
-            [
-              @coordinates_of_the_previous_day.declination.degrees,
-              @coordinates_of_the_day.declination.degrees,
-              @coordinates_of_the_next_day.declination.degrees
-            ]
-          ),
-          @initial_setting
-        )
-      )
-    end
-  end
-end

data/lib/astronoby/events/twilight_events.rb

@@ -1,121 +0,0 @@
-# frozen_string_literal: true
-
-module Astronoby
-  module Events
-    class TwilightEvents
-      TWILIGHTS = [
-        CIVIL = :civil,
-        NAUTICAL = :nautical,
-        ASTRONOMICAL = :astronomical
-      ].freeze
-
-      TWILIGHT_ANGLES = {
-        CIVIL => Angle.from_degrees(96),
-        NAUTICAL => Angle.from_degrees(102),
-        ASTRONOMICAL => Angle.from_degrees(108)
-      }.freeze
-
-      PERIODS_OF_THE_DAY = [
-        MORNING = :morning,
-        EVENING = :evening
-      ].freeze
-
-      attr_reader :morning_civil_twilight_time,
-        :evening_civil_twilight_time,
-        :morning_nautical_twilight_time,
-        :evening_nautical_twilight_time,
-        :morning_astronomical_twilight_time,
-        :evening_astronomical_twilight_time
-
-      def initialize(observer:, sun:)
-        @observer = observer
-        @sun = sun
-        PERIODS_OF_THE_DAY.each do |period_of_the_day|
-          TWILIGHT_ANGLES.each do |twilight, _|
-            zenith_angle = TWILIGHT_ANGLES[twilight]
-            instance_variable_set(
-              :"@#{period_of_the_day}_#{twilight}_twilight_time",
-              compute(period_of_the_day, zenith_angle)
-            )
-          end
-        end
-      end
-
-      # @param period_of_the_day [Symbol] :morning or :evening
-      # @param zenith_angle [Angle] The zenith angle of the twilight
-      def time_for_zenith_angle(period_of_the_day:, zenith_angle:)
-        unless PERIODS_OF_THE_DAY.include?(period_of_the_day)
-          raise IncompatibleArgumentsError,
-            "Only #{PERIODS_OF_THE_DAY.join(" or ")} are allowed as period_of_the_day, got #{period_of_the_day}"
-        end
-
-        compute(period_of_the_day, zenith_angle)
-      end
-
-      private
-
-      # Source:
-      #  Title: Practical Astronomy with your Calculator or Spreadsheet
-      #  Authors: Peter Duffett-Smith and Jonathan Zwart
-      #  Edition: Cambridge University Press
-      #  Chapter: 50 - Twilight
-      def compute(period_of_the_day, zenith_angle)
-        period_time = if period_of_the_day == MORNING
-          observation_events.rising_time
-        else
-          observation_events.setting_time
-        end
-
-        hour_angle_at_period = equatorial_coordinates_at_midday
-          .compute_hour_angle(time: period_time, longitude: @observer.longitude)
-
-        term1 = zenith_angle.cos -
-          @observer.latitude.sin *
-            @equatorial_coordinates_at_midday.declination.sin
-        term2 = @observer.latitude.cos *
-          equatorial_coordinates_at_midday.declination.cos
-        hour_angle_ratio_at_twilight = term1 / term2
-
-        return unless hour_angle_ratio_at_twilight.between?(-1, 1)
-
-        hour_angle_at_twilight = Angle.acos(hour_angle_ratio_at_twilight)
-        time_sign = -1
-
-        if period_of_the_day == MORNING
-          hour_angle_at_twilight = Angle.from_degrees(
-            Constants::DEGREES_PER_CIRCLE - hour_angle_at_twilight.degrees
-          )
-          time_sign = 1
-        end
-
-        twilight_in_hours =
-          time_sign * (hour_angle_at_twilight - hour_angle_at_period).hours *
-          GreenwichSiderealTime::SIDEREAL_MINUTE_IN_UT_MINUTE
-        twilight_in_seconds = time_sign *
-          twilight_in_hours *
-          Constants::SECONDS_PER_HOUR
-
-        (period_time + twilight_in_seconds).round
-      end
-
-      def observation_events
-        @observation_events ||= @sun.observation_events(observer: @observer)
-      end
-
-      def midday
-        date = @sun.time.to_date
-        Time.utc(date.year, date.month, date.day, 12)
-      end
-
-      def sun_at_midday
-        Sun.new(time: midday)
-      end
-
-      def equatorial_coordinates_at_midday
-        @equatorial_coordinates_at_midday ||= sun_at_midday
-          .apparent_ecliptic_coordinates
-          .to_apparent_equatorial(epoch: Epoch.from_time(midday))
-      end
-    end
-  end
-end

data/lib/astronoby/util/astrodynamics.rb

@@ -1,60 +0,0 @@
-# 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.from_radians(solution)
-          end
-
-          eccentric_anomaly_newton_raphson(
-            mean_anomaly,
-            orbital_eccentricity,
-            precision,
-            maximum_iteration_count,
-            current_iteration + 1,
-            Angle.from_radians(solution)
-          )
-        end
-      end
-    end
-  end
-end