B Qa4@sddlZddlZddlZddlZddlZddlZddlZedZ e ej ddl m Z mZmZmZdZddZdddZd d Zdd dZddddddgdfddZdddZedkreejd ddS)Nnumpy)readfile writefileutils1 attributezexample: multilook.py velocity.h5 -r 15 -a 15 multilook.py srtm30m.dem -r 10 -a 10 -o srtm30m_300m.dem # Ignore / skip marginal pixels multilook.py ../../geom_reference/hgt.rdr.full -r 300 -a 100 --margin 58 58 58 58 -o hgt.rdr c Cstjdtjtd}|jdddd|jddd d td d d |jddddtd dd |jdddd|jdddtdddgdd|jddtd d!d"d"d"d"gd#d$|S)%Nz Multilook.) descriptionformatter_classepilogfile+zFile(s) to multilook)nargshelpz-rz--rangez-xlks_xzEnumber of multilooking in range /x direction (default: %(default)s).)desttypedefaultr z-az --azimuthz-ylks_yzEnumber of multilooking in azimuth/y direction (default: %(default)s).z-oz --outfilez7Output file name. Disabled when more than 1 input files)r z-mz--methodmethodaveragenearestz_downsampling method (default: %(default)s) e.g. nearest for geometry, average for observations)rrrchoicesr z--marginmargin)ZTOPZBOTTOMLEFTRIGHTrz?number of pixels on the margin to skip, (default: %(default)s).)rrr metavarrr )argparseArgumentParserRawTextHelpFormatterEXAMPLE add_argumentintstr)parserr%F/home/centos/operations/rsmas_insar/sources/MintPy/mintpy/multilook.py create_parser'sr'cCsbt}|j|d}t|j|_|jdkr<|jdkr nanmedian)datarrrZ method_listmsgrGr5ysizexsize new_shapeZ crop_datatempZ coarse_datar%r%r&multilook_datadsl          $       rYrrc Cst|}t|}t|}t|dt|d}} |d} td|d| |td|td|td||d d d d gk r|d |d | |d ||d f} td |n d d | |f} tj|d } |s:ttj tj |kr.tj|d dt |dt |d| }n tj |}t j|||| d}| dkrftj|||dt|} tdd| D}t}x`| D]V}tdj||tj |d| dkr,t|d}||}t|jd}WdQRXtt|d|d}tt||dd}t|d}n| d | d }tt| d | d ||}xt|D]v}| d |||}| d |||d }t|| d }| d || d |f}t| d | d |t|| d |t| d | d |t|| d |f}td||dkrFtj|||||dd d }n"tj|||dd!d }t|||}|jd krd |jd |d |d |d |d g}n|d |d |d |d g}| dkrtj ||||dd"n|||<qjWqWt!|"d krLtj|d dkrL|#d#d$$d$krLtd%|d#td&d$|d#<| dkrt!|"|d'<tj%||||d(tj&|d)stj&|d*rt'|||S)+aA Multilook input file Parameters: infile - str, path of input file to be multilooked. lks_y - int, number of looks in y / row direction. lks_x - int, number of looks in x / column direction. margin - list of 4 int, number of pixels to be skipped during multilooking. useful for offset product, where the marginal pixels are ignored during cross correlation matching. outfile - str, path of output file Returns: outfile - str, path of output file LENGTHWIDTH FILE_TYPEzmultilooking {} {} file: {} PROCESSORz,number of looks in y / azimuth direction: %dz,number of looks in x / range direction: %dzmultilook method: {}rrHrIrz@number of pixels to skip in top/bottom/left/right boundaries: {}_Zalks_Zrlks)box)z.h5z.he5)metadataref_filecSsg|] }t|qSr%)r-).0ir%r%r& sz"multilook_file..z"multilooking {d:<{w}} from {f} ...)dwfrCrNi@ zbox: {}rF) datasetNamer_xstepystep print_msg)rir_rl)rSriblockrl INTERLEAVEBILz=the input binary file has 2 bands with band interleave as: {}zIfor the output binary file, change the band interleave to BIL as default.BANDS)out_filer`raz.xmlz.aux.xml)(r"rread_attributer/rLospathsplitextgetcwddirnameabspathr#basenameattrupdate_attribute4multilookr layout_hdf5get_dataset_listmaxdicth5pyFiler6prodr5ceilrintr=minreadrYndimwrite_hdf5_blockr-keysgetupperwriteisfilewrite_isce_xml)infilerrr.rr max_memoryatrlengthwidthkr_extdsNamesmaxDigitdsDictdsNamergdsds_sizenum_steprow_steprcr0r1Zbox_iZbox_orSrmr%r%r&multilook_files    $ 2     "          $ rc CsBt|}x,|jD]"}t||j|j|j|j|jdqWtddS)N)rrr.rrzDone.) r2r rrrr.rrr/)r0r1rr%r%r&mainCs r__main__)N)rrrE)N)rssysrr9rrr6logging getLoggerZ np_loggersetLevelWARNING mintpy.utilsrrrr*rrzr r'r2rDrYrr__name__argvr%r%r%r& s&    Y