[c@`sdZddlmZmZmZddlZddlZddlZ ddl m Z m Z ddl mZddl mZdd lmZmZmZdd lmZmZmZd ejfd YZd ejfdYZdejfdYZdejfdYZdS(s An extensible ASCII table reader and writer. Classes to read DAOphot table format :Copyright: Smithsonian Astrophysical Observatory (2011) :Author: Tom Aldcroft (aldcroft@head.cfa.harvard.edu) i(tabsolute_importtdivisiontprint_functionN(t defaultdictt OrderedDicti(tcore(t fixedwidthi(tziptmaptrange(tfirst_true_indextfirst_false_indext groupmoret DaophotHeadercB`skeZdZdZejdZejdejZd Z dZ dZ dZ dZ dZRS( sK Read the header from a file produced by the IRAF DAOphot routine. s\s*#Ks%-?(\d+)\.?\d?[sdfg]s,[#]K\s+ (?P \w+)\s* = (?P .+) $cC`stjj|dS(N(Rt BaseHeadert__init__(tself((s7lib/python2.7/site-packages/astropy/io/ascii/daophot.pyR'sc`s4d }tt}d}x`tt|j|D]F}x=t|t||D]#\}}||j|jqSWq4Wjjrnttt|j|d\}} } t jj } xt j d| dj dD]} ttdjt|tj| } | d kr4| |dd7siitU2tt1c`s tjj|jdS(Ni(tintt re_formattsearchtgroups(R(R(s7lib/python2.7/site-packages/astropy/io/ascii/daophot.pyRNscs`s*|] \}}|t|gfVqdS(N(tsum(t.0tktv((s7lib/python2.7/site-packages/astropy/io/ascii/daophot.pys [s(s#Ns#Us#F(RtlistRRtgettappendtsplittdatat is_multilinetlent first_blocktnptarangetastypetjointitttrepeattfromiterRRtDaophott table_widtht col_widthstdicttitems(Rtgrouped_lines_dicttline_idst coldef_dicttstrippertdefblocktkeytlinet last_namest last_unitst last_formatst N_multilinetitextended_namest get_col_widthtformatstfR0t row_widthst row_shorttwtr((Rs7lib/python2.7/site-packages/astropy/io/ascii/daophot.pytparse_col_defs*s4   % +&* 3  cC`s+|d}t|j}|dkr'd}tt||jt|\}}t|\}} tt||} d| krtt|j| d} | |d Extract info from a header keyword line (#K) tstuffiitunitsiRLiRtvalueRXN(tre_header_keywordtmatchtgroupRtrsplitRIR%(RR9tmtvalst keyword_dict((s7lib/python2.7/site-packages/astropy/io/ascii/daophot.pyROs! &c `s|jstjdn|j|jd}x`|jD]Ut|jjd \}}|dkry|_ n|dkr<|_ q<q<Wt |j g}t j|}||}xt|jD]\}||||__jj_tdrtfddDrDtj_qdj krbtj_qd j krtj_qqqW|jjjdd S(s Initialize the header Column objects from the table ``lines`` for a DAOphot header. The DAOphot header is specialized so that we just copy the entire BaseHeader get_cols routine and modify as needed. Parameters ---------- lines : list List of table lines Returns ---------- col : list List of table Columns s'No column names found in DAOphot headerRMRKRLs##c3`s|]}|jkVqdS(N(RL(Rtx(tcol(s7lib/python2.7/site-packages/astropy/io/ascii/daophot.pys stfgtdRtINDEFt0N(RKRL(RiRj(RQRtInconsistentTableErrort_set_cols_from_namesRPRMRRXR RKRLRR0R'tcumsumt enumeratetstarttendtspanthasattrtanyt FloatTypettypetIntTypetStrTypeR#t fill_valuesR!( RRNtcoldefsRKRYt col_widthtendststartsR>((Rfs7lib/python2.7/site-packages/astropy/io/ascii/daophot.pytget_colss0       ((t__name__t __module__t__doc__tcommenttretcompileRtVERBOSER^taperture_valuesRRGRZROR}(((s7lib/python2.7/site-packages/astropy/io/ascii/daophot.pyR s    5 ( t DaophotDatacB`s/eZejZdZdZdZdZRS(is\s*#cC`stjj|t|_dS(N(RtBaseDataRtFalseR$(R((s7lib/python2.7/site-packages/astropy/io/ascii/daophot.pyRscC`s_|jrHttttj|j}ttt||j _ nt j j ||dS(N(R$tnextRRtstrR"R&ttupletfloattheaderRRRtget_data_lines(RRNtaplist((s7lib/python2.7/site-packages/astropy/io/ascii/daophot.pyRs !( R~RRtFixedWidthSplittertsplitter_classt start_lineRRR(((s7lib/python2.7/site-packages/astropy/io/ascii/daophot.pyRs   tDaophotInputtercB`s>eZdZdZdZejdZddZdZ RS(s\t*t s(#?)[^\\*#]*(\*?)(\\*) ?$ic `stfd|| D\}}}t|}|dkrRdd|| fS|| }t|||!}|dkrdd|fSt||||!} tj||| g} || d| d!} | | |fS(s Search lines for special continuation character to determine number of continued rows in a datablock. For efficiency, depth gives the upper limit of lines to search. c3`s'|]}jj|jVqdS(N(t re_multilineRR(Rtl(R(s7lib/python2.7/site-packages/astropy/io/ascii/daophot.pys siiN(RR RIR R'Rm( RRNtdepthRtspecialtcontt data_startt header_linest first_specialt last_specialtmarkerstmultiline_block((Rs7lib/python2.7/site-packages/astropy/io/ascii/daophot.pytsearch_multilines      cC`s|j|\}}}|dk |j_||j_||j_||jj_|dk rn||d}n|j}|j }|j }g}g} xt |D]\} } |j j | } | ri| j\} }}| s|r| j||} n|r| j||} n|r4| r4|j| n|s|j| | jdj|g}qqtdj| | qW| S(NiRs(multiline re could not match line {}: {}(RRIR#R$RR&RRNtcontinuation_chartmultiline_chart replace_charRnRRRtreplaceR!R*t ValueErrorRL(RRNRtblockRRRRtpartstoutlinesR>R9tmoRRR((s7lib/python2.7/site-packages/astropy/io/ascii/daophot.pyt process_liness:           ( R~RRRRRRRRR(((s7lib/python2.7/site-packages/astropy/io/ascii/daophot.pyRs  0R.cB`sVeZdZdZdgZeZdZeZ e Z e Z dZdZddZRS(s^ Read a DAOphot file. Example:: #K MERGERAD = INDEF scaleunit %-23.7g #K IRAF = NOAO/IRAFV2.10EXPORT version %-23s #K USER = davis name %-23s #K HOST = tucana computer %-23s # #N ID XCENTER YCENTER MAG MERR MSKY NITER \ #U ## pixels pixels magnitudes magnitudes counts ## \ #F %-9d %-10.3f %-10.3f %-12.3f %-14.3f %-15.7g %-6d # #N SHARPNESS CHI PIER PERROR \ #U ## ## ## perrors \ #F %-23.3f %-12.3f %-6d %-13s # 14 138.538 INDEF 15.461 0.003 34.85955 4 \ -0.032 0.802 0 No_error The keywords defined in the #K records are available via the output table ``meta`` attribute:: >>> import os >>> from astropy.io import ascii >>> filename = os.path.join(ascii.__path__[0], 'tests/t/daophot.dat') >>> data = ascii.read(filename) >>> for name, keyword in data.meta['keywords'].items(): ... print(name, keyword['value'], keyword['units'], keyword['format']) ... MERGERAD INDEF scaleunit %-23.7g IRAF NOAO/IRAFV2.10EXPORT version %-23s USER name %-23s ... The unit and formats are available in the output table columns:: >>> for colname in data.colnames: ... col = data[colname] ... print(colname, col.unit, col.format) ... ID None %-9d XCENTER pixels %-10.3f YCENTER pixels %-10.3f ... Any column values of INDEF are interpreted as a missing value and will be masked out in the resultant table. In case of multi-aperture daophot files containing repeated entries for the last row of fields, extra unique column names will be created by suffixing corresponding field names with numbers starting from 2 to N (where N is the total number of apertures). For example, first aperture radius will be RAPERT and corresponding magnitude will be MAG, second aperture radius will be RAPERT2 and corresponding magnitude will be MAG2, third aperture radius will be RAPERT3 and corresponding magnitude will be MAG3, and so on. tdaophotsIRAF DAOphot format tableiPcC`s#tjj||j|j_dS(N(Rt BaseReaderRR#tinputter(R((s7lib/python2.7/site-packages/astropy/io/ascii/daophot.pyRscC`s tdS(N(tNotImplementedError(RRH((s7lib/python2.7/site-packages/astropy/io/ascii/daophot.pytwritesN(R~RRt _format_namet_io_registry_format_aliasesRt_io_registry_can_writet _descriptionR t header_classRt data_classRtinputter_classR/RRIR(((s7lib/python2.7/site-packages/astropy/io/ascii/daophot.pyR.=s<  ( Rt __future__RRRRtnumpyR't itertoolsR+t collectionsRRRRRtextern.six.movesRRR tmiscR R R RR RRtContinuationLinesInputterRRR.(((s7lib/python2.7/site-packages/astropy/io/ascii/daophot.pyt s   ^