HR-e<_~ddlZddlZddlmZmZddlZddlmZ ddl m Z ddl m Z mZddlmZddlmZmZddlmZdd lmZdd lmZd d lmZmZejd ZejdZgdZ dZ!dZ"de"ze!zZ#dZ$dd%e$dZ&dZ'ddddZ(dZ)dZ*dZ+dZ,d Z-d!Z.d"Z/d#Z0dS)$N) OrderedDict defaultdict)units) EarthLocation)Column MaskedColumn)col_copy)Time TimeDelta)BARYCENTRIC_SCALES)FITS_DEPRECATED_SCALES)AstropyUserWarning)CardHeaderz(?P[A-Z]+)[-]+z(?P[A-Z]+)\s*) sdacyminhyrtaBa)zOBSGEO-XzOBSGEO-YzOBSGEO-Z)zOBSGEO-LzOBSGEO-BzOBSGEO-H)DATEzDATE-AVGzDATE-BEGzDATE-ENDzDATE-OBSDATEREFJDREFzMJD-AVGzMJD-BEGzMJD-ENDzMJD-OBSMJDREFTIMEOFFSTIMESYSTIMEUNITTREFDIRTREFPOS)TCTYPTCUNITRPOS(|z)[0-9]+c:tjt|duS)z Check if the FITS header keyword is a time column-specific keyword. Parameters ---------- keyword : str FITS keyword. N)rematchCOLUMN_TIME_KEYWORD_REGEXP)keywords 8lib/python3.11/site-packages/astropy/io/fits/fitstime.pyis_time_column_keywordr/Fs 8. 8 8 DD)UTCzDefault time scale)gzTime columns are jd = jd1 + jd2) TOPOCENTERzTime reference position)r rr#ctjddd<dtjvrddkr$t jdtdd<dd<ddkr%t jdtd d<d d<n-td dd d<fd tD}t|d krtj |dtjid<n|r-t jd |tfdt D}t|d krtj|d<n4|r-t jd |td d<dD]\}}|vr||dd<d Sdddd<d S)z Given the global time reference frame information, verify that each global time coordinate attribute will be given a valid value. Parameters ---------- global_info : dict Global time reference frame information. r scalegpszGlobal time scale (TIMESYS) has a FITS recognized time scale value "GPS". In Astropy, "GPS" is a time from epoch format which runs synchronously with TAI; GPS is approximately 19 s ahead of TAI. Hence, this format will be used.taiformatlocala+Global time scale (TIMESYS) has a FITS recognized time scale value "LOCAL". However, the standard states that "LOCAL" should be tied to one of the existing scales because it is intrinsically unreliable and/or ill-defined. Astropy will thus use the default global time scale "UTC" instead of "LOCAL".utcNzGlobal time scale (TIMESYS) should have a FITS recognized time scale value (got {!r}). The FITS standard states that the use of local time scales should be restricted to alternate coordinates.c(g|]}|v|Sr;.0attr global_infos r. z'_verify_global_info..s(OOOT4;;N;N{4 ;N;N;Nr0unitlocationzaThe geocentric observatory location {} is not completely specified (X, Y, Z) and will be ignored.c(g|]}|v|Sr;r;r<s r.r@z'_verify_global_info..s(SSSt{?R?R;t$?R?R?Rr0zeThe geodetic observatory location {} is not completely specified (lon, lat, alt) and will be ignored.))rmjd)rjd)rfits)valr7ref_timerrE)r getlowerr SCALESwarningswarnrAssertionErrorr7 OBSGEO_XYZlenrfrom_geocentricum OBSGEO_LBH from_geodetic)r?obs_geokeyformat_s` r._verify_global_inforZZso25I I 6 < < > >K 74;.. w 5 ( ( MA#    $)K $)K ! w 7 * * M> #    $)K $(K ! !! &vk)&<==  !% HPOOOZOOOG 7||q"/"?"Sqs"S"S J   M;;A6'??"   TSSSSSS w<<1  &3&A7&KK # #   EEKVG__& '+K #R>> W +  .9#.>'&R&RK # EE  +,u"="= Jr0c|dd}|dd}|dd}|Ut|ddr$t|ddrdS|t jvr||d<d|d<n|tjvrt||d<d|d<n|d kr|d|d<|d|d<n|d kr8tj d |td |d<d |d<nr|dkrDtj d |t|d|d<|d|d<n(dS| |tvs||d|d<|d|d<ndS|:|dkr.|d|d<|dtj dtn?d|d<n9|ddkr(|dtj dtd|d<nd|d<|d|d<dS)a Given the column-specific time reference frame information, verify that each column-specific time coordinate attribute has a valid value. Return True if the coordinate column is time, or else return False. Parameters ---------- global_info : dict Global time reference frame information. column_info : dict Column-specific time reference frame override information. r$Nr%r&Fr4r7TIMEGPSzTable column "{}" has a FITS recognized time scale value "GPS". In Astropy, "GPS" is a time from epoch format which runs synchronously with TAI; GPS runs ahead of TAI approximately by 19 s. Hence, this format will be used.r6r5LOCALaTable column "{}" has a FITS recognized time scale value "LOCAL". However, the standard states that "LOCAL" should be tied to one of the existing scales because it is intrinsically unreliable and/or ill-defined. Astropy will thus use the global time scale (TIMESYS) as the default.r2rCzTime column reference position "TRPOSn" value is "TOPOCENTER". However, the observatory position is not properly specified. The FITS standard does not support this and hence reference position will be ignored.r#zTime column reference position "TRPOSn" is not specified. The default value for it is "TOPOCENTER", and the observatory position has been specified. However, for supporting column-specific location, reference position will be ignored for this column.rIT)rJ TCTYP_RE_TYPEr+ TCTYP_RE_ALGOrKr rLr keysrMrNr7rFITS_TIME_UNIT) column_infor?r4rBrCs r._verify_column_infores OOGT * *E ??7D ) )Dw--H    uRaRy ) )+ m.A.A%).L.L+ 5 [[]]dk ) )#(;;==K $(K ! ! ,133 3 3#9%#@K $(K ! !f__#.w#7K $/$9K ! ! e^^ M<=CF; * , ,C 2122 677"8H>Hctddd}tt}|d}|jD]\}}}|tvr|||<||*t |rXtjd| \}} ||t| |<|||tvr tjd|r ||||<t|t|||rS|D]>\} } t| |r)|j| dz } t#|| | || <?t%|jD]=\} } | dz|vr1t'|| |} | rt#|| | || <>|S) as Read FITS binary table time columns as `~astropy.time.Time`. This method reads the metadata associated with time coordinates, as stored in a FITS binary table header, converts time columns into `~astropy.time.Time` columns and reads global reference times as `~astropy.time.Time` instances. Parameters ---------- hdr : `~astropy.io.fits.header.Header` FITS Header table : `~astropy.table.Table` The table whose time columns are to be read as Time Returns ------- hdr : `~astropy.io.fits.header.Header` Modified FITS Header (time metadata removed) r1r2)r r#T)stripz([A-Z]+)([0-9]+)z TTYPE[0-9]+r)rrcopycards TIME_KEYWORDSremover/r*r+groupsrrPrZrsrmrecolnamesr enumeraterk) hdrrqr? time_columnshcopyrXrrcommentbaseidxrdcolnames r. fits_to_timers,$ ==K{++L HH4H E #y . .UG -  $K  LL     #C ( ( .!4c::AACCID#+0LS "4 ( LL     j RXmS%A%A !&uK  $$$ ,,,S , 2 2 4 4 S S C"; << S.q1!5eGnk!R!Rg"%.11SS W !GL ( (25>;OOK S!5eGnk!R!Rg Lr0c v g}|D]b}t|tr|d}n|jjrt |dn|}||c||d}|j |_ tdt D}tt|j d<|j d}|jt }d |D]}|jrt$nt}|jjdd kr/||jj||j]t/j|j|jgd } ||jj|| d |j||jjd<d ||jjd<|jPd||jjd< 7t=jd |jjtBBd||jjd<|jtDvrOt=jd |jj|jtB |j j#dkrkdD]f} |$ttK | &tNj(d| g3|) fddDVt/j* |jkr(tWd |j||fS)as Replace Time columns in a Table with non-mixin columns containing each element as a vector of two doubles (jd1, jd2) and return a FITS header with appropriate time coordinate keywords. jd = jd1 + jd2 represents time in the Julian Date format with high-precision. Parameters ---------- table : `~astropy.table.Table` The table whose Time columns are to be replaced. Returns ------- table : `~astropy.table.Table` The table with replaced Time columns hdr : `~astropy.io.fits.header.Header` Header containing global time reference frame FITS keywords F) copy_data) copy_indices)rcPg|]#\}}t||d|d$S)rr)r-rrr)r)r=rXrHs r.r@z time_to_fits..<sC   S CFCF ; ; ;   r0__coordinate_columns__NrGformatted_valuerv)axisr)rB coord_type coord_unit time_ref_poszTime Column "{}" has no specified location, but global Time Position is present, which will be the default for this column in FITS specification.r2zQEarth Location "TOPOCENTER" for Time Column "{}" is incompatible with scale "{}".r)xyzOBSGEO-)rhc g|]S}td|t|tjTS)r)r-rr)rrigetattrto_valuerSrT)r=dimrCs r.r@z time_to_fits..~sj !$ !(?#))++(?(?&-h&<&<&E&Eac&J&Jr0zMultiple Time Columns with different geocentric observatory locations ({}, {}) encountered.This is not supported by the FITS standard.),itercols isinstancerrrgindicesr append __class__rnrGLOBAL_TIME_INFOrmrrcolumnsr maskedrserialize_methodreplace_columnrhrrnpstackjd1jd2r4rirCrMrNr7rr size add_columnrrrSrTextendanyrp) rqnew_colsrjnew_colnewtabler coord_meta time_colscol_clsjd12rrCs @r. time_to_fitsrs2H~~!! c6 " " Uhhh//GG;>8;KThs7777QTG    xe44HJHM   ,2244     C/:+.F.FHM*+78J ((..IHBB#&*8,,& 8 $V ,0A A A  # #CHM77393E3E F F F x#'*444  wwt#/F/F/FGGG36)//2C2C 38=!,/25 38=!,/ < 8rs> 00000000------........))))))((((((((000000777777777777 233  122 EDD2 1 &'()43I*>!?!?HHH E E E- 58d>d>d>Nmmm`@,CCC@IIIXAAAHxxxxxr0