B {dc @s8dZddlmZddlZddlZedZddlZddlm m Z ddl Z ddlmZddlmZddlmZdd lmZdd lmZdd lmZd d dgd d dgddddddgdZeejfZeddGdddeZddZd/ddZd0d!d"Z ed#ddd1d$d%Z!d2d'd(Z"d3d)d*Z#d4d+d,Z$d-d.Z%dS)5zCProvides input and output functions for Healpix maps, alm, and cl. )divisionNhealpy)deprecated_renamed_argument) pixelfunc)Alm)UNSEEN)cookbook) deprecated TEMPERATUREZQ_POLARISATIONZU_POLARISATIONZIIZIQZIUZQQZQUZUU)rz1.15.0)Zsincec@s eZdZdS)HealpixFitsWarningN)__name__ __module__ __qualname__rr\/work/yifan.wang/ringdown/master-ringdown-env/lib/python3.7/site-packages/healpy/fitsfunc.pyr2srcsrd}t|trt|}d}t|ddtfddttj D}|rV| t|dkrj|dS|SdS) zReads Cl from a healpix file, as IDL fits2cl. Parameters ---------- filename : str or HDUList or HDU or pathlib.Path instance the fits file name Returns ------- cl : array the cl array FTr)hducsg|]}j|qSr)datafield).0n)fits_hdurr Iszread_cl..rN) isinstance allowed_pathspfopen_get_hdunparrayrangelencolumnsclose)filename opened_fileclr)rrread_cl7s    " r*Fcs|dkr$t|tjr|jn|dj}t|ddddddg}tt|d krtfd d t|dt||D}n2tt|d krtj dd |dg}nt dtj |}d|j d<|jt||ddS)aEWrites Cl into a healpix file, as IDL cl2fits. Parameters ---------- filename : str the fits file name cl : array the cl array to write to file dtype : np.dtype (optional) The datatype in which the columns will be stored. If not supplied, the dtype of the input cl will be used. This changed in `healpy` 1.15.0, in previous versions, cl by default were saved in `float64`. overwrite : bool, optional If True, existing file is silently overwritten. Otherwise trying to write an existing file raises an OSError. Nrr ZGRADIENTZCURLzG-TzC-TzC-Gr cs$g|]\}}tj|d|dqS)z%s)nameformatr")rColumn)rZ column_nameZ column_cl) fitsformatrrrkszwrite_cl..rz%s)r+r,r"z6write_cl: Expected one or more vectors of equal lengthrZCREATOR) overwrite)rr!Zndarraydtype getformatr$shapeziprr- RuntimeError BinTableHDU from_columnsheaderwritetostr)r'r)r0r/ column_namescolstbhdur)r.rwrite_clRs   r=Trc  s^t|dstdt|}t|dkr0|g}|dkrXdd|D}tdt|y$g} x|D]} | t | qdWWn&tk rt |gt |} YnX|dkrt t |ddt d t |d D}nt |t |kstd |dkst|tr|gt |}t ttt |d ks2td tt |d} | dkrVtd g}|rd }t|dtd}t |dkrtdfdd|D}t t| }|dkrd}|tjd||ddxt|||| D]t\}}}}t |dkrH|rHt|}|tj|d|||jdd|dn|tj|d|||dqWtj |}d|j!d<|rd}nd}|df|j!d<|r|df|j!d<d|j!d<| d f|j!d!<|sd"|j!d#<t"| d d$f|j!d%<|rd&nd'd(f|j!d)<|rd*nd+d,f|j!d-<x$| D]}|d d|j!|d<q(W|j#t|| d.dS)/a Writes a healpix map into a healpix FITS file. .. warning:: Starting from healpy 1.15.0, if you do not specify `dtype`, the map will be written to disk with the same precision it is stored in memory. Previously, by default `healpy` wrote maps in `float32`. To reproduce the same behaviour of `healpy` 1.14.0 and below, set `dtype=np.float32`. Parameters ---------- filename : str the fits file name m : array or sequence of 3 arrays the map to write. Possibly a sequence of 3 maps of same size. They will be considered as I, Q, U maps. Supports masked maps, see the `ma` function. nest : bool, optional If True, ordering scheme is assumed to be NESTED, otherwise, RING. Default: RING. The map ordering is not modified by this function, the input map array should already be in the desired ordering (run `ud_grade` beforehand). fits_IDL : bool, optional If True, reshapes columns in rows of 1024, otherwise all the data will go in one column. Default: True coord : str The coordinate system, typically 'E' for Ecliptic, 'G' for Galactic or 'C' for Celestial (equatorial) partial : bool, optional If True, fits file is written as a partial-sky file with explicit indexing. Otherwise, implicit indexing is used. Default: False. column_names : str or list Column name or list of column names, if None here the default column names based on the number of columns: 1 : "TEMPERATURE", 2 : ["Q_POLARISATION", "U_POLARISATION"], 3 : ["TEMPERATURE", "Q_POLARISATION", "U_POLARISATION"], 6 : ["II", "IQ", "IU", "QQ", "QU", "UU"] COLUMN_1, COLUMN_2... otherwise (FITS is 1-based) column_units : str or list Units for each column, or same units for all columns. extra_header : list Extra records to add to FITS header. dtype: np.dtype or list of np.dtypes, optional The datatype in which the columns will be stored. Will be converted internally from the numpy datatype to the fits convention. If a list, the length must correspond to the number of map arrays. Default: use the data type of the input array(s) .. note:: this changed in 1.15.0, previous versions saved in float32 by default overwrite : bool, optional If True, existing file is silently overwritten. Otherwise trying to write an existing file raises an OSError (IOError for Python 2). __len__zThe map must be a sequencerNcSsg|] }|jqSr)r0)rxrrrrszwrite_map..z"setting the output map dtype to %scSsg|] }d|qS)z COLUMN_%dr)rrrrrrsrz%Length column_names != number of mapszMaps must have same lengthz+Invalid healpix map : wrong number of pixelFz'Invalid healpix map : empty partial mapcsg|] }|qSrr)rmm)maskrrrsIZPIXEL)r+r,r"unitiz1024%sz%s)ZHEALPIXzHEALPIX pixelisationZPIXTYPENESTEDRINGz,Pixel ordering scheme, either RING or NESTEDORDERINGz,Ecliptic, Galactic or Celestial (equatorial)ZCOORDSYS)Zxtensionz#name of this binary table extensionZEXTNAMEzResolution parameter of HEALPIXNSIDE)rzFirst pixel # (0 based)ZFIRSTPIXzLast pixel # (0 based)ZLASTPIXEXPLICITIMPLICITzIndexing: IMPLICIT or EXPLICITINDXSCHMPARTIALFULLSKYz'Sky coverage, either FULLSKY or PARTIALOBJECT)r/)$hasattr TypeErrorrZ ma_to_arrayZmaptypelogwarningr9appendr1r$standard_column_namesgetr#AssertionErrorrsetmap npix2nside ValueErrorZ mask_goodr!whereZmin_scalar_typemaxrr-r3ZasarrayZreshapesizer5r6r7 nside2npixr8)r'mnestr0Zfits_IDLZcoordpartialr:Z column_unitsZ extra_headerr/r. curr_dtypensider;pixffZcncur@Zcurr_fitsformatZmm2r<orderingargsr)rAr write_map|sB   &                  rhverbosec  sVd} t|tr tj||d}d} t|||d} | jd} | dkrftd| j d} t | j } nt| } td | t | std | jd d } | d kr|rd pd} td| td| t | }g}t| jdd }|d kr|dkrd}n|dkrd}| jdd }|d krt|dkrT|dkrNtdd}n |dkrt|dkrptdd}|d kr|rdpd}td|td||dkrtt| j jd|}t|dst|ts|f}|r|tdd|D}y| j djtdd} Wn^tjk rz}z| dkr>t '| t"j(|j t"j)d"}||}td$n<|sz| d krzt *| t"j(|j t"j)d"}||}td%yt!|t +|<Wnt,k rYnX|-|qW|rg}x&| j.D]\}}|-||fqW| r|/t|dkr|r|d|fS|dSn4t0fd&dDr@t"1|}|rN||fS|SdS)'a Read a healpix map from a fits file. Partial-sky files, if properly identified, are expanded to full size and filled with UNSEEN. .. warning:: Starting from healpy 1.15.0, if you do not specify `dtype`, the map will be read in memory with the same precision it is stored on disk. Previously, by default `healpy` wrote maps in `float32` and then upcast to `float64` when reading to memory. To reproduce the same behaviour of `healpy` 1.14.0 and below, set `dtype=np.float64` in `read_map`. Parameters ---------- filename : str or HDU or HDUList or pathlib.Path instance the fits file name field : int or tuple of int, or None, optional The column to read. Default: 0. By convention 0 is temperature, 1 is Q, 2 is U. Field can be a tuple to read multiple columns (0,1,2) If the fits file is a partial-sky file, field=0 corresponds to the first column after the pixel index column. If None, all columns are read in. dtype : data type or list of data types, optional Force the conversion to some type. Passing a list allows different types for each field. In that case, the length of the list must correspond to the length of the field parameter. If None, keep the dtype of the input FITS file Default: Preserve the data types in the file nest : bool, optional If True return the map in NEST ordering, otherwise in RING ordering; use fits keyword ORDERING to decide whether conversion is needed or not If None, no conversion is performed. partial : bool, optional If True, fits file is assumed to be a partial-sky file with explicit indexing, if the indexing scheme cannot be determined from the header. If False, implicit indexing is assumed. Default: False. A partial sky file is one in which OBJECT=PARTIAL and INDXSCHM=EXPLICIT, and the first column is then assumed to contain pixel indices. A full sky file is one in which OBJECT=FULLSKY and INDXSCHM=IMPLICIT. At least one of these keywords must be set for the indexing scheme to be properly identified. hdu : int, optional the header number to look at (start at 0) h : bool, optional If True, return also the header. Default: False. verbose : bool, deprecated It has no effect memmap : bool, optional Argument passed to astropy.io.fits.open, if True, the map is not read into memory, but only the required pixels are read when needed. Default: False. Returns ------- m | (m0, m1, ...) [, header] : array or a tuple of arrays, optionally with header appended The map(s) read from the file, and the header if *h* is True. F)memmapT)rrjrGNz6No NSIDE in the header file : will use length of arrayrz NSIDE = %dzWrong nside parameter.rFZUNDEFrDrEz.No ORDERING keyword in header file : assume %szORDERING = %s in fits filerMrKrLrJrHzIncompatible INDXSCHM keywordrIz-No INDXSCHM keyword in header file: assume %sz INDXSCHM = %srr>css$|]}t|tr|n|dVqdS)rN)rr9)rfrrr szread_map..)copyz&Trying to fix a badly formatted headerZfixz:The number of dtypes are not equal to the number of fields)r0z#nside={0:d}, sz={1:d}, m.size={2:d}zOrdering converted to NESTzOrdering converted to RINGc3s|]}|dkVqdS)rNr)rdt)r0rrrls)2rrrrr r7rTrPinforrrrXr\intZ isnsideokrYstripr]r9r#r$r%rNtupleastypeZravelZ VerifyErrorrQverifyrUrOr3rr!Zonesr0Zisnpixokr,Z nest2ringZarangeint32Z ring2nestZmask_bad OverflowErrorrRitemsr&allr")r'rr0r_r`rhrirjr(rrbrcrfszretobjZschmerdrar^Zmnewidxr7keyvaluer)r0rread_map,sD                         ,       ,$      rcCst|s|g}tjt|d|d}|dkr@||kr@tdn |dkrL|}|dkrX|}|dkrd|}||krp|}|dkr|djj}t|\}} t ||k| |k} || }| | } tj ||| d} |d|| d} t } x|D]}t jt| d d |fd |fgd }| |d <|j| |d <|j| |d <t jd tt jd|d d}t jd t|d|d d}t jd t|d|d d}t j|||g}| |qW| jt||ddS)aWrite alms to a fits file. In the fits file the alms are written with explicit index scheme, index = l*l + l + m +1, possibly out of order. By default write_alm makes a table with the same precision as the alms. If specified, the lmax and mmax parameters truncate the input data to include only alms for which l <= lmax and m <= mmax. Parameters ---------- filename : str The filename of the output fits file alms : array, complex or list of arrays A complex ndarray holding the alms, index = m*(2*lmax+1-m)/2+l, see Alm.getidx lmax : int, optional The maximum l in the output file mmax : int, optional The maximum m in the output file out_dtype : data type, optional data type in the output file (must be a numpy dtype). Default: *alms*.real.dtype mmax_in : int, optional maximum m in the input array r)mmaxrzToo big lmax in parameterN)lr^r r)indexirealimag)r0rz l*l+l+m+1)r+r,rCr"unknown)r/)cbZ is_seq_of_seqrZgetlmaxr$rYrr0Zgetlmr!rZgetidxrHDUListemptyrr-r1rur5r6rRr8r9)r'almsZ out_dtypelmaxrZmmax_inr/Zl2maxrr^r~Zidx_in_originalrhdulistalmZout_dataZcindexZcrealZcimagr<rrr write_alms\       rcsZg}d}d}d}ttr(td}xt|D]܉fddtdD\}}} tt|d t } || d| d} | d k rt d | } | } |dkr| }| }n"|| ks|| krtd ||d }t| | | }||j|<| |j|<||q4W|r"t|dkr:|d }n t|}|rR|| fS|SdS) aRead alm from a fits file. In the fits file, the alm are written with explicit index scheme, index = l**2+l+m+1, while healpix cxx uses index = m*(2*lmax+1-m)/2+l. The conversion is done in this function. Parameters ---------- filename : str or HDUList or HDU or pathlib.Path instance The name of the fits file to read hdu : int, or tuple of int, optional The header to read. Start at 0. Default: hdu=1 return_mmax : bool, optional If true, both the alms and mmax is returned in a tuple. Default: return_mmax=False Returns ------- alms[, mmax] : complex array or tuple of a complex array and an int The alms read from the file and optionally mmax read from the file NFTcs g|]}tdj|qS))r)r rr)rr)r'rCrrrqszread_alm..r rr rzNegative m value encountered !z?read_alm: harmonic expansion order in {} HDUs {} does not matchy)rrrrr!Z atleast_1dr#floorsqrtrsrpanyrYr[r4r,rrrrrRr&r$r")r'rZ return_mmaxrZlmaxtotZmmaxtotr(r~ZalmrZalmirr^rrrrr)r'rCrread_almQsF        rcCst|tr$tj||d}t||dSt|tjrbt|trX|t|krXtdt|q||}n,t|tj tj tj tj tj fr|}ntd|S)a  Return an HDU from a FITS file Parameters ---------- input_data : str or HDUList or HDU instance The input FITS file, either as a filename, HDU list, or HDU instance. Returns ------- fits_hdu : HDU The extracted HDU )rj)rzAvailable hdu in [0-%d]z^First argument should be a input_data (str or pathlib.Path), HDUList instance, or HDU instance)rrrrr rrpr$rYZ PrimaryHDUZImageHDUr5ZTableHDUZ GroupsHDUrO)Z input_datarrjrrrrrr s    r cCsFttjdttjdttjdttjdttjdttjdttjdttj dttj d i }y||kr||SWn YnXy t||kr|t|SWn YnXy(tt ||kr|tt |SWn YnXy2tt |d |kr&|tt |d SWn YnXy|t krFd SWn YnXyt |t krrd t |SWn YnXy0t |d t krtd d|D}d |SWn YnXy&t|j |kr|t|j SWn YnXy.t|d j |kr |t|d j SWn YnXtdt |dS)a'Get the FITS convention format string of data type t. Parameters ---------- t : data type The data type for which the FITS type is requested Returns ------- fits_type : str or None The FITS string code describing the data type, or None if unknown type. LBrBJKEDCMrAzA%dcss|]}t|VqdS)N)r$)rsrrrrlszgetformat..zthealpy could not understand the equivalent FITS datatype of {}, please open an issue on the healpy Github repositoryN)r!r0Zbool_Zuint8Zint16ruZint64Zfloat32Zfloat64Z complex64Z complex128typer9r$r[rYr,)tconvrrrrr1st           r1)NF) FNTNFNNrF)rNFFrFTF)NrrrF)rF)NN)&__doc__ __future__rsyslogging getLoggerrPpathlibZastropy.io.fitsioZfitsrnumpyr!Zastropy.utils.decoratorsrrZsphtfuncrrr rr rSr9PathrWarningrr*r=rhrrrr r1rrrrs\          - &  B Z D $