1. 数据准备从GAMMA到StaMPS的关键桥梁搞InSAR时序分析的朋友都知道GAMMA和StaMPS这对黄金搭档用好了能出好成果但数据衔接环节经常让人头疼。我去年处理青藏高原某区域的形变监测项目时光是文件格式转换就折腾了两周。这里分享几个血泪教训首先确保GAMMA输出的*.rslc和*.par文件放在英文路径下路径里出现一个中文字符都可能让后续步骤崩掉。其次建议专门建个工作目录按这个结构组织文件/SBAS_Project /Gamma_Output master.rslc slave1.rslc ... /StaMPS_Input /Processing实测发现最容易被忽略的是DEM参数文件。GAMMA生成的demparm.in需要检查第三行的像素间距是否与影像匹配。有次我遇到mt_prep_gamma报维度不匹配错误最后发现是DEM分辨率设置成了30米而SAR影像实际分辨率是15米。这个小细节会导致后续PS点提取全部失败。2. mt_prep_gamma参数配置与排坑指南2.1 核心参数详解mt_prep_gamma这个脚本是把GAMMA数据喂给StaMPS的关键转换器其参数顺序有严格讲究mt_prep_gamma 20230101 /path/to/data 0.6 5 5 100 100第一参数20230101是主影像日期格式必须与文件名严格一致比如文件名用2023-01-01这里却简写就会报错幅度阈值0.6需要根据场景调整城市区域可提高到0.7-0.8植被覆盖区要降到0.4-0.5分块参数建议先小范围测试我通常先用3x3分块重叠像素设50-100大区域处理再调整到5x52.2 高频报错解决方案错误1Selsbc_patch打不开.rslc文件这个问题90%是因为路径问题。GAMMA对相对路径支持很差建议检查所有输入文件是否在同一目录在mt_prep_gamma里改用绝对路径确保路径不含空格和特殊字符错误2缺少pt2geo输出文件这是GAMMA版本兼容性问题需要修改mt_prep_gamma脚本找到233行注释掉end语句247行的type pt2geo dev/null也要注释掉手动运行pt2geo生成缺失的经纬度文件3. StaMPS初始化配置技巧3.1 参数文件setparm.m的黄金配置经过20个项目验证这套参数组合成功率最高setparm(weed_standard_dev,3); % PS点筛选标准差阈值 setparm(unwrap_method,3D); % 城市区域用3D解缠 setparm(scn_wavelength,800); % 大气滤波波长(m) setparm(plot_scatterer_size,30); % 可视化点大小特别注意ref_lat和ref_lon要设置为研究区中心点有次我直接用了默认值结果相位解缠全部错位。3.2 数据加载避坑指南运行stamps(1,1)时常见三个坑索引超出维度修改sb_load_initial_gamma.m第141-142行参照sb_load_initial.m的矩阵读取方式经纬度乱码把132行的float改为float64矩阵过大报错这是高程和经纬度文件存在0值导致的需要hgt(hgt0) NaN; % 先替换为NaN hgt hgt(~isnan(hgt)); % 再删除4. 相位处理全流程实战4.1 噪声估计的隐藏参数stamps(2,2)做噪声估计时ps_est_gamma_quick.m里有个不起眼但关键的参数max_iter 50; % 默认迭代次数 gamma_thresh 0.3; % 相干性阈值对于低相干性区域如植被覆盖区建议把max_iter加到100gamma_thresh降到0.2。但要注意这会显著增加计算时间我处理300景数据时迭代次数调整后耗时从2小时变成6小时。4.2 PS点选取的实用技巧运行stamps(3,3)时ps_select.m有两个实用技巧在186行后添加代码保存中间结果save(temp_gamma.mat,gamma); % 保存相干性矩阵修改207行的选择逻辑结合幅度离差和相干性双重筛选idx find(gamma0.7 amp_var0.2);4.3 相位解缠的3D优化当stamps(6,6)解缠失败时试试这套组合拳先用默认参数跑一次失败后修改ps_unwrap.m的unwrap_methodunwrap_method 3D; % 替换原来的snaphu调整解缠窗口大小unwrap_win_size 32; % 默认64复杂地形改小去年处理喜马拉雅山区数据时这个调整让解缠成功率从40%提升到85%。5. 结果验证与可视化完成所有步骤后建议按这个顺序检查结果先用ps_plot(w)查看缠绕相位运行ps_plot(u)检查解缠效果最后用ps_plot(v)验证形变速率有个容易忽略的细节默认的120m点间距可能掩盖细节。我习惯在setparm.m里加这行setparm(plot_scatterer_size,15); % 更密集的点位显示对于城市沉降监测还会额外输出时间序列动画ps_plot_ts(persistent); % 生成动态形变图
GAMMA预处理衔接StaMPS的SBAS-InSAR实战：从数据准备到相位解缠的完整避坑指南
发布时间：2026/6/11 16:34:18
1. 数据准备从GAMMA到StaMPS的关键桥梁搞InSAR时序分析的朋友都知道GAMMA和StaMPS这对黄金搭档用好了能出好成果但数据衔接环节经常让人头疼。我去年处理青藏高原某区域的形变监测项目时光是文件格式转换就折腾了两周。这里分享几个血泪教训首先确保GAMMA输出的*.rslc和*.par文件放在英文路径下路径里出现一个中文字符都可能让后续步骤崩掉。其次建议专门建个工作目录按这个结构组织文件/SBAS_Project /Gamma_Output master.rslc slave1.rslc ... /StaMPS_Input /Processing实测发现最容易被忽略的是DEM参数文件。GAMMA生成的demparm.in需要检查第三行的像素间距是否与影像匹配。有次我遇到mt_prep_gamma报维度不匹配错误最后发现是DEM分辨率设置成了30米而SAR影像实际分辨率是15米。这个小细节会导致后续PS点提取全部失败。2. mt_prep_gamma参数配置与排坑指南2.1 核心参数详解mt_prep_gamma这个脚本是把GAMMA数据喂给StaMPS的关键转换器其参数顺序有严格讲究mt_prep_gamma 20230101 /path/to/data 0.6 5 5 100 100第一参数20230101是主影像日期格式必须与文件名严格一致比如文件名用2023-01-01这里却简写就会报错幅度阈值0.6需要根据场景调整城市区域可提高到0.7-0.8植被覆盖区要降到0.4-0.5分块参数建议先小范围测试我通常先用3x3分块重叠像素设50-100大区域处理再调整到5x52.2 高频报错解决方案错误1Selsbc_patch打不开.rslc文件这个问题90%是因为路径问题。GAMMA对相对路径支持很差建议检查所有输入文件是否在同一目录在mt_prep_gamma里改用绝对路径确保路径不含空格和特殊字符错误2缺少pt2geo输出文件这是GAMMA版本兼容性问题需要修改mt_prep_gamma脚本找到233行注释掉end语句247行的type pt2geo dev/null也要注释掉手动运行pt2geo生成缺失的经纬度文件3. StaMPS初始化配置技巧3.1 参数文件setparm.m的黄金配置经过20个项目验证这套参数组合成功率最高setparm(weed_standard_dev,3); % PS点筛选标准差阈值 setparm(unwrap_method,3D); % 城市区域用3D解缠 setparm(scn_wavelength,800); % 大气滤波波长(m) setparm(plot_scatterer_size,30); % 可视化点大小特别注意ref_lat和ref_lon要设置为研究区中心点有次我直接用了默认值结果相位解缠全部错位。3.2 数据加载避坑指南运行stamps(1,1)时常见三个坑索引超出维度修改sb_load_initial_gamma.m第141-142行参照sb_load_initial.m的矩阵读取方式经纬度乱码把132行的float改为float64矩阵过大报错这是高程和经纬度文件存在0值导致的需要hgt(hgt0) NaN; % 先替换为NaN hgt hgt(~isnan(hgt)); % 再删除4. 相位处理全流程实战4.1 噪声估计的隐藏参数stamps(2,2)做噪声估计时ps_est_gamma_quick.m里有个不起眼但关键的参数max_iter 50; % 默认迭代次数 gamma_thresh 0.3; % 相干性阈值对于低相干性区域如植被覆盖区建议把max_iter加到100gamma_thresh降到0.2。但要注意这会显著增加计算时间我处理300景数据时迭代次数调整后耗时从2小时变成6小时。4.2 PS点选取的实用技巧运行stamps(3,3)时ps_select.m有两个实用技巧在186行后添加代码保存中间结果save(temp_gamma.mat,gamma); % 保存相干性矩阵修改207行的选择逻辑结合幅度离差和相干性双重筛选idx find(gamma0.7 amp_var0.2);4.3 相位解缠的3D优化当stamps(6,6)解缠失败时试试这套组合拳先用默认参数跑一次失败后修改ps_unwrap.m的unwrap_methodunwrap_method 3D; % 替换原来的snaphu调整解缠窗口大小unwrap_win_size 32; % 默认64复杂地形改小去年处理喜马拉雅山区数据时这个调整让解缠成功率从40%提升到85%。5. 结果验证与可视化完成所有步骤后建议按这个顺序检查结果先用ps_plot(w)查看缠绕相位运行ps_plot(u)检查解缠效果最后用ps_plot(v)验证形变速率有个容易忽略的细节默认的120m点间距可能掩盖细节。我习惯在setparm.m里加这行setparm(plot_scatterer_size,15); % 更密集的点位显示对于城市沉降监测还会额外输出时间序列动画ps_plot_ts(persistent); % 生成动态形变图
