astronoby 0.6.0 → 0.8.0

data/benchmark/benchmark.rb

@@ -1,259 +0,0 @@
-require "astronoby"
-require "csv"
-require "zip"
-
-class Source
-  NAMES = [
-    ASTRONOBY = "astronoby",
-    IMCCE = "imcce",
-    SUN_CALC = "sun_calc"
-  ].freeze
-
-  attr_accessor :name,
-    :sun_rising_time,
-    :sun_transit_time,
-    :sun_setting_time,
-    :moon_rising_time,
-    :moon_transit_time,
-    :moon_setting_time,
-    :moon_illuminated_fraction
-end
-
-class Comparison
-  SUN_CALC = "sun_calc"
-  ASTRONOBY = "astronoby"
-  NON_APPLICABLE = "n/a"
-
-  TOO_FAR_THRESHOLD = 60 * 5 # 5 minutes
-
-  attr_accessor :sources, :truth
-
-  def initialize
-    @sources = []
-  end
-
-  %i[
-    sun_rising_time
-    sun_transit_time
-    sun_setting_time
-    moon_rising_time
-    moon_transit_time
-    moon_setting_time
-  ].each do |attribute|
-    define_method(:"closest_#{attribute}") do
-      compare(attribute)
-    end
-
-    define_method(:"#{attribute}_too_far?") do
-      too_far?(attribute)
-    end
-  end
-
-  def closest_moon_illuminated_fraction
-    compare(:moon_illuminated_fraction)
-  end
-
-  private
-
-  def compare(attribute)
-    unless truth.public_send(attribute) && sources.all? { |source| source.public_send(attribute) }
-      return NON_APPLICABLE
-    end
-
-    closest_source = sources.min_by do |source|
-      (truth.public_send(attribute) - source.public_send(attribute)).abs
-    end
-
-    closest_source.name
-  end
-
-  def too_far?(attribute)
-    truth_attribute = truth.public_send(attribute)
-    astronoby_attribute = sources
-      .find { _1.name == Source::ASTRONOBY }
-      .public_send(attribute)
-
-    return false unless truth_attribute && astronoby_attribute
-
-    (truth_attribute - astronoby_attribute).abs > TOO_FAR_THRESHOLD
-  end
-end
-
-class Result
-  def initialize
-    @sun_rising_time = []
-    @sun_transit_time = []
-    @sun_setting_time = []
-    @moon_rising_time = []
-    @moon_transit_time = []
-    @moon_setting_time = []
-    @illuminated_fraction = []
-    @sun_rising_time_too_far = []
-    @sun_transit_time_too_far = []
-    @sun_setting_time_too_far = []
-    @moon_rising_time_too_far = []
-    @moon_transit_time_too_far = []
-    @moon_setting_time_too_far = []
-  end
-
-  def add_comparison(comparison)
-    @sun_rising_time << comparison.closest_sun_rising_time
-    @sun_transit_time << comparison.closest_sun_transit_time
-    @sun_setting_time << comparison.closest_sun_setting_time
-    @moon_rising_time << comparison.closest_moon_rising_time
-    @moon_transit_time << comparison.closest_moon_transit_time
-    @moon_setting_time << comparison.closest_moon_setting_time
-    @illuminated_fraction << comparison.closest_moon_illuminated_fraction
-    @sun_rising_time_too_far << comparison.sun_rising_time_too_far?
-    @sun_transit_time_too_far << comparison.sun_transit_time_too_far?
-    @sun_setting_time_too_far << comparison.sun_setting_time_too_far?
-    @moon_rising_time_too_far << comparison.moon_rising_time_too_far?
-    @moon_transit_time_too_far << comparison.moon_transit_time_too_far?
-    @moon_setting_time_too_far << comparison.moon_setting_time_too_far?
-  end
-
-  def display
-    puts "Sun rising time:"
-    tally(@sun_rising_time)
-    puts "Sun transit time:"
-    tally(@sun_transit_time)
-    puts "Sun setting time:"
-    tally(@sun_setting_time)
-    puts "Moon rising time:"
-    tally(@moon_rising_time)
-    puts "Moon transit time:"
-    tally(@moon_transit_time)
-    puts "Moon setting time:"
-    tally(@moon_setting_time)
-    puts "Moon illuminated fraction:"
-    tally(@illuminated_fraction)
-    puts "Sun rising time too far:"
-    tally(@sun_rising_time_too_far)
-    puts "Sun transit time too far:"
-    tally(@sun_transit_time_too_far)
-    puts "Sun setting time too far:"
-    tally(@sun_setting_time_too_far)
-    puts "Moon rising time too far:"
-    tally(@moon_rising_time_too_far)
-    puts "Moon transit time too far:"
-    tally(@moon_transit_time_too_far)
-    puts "Moon setting time too far:"
-    tally(@moon_setting_time_too_far)
-  end
-
-  private
-
-  def tally(data)
-    t = data.tally
-    t.sort_by { |_key, value| -value }.each do |key, value|
-      puts "#{key}: #{value} (#{(value.to_f / t.values.sum * 100).round(2)}%)"
-    end
-    puts "\n"
-  end
-end
-
-data = {}
-result = Result.new
-
-sun_calc_zip_file = File.join(File.dirname(__FILE__), "data/sun_calc.csv.zip")
-imcce_zip_file = File.join(File.dirname(__FILE__), "data/imcce.csv.zip")
-
-puts "Unarchiving sun_calc.csv.zip..."
-
-Zip::File.open(sun_calc_zip_file) do |zip_file|
-  puts "Done unarchiving sun_calc.csv.zip."
-
-  csv_file = zip_file.find { |entry| entry.name.end_with?(".csv") }
-  break unless csv_file
-
-  puts "Parsing sun_calc.csv..."
-
-  csv_content = csv_file.get_input_stream.read
-  CSV.parse(csv_content, headers: true) do |row|
-    data[row["date"]] ||= {}
-    data[row["date"]][row["latitude"]] ||= {}
-    data[row["date"]][row["latitude"]][row["longitude"]] = Comparison.new.tap do |comparison|
-      source = Source.new.tap do |source|
-        source.name = Source::SUN_CALC
-        source.sun_rising_time = Time.new(row["sun_rising_time"]) if row["sun_rising_time"]
-        source.sun_transit_time = Time.new(row["sun_transit_time"]) if row["sun_transit_time"]
-        source.sun_setting_time = Time.new(row["sun_setting_time"]) if row["sun_setting_time"]
-        source.moon_rising_time = Time.new(row["moon_rising_time"]) if row["moon_rising_time"]
-        source.moon_transit_time = Time.new(row["moon_transit_time"]) if row["moon_transit_time"]
-        source.moon_setting_time = Time.new(row["moon_setting_time"]) if row["moon_setting_time"]
-        source.moon_illuminated_fraction = row["illuminated_fraction"].to_f
-      end
-      comparison.sources << source
-    end
-  end
-
-  puts "Done parsing sun_calc.csv."
-end
-
-puts "Unarchiving imcce.csv.zip..."
-
-Zip::File.open(imcce_zip_file) do |zip_file|
-  puts "Done unarchiving imcce.csv.zip."
-
-  csv_file = zip_file.find { |entry| entry.name.end_with?(".csv") }
-  break unless csv_file
-
-  puts "Parsing imcce.csv..."
-
-  csv_content = csv_file.get_input_stream.read
-  CSV.parse(csv_content, headers: true) do |row|
-    comparison = data[row["date"]][row["latitude"]][row["longitude"]]
-    comparison.truth = Source.new.tap do |source|
-      source.name = Source::IMCCE
-      source.sun_rising_time = Time.new(row["sun_rising_time"] + " UTC") if row["sun_rising_time"]
-      source.sun_transit_time = Time.new(row["sun_transit_time"] + " UTC") if row["sun_transit_time"]
-      source.sun_setting_time = Time.new(row["sun_setting_time"] + " UTC") if row["sun_setting_time"]
-      source.moon_rising_time = Time.new(row["moon_rising_time"] + " UTC") if row["moon_rising_time"]
-      source.moon_transit_time = Time.new(row["moon_transit_time"] + " UTC") if row["moon_transit_time"]
-      source.moon_setting_time = Time.new(row["moon_setting_time"] + " UTC") if row["moon_setting_time"]
-      source.moon_illuminated_fraction = row["illuminated_fraction"].to_f
-    end
-  end
-
-  puts "Done parsing imcce.csv."
-end
-
-puts "Comparing data..."
-
-data.each do |date, latitudes|
-  latitudes.each do |latitude, longitudes|
-    longitudes.each do |longitude, comparison|
-      noon = Time.new("#{date}T12:00:00Z")
-      observer = Astronoby::Observer.new(
-        latitude: Astronoby::Angle.from_degrees(latitude.to_i),
-        longitude: Astronoby::Angle.from_degrees(longitude.to_i)
-      )
-      sun = Astronoby::Sun.new(time: noon)
-      sun_observation_events = sun.observation_events(observer: observer)
-      moon = Astronoby::Moon.new(time: noon)
-      moon_observation_events = moon.observation_events(observer: observer)
-
-      source = Source.new.tap do |source|
-        source.name = Source::ASTRONOBY
-        source.sun_rising_time = sun_observation_events.rising_time
-        source.sun_transit_time = sun_observation_events.transit_time
-        source.sun_setting_time = sun_observation_events.setting_time
-        source.moon_rising_time = moon_observation_events.rising_time
-        source.moon_transit_time = moon_observation_events.transit_time
-        source.moon_setting_time = moon_observation_events.setting_time
-        source.moon_illuminated_fraction = moon.illuminated_fraction
-      end
-
-      comparison.sources << source
-      result.add_comparison(comparison)
-    end
-  end
-
-  puts "#{date}: Done."
-end
-
-puts "Done comparing data."
-puts
-puts
-
-puts result.display

data/lib/astronoby/astronomical_models/ephemeride_lunaire_parisienne.rb

@@ -1,143 +0,0 @@
-# frozen_string_literal: true
-
-module Astronoby
-  class EphemerideLunaireParisienne
-    # Lunar theory: Éphéméride Lunaire Parisienne
-    # By Jean Chapront, Michelle Chapront-Touzé
-    # https://en.wikipedia.org/wiki/Ephemeride_Lunaire_Parisienne
-
-    DEGREES_UNIT = 10**-6
-
-    # @return [Array] Periodic terms for the Moon's longitude and distance
-    def self.periodic_terms_for_moon_longitude_and_distance
-      [
-        [0, 0, 1, 0, 6288774, -20905355],
-        [2, 0, -1, 0, 1274027, -3699111],
-        [2, 0, 0, 0, 658314, -2955968],
-        [0, 0, 2, 0, 213618, -569925],
-        [0, 1, 0, 0, -185116, 48888],
-        [0, 0, 0, 2, -114332, -3149],
-        [2, 0, -2, 0, 58793, 246158],
-        [2, -1, -1, 0, 57066, -152138],
-        [2, 0, 1, 0, 53322, -170733],
-        [2, -1, 0, 0, 45758, -204586],
-        [0, 1, -1, 0, -40923, -129620],
-        [1, 0, 0, 0, -34720, 108743],
-        [0, 1, 1, 0, -30383, 104755],
-        [2, 0, 0, -2, 15327, 10321],
-        [0, 0, 1, 2, -12528, 0],
-        [0, 0, 1, -2, 10980, 79661],
-        [4, 0, -1, 0, 10675, -34782],
-        [0, 0, 3, 0, 10034, -23210],
-        [4, 0, -2, 0, 8548, -21636],
-        [2, 1, -1, 0, -7888, 24208],
-        [2, 1, 0, 0, -6766, 30824],
-        [1, 0, -1, 0, -5163, -8379],
-        [1, 1, 0, 0, 4987, -16675],
-        [2, -1, 1, 0, 4036, -12831],
-        [2, 0, 2, 0, 3994, -10445],
-        [4, 0, 0, 0, 3861, -11650],
-        [2, 0, -3, 0, 3665, 14403],
-        [0, 1, -2, 0, -2689, -7003],
-        [2, 0, -1, 2, -2602, 0],
-        [2, -1, -2, 0, 2390, 10056],
-        [1, 0, 1, 0, -2348, 6322],
-        [2, -2, 0, 0, 2236, -9884],
-        [0, 1, 2, 0, -2120, 5751],
-        [0, 2, 0, 0, -2069, 0],
-        [2, -2, -1, 0, 2048, -4950],
-        [2, 0, 1, -2, -1773, 4130],
-        [2, 0, 0, 2, -1595, 0],
-        [4, -1, -1, 0, 1215, -3958],
-        [0, 0, 2, 2, -1110, 0],
-        [3, 0, -1, 0, -892, 3258],
-        [2, 1, 1, 0, -810, 2616],
-        [4, -1, -2, 0, 759, -1897],
-        [0, 2, -1, 0, -713, -2117],
-        [2, 2, -1, 0, -700, 2354],
-        [2, 1, -2, 0, 691, 0],
-        [2, -1, 0, -2, 596, 0],
-        [4, 0, 1, 0, 549, -1423],
-        [0, 0, 4, 0, 537, -1117],
-        [4, -1, 0, 0, 520, -1571],
-        [1, 0, -2, 0, -487, -1739],
-        [2, 1, 0, -2, -399, 0],
-        [0, 0, 2, -2, -381, -4421],
-        [1, 1, 1, 0, 351, 0],
-        [3, 0, -2, 0, -340, 0],
-        [4, 0, -3, 0, 330, 0],
-        [2, -1, 2, 0, 327, 0],
-        [0, 2, 1, 0, -323, 1165],
-        [1, 1, -1, 0, 299, 0],
-        [2, 0, 3, 0, 294, 0],
-        [2, 0, -1, -2, 0, 8752]
-      ]
-    end
-
-    # @return [Array] Periodic terms for the Moon's latitude
-    def self.periodic_terms_for_moon_latitude
-      [
-        [0, 0, 0, 1, 5128122],
-        [0, 0, 1, 1, 280602],
-        [0, 0, 1, -1, 277693],
-        [2, 0, 0, -1, 173237],
-        [2, 0, -1, 1, 55413],
-        [2, 0, -1, -1, 46271],
-        [2, 0, 0, 1, 32573],
-        [0, 0, 2, 1, 17198],
-        [2, 0, 1, -1, 9266],
-        [0, 0, 2, -1, 8822],
-        [2, -1, 0, -1, 8216],
-        [2, 0, -2, -1, 4324],
-        [2, 0, 1, 1, 4200],
-        [2, 1, 0, -1, -3359],
-        [2, -1, -1, 1, 2463],
-        [2, -1, 0, 1, 2211],
-        [2, -1, -1, -1, 2065],
-        [0, 1, -1, -1, -1870],
-        [4, 0, -1, -1, 1828],
-        [0, 1, 0, 1, -1794],
-        [0, 0, 0, 3, -1749],
-        [0, 1, -1, 1, -1565],
-        [1, 0, 0, 1, -1491],
-        [0, 1, 1, 1, -1475],
-        [0, 1, 1, -1, -1410],
-        [0, 1, 0, -1, -1344],
-        [1, 0, 0, -1, -1335],
-        [0, 0, 3, 1, 1107],
-        [4, 0, 0, -1, 1021],
-        [4, 0, -1, 1, 833],
-        [0, 0, 1, -3, 777],
-        [4, 0, -2, 1, 671],
-        [2, 0, 0, -3, 607],
-        [2, 0, 2, -1, 596],
-        [2, -1, 1, -1, 491],
-        [2, 0, -2, 1, -451],
-        [0, 0, 3, -1, 439],
-        [2, 0, 2, 1, 422],
-        [2, 0, -3, -1, 421],
-        [2, 1, -1, 1, -366],
-        [2, 1, 0, 1, -351],
-        [4, 0, 0, 1, 331],
-        [2, -1, 1, 1, 315],
-        [2, -2, 0, -1, 302],
-        [0, 0, 1, 3, -283],
-        [2, 1, 1, -1, -229],
-        [1, 1, 0, -1, 223],
-        [1, 1, 0, 1, 223],
-        [0, 1, -2, -1, -220],
-        [2, 1, -1, -1, -220],
-        [1, 0, 1, 1, -185],
-        [2, -1, -2, -1, 181],
-        [0, 1, 2, 1, -177],
-        [4, 0, -2, -1, 176],
-        [4, -1, -1, -1, 166],
-        [1, 0, 1, -1, -164],
-        [4, 0, 1, -1, 132],
-        [1, 0, -1, -1, -119],
-        [4, -1, 0, -1, 115],
-        [2, -2, 0, 1, 107]
-      ]
-    end
-  end
-end

data/lib/astronoby/epoch.rb

@@ -1,22 +0,0 @@
-# frozen_string_literal: true
-
-module Astronoby
-  class Epoch
-    B1900 = 2415020.3135
-    J1900 = 2415020.0
-    B1950 = 2433282.4235
-    J1950 = 2433282.5
-    J2000 = 2451545.0
-
-    DEFAULT_EPOCH = J2000
-    JULIAN_DAY_NUMBER_OFFSET = 0.5
-
-    def self.from_time(time)
-      time.to_datetime.ajd
-    end
-
-    def self.to_utc(epoch)
-      DateTime.jd(epoch + JULIAN_DAY_NUMBER_OFFSET).to_time.utc
-    end
-  end
-end