TerraSAR-X雷达卫星单极化SM模式介绍
TerraSAR—X卫星及其在地球科学中的应用
2007年总第28annalsofshanghaiastronomicalobservatorycasno282007terrasarx卫星及其在地球科学中的应用中国科学院上海天文台上海200030月15日德国成功发射了terrasarx卫星即将提供高达1m分辨率的多模式多极化的星载sar数据其设计寿命为5该文较全面系统地介绍了这一现代高分辨率雷达卫星的有关情况主要涉及它的科学目标结构和设计参数处理算法和基本产品类型以及在地球科学中的应用等方面冀期以此为我国开展相关领域的研究工作提供参考
Ｔ ｒａ Ａ Ｘ 卫星 及 其 在 地 球 科 学 中 的应 用 ｅｒＳ Ｒ—
陈艳玲 黄
提
城
要
冯天厚
（ 中国科学 院上海天 文台 ， 上海 ２ ０ ３ ） ００ ０
２０ 年 ６ ｌ 德 国成功发射 了 Ｔｒ Ｓ ＲＸ卫星 ， ０７ 月 ５日 ｅ ａＡ － ｒ 即将提供高达 ｌ ｍ分辨率 的多模 式 、 多极 化的星载 Ｓ Ｒ数据 ， 计寿命为 ５ 。该 文较全 面 、 Ａ 其设 年 系统地介 绍 了这 一现代 高分辨 率雷达卫 星 的有关 情况 ， 主要涉及它 的科学 目 、 标 结构和设计参数 ， 处理算法 和基本产品类型 ， 以及在地球科 学 中的应用等方面 ， 以此 为我 国井 展相 关领 域的研究工作提供参考 。 冀期
图 １ 图 ２分别 为 Ｔｒ Ｓ ＲＸ卫 星 的外观 及 内部构造 图。可 以看 出 ， 星 主要 由 Ｓ Ｒ天 和 ｅｒ Ａ — ａ 卫 Ａ 线（ 长约 ４ ８ 、Ａ ．ｍ）Ｓ Ｒ单 元 、 进模 块 、 线 装 置 、 作 用 轮 、 推 总 反 ｘ波 段 下 行 数据 传 输 天 线 以及 太 阳能 电池 阵发 电机 等部分 组成 。 Ｔｒ Ｓ Ｒ Ｘ在距地 ５４ ｍ 高 的极地轨 道 上 围绕 地球 运转 ， 收集 高质 量 的 ｘ波段 雷 达数 ｅｒ Ａ — ａ ｋ １ 可
211198598_基于TerraSAR-X
第48卷第2期2023年3月㊀林㊀业㊀调㊀查㊀规㊀划Forest Inventory and PlanningVol.48㊀No.2Mar.2023doi:10.3969/j.issn.1671-3168.2023.02.003基于TerraSAR-X/TanDEM-X数据的森林高度反演算法研究申晨1,岳彩荣2,章皖秋2,朱腾1,张金兰1(1.广东工贸职业技术学院,广东广州510510;2.西南林业大学,云南昆明650224)摘要:基于RVoG模型,对传统三阶段算法及其改进算法进行研究,以云南省西双版纳州勐腊县为研究区,以TerraSAR-X/TanDEM-X数据为数据源进行森林高度反演算法研究,并结合野外实测数据进行结果验证㊂结果表明,三阶段改进算法对森林高度反演精度优于传统三阶段算法㊂三阶段改进算法对天然林高度反演精度较高;三阶段算法对橡胶林高度反演精度较高㊂关键词:森林高度反演算法;TerraSAR-X/TanDEM-X数据;极化干涉SAR;散射机制;平地效应中图分类号:S758.54;P208;TP31;TP753㊀㊀文献标识码:A㊀㊀文章编号:1671-3168(2023)02-0016-10引文格式:申晨,岳彩荣,章皖秋,等.基于TerraSAR-X/TanDEM-X数据的森林高度反演算法研究[J].林业调查规划,2023,48(2):16-25.doi:10.3969/j.issn.1671-3168.2023.02.003SHEN Chen,YUE Cairong,ZHANG Wanqiu,et al.Forest Height Inversion Algorithm Based on TerraSAR-X/TanDEM-X Data[J].Forest Inventory and Planning,2023,48(2):16-25.doi:10.3969/j.issn.1671-3168.2023.02.003Forest Height Inversion Algorithm Based on TerraSAR-X/TanDEM-X DataSHEN Chen1,YUE Cairong2,ZHANG Wanqiu2,ZHU Teng1,ZHANG Jinlan1(1.Guangdong Polytechnic of Industry and Commerce,Guangzhou510510,China;2.Southwest Forestry University,Kunming650224,China) Abstract:Based on the RVoG model,the three-stage algorithm and its improved algorithm were studied. In Mengla County,Xishuangbanna,Yunnan,TerrasAR-X/tandem-X data were used as data sources to conduct forest height inversion algorithm.Finally,the results were verified by field measurement data. The results showed that the accuracy of the improved three-stage algorithm was better than that of the tra-ditional three-stage algorithm.For natural forest and rubber forest,the improved three-stage algorithm had higher accuracy in height inversion of natural forest.The accuracy of the three-stage algorithm for rubber forest height inversion was high.Key words:forest height inversion algorithm;TerraSAR-X/TanDEM-X data;polarization interferometry SAR;scattering mechanism;flat effect1研究背景及现状1.1研究背景森林冠层高度作为森林垂直结构属性,是用来估计林分蓄积量㊁立地生产力和生物量的重要参数㊂当前,森林冠层高度测量主要有野外人工测量和遥感技术反演两种方式㊂在森林地带,遥感手段反演收稿日期:2021-11-11.资助项目:基于PollInSAR数据的森林高度反演算法研究项目(220000025716).第一作者:申晨(1990-),女,山西潞城人,硕士,助理工程师.主要从事地图制图㊁资源环境遥感研究.Email:514935785@申晨,等:基于TerraSAR -X /TanDEM -X 数据的森林高度反演算法研究森林高度比人工测量具有巨大优势㊂极化干涉合成孔径雷达(PolInSAR)是一种新型的遥感技术,既具有极化合成孔径雷达对植被散射体形状和方向敏感的特性,又具有干涉合成孔径雷达对植被体空间分布和高度敏感的特性,不仅可以提高地形测量精度,同时也获得相关散射机制的物理参数,特别是森林参数为树高的高精度反演以及地物其他各种参数的反演提供了可能,在林业资源管理以及林业调查研究方面具有重大意义㊂1.2研究现状合成孔径雷达(SAR)概念首先于1951年由C.Wiley 等提出,其指出多普勒频移现象可以用逻辑合成一个更大的雷达孔径㊂该理论突破了真实孔径雷达由于不可能无限制增大雷达天线孔径大小带来的方位向分辨率低的缺点[1]㊂随着国外PolIn-SAR 技术的发展,国内很多科研院所开始对该技术进行研究㊂2014年付海强等针对传统PCT 方法中 相干相位-幅度联合反演算法 的缺点,采用RVoG模型,利用改进的非线性迭代算法反演植被高度和地表相位[2]㊂2016年郭胜龙等提出一种基于单基线减缩极化干涉SAR 对森林进行林下地形估计及树高反演的方法[3]㊂2017年龙江平等在相干区域边界优化的基础上,对相位信息和幅度信息进行优化,联合干涉相位和相干幅度信息有效地区分出地表散射和森林体散射,提高了森林高度反演的可靠性㊂2020年陈利军等利用X /L /C 段的干涉数据集在青海典型复杂地形实验区成功获取了DEM 数据产品,证明X 波段干涉测量DEM 的精度与细节质量相对较高[4]㊂2森林高度反演算法研究2.1RVoG 散射模型随机体散射体-地表二层(Random Volume overGround,RVoG)模型是目前极化干涉中应用最广泛的模型之一㊂该模型包含散射体和地面两层结构,如图1所示㊂该模型以两层介质模型描述了厚度为h v 随机分布的森林覆盖模型㊂该模型广泛应用于森林高度提取和消光系数估计[5]㊂在RVoG 模型下,森林高度可理解为随机体的厚度㊂图1㊀RVoG 散射模型Fig.1㊀Scattering model of RVoG㊀㊀RVoG 模型中相关公式为:γc (w ң)=e iΦ0γv +m (w ң)1+m (w ң)γv =II 0ʂγv (ω)I =ʏh v 0exp(ik z z ∗)exp 2σz ∗cos θ0()dz ∗I 0=ʏh v 0exp 2σz ∗cos θ0()dz ∗ìîíïïïïk z =m2πλΔθsin θìîíïïïïïïïïïïï(1)式中:γc (w ң)为随机体散射体-地表二层的总复相干;Φ0为地面干涉相位;m 为给定极化状态下的地体幅度比,有效的地体幅度比参数m (w )可以表征电磁波穿过植被体层时引起的衰减,也可以理解为衰减系数和随机体厚度(森林高度)的函数,相位中心的变化可引起地体幅度比的变化;γv 为体散射相干性,体散射相干性复数值γv 与电磁波的极化状态无关;h v 为森林高度;θ为雷达波的入射角;σ为消光系数;k z 为垂直有效波束,对于单站模式(Monos-tatic)的雷达数据,k z 公式中的m =2,对于双站模式(Bistatic)的雷达数据,k z 公式中的m =1;Δθ为两次观测中入射角的差值;λ为电磁波的波长㊂其中已知量为θ㊁k z ㊁λ㊁Δθ,γc (w ң)由极化干涉基本公式估计获得,未知量为h v 和σ㊂2.2三阶段森林高度反演算法三阶段森林高度反演算法是由Cloude 和Patha-㊃71㊃第2期林业调查规划nassiou 提出的一种经典森林高度反演算法,通过直线拟合㊁地表相位求解㊁森林高度参数估计3个步骤完成㊂第一步:直线拟合将RVoG 模型公式变形可得如下变形公式:γc (w ң)=e iΦ0γv +m (w ң)1+m (w ң)=eiΦ0γv +m (w ң)1+m (w ң)(1-γv )()(2)㊀㊀由公式(2)可知,在复平面内,不同极化状态下的复相干系数呈直线分布㊂然而实际情况下,由于森林地带微气候复杂多变㊁植被自身结构千差万别,复相干系数在复平面内的分布并非一条直线,而是一个长椭圆形的复相干区域,如图2所示㊂图2㊀三阶段森林高度反演算法示意Fig.2㊀Three -stage forest height inversion algorithm㊀㊀因而几乎不存在某一条直线,使不同极化状态下的复相干系数全部在该直线上㊂因此,Cloude 和Pathanassiou 提出利用最小二乘法对不同极化状态下的复相干系数在复平面内进行直线拟合㊂㊀㊀第二步:地表相位求解复相干系数的取值在0~1范围㊂由于地面层的高相干性,可以假设地面的复相干系数模值为1.由公式(2)可知,当地体幅度比m (w )=ɕ时,体散射复相干为模为1的地面相干㊂该地相干恰好为复平面内拟合直线与单位圆的两个交点之一㊂一般认为,占据HV 极化通道的复相干γHV 为植被层的体散射复相干,因此,在复平面内距离γHV 最远的一个交点即为地面层的复相干,该交点对应的相位则为三阶段森林高度反演算法中的地表相位Φ0㊂第三步:森林高度参数估计由RVoG 模型公式(1)可知,复相干系数γv 是消光系数σ和森林高度h v 的函数㊂也就是说,每给定一个σ值和一个h v 值,便可对应一个γv 值㊂通过建立γv ㊁σ㊁h v 之间的二维查找表,可得到不同森林高度和消光系数下对应的复相干系数的理论值㊂正因此,三阶段森林高度反演算法求得的森林高度值皆为整数㊂一般认为,占据HV 极化通道的复相干γHV 为植被层的体散射复相干,对应于地体幅度比m (w )=0情况,将m (w )=0带入公式(1)并变形,可得以下公式:γHV =e iΦ0γv ⇒γv =γHV e i (-Φ0)(3)其中求得的γv 为复相干系数的估计值㊂利用最小差异原则比较二维查找表中复相干系数的理论值γv 和由公式(3)求得的复相干系数的估计值,便可确定对应的消光系数σ和森林高度h v ㊂三阶段森林高度反演算法将复相干系数引入反演模型并考虑了地表相位的生成,是一种科学的森林高度反演算法,三阶段算法流程图见图3所示㊂三阶段方法在实际操作中计算量大㊁运算时间长,地表相位的确定受复相干系数分布以及直线拟合方法的影响较大,当体散射去相干存在误差时,往往会导致反演结果出现误差㊂2.3三阶段森林高度反演算法的改进从相干区域优化㊁体散射复相干优化两个方面提高地表相位的准确程度,实现三阶段算法的改进[6]㊂2.3.1相干区域优化为了扩大相干区域范围,提高地表相位的准确程度,本文提出相干性幅度差最大准则和干涉相位差最大准则,实现相干区域的优化㊂相干性幅度差最大准则是通过寻找不同极化状态下的相干性幅度差,选取相干性幅度差最大时对应的极化状态为最优极化状态[6]㊂即满足公式(4):max γ(ωi )-γ(ωj )=Δr opt (4)式中:ωi ㊁ωj 为极化干涉矢量;γ(ωi )㊁γ(ωj )分别为极化干涉矢量ωi ㊁ωj 所对应的复相干系数;γ(ωi )㊁γ(ωj )分别为极化干涉矢量ωi ㊁ωj 所对应的复相干系数的幅度值㊂㊃81㊃第48卷申晨,等:基于TerraSAR-X/TanDEM-X数据的森林高度反演算法研究图3㊀三阶段算法流程Fig.3㊀Flowchart of the three-stage algorithm ㊀㊀相干性幅度差最大的最优相干估计可以看成是选取相干区域中相干性幅度极大值和极小值作为最优相干性,从而使得相干性幅度差最大㊂干涉相位差最大准则是通过寻找不同极化状态下的干涉相位,该准则的本质与PD算法一致[7],也是通过拉格朗日极值法找到相干区域内不同极化状态复相干系数相位的差值,干涉相位差值最大时对应的极化状态为最优相干性估计㊂即满足公式(5): maxφ(ωi)-φ(ωj)=Δφopt(5)式中:ωi㊁ωj为极化干涉矢量;φ(ωi)㊁φ(ωj)分别为极化干涉矢量ωi㊁ωj所对应的干涉相位㊂干涉相位差最大的最后相干性估计可以有效地分离出不同散射体的相位中心,提高了地表相位的估计精度㊂相干区域边界提取见图4㊂总而言之,相干性幅度差最大准则和干涉相位差最大准则有效地扩充了相干区域范围,使地表相位的估计更加准确,有利于森林植被高度反演精度的提高㊂2.3.2体散射复相干优化由极化干涉先验知识可知,极化通道HV所代表图4㊀相干区域边界提取Fig.4㊀Extraction of the boundary in coherence areas的复相干为植被体散射复相干㊂然而,HV复相干中常含有地表分量,并非纯体散射复相干㊂在森林高度反演中,植被冠层常具有较高的相位中心㊂因此,本文提出以不同极化状态中距离地相位最远的复相干为体散射复相干,代替HV的复相干进行三阶段森林高度反演算法研究㊂三阶段改进算法流程图见图5所示㊂3TerraSAR-X/TanDEM-X数据森林高度反演数据试验与分析3.1研究区概况勐腊县位于西双版纳傣族自治州南部,与老挝接壤㊂全境在北回归线以南,总面积约688510hm2,林地面积占总面积的88.31%㊂林地主要分布在坡度10ʎ~30ʎ㊁海拔500~1000m区域,其中天然林约占林地面积的65.56%,橡胶林约占林地面积的34.44%㊂气候属热带湿润季风气候,全年热量丰富,降水充沛,旱雨两季分明㊂勐腊县是我国热带生态系统和森林植被保存较为完整㊁热区生物资源较为富集地区,境内广泛分布有雨林㊁季雨林㊁常绿阔叶林和竹林等植被类型,被称为热带植物王国㊂3.2研究数据介绍为验证森林高度反演算法对真实数据的估计结果,采用德国航空航天中心提供的TerraSAR-X和㊃91㊃第2期林业调查规划图5㊀三阶段改进算法流程图Fig.5㊀Flowchart of the improved three-stage algorithm TanDEM-X卫星的双星全极化SAR影像进行试验㊂其中主轨道为TanDEM-X平台,副轨道为TerraSAR-X平台,拍摄时两卫星平台分别对同一地面发射电磁波,并分别接收各自的回波,具体参数如表1所示㊂由于主辅雷达平台的拍摄时间仅差10秒,可忽略时间影响㊂表1㊀研究数据参数Tab.1㊀Parameters of study data全极化Bistatic2015年12月13431.86~33.34 2.1㊀㊀分别读取主辅影像的单视复数据后生成主辅影像的Pauli彩色合成图(图6)㊂3.3数据预处理对SAR图像的预处理主要包括:主辅影像配准㊁去地平效应㊁生成极化干涉矩阵T6㊁计算复相干系数㊁计算垂直有效波数㊂预处理流程见图7所示㊂图6㊀主辅影像Pauli合成图Fig.6㊀Pauli-basis composite image of masterand salve image图7㊀预处理流程Fig.7㊀Flowchart of preprocessing3.3.1极化SAR图像滤波高分辨率SAR图像更容易受到斑点噪声的影响,因此有必要在解译前进行滤波处理㊂从滤波效果来看,增强Lee滤波可有效抑制斑点噪声,改善SAR图像的视觉效果㊂通常选择3ˑ3的窗口能比较好地保持地物边界信息等细节特征㊂滤波前后的干涉图如图8所示㊂㊃02㊃第48卷申晨,等:基于TerraSAR-X/TanDEM-X数据的森林高度反演算法研究图8㊀滤波前后干涉图Fig.8㊀Interference patterns of HH beforeand after filter3.3.2去平地效应平地效应是由于受InSAR系统空间几何关系的影响,地面上高程相同的点本应在干涉图上保持不变的相位差发生了变化,造成干涉图相位一定程度的偏移[8]㊂平地效应的存在使干涉结果不能真实地反映地面高程,且干涉条纹过密会给相位解缠带来很大困难㊂平地相位目前主要有两种算法:(1)基于卫星轨道数据计算多普勒方程㊁斜距方程和椭圆方程得到平地相位;(2)监测干涉相位频谱中占优势的条纹频率计算平地相位[9]㊂本实验采用去地平方法为干涉相位频谱来除去平地效应,计算得到的平地相位如图9所示㊂图9㊀平地相位Fig.9㊀Flat earth phase ㊀㊀以HH极化通道为例,通过去平地,平地效应引起地面条纹过密的现象得到有效抑制(图10)㊂图10㊀HH极化干涉图Fig.10㊀Interference patterns of HH polarimetric channel 3.3.3垂直有效波数垂直有效波数的计算对于树高估计是必须的参数,数值差异与主辅影像的入射角之差有关㊂能否准确计算垂直有效波数k z,直接关系到后续树高计算结果的准确程度㊂根据k z公式进行运算得到垂直有效波数k z的文件㊂3.3.4极化干涉数据地理编码TerraSAR-X/TanDEM-X原始数据生成的SLC㊁复相干系数以及最后的森林高度数据均无坐标系统㊂不仅如此,SLC数据与真实地貌呈东西镜像,且距离向与方位向的分辨率不同㊂因此,为了方便森林高度结果与光学影像的匹配,提出极化干涉数据地理编码,使雷达数据处理结果具有光学影像的坐标系统㊂3.3.5复相干系数的计算及优化复相干系数包含了极化干涉SAR数据的相位信息和幅度信息,通常情况下,复相干系数以实部(real)和虚部(imaginary)两部分存储,共4个字节㊂极化干涉复相干系数的计算及优化会直接影响到森林高度反演结果的准确性和可靠性㊂不同极化状态下的复相干系数如图11所示,图中像元的灰度值为该点模值的自然度数,其明暗程度也反映了复相干系数模值的大小㊂按不同林地种类,天然林和橡胶林的复相干系数有微小差异㊂将编码后的复相干系数文件与森林资源二类调查数据进行叠加,提取天然林和橡胶林所对应的复相干系数,进行数据统计分析(表2)㊂㊃12㊃第2期林业调查规划图11㊀不同极化状态下的复相干系数Fig.11㊀Complex coherence coefficient in different polarization㊀㊀由表2可知,对于天然林而言,与其他复相干相比,opt1复相干对天然林的相干性较强,复相干系数最大值为0.989654,最小值为0.317398,均值为0.854106,标准差为0.121,复相干系数反演结果较为稳定㊂对于橡胶林而言,与其他复相干相比,opt1复相干对橡胶林的相干性较强,复相干系数最大值为0.991536,最小值为0.417947,均值为0.927743,标准差为0.057,复相干系数反演结果较为稳定; opt3复相干对橡胶林的相干性较弱,复相干系数最大值为0.669315,最小值为0.002098,均值为0.182702,标准差为0.6293,复相干系数反演结果稳定性较差㊂3.4森林高度反演算法结果分析与讨论利用三阶段森林高度反演算法㊁改进的三阶段森林高度反演算法进行森林高度反演,最终得到如下结果㊂1)树高反演结果的数理统计分析分别对三阶段算法㊁三阶段改进算法估测结果的森林区域进行结果统计,得到6种算法的树高直方图(图12)㊂由图12可知,三阶段算法反演的森林高度主要集中在10~20m范围,在树高约5m出现峰值,在树高约10m处出现大的峰值,相对于其他5种方法,三阶段森林高度反演算法的结果区间范围较小㊂即树高值的分布较为集中,但是数值偏低;三阶段改进算法反演的森林高度主要集中在8~50m范围,树高峰值出现在28m处,与传统的三阶段算法相比,树高反演值的区间范围变大,数值大小整体均有提高(表3)㊂㊃22㊃第48卷申晨,等:基于TerraSAR -X /TanDEM -X 数据的森林高度反演算法研究表2㊀不同林种的相干系数Tab.2㊀Correlation coefficient in different forest categoryHV 0.0013260.7188520.1862440.62800.0003370.7582820.1922420.6244VV0.0013780.9757000.7054640.32860.0161760.9838320.8559870.1484HH +VV 0.0003570.7745290.2041910.62820.0018200.7158390.2013320.6200HH -VV 0.0001860.7984360.2049680.62680.0017530.7080080.1978990.6203HV +VH0.0013260.7188520.1862440.62800.0015720.7037020.1937090.6230opt10.3173980.9896540.8541060.12100.4179470.9915360.9277430.0570opt20.0946280.9526570.5791460.27460.1308240.9468640.6913560.1946opt30.0004640.7321280.1649360.63540.0020980.6693150.1827020.6293PDhigh 0.0000150.9999990.1537850.79430.0001220.9999990.1462850.8050Pdlow 0.0000350.9999990.1533860.71830.0000460.9999990.1436610.7006MaxMag 0.1773730.9999990.7284760.19660.2733780.9999990.8225880.1207MaxPha 0.0288340.9999990.6759070.26630.0450000.9999990.7888340.1572MinMag 0.0086450.9999990.6289520.30160.0450000.9999990.7569180.1797MinPha 0.0000800.9999990.6382660.31210.0891350.9999990.7898600.1566LR 0.0003570.7745290.2041910.62820.0018200.7158390.2013320.6200LL 0.0009650.7886760.2077130.62350.0006100.6948910.2025640.6331RR0.0008870.7849780.2070280.62070.0022970.7156240.2017960.6175图12㊀不同算法的树高直方图Fig.12㊀Histograms of forest height of different algorithm表3㊀不同算法的树高统计Tab.3㊀Forest height estimation indifferent algorithmm ㊀三阶段算法05010.90414.275三阶段改进算法5010.85814.9283.5森林高度反演算法结果验证3.5.1野外实测数据调查为验证TerraSAR -X /TanDEM -X 森林高度反演算法的准确程度,于2016年10月底在云南省西双版纳州勐腊县进行野外数据采集㊂利用激光测距仪㊁微波测距仪㊁标杆㊁卷尺㊁GPS㊁角规等仪器测量了82个森林高度样点㊂其中,天然林19个样点,橡㊃32㊃第2期林业调查规划胶林55个样点,其他林种8个样点㊂野外数据采集主要参数有样地经纬度㊁海拔㊁坡度㊁优势树种㊁郁闭度㊁树高㊁平均冠层高等㊂样点广泛分布于勐腊县城周边林区(图13)㊂图13㊀野外实测样点分布Fig.13㊀Distribution of field sample plots3.5.2树高反演结果验证将采用三阶段算法㊁三阶段改进算法的森林高度反演结果与实测数据进行对比㊂在去除异常值后,制作树高反演散点图,结果如图14,图15,表4所示㊂由树高反演散点图和反演结果精度表可知,三阶段改进算法相关性较强,相关系数R 为0.643,均方根误差(RMSE)为5.53,相较于三阶段算法有明显的改进㊂分别对天然林和橡胶林进行实测数据验证,结果如表5所示㊂从天然林反演结果来看,三阶段改进算法的相关性相对较强,相关系数R 为0.571,均方根误差(RMSE)为7.33,说明三阶段改进算法对地相位和冠层相位的估计有较好的改善㊂图14㊀三阶段算法树高反演散点图Fig.14㊀Tree height inversion scatter plot of thethree -stage algorithm图15㊀三阶段改进算法树高反演散点图Fig.15㊀Tree height inversion scatter plot of theimproved three -stage algorithm表4㊀森林树高反演结果精度Tab.4㊀Precision of tree height inversion results三阶段算法0.483三阶段改进算法0.6435.53表5㊀天然林和橡胶林反演结果精度Tab.5㊀Precision of inversion results in naturalforests and rubber plantation三阶段算法0.3369.200.556 5.61三阶段改进算法0.5717.330.4535.17㊃42㊃第48卷申晨,等:基于TerraSAR-X/TanDEM-X数据的森林高度反演算法研究㊀㊀从橡胶林反演结果来看,三阶段算法的相关性较强,相关系数R为0.556,均方根误差(RMSE)为5.61,三阶段改进算法的相关性次之,相关系数R 为0.453,均方根误差(RMSE)为5.17㊂4结㊀论以RVOG模型为基础分析森林高度反演的三阶段算法㊁三阶段改进算法基本原理,利用TerraSAR-X/TanDEM-X数据进行森林高度反演,并结合野外数据进行结果验证㊂结果表明,三阶段改进算法对森林高度反演的精度优于传统三阶段算法㊂对于天然林而言,三阶段改进算法反演精度较高,相关系数R为0.571,均方根误差(RMSE)为7.33;对于橡胶林而言,三阶段算法反演精度较高,相关系数R为0.556,均方根误差(RMSE)5.56㊂原因有以下几点:(1)TerraSAR-X/TanDEM-X数据为X波段的极化干涉数据,穿透性较弱,因而对森林冠层相位的反演较好,对地表相位的反演效果不佳;(2)天然林结构错综复杂,树种丰富,林层通常为复层林,使电磁波的穿透和森林高度反演增加了难度,橡胶林林冠整齐㊁植株株距较大㊁林下干净,导致电磁波更容易穿透到达地面,为地表相位和冠层相位的反演提供了便利;(3)勐腊县雨林地区地势陡峭,地形和坡度对反演结果的精度造成一定影响㊂参考文献:[1]鞠晓洁.基于Multigen Vega的SAR图像生成仿真系统研究[D].青岛:中国海洋大学,2010.[2]付海强,汪长城,朱建军,等.一种改进的PolInSAR PCT方法反演植被垂直结构[J].测绘工程,2014(11): 56-61,66.[3]郭胜龙,李洋,尹嫱,等.基于简缩极化干涉SAR数据的森林垂直参数反演[J].电子与信息学报,2016(1): 71-79.[4]陈利军,张瑞,郭旭东,等.多源星载SAR地形干涉测量精度分析[J].测绘通报,2020(7):13-17. [5]PAPATHANASSIOU K P,CLOUDE S R.Single-baseline polarimetric SAR interferometry[J].IEEE Transactions on Geoscience and Remote Sensing,2001,39(11):2352-2363.[6]谈璐璐,杨立波,杨汝良.合成孔径雷达简缩极化干涉数据的植被高度反演技术研究[J].电子与信息学报, 2010(12):2814-2819.[7]龙江平,刘峰,段祝庚.联合干涉相位和相干性幅度的极化干涉SAR最优相干性估计[J].测绘学报,2017 (1):73-82.[8]张振宇.合成孔径雷达三维成像技术研究[D].南京:南京理工大学,2009.[9]TIGHE M L,KING D,BALZTER H,et parison of X/C-HH InSAR and L-PolInSAR for Canopy Height Esti-mation in a Lodgepole Pine Forest[J].Proceedings of the Fourth International Workshop on Science&Applications of Sar Polarimetry&Polarimetric Interferometry Poiinsar, 2009.责任编辑:许易琦㊃52㊃第2期。
国外相关星载SAR及InSAR卫星的系统参数
本报告主要调研了国际上到目前为止所存在的一些星载InSAR系统的发展情况，总结了各系统的一些技术指标及参数选择。
以下调研系统中，除了TanDEM-X干涉系统之外，其他的星载SAR系统都不是用来专门进行干涉测量使用的，它们基本的任务还是实现二维高分辨成像，因此大多采用的是重复轨道干涉测量模式。
在进行干涉测量之前，首先要估算此次测量的基线数值，如果不满足要求，此次测量数值就不会采用，因此，对于重复轨道干涉测量的基线实际上是针对需要的测高精度筛选出来的。
1、美国Seasat系统1978年6月，美国国家航空航天局发射了海洋卫星（SeaSat），在卫星上首次装在了合km的面积进行了测绘，该卫星在空间飞行100天，采用的成孔径雷达，对地球表面1亿2是重复轨道干涉模式，首次从空间获得地球表面雷达干涉测量数据。
ERS-1和ERS-2雷达卫星为欧洲空间局分别于1991年和1995年发射，携带有多种有效载荷，包括侧视合成孔径雷达和风向散射计等装置。
ERS-1和ERS-2雷达卫星构成对同一地面访问时间相差一天的星对，使得两次取得的SAR数据之间的相干性得到了一定保障，采用太阳同步晨昏轨道，该系统采用的是重复轨道干涉模式，卫星编队形式为跟飞。
获得。
3、日本JERS-1系统JERS-1雷达系统是日本于1992年发射升空的，采用太阳同步晨昏轨道，该卫星采用了重复轨道干涉模式，但其轨道控制方式不太理想，在交轨方向的基线分量不如日本之后发射的ALOS卫星。
表3中的基线长度是对JERS-1持续观测四年（1993年—1994年）期间的基线变化范围。
雷达卫星Radarsat除了有一个地面卫星数据接收站外，卫星上还载有磁带记录器，可覆盖全球。
该卫星除陆地及海洋应用外，其还肩负两个方面的重要任务：一是对南极大陆提供第一个完全的高分辨率卫星覆盖，二是对全球产生多次卫星覆盖。
Radarsat雷达卫星由加拿大于1995年11月4日发射，具有7种模式、25种波束及不同入射角，因而具有多种分辨率、不同幅宽和多种信息特征，使用于全球环境、土地利用和自然资源监测等。
微波遥感_常用SAR数据产品
多视（Multilooking）处理
单视：处理整个多普勒带宽 ，即指只使用一段合成孔 径长度所成的SAR图像，没有和其它SAR图像进行叠 加，通常所用的合成孔径长度比较长，这样方位分辨 率比较高 多视：把整个多普勒分成几段多普勒数据，即指将整 个有效合成孔径长度分成多段分别对同一场景进行成 像，然后将所得的图像求和叠加得到一副SAR图像， 由几幅SAR图像叠加，就是几视的SAR图像 目的：可以提高SAR图像的信噪比，有效抑制斑点噪 声，但代价是会降低方位分辨率
TerraSAR COSMOSkymed AIRSAR
Microwave Remote Sensing
1 Envisat ASAR
在ENVISAT卫星上载有多个传感器，分别对陆地、海洋 、大气进行观测，其中最主要的就是ASAR（Advanced Synthetic Aperture Radar）的合成孔径雷达传感器
数据组织方式: 数据＋头文件 数据文件: 只存储二进制数据，比如一幅图像有m 行，n列，每个像元用一个浮点数表示，那么这 个数据文件就只存m*n个浮点数 头文件: 存储文本信息，又可称之为元数据。里 面包含可能有图像大小，数据组织方式，轨道信 息，坐标信息等等
Microwave Remote Sensing
（2）Level 1B产品
原始数据及Level 0产品可以通过成像算法并利用标定数据生成Level 1B产品。Level 1B产品可 以是单景产品或条带产品。单景产品是按 “景”定购，而条带产品则包含了整个数据段（segment ）的图象，每个条带产品的最大时间可以达到10分钟。 Level 1B产品分为三类：Precision Image、Single Look Complex Image 和 Medium Resolution
SAR雷达卫星影像数据的基本知识
SAR雷达卫星影像数据的基本知识SAR雷达卫星影像数据的基本知识⽤⼀个⼩天线作为单个辐射单元，将此单元沿⼀直线不断移动，在不同位置上接收同⼀地物的回波信号并进⾏相关解调压缩处理。
⼀个⼩天线通过“运动”⽅式就合成⼀个等效“⼤天线”，这样可以得到较⾼的⽅位向分辨率，同时⽅位向分辨率与距离⽆关，这样SAR就可以安装在卫星平台上⽽可以获取较⾼分辨率的SAR图像。
图1 SAR成像原理⽰意图1、⼏个参重要参数为了更好的理解SAR和SAR图像，需要知道⼏个重要的参数。
分辨率SAR图像分辨率包括距离向分辨率（Range Resolution）和⽅位向分辨率（Azimuth Resolution）。
图2 距离向和⽅位向⽰意图距离向分辨率（Range Resolution）垂直飞⾏⽅向上的分辨率，也就是侧视⽅向上的分辨率。
距离向分辨率与雷达系统发射的脉冲信号相关，与脉冲持续时间成正⽐：Res( r) = c*τ/2其中c为光速，τ为脉冲持续时间。
⽅位向分辨率（Azimuth Resolution）沿飞⾏⽅向上的分辨率，也称沿迹分辨率。
如下为推算过程：真实波束宽度：β= λ/ D真实分辨率：ΔL = β*R = Ls (合成孔径长度)合成波束宽度βs = λ /(2* Ls) = D / (2* R)合成分辨率ΔLs = βs* R = D / 2其中λ为波长，D为雷达孔径，R为天线与物体的距离。
从这个公式中可以看到，SAR系统使⽤⼩尺⼨的天线也能得到⾼⽅位向分辨率，⽽且与斜距离⽆关（就是与遥感平台⾼度⽆关）。
图3 ⽅位向分辨率⽰意图极化⽅式雷达发射的能量脉冲的电场⽮量，可以在垂直或⽔平⾯内被偏振。
⽆论哪个波长，雷达信号可以传送⽔平（H）或者垂直（V）电场⽮量。
接收⽔平（H）或者垂直（V）或者两者的返回信号。
雷达遥感系统常⽤四种极化⽅式———HH、VV、HV、VH。
前两者为同向极化，后两者为异向(交叉)极化。
极化是微波的⼀个突出特点，极化⽅式不同返回的图像信息也不同。
SAR遥感数据产品分级规范应用案例
附录A（资料性附录）SAR遥感数据产品分级规范应用案例A.1TerraSAR-X数据产品分级案例TerraSAR-X卫星为德国研制的一颗高分辨率雷达卫星，携带一颗高频率X波段合成孔径雷达传感器，可以聚束式、条带式和推扫式3种模式成像，并拥有多种极化方式。
其数据产品分级案例见表A.1。
表A.1TerraSAR-X数据产品分级方案与本规范分级对照表SAR遥感数据产品分级规范解释TerraSAR-X数据产品分级解释0级L0原始标准景产品————1级L1A单视复数据产品SSC单视斜距复影像数据产品L1B多视功率产品MGD多视地距探测产品：经过噪声抑制和投影到地距的产品L1C多视功率增强产品————2级L2几何校正产品GEC 地理编码椭球校正产品：用WGS84椭球对影像进行了通用横轴墨卡托投影（UTM）或通用极球面投影（UPS）后的产品3级L3几何精校正产品————4级L4正射校正产品EEC进行地理编码且使用SRTM DEM进行地形起伏纠正的产品ORI SAR高精度正射影像：使用高精度DEM进行地形纠正的产品5级L5FUS区域融合产品MC SAR镶嵌产品OI SAR定位影像：镶嵌产品的子集或者一景影像的一个区域ADM SAR升降轨融合产品L5DTP_DEM数字地形产品数字高程模型产品————L5DTP_DSM数字地表模型产品L5DTP_DLG数字线划图产品L5DTP_DEF地表形变产品L5SCP_POL地表覆盖产品极化特征产品表A.1(续)SAR遥感数据产品分级规范解释TerraSAR-X数据产品分级解释5级L5SCP_CLA 地表覆盖产品地物分类产品————6级L6专题产品OM SAR正射图：具有地图结构和图例的正射制图产品CDM SAR变化监测图SUB SAR地表沉陷图A.2COSMO-SkyMed数据产品分级案例COSMO-SkyMed系统是一个由意大利航天局和意大利国防部共同研发的4颗雷达卫星组成的星座，工作于X波段，提供了具有全球覆盖能力、适应各种气候的日夜获取能力及高分辨率、高精度、高干涉极化测量能力的高效便利的产品服务。
TerraSAR-X卫星成像模式介绍
北京揽宇方圆信息技术有限公司TerraSAR-X卫星成像模式介绍一、卫星概要：TerraSAR-X卫星为德国研制的一颗高分辨率雷达卫星，携带一颗高频率的X波段合成孔径雷达传感器，可以聚束式、条带式和推扫式3种模式成像，并拥有多种极化方式。
可全天时、全天候地获取用户要求的任一成像区域的高分辨率影像。
TanDEM-X于2010年6月21日成功发射，这两颗卫星在3年内将反复扫描整个地球表面，最终绘制出高精度的３Ｄ地球数字模型。
二、TerraSAR-X卫星主要特点1、多分辨率（1m/3m/18.5m）和覆盖区域：对于特定目标区域采用高分辨率，对大面积覆盖采用中等分辨率2、任何其他的商业星载传感器都无法比拟的几何精度3、极高的辐射精度4、不受天气影响，对地球上的任何地点，重访周期最长2.5天（95%的地区可达到2天重访）5、独特的敏捷性（成像和极化模式的快速切换）6、可应用于多种领域：防灾：洪水监测、地震监测、火山监测、滑坡、溢油监测农业：农作物分类、农作物长势监测及估产制图：地物提取、变化监测、地图制图林业：森林分类、林业资源评估与监测水文：土壤湿度监测、沼泽地识别海洋：海冰类型识别、冰川监测、极地监测、海洋表面监测地质：岩性构造、城市地面沉降、矿区沉降监测北京揽宇方圆信息技术有限公司是国内的领先遥感卫星数据机构，遥感行业的国家高新技术企业，整合全球200多颗遥感卫星数据资源，遥感卫星影像数据贯穿中国1960年至今的所有商业卫星影像数据，是中国遥感卫星数据资源最多的专业遥感卫星数据服务机构，提供多尺度、多分辨率、全覆盖的遥感卫星影像数据服务，最大限度的保证了遥感影像数据获取的及时性和完整性。
分发不同性能、技术应用上可以互补的多种卫星影像，包括光学、雷达卫星影像、历史遥感影像等各种卫星数据服务，各种专业应用目的的图像处理、解译、顾问服务以及基于卫星影像的各种解决方案等。
公司拥有完全自主知识产权、高性能、满足大规模遥感数据集中处理的空间大数据管理与服务系统。
微波遥感-2
Rr
c
2 cos
方位分辨率：
Ls Ra
D
R
采用合成孔径技术合成后的天线孔径为真实孔径雷达的 波瓣宽度（方位分辨率）即
Rs
Ls
RD
由于电磁波双程位移，最终的方位分辨率还可提高一倍。
D Rs 2
合成孔径雷达得到的原 始数据还不能叫做图像, 只是一组包含强度、位 相、极化、时间延迟和 频移等信息的大矩阵,叫 做(原始)信号数据(Raw Signal Data)。从信号数 据到图像产品,要经过复 杂的步骤。
d sin dd
二、雷达方程
雷达波束是以天线为中心的球面波。包括： 发射-发射增益-传输吸收-到达-回波吸收-接收 增益-接收总功率。
第二节 成像微波传感器
一、微波辐射计
中文名称：
微波辐射计
英文名称：
microwave radiometer
用以收集和测量地物发射来的微波辐射通量的被动式微波遥感探测仪。 测绘学（一级学科）；摄影测量与遥感学（二级学科） 能定量探测目标物的低电平微波辐射的高灵敏度的接收装置。 地理学（一级学科）；遥感应用（二级学科） 主要用于测量海面辐射强度的微波遥感器。 海洋科技（一级学科）；海洋技术（二级学科）；海洋遥感（三级学科）
第三节 天线、雷达方程和灰度方程
一、雷达天线 天线是将发射机产生的电磁波发射出去，或者 接收地物反射回来的电磁波，并传送给接收机 的重要中间环节。 各向同性的天线（各方向辐射能量相同）是不 存在的。天线立体方向图表示天线能量的分布 空间。
平面天线方向图
TerraSAR-X卫星聚束模式数据干涉处理研究和分析
Ｔ ￣ Ｓ Ｒ卫星采用精确的轨道控制技术， ｅ ａＡ 实时双频 Ｇ Ｓ Ｐ 来获取坐 标点息，从而确保在获取地面信 息的同时用同一个角度对该地区进行 拍摄 ， 证 了同一地 区观测 时卫星斜 视角 基本一 致 。 保 重复 轨道干 涉处理 时干涉影像对的多普勒中心频率差异一般小于 １０ ， ０ Ｈｚ基线距小于 ２ ０ ４ 米 ，特别适合影像的配准和干涉处理。以前的卫星干涉处理方法对 图 ３ 干涉 相位 图 Ｔ ￣ Ｓ Ｒ Ｘ卫星干涉也同样适用 ， ｅａＡ — 在本文 中就不再赘述， 只探讨 Ｔ ｒ ｅ’ 总 结 ｒＳ Ｒ Ｘ卫 星数据 干涉过 程中 区别于其 它卫 星数据 的重采样 过 程 。滑 ａＡ — Ｔｒ Ｓ Ｒ Ｘ 卫星 是 目前 世 界 上最 优 秀 的民用 Ｓ Ｒ卫 星之 一 ， ｅ ａＡ — ｒ Ａ 随 动聚束模式在成像过程中方位向多普勒 中心变化较大 ，在干涉处理的 着 Ｔ ｎ Ｅ Ｘ系统组成 ， ａ Ｄ Ｍ— 该卫星将 【ｌ １以根据 ＩＳ Ｉ ｎＡＲ技术获取世 界上精 重采样过程需要做必要的算法修改。 度最高 ， 分辨率也最高的 Ｄ Ｍ， Ｅ 而聚束模式影像的成功干涉 ， 改变以前 聚束模 式数据 和条带 模式数 据一 样 ， 为零 多普勒 坐标 系 ， 不 同的 Ｓ Ｒ卫星 仅仅 能从 宏观上 检测地 面形变 的技 术现状 ，可 以直 接根据 干 而 Ａ 是聚束模式的多普勒中心频率用影像的第一行开始到最后一行结束多 涉处理结果去监测地物的细节变化信息 。目前聚束模式干涉已在我国 普勒变化范围从 ＋ ００ ｚ 一 ０ ０ ６ ０ Ｈ 到 ６０ ，而在条带模式数据多普勒中心仅 矿 沉降和输电线路的走廊沉降监测中得到了成功应用。 仅是随着列的变化而稍微变化 。聚束模式影像干涉处理时这一差异必 参 考文献
引言
TerraSAR-X卫星数据单极化HS模式详解
单极化 HS 模式
类型 成像模式 产品类型 几何投影 极化模式 分辨率模式 极化数据波段数 极化方式 数据接收范围 仰角波束数量 全性能范围 距离向幅宽 方位向幅宽 绝对辐射精度 相对辐射精度 NESZ DTAR PSLR ISLR 入射角记录 每像素占比特 入射角 （20°-45°） 斜距分辨率（m） 地距分辨率（m） 方位向分辨率（m） 距离向像素间隔 （m） 方位向像素间隔 （m） 有效视数 像素定位精度 辐射分辨率 产品尺寸（MB） 2.1 50 3.1 313 1.5 50 1.5 50 275 2.9 1.2 6.5 7.4 2.0 2 0.8 2 2 0.8 3.2 3.2 2 1.3 1.3 0.8 4.8 4.8 2 3.3 3.3 2 1.1 0.9 20 55 1 16 20 55 1.0 -19 dB <-17 dB -20 dB -13 dB 0 32 SE 1 {HH，VV} 15°— 60° 124 20°— 55° 10km 5km 地距影像 {MGD，GEC，EEC} S RE {MGD，GEC， EEC}_SE_HS_S {MGD，GEC， EEC}_RE_HS_S HS 复数 SSC SSC_HS_S
TerraSAR雷达卫星和WorldDEM产品
风险
能见度非常低的情况下 地形或者是障碍带来的危险
Landing zones Ortho Imagery Topographic DEM potential landing zone
Avionic Geoinformation Layers
Synthetic Vision
Landing Zone Detection
Improved NESZ
14
内容
1. WorldDEM™项目 2. WorldDEM™ 产品 3. WorldDEM™销售
15
1
WorldDEMTM
全球统一标准的DEM产品
SAR Capabilities & WorldDEM - Airbus DS Industry briefing @ Geoscience Australia
20
全球高程信息为油汽田和矿井服务
WorldDEMTM 是全球无接缝的DEM产品，无论是多么便宜或者是困难区域
第一阶段 可行性研究 第二阶段 勘探 第三阶段 开发 第四阶段 开采
1. 2. 3.
根据AOI建立地理数据库 基于地理数据库评估油藏潜力 建议进一步管理措施
1. 2. 3. 4.
现场侦察规划勘探活动 地质/地球物理（如地震学）现场 测量/数据处理与解释 测试钻井和井下地球物理测量 投资成本/环境影响评估，决策制 定蓄水池
下一阶段产品 地面高程模型(DTM) 地表裸露的土地，人为建筑物和植被将被移除
DSMDLeabharlann MWorldDEMcore
WorldDEMTM
WorldDEMcore
22
WorldDEMTM
WorldDEMcore 数据包
TerraSAR-X高分辨率雷达卫星数据介绍
hsatellite
带宽Bandwidth 数据收集范围Data collection range 全效率范围Full performance range
高分辨率SpotLight模式
High Resolution SpotLight Mode
10 km x 10 km 1-4 m 取决于入射角和多极化方式 单极化 (VV or HH) 双极化 (HH/VV) 150 and 300 MHz 15° to 60° 20° to 55° 5 x 10 km
StripMap 模式StripMap mode
30km幅宽30 km swath 3m分辨率3 m resolution
高分辨率&SpotLight模式
High Resolution & SpotLightmode
5/10 km x 10 km 1m分辨率1 m resolution
双重天线接收模式Dual receive antenna mode
卫星轨道
Satellite orbit
hsatellite
空间分辨率Spatial resolution 极化方式Polarization
星下轨迹
Nadir track
Θ1=2
0o
数据收集范围Data collection range 全效率区域Full performance range
30 km (单极化) 15 km (双极化) 50 km 4,200 km technical, 1.650 km acc. to product spec. 3 m （单极化）, 6 m（双极化） Single (VV or HH) Dual (HH/VV, HH/HV, VV/VH) Quad (HH/VV/HV/VH) 15° to 60° 20° to 45°
TerraSAR-X雷达卫星数据介绍
北京揽宇方圆信息技术有限公司
TerraSAR-X雷达卫星数据介绍
一、卫星概要：
TerraSAR-X卫星为德国研制的一颗高分辨率雷达卫星，携带一颗高频率的
X波段合成孔径雷达传感器，可以聚束式、条带式和推扫式3种模式成像，
并拥有多种极化方式。
可全天时、全天候地获取用户要求的任一成像区域的
高分辨率影像。
TanDEM-X于2010年6月21日成功发射，这两颗卫星
在3年内将反复扫描整个地球表面，最终绘制出高精度的３Ｄ地球数字模
型。
二、TerraSAR-X卫星主要特点
1、多分辨率（1m/3m/18.5m）和覆盖区域：对于特定目标区域采用高
分辨率，对大面积覆盖采用中等分辨率
2、任何其他的商业星载传感器都无法比拟的几何精度
3、极高的辐射精度
4、不受天气影响，对地球上的任何地点，重访周期最长2.5天（95%的
地区可达到2天重访）
5、独特的敏捷性（成像和极化模式的快速切换）
6、可应用于多种领域：
防灾：洪水监测、地震监测、火山监测、滑坡、溢油监测
农业：农作物分类、农作物长势监测及估产
制图：地物提取、变化监测、地图制图
林业：森林分类、林业资源评估与监测
水文：土壤湿度监测、沼泽地识别
海洋：海冰类型识别、冰川监测、极地监测、海洋表面监测
地质：岩性构造、城市地面沉降、矿区沉降监测
公司地址：北京市丰台区南三环万柳桥宝隆大厦1－1626电话：4006019091010－57113949。
高分三号卫星雷达极化方式和分辨率
北京揽宇方圆信息技术有限公司高分三号卫星雷达极化方式和分辨率高分三号卫星搭载的传感器是C频段多极化合成孔径雷达，是迄今为止世界上成像模式最多的星载合成孔径雷达，该雷达具有全极化电磁波收发功能，并涵盖了诸如条带、聚束、扫描等12种成像模式（表1）。
空间分辨率从1 m到500 m，幅宽10 km到650 km。
不仅能够用于大范围资源环境及生态普查，还能够清晰地分辨出陆地土地覆盖类型和海面目标，现了既可探地，又可观海，到“一星多用”的效果。
下面为高分三号卫星的12中成像模式的相关介绍。
表1 成像模式聚束模式：观测海面溢油现场尺度、岛礁的位置、面积、建筑物、地上交通线、重要水利工程、泥石流；进行城市规划监测、风景名胜区监测、经济普查与经济活动调查、统计重大项目投资监测、人口普查与城市住户调查、边境反恐监测、全球敏感区域监测、城市规划编制。
精细条带1模式：进行海冰表面拓扑、冰山、海面溢油、海岸尺度、洪涝、洪涝淹没范围、农牧林用地、防洪设施、生态格局动态、毒品原植物监测。
精细条带2模式：进行冰凌或海冰、堰塞水体、森林资源相关地类识别、农业普查、海岸带变迁、浅海地形、内波波长、波向、波速、振幅、深度监测。
标准条带模式：进行积雪范围、干旱范围、海冰监测、湖泊藻类、海洋藻类、海冰类型、冰区航道、海面溢油区域尺度、锋面和涡的位置尺度、舰船、海浪监测。
窄幅扫描模式：进行旱情、近海海冰、水体监测。
宽幅扫描模式：进行海冰外缘线、雪覆盖、雪深、极冰监测。
全极化条带1模式：进行农业普查统计、城市建设专题信息提取。
全极化条带2模式：进行积雪范围、干旱范围、海冰、湖泊藻类、海洋藻类监测。
波成像模式：进行海面风场风速、风向、水体监测、干旱、波长、波高、波向监测。
全球观测模式：进行冰融化阶段、内波、土壤水分、海面溢油、干旱、环境应急、极地冰川监测。
扩展低入射角模式：进行船舶、溢油、海冰、海岸带、海洋维权、海洋环境保护和防灾救灾监测。
TerraSAR-X卫星传感器参数大全
地距影像
{MGD，GEC，EEC}
S
SE
RE
1
{HH，VV}
15°— 60°
124
20°— 55°
10km
5km
-19 dB
<-17 dB
-20 dB
-13 dB
1
16
20
55
20
55
3.2
1.3
4.8
3.3
3.2
1.3
4.8
3.3
2
0.8
2
2
2
0.8
2
2
2.9
1.2
6.5
7.4
2.1
3.1
1.5
1.5
133
12
33
复数 SSC
0 32 1.9 2.2 1.4 2.5 2.0 196
5．单极化 SL 模式产品 类型
{MGD，GEC，
{MGD，GEC，
第 4 页·共 7 页
SSC_SL_S
成像模式 产品类型 几何投影 极化模式 分辨率模式 极化数据波段数 极化方式 数据接收范围 仰角波束数量 全性能范围 距离向幅宽 方位向幅宽 绝对辐射精度 相对辐射精度
2.4
1.5
1.5
产品尺寸（MB）
1500
2667
375
667
2．双极化 SM 模式产品
类型
{MGD，GEC，
{MGD，GEC，
EEC}_SE_SM_D
EEC}_RE_SM_D
成像模式
SM
产品类型
地距影像
几何投影
{MGD，GEC，EEC}
极化模式
D
分辨率模式
TerraSAR-X几何校正及辐射定标模型
自动化技术、计算机技术
科学技术与工程
ScienVol.17 No.14 May 2017 © 2017 Sci.Tech.Engrg.
- T e r r a S A R X 几 何 校 正 及 辐 射 定 标 模 型
且 有 多 种 极 化 类 型 ，全 天 时 、全 天 候 的 任 一 成 像 区 域 的 高 分 辨 率 影 像 都 能 够 获 得 。TerraSAR-X卫 星 主
要 特 点 是 多 种 分 辨 率 类 型 、几 何 精 度 超 过 大 多 数 其 他 的 商 业 卫 星 、辐 射 精 度 极 高 、不 容 易 被 天 气 影 响 、
and radiation calibration for T e r r a S A R - X i m a g e [ J ] • Science Technology
and Engineering，2 0 1 7 ，1 7 ( 1 4 ) ：2 2 3 — 2 2 7
S A R 影 像 的 辐 射 定 标 [4] 是 一 个 建 立 图 像 与 地
射 定 标 模 型 [ J ] . 科 学 技 术 与 工 程 ，2 0 1 7 , 1 7 ( 1 4 ) : 2 2 3 — 2 2 7
W a n g I ，Li S o n g ，W e n g Y i n k a n ，et al. A model of geometric correction
第 一 作 者 简 介 :王 吉 （1 9 9 3 — ) , 男 ，硕 士 研 究 生 。研 究 方 向 ：合成
孔 径 雷 达 的 几 何 定 标 和 辐 射 定 标 工 作 。 E-mail: 331395473 @
雷达卫星影像极化方式
北京揽宇方圆信息技术有限公司雷达卫星影像极化方式雷达卫星影像在极化方面，不同的被观测物体对于入射的不同极化波，后向散射不同的极化波。
因此空间遥感可以使用多波段来增加信息含量，也可以用不同的极化来增强，提高识别目标的准确度。
经验表明，对于海洋应用，L波段的HH极化较敏感，而C波段是VV极化比较好；对于低散射率的草地和道路，水平极化使地物之间有较大的差异，所以，地形测绘用的星载SAR都使用水平极化；对粗糙度大于波长的陆地，HH或VV无明显变化。
下图所示，同一目标对于四种不同极化的成像，V表示垂直极化。
经验表明，不同极化下同一地物的回波强弱不同，图像的色调也不一样，增加了识别地物目标的信息。
相同极化（HH，VV）和交叉极化（HV，VH）的信息比较，可以显著地增加雷达图像信息，而且，植被和其他不同地物的极化回波之间的信息差别比不同波段之间的差别更敏感。
所以，多极化工作是SAR卫星发展方向之一。
一、卫星类型(1)光学卫星：worldview1、worldview2、worldview3、worldview4、quickbird、geoeye、ikonos、pleiades、deimos、spot1、kompsat系例、spot2、spot3、spot4、spot5、spot6、spot7、landsat5(tm)、Sentinel-卫星、landsat(etm)、rapideye、alos、kompsat系例卫星、planet卫星、北京二号、高景一号、资源三号、高分一号、高分二号、环境卫星。
(2)雷达卫星：terrasar-x、radarsat-2、alos雷达卫星、高分三号卫星、哨兵卫星(3)侦查卫星：美国锁眼卫星全系例(1960-1980)二、卫星分辨率(1)0.3米：worldview3、worldview4(2)0.4米：worldview3、worldview2、geoeye、kompsat-3A(3)0.5米：worldview3、worldview2、geoeye、worldview1、pleiades(4)0.6米：quickbird、锁眼卫星(5)1米：ikonos、高分二号、kompsat、deimos(6)1.5米：spot6、spot7、锁眼卫星(7)2.5米：spot5、alos、资源三号、高分一号、锁眼卫星(8)5米:spot5、rapideye、锁眼卫星、planet卫星4米(9)10米:spot5、spot4、spot3、spot2、spot1、Sentinel-卫星(10)15米:landsat5(tm)、landsat(etm)三、卫星国籍(1)美国：worldview1、worldview2、worldview3、quickbird、geoeye、ikonos、landsat5(tm)、landsat(etm)、锁眼卫星、planet卫星(2)法国：pleiades、spot1、spot2、spot3、spot4、spot5、spot6(3)中国：资源三号、高分一号、高分二号、高景卫星(4)德国：terrasar-x、rapideye(5)加拿大：radarsat-2四、卫星发射年份(1)1960-1980年：锁眼卫星(0.6米分辨率至10米)(2)1980-1990年：landsat5(tm)、spot1(3)1990-2000年：spot2、spot3、spot4、landsat(etm)、ikonos(4)2000-2010年：quickbird、worldview1、worldview2、spot5、rapideye、radarsat-2、alos(5)2010-：spot6、spot7、资源三号、高分一号、高分二号、worldview3、worldview4、pleiades、高景卫星、planet卫星优势：1：北京揽宇方圆国内老牌卫星数据公司，经营时间久，行业口碑相传，1800个行业用户选择的实力见证。
雷达干涉测量(崔松整理)之欧阳术创编
雷达干涉测量（崔松整理）时间：2021.02.02 创作：欧阳术第二章绪论第三章雷达SAR：使用短天线一段时间内不断收集回波信号，通过信号聚焦处理方法合成一较大的等效天线孔径的雷达。
1.1雷达及雷达遥感发展概况ENVISAT与ERS的SAR传感器相比，Envisat ASAR的优点主要表现在：扫描合成孔径雷达(ScanSAR)可达到500km的幅照宽度；（ERS只有100km）可获得垂直和水平极化信息；（如果发射的是水平极化方式的电磁波，与地物表面发生作用后会使电磁波极化方向产生不同程度的旋转，形成水平和垂直两个分量，用不同极化方式的天线接收，形成ＨＨ和ＨＶ两种极化方式的图像。
若雷达发射的是垂直极化方式的电磁波，同理，会产生ＶＶ和ＶＨ两种极化方式的图像。
）交替极化模式可使目标同时以垂直极化与水平极化方式成像；有不同的空间分辨率和数据率；可提供7个条带，入射角在15°～45°的雷达数据。
RADARSAT多极化、多入射角ALOSALOS采用了先进的陆地观测技术，能够获取全球高分辨率陆地观测数据。
该卫星载有三种传感器：全色立体测图传感器，新型可见光和近红外辐射计、相控阵型L波段合成孔径雷达（PALSAR）。
PALSAR不受云层，天气和昼夜影响，可全天时全天候对地观测，该卫星具有多入射角，多极化，多工作模式及多种分辨率的特性，最高分辨率可达7m。
（ERS、ENVISAT是多入射角吗？）TerraSAR-XTerraSAR-XTerraSAR-X 是固态有源相控阵的X 波段合成孔径雷达（SAR）卫星，具有多极化、多入射角的特性，具备4 种工作方式和4 种不同分辨率的成像模式。
高分辨率聚束式(High Resolution SpotLight(HS))聚束式(SpotLight Mode(SL))宽扫成像模式(ScanSAR Mode(SC))条带成像模式(Stripmap Mode(SM))COSMO-SkyMedCOSMO-SkyMed星座共包括4颗SAR卫星工作在X波段，具有多极化、多入射角的特性，具备3种工作方式和5种分辨率的成像模式，作为全球第1个分辨率高达1 m的雷达成像卫星星座，COSMO-SkyMed系统将以全天候、全天时对地观测的能力、卫星星座特有的高重访周期和l m高分辨率成像1.2InSAR及发展概况SAR的不足：SAR传感器获取的原始资料主要包含两种信息：一是地面目标区域的二维图像，二是地面目标反射回来的相位SAR成像没有利用回波相位信息。
几种典型的合成孔径雷达系统参数特性介绍(二)
几种典型的合成孔径雷达系统参数特性介绍（二）作者：覃金海来源：《大东方》2018年第04期摘要：本文基于文献检索，主要对目前的典型SAR系统进行了介绍，分别对YSAR和YINSAR、SARMapper系统、TerraSAR-X系统进行了研究与总结，介绍了它们的主要原理以及参数特性。
2.2.1 YSAR和YINSAR自1994年以来在美国国家航空航天局（NASA）的资助之下，美国的杨百翰大学电力与计算机系就开始对低成本和微型的SAR进行探索和研究，并且先后制作出了YSAR和YINSAR的系统样机。
YSAR和YINSAR的系统组成比较简单，通常使用飞行高度在600m左右且为4座或者6座的小型搭载飞机平台，其组成部件多是选取现成的通用部件。
YSAR系统未使用运动补偿装置，成像算法选取的是相位梯度聚焦式（PGA），距离分辨率可以达到1.5 m，方位分辨率可以达到0.5m。
YSAR系统重量约为163 kg（有一半是电池电源的重量），共包括数字子系统、线性调频产生子系统、射频子系统、天线子系统四个部分。
YINSAR系统是基于YSAR系统发展起来的，它采用的是干涉式SAR系统且装备了两根接收天线，且选取了GPS计算运动补偿。
YINSAR的工作载频为9.9 GHz，发射机的峰值功率为10W，整个的YINSAR系统总功率为600 W。
距离分辨率改进到了1.5 m，方位分辨率仍是0.5m。
整个系统由RS-232数据线连接，包括运动测量子系统天线、中频/射频电路与电源子系统、计算机子系统三个部分。
YINSAR系统的软件由GUI（Graphical User Interface）与处理模块、导航模块、数据采集模块、控制模块组成。
目前BYU正在把YINSAR系统和YSAR系统进行结合，形成双频、多模式的SAR系统[1]。
它们的基本参数如表1所示。
2.2.2 SARMapper系统我国首套SARMapper（全称“机载多波段多极化干涉SAR测图系统”）系统于2010年研制成功。
SAR雷达卫星的极化方式
雷达卫星极化方式HH、VV、HV、VH1.什么是SAR极化方式？国科创（北京）信息技术有限公司-极化方式(Polarization）: H水平极化；V 垂直极化，即电磁场的振动方向，卫星向地面发射信号时，所采用的无线电波的振动方向可以有多种方式，目前所使用的有：水平极化（H）：水平极化是指卫星向地面发射信号时，其无线电波的振动方向是水平方向。
垂直极化（V）：垂直极化是指卫星向地面发射信号时，其无线电波的振动方向是垂直方向。
使用H和V线性极化的雷达系统用一对符号表示发射和接收极化，因此可以具有以下通道—HH、VV、HV、VH。
雷达遥感系统常用四种极化方式：（1）HH-用于水平发送和水平接收（2）VV-用于垂直发送和垂直接收（3）HV-用于水平发送和垂直接收（4）VH-用于垂直发送和水平接收这些偏振组合中的前两个被称为相似偏振，因为发射和接收偏振是相同的。
最后两个组合称为交叉极化，因为发送和接收极化彼此正交。
2.什么是SAR中单极化、双极化、全极化？电磁波发射分为水平波（H）和垂直波（V），接收也分为H和V；单极化是指（HH）或者（VV），就是水平发射水平接收或垂直发射垂直接收；如果研究的是气象雷达领域那一般都是（HH）。
双极化是指在一种极化模式的同时，加上了另一种极化模式，如（HH）水平发射水平接收和（HV）水平发射垂直接收。
全极化技术难度最高，要求同时发射H和V，也就是（HH）（HV）（VV）（VH）四种极化方式。
雷达系统可以具有不同级别的极化复杂度：（1）单极化-HH或VV或HV或VH（2）双极化-HH和HV，VV和VH或HH和VV（3）四个极化-HH，VV，HV和VH正交极化（即极化）雷达使用这四个极化，并测量通道之间的相位差以及幅度。
一些双极化雷达还测量通道之间的相位差，因为该相位在极化信息提取中起着重要作用。
雷达卫星影像在极化方面，不同的被观测物体对于入射的不同极化波，后向散射不同的极化波。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
2007年6月15日，德国TerraSAR-X雷达卫星搭乘“第聂伯”火箭从拜科努尔人造卫星发射基地升空。
TrraSAR-X卫星为德国研制的一颗高分辨率雷达卫星，携带一颗高频率的X波段合成孔径雷达传感器，可以聚束式、条带式和推扫式3种模式成像，并拥有多种极化方式。
可全天时、全天候地获取用户要求的任一成像区域的高分辨率影像。
TerraSAR-X/PAZ雷达卫星星座
2015年，与TerraSAR-X完全相同的西班牙卫星PAZ（运行商：Hisdesat）将加入TerraSAR-X的轨道。
两颗卫星将形成一个在覆盖和获取能力方面具有显著提升的卫星星座。
新的星座可支持需要大量数据以及紧急事件等方面的应用，如危机管理、公共安全等。
TerraSAR-X与TanDEM-X采集全球高程模型
与几乎完全相同的双子卫星TanDEM-X一起，为WorldDEMTM 的生产获取基础数据。
WorldDEMTM是一个全球覆盖的、高质量、高精度和大覆盖的数字地形模型。
用于大面积海事监控的宽幅ScanSAR模式
最新的宽幅ScanSAR采集模式，一景数据的最大尺寸可达到405000km²，专为海事监控应用而设计。
宽幅ScanSAR影像可满足针对海上交通，海冰监控及油膜定期探测等应用的大范围持续监测。
用于关键细节的凝视SpotLight
最新的凝视SpotLight采集模式可满足需要物体细节信息的需求。
分辨率最高可以到0.25米，经过进一步优化的辐射质量和非凡的几何传感器精度，可为影像情报分析以及地理信息应用提供精确的细节信息。
详细参数说明：
更多模式访问官网。