astronoby 0.7.0 → 0.8.0

data/lib/astronoby.rb

@@ -1,12 +1,18 @@
 # frozen_string_literal: true
 
+require "astronoby/cache"
+require "astronoby/configuration"
 require "astronoby/constants"
 require "astronoby/angle"
 require "astronoby/angles/dms"
 require "astronoby/angles/hms"
+require "astronoby/constellation"
+require "astronoby/constellations/repository"
+require "astronoby/constellations/data"
+require "astronoby/constellations/finder"
 require "astronoby/distance"
 require "astronoby/ephem"
-require "astronoby/epoch"
+require "astronoby/julian_date"
 require "astronoby/instant"
 require "astronoby/vector"
 require "astronoby/velocity"
@@ -35,6 +41,7 @@ require "astronoby/events/rise_transit_set_event"
 require "astronoby/events/rise_transit_set_events"
 require "astronoby/events/rise_transit_set_calculator"
 require "astronoby/events/twilight_event"
+require "astronoby/events/twilight_events"
 require "astronoby/events/twilight_calculator"
 require "astronoby/geocentric_parallax"
 require "astronoby/mean_obliquity"

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: astronoby
 version: !ruby/object:Gem::Version
-  version: 0.7.0
+  version: 0.8.0
 platform: ruby
 authors:
 - Rémy Hannequin
@@ -37,6 +37,34 @@ dependencies:
     - - "~>"
       - !ruby/object:Gem::Version
         version: 0.4.2
+- !ruby/object:Gem::Dependency
+  name: benchmark
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '0.4'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '0.4'
+- !ruby/object:Gem::Dependency
+  name: irb
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.14'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.14'
 - !ruby/object:Gem::Dependency
   name: rake
   requirement: !ruby/object:Gem::Requirement
@@ -71,14 +99,14 @@ dependencies:
     requirements:
     - - "~>"
       - !ruby/object:Gem::Version
-        version: '2.3'
+        version: '3.0'
   type: :development
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
     - - "~>"
       - !ruby/object:Gem::Version
-        version: '2.3'
+        version: '3.0'
 - !ruby/object:Gem::Dependency
   name: simplecov
   requirement: !ruby/object:Gem::Requirement
@@ -120,16 +148,24 @@ files:
 - CHANGELOG.md
 - CODE_OF_CONDUCT.md
 - CONTRIBUTING.md
-- Gemfile
-- Gemfile.lock
 - LICENSE.txt
 - README.md
 - Rakefile
 - UPGRADING.md
-- benchmark/README.md
-- benchmark/benchmark.rb
-- benchmark/data/imcce.csv.zip
-- benchmark/data/sun_calc.csv.zip
+- docs/README.md
+- docs/angles.md
+- docs/celestial_bodies.md
+- docs/configuration.md
+- docs/coordinates.md
+- docs/ephem.md
+- docs/equinoxes_solstices_times.md
+- docs/glossary.md
+- docs/instant.md
+- docs/moon_phases.md
+- docs/observer.md
+- docs/reference_frames.md
+- docs/rise_transit_set_times.md
+- docs/twilight_times.md
 - lib/astronoby.rb
 - lib/astronoby/aberration.rb
 - lib/astronoby/angle.rb
@@ -147,15 +183,25 @@ files:
 - lib/astronoby/bodies/sun.rb
 - lib/astronoby/bodies/uranus.rb
 - lib/astronoby/bodies/venus.rb
+- lib/astronoby/cache.rb
+- lib/astronoby/configuration.rb
 - lib/astronoby/constants.rb
+- lib/astronoby/constellation.rb
+- lib/astronoby/constellations/data.rb
+- lib/astronoby/constellations/finder.rb
+- lib/astronoby/constellations/repository.rb
 - lib/astronoby/coordinates/ecliptic.rb
 - lib/astronoby/coordinates/equatorial.rb
 - lib/astronoby/coordinates/horizontal.rb
 - lib/astronoby/corrections/light_time_delay.rb
+- lib/astronoby/data/constellations/constellation_names.dat
+- lib/astronoby/data/constellations/indexed_abbreviations.dat
+- lib/astronoby/data/constellations/radec_to_index.dat
+- lib/astronoby/data/constellations/sorted_declinations.dat
+- lib/astronoby/data/constellations/sorted_right_ascensions.dat
 - lib/astronoby/deflection.rb
 - lib/astronoby/distance.rb
 - lib/astronoby/ephem.rb
-- lib/astronoby/epoch.rb
 - lib/astronoby/equinox_solstice.rb
 - lib/astronoby/errors.rb
 - lib/astronoby/events/moon_phases.rb
@@ -164,8 +210,10 @@ files:
 - lib/astronoby/events/rise_transit_set_events.rb
 - lib/astronoby/events/twilight_calculator.rb
 - lib/astronoby/events/twilight_event.rb
+- lib/astronoby/events/twilight_events.rb
 - lib/astronoby/geocentric_parallax.rb
 - lib/astronoby/instant.rb
+- lib/astronoby/julian_date.rb
 - lib/astronoby/mean_obliquity.rb
 - lib/astronoby/moon_phase.rb
 - lib/astronoby/nutation.rb
@@ -208,7 +256,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubygems_version: 3.6.7
+rubygems_version: 3.7.1
 specification_version: 4
 summary: Astronomical calculations
 test_files: []

data/Gemfile

@@ -1,5 +0,0 @@
-# frozen_string_literal: true
-
-source "https://rubygems.org"
-
-gemspec

data/Gemfile.lock

@@ -1,102 +0,0 @@
-PATH
-  remote: .
-  specs:
-    astronoby (0.7.0)
-      ephem (~> 0.3)
-      matrix (~> 0.4.2)
-
-GEM
-  remote: https://rubygems.org/
-  specs:
-    ast (2.4.3)
-    diff-lcs (1.5.1)
-    docile (1.4.1)
-    ephem (0.3.0)
-      minitar (~> 0.12)
-      numo-narray (~> 0.9.2.1)
-      zlib (~> 3.2)
-    json (2.10.2)
-    language_server-protocol (3.17.0.4)
-    lint_roller (1.1.0)
-    matrix (0.4.2)
-    minitar (0.12.1)
-    numo-narray (0.9.2.1)
-    parallel (1.26.3)
-    parser (3.3.7.4)
-      ast (~> 2.4.1)
-      racc
-    prism (1.4.0)
-    racc (1.8.1)
-    rainbow (3.1.1)
-    rake (13.2.1)
-    regexp_parser (2.10.0)
-    rspec (3.13.0)
-      rspec-core (~> 3.13.0)
-      rspec-expectations (~> 3.13.0)
-      rspec-mocks (~> 3.13.0)
-    rspec-core (3.13.2)
-      rspec-support (~> 3.13.0)
-    rspec-expectations (3.13.3)
-      diff-lcs (>= 1.2.0, < 2.0)
-      rspec-support (~> 3.13.0)
-    rspec-mocks (3.13.2)
-      diff-lcs (>= 1.2.0, < 2.0)
-      rspec-support (~> 3.13.0)
-    rspec-support (3.13.2)
-    rubocop (1.75.2)
-      json (~> 2.3)
-      language_server-protocol (~> 3.17.0.2)
-      lint_roller (~> 1.1.0)
-      parallel (~> 1.10)
-      parser (>= 3.3.0.2)
-      rainbow (>= 2.2.2, < 4.0)
-      regexp_parser (>= 2.9.3, < 3.0)
-      rubocop-ast (>= 1.44.0, < 2.0)
-      ruby-progressbar (~> 1.7)
-      unicode-display_width (>= 2.4.0, < 4.0)
-    rubocop-ast (1.44.0)
-      parser (>= 3.3.7.2)
-      prism (~> 1.4)
-    rubocop-performance (1.25.0)
-      lint_roller (~> 1.1)
-      rubocop (>= 1.75.0, < 2.0)
-      rubocop-ast (>= 1.38.0, < 2.0)
-    ruby-progressbar (1.13.0)
-    rubyzip (2.4.1)
-    simplecov (0.22.0)
-      docile (~> 1.1)
-      simplecov-html (~> 0.11)
-      simplecov_json_formatter (~> 0.1)
-    simplecov-html (0.13.1)
-    simplecov_json_formatter (0.1.4)
-    standard (1.49.0)
-      language_server-protocol (~> 3.17.0.2)
-      lint_roller (~> 1.0)
-      rubocop (~> 1.75.2)
-      standard-custom (~> 1.0.0)
-      standard-performance (~> 1.8)
-    standard-custom (1.0.2)
-      lint_roller (~> 1.0)
-      rubocop (~> 1.50)
-    standard-performance (1.8.0)
-      lint_roller (~> 1.1)
-      rubocop-performance (~> 1.25.0)
-    unicode-display_width (3.1.4)
-      unicode-emoji (~> 4.0, >= 4.0.4)
-    unicode-emoji (4.0.4)
-    zlib (3.2.1)
-
-PLATFORMS
-  ruby
-  x86_64-linux
-
-DEPENDENCIES
-  astronoby!
-  rake (~> 13.0)
-  rspec (~> 3.0)
-  rubyzip (~> 2.3)
-  simplecov (~> 0.22)
-  standard (~> 1.3)
-
-BUNDLED WITH
-   2.5.7

data/benchmark/README.md

@@ -1,131 +0,0 @@
-# Benchmark
-
-This is a first attempt to benchmark the accuracy of the library. It is not
-very scientific, but it gives a rough idea.
-
-## Method
-
-The goal is to answer these two questions:
-- Is the library accurate enough compared to a source of truth?
-- Is the library accurate enough compared with other Ruby libraries?
-
-The source of truth is the <abbr title="Institut de Mécanique
-Céleste et de Calcul des Éphémérides">IMCCE</abbr>, a French public
-institude attached to the Paris Observatory. Their ephemerides are used by
-governements and public institutions in multiple European countries, their
-precesion is among the highest in the world.
-They also provide web services to easily access all their data. Many thanks
-for providing such high accuracy data for free.
-
-The other Ruby library is [sun_calc](https://github.com/fishbrain/sun_calc).
-
-474,336 combinations of dates, latitudes and longitudes have been used to
-produce time predictions for the following events:
-- sunrise
-- sun's highest point
-- sunset
-- moonrise
-- moon's highest point
-- moonset
-
-For each combination, we first find out which of SunCalc or Astronoby is the
-closest to the IMCCE. Then we calculate the difference between Astronoby and
-the IMCCE, to discover if the difference is larger than the defined
-threshold of 5 minutes.
-
-## Results
-
-The following output has been generated using what will be part of version 0.6.
-
-```
-Unarchiving sun_calc.csv.zip...
-Done unarchiving sun_calc.csv.zip.
-Parsing sun_calc.csv...
-Done parsing sun_calc.csv.
-Unarchiving imcce.csv.zip...
-Done unarchiving imcce.csv.zip.
-Parsing imcce.csv...
-Done parsing imcce.csv.
-Comparing data...
-2024-01-01: Done.
-2024-01-02: Done.
-...
-2024-12-30: Done.
-2024-12-31: Done.
-Done comparing data.
-
-
-Sun rising time:
-astronoby: 295395 (62.28%)
-sun_calc: 99710 (21.02%)
-n/a: 79231 (16.7%)
-
-Sun transit time:
-astronoby: 434495 (91.6%)
-sun_calc: 39769 (8.38%)
-n/a: 72 (0.02%)
-
-Sun setting time:
-astronoby: 358428 (75.56%)
-n/a: 79231 (16.7%)
-sun_calc: 36677 (7.73%)
-
-Moon rising time:
-astronoby: 290866 (61.32%)
-n/a: 113815 (23.99%)
-sun_calc: 69655 (14.68%)
-
-Moon transit time:
-astronoby: 341902 (72.08%)
-n/a: 101916 (21.49%)
-sun_calc: 30518 (6.43%)
-
-Moon setting time:
-astronoby: 327099 (68.96%)
-n/a: 114308 (24.1%)
-sun_calc: 32929 (6.94%)
-
-Moon illuminated fraction:
-astronoby: 474336 (100.0%)
-
-Sun rising time too far:
-false: 452887 (95.48%)
-true: 21449 (4.52%)
-
-Sun transit time too far:
-false: 396617 (83.62%)
-true: 77719 (16.38%)
-
-Sun setting time too far:
-false: 453208 (95.55%)
-true: 21128 (4.45%)
-
-Moon rising time too far:
-false: 459044 (96.78%)
-true: 15292 (3.22%)
-
-Moon transit time too far:
-false: 384516 (81.06%)
-true: 89820 (18.94%)
-
-Moon setting time too far:
-false: 459222 (96.81%)
-true: 15114 (3.19%)
-```
-
-## Conclusion
-
-As we can see, Astronoby is more accurate than SunCalc in a vast majority of
-cases. When it comes to the Moon's illuminated fraction, Astronoby is always
-more accurate than SunCalc.
-
-`n/a` values means that at least one of the three sources don't have a value
-for the combination of date, latitude and longitude. This happens because
-the Moon and the Sun cannot always rise, transit and set everywhere on Earth
-every day of the year. Latitudes close to the poles are more likely to miss
-data.
-
-Astronoby can be considered "good enough" for more around 90% of the cases,
-which means there is still work to do if we want to always be less than 5
-minutes away from the what the IMCCE provides. We can notice that transit
-times are those that experience the most significant differences.

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/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