一、SAR载荷成像几何与地形畸变机理
1. 斜距投影与成像坐标系
SAR以侧视方式发射并接收电磁波，图像中每个像元由两个参数唯一确定：方位向采样时刻与斜距。方位向对应雷达沿轨飞行的位置，斜距则是天线相位中心到地面目标沿视线方向的直线距离。因此原始SAR图像本质上是地物在斜距方向的投影，与光学影像的中心投影、正射投影均不相同，其坐标系统称为斜距—方位坐标系。在这一坐标系下，地面等距分布的像元在斜距上并不等距，近距端地面分辨率低于远距端，图像几何随地形与轨道姿态剧烈变化。2. 地形起伏引发的几何畸变
当观测区域存在高程起伏时，斜距投影会产生三类典型几何畸变：透视收缩、叠掩与阴影。（1）透视收缩指面向雷达的斜坡上坡面在斜距方向被压缩，坡度越大收缩越严重；
（2）叠掩指高大地物或陡坡顶部先于底部到达雷达，使顶部像元"倒挂"在底部之前；
（3）阴影则指背向雷达的坡面收不到信号，形成无回波区域。
这三类畸变使SAR图像无法直接与地图坐标对应，若仅按椭球面高程为零处理，山区像元的平面位置偏差可达数百米甚至数千米，这是几何校正必须引入DEM的根本原因。
3. 高程误差与定位误差的定量关系
高程对定位误差的影响可用斜距几何定量描述。设局部入射角为θ，目标在斜距上的位置随高程变化的关系近似为δR≈δh·cosθ，对应地距方向的偏移约为δx≈δh/tanθ。以典型入射角30°为例，10米的高程误差将在地距方向引入约17米的定位偏差；入射角增大到45°时，同样的高程误差仍会造成约10米的偏差。可见入射角越大、地形越陡，高程误差对定位精度的放大作用越强，高精度几何校正对DEM高程精度的需求也随之提高。二、几何校正的目标与总体框架
几何校正的目标是将斜距—方位坐标系下的SAR载荷图像变换到指定地图坐标系（如WGS84下的UTM投影或高斯投影），使每个像元获得准确的平面坐标，并恢复地物真实形状。总体框架通常包括四个环节：轨道与成像参数解算、定位模型的建立、DEM辅助的坐标解算与灰度重采样、以及精度评定。定位模型是核心，其作用是根据成像时刻的轨道状态与成像几何，建立像元坐标与地面三维坐标之间的映射关系。DEM在其中扮演高程约束的角色：当定位方程缺少独立高程信息时，DEM提供地表高程作为先验，将二维定位问题转化为在DEM曲面上的三维解算，从而闭合方程组，消除地形引起的系统性误差。根据是否有地面控制点参与，校正又分为有控模式与无控模式，前者依赖地面控制点精化系统参数，后者则完全依赖精密轨道与DEM实现全球定位，是当前遥感业务化处理的主流方向。
三、DEM融合方法的技术路线
1. 距离—多普勒定位模型与DEM迭代解算
距离—多普勒（Range-Doppler, RD）模型是SAR几何校正最经典的定位模型，由三个方程构成：（1）斜距方程描述像元斜距等于天线相位中心与地面目标的空间距离；
（2）多普勒方程描述回波多普勒频率与相对运动速度的关系；
（3）第三个方程描述目标位于某一参考面。
当采用DEM融合时，第三个方程不再使用参考椭球，而是将地面目标约束在DEM确定的高程曲面上。求解时，先以轨道参数与初始高程估计目标的三维坐标，计算斜距与多普勒残差，再用残差修正目标位置并重新读取DEM高程，反复迭代直至残差小于阈值，通常数次迭代即可收敛。该方法的优势在于物理意义明确、模型严密，在精密轨道支持下可实现米级至亚米级的定位精度。
2. 高程约束下的地距重采样与像元灰度重组
定位模型求解得到每个地面网格点的SAR载荷图像坐标后，需要将原始图像的灰度重采样到规则地距网格，生成正射影像。常用做法是反向映射：对输出网格的每个像元，利用DEM高程与定位模型反算其在原始斜距图像中的位置，再通过插值提取灰度。插值方法包括最近邻、双线性与三次卷积，其中三次卷积对边缘保持和纹理保真效果最好，但计算量最大。重采样同时要考虑叠掩与阴影区域的处理：叠掩区域应选择最先到达的信号，阴影区域则标记为无效值，避免伪信息进入正射影像。DEM网格密度与SAR分辨率匹配时，重采样后的影像地形变形得以消除，地物几何关系与真实地表一致。3. 干涉测量中的地形相位补偿
在重轨干涉SAR（InSAR）形变监测中，DEM融合体现为地形相位去除。干涉相位由平地相位、地形相位、形变相位与大气相位叠加而成，其中地形相位仅由观测几何与高程决定。利用外部DEM计算模拟地形相位并从干涉图中减去，即可生成仅含形变、大气与噪声的差分干涉图。地形相位的敏感程度由高程模糊度刻画，其表达式为h_a=λR·sinθ/(2B_⊥)，其中λ为波长，R为斜距，θ为入射角，B_⊥为垂直基线。以C波段（波长约5.6厘米）、斜距850公里、入射角30°、垂直基线100米为例，高程模糊度约为119米，意味着100米高程变化对应一整周干涉相位；若所用DEM存在10米高程误差，将产生约0.08周的残余相位误差，在形变解算中对应约数毫米的误差贡献。因此InSAR应用中，DEM高程精度直接决定差分干涉的可靠性与形变监测的精度下限，实践中常选用TanDEM-X等高精度DEM作为地形相位参考。4. 有理多项式模型与DEM联合精化
有理多项式（RPC）模型以比值多项式近似描述像点坐标与地面坐标的映射关系，因形式统一、不依赖具体传感器物理参数，被广泛应用于星载SAR影像的几何处理。但RPC模型在无控制点条件下受轨道误差影响较大，单独使用难以满足高精度需求。将DEM融入RPC模型的常用做法是：以DEM高程作为第三维约束，解算RPC附加的仿射变换参数，通过区域网平差同时消除系统性的平移、旋转与尺度误差。此方法将物理模型与统计模型结合，既保留了RPC的通用性，又利用DEM约束了高程维的退化问题，尤其适合缺乏精密轨道的国产SAR卫星数据的业务化校正。5. 辐射地形校正与局部入射角补偿
除几何位置外，地形起伏还通过改变局部入射角影响雷达后向散射强度，造成同一地物在不同坡向上的辐射差异。DEM融合的延伸应用是辐射地形校正：利用DEM计算每个像元的局部入射角与投影面积，对后向散射系数进行归一化处理，削弱透视收缩区的辐射增强和背坡区的辐射减弱，使影像辐射特征更接近地物真实散射属性，提升后续分类与定量反演的精度。四、DEM数据源与精度要求
DEM的质量直接决定融合校正的精度上限。目前全球公开与商业DEM主要数据源包括：
1. SRTM全球数字高程模型，其1弧秒版本分辨率约30米、3弧秒版本约90米，绝对高程精度在开阔平原约±6米、整体约±16米，覆盖范围为南北纬60度之间；
2. ALOS World 3D（AW3D30）分辨率30米，由PRISM立体像对生成，标称高程精度约±5米，覆盖范围扩展至南北纬82度；
3. TanDEM-X WorldDEM由德国双星干涉测量生成，分辨率12米，无冰区相对高程精度约0.8至2米、绝对高程精度约3至4米，是当前全球最高精度的商业DEM；此外，ICESat-2与GEDI激光测高数据常作为DEM高程精度验证的独立基准。
工程实践中，高精度几何校正通常要求DEM高程误差控制在数米以内，分辨率不劣于SAR载荷影像地距分辨率的数倍；在陡峭山区，还应优先选择相对精度更高、空洞更少的TanDEM-X或局部高精度DEM，以避免地形细节丢失造成残差。DEM与SAR影像的坐标系、基准面不一致时，还需先进行椭球基准与投影变换的统一处理。
五、精度评估与典型结果
几何校正精度通常以地面控制点（如GPS实测点、人工角反射器、影像匹配同名点）为基准，统计平面定位误差的均方根误差（RMSE）与圆误差概率（CE90、CE95），高程方向采用LE90。典型业务结果表明：在精密轨道与高精度DEM支持下，星载SAR无控制点定位精度可从原始图像的数百米量级提升至10至20米甚至米级；加入地面控制点的有控模式可将精度进一步压缩至数米以内。精度随地形变化明显，平原地区高程误差影响小，精度主要受轨道与姿态误差制约；山区则受DEM精度、叠掩阴影分布与重采样插值共同影响，误差通常大于平原。评估时应按地形类别分层统计，并单独报告叠掩与阴影区域的无效像元占比，避免总体指标掩盖局部退化。
DEM融合方法通过将地表高程作为先验几何信息嵌入距离—多普勒、RPC等定位模型，从根本上消除了地形起伏对SAR载荷影像定位的系统性影响，是星载与机载SAR几何校正的关键环节。从成像机理出发理解畸变来源，依据任务需求匹配DEM数据源，并严格评估分层精度，是该方法工程化应用的核心要义。
MiniSAR聚焦于微型合成孔径雷达（SAR）制造研发，为用户提供定制化机载SAR、轻型MiniSAR、无人机载MiniSAR、SAR数据采集服务、SAR飞行服务等。如您有相关业务需求，欢迎联系！