R1iXddlmZddlZddlZddlZddlZddlmZddlmZddl m Z ddl Z ddl Z ddl mZddlmZdd lmZdd lmZmZmZmZdd lmZdd lmZdd lmZ ddlmZejBdk\r ejDZ#n ejHZ# ddl%m&Z& ddl(m)Z)ddl*m+Z+ e+ddej`jbjddk(Z3dZ4dZ5dZ6dZ7dZ8dZ9dZ:Gd d!e;Zd$Z?d*d%Z@d+d&ZAd'ZBGd(d)e;ZCy#ddl m ZYxYw#e'$rdZ&YwxYw#ddl)m)Z)ddl,m+Z+YxYw#e-$rdZ.Ye/$rdZ.YwxYw),)print_functionN)fnmatch)get_data)time)iers)Splines)load_bundled_npy)_build_legacy_dataparse_deltat_dataparse_deltat_predsparse_leap_seconds) SpiceKernel)EarthSatellite) Timescale)BytesIO)StringIO)r)create_default_context)urlparse)urlopencafileFTPseudoOutputFilecPt|jjddS)z(Return the last path component of a url./)rpathsplit)urls R/home/cursorai/projects/iching/venv/lib/python3.12/site-packages/skyfield/iokit.py _filename_ofr#:s" C=   # #C ( ,,z&https://hpiers.obspm.fr/iers/bul/bulc/z/ftp://ftp.iers.org/products/eop/rapid/standard/z-https://ssd.jpl.nasa.gov/ftp/eph/planets/bsp/z3https://naif.jpl.nasa.gov/pub/naif/generic_kernels/zBhttps://naif.jpl.nasa.gov/pub/naif/generic_kernels/spk/satellites/ct|dS)Nrb)mode)open)rs r" _open_binaryr)Ds 4  r$ceZdZdZddZdZdZdZddZdZ dZ d Z dd Z dd Z dd ZddZddZedZy )LoaderaXA tool for downloading and opening astronomical data files. A default `Loader` that saves data files to the current working directory can be imported directly from the Skyfield API:: from skyfield.api import load But users can also create a `Loader` of their own, if there is another directory they want data files saved to, or if they want to specify different options. The directory is created automatically if it does not yet exist:: from skyfield.api import Loader load = Loader('~/skyfield-data') The options are: ``verbose`` If set to ``False``, then the loader will not print a progress bar to the screen each time it downloads a file. (If the standard output is not a TTY, then no progress bar is printed anyway.) ``expire`` (This option is no longer supported.) Once a `Loader` is created, it can be called like a function to open, or else to download and open, a file whose name it recognizes:: planets = load('de405.bsp') Each loader also supports an attribute and a few methods. cRtjj||_||_g|_ tj |jtttdztdztdzdtfdtfgd|_t"t$t&d|_dt*fgdt,fgdt,fgdt,fgd |_y#t$rQ}|jtjk7r*tjj|jsYd}~d}~wwxYw) Nzfk/satellites/zpck/zpck/a_old_versions/zjup*.bspz*.bsp)finals2000A.allLeap_Second.datzmoon_080317.tfzmoon_pa_de421_1900-2050.bpcz pck00008.tpc.bsp) deltat.data deltat.predsr.*)r/z.bpcz.tpcz.tf)osr expanduser directoryverboseeventsmakedirsOSErrorerrnoEEXISTisdir_IERS2_IERS _NAIF_KERNELS_NAIF_JPLurlsr r rparsersrr)openers)selfr5r6expirees r"__init__zLoader.__init__is++I6    KK ' &$+.>>+86+A),AAU#$   -.1  +& <()<(),'(   1 ww%,,&rww}}T^^/L sC D&AD!!D&cx|jdk(r|Stjj|j|S)z;Return the path to ``filename`` in this loader's directory..)r5r3rjoinrEfilenames r"path_tozLoader.path_tos, >>S Oww||DNNH55r$ctj|j|j}t |z }|dz S)z@Return how recently ``filename`` was modified, measured in days.g@)r3statrNst_mtimer)rErMmtimesecondss r"days_oldzLoader.days_olds5 X./88&5.  r$c^tjj|j|SN)r3rexistsrNrLs r"_existszLoader._existssww~~dll8455r$cd|vr*|}t|jjdd}nt|j|}|r||z }t|j |}t|j |}||tdj||r:|jd||jttd|}||S|j||||} |:|jd|jt| d 5}||cdddS|jd |j|| S#1swYyxYw) z7Open the given file, downloading it first if necessary.://rrNz5Skyfield does not know how to open a file named {0!r}z%{0} Parsing builtin file with {1}()z skyfield.dataz Parsing with {0}()r&z Opening with {0})rrr _searchrBrCrD ValueErrorformat_log__name__rr_assurer() rErMreloadbackupbuiltinr!parseropenerfrs r"__call__zLoader.__call__s9 H C}))//4R8H$))X.Cxx0x0 N,,2F8,<> >  II> 1(;dD!Qay"! II*FOO <$<  "!s EEcD|j|}tjj|}|j ||r|j d|r|rG|t dj ||j d|t|||j||S)Nz File already existsz-Skyfield does not know where to download {!r}z Downloading {0}rb) rNr3rrWr^r\r]downloadr6)rEr!rMrarbrrWs r"r`zLoader._assures||H%% $  II- .6{ !P"(&"244 II)3 / S$ V < r$cT|jj|j|yrV)r7appendr])rEmessageargss r"r^z Loader._logs >7>>401r$cpt|j|}|r||zStdj|)zReturn the URL Skyfield will try downloading for a given filename. Raises ``ValueError`` if Skyfield doesn't know where to get the file based on its name. z&Skyfield doesn't know the URL of {0!r})r[rBr\r])rErMbases r" build_urlzLoader.build_urls;tyy(+ (? "A &*, ,r$Nci}|j|||5}t|D]*\}}|||jj<|D]}|||< , ddd|S#1swY|SxYw)aLoad and parse a satellite TLE file. DEPRECATED: in a misguided attempt to be overly convenient, this routine builds an unweildy dictionary of satellites with keys of two different Python types: integer keys for satellite numbers, and string keys for satellite names. It even lists satellites like ``ISS (ZARYA)`` twice, in case the user wants to look them up by a single name like ``ZARYA``. What a mess. Users should instead call the simple ``tle_file()`` method, and themselves build any dictionaries they need. See the :meth:`~skyfield.iokit.Loader.open()` documentation for the meaning of the ``reload`` and ``filename`` parameters. rarMN)r( parse_tlemodelsatnum) rEr!rarMdrfnamessatnames r"tlez Loader.tleso  YYs6HY ='l s&)#))""#!D!AdG"+>  > s 9AA%c|j|||5}tt|||cdddS#1swYyxYw)u?Load and parse a TLE file, returning a list of Earth satellites. Given a URL or local path to an ASCII text file, this loads a series of TLE “Two-Line Element” sets and returns a list of :class:`~skyfield.sgp4lib.EarthSatellite` objects for them. See :doc:`earth-satellites`. See the :meth:`~skyfield.iokit.Loader.open()` method for the meaning of the ``reload`` and ``filename`` parameters. See the :meth:`parse_tle_file()` function for the meaning of the ``ts`` and ``skip_names`` parameters. rsN)r(listparse_tle_file)rEr!rarMts skip_namesrfs r"tle_filezLoader.tle_files5 YYs6HY =q"j9:> = =s5>cd|vr|j|}|'t|jjdd}|j |}t |||j ||S)uDownload a file, even if it’s already on disk; return its path. You can specify the local ``filename`` to which the file will be saved; the default is to use the final component of ``url``. Set ``backup`` to ``True`` if you want an already-existing file moved out of the way instead of overwritten. Your operating system may raise any of several errors during a download: hostname lookup failure (this is the usual symptom if you are disconnected from the Internet); the server refusing the connection; and the connection closing mid-download. Skyfield makes no attempt to intercept or interpret these errors — which vary by operating system — so your application itself should catch network errors if it needs to avoid printing raw Python exceptions, or if you want to retry failed downloads. rZrrri)rqrrr rNrjr6)rEr!rMrbrs r"rjzLoader.downloadsd$  ..%C  }))//4R8H||H%dDLL8 r$c d|vr8|}tjj|j|}t ||S|'t |jj dd}|j||||}t ||S)a,Open a file, downloading it first if it does not yet exist. Unlike when you call a loader directly like ``my_loader()``, this ``my_loader.open()`` method does not attempt to parse or interpret the file; it simply returns an open file object. The ``url`` can be either an external URL, or else the path to a file on the current filesystem. A relative path will be assumed to be relative to the base directory of this loader object. If a URL was provided and the ``reload`` parameter is true, then any existing file will be removed before the download starts. The ``filename`` parameter lets you specify an alternative local filename instead of having the filename extracted from the final component of the URL. rZrr)r3rrKr5r(rr r`)rEr!r'rarMrbpath_that_might_be_relativers r"r(z Loader.open.s~&  *- '77<<0KLDd# #  }))//4R8H||C66:D$r$cl|j}|rWtd}|d}|tjt |z }|ddz j d}||f}|d}|d} n|dr=|d r5|d r-|d} |d } |d \} } t | | | \}}} nb|jd }|j|5}tj|\}}d d d tj\}}}} ||f}|td d|g}t||| S#1swYExYw)u(Return a `Timescale` built using official Earth rotation data. ``delta_t`` — Lets you override the standard ∆T tables by providing your own ∆T offset in seconds. For details, see :ref:`custom-delta-t`. ``builtin`` — By default, Skyfield uses ∆T and leap second tables that it carries internally; to instead load this data from files, set this option to ``False``. For compatibility with Skyfield ≤ 1.30, if you have on disk the three files ``deltat.data``, ``deltat.preds``, and ``Leap_Second.dat``, then Skyfield will load them. Otherwise, Skyfield will download and use ``finals2000A.all`` from the International Earth Rotation Service. For details, see :ref:`downloading-timescale-files`. ziers.npztt_jd_minus_arange delta_t_1e7gcA leap_dates leap_offsetsr0r1r.r-Nrr)rXr nparangelenroundr rqr(rparse_dut1_from_finals_allbuild_timescale_arraysr r)rEdelta_trcrGarraysdaily_tt daily_delta_tdelta_t_recentrr deltat_data deltat_preds_leap_second_datr!rfutc_mjddut1s r" timescalezLoader.timescaleLsV" LL %j1F23H  #h-0 0H#M2S8??BM%}4N -J!.1L } !N"3:K8L }-K/L!%&7!8 A7I\?8< 4NJ ..!23C31 $ ? ? B  ++GT: >HmZ%}4N  $aG_5N\BB s D**D3c8dj|jS)N )rKr7)rEs r"logz Loader.log~syy%%r$)TF)FFF)FN)FNNFNF)r&FNF)NT)r_ __module__ __qualname____doc__rHrNrTrXrgr`r^rqr{rrjr(rpropertyrr$r"r+r+Gsi B* X6 ! 6% N 2 ,048%*;&4 <0Cd&&r$r+c|j|}||Stjj|\}}|j|}||D]\}}t ||s|cSy)z.Search a Loader data structure for a filename.N)getr3rsplitextr)mappingrMresultrzextpatternresult2s r"r[r[sj [[ "F    *ID# [[ F  & GWx)!' r$ctjj|}tjj|\}}|dk(r t |St dj |)a8Open a file on your local drive, using its extension to guess its type. This routine only works on ``.bsp`` ephemeris files right now, but will gain support for additional file types in the future. :: from skyfield.api import load_file planets = load_file('~/Downloads/de421.bsp') r/zunrecognized file extension: {})r3rr4rrr\r])rrprs r" load_filersY 77  d #D  &ID# f}4  6==dC DDr$c#Kdx}}|D]G}|jdr.t|dk\r|jdr t|dk\r|jd}t|dk(r#|jdj}|g}n;|jdr&|d d jdj}|g}nd }d }|jd}|jd}t |||}|rKd |vrG|j d \} } | jd } |j | |j | ||f|}|}Jy w)a?Parse a file of TLE satellite element sets. DEPRECATED: this routine is overly complicated, doing extra work to try to guess several ways in which the user might want to look up satellites by name. Use ``parse_tle_file()`` instead. TODO: convert this into a wrapper around ``parse_tle_file()``. r$1 E2 s ascii0 Nrz ()) startswithrrstripdecoderr rl) fileobjb0b1b2rzrxline1line2ry short_namesecondary_names r"rtrts<MB MM% SW] MM% SW]7#B2w"}yy)002u%!"v}}W-446IIg&EIIg&E t4C .2ZZ-=* N!/!6!6s!; Z( ^,*   ?sEEc#Kdx}}|D]}|jdrt|dk\r|jdrt|dk\rv|s;|r9|jd}|jdr|dd}|jd }nd}|jd }|jd }t ||||dx}}|}|}yw) uParse lines of TLE satellite data, yielding a sequence of satellites. Given a sequence ``lines`` of byte strings (which can be an open binary file, which acts like a sequence of lines in Python), this routine yields an :class:`~skyfield.sgp4lib.EarthSatellite` for each pair of adjacent lines that start with ``"1 "`` and ``"2 "`` and have 69 or more characters each. If the line preceding a TLE is not part of another TLE, it is used as the satellite’s ``.name``. If you pass a ``ts`` timescale, Skyfield will use it to build the ``.epoch`` date attribute on each satellite; otherwise a timescale derived from Skyfield’s built-in leap second files will be used. If for a particular file you see random lines of text being interpreted as satellite names, set ``skip_names`` to ``True`` and Skyfield will not try to store satellite names. See :doc:`earth-satellites` for details. An exception is raised if the attempt to parse a pair of candidate lines as TLE lines fails. r$rrrs rrNr)rrrrr) linesrrrrrrzrrs r"r~r~s,MB MM% SW] MM% SW]"YYx(=='ABByy)IIg&EIIg&E tR8 8MBBB'sCCc t,ttj}t||}n1tr t|tj}n t|}|tjj}d} |r~trHtdjtjj!|tjn0t#|} t%|j&j)dd} |d zx} } t+td d tj,ztj.ztj0z} d } tj2| | d } tj:|d5} d } |j=|}|sn4|j?||tA|z }| | jC| H|jE dddtjjG|r |r tI| tK| |y#t $r)}t dj||}d|_|d}~wwxYw#t4$r?}|j6t6j8k7r|d z }dj| |} Yd}~nd}~wwxYwQ#t $r }t dj||d}~wwxYw#1swYxYw#t $r!}t dj| ||d}~wwxYw)aDownload a file from a URL, possibly displaying a progress bar. Saves the output to the file named by `path`. If the URL cannot be downloaded or the file cannot be written, an ``IOError`` is raised. Normally, if the standard error output is a terminal, then a progress bar is displayed to keep the user entertained. Specify `verbose=True` or `verbose=False` to override this behavior. Nr)contextzcannot download {0} because {1}zDownloading {0} ...)filezcontent-lengthrz .downloadO_BINARYrriz{0}{1}wbzerror getting {0} - {1}zerror renaming {0} to {1} - {2})&rcertifiwherer_supports_cafile_argument ExceptionIOErrorr] __cause__sysstderrisatty _running_IDLEprintr3rbasename ProgressBarintheadersrgetattrO_CREATO_EXCLO_RDWRr(r9r:r;fdopenreadwriterreportflushrW_rename_original_replace)r!rr6 blocksizerb ssl_context connectionrGe2barcontent_lengthtempbasetempnameflagsifdwlengthdatas r"rjrjs  ! -0 HK k:J & W]]_=J J **##% C  '..rww/?/?/EFzz #d#C !3!3!7!78H"!MNN,,Hx B A & 3bii ?")) KE A  5%0B  2t  DF!y1 #d)#?JJv~6 GGI " ww~~d 4 u  6==c1E F 2 4ww%,,& FAx3H  4 2 D3::3BC C D  0  7>> dA   ssA#H$IKAJ& K$ I-$II J!"5JJ!& K/K  KKK L'LLctjdD]K}|jdd\}}dj|||}tj j |rKnt j|y)NrrJz{0}.old{1}.{2}) itertoolscountrsplitr]r3rrWrename)rnprefixsuffix backup_paths r"rrIs` __Q S!,&--fa@ ww~~k*   IIdK r$ceZdZdZdZy)rcZtjj||_d|_y)Nr)r3rrrMt0)rErs r"rHzProgressBar.__init__Rs((. r$cF|dkryd|z|z}|dk7rt|jz dkryt|_d|dzz}tdj|||j|dk(rdndt j t j jy) Nrdg?#rz [{0:33}] {1:3}% {2}rr)endr)rrrr]rMrrr)rE bytes_so_far bytes_totalpercentrs r"rzProgressBar.reportVs ?  $ 3 sN$''!1C!7 &W\" %,,S'4==I"c>$ E r$N)r_rrrHrrr$r"rrQs  r$rr)NiF)D __future__rrr3r:rrpkgutilrrrnumpyrrrcurvelibr functionsr io_timescaler r r rjpllibrsgp4librtimelibriorr version_inforeplacerrsslr ImportError urllib.parserurllib.requestrurllib2 TypeErrorrr\r __class__r_rr#r>r=rAr?r@r)objectr+r[rrtr~rjrrrr$r"rs^% ' #-vzzHyyH"* %& % Bt%%..2DD - 1 :6E L!y&Vy&v E *X*XP d!&m-,"!"  !& %% $%sB"D D D'% D: DD$#D$'D7:E E  E