U bb@sddlZddlZddlZddlZddlZddlZddlZddlZddl Z ddl m Z m Z m Z ddlmZmZmZddlmZddlZ deddeddeddZdZd Zd d Zdd d ZdddZddZGdddZdddZe dkr eej!dddS)N)sensorcluster RAMP_LIST)readfile writefileutils) STEP_LISTa=Command line options for steps processing with names are chosen from the following list: {} {} {} In order to use either --start or --dostep, it is necessary that a previous run was done using one of the steps options to process at least through the step immediately preceding the starting step of the current run.  a=example: smallbaselineApp.py #run with default template 'smallbaselineApp.cfg' smallbaselineApp.py #run with default and custom templates smallbaselineApp.py -h / --help #help smallbaselineApp.py -H #print default template options smallbaselineApp.py -g #generate default template if it does not exist smallbaselineApp.py -g #generate/update default template based on custom template smallbaselineApp.py --plot #plot results without run # Run with --start/stop/dostep options smallbaselineApp.py GalapagosSenDT128.template --dostep velocity #run at step 'velocity' only smallbaselineApp.py GalapagosSenDT128.template --end load_data #end after step 'load_data' zreference: Yunjun, Z., H. Fattahi, and F. Amelung (2019), Small baseline InSAR time series analysis: Unwrapping error correction and noise reduction, Computers & Geosciences, 133, 104331, doi:10.1016/j.cageo.2019.104331. cCstjdtjtdtd}|jdddd|jdd d d d d |jddddd|jddddd|jddddd|jddddd|dt}|jddd td!d"d#|jd$d%d&d td'd(d#|jd)d*d d+d,|S)-Nz;Routine Time Series Analysis for Small Baseline InSAR Stack ) descriptionformatter_classepilogcustomTemplateFile?z[custom template with option settings. ignored if the default smallbaselineApp.cfg is input.)nargshelp--dirz --work-dirworkDir./z'work directory, (default: %(default)s).)destdefaultr-ggenerate_template store_truez:generate default template (if it does not exist) and exit.)ractionrz-Hprint_templatez)print the default template file and exit.-vz --versionzprint software version and exit)rr--plotplotz5plot results [only] without running smallbaselineApp.z#steps processing (start/end/dostep)z--start startStepSTEPrz:start processing at the named step (default: %(default)s).)rmetavarrrz--endz--stopendStepz7end processing at the named step (default: %(default)s)z--dostepdoStepz%run processing at the named step only)rr"r) argparseArgumentParserRawTextHelpFormatter REFERENCEEXAMPLE add_argumentadd_argument_group STEP_HELPr)parserstepr0M/home/centos/operations/rsmas_insar/sources/MintPy/mintpy/smallbaselineApp.py create_parser<s>       r2c Cst}|j|d}|r|n tjdd|_tjtjtj d}|j rt |d}| }W5QRXz6ddl m}ddlm}|}|||d d d Wntk rt|YnXtd|jrttjjtd|js6tjtj|s6|js6|ttd } | d 7} | d7} td| |jrtj|jsXt|jtj|j|_tj|jtj|krd|_|jdgkrd} tdnd} t|t | d|_!|S)zCommand line parser.)argsNdefaults/smallbaselineApp.cfgrr)Console)Syntaxcfgr)Zbackground_colorz$no template file found! It requires:z& 1) input a custom template file, ORzN 2) there is a default template 'smallbaselineApp.cfg' in current directory.zERROR: rTz*plot smallbaselineApp results without run.F) step_list plot_only)"r2 parse_argssysargvospathjoindirnamemintpy__file__ropenreadZ rich.consoler7Z rich.syntaxr8print ImportErrorexitversionversion_descriptionrisfilebasenamer print_usager* SystemExitFileNotFoundErrorabspathread_inps2run_stepsrrunSteps) iargsr.inps template_fileflinesr7r8Zconsolemsgr;r0r0r1cmd_line_parse[sR          rZFc CsFdD]<}t||}|r||krd|}|d|7}t|q|jrX|j|_|j|_||j}||j}||krd|j|j}t||||d}|js|rg}t|dkr:t|dkrt t j j n t t j j t dtjt dtt d tjt|t d ||dd t d |S) z(read/get run_steps from input arguments.)r r#r%zInput step not found: {}z Available steps: {}z2input start step "{}" is AFTER input end step "{}"r4rz --RUN-at-{}--zCurrent directory: {}z+Run routine processing with {} on steps: {}zRemaining steps: {}Nz2--------------------------------------------------)varsformat ValueErrorr%r r#indexrlenrGrCrJrKlogodatetimenowr?getcwdr@rMrD) rUr:r;keyvaluerYidx0idx1 run_stepsr0r0r1rRs6         rRcstjtjtjd}tj|d}tj|sNtd|t ||nNt |}t |t fdd|Drtdt ||t||S)zGet the latest version of default template file. If an obsolete file exists in the working directory, the existing option values are kept. r5zsmallbaselineApp.cfgz/copy default template file {} to work directorycsg|]}|kqSr0)keys).0rdcdictr0r1 sz8get_the_latest_default_template_file..zAobsolete default template detected, update to the latest version.)r?r@rArBrCrDrLrGr\shutilcopy2r read_templateanyriutupdate_template_file)work_dirZlfilecfileldictr0rkr1$get_the_latest_default_template_files     rwc@seZdZdZd;ddZddZddZd d Zd d Zd dZ ddZ ddZ ddZ ddZ ddZddZedd9d:ZdS)?TimeSeriesAnalysiszY Routine processing workflow for time series analysis of small baseline InSAR stacks NcCs*||_tj||_tjt|_dSN)rr?r@rQrrccwd)selfrrr0r0r1__init__szTimeSeriesAnalysis.__init__cCstd|_|jr4tjtj|jd|_td|jtj|jddt |jtd|jt |j|_ | dS)zThe starting point of the workflow. It runs everytime. It 1) grab project name if given 2) go to work directory 3) get and read template(s) options Nrz Project name:Texist_okzGo to work directory:) projectNamerr?r@splitextrMrGmakedirsrchdirrw templateFile_read_templater{r0r0r1rEs    zTimeSeriesAnalysis.openc Csd|_|jrtd|jt|j}ddddddddd}|D] \}}||krB||||<qBdD](}||krh||dd||<qhd |kr|d |d <d D]}||kr| |qt ||_td t |j |j|_ d D]}tj|j|}tj|dd|j|j fD]b}|r t jtj|tj||ddddkr t||tdtj|tj|q qtd|j t|j |_t |j|_|jddkrdD]4}|j|dkrd|jd<td|qqdS)Nzread custom template file:autoyesno)defryontruenofffalse) mintpy.derampmintpy.troposphericDelay.method-_ processorzmintpy.load.processor) SUBSET_XMIN SUBSET_YMINz6update default template based on input custom template)inputspicTr}F)out_filein_filecheck_readable print_msgrunz© {} to {:<8} directory for backup.zread default template file:mintpy.geocode)mintpy.save.hdfEos5mintpy.save.kmzz*Turn ON mintpy.geocode in order to run {}.)customTemplaterrGrrpitemsrilowerreplacepopdictrrrsrr?r@rArr run_or_skiprMrnror\templatecheck_template_auto_value)r{rlZstandardValuesrdrebackup_dirname backup_dirtfiler0r0r1rsf              z!TimeSeriesAnalysis._read_templatecCs|d|jg}|jr$||jg7}|jr8|d|jg7}tdd|tj|t |j t j |j dddd\}}}|jrtd t j|t j|t j|jt ||jt ||j|s|jd d |jd d td |dS)aHLoad InSAR stacks into HDF5 files in ./inputs folder. It 1) copy auxiliary files into work directory (for Unvi of Miami only) 2) load all interferograms stack files into mintpy/inputs directory. 3) check loading result 4) add custom metadata (optional, for HDF-EOS5 format only) --templatez --projectz load_data.py Trrz:updating {}, {} metadata based on custom template file: {}F) print_aux) normal_endz.step {}: NOT all required dataset found, exit.N)_copy_aux_filerrrrGrArC load_datamainr?rrrrcheck_loaded_datasetr\r@rM add_attributer plot_resultclose RuntimeError)r{ step_namerT load_complete stack_file geom_filer0r0r1 run_load_dataGs.         z TimeSeriesAnalysis.run_load_datacstdr|jrtjtd|jdddg}tjfdd|Ddd}|D]D}tjtj||d d d krPt ||j t d tj|qPdS) NZ SCRATCHDIRzPROCESS/unavco_attributes.txtzPROCESS/bl_list.txtzSLC/summary*slc.jpgcsg|]}tj|qSr0)r?r@rArjiproj_dirr0r1rmsz5TimeSeriesAnalysis._copy_aux_file..T)rQFrrrrzcopy {} to work directory)r?getenvrr@rArr get_file_listrrMrnrorrGr\)r{flistfnamer0rr1rxsz!TimeSeriesAnalysis._copy_aux_filecstjjdddd\}}tjjd}fdddD}zd d|Dd }Wnd }YnXtjjd }d t|krtd ||tj ||ddkrtj |d d\}}t|} dD]8} | | krtd | tj || dd } d || d k<qd |d<t j |||d|djg} tdd| tj| |djdg} t|} tdd| Dr| ddddd g7} n&td!d| Dr| dd"dd#d$g7} td%d| jd&rtj |||jgdd'dkrtj| ntd(d S))z>Modify network of interferograms before the network inversion.Frr4rzcoherenceSpatialAvg.txtcsg|]}tjj|dqS)z network.pdf)r?r@rArrrr0r1rmsz?TimeSeriesAnalysis.run_network_modification..)rcSsg|]}tj|r|qSr0r?r@rLrr0r0r1rms rN waterMask.h5 waterMaskz#generate {} from {} for conveniencyrrr) datasetName)latitude longitudez)set pixels with 0 in {} to 0 in waterMask FILE_TYPE)rmetadata-tz modify_network.pyr --nodisplaycss|]}d|kVqdSphaseNrrr0r0r1 sz>TimeSeriesAnalysis.run_network_modification..z-d coherencerz0.2z1.0css|]}d|kVqdSoffsetNrrr0r0r1rs offsetSNR020z plot_network.py mintpy.plotrz1mintpy.plot is turned OFF, skip plotting network.)rrrrr?r@rArget_dataset_listrGr\rrFrwriterrCmodify_networkrrqrZ plot_network)r{rrrZcoh_txtZnet_figwater_mask_fileZ water_maskatrZ ds_name_listds_namedsrTdsNamesr0rr1run_network_modificationsP       z+TimeSeriesAnalysis.run_network_modificationcCstj|jddd}t|}tj|jd}tj|jd}tj|jd}tdd|Dr|d d |d g}t d d |t j |tdd|Dr|ddd |d g}n"tdd|Dr|ddd |d g}t dd |t j |dS)z.Generate auxiliary files from ifgramStack fileFrr4maskConnComp.h5avgSpatialCoh.h5zavgSpatialSNR.h5css|]}d|kVqdSrrrr0r0r1rsz>TimeSeriesAnalysis.generate_ifgram_aux_file..z --nonzero-o--update generate_mask.pyrcss|]}d|kVqdSrrrr0r0r1rs --datasetrcss|]}d|kVqdSrrrr0r0r1rsr temporal_average.pyN)rrrrrrr?r@rArqrGrC generate_maskrtemporal_average)r{rr mask_filecoh_fileZsnr_filerTr0r0r1generate_ifgram_aux_files    z+TimeSeriesAnalysis.generate_ifgram_aux_filecCsz|tj|jdddd\}}}tj|jd}|d|jd|g}|dk rZ|d |g7}td d |t j |dS) zSelect reference point. It 1) generate mask file from common conn comp 2) generate average spatial coherence and its mask 3) add REF_X/Y and/or REF_LAT/LON attribute to stack file Frr4rr-cNz--lookupz reference_point.pyr) rrrrrr?r@rArrGrCreference_pointr)r{rrr lookup_filerrTr0r0r1run_reference_points  z&TimeSeriesAnalysis.run_reference_pointcCstj|jddd}tj|jd}|ddd|dg}td d |tj |tj|jd }|d |d ddg}tdd |tj |dS)a A quick overview on the interferogram stack for: 1) avgPhaseVelocity.h5: possible ground deformation through interferogram stacking 2) numTriNonzeroIntAmbiguity.h5: phase unwrapping errors through the integer ambiguity of phase closure Frr4avgPhaseVelocity.h5r unwrapPhaserrrrrz --water-maskz--actionZ calculate unwrap_error_phase_closure.pyN) rrrrr?r@rArGrCrrunwrap_error_phase_closure)r{rrZ pha_vel_filerTrr0r0r1run_quick_overviews  z%TimeSeriesAnalysis.run_quick_overviewcCs|jd}|stddStj|jddd}tj|jd}|d|jd g}|d |g}|d kr~td d |t j |n|dkrtdd |t j |nr|dkr|ddddg7}td d |t j ||ddddg7}tdd |t j |nt d|dS)zCorrect phase-unwrapping errorszmintpy.unwrapError.methodz)phase-unwrapping error correction is OFF.NFrr4rrrz --cc-maskZbridgingz unwrap_error_bridging.pyrZ phase_closurerzbridging+phase_closurez-irrunwrapPhase_bridging!unwrapPhase_bridging_phaseClosurezun-recognized method: {})rrGrrrrr?r@rArrCZunwrap_error_bridgingrrr]r\)r{rmethodrrZ iargs_bridgeZ iargs_closurer0r0r1run_unwrap_error_corrections:    z.TimeSeriesAnalysis.run_unwrap_error_correctioncCsJtj|jddd}|d|jdg}tdd|tj|| dS) zInvert network of interferograms for raw phase time-series. 1) network inversion --> timeseries.h5, temporalCoherence.h5, numInvIfgram.h5 2) temporalCoherence.h5 --> maskTempCoh.h5 Frr4rrz ifgram_inversion.pyrN) rrrrrrGrArCZifgram_inversionr generate_temporal_coherence_mask)r{rrrTr0r0r1run_network_inversion&s  z(TimeSeriesAnalysis.run_network_inversionc stjjddd}tjjd}tjjd}jd}|d|d|g}jd d kr|d t|kr||d |d d ddg7}t dd|dddD}t dd}tj ||dddkrd}nVt d ||t |t fdd|Drd}t d |nt d |t d ||dkrdtj|i|D]}j||<qBt|tt|dd k} t d! | tjd"} | | krd# t| } | d$7} | d%7} | d&7} t| d'S)(z4Generate reliable pixel mask from temporal coherenceFrtemporalCoherence.h5maskTempCoh.h5"mintpy.networkInversion.minTempCoh-mr"mintpy.networkInversion.shadowMaskT shadowMask--base--base-dataset --base-value1rrcSsg|]}d|qS)zmintpy.networkInversion.{})r\rr0r0r1rmJszGTimeSeriesAnalysis.generate_temporal_coherence_mask..)Z minTempCohr zupdate mode: ONskip)rrrrz?1) output file: {} already exists and newer than input file: {}c3s(|] }tj||dkVqdS)FalseN)strrgetrrr{r0r1rRszFTimeSeriesAnalysis.generate_temporal_coherence_mask..z82) NOT all key configuration parameters are the same: {}z42) all key configuration parameters are the same: {}zrun or skip: {}rgznumber of reliable pixels: {}z#mintpy.networkInversion.minNumPixelz,Not enough reliable pixels (minimum of {}). zTry the following: zT1) Check the reference pixel and make sure it's not in areas with unwrapping errors zG2) Check the network and make sure it's fully connected without subsetsN)rrrrr?r@rArrrrGrr\read_attributerqrCrrrnpsumrFfloatintr) r{r tcoh_filertcoh_minrTZ config_keysflagrd num_pixel min_num_pixelrYr0rr1r8sL       z3TimeSeriesAnalysis.generate_temporal_coherence_maskrcstj|}tj|d}tj|d}t|}dddddg}t|D]}|}|dkr|ddrd tj |d }nR|dkr|d }|rd tj |d }n$|dkrP|d }|d}|r|dkrd tj |d }nX|dkrd tj |d }n6|dkrTimeSeriesAnalysis.get_timeseries_filename..csg|]}|dqS)r+r0r-r/r0r1rmsreference_dater$velocitygeocodeY_FIRST geo/geo_{}hdfeos5)r?r@rQrArrrr startswithr\rr]rsortedlistsetrirBrM)rrtZfname0Zfname1rZphase_correction_stepsr.rmodelrYr/fnamesfdirfbaser0r/r1get_timeseries_filenameos              z*TimeSeriesAnalysis.get_timeseries_filenamec Cstj|jddd}||j|j|}|d}|d}||kr|||d|g}tdd|tj||d d krtj |n$t |}| d d }td |d S)zCorrect local oscillator drift (LOD). Automatically applied for Envisat data. Automatically skipped for all the other data. Frrr+r,rz local_oscilator_drift.pyrrrr%Nz6No local oscillator drift correction is needed for {}.)rrrrr?rrGrArrCZlocal_oscilator_driftrrrrr\) r{rrr<rrrTrsatr0r0r1%run_local_oscillator_drift_corrections   z8TimeSeriesAnalysis.run_local_oscillator_drift_correctioncCstj|jddd}||j|j|}|d}|d}||kr|d|d|dg}td d |tj||d d krd dlm }| |ntddS)z Correct solid Earth tides (SET).Frrr+r,rrrz solid_earth_tides.pyrrrr)solid_earth_tidesz No solid Earth tides correction.N) rrrrr?rrGrArrCrBr)r{rrr<rrrTrBr0r0r1 run_solid_earth_tides_corrections  z3TimeSeriesAnalysis.run_solid_earth_tides_correctionc Cstj|jddd}tj|jd}||j|j|}|d}|d}||kr|jd}|jd}|jd } |jd } d d } | d kr|jd} |jd} |d|d|d|d|d| d| g }t dt dd|tj ||ddkrt j |q| dkr^|jd}d|d|d|d|g}t dt d d|tj ||ddkrt j |q| d!krd|d"|d|d#| g}t d$t d%|d&|}tj |||gddkrtj|r| || |kr||d|d'g}t d(t d)|t d*d|t j |nt|d+krDd,d-l m}t d.d|| |nB|d/krvd,d0l m}t d1d|| |ntd2|d3nt d4d5S)6zCorrect tropospheric delays.Frrrr+r,z"mintpy.troposphericDelay.polyOrderr'z#mintpy.troposphericDelay.weatherDirrcSst|}|d|dfS)NLENGTHWIDTH)rr)rrr0r0r1get_dataset_sizes zNTimeSeriesAnalysis.run_tropospheric_delay_correction..get_dataset_sizer(zmintpy.troposphericDelay.looksz'mintpy.troposphericDelay.minCorrelationrz-pr r-lrz>tropospheric delay correction with height-correlation approachz tropo_phase_elevation.pyrrrr)z!mintpy.troposphericDelay.gacosDirz-frz1tropospheric delay correction with gacos approachz tropo_gacos.pyr*z--modelz-wzVAtmospheric correction using Weather Re-analysis dataset (PyAPS, Jolivet et al., 2011)zWeather Re-analysis dataset:z./inputs/{}.h5z--force,--------------------------------------------z'Use existed tropospheric delay file: {}z diff.py)ZERA5r) tropo_pyaps3z tropo_pyaps3.py)ZMERRAZNARR) tropo_pyapsz tropo_pyaps.pyzun-recognized dataset name: .z!No tropospheric delay correction.N)rrrrr?r@rAr?rupperrGrrCZtropo_phase_elevationrZ tropo_gacosr\rLdiffrIrJr])r{rrrr<rrZ poly_order tropo_modelZ weather_dirrrFZ tropo_lookZ tropo_min_corrTZ GACOS_dir tropo_filerIrJr0r0r1!run_tropospheric_delay_corrections|                  z4TimeSeriesAnalysis.run_tropospheric_delay_correctioncCs|jd}|jd}||j|j|}|d}|d}||krtd||d|d|d|d g}td d |tj|ntd d S)z5Estimate and remove phase ramp from each acquisition.zmintpy.deramp.maskFilerr+r,z,Remove for each acquisition a phase ramp: {}z-sr rrz remove_ramp.pyrzNo phase ramp removal.N) rr?rrGr\rArCZ remove_rampr)r{rrrr<rrrTr0r0r1run_phase_deramping?s  z&TimeSeriesAnalysis.run_phase_derampingcCstj|jddd}||j|j|}|d}|d}||kr|d|jd|dg}|jd rh|d |g7}td d |tj |ntd dS)zbstep - correct_topography Topographic residual (DEM error) correction (optional). Frrr+r,rrrz,mintpy.topographicResidual.pixelwiseGeometryrz dem_error.pyrz#No topographic residual correction.N) rrrrr?rrrGrArCZ dem_errorr)r{rrr<rrrTr0r0r1#run_topographic_residual_correctionQs  z6TimeSeriesAnalysis.run_topographic_residual_correctioncCsFd}tj|r:|d|jg}tdd|tj|ntddS)z-Noise evaluation based on the phase residual.ztimeseriesResidual.h5rz timeseries_rms.pyrz9No residual phase file found! Skip residual RMS analysis.N) r?r@rLrrGrArCtimeseries_rmsr)r{rZres_filerTr0r0r1run_residual_phase_rmsfs  z)TimeSeriesAnalysis.run_residual_phase_rmscCsj|jdr^d|jg}||j|j|d}|D]}||g7}q0tdd|tj|ntddS)z;Change reference date for all time-series files (optional).zmintpy.reference.daterr+z reference_date.pyrzNo reference date change.N) rrr?rrGrArCr1r)r{rrTin_filesrr0r0r1run_reference_daters   z%TimeSeriesAnalysis.run_reference_datec Cs"||j|j|d}tj|jd}|d|jd|dg}tdd|tj ||jd }tj|jd |}tj |rtjtj|d }d tj|d |}tj|j|}|d|jd|dg}t|} |d | d d| d| dg7}tdd|tj |dS)z7Estimate average velocity from displacement time-seriesr+ velocity.h5rrrz timeseries2velocity.pyrr'z inputs/{}.h5rz{}{}.h5z --ref-dateREF_DATEz--ref-yxREF_YREF_XN)r?rrr?r@rArrGrCZtimeseries2velocityrrLr\rLrrMrr) r{rts_filevel_filerTrNrOsuffixZtropo_vel_filerr0r0r1run_timeseries2velocitys"   z*TimeSeriesAnalysis.run_timeseries2velocityc Cs|jdrx||j|j|d}t|}d|krntj|jd}tj |ddt j |jddd d \}}|d d |d g}|d|d|j d|dg}t dd|tj|tj|dtj|}tj|d} tj|d} |jd} | d| d| g}|jddkrcreating Google Earth KMZ file for geocoded velocity file: ...rWr4zgeo/geo_velocity.h5z{}.kmzrrz save_kmz.pyrcSsg|]}tj|r|qSr0rrr0r0r1rms z.z./geo/{}z./pic/{}NFrrz,save velocity to Google Earth format is OFF.)rrGr?r@rArrrrir\rrMr_rrrrCZsave_kmzr)r{rr\rZkmz_filerTr>Z kmz_filesr0r0r1run_save2google_earths&    z(TimeSeriesAnalysis.run_save2google_earthc Cs|jddkrx||j|j|d}|jdk r@t||jtj|jd}tj|jd}tj|jd}tj |jdd d }d |krtj|jd }tj|jd }tj|jd}tj|jd tj |}t d|d|d|d|d|d|j g }t dd|t|}t|} td | } t| dkrH| dnd} tj| |||||gddkrtj|nt ddS)zUSave displacement time-series and its aux data in geo coordinate into HDF-EOS5 formatrTr+NrrrFrrr_zgeo/geo_temporalCoherence.h5zgeo/geo_avgSpatialCoh.h5zgeo/geo_maskTempCoh.h5r5rHz--tcz--ascr rrz save_hdfeos5.pyrz{}_*.he5rrrz+save time-series to HDF-EOS5 format is OFF.)rr?rrrrrr?r@rArr\rMrGrrrrget_unavco_mission_namerr_rrCZ save_hdfeos5r) r{rr[rZ scoh_filerrrTrZSATZ hdfeos5_filesZ hdfeos5_filer0r0r1run_save2hdfeos5sF   z#TimeSeriesAnalysis.run_save2hdfeos5cCsz|D]n}td||dkr,||q|dkr@||q|dkrT||q|dkrh||q|dkr|||q|dkr||q|dkr||q|d kr| |q|d kr| |q|d kr| |q|d kr| |q|d kr | |q|dkr ||q|dkr6||q|dkrL||q|dkrb||q|dkr||qdS)zrun the chosen steps.z5 ******************** step - {} ********************rrrZquick_overviewcorrect_unwrap_errorZinvert_networkr r!r"r#r$Z residual_RMSr1r2r3Z google_earthr6N)rGr\rrrrrrrArCrPrQrRrTrVr^rcrdrf)r{r0r.r0r0r1r sJ                       zTimeSeriesAnalysis.runTc std|jd}tt|jd}tt|jd}t|jd}tj|jdddd \}}t j |jd }t j |jd } t j |jd } rt j n|rt j |n|}|rt j |n|}|rt j |n|}| rt j | n| } d ddddddg} dd|d|ddgddddddgd dddddg|g|gdddddgddgd dd!dddgd"dd!gd#dgd$dgd%dgd&gd'dddddgd(dddddgd)g| d*g| t j | d+ddgt j | d,ddgt j | d-ddgt j | d.d/gt j | d0g| d/|d1dd!gd2dd!gg} g} | D]~}|d3|dd }}|sXq4d4|krt t|}t|d3kr|D]}| |g|qn|| kr4| |q4d5d6| D} td7d6D7fd8d6| D} d9t|d:d;dd6| D} t}d}|jd?}|r|jd@}tj||}tjt| d|dA\}}}dB|krdCkrnndDn dD|dC}t|tt||d}|r|dkrtdE|||dFfdGdH| Dn| D]}tj|qtdItdJ}|D]}t || qtdKtdL}|tdM7}|tdN7}|tt j | dN7}|D]$}t!|t j | t j "|qzt#t|dO\}}tdP||dQ}|rt|d S)Rz1Plot data files and save to figures in pic folderz= ******************** plot & save to pic ********************r'zmintpy.plot.maxMemoryzmintpy.compute.maxMemoryzmintpy.plot.dpiFrr4Nrr_rz--noaxisz-ucmz--wrapz --wrap-rangez-55rWz--demz--maskrrgrayrrrz unwrapPhase-z --zero-maskcmyz coherence-rzconnectComponent-zunwrapPhase_bridging-zunwrapPhase_phaseClosure-z"unwrapPhase_bridging_phaseClosure-rrrrztimeseries_*.h5rar`zgeo_avgSpatialCoh.h5zgeo_velocity.h5r2zgeo_timeseries*.h5z.h5znumInvIfgram.h5r*cSs g|]}tj|dr|qS)rrrjrTr0r0r1rmsz2TimeSeriesAnalysis.plot_result..cSsg|]}|dqS)rr0)rjxr0r0r1rmscs4g|],}|dks,|dkr|dkr|qS)rr4r0rm)stack_dset_listrr0r1rms    z--dpiz --noverboserrz--memorycsg|] }|qSr0r0rm) opt_commonr0r1rmszmintpy.compute.clusterzmintpy.compute.numWorker)rmaxParallelNumg@g@g?z¶llel processing using {} cores ...)Zn_jobsc3s|]}tjj|VqdSry)rCviewrrm)delayedr0r1rsz1TimeSeriesAnalysis.plot_result..z© *.txt files into ./pic directory.z*.txtz,move *.png/pdf/kmz files to ./pic directory.z*.pngz*.pdfz*.kmz<z'time used: {:02.0f} mins {:02.1f} secs.aExplore more info & visualization options with the following scripts: info.py #check HDF5 file structure and metadata view.py #2D map view tsview.py #1D point time-series (interactive) transect.py #1D profile (interactive) plot_coherence_matrix.py #plot coherence matrix for one pixel (interactive) plot_network.py #plot network configuration of the dataset plot_transection.py #plot 1D profile along a line of a 2D matrix (interactive) save_kmz.py #generate Google Earth KMZ file in raster image save_kmz_timeseries.py #generate Goodle Earth KMZ file in points for time-series (interactive) )$rGrrLabsrrrrrrr?r@rArelpathr8globr_appendrrrtimerZ DaskClusterZformat_num_workercheck_parallelminmaxr\rCrrrrnromoverMdivmod) r{rrNZmax_plot_memory max_memoryZfig_dpirrrZgeo_dirZpic_dirZopt4tsZ iargs_list0Z iargs_listrTrr3r< start_timeZ run_parallelZ cluster_type num_workers num_coresParallelZ plot_memoryZtfilesrZpfilesZpfilemsrYr0)rsrprorr1rDs +     ,  " zTimeSeriesAnalysis.plot_resultcCs<td|jt|j|r8d}|d7}|d7}t|dS)NzGo back to directory:z1 ################################################z. Normal end of smallbaselineApp processing!)rGrzr?r)r{rrYr0r0r1rs  zTimeSeriesAnalysis.close)NN)r)T)T) __name__ __module__ __qualname____doc__r|rErrrrrrrrrr staticmethodr?rArCrPrQrRrTrVr^rcrdrfrrrr0r0r0r1rxs< D16&7 [M &+: $rxcCst}t|}t|j|j}||j|jd|jsR|j drZt |jdkrZ| | t t|d\}}td||dS)Nr/rr4rtz'Time used: {:02.0f} mins {:02.1f} secs )ryrZrxrrrErrSrrr_rrr~rGr\)rTrrUapprrr0r0r1rsr__main__r4)N)F)N)"r?r=rwryrarnr&numpyrrCmintpy.objectsrrr mintpy.utilsrrrrrZmintpy.defaults.templaterZmintpy.workflowr\r-r*r)r2rZrRrwrxrrr>r0r0r0r1 sF    A -