o Concatenating two different tracks based on the overlay region) descriptionformatter_classepilogdominantzdominant track. )nargshelp affiliatezaffiliate track. z-bz--bboxSNWE)SNWEzBounding box of area to be geocoded. Include the uppler left corner of the first pixel and the lower right corner of the last pixel)desttypermetavarrz--rewrite_affiliate store_trueFz:whether rewrite affiliate *.h5 file after add the offset. )actiondefaultrz8whether plot profiles for overlay regions of both tracks)titlez --plotpairz;whether plot profiles for overlay regions of both tracks. z --pro_numpro_numzprofile numbers )rrrz--llLONLATzglongitude and latitude of point that you want to profile pass through.use this option when pro_num = 1.)rrrrz-gz --geometrygeometryz geoemtry data of dominant track.zmosaic optionsz-otz--overlapTreatmentoverlapTreatment?aprocess method of incidence and azimuth angle for overlapping region,average: using average value of dominant and affiliate track for the overlapping region;dominant: using dominant swath value for the overlapping region;affiliate: using affiliate swath value for the overlapping region)rrrrz--outputoutputz output name)rrrz--outdiroutdir) argparseArgumentParserRawTextHelpFormatterEXAMPLE add_argumentfloatadd_argument_groupintstr)parserZ plotpair_optmosaicr3N/home/exouser/operations/rsmas_insar/tools/MimtPy/mimtpy/concatenate_offset.py create_parser3s0   r5cCst}|j|d}|S)N)args)r5 parse_args)iargsr1inpsr3r3r4cmd_line_parse\s r:c stdtd}}dvrsz$get_bounding_box..)rr"rcrA)z LON_REF{}rBrDrGr3r4rItrJrr"rrK)r/keysr-npmean) rHlengthwidthlat0lon0lat_steplon_steplat1lon1latslonsr3rGr4get_bounding_boxbs      rYcCs0t|||d}t|||d}||fS)N?r/)rQrRZlat_nZlon_nrSrTrowcolmr3r3r4 calculate_rc{sr^cCs8t|||dd}t|||dd}||fS)NrZrr[)rQrRrUrVrSrTrowscolmsr3r3r4calculate_roco_overlaysracCs||kr|S|SNr3abr3r3r4maxrfcCs||kr|S|Srbr3rcr3r3r4minrgrhc) Cs0t}t|\}}} } t|\} } } }t|| }t|| }t| | }t| |}|jdkr5|}|}|}|}nHt|jd}t|jd}t|jd}t|jd}||kr]|td||kri|tdt||}t||}t||}t||}t||||t|dt|d\}}t ||||t|dt|d\}}t | | ||t|dt|d\}} |dks|dkr|||||||f}!||||| | |f}"|!|"}#t |#}$d | vr|d }%|d }&|d }'|d }(n|d }'|d }(|$|||| ||||f S|||| ||||fS)Nrr"rrKz/input bounding box error! Wrong latitude order!z0input bounding box error! Wrong longitude order!r?r@r=REF_LATREF_LONREF_XREF_Y) r5rYrhrfrr- print_usage Exceptionrar^rM nanmedianrL))r9m_atrm_datas_atrs_datatypeflagr1m_lat0m_lon0m_lat1m_lon1s_lat0s_lon0s_lat1s_lon1Z over_lat0_tmpZ over_lon0_tmpZ over_lat1_tmpZ over_lon1_tmp over_lat0 over_lon0 over_lat1 over_lon1 user_lat0 user_lon0 user_lat1 user_lon1 overlay_rows overlay_colmsm_row0m_colm0s_row0s_colm0 m_overlay s_overlayZoffset_overlayoffsetref_latref_lonref_xref_yr3r3r4calculate_overlaysV         &&&    rc Cs||}dtjd|jddddtjd|jddd}|jdkr3|jd}||}|jrBtj |||d|S)N{}{}{}r.rz_offset.out_filemetadata) rCospathsplitjoinrr'rewrite_affiliaterwrite)r9rrrrs s_data_offsetrout_dir out_dir_filer3r3r4rsJ  rcCs|j\}}tj||ftjdtj}t|D]B}t|D];}|||ftjkrD|||ftjkrDt|||f|||fg|||f<qt|||f|||fg|||f<qq|S)Ndtype)shaperMzerosfloat32nanrangesumnansum)Zmatrix1Zmatrix2r_r`Z matrix_sumr\r]r3r3r4sum_matrix_nans   $((rcR Cs t|}|j\}}|}|j\}}|||| ||| f}|||| ||| f}tt||}tj| | ftjd}t|}t|}d||<d||<||@}||B}tj ||<d||<||}| r|j dkr~||||<d||<n|j dkr||||<d||<||}||| }||| }tj ||ftjdtj } t |\}!}"}#}$t |\}%}&}'}(t |!|%})t|"|&}*t|#|'}+t |$|(},t|d}-t|d}.t|)|*|!|"|-|.\}/}0t|)|*|%|&|-|.\}1}2t|)|*| | |-|.\}3}4|| |/|/||0|0|f<|| |1|1||2|2|f<|| |3|3| |4|4| f<|}5||5d <||5d <|*|5d <|)|5d <|)|5d <|)|5d <|)|5d <|d dkr~|d|5d<|d|5d<|d|5d<|d|5d<|d|5d<|d|5d<|d|5d<|d|5d<n0|d|5d<|d|5d<|d|5d<|d|5d<|d|5d<|d|5d<|d|5d<|d|5d<|jr|jd}6|j}7tj|7ddd}8|8||| ||| f}9| | d|-}:| | d|.};|d}t|=d}?t|=d}@t|=d}At|@|>|A|?}B|jdkrVt|jd}Ct|jd}Dtd|B|C|D||| | |||6 }E|E|9|Ej|Ejf}Ft !|Ej"|Ej#|Ej$|Ej%|Ej&|Ej'|Ej(|Ej)|F|6 | |5fSt *|j| | |:|;||}Gg}Hg}I|GdddfD]s}Jt+}Ktt,|J|B|Gt,|Jddf|Gt,|Jddf||| | |||6 }E|Et,|J|Kd<t-|Ej"|Ej#g|Kd<t-|Ej$|Ej%g|Kd<|Ej&|Kd <|Ej'|Kd!<|H.|Kt+}Lt,|J|Ld<|9|Ej|Ejf|Ld"<|I.|L|Ej(}M|Ej)}Nqnt /|j|H|I\}O}P|d#d$krd$}Qnd%}Qt 0|O|P|M|N|Q|6| |5fS)&zmosaic two tracksrrr"rrrr?r@r;r<r>r=ORBIT_DIRECTIONZ DESCENDINGLAT_REF1LON_REF1LAT_REF3LON_REF3LAT_REF2LON_REF2LAT_REF4LON_REF4height datasetNameZscene_footprintz ([\d+\.]+)rKNNOZp_startZp_endrqrsvalue FILE_TYPEvelocity displacement)1copydeepcopyrrMnanmeanabsonesrisnanrr$rrYrfrhr-r^Zplotpairr'r#rreadrefindallmathatanr r!r Zprofile_extractZrow_noZcolm_noprofilesZ profile_plot lon_start lat_startZlon_endZlat_endZ m_profileZ s_profilem_names_nameZsearch_profilesdictr/arrayappendZprofile_averageZ profiles_plot)Rr9rprqrrrrrrrr}r~rrdsliceZmm_atrm_rowsm_colmsrss_rowss_colmsrrZ ms_differenceZ mosaic_freqZ m_over_posZ s_over_posboth_nanZnone_nanZmosaic_overlay_tmpZmosaic_overlay mosaic_rows mosaic_colms mosaic_datarurvrwrxryrzr{r|Z mosaic_lat0Z mosaic_lon0Z mosaic_lat1Z mosaic_lon1Zmosaic_lat_stepZmosaic_lon_stepZ m_row_loc0Z m_colm_loc0Z s_row_loc0Z s_colm_loc0Zoverlay_row_loc0Zoverlay_colm_loc0 mosaic_atrr'Z geometry_fileZ topographyZ overlap_demrrpolygonZlonlatsZlon_urZlat_urZlon_lrZlat_lrZ angle_radZ point_lonZ point_latZpro_objZ dem_profileZ pro_catalogZprofile_dict_listZprofile_dem_listZpro_NOZ profile_dictZ profile_demrrZprofile_dict_finalZprofile_dem_finalfiletyper3r3r4 mosaic_trackss                                         ,!@       rc Cs|jdkr+dtjd|jddddtjd|jddd}ndd|jdtjd|jddd}|jdkrO|jd}||}tj |||ddS) zwrite mosaic dataNrrrrrz_mosaic.r) r&rCrrrrrr'rr)r9rrrrrr3r3r4 write_mosaics L4  rcCsF|d}|dkr d}|S|dkrd}|S|ddkrd}|Sd }|S) zOjudage dominant/affiliate image is data(velocity or .unw) or dataset (geometry)rr#rrr" coherencerKr)find)rp file_typertr3r3r4judge_data_datasetssrcCs:t|d}t|d}g}g} g} ||krXt|D]7} || } tt|| } | | dk}|jdkrT||| | || | t|| dkdq|}n?t|D]7} || } tt|| } | | dk}|jdkr| || | || |t|| dkdq]| }||| | fS)z;match the date in dominant and affiliate timeseries datasetrg{Gz?r)rdate_list2vectorrMarangerrsizer) m_dateList s_dateListm_dims_dimm_bperps_bperpZ m_datevectorZ s_datevectorm_Dates_DatebperprFdateZdate_subZdate_obj date_finalr3r3r4 date_matchs6     rc1Cst|}td|j}td|j}t|}|dkrgtd|jd}td|jd}t||||||\ }}} } } } } }}t ||||}t ||||||| | | | | || \}}t |||dS|dkrgd}d|j}d|j}t }|D]X}tj||dd}tj||dd}t||||||\}} } } } } }}|dks|dkrt ||||||| | | | | |||d\}}nt ||||||| | | | | || \}}|||<qt |||dS|d kr d|j}d|j}tj|d dd}tj|d dd}t |d }|d }t|}t |d }|d }t|}|jd}|jdd \}} |jd}!|jdd \}"}#t }t||||!||\}$}%}&}'t|%}(t|||dddddf||dddddf|\ }}} } } } } }}||"|})| |#|}*tj|(|)|*fd}+d},t|%|&D]O\}-}.tj||-dd}/tj||.dd}0t|||/||0|\ }}} } } } } }}t ||||0}t |||/|||| | | | | || \|+|,ddddf<}|,d},qtj|'tjd|d <tj|$tjd|d<|+|d <t|$d|d<t |||dS|d krYd|j}d|j}t }t|d}t|d}t||||||\}} } } } } }}t ||||||| | | | | || \}}t |||dSdS)zsimulate LOS displacementrrr) azimuthAnglerincidenceAnglelatitude longitudeslantRangeDistancerrr)rr"rrrrKN)rrrREF_DATE)r:rread_attributerrrrrrrrrrh5pyFiler get_date_listrrlenrMemptyziprrstring_r0)1r8r9rprrrtrqrsrrrrrr}r~rrrrrZ datasliceZm_fileZs_fileZmosaic_datasetrZ m_bperp_daterrZ s_bperp_daterrrrrrrrrrrrZ mosaic_dimrrZmosaic_timeseriesrFm_dates_dateZ dominant_dataZaffiliate_datar3r3r4mains$$  "*$          H  $6    "$r__main__rb),rsysr(rnumpyrMjsonrrrmatplotlib.pyplotpyplotplt mintpy.utilsrrrmintpy.objectsrrrrr Zmimtpy.objects.profilesr objectsrr+r5r:rYr^rarfrhrrrrrrrr__name__r3r3r3r4sD   )K 0 "s