1. DEM配准与地理编码的核心概念
第一次接触GAMMA处理SAR数据时,最让我头疼的就是DEM配准和地理编码这两个环节。简单来说,DEM配准就是把外部数字高程模型(比如SRTM或ASTER数据)与SAR影像对齐的过程,而地理编码则是把斜距几何的SAR数据转换到地理坐标系下。听起来简单,但实际操作中会遇到各种"坑"。
举个例子,去年处理一批Sentinel-1数据时,我发现直接用原始DEM会导致配准偏移200多米。后来排查发现是没考虑大地水准面转换——DEM用的是正高系统(以海平面为基准),而SAR测量的是大地高(以椭球面为基准)。这个细节在GAMMA的dem_import命令中通过geoid参数体现,选择EGM2008模型后误差立刻降到了10米内。
需要模型API调用? 免费领10W Token,多模型网关一键接入 Claude、DeepSeek 等主流模型。
2. DEM数据准备与导入
2.1 DEM格式选择与预处理
GAMMA支持多种DEM格式,但最常用的是GeoTIFF。我习惯先用QGIS检查DEM质量:
- 检查是否有空洞(NODATA值)
- 确认覆盖范围大于研究区
- 分辨率建议比SAR数据高2-3倍(比如Sentinel-1用30米DEM)
bash复制# 用gdal检查DEM信息
gdalinfo input_dem.tif
2.2 关键参数解析
执行dem_import时,这几个参数最容易出错:
- geoid:必须选对模型(EGM96或EGM2008)
- dem_proj:DEM的投影参数
- region_cut:可以裁剪研究区域
bash复制dem_import input_dem.tif dem.par dem - - - EGM2008
注意:如果后续出现配准偏差,90%的问题都出在DEM导入环节,务必反复检查投影参数。
3. 初始查找表生成与优化
3.1 gc_map2的实战技巧
生成初始查找表时,lat_ovr和lon_ovr这两个过采样因子很关键。我的经验值是:
- 平坦地区用1.5
- 山区用2.0-3.0
- 超高精度需求可以到5.0
bash复制gc_map2 20210311.vv.mli.par dem.par dem dem_seg.par dem_seg lookup_table 2.0 2.0 -
