I found it hard to compare the results on this library with different sources like the IMCCE, NASA JPL, Heavens Above, etc.

The main reason is the difference of meaning of coordinates. Regarding the Sun, for example, books like Practical Astronomy or Celestial Calculations provide algorithms to calculate coordinates to the date (current epoch). There are in fact different definitions of coordinates:

• Mean coordinates for standard epoch J2000
• Apparent coordinates for standard epoch J2000
• Mean coordinates for the date's epoch
• Apparent coordinates for the date's epoch

While we already have a way to convert equatorial coordinates from one epoch to another, the current algorithms I have access to will provide only Sun ecliptic coordinates for the date's epoch.

I can still manage to make the difference between true and apparent coordinates, which is the purpose of this change.

This change removes Astronoby::Sun#ecliptic_coordinates and Astronoby::Coordinates::Ecliptic#to_equatorial in favour of the following:

• Astronoby::Sun#true_ecliptic_coordinates: true coordinates
• Astronoby::Sun#apparent_ecliptic_coordinates: include corrections for nutation and aberration
• Astronoby::Coordinates::Ecliptic#to_true_equatorial: true equatorial coordinates for the date
• Astronoby::Coordinates::Ecliptic#to_apparent_equatorial: include corrections for true obliquity

time = Time.utc(1988, 7, 27, 0, 0, 0)
epoch = Astronoby::Epoch.from_time(time)

true_ecliptic_coordinates = Astronoby::Sun.new(epoch: epoch).true_ecliptic_coordinates
true_equatorial_coordinates = ecliptic_coordinates.to_true_equatorial(epoch: epoch)

true_equatorial_coordinates.right_ascension.str(:hms)
# => "8h 26m 3.6131s"

true_equatorial_coordinates.declination.str(:dms)
=> "+19° 12′ 43.1836″"

apparent_ecliptic_coordinates = Astronoby::Sun.new(epoch: epoch).apparent_ecliptic_coordinates
apparent_equatorial_coordinates = ecliptic_coordinates.to_apparent_equatorial(epoch: epoch)

apparent_equatorial_coordinates.right_ascension.str(:hms)
# => "8h 26m 3.7342s"

apparent_equatorial_coordinates.declination.str(:dms)
# => "+19° 12′ 50.2238″"

As you can see the difference is very subtle, but my goal is to provide accurate data and more importantly accurate concepts. I know this library still suffers from some approximations, but I aim to correct this problems over time.