R1i%dZddlZddlmZddlmZddlmZddlm Z ddl m Z m Z dd l mZdd lmZdd lmZmZmZed e DZGd deZGddeZGddeZGddeZdZdZdZy)z2An interface between JPL ephemerides and Skyfield.N) defaultdict)OutOfRangeError)SPK)target_name_pairs)AU_KMDAY_S)EphemerisRangeError)compute_calendar_date)VectorFunction VectorSum_jpl_code_name_dictc#*K|] \}}||f ywN).0targetnames S/home/cursorai/projects/iching/venv/lib/python3.12/site-packages/skyfield/jpllib.py rs):~T6N):scFeZdZdZdZdZdZdZdZdZ dZ d Z d Z y ) SpiceKernelaEphemeris file in NASA .bsp format. A "Spacecraft and Planet Kernel" (SPK) file from NASA provides |xyz| coordinates for bodies in the Solar System like the Sun, planets, moons, and spacecraft. You can download a .bsp file yourself and use this class to open it, or use the Skyfield ``load()`` function to automatically download a popular ephemeris. Once loaded, you can print this object to the screen to see a report on the segments that it includes: >>> planets = load('de421.bsp') >>> print(planets) SPICE kernel file 'de421.bsp' has 15 segments JD 2414864.50 - JD 2471184.50 (1899-07-28 through 2053-10-08) 0 -> 1 SOLAR SYSTEM BARYCENTER -> MERCURY BARYCENTER 0 -> 2 SOLAR SYSTEM BARYCENTER -> VENUS BARYCENTER 0 -> 3 SOLAR SYSTEM BARYCENTER -> EARTH BARYCENTER 0 -> 4 SOLAR SYSTEM BARYCENTER -> MARS BARYCENTER 0 -> 5 SOLAR SYSTEM BARYCENTER -> JUPITER BARYCENTER 0 -> 6 SOLAR SYSTEM BARYCENTER -> SATURN BARYCENTER 0 -> 7 SOLAR SYSTEM BARYCENTER -> URANUS BARYCENTER 0 -> 8 SOLAR SYSTEM BARYCENTER -> NEPTUNE BARYCENTER 0 -> 9 SOLAR SYSTEM BARYCENTER -> PLUTO BARYCENTER 0 -> 10 SOLAR SYSTEM BARYCENTER -> SUN 3 -> 301 EARTH BARYCENTER -> MOON 3 -> 399 EARTH BARYCENTER -> EARTH 1 -> 199 MERCURY BARYCENTER -> MERCURY 2 -> 299 VENUS BARYCENTER -> VENUS 4 -> 499 MARS BARYCENTER -> MARS To retrieve the one or more vectors necessary to compute the position of a body relative to the Solar System barycenter, look up the body by its name or official SPICE identifying integer: >>> planets['earth'] 3 EARTH BARYCENTER 'de421.bsp' segment 3 EARTH BARYCENTER -> 399 EARTH> >>> planets[499] 4 MARS BARYCENTER 'de421.bsp' segment 4 MARS BARYCENTER -> 499 MARS> The result will be a :class:`~skyfield.vectorlib.VectorFunction` instance that you can ask for a position at a given input time. c||_tjj||_t j ||_|j jDcgc]}t||c}|_td|jDjd|jD|_ ycc}w)Nc34K|]}|jywrcenterrss rrz'SpiceKernel.__init__..Is9=a=c34K|]}|jywrrrs rrz'SpiceKernel.__init__..Is@:+8a=r) pathosbasenamefilenameropenspksegments SPICESegmentsetunioncodes)selfr"rs r__init__zSpiceKernel.__init__Ds ((. 88D>8<8I8IJ8I1dA.8IJ 94==99??@:+/==@:: KsB;c`djt|j|jS)Nz <{0} {1!r}>)formattype__name__r"r-s r__repr__zSpiceKernel.__repr__Ls###DJ$7$7CCc  |jj}dj|jt |g}dj}dx}}|D]}||j k7s||j k7ri|j |j }}|tt|}|tt|}|jdj|||||jt|dj|S)Nz(SPICE kernel file {0!r} has {1} segmentsz{0}-{1:02}-{2:02}z, JD {0:.2f} - JD {1:.2f} ({2} through {3}) ) r'r(r0r%lenstart_jdend_jdr intappend_format_segmentjoin) r-r(lines format_datestartendrstartsendss r__str__zSpiceKernel.__str__Os88$$;&H 68)00 A "cQXXoZZs$&;CJ&GH"$9#c($CD K$fUC>@ LL+ ,yyr5c|jj|jdd=|jj y)zClose this ephemeris file.N)r'closer(r,clearr3s rrGzSpiceKernel.close_s.  MM!  r5c6|jjS)aReturn the comments string of this kernel. The resulting string often contains embedded newlines, and is formatted for a human reader. >>> print(planets.comments()) ; de421.bsp LOG FILE ; ; Created 2008-02-12/11:33:34.00. ... LEAPSECONDS_FILE = naif0007.tls SPK_FILE = de421.bsp ... )r'commentsr3s rrJzSpiceKernel.commentshs xx  ""r5ctt}tD](\}}||jvs||j |*t |S)a*Return all target names that are valid with this kernel. >>> pprint(planets.names()) {0: ['SOLAR_SYSTEM_BARYCENTER', 'SSB', 'SOLAR SYSTEM BARYCENTER'], 1: ['MERCURY_BARYCENTER', 'MERCURY BARYCENTER'], 2: ['VENUS_BARYCENTER', 'VENUS BARYCENTER'], 3: ['EARTH_BARYCENTER', 'EMB', ... The result is a dictionary with target code keys and name lists as values. The last name in each list is the one that Skyfield uses when printing information about a body. )rlistrr,r<dict)r-dcoders rnameszSpiceKernel.nameszsC  +JD$tzz!$t$,Awr5cZt|tr|}nA|j}tj |}|t dj |||jvrHdjd|jD}tdj |j|||S)aoTranslate a target name into its integer code. >>> planets.decode('Venus') 299 Raises ``ValueError`` if you supply an unknown name, or ``KeyError`` if the target is missing from this kernel. You can supply an integer code if you already have one and just want to check whether it is present in this kernel. zunknown SPICE target {0!r}z, c32K|]}t|ywr)_format_code_and_name)rcs rrz%SpiceKernel.decode..sM*Q 5a 8*sz@kernel {0!r} is missing {1!r} - the targets it supports are: {2}) isinstancer;upper_jpl_name_code_dictget ValueErrorr0r,r>KeyErrorr%)r-rrOtargetss rdecodezSpiceKernel.decodes dC D::.sNX'W^^W5XsrrN) r\r(rMtuple_centerr8rrr )r-rr( segment_dictchainrs r __getitem__zSpiceKernel.__getitem__sV$==NXNN gfl34 u:?8Odd qr!!//r5ct|tr|}n#tj|j }||j vSr)rUr;rWrXrVr,)r- name_or_coderOs r __contains__zSpiceKernel.__contains__s9 lC (D&**<+=+=+?@Dtzz!!r5N) r2 __module__ __qualname____doc__r.r4rErGrJrPr\rerhrr5rrrs6/`:D #$,4 0"r5rc.eZdZdZdZedZdZy)r)c|jdk(rtjtS|jdk(rtjtSt dj |j)Nz#SPK data type {0} not yet supported) data_typeobject__new__ChebyshevPositionChebyshevPositionVelocityrYr0)cls ephemeris spk_segments rrrzSPICESegment.__new__s^  A %>>"34 4  A %>>";< <> &!6!679 9r5cd||_|j|_|j|_||_yr)rvrrrw)r-rvrws rr.zSPICESegment.__init__s+"!(( !(( &r5cLdj|jjS)Nz {0!r} segment)r0rvr"r3s r vector_namezSPICESegment.vector_names%%dnn&9&9::r5c|j}|j|j|j|jfSr)rwtdb_jdr9r:)r-tsrs r time_rangezSPICESegment.time_ranges2   yy$bii&999r5N)r2rirjrrr.propertyrzr~rr5rr)r)s%9' ;;:r5r)ceZdZdZy)rsc|j} |j|j|j\}}|tz |tz ddfS#t$r}|j |j \}}d|jddz}d|jddz}d|d|} |j} |j}t| ||| |}d|_ |d}~wwxYw)Nz%04d-%02d-%02droz$ephemeris segment only covers dates z through ) rwcompute_and_differentiatewhole tdb_fractionrr~r} tdb_calendarout_of_range_timesr __cause__r) r-tr_positionvelocitye start_timeend_timertextmasks r_atzChebyshevPosition._ats"" !(!B!B") Hh%E!14== #'??144#8 J :#:#:#>AJKQOD''D&&G#D*hgNAAKG s)A CBCCNr2rirjrrr5rrsrss>r5rsceZdZdZy)rtc|jj|j|j}|ddtz |ddt ztz ddfS)Nro)rwcomputerrrr )r-rpvs rrzChebyshevPositionVelocity._atsJ    % %aggq~~ >"1v~r!"v~5tTAAr5Nrrr5rrtrtsBr5rtc#JK||vr||}||j}||vryyw)z {1:<3} {2} -> {3})rrXrrr0)r_cnametnames rr=r=sR  # #GNNI >E  # #GNNI >E , 3 3u 66r5)rkr# collectionsrjplephem.exceptionsr jplephem.spkrjplephem.namesr constantsrr errorsr timelibr vectorlibr r rrMrWrqrr)rsrtrbrSr=rr5rrs8 #/,#'*EE):i"&i"V:>:0> >(B B ( 6r5