一、MiniSAR基础成像原理与原始数据预处理
1. MiniSAR成像基本原理
SAR属于主动侧视微波成像传感器,通过雷达平台沿航迹方向运动,合成虚拟大孔径,实现距离向与方位向二维高分辨率成像。MiniSAR作为轻量化版本,雷达天线向侧下方发射线性调频微波脉冲,脉冲接触地表目标后产生后向散射回波,载荷同步记录回波幅度、相位、时间信息。距离向分辨率依靠脉冲带宽调制实现,方位向分辨率依靠平台移动形成合成孔径获得。雷达回波记录的相位信息包含两部分:目标与雷达之间的斜距相位,以及地表目标本身散射引起的固有相位。单视SAR影像仅利用回波幅度生成二维强度图像;而干涉SAR(InSAR)利用同一区域两次观测的回波相位差,反演地表高程,是生成DSM与三维点云的核心理论基础。MiniSAR获取干涉数据一般采用单航过双天线模式或者重复航过模式,单航过双天线方案更适合无人机平台,可消除时间去相干影响,更适用于快速地形测绘。
2. 原始数据预处理流程
原始MiniSAR回波数据不能直接用于干涉处理与产品生成,预处理是整个数据链路的前置基础,处理步骤依次为原始回波解码、运动补偿、距离多普勒成像、辐射定标、多视处理、配准。(1)原始回波解码与数据解析:提取雷达射频原始采样数据,同步读取POS(定位定姿)数据,包含平台GPS三维位置、航向、俯仰、横滚角、角速度等姿态参数。POS数据用于记录载荷瞬时空间状态,是后续几何校正、运动补偿的基础。MiniSAR无人机平台飞行抖动明显,POS采样频率需要与雷达脉冲采样频率匹配。
(2)运动补偿:理想SAR成像假设平台匀速直线飞行,但无人机存在气流扰动,实际轨迹为不规则曲线,斜距误差会带来方位向散焦、相位畸变。运动补偿利用POS数据修正每个脉冲对应的雷达瞬时位置,补偿航迹偏离理想航线带来的斜距误差,提升图像聚焦质量,抑制相位噪声。
(3)距离多普勒成像:对回波信号进行脉冲压缩,完成距离向聚焦;再利用方位向多普勒历史进行方位向聚焦,生成单视复数影像(SLC)。SLC影像每个像素同时保存复数幅度与相位信息,是干涉处理的核心基础数据。
(4)辐射定标:将SLC影像像素幅度转换为地表后向散射系数σ⁰,消除雷达发射功率、天线增益、距离衰减等系统因素带来的幅度偏差,使影像灰度值能够真实反映地物微波散射特性,为DOM产品提供物理量纲。
(5)多视处理:对SLC影像进行方位向、距离向多视平均,抑制SAR图像固有的斑点噪声(Speckle),代价是空间分辨率降低。多视参数需要结合应用需求权衡,应急测绘场景可采用少视模式保留高分辨率,地形高程解算场景采用多视降低噪声。
(6)影像配准:针对两幅SLC干涉影像,以主影像为基准,将副影像在距离、方位向进行亚像素级配准,保证同名像素对应同一地表目标。配准残差是影响干涉相位质量的关键指标,通常要求配准精度优于0.1个像素,否则干涉相位相干性快速下降。
二、DSM(数字地表模型)生成技术
DSM描述包含建筑物、树木等地物在内的地表顶面高程,区别于剔除地物的DEM数字高程模型,MiniSAR生成的高程产品天然为DSM。MiniSAR-DSM基于InSAR干涉相位解算高程,完整处理流程包含干涉图生成、相干系数计算、相位解缠、轨道基线精化、高程反演、地理编码、格网采样。
1. 干涉图与相干系数生成
配准后的主、副SLC影像共轭相乘,得到干涉相位图。干涉相位由平地相位、地形相位、形变相位、噪声相位叠加组成。平地相位由地球曲面和基线几何带来,需要基于轨道参数去除平地相位,得到仅包含地形起伏的残余相位。相干系数衡量两幅影像相位一致性,取值0~1。高相干区域相位稳定,适合高程解算;植被茂密区、水体区域相干系数低,相位噪声大,高程可靠性差。相干系数可作为DSM质量掩膜,用于后续剔除不可靠高程点。
2. 相位解缠
干涉相位观测值被包裹在(-π,π]区间,而真实相位是连续的,相位解缠就是恢复真实连续相位的关键步骤,也是InSAR生成DSM最核心的难点。常用算法包括枝切法、最小费用流法、区域增长法。枝切法适合噪声较少区域,通过识别相位不连续残差点,构建枝切线阻止相位错误传播;最小费用流算法抗噪能力更强,在MiniSAR低空数据中应用更广,能够处理部分不连续相位区域。在叠掩、阴影边界,相位突变严重,相位解缠容易出现跳变错误,需要结合相干系数进行掩膜处理,标记失效区域。
3. 基线精化与高程反演
基线是两次观测天线之间的空间矢量,基线参数误差会直接带来高程系统性偏差。利用地面控制点GCP,对基线矢量进行精化校正,降低基线估计误差。在基线参数、解缠相位确定后,根据InSAR几何模型,相位差与高程建立数学关系,逐像素反演每个像元对应的地表高程。4. 地理编码与DSM格网生成
雷达SLC影像为斜距坐标系,高程结果需要转换到大地坐标系(如CGCS2000)。基于雷达几何模型,将斜距坐标下的高程点投影至地理平面,按固定分辨率规则格网采样,生成栅格格式DSM产品。对于解缠失效的低相干像素,采用邻域插值填充,同时在产品元数据中标记插值区域。完成后进行高程异常检测,剔除明显跳变的粗差,得到最终DSM栅格产品。三、雷达DOM(数字正射影像)生成技术
MiniSAR的DOM是雷达后向散射强度正射影像,不同于光学DOM,灰度代表地表微波后向散射系数,反映地物介电常数、表面粗糙度、几何形态。雷达DOM生成不依赖干涉相位,基于单幅SLC强度影像开展辐射校正、几何正射校正、镶嵌裁切、灰度均衡,生成地理编码的栅格影像产品。
1. 强度影像预处理
提取SLC影像幅度,经过辐射定标得到后向散射系数σ⁰影像,再进行斑点噪声滤波。常用滤波算法:Lee滤波、Frost滤波、Gamma MAP滤波。滤波目的在保留边缘特征前提下抑制斑点噪声,避免地物边界模糊。2. 雷达几何正射校正
SAR侧视成像存在透视收缩、叠掩、阴影几何畸变,无法直接使用光学正射校正模型。MiniSAR雷达正射校正采用距离多普勒(RD)模型,利用POS轨道参数与DSM高程数据,建立雷达斜距坐标与大地坐标之间映射关系。正射校正的核心思路:对于地理坐标系下每个输出格网点,利用DSM获取该点高程,代入RD模型反向计算该点在原始SAR斜距影像中的对应像素坐标,重采样灰度值。重采样一般采用双线性插值,兼顾运算效率与图像平滑度。该步骤必须使用同批次生成的DSM作为高程输入,补偿地形起伏带来的几何畸变,消除透视收缩,得到符合平面坐标的正射影像。
3. 影像灰度均衡与产品输出
单航带雷达影像存在距离向辐射梯度,近距端回波强度普遍高于远距端。采用直方图匹配、分块均衡消除航带内灰度不均匀;多航带作业时,进行航带间辐射匀色与影像镶嵌,接边区域平滑过渡。完成镶嵌后按标准分幅裁切,生成GeoTIFF格式雷达DOM,附带坐标信息、分辨率、辐射定标参数等元数据。雷达DOM优势在于全天候获取地表散射特征,水体、裸土、建筑、植被在雷达图像上呈现明显灰度差异,可用于地物分类、地表覆盖调查;局限在于无法直接表达高程信息,叠掩区域仍然存在地物错位,需要结合DSM联合判读。
四、MiniSAR三维点云生成技术
MiniSAR三维点云为离散三维坐标集合,每个点包含X/Y/Z地理坐标与后向散射强度属性,可直接导入点云软件进行三维可视化、剖面提取、地物建模。MiniSAR点云有两条生成技术路线:基于InSAR栅格DSM采样点云、基于方位距离坐标直接反演的稀疏点云,工程中以DSM栅格采样方案为主。
1. 基于DSM栅格生成密集点云
DSM是规则格网栅格,每个栅格单元对应一组平面坐标与高程值,遍历所有有效栅格,将栅格中心转换为三维坐标(X,Y,Z),并同步赋予该位置DOM对应的后向散射强度,形成密集三维点云。低相干、解缠失效区域的栅格直接剔除,不生成点,减少粗差点。该方法实现简单,点云密度由DSM分辨率控制,处理效率高,适合常规测绘生产。该路线的局限性:点云严格继承DSM的高程误差,相位解缠错误会直接转化为点云高程粗差;格网采样点呈规则排列,缺乏真实地物边缘的不规则采样特性。
2. 基于单像素雷达几何直接反演稀疏点云
利用SLC影像每个像素的斜距、多普勒中心、平台POS参数,结合干涉解缠得到的相位,直接求解每个散射中心三维坐标,不经过栅格DSM中间产物。该方式得到的点云采样位置为SAR原始成像网格,并非均匀地理格网,更贴近雷达散射中心真实位置,对陡峭山体、建筑立面散射点还原效果更好,但点分布不均匀,存在大量噪声点,需要进行点云滤波。3. 点云后处理与质量优化
原始MiniSAR点云含有大量噪声,需要多级滤波处理。(1)采用统计离群点滤波,剔除邻域距离异常的孤立粗差点;
(2)采用体素下采样,在保留地物形态前提下降低点云数据量;
(3)针对植被区域,可根据相干系数阈值分离地表散射点与冠层散射点。
完成滤波后的点云可输出LAS/LAZ标准点云格式,点属性包含三维坐标、强度、相干系数、点分类标记。MiniSAR点云可用于地形剖面提取、土石方量估算、建筑物轮廓三维可视化。与激光雷达点云对比,MiniSAR点云高程精度相对偏低,但不受光照云雾限制,适合应急场景快速获取三维点云成果。
五、三类产品特征对比与质量控制
1. 产品信息特征对比
(1)DSM为规则栅格高程产品,核心承载高程信息,适合等高线生成、地形分析、土方计算;(2)雷达DOM为栅格强度影像,承载地物微波散射特征,用于地物识别、边界判读;
(3)三维点云为离散矢量三维数据,适合三维可视化、局部细节剖面分析。
三类产品由同一MiniSAR观测数据派生,在工程测绘中经常配套使用,DOM用于平面地物判读,DSM提供高程背景,点云用于三维细节展示。
2. 质量控制要点
(1)DSM质量:检查高程粗差、相位解缠跳变、叠掩区域高程失真;利用地面控制点统计高程中误差;使用相干系数图标记低可靠区域。(2)DOM质量:检查几何接边误差、辐射不均匀、斑点噪声残留;验证地物平面定位精度,叠掩区域单独标注。
(3)三维点云质量:统计点云粗差率、点密度,剔除离群噪声点,评估地物轮廓完整性。
3. 典型误差来源
MiniSAR低空航飞的核心误差来源:(1)平台姿态抖动引入的运动误差;
(2)基线估计偏差带来高程系统性偏移;
(3)地形叠掩、阴影区域相位失效;
(4)植被体散射造成时间去相干;斑点噪声降低相位解缠稳定性。
在工程作业时,可通过优化航飞方案、增加地面控制点、多航带观测、合理选择滤波参数降低误差影响。
MiniSAR凭借轻量化、全天候遥感能力,能够在低空平台快速获取微波回波数据,经过信号预处理、干涉解算、几何校正,可生成DSM、雷达DOM、三维点云三类标准化测绘产品。DSM依靠InSAR干涉相位实现地表高程重建;雷达DOM基于单幅强度影像结合DSM开展距离多普勒正射校正,消除SAR侧视几何畸变;三维点云可从DSM栅格采样或者雷达散射中心直接反演得到,实现地表三维离散表达。
MiniSAR聚焦于微型合成孔径雷达(SAR)制造研发,为用户提供定制化机载SAR、轻型MiniSAR、无人机载MiniSAR、SAR数据采集服务、SAR飞行服务等。如您有相关业务需求,欢迎联系!