astronoby 0.3.0 → 0.4.0

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 2ab8b9191b1b2452c678389fbc4c1750c4080c8f3cddf55cf7aebad2d99a0da7
-  data.tar.gz: f823456f556792e1d2ea3747f3665c84442a907c3a08271df5f871204bd09f10
+  metadata.gz: 3c65785b0c459ba8add86fcceb906205caa79a80f47b7f33c0b7ee008173fd14
+  data.tar.gz: 66478c6055d84e34f560490f225ae3e3d5020972e729b37e8f6f0a6566ae173e
 SHA512:
-  metadata.gz: b14d2a19972991fecc5b0a12c15fe18f3a9a07d0e9fc705fe2d984af14e45ae82962332361e65dace20c27b5f6e05e19483385731770abe463a3bf28adcfb354
-  data.tar.gz: 725810244872e539789b23d2b28ce656781a1f9121fabc76c3e2ebc5c9e6616900bff852290cf8719a10feb7a06afa42f7d02f4e00ac3ba2dadc4dea47d6a24c
+  metadata.gz: 4bae651aabaec62d258a768952eb3cd830d99639acf33f02b9f0ab6fd2eeff8215d9a30087a5818736b777c1a86344918f98f33f987f65829a8712cb20477c55
+  data.tar.gz: 401659e2d49599246ad2ea1a455a6d60b9797c49e730cacf8a614f4235929263c9318644373a96ea496969fcb814cc93144a8dbfffd0a0318e1628500bb1c8e8

data/CHANGELOG.md

@@ -1,5 +1,77 @@
 # Changelog
 
+## 0.4.0 - 2024-04-29
+
+_If you are upgrading: please see [UPGRADING.md]._
+
+[UPGRADING.md]: https://github.com/rhannequin/astronoby/blob/main/UPGRADING.md
+
+### Bug fixes
+
+* Fix ecliptic to equatorial epoch ([#56])
+
+[#56]: https://github.com/rhannequin/astronoby/pull/56
+
+### Features
+
+* Add twilight times ([#49])
+* Add interpolation method ([#52])
+* Add decimal_hour_to_time util ([#53])
+* Calculate leap seconds for an instant ([#54])
+* Add `Angle#-@` ([#55])
+* Enable equivalence and hash equality to `Observer` ([#57])
+* Twilight events dedicated class ([#61])
+
+[#49]: https://github.com/rhannequin/astronoby/pull/49
+[#52]: https://github.com/rhannequin/astronoby/pull/52
+[#53]: https://github.com/rhannequin/astronoby/pull/53
+[#54]: https://github.com/rhannequin/astronoby/pull/54
+[#55]: https://github.com/rhannequin/astronoby/pull/55
+[#57]: https://github.com/rhannequin/astronoby/pull/57
+[#61]: https://github.com/rhannequin/astronoby/pull/61
+
+### Improvements
+
+* Upgrade bundler from 2.3.11 to 2.5.7 by @dorianmariecom ([#45])
+* Drop `BigDecimal` ([#46])
+* Bump rake from 13.1.0 to 13.2.0 ([#47])
+* Increase Ruby versions support ([#48])
+* Bump rake from 13.2.0 to 13.2.1 ([#51])
+* Dedicated constants class ([#62])
+* Improve accuracy of equation of time ([#63])
+* Twilight times better accuracy ([#65])
+* Update UPGRADING.md ([#66])
+* release: Bump version to 0.4.0 ([#67])
+
+[#45]: https://github.com/rhannequin/astronoby/pull/45
+[#46]: https://github.com/rhannequin/astronoby/pull/46
+[#47]: https://github.com/rhannequin/astronoby/pull/47
+[#48]: https://github.com/rhannequin/astronoby/pull/48
+[#51]: https://github.com/rhannequin/astronoby/pull/51
+[#62]: https://github.com/rhannequin/astronoby/pull/62
+[#63]: https://github.com/rhannequin/astronoby/pull/63
+[#65]: https://github.com/rhannequin/astronoby/pull/65
+[#66]: https://github.com/rhannequin/astronoby/pull/66
+[#67]: https://github.com/rhannequin/astronoby/pull/67
+
+### Backward-incompatible changes
+
+* More accurate rising, transit and setting times ([#50])
+* Observation events dedicated and centralized class ([#60])
+* Change `Astronoby::Sun` constructor ([#64])
+
+[#50]: https://github.com/rhannequin/astronoby/pull/50
+[#60]: https://github.com/rhannequin/astronoby/pull/60
+[#64]: https://github.com/rhannequin/astronoby/pull/64
+
+### New Contributors
+
+* @dorianmariecom made their first contribution in [#45]
+
+[#45]: https://github.com/rhannequin/astronoby/pull/45
+
+**Full Changelog**: https://github.com/rhannequin/astronoby/compare/v0.3.0...v0.4.0
+
 ## 0.3.0 - 2024-03-29
 
 _If you are upgrading: please see [`UPGRADING.md`]._

data/Gemfile.lock

@@ -1,7 +1,7 @@
 PATH
   remote: .
   specs:
-    astronoby (0.3.0)
+    astronoby (0.4.0)
       matrix (~> 0.4.2)
       rake (~> 13.0)
       rspec (~> 3.0)
@@ -21,7 +21,7 @@ GEM
       racc
     racc (1.7.3)
     rainbow (3.1.1)
-    rake (13.1.0)
+    rake (13.2.1)
     regexp_parser (2.9.0)
     rexml (3.2.6)
     rspec (3.13.0)
@@ -77,4 +77,4 @@ DEPENDENCIES
   standard (~> 1.3)
 
 BUNDLED WITH
-   2.3.11
+   2.5.7

data/README.md

@@ -32,7 +32,7 @@ major one is released.
 
 ```rb
 angle1 = Astronoby::Angle.from_degrees(90)
-angle2 = Astronoby::Angle.from_radians(Astronoby::Angle::PI / 2)
+angle2 = Astronoby::Angle.from_radians(Math::PI / 2)
 angle3 = Astronoby::Angle.from_hours(12)
 
 angle1 == angle2
@@ -71,47 +71,82 @@ horizontal.azimuth.str(:dms)
 
 ```rb
 time = Time.utc(2023, 2, 17, 11, 0, 0)
-epoch = Astronoby::Epoch.from_time(time)
 
 latitude = Astronoby::Angle.from_degrees(48.8566)
 longitude = Astronoby::Angle.from_degrees(2.3522)
 
-sun = Astronoby::Sun.new(epoch: epoch)
+sun = Astronoby::Sun.new(time: time)
 
 horizontal_coordinates = sun.horizontal_coordinates(
   latitude: latitude,
   longitude: longitude
 )
 
-horizontal_coordinates.altitude.degrees.to_f
-# => 27.50008242057459
+horizontal_coordinates.altitude.degrees
+# => 27.500082409271247
 
 horizontal_coordinates.altitude.str(:dms)
-# => "+27° 30′ 0.2967″"
+# => "+27° 30′ 0.2966″"
 ```
 
 ### Sunrise and sunset times and azimuths
 
 ```rb
-date = Date.new(2015, 2, 5)
-epoch = Astronoby::Epoch.from_time(date)
+time = Time.new(2015, 2, 5)
 observer = Astronoby::Observer.new(
   latitude: Astronoby::Angle.from_degrees(38),
   longitude: Astronoby::Angle.from_degrees(-78)
 )
-sun = Astronoby::Sun.new(epoch: epoch)
+sun = Astronoby::Sun.new(time: time)
+observation_events = sun.observation_events(observer: observer)
 
-sun.rising_time(observer: observer)
-# => 2015-02-05 12:13:26 UTC
+observation_events.rising_time
+# => 2015-02-05 12:12:59 UTC
 
-sun.rising_azimuth(observer: observer).str(:dms)
-# => "+109° 41′ 22.2585″"
+observation_events.rising_azimuth.str(:dms)
+# => "+109° 46′ 43.145″"
 
-sun.setting_time(observer: observer)
-# => 2015-02-05 22:35:12 UTC
+observation_events.transit_time
+# => 2015-02-05 17:25:59 UTC
 
-sun.setting_azimuth(observer: observer).str(:dms)
-# => "+250° 18′ 37.7414″"
+observation_events.transit_altitude.str(:dms)
+# => "+36° 8′ 15.7638″"
+
+observation_events.setting_time
+# => 2015-02-05 22:39:27 UTC
+
+observation_events.setting_azimuth.str(:dms)
+# => "+250° 23′ 33.614″"
+```
+
+### Twilight times
+
+```rb
+time = Time.new(2024, 1, 1)
+sun = Astronoby::Sun.new(time: time)
+observer = Astronoby::Observer.new(
+  latitude: Astronoby::Angle.from_degrees(48.8566),
+  longitude: Astronoby::Angle.from_degrees(2.3522)
+)
+twilight_events = sun.twilight_events(observer: observer)
+
+twilight_events.morning_astronomical_twilight_time
+# => 2024-01-01 05:49:25 UTC
+
+twilight_events.morning_nautical_twilight_time
+# => 2024-01-01 06:27:42 UTC
+
+twilight_events.morning_civil_twilight_time
+# => 2024-01-01 07:07:50 UTC
+
+twilight_events.evening_civil_twilight_time
+# => 2024-01-01 16:40:01 UTC
+
+twilight_events.evening_nautical_twilight_time
+# => 2024-01-01 17:20:10 UTC
+
+twilight_events.evening_astronomical_twilight_time
+# => 2024-01-01 17:58:26 UTC
 ```
 
 ### Solstice and Equinox times

data/UPGRADING.md

@@ -7,6 +7,30 @@ changes to it as long as a major version has not been released.
 If you are already using Astronoby and wish to follow the changes to its
 public API, please read the upgrading notes for each release.
 
+## Upgrading from 0.3.0 to 0.4.0
+
+### `Body` class removed ([#50])
+
+The different behaviors from `Body` have been moved to other classes like
+`Events::ObservationEvents`.
+
+[#50]: https://github.com/rhannequin/astronoby/pull/50
+
+### Rising and setting times and azimuths removed from `Sun` ([#60])
+
+`#rising_time`, `#rising_azimuth`, `#setting_time` and `#setting_azimuth`
+have been removed from `Astronoby::Sun` and moved to
+`Astronoby::Events::ObservationEvents`.
+
+[#60]: https://github.com/rhannequin/astronoby/pull/60
+
+### `Sun` constructor changed ([#64])
+
+The `Sun` constructor now doesn't accept the `epoch` key argument anymore,
+but only a new `time` key argument.
+
+[#64]: https://github.com/rhannequin/astronoby/pull/64
+
 ## Upgrading from 0.2.0 to 0.3.0
 
 ### `Sun#ecliptic_coordinates` method removed (#41)

data/lib/astronoby/aberration.rb

@@ -2,6 +2,8 @@
 
 module Astronoby
   class Aberration
+    MAXIMUM_SHIFT = Angle.from_degrees(20.5)
+
     def self.for_ecliptic_coordinates(coordinates:, epoch:)
       new(coordinates, epoch).apply
     end
@@ -18,15 +20,15 @@ module Astronoby
     #  Chapter: 36 - Aberration
     def apply
       delta_longitude = Angle.from_degrees(
-        -20.5 * (
+        -MAXIMUM_SHIFT.degrees * (
           sun_longitude - @coordinates.longitude
-        ).cos / @coordinates.latitude.cos / 3600
+        ).cos / @coordinates.latitude.cos / Constants::SECONDS_PER_DEGREE
       )
 
       delta_latitude = Angle.from_degrees(
-        -20.5 *
+        -MAXIMUM_SHIFT.degrees *
         (sun_longitude - @coordinates.longitude).sin *
-        @coordinates.latitude.sin / 3600
+        @coordinates.latitude.sin / Constants::SECONDS_PER_DEGREE
       )
 
       Coordinates::Ecliptic.new(
@@ -37,7 +39,7 @@ module Astronoby
 
     def sun_longitude
       @_sun_longitude ||= Sun
-        .new(epoch: @epoch)
+        .new(time: Epoch.to_utc(@epoch))
         .true_ecliptic_coordinates
         .longitude
     end

data/lib/astronoby/angle.rb

@@ -1,23 +1,10 @@
 # frozen_string_literal: true
 
-require "bigdecimal/math"
-
 module Astronoby
   class Angle
     include Comparable
 
-    PRECISION = 14
-    PI = BigMath.PI(PRECISION)
-    PI_IN_DEGREES = BigDecimal("180")
-
-    FULL_CIRCLE_IN_RADIANS = (2 * PI)
-
-    RADIAN_PER_HOUR = PI / BigDecimal("12")
-    MINUTES_PER_DEGREE = BigDecimal("60")
-    MINUTES_PER_HOUR = BigDecimal("60")
-    SECONDS_PER_MINUTE = BigDecimal("60")
-    SECONDS_PER_HOUR = MINUTES_PER_HOUR * SECONDS_PER_MINUTE
-
+    MIN_PRECISION = 10
     FORMATS = %i[dms hms].freeze
 
     class << self
@@ -26,28 +13,32 @@ module Astronoby
       end
 
       def from_radians(radians)
-        normalized_radians = radians.remainder(FULL_CIRCLE_IN_RADIANS)
+        normalized_radians = radians.remainder(Constants::RADIANS_PER_CIRCLE)
         new(normalized_radians)
       end
 
       def from_degrees(degrees)
-        radians = degrees / PI_IN_DEGREES * PI
+        radians = degrees / Constants::PI_IN_DEGREES * Math::PI
         from_radians(radians)
       end
 
       def from_hours(hours)
-        radians = hours * RADIAN_PER_HOUR
+        radians = hours * Constants::RADIAN_PER_HOUR
         from_radians(radians)
       end
 
       def from_hms(hour, minute, second)
-        hours = hour + minute / MINUTES_PER_HOUR + second / SECONDS_PER_HOUR
+        hours = hour +
+          minute / Constants::MINUTES_PER_HOUR +
+          second / Constants::SECONDS_PER_HOUR
         from_hours(hours)
       end
 
       def from_dms(degree, minute, second)
         sign = degree.negative? ? -1 : 1
-        degrees = degree.abs + minute / MINUTES_PER_HOUR + second / SECONDS_PER_HOUR
+        degrees = degree.abs +
+          minute / Constants::MINUTES_PER_HOUR +
+          second / Constants::SECONDS_PER_HOUR
         from_degrees(sign * degrees)
       end
 
@@ -70,20 +61,16 @@ module Astronoby
     attr_reader :radians
 
     def initialize(radians)
-      @radians = if radians.is_a?(Integer) || radians.is_a?(BigDecimal)
-        BigDecimal(radians)
-      else
-        BigDecimal(radians, PRECISION)
-      end
+      @radians = radians
       freeze
     end
 
     def degrees
-      @radians * PI_IN_DEGREES / PI
+      @radians * Constants::PI_IN_DEGREES / Math::PI
     end
 
     def hours
-      @radians / RADIAN_PER_HOUR
+      @radians / Constants::RADIAN_PER_HOUR
     end
 
     def +(other)
@@ -94,6 +81,10 @@ module Astronoby
       self.class.from_radians(@radians - other.radians)
     end
 
+    def -@
+      self.class.from_radians(-@radians)
+    end
+
     def sin
       Math.sin(radians)
     end
@@ -144,31 +135,32 @@ module Astronoby
       sign = deg.negative? ? "-" : "+"
       absolute_degrees = deg.abs
       degrees = absolute_degrees.floor
-      decimal_minutes = MINUTES_PER_DEGREE * (absolute_degrees - degrees)
+      decimal_minutes = Constants::MINUTES_PER_DEGREE *
+        (absolute_degrees - degrees)
       absolute_decimal_minutes = (
-        MINUTES_PER_DEGREE * (absolute_degrees - degrees)
+        Constants::MINUTES_PER_DEGREE * (absolute_degrees - degrees)
       ).abs
       minutes = decimal_minutes.floor
-      seconds = SECONDS_PER_MINUTE * (
+      seconds = Constants::SECONDS_PER_MINUTE * (
         absolute_decimal_minutes - absolute_decimal_minutes.floor
       )
 
-      Dms.new(sign, degrees, minutes, seconds.to_f.floor(4))
+      Dms.new(sign, degrees, minutes, seconds.floor(4))
     end
 
     def to_hms(hrs)
       absolute_hours = hrs.abs
       hours = absolute_hours.floor
-      decimal_minutes = MINUTES_PER_HOUR * (absolute_hours - hours)
+      decimal_minutes = Constants::MINUTES_PER_HOUR * (absolute_hours - hours)
       absolute_decimal_minutes = (
-        MINUTES_PER_HOUR * (absolute_hours - hours)
+        Constants::MINUTES_PER_HOUR * (absolute_hours - hours)
       ).abs
       minutes = decimal_minutes.floor
-      seconds = SECONDS_PER_MINUTE * (
+      seconds = Constants::SECONDS_PER_MINUTE * (
         absolute_decimal_minutes - absolute_decimal_minutes.floor
       )
 
-      Hms.new(hours, minutes, seconds.to_f.floor(4))
+      Hms.new(hours, minutes, seconds.floor(4))
     end
   end
 end