stamps(1,1)
stamps(2,2)
%%修改这个参数设置为n可以大幅缩短运行时间%%
setparm(select_reest_gamma_flag,n)
%%如果你的研究区包含海洋等想去除的部分，将下面这个参数改为y%%
weed_zero_elevation
stamps(3,3)
%%修改这个参数可以将同相邻点的位移量，如果是取自同一个像素中心的情况去除%%
setparm(weed_neighbours, y)
stamps(4,4)
stamps(5,5)
%%以下是使用GACOS进行大气对流层校正(gacos_datapath和UTC_sat参数根据实际修改)%%
getparm_aps
load(parms.mat)
save(parms_aps.mat, heading, lambda, -append)
setparm_aps(gacos_datapath, ./gacos)
setparm_aps(UTC_sat, 10:43)
%%以下代码是为了解决：the error la2.mat is empty with no value.%%
load inc2.mat
load la2.mat
la=inc
save(la2.mat, la)
aps_weather_model(gacos, 1,3)
stamps(6,6)
setparm(scla_deramp, y)
stamps(7,7)
%%以下步骤是保存将研究区的PS点位移结果输出到stamps_visualizer的csv文件%%
ps_plot(v-dao, a_gacos, ts)
load parms.mat
corr_method = a_gacos
ps_plot(v-dao,corr_method, -1)
load ps_plot_v-dao.mat
lon2_str = cellstr(num2str(lon2))
lat2_str = cellstr(num2str(lat2))
lonlat2_str = strcat(lon2_str, lat2_str)
lonlat_str = strcat(cellstr(num2str(lonlat(:,1))), cellstr(num2str(lonlat(:,2))))
ind = ismember(lonlat_str, lonlat2_str)
disp = ph_disp(ind)
disp_ts = ph_mm(ind,:)
export_res = [lon2 lat2 disp disp_ts]
metarow = [ref_centre_lonlat NaN transpose(day)-1]
k = 0
export_res = [export_res(1:k,:); metarow; export_res(k+1:end,:)]
export_res = table(export_res)
writetable(export_res,research_gacos.csv) %%这后面是输出csv文件的名字
