Kosmorro
/
lib
miroir de https://github.com/Kosmorro/lib
1
0
Bifurcation 0
Code Tickets Versions 18 Activité
A library that computes the ephemerides.
Vous ne pouvez pas sélectionner plus de 25 sujets Les noms de sujets doivent commencer par une lettre ou un nombre, peuvent contenir des tirets ('-') et peuvent comporter jusqu'à 35 caractères.
149 Révisions
7 Branches
1.2 MiB
 
 
Aborescence: 2772edba16
Branches Tags
${ item.name }
Créer la branche ${ searchTerm }
de '2772edba16'
${ noResults }
lib/kosmorrolib/ephemerides.py

358 lignes
21 KiB
Brut Annotations Historique 

  1. #!/usr/bin/env python3

  2. 

  3. #    Kosmorrolib - The Library To Compute Your Ephemerides

  4. #    Copyright (C) 2021  Jérôme Deuchnord <jerome@deuchnord.fr>

  5. #

  6. #    This program is free software: you can redistribute it and/or modify

  7. #    it under the terms of the GNU Affero General Public License as

  8. #    published by the Free Software Foundation, either version 3 of the

  9. #    License, or (at your option) any later version.

  10. #

  11. #    This program is distributed in the hope that it will be useful,

  12. #    but WITHOUT ANY WARRANTY; without even the implied warranty of

  13. #    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

  14. #    GNU Affero General Public License for more details.

  15. #

  16. #    You should have received a copy of the GNU Affero General Public License

  17. #    along with this program.  If not, see <https://www.gnu.org/licenses/>.

  18. 

  19. from datetime import date, datetime, timedelta

  20. from typing import Union

  21. 

  22. from skyfield.searchlib import find_discrete, find_maxima

  23. from skyfield.timelib import Time

  24. from skyfield.constants import tau

  25. from skyfield.errors import EphemerisRangeError

  26. 

  27. from .model import Position, AsterEphemerides, MoonPhase, Object, ASTERS

  28. from .dateutil import translate_to_utc_offset, normalize_datetime

  29. from .core import (

  30.     get_skf_objects,

  31.     get_timescale,

  32.     get_iau2000b,

  33.     deprecated,

  34.     alert_deprecation,

  35. )

  36. from .enum import MoonPhaseType

  37. from .exceptions import OutOfRangeDateError

  38. 

  39. RISEN_ANGLE = -0.8333

  40. 

  41. 

  42. def _get_skyfield_to_moon_phase(

  43.     times: [Time], vals: [int], now: Time, utc_offset: Union[int, float]

  44. ) -> Union[MoonPhase, None]:

  45.     tomorrow = get_timescale().utc(

  46.         now.utc_datetime().year, now.utc_datetime().month, now.utc_datetime().day + 1

  47.     )

  48. 

  49.     next_phase_time = None

  50.     i = 0

  51. 

  52.     # Find the current moon phase:

  53.     for i, time in enumerate(times):

  54.         if now.utc_datetime() <= time.utc_datetime():

  55.             if time.utc_datetime() >= tomorrow.utc_datetime():

  56.                 i -= 1

  57.             break

  58. 

  59.     current_phase = MoonPhaseType(vals[i])

  60. 

  61.     if current_phase in [

  62.         MoonPhaseType.NEW_MOON,

  63.         MoonPhaseType.FIRST_QUARTER,

  64.         MoonPhaseType.FULL_MOON,

  65.         MoonPhaseType.LAST_QUARTER,

  66.     ]:

  67.         current_phase_time = translate_to_utc_offset(

  68.             times[i].utc_datetime(), utc_offset

  69.         )

  70.     else:

  71.         current_phase_time = None

  72. 

  73.     # Find the next moon phase

  74.     for j in range(i + 1, len(times)):

  75.         if vals[j] in [0, 2, 4, 6]:

  76.             next_phase_time = translate_to_utc_offset(

  77.                 times[j].utc_datetime(), utc_offset

  78.             )

  79.             break

  80. 

  81.     return MoonPhase(current_phase, current_phase_time, next_phase_time)

  82. 

  83. 

  84. def get_moon_phase(

  85.     for_date: date = date.today(), utc_offset: Union[int, float] = 0, **argv

  86. ) -> MoonPhase:

  87.     """Calculate and return the moon phase for the given date, adjusted to the given UTC offset if any.

  88. 

  89.     Get the moon phase for the 27 March, 2021:

  90. 

  91.     >>> get_moon_phase(date(2021, 3, 27))

  92.     <MoonPhase phase_type=MoonPhaseType.WAXING_GIBBOUS time=None next_phase_date=2021-03-28 18:48:10.902298+00:00>

  93. 

  94.     When the moon phase is a new moon, a first quarter, a full moon or a last quarter, you get the exact time

  95.     of its happening too:

  96. 

  97.     >>> get_moon_phase(datetime(2021, 3, 28))

  98.     <MoonPhase phase_type=MoonPhaseType.FULL_MOON time=2021-03-28 18:48:10.902298+00:00 next_phase_date=2021-04-04 10:02:27.393689+00:00>

  99. 

  100.     Get the moon phase for the 27 March, 2021, in UTC+2:

  101. 

  102.     >>> get_moon_phase(date(2021, 3, 27), utc_offset=2)

  103.     <MoonPhase phase_type=MoonPhaseType.WAXING_GIBBOUS time=None next_phase_date=2021-03-28 20:48:10.902298+02:00>

  104. 

  105.     Note that the `utc_offset` argument was named `timezone` before version 1.1. The old name still works, but will be dropped in the future.

  106. 

  107.     >>> get_moon_phase(date(2021, 3, 27), timezone=2)

  108.     <MoonPhase phase_type=MoonPhaseType.WAXING_GIBBOUS time=None next_phase_date=2021-03-28 20:48:10.902298+02:00>

  109. 

  110.     Note that the moon phase can only be computed for a date range.

  111.     Asking for the moon phase with an out of range date will result in an exception:

  112. 

  113.     >>> get_moon_phase(date(1000, 1, 1))

  114.     Traceback (most recent call last):

  115.         ...

  116.     kosmorrolib.exceptions.OutOfRangeDateError: The date must be between 1899-08-09 and 2053-09-26

  117.     """

  118. 

  119.     if argv.get("timezone") is not None:

  120.         alert_deprecation(

  121.             "'timezone' argument of the get_moon_phase() function is deprecated. Use utc_offset instead."

  122.         )

  123.         utc_offset = argv.get("timezone")

  124. 

  125.     earth = get_skf_objects()["earth"]

  126.     moon = get_skf_objects()["moon"]

  127.     sun = get_skf_objects()["sun"]

  128. 

  129.     def moon_phase_at(time: Time):

  130.         time._nutation_angles = get_iau2000b(time)

  131.         current_earth = earth.at(time)

  132.         _, mlon, _ = current_earth.observe(moon).apparent().ecliptic_latlon("date")

  133.         _, slon, _ = current_earth.observe(sun).apparent().ecliptic_latlon("date")

  134.         return (((mlon.radians - slon.radians) // (tau / 8)) % 8).astype(int)

  135. 

  136.     moon_phase_at.rough_period = 7.0  # one lunar phase per week

  137. 

  138.     today = get_timescale().utc(for_date.year, for_date.month, for_date.day)

  139.     start_time = get_timescale().utc(for_date.year, for_date.month, for_date.day - 10)

  140.     end_time = get_timescale().utc(for_date.year, for_date.month, for_date.day + 10)

  141. 

  142.     try:

  143.         times, phases = find_discrete(start_time, end_time, moon_phase_at)

  144.         return _get_skyfield_to_moon_phase(times, phases, today, utc_offset)

  145. 

  146.     except EphemerisRangeError as error:

  147.         start = translate_to_utc_offset(error.start_time.utc_datetime(), utc_offset)

  148.         end = translate_to_utc_offset(error.end_time.utc_datetime(), utc_offset)

  149. 

  150.         start = date(start.year, start.month, start.day) + timedelta(days=12)

  151.         end = date(end.year, end.month, end.day) - timedelta(days=12)

  152. 

  153.         raise OutOfRangeDateError(start, end) from error

  154. 

  155. 

  156. def get_ephemerides(

  157.     position: Position,

  158.     for_date: date = date.today(),

  159.     utc_offset: Union[int, float] = 0,

  160.     **argv

  161. ) -> [AsterEphemerides]:

  162.     """Compute and return the ephemerides for the given position and date, adjusted to the given UTC offset if any.

  163. 

  164.     Compute the ephemerides for July 7th, 2022:

  165. 

  166.     >>> get_ephemerides(Position(36.6794, 4.8555), date(2022, 7, 7))

  167.     [<AsterEphemerides rise_time=2022-07-07 04:29:00 culmination_time=2022-07-07 11:46:00 set_time=2022-07-07 19:02:00 aster=<Object type=STAR name=SUN />>,

  168.     <AsterEphemerides rise_time=2022-07-07 12:16:00 culmination_time=2022-07-07 18:06:00 set_time=2022-07-07 23:54:00 aster=<Object type=SATELLITE name=MOON />>,

  169.     <AsterEphemerides rise_time=2022-07-07 03:36:00 culmination_time=2022-07-07 10:58:00 set_time=2022-07-07 18:20:00 aster=<Object type=PLANET name=MERCURY />>,

  170.     <AsterEphemerides rise_time=2022-07-07 02:30:00 culmination_time=2022-07-07 09:44:00 set_time=2022-07-07 16:58:00 aster=<Object type=PLANET name=VENUS />>,

  171.     <AsterEphemerides rise_time=2022-07-07 00:05:00 culmination_time=2022-07-07 06:39:00 set_time=2022-07-07 13:14:00 aster=<Object type=PLANET name=MARS />>,

  172.     <AsterEphemerides rise_time=2022-07-07 22:59:00 culmination_time=2022-07-07 05:11:00 set_time=2022-07-07 11:20:00 aster=<Object type=PLANET name=JUPITER />>,

  173.     <AsterEphemerides rise_time=2022-07-07 21:06:00 culmination_time=2022-07-07 02:29:00 set_time=2022-07-07 07:48:00 aster=<Object type=PLANET name=SATURN />>,

  174.     <AsterEphemerides rise_time=2022-07-07 00:47:00 culmination_time=2022-07-07 07:42:00 set_time=2022-07-07 14:38:00 aster=<Object type=PLANET name=URANUS />>,

  175.     <AsterEphemerides rise_time=2022-07-07 22:27:00 culmination_time=2022-07-07 04:25:00 set_time=2022-07-07 10:20:00 aster=<Object type=PLANET name=NEPTUNE />>,

  176.     <AsterEphemerides rise_time=2022-07-07 19:46:00 culmination_time=2022-07-07 00:41:00 set_time=2022-07-07 05:33:00 aster=<Object type=PLANET name=PLUTO />>]

  177. 

  178.     UTC offset can be optionnaly set to adapt the hours to your location:

  179. 

  180.     >>> get_ephemerides(Position(36.6794, 4.8555), date(2022, 7, 7), utc_offset=2)

  181.     [<AsterEphemerides rise_time=2022-07-07 06:29:00 culmination_time=2022-07-07 13:46:00 set_time=2022-07-07 21:02:00 aster=<Object type=STAR name=SUN />>,

  182.     <AsterEphemerides rise_time=2022-07-07 14:16:00 culmination_time=2022-07-07 20:06:00 set_time=2022-07-07 01:27:00 aster=<Object type=SATELLITE name=MOON />>,

  183.     <AsterEphemerides rise_time=2022-07-07 05:36:00 culmination_time=2022-07-07 12:58:00 set_time=2022-07-07 20:20:00 aster=<Object type=PLANET name=MERCURY />>,

  184.     <AsterEphemerides rise_time=2022-07-07 04:30:00 culmination_time=2022-07-07 11:44:00 set_time=2022-07-07 18:58:00 aster=<Object type=PLANET name=VENUS />>,

  185.     <AsterEphemerides rise_time=2022-07-07 02:05:00 culmination_time=2022-07-07 08:39:00 set_time=2022-07-07 15:14:00 aster=<Object type=PLANET name=MARS />>,

  186.     <AsterEphemerides rise_time=2022-07-07 01:02:00 culmination_time=2022-07-07 07:11:00 set_time=2022-07-07 13:20:00 aster=<Object type=PLANET name=JUPITER />>,

  187.     <AsterEphemerides rise_time=2022-07-07 23:06:00 culmination_time=2022-07-07 04:29:00 set_time=2022-07-07 09:48:00 aster=<Object type=PLANET name=SATURN />>,

  188.     <AsterEphemerides rise_time=2022-07-07 02:47:00 culmination_time=2022-07-07 09:42:00 set_time=2022-07-07 16:38:00 aster=<Object type=PLANET name=URANUS />>,

  189.     <AsterEphemerides rise_time=2022-07-07 00:31:00 culmination_time=2022-07-07 06:25:00 set_time=2022-07-07 12:20:00 aster=<Object type=PLANET name=NEPTUNE />>,

  190.     <AsterEphemerides rise_time=2022-07-07 21:46:00 culmination_time=2022-07-07 02:41:00 set_time=2022-07-07 07:33:00 aster=<Object type=PLANET name=PLUTO />>]

  191. 

  192. 

  193.     Note that the `utc_offset` argument was named `timezone` before version 1.1. The old name still works, but will be dropped in the future.

  194. 

  195.     >>> get_ephemerides(Position(36.6794, 4.8555), date(2022, 7, 7), timezone=2)

  196.     [<AsterEphemerides rise_time=2022-07-07 06:29:00 culmination_time=2022-07-07 13:46:00 set_time=2022-07-07 21:02:00 aster=<Object type=STAR name=SUN />>,

  197.     <AsterEphemerides rise_time=2022-07-07 14:16:00 culmination_time=2022-07-07 20:06:00 set_time=2022-07-07 01:27:00 aster=<Object type=SATELLITE name=MOON />>,

  198.     <AsterEphemerides rise_time=2022-07-07 05:36:00 culmination_time=2022-07-07 12:58:00 set_time=2022-07-07 20:20:00 aster=<Object type=PLANET name=MERCURY />>,

  199.     <AsterEphemerides rise_time=2022-07-07 04:30:00 culmination_time=2022-07-07 11:44:00 set_time=2022-07-07 18:58:00 aster=<Object type=PLANET name=VENUS />>,

  200.     <AsterEphemerides rise_time=2022-07-07 02:05:00 culmination_time=2022-07-07 08:39:00 set_time=2022-07-07 15:14:00 aster=<Object type=PLANET name=MARS />>,

  201.     <AsterEphemerides rise_time=2022-07-07 01:02:00 culmination_time=2022-07-07 07:11:00 set_time=2022-07-07 13:20:00 aster=<Object type=PLANET name=JUPITER />>,

  202.     <AsterEphemerides rise_time=2022-07-07 23:06:00 culmination_time=2022-07-07 04:29:00 set_time=2022-07-07 09:48:00 aster=<Object type=PLANET name=SATURN />>,

  203.     <AsterEphemerides rise_time=2022-07-07 02:47:00 culmination_time=2022-07-07 09:42:00 set_time=2022-07-07 16:38:00 aster=<Object type=PLANET name=URANUS />>,

  204.     <AsterEphemerides rise_time=2022-07-07 00:31:00 culmination_time=2022-07-07 06:25:00 set_time=2022-07-07 12:20:00 aster=<Object type=PLANET name=NEPTUNE />>,

  205.     <AsterEphemerides rise_time=2022-07-07 21:46:00 culmination_time=2022-07-07 02:41:00 set_time=2022-07-07 07:33:00 aster=<Object type=PLANET name=PLUTO />>]

  206. 

  207.     Objects may not rise or set on the given date (e.g. they rise the previous day or set the next day).

  208.     In this case, you will get `None` values on the rise or set time.

  209. 

  210.     If an objet does not rise nor set due to your latitude, then both rise and set will be `None`:

  211. 

  212.     >>> north_pole = Position(70, 20)

  213.     >>> south_pole = Position(-70, 20)

  214.     >>> get_ephemerides(north_pole, date(2021, 6, 20))

  215.     [<AsterEphemerides rise_time=None culmination_time=2021-06-20 10:42:00 set_time=None aster=<Object type=STAR name=SUN />>,

  216.     <AsterEphemerides rise_time=2021-06-20 14:30:00 culmination_time=2021-06-20 18:44:00 set_time=2021-06-20 22:53:00 aster=<Object type=SATELLITE name=MOON />>,

  217.     <AsterEphemerides rise_time=2021-06-20 22:56:00 culmination_time=2021-06-20 09:47:00 set_time=2021-06-20 20:34:00 aster=<Object type=PLANET name=MERCURY />>,

  218.     <AsterEphemerides rise_time=None culmination_time=2021-06-20 12:20:00 set_time=None aster=<Object type=PLANET name=VENUS />>,

  219.     <AsterEphemerides rise_time=None culmination_time=2021-06-20 13:17:00 set_time=None aster=<Object type=PLANET name=MARS />>,

  220.     <AsterEphemerides rise_time=2021-06-20 23:06:00 culmination_time=2021-06-20 03:04:00 set_time=2021-06-20 06:58:00 aster=<Object type=PLANET name=JUPITER />>,

  221.     <AsterEphemerides rise_time=2021-06-20 23:28:00 culmination_time=2021-06-20 01:48:00 set_time=2021-06-20 04:05:00 aster=<Object type=PLANET name=SATURN />>,

  222.     <AsterEphemerides rise_time=2021-06-20 21:53:00 culmination_time=2021-06-20 07:29:00 set_time=2021-06-20 17:02:00 aster=<Object type=PLANET name=URANUS />>,

  223.     <AsterEphemerides rise_time=2021-06-20 22:51:00 culmination_time=2021-06-20 04:22:00 set_time=2021-06-20 09:50:00 aster=<Object type=PLANET name=NEPTUNE />>,

  224.     <AsterEphemerides rise_time=None culmination_time=2021-06-20 00:40:00 set_time=None aster=<Object type=PLANET name=PLUTO />>]

  225. 

  226.     >>> get_ephemerides(north_pole, date(2021, 12, 21))

  227.     [<AsterEphemerides rise_time=None culmination_time=2021-12-21 10:38:00 set_time=None aster=<Object type=STAR name=SUN />>,

  228.     <AsterEphemerides rise_time=None culmination_time=2021-12-21 00:04:00 set_time=None aster=<Object type=SATELLITE name=MOON />>,

  229.     <AsterEphemerides rise_time=None culmination_time=2021-12-21 11:33:00 set_time=None aster=<Object type=PLANET name=MERCURY />>,

  230.     <AsterEphemerides rise_time=2021-12-21 11:58:00 culmination_time=2021-12-21 12:33:00 set_time=2021-12-21 13:08:00 aster=<Object type=PLANET name=VENUS />>,

  231.     <AsterEphemerides rise_time=None culmination_time=2021-12-21 08:54:00 set_time=None aster=<Object type=PLANET name=MARS />>,

  232.     <AsterEphemerides rise_time=2021-12-21 11:07:00 culmination_time=2021-12-21 14:43:00 set_time=2021-12-21 18:19:00 aster=<Object type=PLANET name=JUPITER />>,

  233.     <AsterEphemerides rise_time=2021-12-21 11:32:00 culmination_time=2021-12-21 13:33:00 set_time=2021-12-21 15:33:00 aster=<Object type=PLANET name=SATURN />>,

  234.     <AsterEphemerides rise_time=2021-12-21 09:54:00 culmination_time=2021-12-21 19:13:00 set_time=2021-12-21 04:37:00 aster=<Object type=PLANET name=URANUS />>,

  235.     <AsterEphemerides rise_time=2021-12-21 10:49:00 culmination_time=2021-12-21 16:05:00 set_time=2021-12-21 21:21:00 aster=<Object type=PLANET name=NEPTUNE />>,

  236.     <AsterEphemerides rise_time=None culmination_time=2021-12-21 12:31:00 set_time=None aster=<Object type=PLANET name=PLUTO />>]

  237. 

  238.     >>> get_ephemerides(south_pole, date(2021, 6, 20))

  239.     [<AsterEphemerides rise_time=None culmination_time=2021-06-20 10:42:00 set_time=None aster=<Object type=STAR name=SUN />>,

  240.     <AsterEphemerides rise_time=2021-06-20 11:10:00 culmination_time=2021-06-20 19:06:00 set_time=2021-06-20 01:20:00 aster=<Object type=SATELLITE name=MOON />>,

  241.     <AsterEphemerides rise_time=2021-06-20 07:47:00 culmination_time=2021-06-20 09:47:00 set_time=2021-06-20 11:48:00 aster=<Object type=PLANET name=MERCURY />>,

  242.     <AsterEphemerides rise_time=None culmination_time=2021-06-20 12:20:00 set_time=None aster=<Object type=PLANET name=VENUS />>,

  243.     <AsterEphemerides rise_time=2021-06-20 12:14:00 culmination_time=2021-06-20 13:17:00 set_time=2021-06-20 14:21:00 aster=<Object type=PLANET name=MARS />>,

  244.     <AsterEphemerides rise_time=2021-06-20 18:32:00 culmination_time=2021-06-20 03:04:00 set_time=2021-06-20 11:32:00 aster=<Object type=PLANET name=JUPITER />>,

  245.     <AsterEphemerides rise_time=2021-06-20 15:20:00 culmination_time=2021-06-20 01:48:00 set_time=2021-06-20 12:12:00 aster=<Object type=PLANET name=SATURN />>,

  246.     <AsterEphemerides rise_time=2021-06-20 04:32:00 culmination_time=2021-06-20 07:29:00 set_time=2021-06-20 10:26:00 aster=<Object type=PLANET name=URANUS />>,

  247.     <AsterEphemerides rise_time=2021-06-20 21:28:00 culmination_time=2021-06-20 04:22:00 set_time=2021-06-20 11:13:00 aster=<Object type=PLANET name=NEPTUNE />>,

  248.     <AsterEphemerides rise_time=None culmination_time=2021-06-20 00:40:00 set_time=None aster=<Object type=PLANET name=PLUTO />>]

  249. 

  250.     >>> get_ephemerides(south_pole, date(2021, 12, 22))

  251.     [<AsterEphemerides rise_time=None culmination_time=2021-12-22 10:39:00 set_time=None aster=<Object type=STAR name=SUN />>,

  252.     <AsterEphemerides rise_time=None culmination_time=2021-12-22 01:01:00 set_time=None aster=<Object type=SATELLITE name=MOON />>,

  253.     <AsterEphemerides rise_time=None culmination_time=2021-12-22 11:35:00 set_time=None aster=<Object type=PLANET name=MERCURY />>,

  254.     <AsterEphemerides rise_time=None culmination_time=2021-12-22 12:27:00 set_time=None aster=<Object type=PLANET name=VENUS />>,

  255.     <AsterEphemerides rise_time=None culmination_time=2021-12-22 08:53:00 set_time=None aster=<Object type=PLANET name=MARS />>,

  256.     <AsterEphemerides rise_time=2021-12-22 05:52:00 culmination_time=2021-12-22 14:40:00 set_time=2021-12-22 23:26:00 aster=<Object type=PLANET name=JUPITER />>,

  257.     <AsterEphemerides rise_time=2021-12-22 02:41:00 culmination_time=2021-12-22 13:29:00 set_time=2021-12-22 00:21:00 aster=<Object type=PLANET name=SATURN />>,

  258.     <AsterEphemerides rise_time=2021-12-22 16:01:00 culmination_time=2021-12-22 19:09:00 set_time=2021-12-22 22:17:00 aster=<Object type=PLANET name=URANUS />>,

  259.     <AsterEphemerides rise_time=2021-12-22 08:59:00 culmination_time=2021-12-22 16:01:00 set_time=2021-12-22 23:04:00 aster=<Object type=PLANET name=NEPTUNE />>,

  260.     <AsterEphemerides rise_time=None culmination_time=2021-12-22 12:27:00 set_time=None aster=<Object type=PLANET name=PLUTO />>]

  261. 

  262.     Please note:

  263.     - The ephemerides can only be computed for a date range. Asking for the ephemerides with an out of range date will result in an exception:

  264. 

  265.       >>> get_ephemerides(Position(50.5824, 3.0624), date(1000, 1, 1))

  266.       Traceback (most recent call last):

  267.           ...

  268.       kosmorrolib.exceptions.OutOfRangeDateError: The date must be between 1899-07-29 and 2053-10-07

  269. 

  270.     - The date given in parameter is considered as being given from the point of view of the given UTC offset.

  271.       Using a UTC offset that does not correspond to the place's actual one can impact the returned times.

  272.     """

  273. 

  274.     if argv.get("timezone") is not None:

  275.         alert_deprecation(

  276.             "'timezone' argument of the get_ephemerides() function is deprecated. Use utc_offset instead."

  277.         )

  278.         utc_offset = argv.get("timezone")

  279. 

  280.     def get_angle(for_aster: Object):

  281.         def fun(time: Time) -> float:

  282.             return (

  283.                 position.get_planet_topos()

  284.                 .at(time)

  285.                 .observe(for_aster.skyfield_object)

  286.                 .apparent()

  287.                 .altaz()[0]

  288.                 .degrees

  289.             )

  290. 

  291.         fun.rough_period = 1.0

  292.         return fun

  293. 

  294.     def is_risen(for_aster: Object):

  295.         def fun(time: Time) -> bool:

  296.             return get_angle(for_aster)(time) > RISEN_ANGLE

  297. 

  298.         fun.rough_period = 0.5

  299.         return fun

  300. 

  301.     # The date given in argument is supposed to be given in the given UTC offset (more natural for a human),

  302.     # but we need it in UTC. Subtracting the offset to get it in UTC.

  303. 

  304.     start_time = get_timescale().utc(

  305.         for_date.year, for_date.month, for_date.day, -utc_offset

  306.     )

  307. 

  308.     end_time = get_timescale().utc(

  309.         for_date.year, for_date.month, for_date.day + 1, -utc_offset

  310.     )

  311. 

  312.     ephemerides = []

  313. 

  314.     try:

  315.         for aster in ASTERS:

  316.             times, risen_info = find_discrete(start_time, end_time, is_risen(aster))

  317.             culmination_time, _ = find_maxima(start_time, end_time, get_angle(aster))

  318. 

  319.             rise_time, set_time = None, None

  320.             culmination_time = (

  321.                 culmination_time[0] if len(culmination_time) == 1 else None

  322.             )

  323. 

  324.             for i, time in enumerate(times):

  325.                 time_dt = normalize_datetime(

  326.                     translate_to_utc_offset(time.utc_datetime(), to_tz=utc_offset)

  327.                 )

  328. 

  329.                 if time_dt is not None and time_dt.day != for_date.day:

  330.                     continue

  331. 

  332.                 if risen_info[i]:

  333.                     rise_time = time_dt

  334.                 else:

  335.                     set_time = time_dt

  336. 

  337.             if culmination_time is not None:

  338.                 culmination_time = normalize_datetime(

  339.                     translate_to_utc_offset(

  340.                         culmination_time.utc_datetime(),

  341.                         to_tz=utc_offset,

  342.                     )

  343.                 )

  344. 

  345.             ephemerides.append(

  346.                 AsterEphemerides(rise_time, culmination_time, set_time, aster=aster)

  347.             )

  348.     except EphemerisRangeError as error:

  349.         start = translate_to_utc_offset(error.start_time.utc_datetime(), utc_offset)

  350.         end = translate_to_utc_offset(error.end_time.utc_datetime(), utc_offset)

  351. 

  352.         start = date(start.year, start.month, start.day + 1)

  353.         end = date(end.year, end.month, end.day - 1)

  354. 

  355.         raise OutOfRangeDateError(start, end) from error

  356. 

  357.     return ephemerides