U bb@sddlZddlZddlZddlZddlZddlZddlm Z ddl m Z ddl Z e jdeddddlmZddlmZmZmZmZmZmZddlmZddlmZmZmZm Z!m"Z#dd l$m%Z%dd l&m'Z'm(Z(d Z)d Z*d dZ+d7ddZ,ddZ-ddZ.ddZ/d8ddZ0ddZ1d9ddZ2ddZ3ggdfd d!Z4d"d#Z5d$d%Z6d&d'Z7d(d)Z8d*d+Z9d,d-Z:d:d/d0Z;Gd1d2d2Zd5kre=ej?d6ddS)<N)make_axes_locatableignore matplotlib)categorymodule)giantIfgramStackgeometryDatasetNamesifgramDatasetNames ifgramStacktimeseriesKeyNames timeseries)GPS) arg_groupptimereadfileutilsplot)multilook_data)subsetversionaexample: view.py velocity.h5 view.py velocity.h5 velocity --wrap --wrap-range -2 2 -c cmy --lalo-label view.py velocity.h5 --ref-yx 210 566 #change reference pixel for display view.py velocity.h5 --sub-lat 31.05 31.10 --sub-lon 130.05 130.10 #subset in lalo / yx view.py timeseries.h5 view.py timeseries.h5 -m no #do not use auto mask view.py timeseries.h5 --ref-date 20101120 #change reference date view.py timeseries.h5 --ex drop_date.txt #exclude dates to plot view.py timeseries.h5 '*2017*' '*2018*' #all acquisitions in 2017 and 2018 view.py timeseries.h5 20200616_20200908 #reconstruct interferogram on the fly view.py ifgramStack.h5 coherence view.py ifgramStack.h5 unwrapPhase- #unwrapPhase only in the presence of unwrapPhase_bridging view.py ifgramStack.h5 -n 6 #the 6th slice view.py ifgramStack.h5 20171010_20171115 #all data related with 20171010_20171115 view.py ifgramStack.h5 'coherence*20171010*' #all coherence related with 20171010 view.py ifgramStack.h5 unwrapPhase-20070927_20100217 --zero-mask --wrap #wrapped phase view.py ifgramStack.h5 unwrapPhase-20070927_20100217 --mask ifgramStack.h5 #mask using connected components # GPS (for one subplot in geo-coordinates only) view.py geo_velocity_msk.h5 velocity --show-gps --gps-label #show locations of available GPS view.py geo_velocity_msk.h5 velocity --show-gps --gps-comp enu2los --ref-gps GV01 view.py geo_timeseries_ERA5_ramp_demErr.h5 20180619 --ref-date 20141213 --show-gps --gps-comp enu2los --ref-gps GV01 # Save and Output view.py velocity.h5 --save view.py velocity.h5 --nodisplay view.py geo_velocity.h5 velocity --nowhitespace a\Plot Setting: plot.name = 'Yunjun et al., 2016, AGU, Fig 4f' plot.type = LOS_VELOCITY plot.startDate = plot.endDate = plot.displayUnit = cm/yr plot.displayMin = -2 plot.displayMax = 2 plot.colormap = jet plot.subset.lalo = 33.05:33.15, 131.15:131.27 plot.seed.lalo = 33.0651, 131.2076 c CsPtjdtjtd}|dd}|jdtdd|jdtd gd d |jd d ddtd gdd|jddddd|jddddd gdd|jddddd d|jd!d"d#td$|jd%d&dd'd|jd(d)d*d+d,d-d.d/hd0d1t |}t |}t |}t |}t |}t |}t |}t |}t |}t |}t |}|S)2NzPlot InSAR Product in 2D) descriptionformatter_classepilogz Input FilezFile/Dataset to displayfilezfile for display)typehelpdset*z8optional - dataset(s) to display (default: %(default)s).)rnargsdefaultrz-nz --dset-num dsetNumListNUMzMoptional - order number of date/dataset(s) to display (default: %(default)s).)destmetavarrrrrz --nosearch search_dset store_falsez#Disable glob search for input dset.)r"actionrz--exz --exclude exDsetListZDsetz3dates will not be displayed (default: %(default)s).)r"r#rrrz --show-keptz--show-kept-ifgramplot_drop_ifgramz@display kept interferograms only, without dropped interferogramsz--plot-settingdisp_setting_filez!Template file with plot setting. )r"rz --noverbose print_msgzz cmd_line_parse..zview.py Agg)*rK parse_argssysargvvarsospathisfileFileNotFoundErrorformatr'rread_date_listoutfiledisp_figsave_fig lalo_step lalo_label zero_maskrOdisp_whitespace disp_axis disp_title disp_cbar disp_tickflip_lrflip_ud auto_fliplistsetrread_attributerkeysprintr*vprintjoinr)%update_inps_with_display_setting_filepltswitch_backend)iargsrGinpskeyfnameZ geo_opt_namesopt_namerIrIrJcmd_line_parsesJ      rc Cstdd}g}|jD]V}|tjtj|dtj|g}dd|D}t|dkrh||dqd}q|dkrtj |j }t dd|D}||krd}ntd ||j |S) Nzupdate mode: ONskippiccSsg|]}tj|r|qSrI)rcrdre.0irIrIrJ s zrun_or_skip..rruncSsg|]}tj|qSrI)rcrdgetmtimerrIrIrJrsz4{} exist and are newer than input file: {} --> skip.) r|rircrdr}dirnamebasenamelenappendrrminrg)rflagZoutfilesrfnamestitorIrIrJ run_or_skips  $ rcCs0t|}|js&d|kr&|d|_|jsFd|krFt|d|_|jsfd|krft|d|_|jsd|kr|d|_|jsd|krdd|d dd  d d D|_|j sd|krd d|d dd  d d D|_ |j s,d|kr,dd|d dd  D|_ |S)z2Update inps using values from display setting filezplot.displayUnitzplot.displayMinzplot.displayMaxz plot.colormapzplot.subset.lalocSsg|] }t|qSrIfloatrnrIrIrJrsz9update_inps_with_display_setting_file..,r]rcSsg|] }t|qSrIrrrIrIrJrszplot.seed.lalocSsg|] }t|qSrIrrrIrIrJrszplot.referenceLalo) r read_template disp_unitrzdisp_minrdisp_maxcolormap subset_latreplacesplit subset_lonref_lalo)rZ disp_set_fileZ disp_set_dictrIrIrJr~s    (( r~csJt|}tt|\__|j_|j_ddj j f_ t dt j t dt jt dt |j t dt jt dtjdjdjd jdd krd _jrt d jdkrd _tj|jjdjd_jrt djjsV|jsVtdtdg_n*|jdjd_t djjrjr|jjddd|jjdddg_n0d|krt |dt |dg_nd_t!|tj"|j#j$j%jd\_#_$t&|jdj'st(fdddDsj)dkrrjd*dddksjdd krd!d"g_'nj)d#ksj$rj%_'j+sЈjrʈj,rd$_+nd"_+j-rtj.|jd\_/_0j1stj2|d%jtd&_1t d'j1j3s>d(j1j4g_3t5S))Nrzdata coverage in y/x: zsubset coverage in y/x: zdata coverage in lat/lon: zsubset coverage in lat/lon: zH------------------------------------------------------------------------rMrg@@Fz=area exceed 2000 pixels, turn off default DEM contour displayT) datasetNamer*z$input reference point in lat/lon: {}zbWARNING: --ref-lalo is NOT supported when 1) file is radar-coded AND 2) no lookup table file foundz1 --> ignore the --ref-lalo input and continue.z$input reference point in y /x : {}y coord_typexREF_LATREF_LON)rwrap wrap_ranger*r*c3s|]}|jkVqdSrX)rlowerrrrIrJ @sz1update_inps_with_file_metadata..) coherencez.corr -)rcmask?)z.intg? FILE_PATH) datasetNames inps_dictzfigure title: {}{})6ut coordinatersubset_input_dict2boxrbpix_boxgeo_boxcheck_box_within_data_coverage box_pixel2geowidthlengthdata_boxr|r9maxdisp_dem_contourrN multilook_num multilookppauto_colormap_namerrr*rrg lookup_filer{ geo2radarref_yxyx2lalorzr read_pts2inpscheck_disp_unit_and_wraprrrcheck_map_projectionvlimanyrr transparencydisp_dem_shadervauto_flip_directionrtru fig_titleauto_figure_titlerifig_extupdate_figure_setting)rmetadatacoordrIrrJupdate_inps_with_file_metadatas         rTcCsd|_|dd|_|jr|jdr|jdkr|js@|jrd}|jdrlt |_|rt |dnp|jd rd | kr|d }t ||_|rt |d ||jrtd nt d |jt dd|_|S)a_Check/initiate the map projection object via cartopy based on inps and metadata. Cartopy requires that: 1. file is geocoded AND 2. file coordinates are in the unit of degrees / meters AND 3. set to display in geo-coordinates Use cartopy (by initiating inps.map_proj_obj) ONLY IF: 1. show fancy lat/lon label via --lalo-label OR 2. show coastline via --coastline This function will update the following variables: inps.map_proj_obj # cartopy.crs.* object or None inps.coord_unit # degree or meter inps.fig_coord # geo or radar NY_UNITdegreesdegmetergeoz!initiate cartopy map projection: r PlateCarreerUTM_ZONEz UTM zone z1--lalo-label is NOT supported for projection: UTMz*WARNING: Un-recognized coordinate unit: {}z1 Switch to the native Y/X and continue to plotradar) map_proj_objgetr coord_unitr startswith fig_coordrm coastlineccrsrr{rzUTM ValueErrorrg)rrr*msgutm_zonerIrIrJrbs8       rcCs|jrJz,|jd|jd}|jd|jd}Wn YnX|jdkrt|t|jd|j}t|t|jd|j}t|jdkr|||f}tj|st |s||8}t d |jn"d ||}|d7}t |d|_qPt|jdkrPd|kr"|jdkrbnn skip re-referencing.rr Frboxr*z1input reference point {} is out of data coverage!)rrrrr*zApply math operation: {}r,r-r.gr/rr0gf@r1z un-recognized math operation: {}rdata range: {} {}display range: {} {})(rrrr:rshapenpma is_maskedisnanr|rgr{squeezerrreadrrralltilereshaperscale_data4disp_unit_and_rewraprrrr* disp_scalerr+r,r-pinanminnanmaxdlim'auto_adjust_colormap_lut_and_disp_limitcmap_lut)datarrref_yref_xref_valrref_boxrIrIrJupdate_data_with_plot_inpss    D  ",           rc sljr tnddas*tdgt|tjtrPtj jj j dj_j rvtj j jjjd\}tdtjj\ jrʈjdrʈjdkrtd jrtd j|jjd j rtd tj|jjd t|tdjrjrjrtjdjdd} |d|ddd\}}|jd|jdf}|8tdjd|d|ddj!d_"jdjdjdjdf} |j#jj$dj$d| ddjj%dd } jdr<jdjddkrd)k r>|j3j>d)d)dfj>d)d)dfj7j5d9d+d,td-d}|rt?@d:d;gdgd?d@gd:d;gg}|j3|d)d)df|d)d)dfdAj5d9d+d,|A| dd|B| ddBtCjDEtCjDF|dCdDtCjDGrVz4tHjIdEjddFd tHjIdGjddFd Wn2dHtCjDJ}|dI7}t|d)YnXnd) fdJd/} | |_:d)}jKrtL|}|jMjNjOjOtPjQdK}tR| |\}jSr|jTjUj*jVdLjWrtdM|XjYrtdN|Zj[s"|\dOtdPj]rD|j^_d%|j`_d%n|abg|cbg|| |fS)Qa?Plot one slice of matrix Parameters: ax : matplot.pyplot axes object data : 2D np.array, metadata : dictionary, attributes of data inps : Namespace, optional, input options for display Returns: ax : matplot.pyplot axes object inps : Namespace for input options im : matplotlib.image.AxesImage object cbar : matplotlib.colorbar.Colorbar object Example: import matplotlib.pyplot as plt import mintpy.utils.readfile as readfile import mintpy.view as pv data, atr = readfile.read('velocity.h5') fig = plt.figure() ax = fig.add_axes([0.1,0.1,0.8,0.8]) ax = pv.plot_slice(ax, data, atr)[0] plt.show() c_sdSrXrIrYrIrIrJr[r\zplot_slice..rvlist)rrr*zdisplay data in transparency: rrzplot in geo-coordinatez#draw coast line with resolution: {}) resolutionzplotting DEM background ...)axrdemrr*zplotting image ...)siteFrrrMrz.format_coordrzplotting in Y/X coordinate ...zplotting Data ...?)rr r!r"r$r%r')r4NNREF_YREF_Xblackiiiiiiirrrzinputs/geometryRadar.h5latituder longitudez4WARNING: no latitude / longitude found in file: {}, z'skip showing lat/lon in the status bar.csd||}tt|jd}tt|jd}d|krPkrnnd|krhkrnn||f}t|stj|r|d7}n|d|7}jrЈ||f}t|s|d|7}r|d||f||f7}|S)Nzx={:.1f}, y={:.1f}rrMr@rAz , h={:.0f} mz, lat={:.4f}, lon={:.4f}) rgr:rrintrrrrrN)rrrrCrDrErF)rr geom_filerlatslonsrKrLrIrJrMs 0     )pad axes_class)fontsizecolorzflip figure left and rightzflip figure up and downoffzturn off axis display)dr*r{r|rr isinstancerr9r ColormapExtr cmap_vlistrNread_demrrrrrrrrrg coastlinesplot_dem_backgroundrrdisp_gps gps_component ref_gps_siter Zget_stat_lat_lonrrdisp_ref_pixelimshowr animation disp_scalebarZscalebar draw_scalebar scalebar_padr,rmdraw_lalo_labelrlr-r.r tick_paramsrr ref_markerref_marker_sizepts_lalo pts_markerpts_marker_sizeplot_gpsrMrrzr:pts_yxrarrayset_xlimset_ylimrcrdr}rrerrrrrr append_axescbar_loc cbar_sizerAxes plot_colorbarrq set_titler font_colorrt invert_xaxisru invert_yaxisrpaxisrsxaxis set_visibleyaxis get_xaxis set_ticks get_yaxis)rrrr dem_metadataZ ref_site_lalorrref_datar"imSNWErMrrZ temp_testrurcbardividercaxrI) rrGrrrHrIrJrVrrWrXrKrLrJ plot_slices~           &(   $       $0    rcCst|j}d}|jds.|d|d7}|d|dtj|j7}d|krj|d|d7}t d tj t t j t |t|d |_t|d |_|d|_tjtj |jd |_tj|jd |_t d|j|jftj|jdd|_t|\}}||fS)Nzinput file is )z.h5z.he5z{} PROCESSORz {} file: {} FILE_TYPE DATA_TYPEz in {} formatz run {} in {}LENGTHWIDTHrrMzfile size in y/x: {}T) no_complex)rryrendswithrgrcrdabspathrzr|r__file__rversion_descriptionr:rrrsplitextZfileBaseZfileExtget_slice_list sliceListread_dataset_input)ratrrrIrIrJread_input_file_infos$     rcsdd|D}|rt|tr"|g}g}|rj|D]8dkrDddd|fddD7}q.n|fdd|D7}ttt|}|fdd|D7}ttt|}fd d|D}||fS) z0Get dataset(es) from input dataset / dataset_numcSsg|]}|qSrIrIrrrIrIrJr5sz(search_dataset_input..rz*{}*z.*cs g|]}t|dk r|qSrX)rematchre)dsrIrJrEscsg|]}|kr|qSrIrIrallListrIrJrIscsg|]}|qSrI)indexrrrIrJrKscsg|] }|qSrIrIrrrIrJrQs)r^r9rgrsortedrwrx)rZinListZ inNumListr$tempListZ outNumListZoutListrI)rrrJsearch_dataset_input2s"   rc stjdkstjdkrtjdkr:tdjjdkr\d_tdjnHjdkrtjdkrdjdkrjdd\}}|g_|_t j jjjd_njd krd g_njd krd g_ntjd krd_nbjdkr8t j }|j dd|j dddg_n&j _jdkr^ddjD_t j jjjd_t j jgj\__g}tj }js|dd krtdtj fddj D}fdd|D}ntd|dd_tttjtjt|_fddjD_tj_jrjtkrtd_t j jggjdd}|stdjd_n|_jdkr tdtjj tj td tjtd!tjnPtdtjj tj td"tjjtd#tjjjrjtkrtd$jjdkrd%}|d&j 7}t|tjj jdddd'}|fS)(zFCheck input / exclude / reference dataset input with file dataset listrzinput dataset: "{}"velocityFz#turning glob search OFF for {} filer rM_r unwrapPhaseHDFEOS displacementgiantTimeseriesreconsrrrgeometrycSsg|]}|ds|qSbperp)rrrIrIrJr}s z&read_dataset_input..rz&do not show the dropped interferogramscs"g|]}|ddkr|qS)rrMrr)date12_drop_listrIrJrscsg|]}j|qSrI)rrrrrIrJrszH--show-kept option does not apply to file type: {}, ignore and continue.Tcsg|]}j|qSrI)rrrrIrJrsNzVWARNING: input reference date {} is not included in input file! Ignore it and continuer znum of datasets in file {}: {}znum of datasets to exclude: {}znum of datasets to display: {}zdatasets to exclude ({}): {}zdatasets to display ({}): {}zinput reference date: {}zNo input dataset found!z available datasets: {}r)rrr r|rgrr$rref_daterrrropenr'Z exDsetNumListrryr(r get_drop_date12_listrrwrxdsetNumr r{rcrdr Exception) rdate1date2objZ drop_num_listrZdrop_slice_listrrrI)rrrJrVs &           &       rcstjdjd}tjdjd}jdkrjsFd_jsjrjdkrtjdjd}tjdjd}tj ||fj j d_njstj _t djdjdjsd_jtjd krjd7_qjdkrxjdkrxtfd d d Drx|d |g}jddjdg}tj||j\__ttjtjjt_t dtjt dtjt dtjt dtjjsd_jjdkrd_jrZtjjd_tjtjjd\_ _!j!"_!ntjj#_tjtjj#d_ jdjdjdjd||krƈj d7_ j$rj d7_ j%dj%ddtj&krj tj%dj%d7_ j'r:j dj'7_ j(rPj d7_ jdkrrdj j!g_n fddt)djdD_fddjD_S)zUpdate figure setting based on number of subplots/datasets 1) fig_size and font_size 2) for multi: figure/row/column number 3) for multi: output file name rrMrrr)ds_shaperrr*zfigure size : [{:.2f}, {:.2f}]gd@c3s|]}|jkVqdSrXrarrrIrJrsz(update_figure_setting..)z--nrowsz--ncols皙?gffffff?zdataset number: zrow number: zcolumn number: zfigure number: 2_sub_wrapZ_refZ_exrcs"g|]}djt|jqS)z{}_{}{})rg outfile_baser9rrjrrIrJrsz)update_figure_setting..csg|]}tjj|qSrI)rcrdr}outdir)rrirrIrJrs)*rrrr,fig_sizerrabsrauto_figure_sizerrr*default_figsize_multir|rgfig_num fig_row_num fig_col_numrauto_row_col_numrceilastyper:r9rircrdrrrrrrrrrrr rr'range)rrr data_shapeZ fig_size4plotrIrrJrs     "2 "  rc stj||tjdjdjtjdjdjftjd}tjdkrjdkrt dfdd t ||D}t |jjjj d }| d d d stjjf|d|dd<n*t dttjjf|d|dd<jdkrBjddrBjrBj\}}|||d|df} tjj|| ddd} t |jdD]8} || ddddfdk} || | f| | 8<qhnt dtj||j d} t ||D]n} tjjj| jjjdd d}||| |ddddf<| j| |dj| dddq| jrt djtjjjjjjdd d} || 8}tjdksjdkstfdddDsjdd d!gksjd"d t tjDkrd#}nd}|r t||\}n&tfd$dd%Dr,td&d_ j!rTt d'tj"#|dk|}j$dk rt d(t%j$|jdddf}tj"#|dk|}t&|t'|g_(|r tfd)dd*Dr jdd+rjr jdd,kr t)j*|d-dd.\_+_,|S)/zIRead multiple datasets for one figure into 3D matrix based on i_start/endrrMrrdtype)r rr rrzreading data as a 3D matrix ...csg|]}j|qSrI)rrrrIrJr"sz$read_data4figure..rrxstepystepr*rfloat32complexNz;input data is complex, calculate its amplitude and continuer rFrrz)reading data as a list of 2D matrices ...maxValuer*/rsuffixzconsider input reference date: )r Z inversionc3s|]}|jkVqdSrX)dsetFamilyList)rdrrIrJr_sz#read_data4figure..) horizontalverticalrmodelresidualcSsg|]}d|dqS)zband{}rMrgrrIrIrJrasTc3s|]}|jkVqdSrXrrrrIrJrjs)z-uz--unitzeWARNING: -u/--unit option is disabled for multi-subplots with different units! Ignore it and continuemasking pixels with zero value masking datac3s|]}|jkVqdSrXr)rargrrIrJrzs)z-vz--vlimz--wrapunwrapr ) num_multilookr*)-rzerosr:rrrrrrr|rdictr*rrrrrr file_ref_yxrr progressBarrupdatercloserrrrr{rrnr masked_wheremskrr r r rrrr)i_starti_endrrr dset_listrZrrrrrmaskprog_barrZsame_unit4all_subplotsrrIrrJread_data4figures "  *  rc Cs|jr&tj|d|j|j|j||jd|j}|dkrHt |t |g}|j dd|j dd|j dd|j ddf}|j ||j |d|dd|j|dd }|jrd \}} |jr|jd|jd}} n&d |krt|d t|d }} |r| r|j| ||j|jd ||dd||dd|jrV|j|jdkrr|g|gdd} | |_|jrd} |j t!ks|j"d#dr|j"|$dd%dd} zt&j'(| d)dd} Wn| } YnXn|j"|} |j*+| }t,|j-dkrRd| krR| %dd} |j.ddkrR|j/+| }|j|j}|dkrn| } nfd|krdkrnnd0|| } ns|?d,|S)-zPlot one subplot for one 3D array 1) Plot DEM, data and reference pixel 2) axes setting: tick, ticklabel, title, axis etc. N)rr dem_shade dem_contourdem_contour_seqrr*rrNrrrMr)rr r!r$r%r"r'rOrPrQr9rrcSs d||S)Nzx={:.1f}, y={:.1f}, v =r)rrrIrIrJrMsz)plot_subplot4figure..format_coordrrrz%Y%m%drr z{} {}rr z{})size)r*prop)r[crimsonr\bold fontweightPLATFORMSenZSENTINEL1_SENSORAZ SENTINEL1_IPFz : {}r])@rNrrcrrrr*rrr r rrhrrrgrrzr:rrorprwrxrsrrrrrrMrqrr rrrrdtdatetimestrptime isoformatrrrrr date12Listrg fig_title_inrr,add_inner_titledropDatasetListrdisp_title4sentinel1 mplColorsr~rtrrurrpr)rrrrrrr"rrrrMZ subplot_titledate_strZ title_strZ title_indZ num_subplotr kwargZ s1_sensorZs1_IPFrIrIrJplot_subplot4figures               rc CsPdt||j|d}tdt|tj||j|j|jddd\}}|j j || }|d|j|j}t |j||j|jg}t||||}t|jtrtj|j|j|jdj|_tdtj|||jd} t||D]} | |} t| ||| || d d d d f|d } |jrt|jst|j| || d d d dd} tj|j |j!| d"dd }|rt| #|| j$| d|j!| %dddq| &~t||t'|D]} |(|| qt)*|j+d|j,dgt)-|j+d|j,dgg|_+td|j,|j.|jr0td|j|j.dd}|jr|jstdtd|/|jdd|jd|||nf|||d}|jddkr|d}td |j0d!d"|1d#d$|dd%|g}t2|| |\}} n ||||j3rLtd&t4j56|j|d|j7|j8|j|dd'd|j7d(|j9sL|:d S))zPlot one figure with multiple subplots 1) create figure 2) read all data into 3D array 3) loop to plot each subplot using plot_subplot4figure() 4) common colorbar and save zFigure {} - {}rM(----------------------------------------T)numfigsizenrowsncolssharexshareyrz plotting ...rN)rrrgQ?rNr(r)rrYshrinkaspect orientationrUNITrrrrrrcSsH|jddddddd|js"|jr6|j|j|jdn|jrD|dS)Ng{Gz?g\(\?g?)r5r6r8r7wspacehspace)r*r))subplots_adjust fig_wid_space fig_hei_spacefig_tight_layout tight_layout)figrrIrIrJadjust_subplots_layout?s z+plot_figure..adjust_subplots_layoutz-Note: different color scale for EACH subplot!z0Adjust figsize for the colorbar of each subplot.rg?g @rz show colorbarg(\?)r6gGz?rg{Gzt?save figure to {} with dpi={}tight) bbox_inches transparentdpi);rgr9rir|rsubplotsrrrcanvasmanagerset_window_titleflattenrrrr^rrr_rr`rrr*rrrrrcolorbarrryrrr set_labelrrrrdelaxesrr dlim_allr r rset_size_inchesr+add_axesr}rkrcrdrfig_dpisavefigrjclf)rrrrr0axsrrrrridxrr data_unitr1Z cbar_lengthrrIrIrJ plot_figures     $     $rHc stttdd|jD|_tttdd|jD|_t|jdkrv|jddkrvtttdd|jD|_|j dkr|j rt d |j n2|j r|j dkrtj|j|j|j|j|j d |_ |jd k r|j dkrt|j|j |j d d |_d |_|jd krdkrtdtd}}tdtd}}d|krb|krnn6d|kr~|krnn||g|_td |j|jdkr|jd_n|jdkr|jd_g|_|jdkr<|jr.css|]}|ddVqdS)StdrN)rrrIrIrJrrsrMrr css|]}|ddVqdS)rrMNrrrIrIrJrtsz/multilook {0} by {0} with nearest interpolation) max_memoryr*Nr)methodrrPrQrrrz#consider reference pixel in y/x: {}rg$@dg4@Frcsg|]}d|qS)z{}-{}rr)rrIrJrsz*prepare4multi_subplots..z@mark interferograms with 'dropIfgram=False' in red colored titlec3s|]}||kVqdSrXrIr)rrrIrJrs)rrzreading DEM: {} ... heightr)rrr*rz;WARNING: DEM file has a different size from the data file. zNThis feature is only supported for single subplot, and not for multi-subplots.z --> Ignore it and continue.)0rrwrxrrrrrrrr*r{rgrrauto_multilook_numrrr maxMemoryrrrrrzr:r|rrprrqr rrrrNrryrrcrdrrprepare_dem_backgroundrrrr) rrrrrrrZdropDate12ListrrrI)rrrrJprepare4multi_subplotsjs8         rQFc CsRt|dd}t|t|\}}t||}tj|j|j|j d|j |j |j |j d\|_|_tj|j|j d|j |j d\}}|jr|tj|j|j|j ddd8}|jdkr>|j dddd r>d |kr>t|d t|d }}tj|j|j d|||d|dfddd}||d k|8<|jrXtj|d k|}|jdk rztj|jd k|}ntj|jtjd |_tj|r|j|j9_|jt|j 9_n|jt|9_t!|||\}}|dkrDdd|j"D}|j#dk r t$|j#dni} t%j&|d| d\} } |||| fS|||fSdS)aPrepare data from command line as input, for easy call plot_slice() externally Parameters: cmd - string, command to be run in terminal Returns: data - 2D np.ndarray, data to be plotted atr - dict, metadata inps - namespace, input argument for plot setup Example: subplot_kw = dict(projection=ccrs.PlateCarree()) fig, ax = plt.subplots(figsize=[4, 3], subplot_kw=subplot_kw) W, N, E, S = (-91.670, -0.255, -91.370, -0.515) # geo_box cmd = 'view.py geo_velocity.h5 velocity --mask geo_maskTempCoh.h5 --dem srtm1.dem --dem-nocontour ' cmd += f'--sub-lon {W} {E} --sub-lat {S} {N} -c jet -v -3 10 ' cmd += '--cbar-loc bottom --cbar-nbins 3 --cbar-ext both --cbar-size 5% ' cmd += '--lalo-step 0.2 --lalo-loc 1 0 1 0 --scalebar 0.3 0.80 0.05 --notitle' data, atr ,inps = view.prep_slice(cmd) ax, inps, im, cbar = view.plot_slice(ax, data, atr, inps) plt.show() rMNrrOrrr r!r*rFrrrrPrQrrTcSsg|] }|dqS)g@rIrrIrIrJrszprep_slice..r/Figurer r subplot_kw)'rrr|rrr read_maskrrOrr mask_vmin mask_vmaxr*rrrrrrrzr:rnrrronesrbool_rrrrrrrrrr7) cmdZauto_figrrrrrrr rVr0rrIrIrJ prep_slicesn          r]c@s*eZdZdZd ddZddZddZdS) viewerzClass for view.py Example: import matplotlib.pyplot as plt from mintpy.view import viewer cmd = 'view.py timeseries.h5' obj = viewer(cmd) obj.configure() obj.plot() NcCs$|r|dd}||_||_dS)NrM)rr\r)selfr\rrIrIrJ__init__s zviewer.__init__c Cst|j}t|\}|_t||j}d|_|jrH|jsHt|dkrHd|_|j D]\}}t |||qRt j |j|j|jd|j|j|j|jd\|_|_|jS)NrrrrR)rrrrrr update_moderjr__dict__itemssetattrrrWrrOrrrXrYr*r)r_rrvaluerIrIrJ configure"s"  zviewer.configurec Cs|jdkrtdtj|j|jd|jdd\}|_|jr\|tj|j|j|jddd8}|j dkr|jd dd drd |j krt |jd t |jd }}tj|j|jd|||d|dfddd}||d k|8<|jrtd tj|d k|}|jdk r:td tj|jd k|}ntj|jtjd|_t|j}|jdk rtd|jtj||jk|}n6|dk rt|std|tj||k|}tj|r|j|j9_|jt|j9_n|jt|9_t||j|\}}|jdk r sR       / @i 5\ #$iWnxnT Q