Zf:ddlZddlZddlZddlZej dedddlmZ ddl m Z ddl ZddlmZddlmZmZddlmZddlmZmZmZmZdd lmZmZmZm Z!d Z"d Z#dd Z$dZ%ddZ&dZ'dZ(ggd fdZ)dZ*dZ+dZ,dZ-dZ.dZ/GddZ0y)Nignore matplotlib)categorymodule)make_axes_locatable)subsetversion)multilook_data)TIMESERIES_KEY_NAMESgiantIfgramStackgnss ifgramStack)plotptimereadfileutilsc tdd}|jrd}|Sg}|jD]}|tjj tjj |dtjj|g}|Dcgc]$}tjj|s#|&}}t|dkDr|j|dd}|dk(r_tjj|j}td|D}||kDrd}|St|d|jd|Scc}w) Nzupdate mode: ONskiprunpicrc3ZK|]#}tjj|%ywN)ospathgetmtime).0is D/home/exouser/operations/rsmas_insar/tools/MintPy/src/mintpy/view.py zrun_or_skip..7s7!!!$7s)+z& exist and are newer than input file: z --> skip.)vprintdisp_figoutfilerrjoindirnamebasenameisfilelenappendrfilemin)inpsflagoutfilesfnamefnamesrtitos r run_or_skipr2"s"  D }} Hbggooe& ignore the --ref-lalo input and continue.z"input reference point in y /x : REF_LATREF_LON) disp_unitwrap wrap_ranger8r8c3VK|] }|jjv"ywr)keylowerrrr+s rrz1update_inps_with_file_metadata..s!E!TXX^^%%Es&)) coherencez.corr-)rCcmask?)z.intg? FILE_PATH) datasetNames inps_dictzfigure title:  _nws)7ut coordinatersubset_input_dict2boxvarspix_boxgeo_boxcheck_box_within_data_coverage box_pixel2geowidthlengthdata_boxr strmaxdisp_dem_contourdem_file multilook_num multilookppauto_colormap_namecolormapdsetr8ref_lalo lookup_fileprint geo2radarref_yx radar2geokeysfloat read_pts2inpscheck_disp_unit_and_wrapr;r<r=check_map_projectionvlimanyr@split transparencydisp_dem_shade auto_flipauto_flip_directionflip_lrflip_ud fig_titleauto_figure_titler"replacedisp_whitespacefig_extupdate_figure_setting)r+metadatacoordfbases` rupdate_inps_with_file_metadatarAsH MM( #E!'!=!=d4j(!SDL$,77 EDL&&t||4DL4::t{{3DM %c$--&8 89 %c$,,&7 78 )#e.A.A$--.P*Q QR )#dll*; ;< UV 4<<?T\\!_ , <<?T\\!_ ,.034 % == R S A)) IIaL.. DM }}3DMM?CD||E$5$5 v w E FDM //$--*:DMM!>-- [ !4j   DNN *+ <<&&sB/## VOE '$,,01  &D KTcd|_|jddj|_|jr)|jj dr |j dk(r|js |jrd}|jj dr)tj|_t|dz|S|jj d rd |jvrI|d }tj||_t|d |z|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: 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 plotyx) map_proj_objgetrA coord_unitrS startswith fig_coord lalo_label coastlineccrsrr riUTM ValueErrorre)r+r}r8msgutm_zones rrmrms1$Dll8Y7==?DO **+;<%'DNN2 ?? % %e , $ 0 0 2D  3& '" K__ ' ' 0X]]_,#J/$(HHX$6!sy 334??$%XYY K @@QRSIJ!% Krc |jr|jd|jd|jddz|jddzg|_|jd|jdz }|jd|jdz }|jdkDrbt |t |jdz z |jz }t |t |jdz z |jz }t |j dk(rd|cxkr|j dkr&nn#d|cxkr|j dkr nn|||f}n=tj|j|jd|jdd}tjj|s4tj|s||z}td|jnd |d |d }|d z }t!|d|_nt |j d k(rd|cxkr|j dkrrtd|j>|j>dk(rtj@|}n|j>dk(rtjB|}n|j>dk(r|dz}ny|j>dk(rd|z }nd|j>dk(r|dtjDz z}n>|j>dk(r|tjDdz z}nt'd|j>tjF|tjH|g|_%|j<s4t/jL||j8\|_'|_|_(td|jJd|j2td |j<d|j2||fS)!Nr5rr6Fr7boxr8zset reference pixel to: z WARNING: input reference pixel (, z!) has either masked or NaN value!z -> skip re-referencing.r4 timeserieszinput reference point z is out of data coverage!)r}r;r<r=r8zApply math operation: squaresqrtreverseginverserGrad2deggf@deg2radzun-recognized math operation: r>data range: rKdisplay range: ))rgref_boxrRr]intr'shaperreadr)rbnpma is_maskedisnanr resqueezer@ralltilereshaper_scale_data4disp_unit_and_rewrapr;r<r=r8 disp_scalernmath_operationrrpinanminnanmaxdlim'auto_adjust_colormap_lut_and_disp_limitcmap_lut unique_values)datar}r+ref_yref_xref_valrs rupdate_data_with_plot_inpsrs {{ ADKKN A*DKKNQ,>@  Aa0 Aa0    !T%7%7!%;!<<@R@RRSET%7%7!%;!<<@R@RRSE tzz?a E*DJJrN*qE/JDJJrN/Jue|,"--II $ !  #   55??7+BHHW4E1$++?@8r%Hij11c " _ !E*DJJrN*qE/JDJJrN/J**T!UE/%:;\)"--II $  #   !#9$++F_!`aa55??7+8I7I0JAq 9Atzz!}djjYZm;\]]1$++?@8r%Hij11c "  44 .. YY??.. T ^ _ Y yyOO  '(;(;'<=>  H ,ryy  F *rwwt}  I -   I -rDy  I -ruu 0D  I - 0D!?@S@S?TUVV4"))D/2DI 99 "II$Z^ZhZhi    _TYYKq(8 9: _TYYKq(8 9: :rc  ddlm}|jdd}td|Drt |}g}t j |Dcgc]}|jdc}d}t j |Dcgc]}|jdc}d}t|D]\}tfdt||DrZt||D]J\} } | k(r|dd} | cxkr| kr nn  d|zz } (| k(s. d|dd zz } |j| |~|j|||} | jrtnd at| \} } t!| | } t#j$| j&| j(| j*d| j,| j.| j0| j \| _| _t5j6| j&| j*d| j,| j \}} | j8r=|t5j6| j&| j8| j,d dz}| j:dvr| j*djddjdr{d| j=vrit?| dt?| d}}t5j6| j&| j*d|||dz|dzfd d}||dk7xx|zcc<| j@r#t jBjE|dk(|}| j2.t jBjE| j2dk(|}n4t jF|jHt jJ| _t jBjM|rT| xj2|jNzc_| xj2t jP|jRzc_n)| xj2t jP|zc_tU|| | \}} |r`| jVDcgc]}|dz  }}| jXt[| jXni}t]j^||d}|| | |fS|| | fScc}wcc}wcc}w)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_parser5Nc3>K|]}|jdyw)"'N)r)rxs rrzprep_slice..ms 41<< $ 4rc3DK|]\}}|cxkxr|kncywr)rind0ind1rs rrzprep_slice..vs"It41$$$Is rKrcyrrargskwargss rzprep_slice..r mask_filer7rvminvmaxr8rFrrD unwrapPhaseREF_YREF_XrFdtypeg@ projectionfigsize subplot_kw)0mintpy.cli.viewrrprolistrwhererendswith enumeratezipr(r8rer read_input_file_inforr_ read_maskr)rrbrR mask_vmin mask_vmaxmskrrref_dater@rir zero_maskr masked_whereonesrbool_rmaskrrrfig_sizerdictpltsubplots)cmdauto_figriargs temp_iargsrind0sind1s temp_iargrrtempr+atrrrrref_datarrraxs ` r prep_slicerLs/6/ IIKOE  4e 44%[ ZH!,,{3HI!L:Fa!**[1FGJ%j1 (LAyIs5%7HII"%eU"3JD$Dy(}Di/din 44 T* Y' (" % DnnU*FF%T*ID# )$ 4D!|| ..IIaL LL ^^ ^^.. DHdn  IIaL LL.. ID# }}  II       O# ! ""3'*55mD388:%3w<(#c'l*;u== II ! uQwa0    TRZH$ ~~uu!!$"*d3 xxuu!!$((b.$7774::RXX6 uut TYYJ RXXdii((( RXXd^O#+D#trj@rԉjBrtEjFjHjB}|jK} jM| d | d d d\} } | jd z | jd z f} | ztdjBd| d dd| d dd| ddjN d_(jRr&tjT|j} njVdk(rPtd jVd!|jXfj2d"jZj\d#|} njVd$k(r{td jVd%|j2j"\}}|j^||fjad&jbd'|} |jed(ntgd)jVd*jhrCtjj|jh|jljnj+j(j/d,r)jdjd z d-kDrd_8jprptd.jrtjt|jj(jrjvjxjz/j|rtj~|jjjjjjr jn jxjj0 n"|jd1d2jxd3d3d3d34jPrZjrN|jjd jd jj5td6jX|jjd7d7d fjd7d7d fjjd8d9:td;j>r$tj|||j}jr$t%j<j"\"%&f d<}||_Onyd=_td>jr/tdtj:|d7jjd d?z jdd?z jdd?z jd d?z f_jRr%tjT|j} nډjVdk(r6td@|jXfj2jZdA|} njVd$k(rmtdB|j2j"\}}|j^||fjad&jbd'|} |jed(ntgd)jVd*|jjxCjPrdD\}}jrjd jd }}n.dE|jvrt|dEt|dF}}|r6|r4|j||jj5td6jX|jjd7d7d fjd7d7d fjjdGd9:td;d}|rStjdHdIgdJdKgdJdLgdMdNgdHdIgg}|j|d7d7d f|d7d7d fdOjdGd9:|jj2d d|jj2ddPtjjtjj|dQdR!tjj!rtj! t fdSdTDrMtj!dUjdVd #tj!dWjdVd $n8dXtjj!d}|dYz }t|d7!nd7!!"#$%&f dZ}||_Od7}jrft|}|jjjjtj[}tj| |\}jr2|jjjxj\jrtd]|jjrtd^|jjs|jed_td`jr[jdak(r|jdd3dd3b|jjd3|jjd3n>|jjg|jjgjrjt+jc}jddzd k(rde|df<tj|jfi|tdgjdh|| |fS)iaPlot one slice of matrix Parameters: ax : matplot.pyplot axes object data : 2D np.ndarray, metadata : dictionary, attributes of data inps : Namespace, 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: See prep_slice() above for example usage. cyrrrs rrzplot_slice..rrextentds_shapec|\}}tj|d|d|}tj|d|d|ddd}tj||\}}|j|jfS)a+Get mesh grid coordinates for a given extent and shape. Parameters: extent - tuple of float for (left, right, bottom, top) in data coordinates shape - list of int for [length, width] of the data Returns: xx/yy - 1D np.ndarray of the data coordinates rr5r6r4Nr)rlinspacemeshgridflatten)rrheightrVryxxyys rextent2meshgridz#plot_slice..extent2meshgridsu !  KKq 6!9e 4 KKq 6!9f 5dd ;Q"Bzz|RZZ\))rrvlist)rRrSr8zdisplay data in transparency: )rrr5)cmaprralphazorderrrzplot in geo-coordinater6r4z!draw coast line with resolution: ) resolution linewidthzplotting DEM background ...)rrSdemr+r8)sitez=referencing InSAR data to the pixel nearest to GNSS station: z at [z.6frz] by substrating .3frKFr>imagezplotting data as z! via matplotlib.pyplot.imshow ...upper)rorigin interpolationanimatedscatterz7 via matplotlib.pyplot.scatter (can take some time) ...o)cmarkersequalzUn-recognized plotting style: !)rfaultline_fileSNWErmin_distr8rzplot scale bar: )rrSunitloclabelpad font_sizer) rrS lalo_steplalo_loc lalo_max_num lalo_offsetr.rr8bothinT)which direction labelsizeleftrighttopbottommszplot reference pointNkrG)r=mecmewzplot points of interestc d|ddd|dd} j||\}}|jdz}|jdz}d|cxkrkrn|Sd|cxkrkrn|S ||f}|tj|stjj |rdnd|d z }j rh j||\}}|dz}|dz}d|cxkrkr8nn5d|cxkr kr*nn' ||f}|tj|rd nd |d z }|d |jdzdz }|d|jdzdz }|S)NzE=.f, N=r5r, v=[], v=r, h=[], h=.1fz, x=z.0f, y=)lalo2yxrRrrrrr\)rr rrowcolvdem_rowdem_colhr~ coord_demrrdem_len dem_pix_boxdem_widr+ lalo_digitnum_colnum_rows r format_coordz plot_slice..format_coordsq:,a(Qq Ao,>?C}}Q*HC 4<<? "C 4<<? "CC!'!J'(3&8&8JcN288A;"%%//!2DxDQRSVPW.X=='0'8'8A'>$GW{1~-G{1~-GG-g-!w2H2H 01288A;xd1S'NJc$,,q/1#677c$,,q/1#677Jrrzplotting in Y/X coordinate ...?z.plotting data via matplotlib.pyplot.imshow ...)rrzDplotting data via matplotlib.pyplot.scatter (can take some time) ...)r7NNrrblackiiiiiiirDrHzinputs/geometryRadar.h5c3&K|]}|v ywrr)rr geom_ds_lists rrzplot_slice..sH1 $Hs)latitude longituder`rraz0WARNING: no latitude / longitude found in file: z'skip showing lat/lon in the status bar.cP d|dd|d}ttj| jdz }ttj| jdz }d|cxkrkrn|Sd|cxkrkrn|S||f}|tj|stj j |rdnd|dz } jr(||f}|tj|rd nd |dd z } r$|d ||fd dz }|d ||fd dz }|S)Nx=rIrJrr5rErFrrGrHz mz, E=rBrCrD)rrrintrRrrrr\)rr rrMrLrNrQrr geom_filer+rVlatslonsrWrXs rrYz plot_slice..format_coordsGqgT!C)Cbgga$,,q/123Cbgga$,,q/123CC!'!J'(3&8&8JcN288A;"%%//!2DxDQRSVPW.X==CH Arxx{8$qgR8HHCT$sCx.:,a!@AACT$sCx.:,a!@AACJr)pad axes_class)fontsizecolorzflip figure left and rightzflip figure up and downoffzturn off axis displayr;)r;r: labelbottomlabeltop)rotationZcentervazrotate Y-axis tick labels by z deg)}r8rer tupler isinstancerarYr_ ColormapExtr cmap_vlistr\read_demrRrSrqrrNget_lalo_digit4displayrrrnrrrr coastlinescoastline_linewidthplot_dem_backgroundrO disp_gnssgnss_component ref_gnss_siter get_gnss_class gnss_sourceget_site_lat_lonrfr;disp_ref_pixeldisp_dem_blendplot_blend_imagestyleimshowr animationr r scatter_marker_sizeaxisrr'plot_faultlinefaultline_linewidthfaultline_min_dist disp_scalebarscalebar draw_scalebar scalebar_padr.scalebar_linewidthrdraw_lalo_labelr/r0r1r2lalo_font_sizer tick_paramsrcr ref_markerref_marker_sizepts_lalo pts_markerpts_marker_size plot_gnssrYrgrirpts_yxrarrayset_xlimset_ylimrrr#r$r&rget_dataset_listrrr% disp_cbarr append_axescbar_loc cbar_sizerAxes plot_colorbar disp_title set_titlerw font_colorru invert_xaxisrv invert_yaxis disp_axis disp_tickxaxis set_visibleyaxis get_xaxis set_ticks get_yaxis ylabel_rotsetpget_yticklabels)'rrr}r+r dem_metadatarr(gnss_obj ref_site_lalor rrimrrrYrr temp_testrrcbardividercax tick_kwargsr~rRrrSrTrUr_rerVrfrgrWrXs' ` ` @@@@@@@@@@@@@r plot_slicersL nnU*FF * * *$--% MM]]//  (  }})+ MMLLLLnn * &\; +C0A0A,B BCzzGW**8PJt}}499Q> 6t~~6FG H MMT^^t?W?WM X == 0 1  " " ..   h' >>d11d6H6H GSWSeSe$(NN>6<>BZZ9 $ &tzzl2ij k$T[[$**=FBBgdllnSDD\D\g`fgB GGG =djj\KL L      #222200..   ?? % %e ,$,,q/DLLQRO2SWY1Y!&D     %dmm_5 6    __MM**..11  ??    ..!.. ,,151D1D$--$..,,.. * NN44>> $Dd4  I   4== GGDMM!$dmmA&6OO(<(<  > ) * == $ GGDMM!Q$'q!t)<OO(<(<  % , - >>b$h$..QB == l3I"yy GW  .'/0 == 0 1  " "..  ||As*DLLOC,?||As*DLLOC,?A    $$RsDDNNSB ZZ7 " C D4` 4CUCU`Y_`B ZZ9 $ Y Z$T[[$**=FBBgdllnSDD\D\g`fgB GGG =djj\KL L 0   %LE5{{#{{1~t{{1~uHMMO+"8G#45s8G;L7Muudoo$:N:NO-. ;; " GGDKK1%t{{1a4'8OO(<(<R  ) , - XXT T T T T F GGF1a4L&A,8L8LRY_aG b DKK!$% DKK!$% GGLL+1F!GIbc 77>>) $#44Y?LH.GHH}}YJT\\ejklmn}}YKT\\ejklmnHIYIYZcIdHeegh@@s  I  &' D ~~%b)!!$--T^^`c`h`h!i%%dB4 d  T^^dnnDOO T ||+,  ||()  >> &' ~~ ==H $ NN%TutN T T" T"   $   $ DOO4 ??R 1 $ (K  ##%55.t.?tDE tR rc$tj|j}d}|jjds|dj |dz }|dj |dt j j|jz }d|jvr|dj |dz }td t j jtd tjt|t|d |_t|d |_|d|_t j j%t j j|jd |_t j j%|jd|_td|j|j ftj*|jd|_t/|\}}||fS)Nzinput file is )z.h5z.he5z{} PROCESSORz {} file: {} FILE_TYPE DATA_TYPEz in {} formatzrun z in LENGTHWIDTHrr5zfile size in y/x: T) no_complex)rread_attributer)rformatrrabspathrir r%__file__r version_descriptionrrWrVr@splitextfileBasefileExtget_slice_list sliceListread_dataset_input)r+rrs rrrDs  ! !$)) ,C C 99  o . u||C ,--=  K 0"''//$))2L MMCchhj  %%c+&677 T"''""8,-T'2M2M1N OP 3Kc(m$DKS\"DJ;DHGG$$RWW%5%5dii%@A!DDM77##DII.q1DL djj 9: ;<,,TYY4HDN#4(ID# 9rc 4|Dcgc]}|}}|rt|tr|g}g}|rN|D]H}d|vrd|d}|jdd}||Dcgc]}tj|||c}z }Jn||Dcgc] }||vs| c}z }t t t|}||Dcgc]}|j|c}z }t t t|} | Dcgc]}|| } }| | fScc}wcc}wcc}wcc}wcc}w)z0Get dataset(es) from input dataset / dataset_num*z.*) rtrYryrematchsortedrsetindex) all_listin_list in_num_list search_dsetrtempListdser outNumListoutLists rsearch_dataset_inputrbs0**1*K* gs #iG  Qb=RDBZZD)P1BHHRO4OQPP Q G=qqH}= =H$s8}-.8< 0s) DD+D< D D 1D/ Dc4 t|jdkDst|jdkDr t|jdkDrtd|jd|jdk(r!d|_td|jdnh|jdk(rYt|jd k(rAd |jdvr0|jdj d \}}|g|_||_t|j|j|j|j d |_n5|jd k(r d g|_n|jdk(r dg|_n|jdk(rd|_n|jdk(rOt|j}|jd|jdj ddg|_nQ|j|_|jdk(r1|jDcgc]}|jdr|c}|_t|j|j|j|j d |_t|j|jg|j \|_|_g}t!j"|jd}|j$s|d k(rtdt'|jj)}|jDcgc]}|j dd |vs|}}|Dcgc]}|jj+|}}nt-d|dd|_t/t1t3|jt3|jz t3|z |_|jD cgc]} |j| c} |_t|j|_|jr~|jt6vrd|_t|j|jgg|j dd} | s(d|jd} | dz } t-| d|_n| |_|jd vrtd!t8j:j=|jd"t|jtd#t|jtd$t|jntd!t8j:j=|jd"t|jtd%t|jd&|jtd't|jd&|j|jr*|jt6vrtd(|j|j4dk(rd)} | d*|jz } t?| t!j"|j|jdj dd+} || fScc}wcc}wcc}wcc} w),zFCheck input / exclude / reference dataset input with file dataset listrzinput dataset: "rvelocityFzturning glob search OFF for z filerr5_)rrrrrrHDFEOS displacementgiantTimeseriesreconsr r>rDgeometrybperprz&do not show the dropped interferogramsz9WARNING: --show-kept option does not apply to file type: z, ignore and continue.TNzWARNING: 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: r7) r'rb dsetNumListr r@rrprrrr r)openr exDsetList exDsetNumListrrplot_drop_ifgramrget_drop_date12_listrrerrrdsetNumr rrr% Exception) r+date1date2objr drop_num_listftypedate12_drop_listdrop_slice_listrrrrs rrrs| 499~S!1!12Q6 tyy>A  %dii[2 3 88z !$D  1$((5A B XX %#dii.A*=#STBU99Q<--c2LE5DI!DM0^^II(((( *+, - 88} $&DI XX !'(DI XX* * DI XX+ +"499-C HHuH %q)//4Q78DIDIxx:%(, O1g9NQO /^^II(((( *+, -+?$$ +&'DOT'M  # #DII .{ ;E  M ! ; <*4995JJL *...`QAGGCLOO_<_q`O`>MNT^^11!4NMN MeWTjk l$(D !d3t'7'7#83t?Q?Q;R#RUXYfUg#ghiD,0,<,<=q"=DItyy>DL }} 88/ / DM'^^]]O(( *+, -./ 0 24==/A`aC + +C #J DM$DM xx?")"''*:*:499*E)FbT^^I\H]^_-c$//.B-CDE-c$))n-=>?)"''*:*:499*E)FbT^^I\H]^_&s4??';&s$Z&ZZ 8Z "ZZctjdjdz }tjdjdz }jdk(rjsd_jsj r_j dk(rPtj dj dz }tj dj dz }tj||fjj_Sjstj_tdjdd d jdd d jsgd_jtjz d kDr;xjdz c_jtjz d kDr;jdk(rj dk(rut#fd dDra|dz|g}jddzjdg}tj$j||j\__tjtjj zz _t'j(jj+t,_tdt/jztdt/jztdt/j ztdt/jzjdk(r\dj0vrNtj||j zfjd_tdjjs*d_jj zdkDrd_j2rt4j6j9j2d_t4j6j=t4j6j?j2d\_ _!jBjE_!nt4j6j9jF_t4j6j=t4j6j?jFd_ jdjdz jdjdz z||zkrxj@dz c_ jHrxj@dz c_ jJdjJdz dt&jLzk7r;xj@t/jJdjJdz z c_ jNr"xj@djNzz c_ jPrxj@dz c_ jdk(r!j@jBg_nPtSdjdzDcgc](}j@dt/|jB*c}_j2Dcgc],}t4j6jUj:|.c}_Scc}wcc}w) 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 r4r5r6rr)rrr8zfigure size : [z.2fr]gd@c3:K|]}|jvywrargvrBs rrz(update_figure_setting..sKq*K)z--nrowsz--ncols皙?gffffff?zdataset number: zrow number: zcolumn number: zfigure number: z --figsizeFz'row number is 1, adjust figure size to 2_sub_wrap_ref_exr)+rjrRrr.rrSrabsr_auto_figure_sizerr8default_figsize_multir fig_num fig_row_num fig_col_numrauto_row_col_numrceilastyperrYrr"rrr$outdirrr% outfile_baser{rAr)r<r=rrrranger#)r+rWrV data_shape fig_size4plotjr"s` rr|r|s 4<<?4<<?2 3F $,,q/$,,q/1 2E ||q~~DN}}||% 7T\\!_t||A>?T\\!_t||A>?// %....*DM\ KO}}44DMq!1# 6bq9I#8NaPQ||DL,,uT\\22U: ! ,,uT\\22U:    !d&6&6!&;K4JKK!*e,J!]]1-d2DMM!4DEM131D1D  2 .D d. T\\*U43C3CdFVFV3V-WW wwt||,33C8 !#dll"334!#d&6&6"778!#d&6&6"778!#dll"334   q [ %A// %(8(8"89..!DM .>rww?O?OPTP\P\]^P_?`.a +D t|<<--/DL''//$))4DK " 0 01A1A$))1L Ma PD  Q Q/$,,q/$,,q/2QRUZ[aUaa!!V+!yy!!W,!OOA&);;"%%G%%T__Q-?$//RSBT-T)UU%}}!!VDMM%99!!!U*! <<1 #001$,,@ADL&+1dll1n%=? ! $0013q6(4<<.I?DLJN,,Ww T[[':W K ?Ws -\?1]c  tj||z tjdjdz jz tjdjdz jz ftj }t jjjj}jr&j\}}t |||dz|dzfd}tjdk(rajd vrRtd t||D cgc]} j| } } |jd d j!d s)t#j$j&fd| i|d|ddnFtdtj(t#j$j&fd| i|d|ddjdj!drjrt#j$j&fd| id} t|j*dD]#} || ddddfdk7} || | fxx| | zcc<%n:tdd|d<t-j.||z j} t||D]} | j1| |z dzj| j3ddt#j$j&fdj| i|d}j| j!drNjrBt#j$j&fdj| id}||dk7xx|zcc<||| |z ddddf<| j5j6rLtdj6zt#j$j&fdj6i|d} || z}tjdk(sjdvst9fddDsjgdk(swjttjD cgc] } d| dz c} k(s8t9djDst9d jDrd!}nd}|rt;||\}n-t=fd"d#Drd$}|d%z }t?|d_ jBr.td&tjDjG|dk(|}jH^td'tjJjH|j*dddf}tjDjG|dk(|}tjL|tjN|g_(|rzt9fd(d)Drfjdj!d*r jr<jdd+vr+tSjT|d,d-\_+_,_-|Scc} wcc} w).zIRead multiple datasets for one figure into 3D matrix based on i_start/endr4r5r6rr)rxstepystepr8F)rr8)rrrrrzreading data as a 3D matrix ...rfloat32complexr7Nz;input data is complex, calculate its amplitude and continuerrFz)reading data as a list of 2D matrices ...r8maxValuer8/rsuffixzconsider input reference date: )r inversionc3:K|]}|jvywr)dsetFamilyList)rdr+s rrz#read_data4figure..sP1+++Pr) horizontalvertical)rmodelresidualbandc3>K|]}|jdyw) AmplitudeNrrr's rrz#read_data4figure..sHq1::k*Hrc3>K|]}|jdyw)PhaseNr/r0s rrz#read_data4figure..sD11::g&DrTc3:K|]}|jvywrr)rrr+s rrz#read_data4figure..s8!qDII~8r)z-uz--unitzNWARNING: -u/--unit option is disabled for multi-subplots with different units!zIgnore it and continue.masking pixels with zero value masking datac3:K|]}|jvywrr)rargr+s rrz#read_data4figure..sOSCtyy(Or)z-vz--vlimz--wrapunwrap)r ) num_multilookr8).rzerosrrRr]rrr8 file_ref_yxr'r&r@r rrbrrrrr)r rr progressBarupdaterpcloserrrrorer;rrrrrrrrr_rrrnr)i_starti_endr+r}rrrr ref_kwargsr dset_listrrprog_barr'ref_dsame_unit4all_subplotsrrs ` rread_data4figurerG[s 88  T\\!_t||A .$2D2D DE T\\!_t||A .$2D2D DEZZ  D LL    .. F '' uueU1WeAg>%P  D   A%bb01+0%+@AaTYYq\A A||K3>>yImmDIIO9OOPQRDG P QffX]]499V)VvVWXYZDG   q ! , ,] ;@P@P}}TYYTITTUVWH4::a=) -Aq!G}*QW !,  - :;#{$$eGmt~~Vw& (A OOAK!ODIIaL4F4Fs4KB4OO P diiLTYYq\LVLQOAyy|&&}5$:J:J diiXTYYq\XZXYZ[!Q$5 '(DWa" # (  }}0>?==P PPQRS  D   A%xx66P5OPP""&AA""U3tGZGZC[=\&]acU|&]]HD4G4GHHD0C0CDD!%!&/hE d 8'78 8bC , ,C #J ~~/0uu!!$"*d3 xx~ggdhhA1 56uu!!#)T24"))D/2DIO4NOO((+66x@IYIY##A&i7797a7a 8 4 ty$"4 KuBf'^s =Y >Y%c  |jrDtj|d|j|j|j ||j |j |jn)tj|tj|g}|jddz |jddz |jddz |jddz f|_ |j||j|d|d|j|j |jd}|j"rd \}}|j$r|j$d|j$d}}n.d |j'vrt)|d t)|d }}|r+|r)|j+|||j,|j. |j1|jdd|j3|jdd |j4r|j6|j8zdkDr>|j;j=g|j?j=gd} | |_ |jBrvd} |jDtFvs|jHdjKdrd|jH|jMddjOdd} tPjRjU| djWdd} n|jH|} |jXj[| } t]|j^dk(rEd| vrA| jOdd} |j`ddk(r|jbj[| } |j6|j8z}|dkr| } nJd|cxkrdkr nn| d| } n4d|cxkrdkr&nn#| jMddjMdd} n| } | r|jdr0tg|jh}tjj|| d|ntg|jh}|jH||jlvr d|d <d!|d"<na|jod#dd$k(rL|jpr@|d%jO|}|d&k(rd'nd(|d <|d)jO|}| d*|z } |jr| fi||jtr|jw|jxr|j{|j|s|jd+|S#| } YLxYw),zPlot one subplot for one 3D array 1) Plot DEM, data and reference pixel 2) axes setting: tick, ticklabel, title, axis etc. N)rrS dem_shade dem_contourdem_contour_seqr+r8rrZr6r4r5)rrrrrrrr[rrr<r]r9cd|dd|ddS)NrcrIrJz, v =r)rr s rrYz)plot_subplot4figure..format_coordsAc7$qgU++rrr!rDrz%Y%m%drr rr)size)r,prop)rjcrimsonrkbold fontweightPLATFORMSenSENTINEL1_SENSORAC0C1 SENTINEL1_IPFz : rl)@r\r_r{rIrJrKr8rnrrrrRrrrarrqrrgrirrrrrrrrrrrrrYrr@r rbrryrpdtdatetimestrptime isoformatrrr'r&r date12List fig_title_inrr.add_inner_titledropDatasetListrdisp_title4sentinel1rrurrvrrr)rr+rrr}rnrrrrY subplot_titledate_str title_str title_ind num_subplotrQkwarg s1_sensors1_IPFs rplot_subplot4figurerms  }} nn(( 00nn & -499BIIdORYYt_3UD<<?3& Q(;<<?3& Q(;=DK 4dmm$q'Q!%!3!34;L;L++a  1B  ! u ;;;;q>4;;q>5E   'x013x7H3I5E U GGE5$//d6J6JG KKK Aa !KK Aa ! >>T--0@0@@2E   $   $,"BO  88+ +tyy|/F/Fw/Oyy|++C4::3?CH ) " 4 4Xx H R R TUVWY Z  ! I,,Y7I4&&'1, 1A%OOC03 88K(M9 $ 5 5i @I**T-=-==Ka ) [&B&#,+R { ; k'R' ) 1 1#d ; C CC M #,+    0""2}!$Gdnn599Q<4#7#77%.E'N*0E,' ||J5>4C\C\$,-?$@$F$F$H$K 1:c1Atg!)/!:!@!@!B1!E%3vh7  ]4e4 ||  ||  >>  Is ) ( s 1UU c  dt|d|j|dz }tdt|tj||j |j |jdd\}}|jjj||j}|dz |j z|jz}t|j||j |jzzg}t||||}t|j trEt#j$|j |j&|j(j |_tdt+j,||z |j. } t1||D]} | |z } | j3| dz|j4| j7d d  t9| ||| || d d d d f|} |j:sg|j<rt|j?| || dddd} tAjB|jD|j4| jGdd }|s| jI|| jK~t1||z tM|D]} |jO|| tQjR|jTd|jVdgtQjX|jTd|jVdgg|_*td|jVd|jZ|j<r%td|j<d|jZd}|j:r|j<sOtdtd|j]|j ddz|j d|||n~|||d}|j ddkDr|d z}td!|j_d"#|jad$d%|z d z d&|g}t#jb| |\}} n ||||jdrtd'tfjhjk|j|dz d(|jl|jo|j|dz d)d|jl*|jps|jsy )+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 - r5(----------------------------------------T)numrnrowsncolssharexshareyrz plotting ...rr!rr"N)rrr}gQ?rZr*r))rrhshrinkaspect orientationrUNITrrrKrc|jdddddd|js |jr(|j|j|jy|jr|j y)Ng{Gz?g\(\?g?)r8r9r;r:wspacehspace)r{rz)subplots_adjust fig_wid_space fig_hei_spacefig_tight_layout tight_layout)figr+s radjust_subplots_layoutz+plot_figure..adjust_subplots_layoutsz TT    !3!3   t'9'9'+'9'9  ;  " "    rz-Note: different color scale for EACH subplot!z0Adjust figsize for the colorbar of each subplot.rg?g @r6z show colorbarg(\?)r9gGz?rGg{Gzt?save figure to with dpi=tight) bbox_inches transparentdpi):rYr"r rrrrrcanvasmanagerset_window_titler r*rrGrtrar_rurrvrr=r8rr>rbrprmrrncolorbarrrr)r set_labelr?r'delaxesrrdlim_allrrr;set_size_inchesr|add_axesrsave_figrrrfig_dpisavefigr!clf)rr+r}rwraxsr@rArrDridxrr data_unitr cbar_lengthrs r plot_figurerJsT#a&T\\!A#%6$78I 56 9|| t}}D HC JJ'' 2 ++-C1u(((4+;+;;G w)9)9D  %-4>>RH 7E "*'kAdiil&8&8&=b&AB! t3xc1ai  >>$))<<s3xT#b^h<iD!// tyyQR|TXXY_aefIy)%*& NN 5=#c( + CFYY a 0$))A,?@YY a 0$))A,?@BDM _TYYKq(8 9: yy 1T^^,<=>  ~~yy B C E F    a 03 6 a8H I "3 - #3 -K}}Q#%q ? #   d  +,,s;&95+NOC))$C8JD$ sD) }}ac1B!C DJt||n]^ DLL1%7RVR^R^ _}} GGI rc @  tt|jDchc]}|jddc}|_tt|jDchc]}|j ddc}|_t |jdk(rU|jddk(rCtt|jDchc]}|jddc}|_ |jdkDr1|jr%tdj|jnn|jrb|jdk(rStj |j"|j$|j&z|j(|j |_ |j*A|jdkDr2t-|j*|j|jd |_d |_|j0d vrdj3vrut5dt5d}}t5dt5d}}d|cxkr|kr2nn/d|cxkr|kr$nn!||g|_t7d|j.|j8dkDr|xj:dzc_n$|j8dkDr|xj:dzc_g|_|j0dk(r|j>r{tA|jB}|jEd|jG}|jD]*} |xj<|D cgc] } | d|  c} z c_,t7d|jHr tKjL|jH tO fddDrt7dtPjRjU|jHdtKjV|jHd|j"|j|jdd} tjX| ||j \|_-|_.|_/|Sd |_$d!|_0d"} | d#z } | d$z } t| |Scc}wcc}wcc}wcc} w)%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 rDrStdrLr5rrz/multilook {0} by {0} with nearest interpolation) max_memoryr8Nnearest)methodrrrrrz!consider reference pixel in y/x: r9g$@dg4@Fr>z@mark interferograms with 'dropIfgram=False' in red colored titlec34K|]}||k(ywrr)rrrr}s rrz)prepare4multi_subplots..sK!|A(1+-Ks)rrz reading DEM: z ... r )r7rrrr8)rr+r8rGz;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.)1rrrbrpr&ryr'rrr`r]r8rerr^r_auto_multilook_numrRrr maxMemoryrr r<r@rirr rrrcrrr)rrr\rrrrrr%rprep_dem_backgroundrIrJrKrq)r+r}rrrrWrVrdropDate12Listrrrrrs ` @rprepare4multi_subplotsrs! &J1qwws|A&J!KLD DDWDW&Xqqyyr':&X!YZD 4  1$+)>-)O &O1qwws|A&O!PQ A$.. ?FFtGYGYZ[ D..!3  22 LL$**T-=-==~~nn  xx 2 2Q 6! HH        D xx?"w(--/'A8G,-s8G3D/EuHX./Xg5F1G   1#5#5 %u~D  6t7G7G6HI J ||b #   #D xx= T__$))$ 5!113$$ IA  $HAs!A3Z$H H  IQR }}..t}}= K7JK K ]277#3#3DMM#B"C5I J-- $LL(((( CFHE[E[..F BDND,d.B K !DM #D OC c cC 6 6C #J Ka'K&X&O^%IsR R R R c$eZdZdZddZdZdZy)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() NcJ|r|jdd}||_||_y)Nr5)rprr)selfrrs r__init__zviewer.__init__#s$ IIKOE rc j|jrtnda|jst j dt |\}|_t||j}d|_ |jrt|dk(rd|_ |jjD]\}}t|||tj |j"|j$|j&d|j(|j*|j,|j\|_|_|jS)Ncyrrrs rrz"viewer.configure..+rrAggrrrr)r8rer r!rswitch_backendrrrr, update_moder2__dict__itemssetattrr_rr)rrbrRrrr)rr+r@values r configurezviewer.configure)s...J}}   u %-d3dh-dDHH=    D 1V ;DI----/ &JC D#u % &$&<< IInn !  nn$& $.yyrc , |jdk(rtdtj|j|j d|j d\}|_|jr=|tj|j|j|j ddz}|jdvr|j djddjdrd |jjvr}t|jd t|jd }}tj|j|j d|||dz|dzfdd}||d k7xx|zcc<|jr.td tj j#|d k(|}|j$9td tj j#|j$d k(|}n4tj&|j(tj*|_tj,|j}|j.9td|j.tj0|||j.k(<n;|9tj2|s$td|tj0|||k(<tj j5|rT|xj$|j6zc_|xj$tj2|j8zc_n)|xj$tj2|zc_t;||j|\}}|j<t?|j<ni}tAjB|jD|\}}|jFs|jIddddtK|||j|d}|jLrtdtNjPjS|jTdd|jV|jFs-|jY|jTdd|jVd n,|jY|jTdd|jVd|jZs|j]n gd} t_ta| ta|jbz} | D]} ted| dtg||j}d d g|_4tkd|jldzD]} to| ||j|jldkDratdtd|jhd|jp|jrr%td |jrd|jp|jZrtd!tAjty)"Nr5zreading data ...rFrrrDrrrrFr4r5rz%masking pixels with NO_DATA_VALUE of rr)r8r9r;r:rrT)rr pad_inchesr)rrr)z --show-gnssz --coastlinez --lalo-labelz --lalo-stepz --scalebarz--pts-yxz --pts-laloz --pts-filez WARNING: z= is NOT supported for multi-subplots, ignore it and continue.)r}rozall data range: rKzdisplay range: z showing ...);rr rrr)rbrRrrr@rprrirrrrrrrrrget_no_data_value no_data_valuenanrrrrrrrrrrrzr|rrrrrr"rrr!rrrrrerrrrrr;rnshow) rrrrr no_data_valrrr opt_namesopt_namers rrz viewer.plotJs <<1  % &%]]  IIaLLL !ND$(}} II $  # %&' ((O+ ! **3/2==mL488==?2"488G#45s488G;L7Mu#==II $ ! uQwa8# %&' ( TRZ H, ~~78uu))$"*d;xx#~&uu))$((b.$?774::RXX>#44TYY?K!!->t?Q?Q>RST3566TT///0(+1F>{mLM,.FFT[()uut$TYYJ&RXXdii000RXXd^O+4D$((DIJD$@D?P?P?\):):;bdJll4==ZPGC''##!#Bb$$7:D}}a)I(J*UYUaUaTbcd++KK QTt||`cKdKK QTt||ahKi}}GGI AIS^c$))n<=I% k (+hij k*$BD HDM1dllQ./ 8Atdhh7 8||aAB)$--$..9IJK99-dii[$..9IJK == = ! HHJrr[)__name__ __module__ __qualname____doc__rrrrrrrrs  Bkrr)T)F)1r]r\rrwarningsfilterwarnings UserWarning cartopy.crscrsrmatplotlib.pyplotpyplotrnumpyrmpl_toolkits.axes_grid1rmintpyrr mintpy.multilookr mintpy.objectsr r r r mintpy.utilsrr_rrrrNr2rrmrrrrrrr|rGrmrrrrrrrs ;|L7"+ BA>od/fbLDPl` <,.24 FodaHvrslo dWv^^r