Project

General

Profile

display displacement rate in Google Earth

Added by Gy Lona almost 7 years ago

Dear Gamma users,

I want to display a deformation time series in Google Earth. After determining the linear LOS deformation rate and calculate the coherence I also generated an ascii file with disp_prt. My problem is that the LOS displacement values (written into the kml file) are not the same if I compare them with the vu_disp command (LOS rate from kml and slope from vu_disp differ - see attached png). Could you help me please where did I make a mistake?
The following commands are executed:

#generate deformation rate
def_mod_pt pt - pSLC_par - itab_ts pbase_ts 0 pdiff_ts_Ba 0 $REF_POINT pres_B pdh_B pdef_B punw_B psigma_B pmask_B 60. -0.2 0.2 5.0 5 pdh_err_B pdef_err_B
#display result
vu_disp pt pmask_B pSLC_par itab_ts pdisp_ts_B pdef_B phgt2 psigma_B pdh_err_B pdef_err_B - pdef_B.ras -0.16 0.04 2 128

#generate kml file
pt2geo pt pmask_B 20170502.rslc.par - phgt2 ../DEM/EQA.dem_par ../DEM/20170502.diff_par 5 1 pt_map pmap pmapll
pt2d pt_map pmask_B pdef_B 1 P.def1 400 600 10. 10. 1 1 2 3 3 1
ras8_float P.def1 - 400 P.def1.bmp 1 0. 240. - - - 0 -0.025 0.025 0 0 1. .35 1 1 0 1 1 1
ras_dB EQA.20170502.mli 462 1 0 1 1 -22. 3.5 0
comb_hsi P.def1.bmp EQA.20170502.mli.bmp P.EQA.def1.map.bmp
ras_pt pt_map - P.EQA.def1.map.bmp P.EQA.PALSAR.def1_025.bmp 1 1 255 0 0 15 0 2 4004

disp_prt pt_map pmask_B - pSLC_par itab_ts pmap phgt2 pdef_B pcct_B pdh_err_B pdef_err_B pdisp_ts_B 4004 P_S1_items.txt P_S1_disp_tab.txt

kml_pt P_S1_disp_tab.txt 5 4 7 "LOS rate [mm/year]" 6 "height" 8 coherence 1 P_S1_disp.tmp.kml button_master.png gamma_logo.png - 1 0. 240. 0.7 0.7 1.0 1.0 0 -25.0 25.0 0
sed 's/>0.20.4 P_S1_disp_025.kml
zip P_S1_disp_025.kmz ./P_S1_disp_025.kml button_master.png gamma_logo.png pdef_final.colors.jpg