o Ce@s>ddlZddlZddlZddlZejdeddddlmZ ddl m Z ddl ZddlmZddlmZmZddlmZddlmZmZmZdd lmZdd lmZmZm Z m!Z"d d Z#d dZ$d+ddZ%ddZ&d,ddZ'ddZ(ddZ)ggdfddZ*ddZ+dd Z,d!d"Z-d#d$Z.d%d&Z/d'd(Z0Gd)d*d*Z1dS)-NignoreZ matplotlib)categorymodule)make_axes_locatable)subsetversion)multilook_data)TIMESERIES_KEY_NAMESgiantIfgramStack ifgramStack)GPS)plotptimereadfileutilsc Cstdd}|jr d}|Sg}|jD]+}|tjtj|dtj|g}dd|D}t|dkr;| |dqd}q|dkretj |j }t dd |D}||krZd}|St|d |j d |S) Nzupdate mode: ONskiprunZpiccSsg|] }tj|r|qS)ospathisfile.0irrD/home/exouser/operations/rsmas_insar/tools/MintPy/src/mintpy/view.py +zrun_or_skip..rcss|] }tj|VqdSN)rrgetmtimerrrr 3szrun_or_skip..z& exist and are newer than input file: z --> skip.) vprintdisp_figoutfilerrjoindirnamebasenamelenappendrfilemin)inpsflagZoutfilesfnamefnamesZtitorrr run_or_skips( $ r/csZt|}tt|\__|j_|j_ddj j f_ t dt j t dt jt dt |j t dt jt dtjdjdjd jdd krtd _jrtt d jdkr|d _tj|jjdjd_jrt djjs|jstdtdg_n|jdjddd _t djjrوjr|jdjddd _nd|vrt|dt|dg_nd_t |tj!|j"j#j$jd\_"_#t%|jdj&sOt'fdddDs;j(dkr3jd)dddvs;jddkrAddg_&nj(d vsKj#rOj$_&j*sbjr_j+r_d!_*nd_*j,rrtj-|jd\_._/j0stj1|d"jtd#_0t d$j0j2sj03d%d&}j4s|d'7}|j5g_2t6S)(Nrzdata coverage in y/x: zsubset coverage in y/x: zdata coverage in lat/lon: zsubset coverage in lat/lon: zH------------------------------------------------------------------------g@@Fz=area exceed 2000 pixels, turn off default DEM contour displayT) datasetName print_msgz"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 : REF_LATREF_LON) disp_unitwrap wrap_ranger4r4c3s|] }|jvVqdSr)keylowerrr*rrrsz1update_inps_with_file_metadata..) coherencez.corr -)r>cmask?)z.intg? FILE_PATH) datasetNamesZ inps_dictzfigure title:  Z_nws)7ut coordinaterZsubset_input_dict2boxvarspix_boxgeo_boxZcheck_box_within_data_coverageZ box_pixel2geowidthlengthZdata_boxr strmaxdisp_dem_contourdem_file multilook_num multilookppZauto_colormap_namecolormapdsetr4ref_laloZ lookup_fileprint geo2radarref_yxZ radar2geokeysfloatZ read_pts2inpsZcheck_disp_unit_and_wrapr7r8r9check_map_projectionvlimanyr;split transparencydisp_dem_shade auto_flipZauto_flip_directionflip_lrflip_ud fig_titleZauto_figure_titler"replacedisp_whitespacefig_extupdate_figure_setting)r*metadatacoordfbaserr=rupdate_inps_with_file_metadata=s    " $   & rnTcCsd|_|dd|_|jrl|jdrl|jdkrl|js |jrld}|jdr5t |_t |d|S|jd rld | vr]|d }t ||_t |d ||jr[td |Std |jtdd|_|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 yx NY_UNITdegreesdegmetergeoz!initiate cartopy map projection: rr PlateCarreersUTM_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 plotyx) map_proj_objgetr< coord_unitrK startswith fig_coord lalo_label coastlineccrsrur r[UTM ValueErrorrX)r*rkr4msgutm_zonerrrr]s6       r]cCs|jrO|jd|jd|jdd|jddg|_|jd|jd}|jd|jd}|jdkrPt|t|jd|j}t|t|jd|j}t|jdkrd|krd|jdkr|nnd|krs|jdkr|nn|||f}ntj|j |j d|jddd}t j |st |s||8}td|jnd |d |d }|d 7}t|d|_nt|jd krNd|kr|jdkrnnd|kr|jdkrnn t |dd||f}n|jdkrtj|j |j |jddd}n td|jdt j |s:t t |r:|t |dddd|jd|jdf8}td|jnd |d |d }|d 7}t|d|_nd|_tj|||j|j|j|jd\}|_|_|_|jrp|j|_|jrtd|j|jdkrt |}nF|jdkrt !|}n:|jdkr|d9}n/|jdkrd|}n$|jdkr|dt j"9}n|jdkr|t j"d9}ntd|jt #|t $|g|_%|jstj&||jd\|_'|_|_(td|j%d|jtd |jd|j||fS)!Nr1rr2Fr3boxr4zset reference pixel to: z WARNING: input reference pixel (, z!) has either masked or NaN value!z -> skip re-referencing.r0 timeserieszinput reference point z is out of data coverage!)rkr7r8r9r4zApply math operation: squaresqrtreverseginverserBrad2deggf@deg2radzun-recognized math operation: r:data range: rEdisplay range: ))rZZref_boxrJrRintr&shaperreadr(rVnpma is_maskedisnanr rXsqueezer;ralltilereshaperTZscale_data4disp_unit_and_rewrapr7r8r9r4Z disp_scaler^math_operationrrpinanminnanmaxdlim'auto_adjust_colormap_lut_and_disp_limitcmap_lut unique_values)datarkr*ref_yref_xref_valrrrrupdate_data_with_plot_inpss << ",rFc sdddlm}|dd}tdd|Drt|}g}tdd|Dd}td d|Dd}t|D]U\}tfd dt||Drt||D]7\}} |kr_|dd} qP|kri| krrnn| d |7} qP| kr| d |dd 7} | | nqPq9| |q9||} | j rt nd da t | \} } t| | } tj| j| j| jd| j| j| j| j d\| _| _tj| j| jd| j| j d\} } | jr| tj| j| j| jddd8} | jdvr-| jddddr-d| vr-t| dt| d}}tj| j| jd|||d|dfddd}| | dk|8<| jr:tj !| dk| } | jdurKtj !| jdk| } n tj"| j#tj$d| _tj %| rr| j| j&9_| jt'| j(9_n | jt'| 9_t)| | | \} } |rdd| j*D}| j+durt,| j+dni}t-j.||dd}| | | |fS| | | fS)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: from cartopy import crs as ccrs from mintpy.view import prep_slice, plot_slice # initiate matplotlib figure/axes subplot_kw = dict(projection=ccrs.PlateCarree()) fig, ax = plt.subplots(figsize=[4, 3], subplot_kw=subplot_kw) # compose view.py command cmd = 'view.py geo_velocity.h5 velocity --mask geo_maskTempCoh.h5 --dem srtm1.dem --dem-nocontour ' cmd += f'--sub-lon -91.7 -91.4 --sub-lat -0.5 -0.3 -c jet -v -3 10 ' cmd += '--cbar-loc bottom --cbar-nbins 3 --cbar-ext both --cbar-label "LOS velocity [cm/year]" ' cmd += '--lalo-step 0.2 --lalo-loc 1 0 1 0 --scalebar 0.3 0.80 0.05 --notitle' # call prep/plot_slice() data, atr, inps = prep_slice(cmd) ax, inps, im, cbar = plot_slice(ax, data, atr, inps) plt.show() r)cmd_line_parser1Ncs|]}|dVqdS)"'Nr{rxrrrrizprep_slice..cSg|]}|dqSrrrrrrrozprep_slice..cSrrendswithrrrrrprc3s,|]\}}|ko|knVqdSrr)rind0ind1rrrrrs*rErc_dSrrargskwargsrrrzprep_slice.. mask_filer3rvminvmaxr4rFr r? unwrapPhaseREF_YREF_XrAdtypecSsg|]}|dqS)g@rrrrrr projectionfigsize subplot_kw)/Zmintpy.cli.viewrr`r_listrwhere enumeratezipr'r4rXr read_input_file_infornrT read_maskr(rrVrJ mask_vmin mask_vmaxmskrrref_dater;r{r[r zero_maskr masked_whereonesrbool_rmaskrrrfig_sizerxdictpltsubplots)cmdZauto_figriargsZ temp_iargsZind0sZind1sZ temp_iargrrtempr*atrrrrref_datarraxrrr prep_sliceHs           rc s jrtnddat jtrtj j j j dj _ j r/tj j j j jd\}tdt jj\ tj| jd j r jdr jdkrtd  j d  j d  j d  j d f _ j d  j d  j d  j d f} jrtd j|j j jd j rtdtj| j  jdt| jr jr jrt jdjdd} |d |d d d \}}| j d | j d f} | 8td jd|d dd|d dd| dd j! d _" j#rtj$| jd} ntd|j% j j&d  j&d  jd j' j j(d d } j)rEtj*| j)| j+ j, jd jd r\ j d  j d d!kr\d _- j-rztd" j.tj/| j j j. j0 j1 j2d# j3rtj4| j j5 j6 j7 j8 j9r j9n j1 j: jd$ n |j;d%d& j1d'd'd'd'd( j"rƈ j<r|j= j j?d)td* j@d+ur|j= j@d+d+d f j@d+d+d f jA jBd,d-d.td/ jrtjC|| | jd} j rt|j\ f d0d1} | |_Dnxd2 _td3 j r7tdtj|d+ jd j d d4 j d d4 j d d4 j d d4f _ j#rbtj$| jd} ntd|j% j j&d  j&d  j j' jd d5} |j; j1d6 j"rd7\} } jEr jEd  jEd } } nd8|FvrtG|d8tG|d9} } | r| r|j=| | j> j?d)td* jHd+ur|j= jHd+d+d f jHd+d+d f jA j?d:d-d.td/d}|rtIJd;dgd=d?gd@dAgd;d.rZvlist)rJrKr4zdisplay data in transparency: )rzrqrtzplot in geo-coordinaterr2r0r1z!draw coast line with resolution: ) resolution linewidthzplotting DEM background ...)rrKdemr*r4)siteFr:z=referencing InSAR data to the pixel nearest to GNSS station: z at [z.6frz] by substrating .3frEzplotting data ...upper) cmaprrextentorigin interpolationalphaZanimatedzorder)rfaultline_fileSNWErmin_distr4rrzplot scale bar: )rrKunitlocZlabelpad font_sizer) rrK lalo_steplalo_loc lalo_max_num lalo_offsetrrr4bothinT)whichZ direction labelsizeleftrighttopbottommszplot reference pointNkrB)r ZmecZmewzplot points of interestc sd|ddd|dd}j|ddjd}j|ddjd }d|kr4 krn|Sd|krA krn|S||f}|t|sVtj|rXd nd |d 7}jrj|ddd}j|ddd }d|krkrnn!d|krkrnn||f}|t|rd nd|d7}|d|jdd7}|d|jd d7}|S)NzE=.f, N=lon)Z coord_typerlatr1, v=[], v=r, h=[], h=.1fz, x=z.0f, y=)Zlalo2yxrJrrrrrQ) ryrcolrowvZdem_colZdem_rowh) rl coord_demrrdem_len dem_pix_boxdem_widr* lalo_digitnum_colnum_rowrr format_coordfs$$  *0 z plot_slice..format_coordrwzplotting in Y/X coordinate ...?)rrrrrrr)rNNrrZblackiiiiiiir?rCzinputs/geometryRadar.h5c3s|]}|vVqdSrrr) geom_ds_listrrrszplot_slice..)latitude longituder*rr+z0WARNING: no latitude / longitude found in file: z'skip showing lat/lon in the status bar.cs(d|dd|d}tt|jd}tt|jd}d|kr,krn|Sd|kr9krn|S||f}|t|sNtj|rPdnd|d7}jrp||f}|t|rhd nd |dd 7}r|d ||fd d7}|d||fd d7}|S)Nx=rrrr1rrrrrz mz, E=rrr)rrrintrJrrrrQ)rrrrrrr) rr geom_filer*r"latslonsr#r$rrr%s   *    )padZ axes_class)fontsizecolorzflip figure left and rightzflip figure up and downoffzturn off axis displayr )r r Z labelbottomZlabeltop)ZrotationZcentervazrotate Y-axis tick labels by z deg)tr4rXr isinstancerUrNrT ColormapExtr cmap_vlistrQZread_demrJrKrarrGget_lalo_digit4displayrzr{r|rr~Z coastlinescoastline_linewidthplot_dem_backgroundrHdisp_gps gps_component ref_gps_siter Zget_stat_lat_lonrYr7disp_ref_pixeldisp_dem_blendZplot_blend_imageimshowr^r animationrZplot_faultlinefaultline_linewidthfaultline_min_dist disp_scalebarscalebarZ draw_scalebar scalebar_padrscalebar_linewidthr}Zdraw_lalo_labelrrrrlalo_font_sizerxZ tick_paramsrWr ref_markerref_marker_sizepts_lalo pts_markerpts_marker_sizeZplot_gpsr%rZr[rpts_yxrarrayset_xlimset_ylimrrr#r$rrget_dataset_listrrr% disp_cbarrZ append_axescbar_loc cbar_sizerZAxes plot_colorbar disp_title set_titlerf font_colorrd invert_xaxisre invert_yaxis disp_axisaxis disp_tickZxaxisZ set_visibleZyaxis get_xaxis set_ticks get_yaxis ylabel_rotrZsetpZget_yticklabels)rrrkr* dem_metadatarZ ref_site_lalorrrimr%rrZ temp_testrQrcbarZdividercaxrr)rlrrrrr r!r)r.r*r"r/r0r#r$r plot_slices     &$     (    $       $0      rjcCs"t|j}d}|jds|d|d7}|d|dtj|j7}d|vr5|d|d7}t d tj t d 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 {} formatzrun z in LENGTHWIDTHrr1zfile size in y/x: T) no_complex)rread_attributer(rformatrrabspathr[r r%__file__rversion_descriptionrrMrLr;splitextZfileBaseZfileExtget_slice_list sliceListread_dataset_input)r*rrrrrrs$     rcsdd|D}|rTt|tr|g}g}|r6|D]dvr#dddd|fddD7}qn |fdd|D7}ttt|}|fdd|D7}ttt|}fdd|D}||fS) z0Get dataset(es) from input dataset / dataset_numcSsg|]}|qSrrrrrrr>sz(search_dataset_input..*z.*cs g|] }t|dur|qSr)rematchre)dsrrrNs csg|]}|vr|qSrrrall_listrrrQcsg|]}|qSr)indexr}rrrrSrcsg|]}|qSrrrrrrrYr)r8rNrgsortedrset)rin_list in_num_list search_dsetZtempListZ outNumListZoutListr)rrrsearch_dataset_input;s$    rc stjdkstjdkrdtjdkrtdjdjdkr0d_tdjdn$jdkrTtjd krTd jdvrTjdd \}}|g_|_tj jjjd d _n[jd krnd g_nBjdkrxdg_n8jdkrd_n/jdkrt j }|j dd|j dddg_nj _jdkrddjD_tj jjjd d _tj j gjd \_ _g}tj d}js |d krtdtj fddj D}fdd|D}n td|dd_tttjtjt|_fddjD_tj_jrdjtvr;d _tj jggjd dd}|sad!jd"}|d#7}t|d _n|_jd$vrtd%tjj d&tj td'tj td(tjn.td%tjj d&tj td)tj d*j td+tjd*jjrшjtvrtd,jjdkrd-}|d.j 7}t|tjj jdddd/} | fS)0zFCheck input / exclude / reference dataset input with file dataset listrzinput dataset: "rvelocityFzturning glob search OFF for z filerr1_)rrrrr rHDFEOS displacementgiantTimeseriesreconsr r:r?geometrycSsg|] }|ds|qSbperprrrrrrz&read_dataset_input..rlz&do not show the dropped interferogramscs"g|] }|ddvr|qSr?r1r`r)date12_drop_listrrrs"csg|]}j|qSr)rxrrr=rrrrz9WARNING: --show-kept option does not apply to file type: z, ignore and continue.Tcg|]}j|qSr)rxrr=rrrrNzWARNING: input reference date z is not included in input file!zIgnore it and continuerznum of datasets in file z: znum of datasets to exclude: znum of datasets to display: zdatasets to exclude (z): zdatasets to display (zinput reference date: zNo input dataset found!z available datasets: r3)r&rV dsetNumListr r;rr`rrrxr r(open exDsetListZ exDsetNumListrrqplot_drop_ifgramr get_drop_date12_listrXrrrdsetNumr rrr% Exception) r*date1date2objZ drop_num_listftypeZdrop_slice_listrrrr)rr*rry^s &         &   $$  rycs tjdjd}tjdjd}jdkrUjs#d_jsSjrFjdkrFtjdjd}tjdjd}tj ||fj j d_Sjs\tj _t djdd d jdd d jsd_jtjd krjd7_jtjd ksjdkrňjdkrtfd ddDr|d|g}jddjdg}tj||j\__tjtjj_tjt_t dtjt dtjt dtjt dtjjdkr$djvr$tj ||jfj dd_t djjs7d_jjdkr7d_jr]tjjd_tjtjjd\_ _!j!"_!n{tjj#_tjtjj#d_ jdjdjdjd||krj d7_ j$rj d7_ j%dj%ddtj&krj tj%dj%d7_ j'r͈j dj'7_ j(r؈j d7_ jdkrj j!g_nfd d!t)djdD_fd"d!jD_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 r0r1r2rrt)Zds_shaperVr4zfigure size : [z.2fr]gd@c3|]}|jvVqdSrargvrr=rrrrz(update_figure_setting..)z--nrowsz--ncols皙?gffffff?zdataset number: zrow number: zcolumn number: zfigure number: z --figsizeFz'row number is 1, adjust figure size to 2_sub_wrapZ_refZ_excs&g|]}jdt|jqS)r) outfile_baserNrirjr=rrr-sz)update_figure_setting..csg|] }tjj|qSr)rrr#outdir)rr"r=rrr/r)*r\rJrrrrKr|absrTZauto_figure_sizerVr4Zdefault_figsize_multir fig_num fig_row_num fig_col_numrZauto_row_col_numrceilastyperrNrr"rrr$rrvr%rrir<r(r8r9rrrrange)r*rMrLZ data_shapeZ fig_size4plotrr=rrjs N&   "2 "  rjcstj||tjdjdjtjdjdjftjd}tjjjjd}jrIj\}}t|||d|dfdd}t j dkr׈j d vrt d fd d t ||D} |d ddstjjfd| i|d|dd<nt dttjjfd| i|d|dd<j ddrֈjrtjjfd| i|d} t |jdD]} || ddddfdk} || | f| | 8<qnwt dd|d<tj||jd} t ||D][} | j| |dj| dddtjjfdj| i|d}j| dr<jr.rmfloat32complexr3Nz;input data is complex, calculate its amplitude and continuerrAz)reading data as a list of 2D matrices ...r4maxValuer4/rsuffixzconsider input reference date: )rZ inversionc3|]}|jvVqdSr)dsetFamilyListrdr=rrrrz#read_data4figure..) horizontalvertical)rmodelresidualcSsg|] }d|dqS)bandr1rrrrrrrcsr)Z AmplitudeNrrrrrrrcsr)ZPhaseNrrrrrrrTc3rrrrr=rrrr)z-uz--unitzNWARNING: -u/--unit option is disabled for multi-subplots with different units!zIgnore it and continue.masking pixels with zero value masking datac3rrr)rargr=rrrr)z-vz--vlimz--wrapunwrapr )Z num_multilookr4).rzerosrrJrRrrr4 file_ref_yxr&rr;r rryr{rrr(rrr progressBarupdaterVr`closerrrr_rXr7rrrrrrrrrTrrr^r)i_starti_endr*rkrrrrZ ref_kwargsZ dset_listrrrprog_barrZref_dZsame_unit4all_subplotsrrrr=rread_data4figure4s  &*&""   rc Cs|jrtj|d|j|j|j||jd|jdur|jn t |t |g}|j dd|j dd|j dd|j ddf|_ |j ||j|d|d|j|j|j dd}|jrd \}}|jrk|jd|jd}}nd |vr~t|d t|d }}|r|r|j|||j|jd ||j dd||j dd |jr|j|jdkr|g|gdd} | |_ |j!rd} |j"t#vs|j$d%dr|j$|&dd'dd} zt(j)*| d+dd} Wnt| } Ynm|j$|} |j,-| } t.|j/dkr&d| vr&| 'dd} |j0ddkr&|j1-| } |j|j}|dkr4| } n2d|kr@dkrJnn| d| } nd|krVdkrcnn | &dd&dd} n| } | r|j2r}t3|j4d}tj5|| d|dnPt3|j4d}|j$||j6vrd |d!<d"|d#<n/|7d$dd%kr|j8r|d&'|}|d'krd(nd)|d!<|d*'|}| d+|7} |j9| fi||j:r|;|j<r|=|j>s|?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)rrK dem_shade dem_contourdem_contour_seqr*r4rr&r2r0r1)rrrrrrrr'rrr r(rcSsd|dd|ddS)Nr,rrz, v =r)rrrrrr%sz)plot_subplot4figure..format_coordrrr?rz%Y%m%drlr  rr)size)rprop)r2Zcrimsonr3ZboldZ fontweightPLATFORMZSenZSENTINEL1_SENSORAC0C1Z SENTINEL1_IPFz : r4)@rQrTr=rrrr4r^rrrrJrrCrUrrarArZr[rr rLrMrSrTrarrrbrcrdr%rZr;r rVr{rgr`dtdatetimestrptime isoformatrxrr&rr date12List fig_title_inrrZadd_inner_titledropDatasetListrydisp_title4sentinel1r[rdr]rer^r_r`)rr*rrrkr^rgrrr%Z subplot_titledate_strZ title_strZ title_indZ num_subplotrkwargZ s1_sensorZs1_IPFrrrplot_subplot4figures$           rc CsXdt|d|j|d}tdt|tj||j|j|jddd\}}|jj || }|d|j|j}t |j ||j|jg}t||||}t|jtrctj|j|j|jdj|_tdtj|||jd } t||D]T} | |} | j| d|j| d d d t| ||| || d d d d f|d} |jr|js|j| || ddddd} t j!|j"|j| d#dd }|r| $|qv| %~t||t&|D] } |'|| qt()|j*d|j+dgt(,|j*d|j+dgg|_*td|j+d|j-|jrtd|jd|j-dd}|jrv|jsBtdtd|.|jdd|jd|||n9|||d}|jdd krU|d!}td"|j/d#d$|0d%d&|d!d'|g}t1|| |\}} n||||j2rtd(t3j45|j|dd)|j6|j7|j|dd*d|j6d+|j8s|9d 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 z - r1(----------------------------------------T)numrnrowsncolsZsharexZshareyrz plotting ...rrrrN)rrrkgQ?r&rr)rr1shrinkZaspectZ orientationrUNITrrrErcSsJ|jddddddd|js|jr|j|j|jddS|jr#|dS)Ng{Gz?g\(\?g?)rrr r wspacehspace)rr)subplots_adjust fig_wid_space fig_hei_spacefig_tight_layoutZ tight_layout)figr*rrradjust_subplots_layoutis z+plot_figure..adjust_subplots_layoutz-Note: different color scale for EACH subplot!z0Adjust figsize for the colorbar of each subplot.rg?g @r2z show colorbarg(\?)rgGz?rBg{Gzt?save figure to with dpi=tight) bbox_inches transparentdpi):rNr"r rrrrrZcanvasmanagerZset_window_titleflattenr)rrr8rUrTr9rr:rrr4rrrVr`rrVr^Zcolorbarrrqr(ry set_labelrr&Zdelaxesrrdlim_allrrr7Zset_size_inchesrZadd_axesrYsave_figrrrsfig_dpisavefigr!clf)rr*rkrfrZaxsrrrrridxrgrhZ data_unitrZ cbar_lengthrirrr plot_figure#s   "      (rc sttdd|jD|_ttdd|jD|_t|jdkr5|jddkr5ttdd|jD|_|jdkrF|j rFt d |jn|j r_|jdkr_t j|j|j|j|j|j d |_|jd uru|jdkrut|j|j|jd d |_d |_|jd vrdvrtdtd}}tdtd}}d|kr|krnnd|kr|krnn ||g|_td|j|jdkr|jd_n |jdkr|jd_g|_|jdkr|jrt|j}|j dd|!}|jD]|jfdd|D7_qtd|j"rit#$|j"t%fdddDrUtd t&j'(|j"d!t#j)|j"d"|j|j|jdd#d}t j*|||j d$\|_+|_,|_-|Sd |_"d%|_.d&} | d'7} | d(7} t | |S))zPrepare for multiple subplots: 1) check multilook to save memory 2) read existed reference pixel info for unwrapPhase 3) read dropIfgram info 4) read and prepare DEM for background cSh|] }|ddqS)r?rrrrrr rz)prepare4multi_subplots..cSsh|]}|ddqS)StdrF)rgrrrrrrr1rlr cSrrrrrrrrrz/multilook {0} by {0} with nearest interpolation) max_memoryr4Nnearest)methodrrrrnrorz!consider reference pixel in y/x: rg$@dg4@Fr:csg|] }d|qS)r?rrrrrrrz*prepare4multi_subplots..z@mark interferograms with 'dropIfgram=False' in red colored titlec3s |] }||kVqdSrrr)rfrkrrrsz)prepare4multi_subplots..)rnroz reading DEM: z ... height)r3rrrr4)rr*r4rBz;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.)/rrrVrr&rrxrrRr4rXrrrSrTZauto_multilook_numrJrr maxMemoryrrrr;r[rr rrMrrZr r(rrrQrrqrrrr%rZprep_dem_backgroundrrrra) r*rkrrrMrLrZdropDate12Listrrr)rfrrkrprepare4multi_subplotss0         rc@s*eZdZdZd ddZddZddZdS) viewerzviewer class definition. Example: cmd = 'view.py timeseries.h5' inps = cmd_line_parse(cmd.split()[1:]) from mintpy.view import viewer obj = viewer() obj.configure(inps) obj.plot() NcCs$|r |dd}||_||_dS)Nr1)r`rr)selfrrrrr__init__s zviewer.__init__c Cs|jrtndda|jstdt|\}|_t||j}d|_ |j r-t |dkr-d|_ |j D] \}}t|||q2tj|j|j|jd|j|j|j|jd\|_|_|j S)Nc_rrrrrrrrrz"viewer.configure..ZAggrrrr)r4rXr r!rZswitch_backendrrrnr+ update_moder/__dict__itemssetattrrTrr(rrVrJrrr)rr*r;valuerrr configures(   zviewer.configurec Cs|jdkrytdtj|j|jd|jdd\}|_|jr.|tj|j|j|jddd8}|j dvrw|jd dd drwd |j vrwt |jd t |jd }}tj|j|jd|||d|dfddd}||d k|8<|jrtd tj|d k|}|jdurtd tj|jd k|}n tj|jtjd|_t|j}|jdurtd|jtj||jk|}n|durt|std|tj||k|}tj|r|j|j9_|jt|j9_n |jt|9_t||j|\}}|jdurt|jdni}tj |j!|d\}}|j"s1|j#dddddt$|||j|d}|j%rxtdt&j'(|j)dd|j*|j"sc|j+|j)dd|j*d dn |j+|j)dd|j*dd|j,sx|-ncgd} t.t/| t/|j0@} | D] } t1d| dqt2||jd}d d g|_3t4d|j5dD] } t6| ||jdq|j5dkrtdtd|j3d|j7|j8rtd |j8d|j7|j,rtd!t9dS)"Nr1zreading data ...rFrrr?rrrrArrrz%masking pixels with NO_DATA_VALUE of rr)rrr r rrT)rrZ pad_inchesr)rrr)z --show-gpsz --coastlinez --lalo-labelz --lalo-stepz --scalebarz--pts-yxz --pts-laloz --pts-filez WARNING: z= is NOT supported for multi-subplots, ignore it and continue.)rkrzall data range: rEzdisplay range: z showing ...):rr rrr(rVrJrrr;r`r{r[rrrrrrrrrget_no_data_value no_data_valuerrrrrrxrrrrrhrrjr rrrsr"r r r!r rrrrXrrrrrr7r^show) rrrrrZ no_data_valrrrZ opt_namesopt_namerrrrr #s       $  z viewer.plotr')__name__ __module__ __qualname____doc__rr!r rrrrrs   !r)T)F)2rrrr{warningsfilterwarnings UserWarningZ cartopy.crscrsrmatplotlib.pyplotpyplotrnumpyrZmpl_toolkits.axes_grid1rmintpyrrZmintpy.multilookrZmintpy.objectsr r r Zmintpy.objects.gpsr mintpy.utilsr rTrrrrGr/rnr]rrrjrrryrjrrrrrrrrrs@       r3 f O#rdyvr[