B Qa %@sddlZddlZddlZddlZddlmZddlZddlZddl Z ddl m Z m Z mZmZmZmZmZdZdgdgdddgd d d gd gd dgdddgddgddgddgddgddddddd gd!gd"d#gd$gd%gd&d'd(gd)d*d+d,gd-d.d/d0d1gd2gd3gd4gd5gd6gd7gd8gd9gd:gd;dgd?d@gdAdBgdCgdDgdEgdFgdG$ZdHdIdJdKdLdMdNdOdPdQ ZdRdIdSdKdTdUdVdWdPdQ ZdXdXdYdZd[d\d]d^d_d`da ZdbdcgZdRdSdTdUdVdWdPdIdKdddedf ZdgdhdidjZdkdldmZddodpZddqdrZddsdtZddvdwZddxdyZ dzd{Z!dd|d}Z"dd~dZ#dddZ$dddZ%dddZ&ddZ'ddZ(ddZ)ddZ*dddZ+ddZ,dddZ-dddZ.dddZ/dddZ0dddZ1dddZ2dddZ3dddZ4dddZ5ddZ6dS)N)datasetUnitDictgeometrygiantIfgramStackgiantTimeseries ifgramStack timeseriesHDFEOSgJxޱAZ azimuth_looksZ range_looksZazimuthPixelSizeZazimuth_pixel_spacingZaz_pixel_spacingZrangePixelSizeZrange_pixel_spacingZrg_pixel_spacingZ center_timeZdataType data_typeZ earthRadiusearth_radius_below_sensor earth_radiusZ HEADING_DEGheadingaltitude SC_heightZ number_bandsZbandsscheme interleavelengthZ FILE_LENGTHlines azimuth_linesnlinesaz_sampinterferogram_azimuth_linesZ passDirectionZspacecraftNamesensorZ polarizationZprfZ startingRangeZnear_range_slcZ near_rangeZ wavelengthZ WavelengthZradarWavelengthZradar_wavelengthwidthZWidthZsamples range_sampinterferogram_widthref_dateZref_latZref_lonref_xref_yZ subset_x0Z subset_x1Z subset_y0Z subset_y1Z corner_lonZ corner_eastZ corner_latZ corner_northZpost_lonZ post_eastZpost_latZ post_northZ swathNumberZfirstFrameNumberZlastFrameNumberZ trackNumber)$ZALOOKSZRLOOKSAZIMUTH_PIXEL_SIZERANGE_PIXEL_SIZECENTER_LINE_UTC DATA_TYPE EARTH_RADIUSZHEADINGHEIGHTBANDS INTERLEAVELENGTHORBIT_DIRECTIONPLATFORMZ POLARIZATIONPRFZSTARTING_RANGE WAVELENGTHWIDTHREF_DATEREF_LATREF_LONREF_XREF_Y SUBSET_XMINZ SUBSET_XMAX SUBSET_YMINZ SUBSET_YMAXX_FIRSTY_FIRSTX_STEPY_STEPZ beam_swathZ first_frameZ last_frameZrelative_orbitBYTEuint16SHORTuint32INTFLOATDOUBLECFLOAT complex128)  uint8int16int32float32float64 complex64r@rArBrCrDrErFrGrH) rIint8r8rJr:rKrLrMrNr?z.tifz.grdint64uint64) 1234569Z12Z13Z14Z15BSQBILBIP)ZBANDZLINEPIXELz little-endianz big-endian)0rRTc Csd}t|tr"|ddd}nt|tr:|dd}t||d}t|dt|d} } |spdd| | f}tjtj |d } | dkrt ||||||d } nt |||||d \} }|rt j| |d } | |fS) adRead one dataset and its attributes from input file. Parameters: fname : str, path of file to read datasetName : str or list of str, slice names box : 4-tuple of int area to read, defined in (x0, y0, x1, y1) in pixel coordinate x/ystep : int, number of pixels to pick/multilook for each output pixel data_type : numpy data type, e.g. np.float32, np.bool_, etc. Returns: data : 2/3-D matrix in numpy.array format, return None if failed atr : dictionary, attributes of data, return None if failed Examples: from mintpy.utils import readfile data, atr = readfile.read('velocity.h5') data, atr = readfile.read('timeseries.h5') data, atr = readfile.read('timeseries.h5', datasetName='timeseries-20161020') data, atr = readfile.read('ifgramStack.h5', datasetName='unwrapPhase') data, atr = readfile.read('ifgramStack.h5', datasetName='unwrapPhase-20161020_20161026') data, atr = readfile.read('ifgramStack.h5', datasetName='coherence', box=(100,1100, 500, 2500)) data, atr = readfile.read('geometryRadar.h5', datasetName='height') data, atr = readfile.read('geometryRadar.h5', datasetName='bperp') data, atr = readfile.read('100120-110214.unw', box=(100,1100, 500, 2500)) Nr-) datasetNamer&r+r@)z.h5z.he5)r_boxxstepystep print_msg)r_r`rarb)dtype) isinstancelistsplitstrread_attributeintospathsplitextbasenamelowerread_hdf5_fileread_binary_filenparray) fnamer`r_rcrarbr Z dsname4atratrrrfextdatarxK/home/centos/operations/rsmas_insar/sources/MintPy/mintpy/utils/readfile.pyreads2      rzc st|}g}x(dd|DD]}||kr||qWdd|Dfdd|Dsf|dgnttrvgtddDrfddDdd dfd dD} d d| D} t|d $fd ddgD} fdd|D} t| dkr$| d} n,t| dkr@| d} nt d || j dkr| |d|d|d|df} ||dkrt |d|d|}t |d|d|}| t |dd|t |dd|f} | d|d|f} n| j dkrt j| jdt jd}| r2| dgkr@d|dd<n:ddfdd|DD}x| D]}d|||<qbW||dkr| ||d|d|d|df} n(t |d|d|}t |d|d|}t |}t |||f| j} t |d}xt|D]}|rNtjdd|dd|dtj| |||d|d|d|df}|t |dd|t |dd|f}|d|d|f| |ddddf<qW|rtdtdd| jDrt | } WdQRX| S)a Parameters: fname : str, name of HDF5 file to read datasetName : str or list of str, dataset name in root level with/without date info 'timeseries' 'timeseries-20150215' 'unwrapPhase' 'unwrapPhase-20150215_20150227' 'HDFEOS/GRIDS/timeseries/observation/displacement' 'recons' 'recons-20150215' ['recons-20150215', 'recons-20150227', ...] '20150215' 'cmask' 'igram-20150215_20150227' ... box : 4-tuple of int area to read, defined in (x0, y0, x1, y1) in pixel coordinate x/ystep : int, number of pixels to pick/multilook for each output pixel Returns: data : 2D/3D array atr : dict, metadata cSsg|]}|ddqS)r^r)rg).0irxrxry sz"read_hdf5_file..cSsg|]}d|kr|qS)r^rx)r{r|rxrxryr}scsg|]}|kr|qSrxrx)r{r|) ds_2d_listrxryr}srcss&|]}|p|ddVqdS)NrC)isdigit)r{xrxrxry (sz!read_hdf5_file..csg|]}dd|qS)z{}-{}r)format)r{r) ds_3d_listrxryr})sr^csg|]}|dqS))replace)r{r)dsFamilyrxryr}-scSs"g|]}|dr|ddqS)r^r@N) startswith)r{rrxrxryr}.srcsg|]}|kr|qSrx)keys)r{r|)frxryr}3scs"g|]}dd|kr|qS)z/{}r)r)r{r|)r_rxryr}4sz%input dataset {} not found in file {}rAr@rBN)rdrTcSsg|]}|dddqS)r^r@)rg)r{r|rxrxryr}Tscsg|]}|r|qSrx)r)r{j)rrxryr}Us zreading slice /z...css|]}|dkVqdS)r@Nrx)r{r|rxrxryr{s)get_slice_listappendrerhallrgh5pyFilelen ValueErrorrndimrjrrzerosshapebool_indexsumrdwhererangesysstdoutwriteflushprintanysqueeze)rtr_r`rarbrc slice_listds_listr|Z inputDateListdsNamesZ dsNamesOlddsrwxsizeysizeZ slice_flag date_listdZ num_sliceindsZd2rx)r_rr~rrryrps      " ,rpc Cs>tjtj|d}t||d}|d}t|d}t|d} |sXdd| |f}|dd} |d d } t|d d } |d d} d}d}|dkr\ddddd}| | kr|| } |d dd}|dkr t d| }|rN|dkrNd}nB|dkrJ|rB|dkr,d}n|dkrkr|d d?} |d d7} d@| kr|d@dAkrd } n|dBkrntdC||dBkrt||||||dD}n"t||| f|| | | | ||||dE }d|kr6| |d<||fS)FaRead data from binary file, such as .unw, .cor, etc. Parameters: fname : str, path/name of binary file datasetName : str, dataset name for file with multiple bands of data e.g.: incidenceAngle, azimuthAngle, rangeCoord, azimuthCoord, ... box : 4-tuple of int area to read, defined in (x0, y0, x1, y1) in pixel coordinate x/ystep : int, number of pixels to pick/multilook for each output pixel Returns: data : 2D array in size of (length, width) in BYTE / int16 / float32 / complex64 / float64 etc. atr : dict, metadata of binary file r@)r_ PROCESSORr&r+rr!rL BYTE_ORDERz little-endianr$rRr%rZphase)iscerOrMrN)bytefloatdoublecfloat FILE_TYPE.r)unwcorrA)Zband1 intensity magnitude)Zslc) amplituderrr)band2rcomplexlos)razhead)ZincLocallocalrZband3rB)magZamp)ranglerealimag)cpxr)roipac)z.unwz.corz.hgtz.msk)z.int)z.amp)z.demz.wgs84rJ)z.flgz.bytr)z.trans)rZazimuthgammaz big-endian)z.unwz.corz.hgt_simz.demz.ampz.ramp)z .utm_to_rdcr[z.slc complex32)z.mlisnaprYz byte orderr])gdalZgmtsarZhyp3ZcosicorrzUnknown InSAR processor: {})r`bandcpx_bandrarb) rr`r byte_ordernum_bandrrrrarb)rkrlrmrnrorirjgetupperrrminrrrrr read_gdal read_binary)rtr_r`rarbrvru processorrrr rrrrr dataTypeDictkrrwrxrxryrqs                                       rqFc stjtj|\}}|}x|dkr.NrF)rc)r)r)r)r)rr&r+cs|t|tjrx|jddfkrx|jdkr8tn@|jdkrdtfddt|jdD7ant d |jdS)NrArBcsg|]}d|dqS)z{}-{}r@)r)r{r|)namerxryr}sz?get_slice_list..get_hdf5_2d_dataset..rz#file has un-defined {}D dataset: {}) rerrrrrrrwarningswarnr)robj)rr)rryget_hdf5_2d_dataset~s "   "z+get_slice_list..get_hdf5_2d_datasetr$rR)z.transz .utm_to_rdc rangeCoord azimuthCoordrrAincidenceAngle azimuthAngle)z.unwrr)z.intz.slcroff)z.bip azimuthOffset rangeOffsetzcov.biprBZazimuthOffsetVarZrangeOffsetVarZ offsetCovarcSsg|]}d|dqS)zband{}r@)r)r{r|rxrxryr}s)rkrlrmrnrorirrrropen sliceListrrrrrrrj visititemsrrendswithr) rtZ no_complexfbaservrurd1_listrrrrx)rrrryrKsr                        " rc s|r |gStj|d}|dkr~t|}t|dt|dfdd}gat|d}| |WdQRXnt |atS) zCGet list of 2D and 3D dataset to facilitate systematic file readingr@)z.h5z.he5r&r+cs0t|tjr,|jddfkr,t|dS)Nr)rerrrrr)rr)rrrxryget_hdf5_datasets"z*get_dataset_list..get_hdf5_datasetrN) rkrlrmrorirjrrrrr)rtr_rvrurrrx)rrryget_dataset_listsrc CsVtj|d}|dkr dSd}t|d}t|d}||j}WdQRX|S)z+Get the compression type of input HDF5 filer@)z.h5z.he5Nrr)rkrlrmrorrr compression)rtextrds_namerrxrxryget_hdf5_compressions rc s tjtj|\}}|}tj|sNd|}|dt7}t||dkr|rj| dd}t |d<t j }fddDfd dDWd QRXd d d g}tfdddDrd}ntfdddDrd}ntfdddDrd}ntfdddDr:d}ndkrJd}ndkrZd}ntfdddDrxd }ntfd!d|Drttt|@d}nPd"|kr|d"}n|D]2\} } y| d)|| <Wn| || <YnXqWt||d"<d } t |d`|r\|kr\|} n$d*d+}ga|trtd} | d k rt| j |d,<Wd QRXd-|kr|d-|d.<d.|krd/|d.<nD|d0kri}d1|d2<d3|d$<d4|d5<d6|d7<d8|d9<d8|d:<t!|d d;t!|d;d }}|d<d=|d><|d?|d@<dA|d.<||d"<dB|d,<dC|dD<dE|dF<dG|dH<tj"||d2|d$dIkrdJ}|dK|d2|d$|d,7}t|nH|dL|dM|dNtj|ddO|dP|dQg}dRd|D}r,fdSd|D}t|dkrZ|t#krL|g}nt$dT|i}|%dUrtdVd|DrdW|d.<ntdXd|DrdY|d.<dZd|D}t|dkr0|t&|dn\td[d|Drd\|d.<n>td]d|Drd.|kr0d^|d.<n|t#kr0d_|d.<d.|krFd/|d.<|d}tj|d`}x |dakr~tj|\}}q`W|s|}|dLkr|t'|d"|kr||d"<n|dMkr|t&|d"|kr||d"<n|dNkr|t(|||d"<n|dOkr|t)||d.dWkrtj|}|*dbrZdc|d"<nL|*ddrpde|d"<n6|*dfrdg|d"<n dh|krdh|d"<n |di|d"<n||d"<n&|dPgt#kr|t+|||d"<djdBdkdldm}|,d,dn}|dnkr||kr|||d,<|r,|-dodp}|d"-dqdp}|dkrr|rh|t.krht.||dH<ndr|dH<nT|r|t.krt.||dH<n2dH|kr|t.krt.||dH<nds|dH<dt|krdu|kr|dt|du<tj/||dt<t0|}|S)vaRead attributes of input file into a dictionary Parameters: fname : str, path/name of data file datasetName : str, name of dataset of interest, for file with multiple datasets e.g. unwrapPhase in ifgramStack.h5 coherence in ifgramStack.h5 height in geometryRadar.h5 latitude in geometryRadar.h5 ... Returns: atr : dict, attributes dictionary zinput file not existed: {} zcurrent directory: )z.h5z.he5r^rrcs g|]}t|tjr|qSrx)rerGroup)r{r|)rrxryr}sz"read_attribute..cs.g|]&}t|tjr|jdkr|qS)rA)rerrr)r{r|)rrxryr}sNZinterferograms coherencewrappedc3s|]}|kVqdS)Nrx)r{r|)rrxryrsz!read_attribute..) unwrapPhaserrrc3s|]}|kVqdS)Nrx)r{r|)rrxryrs)heightlatituderrc3s|]}|kVqdS)Nrx)r{r|)rg1_listrxryrs)rrc3s|]}|kVqdS)Nrx)r{r|)rrxryrs)velocityrrreconsrc3s|]}|kVqdS)Nrx)r{r|)rrxryrs)igramfigramrc3s|]}|kVqdS)Nrx)r{r|)rrxryrsrzunrecognized file type: r+cSs0t|jdkr,d|jkr,tt|jdS)Nrr+)rattrsratr_listrdict)rrrxrxryget_hdf5_attrssz&read_attribute..get_hdf5_attrscSsg|] }t|qSrx)r)r{r|rxrxryr})sz!No attribute WIDTH found in file:utf8cSs$t|tjr |jdkr t|dS)NrA)rerrrrr)rrrxrxryrBsz(read_attribute..get_hdf5_datasetr!ZINSAR_PROCESSORrmintpyz.dehmipr&i(#g2.Z r6g2.Z ?r5degreesY_UNITX_UNITrAg?g?r4gY@r3ZGSIrLZLATLONZ PROJECTIONZWGS84ZGEODETIC_DATUMmZUNITrCzAinput DEHM file size do NOT match with the pre-defined 10m DEHM: z{} * {} in {}!z.rscz.xmlz.parz.hdrz.vrtz.aux.xmlcSsg|]}tj|r|qSrx)rkrlisfile)r{r|rxrxryr}|scsg|]}|r|qSrx)r)r{r|) metafile_extrxryr}sz(No metadata file found for data file: {}z.imgcss|]}|dVqdS)z.hdrN)r)r{r|rxrxryrsrcss|]}|dVqdS))z.xmlz.hdrz.vrtN)r)r{r|rxrxryrsrcSsg|]}|dr|qS)z.xml)r)r{r|rxrxryr}scss|]}|dVqdS)z.parN)r)r{r|rxrxryrsrcss|]}|dVqdS)z.rscN)r)r{r|rxrxryrsrrr@)z.geoz.rdrz.fullz.wgs84z.grdrz.unw)Zcohrz.corZ phase_ifgz.intdemz file typerOrMrN)rrrrnoneZStdrrradianrRZ FILE_PATHZ OG_FILE_PATH)1rkrlrmrnrorrgetcwd Exceptionrgrrrrrrrfsetrrr get_metadatarupdaterrrrargmaxitemsdecoderhrrdrgetsizeGDAL_FILE_EXTSFileNotFoundErrorr read_isce_xmlread_roipac_rscread_gamma_par read_envi_hdrr read_gdal_vrtrrrabspathstandardize_metadata)rtr_rrrvmsgruZpy2_mintpy_stack_filesrrZnum_listkeyvaluerryyxxZ metafilesZ xml_filesZmetafileZmeta_extrr rx)rrrrryrisp                         " "                              ricsdkr ttx t|D]\}}||<q WfddDfddD}xB|D]:}|}fdd|D}t|dkrf|d|<qfWS)zQConvert metadata input ROI_PAC/MintPy format (for metadata with the same values).Ncsg|]}|kr|qSrx)r)r{r|) standardKeysrxryr}sz(standardize_metadata..csg|]}|kr|qSrx)r)r{r|)metaDictrxryr} scsg|]}|kr|qSrxrx)r{r|)in_keysrxryr}sr)standardMetadataKeysriterrrr)Z metaDictInr#rr Zstd_keysZstd_keyZ cand_keysrx)r%r$r#ryrs   r=c Csi}i}d}g}dd}d}tj|rHt|d} | }WdQRXnt|tr\|d}x|D]} | } dd| |d D} t | d ks| d r| d ri}d }n&|r| ds| dsd}| |qb| d} t | d dddd} tj | } tj| } | dkr(| || <|rb|| r@| || <qbd}| || || <qbW|rj| |t |dkrt||d<|S)aReads the template file into a python dictionary structure. Parameters: fname : str full path to the template file delimiter : str string to separate the key and value print_msg : bool print message or not Returns: template_dict : dict file content Examples: tmpl = read_template(KyushuT424F610_640AlosA.template) tmpl = read_template(R1_54014_ST5_L0_F898.000.pi, ':') FcSs.|d}|dkrdS|ddko,t|dkS)NrFrplotr@)rgr)Z attributetokensrxrxryis_plot_attribute.s z(read_template..is_plot_attributeNr cSsg|] }|qSrx)strip)r{r|rxrxryr}@sz!read_template..r@rA)%#!>Tr.r/rrplotAttributes)rkrlrr readlinesrerhrgr-rrrr expanduser expandvarsjsondumps)rt delimiterrcZ template_dictZplotAttributeDictZinsidePlotObjectr2r+rrlinecrr rxrxry read_templatesL                r; c Cst|d}|}WdQRXi}xR|D]J}dd|dd|dD}|d}|dd d }|||<q(Wt|}|S) aaRead ROI_PAC .rsc file into a python dict structure. Parameters: fname : str. File path of .rsc file. Returns: rscDict : dict Dictionary of keys and values in RSC file. Examples: from mintpy.utils import readfile atr = readfile.read_roipac_rsc('filt_101120_110220_c10.unw.rsc') rNcSsg|] }|qSrx)r-)r{r|rxrxryr}~sz#read_roipac_rsc.. r<r@rr,r)rr3r-rrgr) rtr8rrZrscDictr9r:rr rxrxryrms  " r:c Cst|d}||d}WdQRXi}xz|D]r}dd||dD}t|dksd|drjtq0|d}t|dd d d dd} | ||<q0Wt |}t |}|S) aRead GAMMA .par/.off file into a python dict structure. Parameters: fname : str. File path of .par, .off file. delimiter : str, optional String used to separate values. skiprows : int, optional Skip the first skiprows lines. Returns: parDict : dict Attributes dictionary rNcSsg|] }|qSrx)r-)r{r|rxrxryr}sz"read_gamma_par..r@rA)r.r/rr,rr/) rr3r-rgrrnextrhrattribute_gamma2roipacr) rtr8skiprowsrrZparDictr9r:rr rxrxryrs  ( rcst}x t|D]\}||<qWx2|D]&tfdddDr2||d<q2Wx2|D]&tfdddDrf||d<qfWd|krt|}tt||d <d |krt|}t||d <d }||krt||}t|||d <d|kr2||d<d|krt|}d|krhdksnd|krdkrnn d|d<nd|d<d|krt|}d|krdkrnn d|d<nd|d<|S)z2Convert Gamma metadata into ROI_PAC/MintPy format.c3s|]}|VqdS)N)r)r{r|)rrxryrsz)attribute_gamma2roipac..)rrrrr&c3s|]}|VqdS)N)r)r{r|)rrxryrs)rrrr+Zradar_frequencyr*r r"Zsar_to_earth_centerr#rr(r iihiZZ ascendingr'Z descendingZ azimuth_anglerz-1 ANTENNA_SIDErR)rr'rrrrrhSPEED_OF_LIGHT)Z par_dict_inZpar_dictr ReZkey2rx)rryr@sJ        6   r@c CsXt|}i}|jdrx0|dD]"}|d}|dj}|||<q*Wxt ddgddgD]\}}|d |}t |d dj}t |d dj} t |d krdt |d krd||d |<| |d|d |<d|d |<qdWnT|jdkrL|d} x,| dD]}|d}|j}|||<qWt |d|d<t|}|S)z1Read ISCE .xml file into a python dict structure.imagepropertyrr Z coordinate1Z coordinate2XYz./component[@name='{}']z./property[@name='delta']z!./property[@name='startingvalue']g?gHz>z{}_STEPg@z{}_FIRSTrz{}_UNITZ PAMDatasetz./Metadata[@domain='ENVI']ZMDIrr )ETparseZgetroottagrfindallrfindtextziprrabsENVI2NUMPY_DATATYPEr) rtrootZxmlDictchildrr Z coord_nameprefixZv_stepZv_firstmetarxrxryrs0       rcCst|dd}t|dd|d<t|dd|d<d |krd d |d d D}tt|d }tt|dd}t|dkrt|dkrd|d<d|d<t||d<t||d<tt|d|d|d<tt|d|d|d<t |}|S)z0Read ENVI .hdr file into a python dict structurer()r8z data typerUr!z byte orderrRrzmap infocSs$g|]}|ddddqS){r})rr-)r{r|rxrxryr}sz!read_envi_hdr..,rDrEgg?gHz>rrrr5r6rBg@r3rCr4) r;rSrENVI_BYTE_ORDERrrgrRrrhr)rtruZmap_infox_stepy_steprxrxryr s     rcCsyddlm}m}Wntk r0tdYnX|||j}i}|j|d<|j|d<|j|d<| dj }t ||d<| d  d d }t||d <|}|d}|d } t|d} t|d d} | |d<| |d<|| d|d<| | d|d<|j|d} | dd|d<| } | r\d| kr\| dd|d<d|d<d|d<nJt| dkrt| dkrd |d<d |d<|dd!kr|dd"8<t|}|S)#zRead GDAL .vrt file into a python dict structure using gdal Modified from $ISCE_HOME/applications/gdal2isce_xml.gdal2isce_xml() written by David Bekaert. r)rosrzCannot import gdal and osr!r+r&r$r@r!ZIMAGE_STRUCTUREr%r\rBrDgr5r6g@r3r4)ZwktZ AUTHORITYZEPSGZUTMzUTM zoneZUTM_ZONEZmetersrrg?gHz>rgf@gv@)osgeorr^ ImportErrorOpen GA_ReadOnly RasterXSize RasterYSizeZ RasterCount GetRasterBandDataTypeGDAL2ISCE_DATATYPEZ GetMetadatarENVI_BAND_INTERLEAVEZGetGeoTransformrRZSpatialReferenceZ GetProjectionZ GetAttrValueZGetNamergr-r)rtrr^rrur r transformx0Zy0r\r]ZsrsZsrs_namerxrxryr&sH        rc Cst|d}|}WdQRXi}xR|D]J}dd|dd|dD}|d}|dd d }|||<q(Wt|}t|}|S) zRead GMTSAR .prm file into a python dict structure. Parameters: fname : str. File path of .rsc file. Returns: prmDict : dict Dictionary of keys and values in the PRM file. rNcSsg|] }|qSrx)r-)r{r|rxrxryr}psz#read_gmtsar_prm..r=r<r@rr,r)rr3r-rrgattribute_gmtsar2roipacr) rtr8rrZprmDictr9r:rr rxrxryread_gmtsar_prmbs  " rmc Cs*t}x t|D]\}}|||<qWd}||kr^||}|dkrVd|d<nd|d<d}||krt||}t|d}t|d}t|d }|||||}||d <d }||krt||}t|d |d <t|dt|dd } | t| } t| ddd|d<|S)zXConvert GMTSAR metadata into ROI_PAC/MintPy format (for metadata with different values).ZlookdirRz-1rDrRZSC_velr r)rrZ rng_samp_rateg@rZSC_clock_startZ SC_clock_stopg8@gN@r ) rr'rrrrrErjrh) Z prm_dict_inZprm_dictrr ZvelrFr)rZ az_pixel_sizeZ dt_centerZt_centerrxrxryrl{s2            rllrc Cs|\} } |sdd| | f}|dkr`td|\} }t| dkrT| d} tt|d}d| |}|}|dkrtj|||d| |d d | |}||d|d| |d|d| |d|d f}n|d kr>tj|||d| |d d | |}||d|dt |d|d ||df}n~|d krtj|||d| |d| d d | }|| |d|d| |d|d|d|d f}n t d|| dd dr@| dr|j }n\| dr|j}nH| drt|}n0| dr(t|}n| dr6n t d|| | dkrt|d |d| }t|d|d| }|t| d d| t| d d| f}|d|d|f}|S)aRead binary file using np.fromfile. Parameters: fname : str, path/name of data file to read shape : tuple of 2 int in (length, width) box : tuple of 4 int in (x0, y0, x1, y1) data_type : str, data type of stored array, e.g.: bool_ int8, int16, int32 float16, float32, float64 complex64, complex128 byte_order : str, little/big-endian num_band : int, number of bands interleave : str, band interleav type, e.g.: BIP, BIL, BSQ band : int, band of interest, between 1 and num_band. cpx_band : str, e.g.: real, imag, imaginary phase, mag, magnitude cpx x/ystep : int, number of pixels to pick/multilook for each output pixel Returns: data : 2D np.array Examples: # ISCE files atr = read_attribute(fname) shape = (int(atr['LENGTH']), int(atr['WIDTH'])) data = read_binary('filt_fine.unw', shape, num_band=2, band=2) data = read_binary('filt_fine.cor', shape) data = read_binary('filt_fine.int', shape, data_type='complex64', cpx_band='phase') data = read_binary('burst_01.slc', shape, data_type='complex64', cpx_band='mag') data = read_binary('los.rdr', shape, num_band=2, band=1) r)bbigz big-endianzieee-bez (\d+|\D+)r@z>{}{}rZrB)rdcountr_rAr[rYzunrecognized band interleaving:r1rr:rrpharzunrecognized complex band:N)rerNrrjrrrrfromfilereshapearangerrrrrrabsolute)rtrr`r rrrrrrarbrrletterZdigitrwrrrxrxryrs`"    .  & $(       rc Csyddlm}Wntk r,tdYnX|||j}||}|sZdd|j|jf}||d|d|d|df} t |j } | dd d r| d r| j } nT| d r| j} nB| d rt| } n,| d rt| } n| drn td|||dkrt|d|d|} t|d|d|} | t|dd|t|dd|f} | d| d| f} | S)aRead binary data file using gdal. Parameters: fname : str, path/name of data file to read box : tuple of 4 int in (x0, y0, x1, y1) band : int, band of interest, between 1 and num_band. cpx_band : str, e.g.: real, imag, phase, mag, cpx x/ystep : int, number of pixels to pick/multilook for each output pixel Returns: data : 2D np.array r)rzCannot import gdal!r@rBrAr1rr:rrrtr)rrzunrecognized complex band:N)r`rrarbrcrfrdreZ ReadAsArrayrhrgrrrrrrrryrrj) rtr`rrrarbrrZbndrwr rrrxrxryrs:  (         rc Cst|}d|krd|d<tt|d}tt|d}|sLdd||g}d}|dkr\d}tj|||d d |d |d d |}||d |d |d|d f}||d |d ||d||d f} || |fS) aReads roi_pac data (RMG format, interleaved line by line). ROI_PAC file: .unw, .cor, .hgt, .trans, .msk RMG format (named after JPL radar pionner Richard M. Goldstein): made up of real*4 numbers in two arrays side-by-side. The two arrays often show the magnitude of the radar image and the phase, although not always (sometimes the phase is the correlation). The length and width of each array are given as lines in the metadata (.rsc) file. Thus the total width width of the binary file is (2*width) and length is (length), data are stored as: magnitude, magnitude, magnitude, ...,phase, phase, phase, ... magnitude, magnitude, magnitude, ...,phase, phase, phase, ... ...... box : 4-tuple defining the left, upper, right, and lower pixel coordinate. Example: a,p,r = read_float32('100102-100403.unw') a,p,r = read_float32('100102-100403.unw',(100,1200,500,1500)) r!Zrmgr+r&rf4)rprqz big-endianzieee-bez>f4rBrA)rdrsr@)rirrjrrrrvrw) rtr`rrurrr rwrrrxrxry read_float32@s$  $r|cCst|}d|krd|d<tt|d}tt|d}|sLdd||g}d}|dkr\d}tj|||d |d |d |}||d |d |d|d f}||fS) z@Read real float64/double data matrix, i.e. isce lat/lon.rdr r!rMr+r&rf8)rprqz big-endianzieee-bez>f8rB)rdrsr@rA)rirrjrrrrvrw)rtr`rrurrr rwrxrxryread_real_float64ms   r~c Cst|}d|krd|d<tt|d}tt|d}|sLdd||g}d}|dkr\d}tj|||d |d |d |}||d |d |d|d f}|d krt|} n8|dkrt|} n$|dkr|j } n|dkr|j } n|} | |fS)acRead complex float 32 data matrix, i.e. roi_pac int or slc data. old name: read_complex64() ROI_PAC file: .slc, .int, .amp Data is sotred as: real, imaginary, real, imaginary, ... real, imaginary, real, imaginary, ... ... Parameters: fname : str, input file name box : 4-tuple defining (left, upper, right, lower) pixel coordinate. byte_order : str, optional order of reading byte in the file band : str output format, default = phase phase, amplitude, real, imag, complex Returns: data : 2D np.array in complex float32 Example: amp, phase, atr = read_complex_float32('geo_070603-070721_0048_00018.int') data, atr = read_complex_float32('150707.slc', 1) r!rNr+r&rZc8)rprqz big-endianzieee-bez>c8rB)rdrsr@rArrrr) rirrjrrrrvrwrryrr) rtr`rrrurrr rwZdataOutrxrxryread_complex_float32s2    rcCst|}d|krd|d<tt|d}tt|d}|sLdd||g}d}|dkr\d}tj|||d |d |d |}||d |d |d|d f}||fS) aRead real float 32 data matrix, i.e. GAMMA .mli file Parameters: fname : str, path, filename to be read byte_order : str, optional, order of reading byte in the file Returns: data : 2D np.array, data matrix atr : dict, attribute dictionary Usage: data, atr = read_real_float32('20070603.mli') data, atr = read_real_float32('diff_filt_130118-130129_4rlks.unw') r!rLr+r&rr{)rprqz big-endianzieee-bez>f4rB)rdrsr@rA)rirrjrrrrvrw)rtr`rrurrr rwrxrxryread_real_float32s   rcCsTt|}d|krd|d<tt|d}tt|d}|sLdd||g}d}|dkr\d}tj|||d d |d |d d |}||d |d d |dd |d f}td t |d } || d |d |d |d |d} || |d |d |d |d} |r.| | |fSt | | } t | | } | | |fSd S)amRead complex int 16 data matrix, i.e. GAMMA SCOMPLEX file (.slc) Gamma file: .slc Inputs: file: complex data matrix (cpx_int16) box: 4-tuple defining the left, upper, right, and lower pixel coordinate. Example: data,rsc = read_complex_int16('100102.slc') data,rsc = read_complex_int16('100102.slc',(100,1200,500,1500)) r!rr+r&ri2)rprqz big-endianzieee-bez>i2rBrA)rdrsr@N) rirrjrrrrvrwflattenrxrhypotarctan2)rtr`rrrurrr rwZodd_idxrrrrrxrxryread_complex_int16s2   $    rcCst|}d|krd|d<tt|d}tt|d}|sLdd||g}d}|dkr\d}tj|||d |d |d |}||d |d |d|d f}||fS) Nr!rJr+r&rr)rprqz big-endianzieee-bez>i2rB)rdrsr@rA)rirrjrrrrvrw)rtr`rrurrr rwrxrxryread_real_int16s   rcCs|dr|ddd}n|d}t|dd}d|krLd|d<tt|d}tt|d}|s|dd||g}tj|tj|d |d  |d |}||d |d |d|d f}||fS) zuRead binary file with flags, 1-byte values with flags set in bits For ROI_PAC .flg, *_snap_connect.byt file. z_snap_connect.bytrz.unw.rscz.rscr!boolr+r&rB)rdrsr@rA) rrgrirrjrrrrvrrw)rtr`ZrscFilerurrrwrxrxry read_bools    rc Cstj|tddtdddf}g}x|D]}||q,Wtj|dd}|dddf}t|dddf}t|dddf}t|dddf}|||||fS) N)rr@)rdusecolsr)r@rArBrC)rr@rArB)rrloadtxtbytesastyperhrrs) rtyyyymmddZYYYYMMDDyrwdatesZnorthZeastuprxrxry read_GPS_USGS6s r)NNTr@r@N)NNr@r@T)NNr@r@)F)N)NN)N)r(T)r<)r>rB)r() NrLror@rZr@rr@r@)Nr@rr@r@)Nro)Nro)Nror)Nro)NroF)Nro)N)7rkrrurZdefusedxml.ElementTreeZ ElementTreerKrr6numpyrrmintpy.objectsrrrrrrrrEr&rhZGDAL2NUMPY_DATATYPEZNUMPY2GDAL_DATATYPErrSrir[rzrprqrrrrirr;rrr@rrrrmrlrrr|r~rrrrrrrxrxrxry s $     9 ;  K `  ,  W  @%< ) \ > -  :  -