合成孔径雷达(SAR)是一种基于雷达原理的先进遥感技术,通过主动发射电磁波并接收地面反射回波,生成高分辨率的雷达图像。它能够在任何天气条件下工作,不受云层、雾气和黑暗环境的影响,为地球观测提供了独特的能力。本文将详细介绍SAR卫星的波段分类、成像模式、产品层级以及正射影像制备与InSAR技术的相关知识,帮助读者通俗理解这一复杂而重要的遥感技术。
一、SAR卫星波段分类及其穿透力特点
SAR卫星根据使用的电磁波波段不同,可分为X波段、C波段、L波段和P波段。波长越长,穿透力越强,但分辨率相应降低。这一特点使得不同波段的SAR卫星适用于不同场景的观测任务。
1. X波段(约3厘米波长)
- 波长范围:3厘米左右
- 频率范围:9-10 GHz
- 穿透能力:较弱,主要穿透树叶表层
- 分辨率:最高,可达0.3米(如德国TerraSAR-X卫星)
- 典型应用:城市精细监测、建筑结构识别、港口设施观测等
X波段SAR卫星具有最高的空间分辨率,能够清晰显示地面细节。例如,TerraSAR-X卫星能在夜间持续追踪伊拉克装甲部队动向,为空袭提供精准目标坐标。然而,它的穿透力较弱,在茂密森林或厚植被覆盖区域效果不佳 。
2. C波段(约6厘米波长)
- 波长范围:6厘米左右
- 频率范围:5.4-5.6 GHz
- 穿透能力:中等,能部分穿透植被
- 分辨率:中等,约5米(如欧洲Sentinel-1卫星)
- 典型应用:海洋监测、农业物候监测、林地和稀疏森林覆盖区域分析等
C波段是SAR卫星的主流选择,如Sentinel-1系列卫星广泛应用于全球环境监测。它在穿透能力和分辨率之间取得了平衡,能够获取上层林冠的信息,同时具备较高的数据获取效率。Sentinel-1卫星的VH极化数据能揭示雨季洪水淹没范围,精度超过90% 。
3. L波段(约23厘米波长)
- 波长范围:23厘米左右
- 频率范围:1.25-1.65 GHz
- 穿透能力:较强,能穿透较厚植被和部分土壤
- 分辨率:较低,约30米
- 典型应用:森林生物量监测、冻土区观测、地下古河道探测等
L波段SAR卫星具有较强的穿透能力,如日本的ALOS-2卫星和我国的陆地探测一号卫星。L波段雷达能深入森林进行生物量测量,反向散射更多来自树干和树枝,而非树叶。在长白山常绿针叶林监测中,L波段SAR数据能有效评估森林健康状况和监测森林退化。此外,L波段SAR曾成功发现撒哈拉沙漠中的地下古河道,展示了其穿透地表的能力 。
4. P波段(约70厘米波长)
- 波长范围:70厘米左右
- 频率范围:300-400 MHz
- 穿透能力:最强,能深入树干层和干燥土壤
- 分辨率:较低,约100米
- 典型应用:全球森林三维结构测绘、地下工事探测、高纬度地区监测等
P波段SAR卫星代表了穿透能力的最高水平,如欧洲空间局(ESA)的BIOMASS卫星。P波段雷达(波长约70厘米)的穿透能力显著优于C波段和L波段,可深入树冠层观测树干、枝干等木质生物量。BIOMASS卫星的观测任务分为层析成像阶段和干涉成像阶段,前者精度可达15-20米(垂直)、200米(水平),将生成人类首个全球三维森林结构图。在军事应用中,P波段SAR可穿透干燥土壤、冰层,识别地下5-10米深的简易弹药库或隧道入口 。
波段选择的考量因素:在实际应用中,选择SAR波段需综合考虑穿透需求和分辨率要求。例如,监测城市建筑时,X波段的高分辨率更为合适;评估森林碳储量时,L或P波段的强穿透能力则更为重要;而进行全球海洋监测时,C波段的平衡特性最为适合 。
二、SAR成像模式:条带模型与滑动聚束模式
SAR卫星根据不同观测需求,采用了多种成像模式。其中,条带模式和滑动聚束模式是最常用的两种。
1. 条带模式(Stripmap)
条带模式是SAR卫星的基本工作模式,特点如下:
- 成像原理:雷达天线以固定角度照射地面,随着卫星沿轨道前进,对地面连续成像
- 天线指向:固定不变,波束宽度较大
- 覆盖范围:较宽,如Sentinel-1的80公里幅宽
- 分辨率:中等,方位向分辨率约为天线长度的一半
- 数据特点:保留幅度和相位信息,以复数形式存储(L1A级产品)
- 适用场景:大范围地形测绘、灾害评估、海洋监测等
通俗类比:条带模式类似于用广角相机连续拍摄地面上的一条长带,能够覆盖更广阔的区域,但细节表现相对有限。
条带模式通过多普勒波束锐化(DBS)技术实现成像,其核心是利用天线在运动轨迹上的位置变化和回波信号的多普勒频移,合成虚拟长天线,从而提高方位分辨率。例如,欧洲航天局的Sentinel-1卫星使用条带模式对布鲁塞尔地区进行成像,展示了安特卫普港、市区和机场跑道等不同地物的特征 。
2. 滑动聚束模式(Sliding Spotlight)
滑动聚束模式是一种结合条带模式和聚束模式的创新技术,特点如下:
- 成像原理:天线波束在地面的移动速度介于0和卫星飞行速度之间,通过控制波束扫描速度增加相干累积时间
- 天线指向:动态调整,波束宽度较小但扫描范围较大
- 覆盖范围:较条带模式窄但比聚束模式宽
- 分辨率:高,可达0.5米(如大兴国际机场测绘案例)
- 数据特点:保留幅度和相位信息,以复数形式存储(L1A级产品)
- 适用场景:城市区域监测、大型基础设施形变分析等
通俗类比:滑动聚束模式类似于用望远镜在移动中连续扫描,既保持了高分辨率,又扩大了覆盖范围。
滑动聚束模式通过控制天线的扫描速度,使波束在地面的移动速度减慢,从而增加了相干积累时间,提高了方位向分辨率。与定点聚束不同,滑动聚束通过控制波束辐射区在地面的运动速度,平衡了分辨率与测绘带宽的关系 。例如,我国"航天宏图一号"卫星采用滑动聚束模式,分辨率达0.5米级,可用于监测城市扩张、基础设施变形等 。
条带模式与滑动聚束模式的对比:
滑动聚束模式在保持高分辨率的同时,显著扩大了覆盖范围,为需要兼顾精度和面积的观测任务提供了理想选择 。
三、SAR产品层级:从原始数据到应用产品
SAR卫星数据经过不同级别的处理,形成了多层次的产品体系。这些产品根据处理程度和用途的不同,具有不同的特性和应用价值。
1. L0级:原始信号数据
- 定义:未经任何处理的原始雷达回波信号数据
- 特点:包含所有原始信息,但无法直接使用
- 存储形式:复数形式(包含幅度和相位信息)
- 用途:作为数据处理链中的原始数据,用于生成更高级别的数据产品
L0级数据是SAR卫星直接接收的原始信号,包含了丰富的信息,但需要经过复杂的处理才能转化为可用的图像。这类数据通常由专业机构处理,普通用户较少直接使用。
2. L1A级:单视复数影像产品
- 定义:经过成像处理和辐射校正的原始数据
- 特点:保留幅度和相位信息,以复数形式存储,具有斜距或地距投影
- 分辨率:较高,保留原始分辨率
- 用途:干涉测量(如InSAR形变监测)、地表变化分析等专业应用
通俗解释:L1A级产品类似于"未经修饰的高清照片",保留了所有原始细节和信息,但需要专业处理才能分析使用。它特别适合需要进行地形变化监测或地表形变分析的场景,如监测地震后的地面破裂情况。
例如,Sentinel-1卫星的L1A级产品(SLC模式)分辨率达3-5米,可用于北加州地震造成的地面破裂监测 。L1A级数据是干涉测量的基础,能够检测到毫米级的地表形变 。
3. L1B级:多视复数影像产品
- 定义:在L1A基础上进行多视处理的数据产品
- 特点:保留幅度和相位信息,以复数形式存储,但进行了多视处理以降低噪声
- 分辨率:较低,通常为L1A的4视处理(如Sentinel-1的GRD模式分辨率为5米)
- 用途:直接视觉分析、地物分类、变化检测等
通俗解释:L1B级产品可以理解为"去噪版高清照片",通过对同一区域的多次观测进行平均处理,减少了图像中的"雪花"(斑点噪声),使图像更加清晰。这种处理虽然牺牲了部分分辨率,但提高了图像的可读性和分析效率。
L1B级产品特别适合于需要直接进行地物识别和分类的场景。例如,Sentinel-1的GRD模式数据可用于农业监测,区分不同作物类型。多视处理还能够在一定程度上消除地形起伏引起的影像畸变,提高图像的几何一致性。
4. L1C级:系统几何校正产品
- 定义:经过系统几何校正的数据产品,形成具有地图投影的图像
- 特点:已进行几何校正,具备地理坐标参考,但未考虑地形起伏
- 存储形式:实数形式(仅保留幅度信息)
- 分辨率:与L1B相同或略低
- 用途:基础地理测绘、地图制作、地物定位等
通俗解释:L1C级产品相当于"已对齐到地图上的照片",它将雷达影像转换为标准地图坐标系(如UTM),便于与其他地理信息数据进行叠加分析。这种产品消除了卫星姿态变化引起的影像偏移,但尚未考虑地形起伏的影响,因此在复杂地形区域仍存在一定的投影误差。
例如,我国陆地探测一号卫星的L1C级产品可应用于地质灾害监测,如监测冻融滑坡活动性、表碛覆盖型冰川表面动力过程等。这类产品为后续的地形校正和应用分析提供了基础 。
5. L2 GEC级:地形校正产品
- 定义:在L1C基础上使用DEM进行地形校正的数据产品
- 特点:结合DEM消除地形起伏影响,生成正射影像或高程模型
- 存储形式:实数形式(仅保留幅度信息)
- 分辨率:与L1C相同或略低
- 用途:高精度地形分析、工程规划、灾害预警等
通俗解释:L2 GEC级产品可以视为"地形修正版地图照片",它利用数字高程模型(DEM)对L1C级产品进行地形校正,消除了因地形起伏导致的影像畸变,使图像更加贴近真实地表形态。这种产品特别适合于需要精确地形信息的场景,如城市规划、灾害评估等。
例如,我国"宏图一号"卫星搭载的多基线干涉SAR系统可提供L2级地形校正产品,高程精度优于2米,具备1:5万比例尺测绘能力 。在三峡库区滑坡监测中,InSAR技术结合DEM数据生成的地形校正产品能够有效识别滑坡区域的形变情况。
6. L3级:几何精校正产品
- 定义:在L2基础上使用地面控制点进行精校正的产品
- 特点:具有原始空间分辨率且空间位置准确,精度最高
- 存储形式:实数形式(仅保留幅度信息)
- 分辨率:与L1C相同或略低
- 用途:高精度工程应用、地形测绘、地图更新等
通俗解释:L3级产品是"经过专业校准的地图照片",通过在L2级产品基础上加入地面控制点进行精校正,进一步提高了图像的几何精度。这种产品适用于对空间精度要求极高的场景,如工程建设、土地规划等。
7. L4级:镶嵌图像产品
- 定义:将多个L3级产品无缝拼接形成的连续大区域影像
- 特点:覆盖范围大,拼接线平滑,色调一致
- 存储形式:实数形式(仅保留幅度信息)
- 分辨率:与L3相同或略低
- 用途:大区域监测、地图底图、环境评估等
通俗解释:L4级产品类似于"无缝拼接的全景照片",它将多个相邻的L3级产品拼接成一个完整的区域影像,便于进行大范围的地表分析和监测。这种产品广泛应用于城市扩张监测、环境变化评估等领域。
SAR产品层级的应用适配:
不同级别的SAR产品适用于不同的应用场景,用户可根据需求选择合适的产品级别:
- 科研与专业分析:通常使用L1A或L1B级产品,因为它们保留了相位信息,便于进行干涉测量和地表变化分析
- 基础地理应用:L1C级产品已具备地理参考,可直接用于地图制作和地物定位
- 高精度地形测绘:L2 GEC或L3级产品结合DEM数据,能够提供更精确的地形信息
- 大区域监测:L4级镶嵌产品适合进行大范围的地表监测和分析
四、正射影像的制备流程
正射影像是通过几何校正消除传感器畸变、地形起伏等因素影响的遥感影像,它兼具地图几何精度与影像特征,广泛应用于土地监测、灾害评估、城市规划等领域。正射影像的制备是一个系统化的过程,主要包括以下四个步骤:
1. 任务规划
任务规划是正射影像制备的首要步骤,涉及:
- 卫星轨道选择:根据观测区域地理位置和范围,选择合适的卫星轨道和过境时间
- 成像时间窗口:考虑地表条件(如植被覆盖、积雪情况)和气象条件(如云层覆盖),选择最佳成像时间
- 覆盖区域划分:确定需要覆盖的地理区域,并规划相应的影像采集方案
- 资源调度:包括卫星数据采集、存储和传输资源的合理配置
通俗解释:任务规划类似于"安排一次高质量的航拍",需要考虑何时拍、拍哪里、用什么设备拍等问题。例如,Sentinel-1卫星在纳米比亚洪灾监测任务中,选择了最佳的观测时间和覆盖区域,以获取最清晰的灾情影像 。
任务规划需要考虑多种因素,如卫星轨道高度、倾角、重访周期等。例如,Sentinel-1卫星在双星组网模式下,重访周期约为6天,能够为区域提供相对稳定的时间序列数据,便于进行地表变化监测 。
2. 辐射校正与几何校正
辐射校正与几何校正是正射影像制备的核心步骤,包括:
辐射校正:消除传感器灵敏度差异、太阳高度角变化、大气散射等影响,使图像灰度值真实反映地表反射特性
- 大气校正:使用6S模型等算法消除大气散射和吸收的影响
- 传感器定标:校正雷达系统本身的辐射特性,确保不同像元之间的灰度一致性
- 地形阴影校正:利用DEM数据消除因地形起伏引起的阴影效应
几何校正:消除传感器内部畸变、外部姿态变化和地球曲率等因素引起的影像畸变
- 系统几何校正:基于卫星姿态和轨道参数,将雷达坐标系转换为地图坐标系
- 地形校正:结合DEM数据,消除地形起伏引起的投影差
- 区域网平差:通过连接点和控制点,优化影像之间的相对位置关系
通俗解释:辐射校正就像"调整相机的曝光和白平衡",使图像能够真实反映地表情况;几何校正就像"修正地图的比例尺和形状",确保图像上的每个像素对应地表的真实位置。
辐射校正中常用的方法包括波段比值法和暗像元法。波段比值法通过将不同波段的图像相除,消除阴影影响,提高地物区分度 。例如,水和沙滩在单波段图像中亮度相近,难以区分,但通过波段比值处理后,它们的差异变得明显。
几何校正中,区域网平差是一个关键步骤。区域网平差通过量测连接点和控制点,在立体或平面区域网内计算卫星影像定向参数,实现多时相影像的联合处理。例如,我国地质调查局应用InSAR技术进行地面沉降监测,通过区域网平差实现了三大地面沉降区的监测成果"一张图",总体精度优于±5毫米 。
3. 区域网构建
区域网构建是正射影像制备的关键环节,主要包括:
- 控制点布设:在观测区域内设置地面控制点,为几何校正提供基准
- 连接点提取:从多个影像中提取共同的地物特征点,作为影像间的连接参考
- 平差计算:基于控制点和连接点,计算影像间的相对位置关系和定向参数
- 立体区域网平差:用于生成数字高程模型(DEM)和数字正射影像(DOM),需量测三维连接点并进行光束法平差
- 平面区域网平差:适用于正射影像生产,仅需平面坐标参与平差计算
通俗解释:区域网构建类似于"制作拼图时的参考线",通过设置控制点和提取连接点,确保各个影像能够准确拼接成一个完整的区域影像。例如,在制作大兴机场的高分辨率正射影像时,需要在机场区域内设置多个控制点,并从多景卫星影像中提取连接点,进行区域网平差处理,确保影像的几何精度。
区域网平差的精度指标因分辨率不同而异。例如,无控制条件下1米分辨率影像平面位置中误差不超过5米,高程中误差7米;优于2.5米分辨率影像的平面中误差不超过7.5米,高程中误差达6米 。
4. 正射影像制备
正射影像制备是整个流程的最终目标,主要包括:
影像镶嵌:将多个校正后的影像拼接成一个完整的区域影像
- 拼接线选取:在影像边缘区域确定最佳的拼接位置,确保拼接后的影像视觉效果一致
- 拼接精度:要求拼接误差不大于2个像元
辐射均衡与色调融合:使不同影像之间的色调和亮度达到一致,提高整体视觉效果
质量验证:检查平面精度是否达到数字正射影像(DOM)制图标准
成果输出:生成标准格式的正射影像产品,如GeoTIFF格式,包含完整的地理坐标信息
通俗解释:正射影像制备就像"将多张照片拼接成一张完整的地图照片",需要确保拼接后的影像没有明显的接缝,色调和亮度一致,能够准确反映地表的真实情况。例如,我国已实现2米分辨率DOM的季度全国覆盖,重点区域优于1米分辨率DOM的年度更新,为国土调查、城市规划提供了重要数据支持 。
正射影像制备过程中,数字微分纠正技术是消除地形起伏影响的核心手段。它将影像分割为微小面积的纠正单元(如像素单元),逐单元进行坐标变换与灰度重采样,最终生成具有统一比例尺的地面正射投影影像 。与传统光学机械纠正相比,数字微分纠正以面元素取代线状缝隙作为纠正单元,通过灰度双线性插值或三次卷积插值算法实现高精度影像重构 。
五、InSAR技术及其产品:DOM、DEM、DSM
InSAR(合成孔径雷达干涉测量)是一种基于SAR影像的相位信息进行地表高程或形变测量的技术。它通过处理同一区域两幅或多幅复数值雷达影像的相位干涉信息,获取地表三维形态或微小形变信息。
1. InSAR技术原理
InSAR技术的基本原理是:
- 数据获取:获取同一地区具有适当基线参数的两幅或多幅SAR影像
- 图像配准:将不同时间或不同视角的SAR影像精确对齐,确保像素对应
- 干涉图生成:通过复共轭相乘两幅复图像,得到干涉图,其中的相位信息反映了地表的高程或形变
- 相位解缠:将干涉图中包裹的相位值(范围[-π, π])转换为连续的相位值,获取真实地表高程或形变信息
- 地形相位提取:通过基线参数和相位信息计算地表高程或形变
通俗解释:InSAR技术类似于"用两把尺子测量同一物体的高度",通过比较两次测量结果之间的差异,计算出物体的实际高度或变化情况。例如,汶川地震后,研究人员利用InSAR技术获取了震后高精度地表形变数据,揭示了地震引起的地表位移场,为评估地震影响和灾后重建提供了重要依据 。
InSAR技术的关键挑战在于大气相位延迟和时间失相干。大气中的水汽和温度变化会影响相位测量精度,而地表变化(如植被生长)会导致图像之间的失相干。近年来,时序InSAR技术通过多时相数据融合和永久散射体分析,显著提高了形变监测的精度 。
2. InSAR主要产品
InSAR技术生成的主要产品包括DOM、DEM和DSM,它们各有特点和用途:
(1) DOM(数字正射影像)
- 定义:利用数字高程模型(DEM)对数字化SAR影像,经逐像元进行投影差改正、镶嵌,按地图图幅范围裁切生成的数字正射影像
- 特点:兼具地图几何精度与影像特征,精度高、信息丰富、直观逼真、获取快捷
- 生成流程:SAR影像→辐射校正→几何校正→结合DEM进行数字微分纠正→镶嵌→DOM
- 应用领域:作为地图分析背景控制信息,评价其他数据的精度、现实性和完整性,从中提取自然资源和社会经济发展信息
通俗解释:DOM就像"一张没有地形畸变的雷达照片",它通过DEM数据对SAR影像进行校正,消除了因地形起伏引起的影像畸变,使图像上的每个像素对应地表的真实位置。例如,互联网卫星地图大多采用DOM形式,既保留了影像的直观特征,又具备地图的几何精度。
(2) DEM(数字高程模型)
- 定义:描述地表高程的数字模型,反映裸露地形的高程特征
- 特点:仅包含地形高程信息,不包括建筑物、植被等地表物体
- 生成流程:SAR影像→配准→干涉图生成→相位解缠→基线校正→DEM
- 应用领域:地形分析、水文模拟、地质勘探、工程建设、无线通讯等
通俗解释:DEM就像"一张地形高低的地图",它用数字表示地表各点的高程,能够生成等高线、坡度图等信息。例如,在工程建设上,DEM可用于土方量计算、通视分析等;在防洪减灾方面,DEM是进行水文分析的基础 。
(3) DSM(数字表面模型)
- 定义:描述地表及地物表面三维形态的数字模型,包含自然地形、建筑物、植被等顶部的高程信息
- 特点:涵盖地表所有要素的高程数据,包括地表物体
- 生成流程:SAR影像→配准→干涉图生成→相位解缠→基线校正→结合多源数据(如光学影像、LiDAR)→DSM
- 应用领域:城市三维建模、遥感影像配准、森林碳储量监测、巡航导弹导航等
通俗解释:DSM可以理解为"一张包含所有地表物体高度的地图",它不仅反映了地面的高低,还包括了建筑物、树木等物体的高度信息。例如,在城市规划中,DSM可用于三维城市建模;在巡航导弹导航中,DSM能帮助导弹"逢山让山,逢森林让森林",实现低空飞行 。
DEM与DSM的区别:
通俗解释:DEM是"裸地地图",只显示地面高低;而DSM是"全要素地图",显示地面以及上面建筑物、树木等物体的高度。当从DSM高程数据中去除植被和人造要素时,会生成一个DEM,即裸地高程模型。
3. InSAR技术在实际应用中的价值
InSAR技术已在多个领域展现出重要价值:
地质灾害监测:如滑坡、地面沉降等。我国地质调查局应用InSAR技术进行地面沉降监测,实现了三大地面沉降区的监测成果"一张图",总体精度优于±5毫米,为灾害预警提供了重要依据 。
地震形变分析:如汶川地震后,研究人员利用InSAR技术获取了震后高精度地表形变数据,揭示了地震引起的地表位移场,为评估地震影响和灾后重建提供了重要依据 。
城市扩张监测:通过多时相InSAR数据,可监测城市区域的地表形变,评估城市扩张对地表的影响 。
森林碳储量监测:BIOMASS卫星的P波段InSAR技术能够直接反演森林主干层的三维结构,提供全球性、系统性的森林生物量变化数据,为碳中和研究提供支持 。
基础设施形变监测:如水电站大坝、铁路、桥梁等。长江科学院应用InSAR技术开发了CK-DSM大坝安全监测云服务系统,为乌东德、白鹤滩、溪洛渡、向家坝等水电工程提供了安全监测管理及分析服务 。
InSAR技术的局限性:
- 大气干扰:大气中的水汽和温度变化会影响相位测量精度,需要通过大气校正模型进行补偿
- 时间失相干:地表变化(如植被生长)会导致图像之间的失相干,影响干涉处理效果
- 处理复杂度:InSAR数据处理涉及多步骤、多参数调整,需要专业软件和算法支持
- 分辨率限制:InSAR技术的空间分辨率通常低于光学遥感,且受基线参数影响较大
随着技术进步,InSAR技术正在向更高精度、更广应用的方向发展。例如,时序InSAR技术通过永久散射体法提升形变反演精度,当前研究聚焦多源数据融合与人工智能辅助处理以克服大气误差及时空失相干等技术瓶颈 。
六、总结与展望
SAR技术作为一种全天候、全天时的主动遥感手段,已在多个领域展现出独特优势。不同波段的SAR卫星具有不同的穿透能力和分辨率,为各种应用场景提供了多样化选择。条带模式和滑动聚束模式则在覆盖范围和分辨率之间取得了不同平衡,满足了从大范围监测到高精度分析的需求。
SAR产品层级从L0到L4,反映了数据处理的逐步深入和应用的逐步扩展。L1A/L1B级产品保留了相位信息,适合专业分析;L1C/L2级产品则提供了地理参考和地形校正,便于直接应用。正射影像的制备流程则确保了遥感影像的几何精度,使其能够作为地图分析的基础。
InSAR技术通过分析SAR影像的相位信息,能够获取地表高程或微小形变信息,生成DOM、DEM和DSM等产品。这些产品在灾害监测、环境评估、工程建设等领域具有广泛的应用价值。虽然InSAR技术面临大气干扰和时间失相干等挑战,但随着多源数据融合和人工智能辅助处理等技术的发展,其应用前景将更加广阔。