Zfi+ddlZddlZddlmZddlmZedZdZdZdZ ddZ dd Z d Z dd Z ed k(re ejd dyy)N)get_template_content) arg_utilsvelocityureferences: Efron, B., and R. Tibshirani (1986), Bootstrap methods for standard errors, confidence intervals, and other measures of statistical accuracy, Statistical Science, 54-75, doi:10.1214/ss/1177013815. Fattahi, H., and F. Amelung (2015), InSAR bias and uncertainty due to the systematic and stochastic tropospheric delay, J. Geophy. Res. Solid Earth, 120(12), 8758-8773, doi:10.1002/2015JB012419. Hetland, E., Musé, P., Simons, M., Lin, Y., Agram, P., & DiCaprio, C. (2012). Multiscale InSAR time series (MInTS) analysis of surface deformation. Journal of Geophysical Research: Solid Earth, 117(B2), doi:10.1029/2011JB008731 aexample: timeseries2velocity.py timeseries_ERA5_demErr.h5 timeseries2velocity.py timeseries_ERA5_demErr_ramp.h5 -t KyushuAlosDT73.txt timeseries2velocity.py timeseries.h5 --start-date 20080201 --end-date 20100508 timeseries2velocity.py timeseries.h5 --ex exclude_date.txt # complex time functions timeseries2velocity.py timeseries_ERA5_demErr.h5 --poly 3 --period 1 0.5 --step 20170910 timeseries2velocity.py timeseries_ERA5_demErr.h5 --poly 1 --exp 20170910 90 timeseries2velocity.py timeseries_ERA5_demErr.h5 --poly 1 --log 20170910 60.4 timeseries2velocity.py timeseries_ERA5_demErr.h5 --poly 1 --log 20170910 60.4 200 --log 20171026 200.7 timeseries2velocity.py timeseries_ERA5_demErr.h5 --poly 1 --polyline 20190101 20200501 # uncertainty quantification of the estimated time functions timeseries2velocity.py timeseries_ERA5_demErr.h5 --uq residue timeseries2velocity.py timeseries_ERA5_demErr.h5 --uq covariance --ts-cov timeseriesCov.h5 timeseries2velocity.py timeseries_ERA5_demErr.h5 --uq bootstrap z-exclude_date.txt: 20040502 20060708 20090103 c d}tdztzdztz}tj dd}t j |||||}|jdd|jd d d d |jdddd |jddddt j|d}|jd}|jdddd |jddd d! |jd"d#d$d%gd&tz'|jd(d)}|jd*d+d,d-d.hd/d01|jd2d3d4d5 |jd6d7d8td9d:;t j|}|jdd?d@ddA|jdBdCdDdEdFGt j|}|S)HNz4Estimate velocity / time functions from time-series. .)synopsis descriptionepilog subparserstimeseries_filez.Time series file for time function estimation.)helpz --templatez-t template_fileztemplate file with options)destrz-oz--outputoutfilezoutput file namez--update update_mode store_truezEnable update mode, and skip estimation if: 1) output file already exists, readable and newer than input file 2) all configuration parameters are the same.)ractionrF)plotzDates of interestz-sz --start-date startDatez'start date for time function estimationz-ez --end-dateendDatez%end date for time function estimationz--exz --ex-date excludeDate+zxdate(s) not included in time function estimation, i.e.: --exclude 20040502 20060708 20090103 --exclude exclude_date.txt )rnargsdefaultrzUncertainty quantification (UQ)z+Estimating the time function parameters STDz--uqz --uncertaintyuncertaintyQuantificationVALresidue>r bootstrap covariancezUncertainty quantification method (default: %(default)s). residue - STD from time series fitting residue (Fattahi & Amelung, 2015) covariance - STD from time series covariance bootstrap - STD from bootstrap resampling (Efron & Tibshirani, 1986))rmetavarrchoicesrz--ts-covz --ts-cov-filetimeSeriesCovFilezB4D time-series (co)variance file for time function STD calculationz--bcz--bootstrap-countbootstrapCountiz>number of iterations for bootstrapping (default: %(default)s).)rtyperrz Residual filez4Save residual displacement time-series to HDF5 file.z --save-resz--save_residualsave_resz8Save the residual displacement time-series to HDF5 file.z --res-filez--residual-fileres_fileztimeseriesResidual.h5zJOutput file name for the residual time-series file (default: %(default)s).)rrr) REFERENCETEMPLATEEXAMPLE__name__splitrcreate_argument_parser add_argumentadd_reference_argumentadd_argument_group DROP_DATE_TXTintadd_timefunc_argumentadd_memory_argument)r r r nameparserdateuqresids W/home/exouser/operations/rsmas_insar/tools/MintPy/src/mintpy/cli/timeseries2velocity.py create_parserr<6sEEH   (4 /' 9F >># r "D  - - xXfQ[]F )0`a  dGcd jy?QR  |MN - -f5 AF  $ $%8 9Dd> DFd$//)::YZ[[ }}A-QRR4==)DM4==)DM <<  !1!1$2F2F!GHK 4 4[[.r2FF!&- Kr=cDtdtjj|zddlm}m}m}t|}|j|jddg}|j|}d}t|jDcgc]}||z|jvs|} }| D]0} ||| z} | s| dvr|j| || <(| d vr$|j| jd || <P| d vr<| j!d d jd D cgc] } t#| c} || <| d vr%| j!d d jd || <| dvrU| j!dd j!dd } | jd D cgc]} | jd c} || <| dvr| || <| dvs#t%| || <3d} | |jvr|| rt#|| |_|Scc}wcc} wcc} w)z.Read input template file into inps.excludeDatez!read options from template file: r)ptimer@rA[]) skip_charszmintpy.timeFunc.)rr)r,)rS;)rT)rUrV/|)rr$)rRr%zmintpy.compute.maxMemory)rQrXrYrbrMror@rAvars read_templatercheck_template_auto_valuelistkeysyyyymmddr-replacefloatr3 maxMemory) rrdror@reiDicttemplate key_prefixikey_listkeyvaluexs r;rPrPs -bgg.>.>}.M MN99 JE%%d&8&8c3Z%PH++H5H$J -QaA1PQHQ(C( .."^^E2c '"^^EKK,<=c  $05 c#0F0L0LS0QR1eAhRc  $"]]3s399#>c & c#.66s3?49KK4DEqaggclEc JJ"c 88 Zc /(2 %C hmmo(3-x}- K=RSFsH)H(H#Hcjt|}ddlm}m}|jr ||dk(ry||y)Nr) run_or_skiprun_timeseries2time_funcskip)rmmintpy.timeseries2velocityrrr)rcrdrrs r;mainrs5 % D  K-7T"r=__main__rF)N)rXsysmintpy.defaults.templaterrMrr*r)r+r2r<rmrPrr,argvr=r;rsm 9" +   & <~EP*\#& z!"r=