StaMPS中进行PSInSAR解析的命令

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文件的名字