拉布拉多海一年平整冰厚度SAR反演算法
刘眉洁;戴永寿;张杰;张晰;孟俊敏
【摘 要】在利用实测一年平整冰厚数据和同步全极化RADARSAT-2 SAR影像、全面分析不同SAR极化参数对海冰厚度敏感度的基础上,发现Alpha角最高;从理论上证实利用两者相关性反演冰厚的可行性,将所得经验方程用于冰厚反演,其结果与实测数据吻合较好;与现有常用方法的反演结果进行对比表明,新方法误差较小;Alpha角可以用于一年平整冰厚度反演,证明全极化SAR海冰信息反演冰厚确有优势.
【期刊名称】《中国石油大学学报（自然科学版）》
【年(卷),期】2014(038)003打电话 英文
【总页数】7页(P186-192)
【关键词】海冰厚度;Alpha角;海冰厚度反演算法;极化SAR;极化散射特性;拉布拉多海
【作 者】刘眉洁;戴永寿;张杰;张晰;孟俊敏
【作者单位】中国石油大学信息与控制工程学院,山东青岛266580;国家海洋局第一海洋研究所,山东青岛266061;青岛大学物理科学学院,山东青岛266071;中国石油大学信息与控制工程学院,山东青岛266580;国家海洋局第一海洋研究所,山东青岛266061;国家海洋局第一海洋研究所,山东青岛266061;国家海洋局第一海洋研究所,山东青岛266061
【正文语种】中 文
【中图分类】P731.32
海冰厚度影响着气候变化和海-气能量与物质交换。一年冰对海-气界面的物理过程影响尤为显著，其厚度提取具有重要研究意义［1-4］。现场勘测、高度计探测和航标测量等现有冰厚获取手段难以同时满足长期、大范围和高分辨率的需求。全极化(SAR)能够提供时空连续的高分辨率海冰散射信息。海冰SAR散射特性主要依赖于海冰表面粗糙度和复介电常数［5-6］，两者又与海冰生长过程密切相关，可将散射特性与冰厚联系起来用于冰厚的提取［7-15］。一年平整冰的散射特性可能受介电常数影响较大，所以需更多关注与介电常数相关的极化参数，如Alpha角(α)。笔者利用拉布拉多海全极化RADARSAT-2 SAR影像和同步现场冰厚数据，分析海冰SAR极化特性对一年平整冰厚度的响应特性，选出敏感度最
高的特性参数，并从理论上证明其有效性，进而得到相关性的经验方程，从而进行冰厚反演和误差分析，并与其他方法的结果进行对比分析。
1 研究区域和数据源
SAR数据为拉布拉多海RADARSAT-2全极化单视复影像，共四景，每景影像分辨率8 m，覆盖范围25 km×25 km(表1)，并都进行了辐射定标和几何校正，如图1所示。
拉布拉多海包括一年冰和多年冰，四景SAR影像以一年平整冰和一年变形冰为主。图2为影像在拉布拉多海覆盖区域示意图。
表1 RADARSAT-2全极化SAR参数Table 1 Properties of RADARSAT-2 SAR获取时间 入射角 /(°)2011年3月19日21:51:02 42 2011年3月19日21:51:05 42 2011年3月20日09:56:00 49 2011年3月19日10:25:08 29
石家庄北大青鸟学校图1 伪彩色合成图像Fig.1 Fal colour composite images
图2 SAR影像在拉布拉多海覆盖区域示意图Fig.2 Map of test sites
sfzSAR影像的同步现场海冰厚度数据来自加拿大贝德福德海洋研究所，是利用机载电磁感应设备(HEM)获取的。HEM通常由激光测距仪和电磁感应设备组成。前者测量海冰表面到接收器的距离，后者获取冰水分界面到接收器的距离，两个距离的差值即为海冰(和雪)厚度［17］。该设备测量平整冰厚度时精度较高，误差±0.1 m;但测量变形冰时误差较大，可达实际冰厚的 50%［1，15］。
因为研究对象是一年平整冰，所以通过影像解译，尽量剔除变形冰的HEM数据，只用平整冰数据，这样也可以保证所用冰厚数据具有较高精度，如图3所示(2011年3月20日09:56:00的SAR影像)。
2 海冰厚度和海冰后向散射特征的相关性分析
2.1 样本的选取
先根据四景影像的经纬度范围，提取该范围内0～2 m的HEM一年平整冰厚度数据;再以HEM数据中每一海冰厚度数据的经纬度值为中心点，利用3×3的窗口从RADARSAT-2影像中获取该点的平均极化特征参数值，并将其与对应的冰厚关联起来。由于同一海冰厚度对
应的每一SAR极化特征参数值可能会在一定范围内变化，直接用上面得到的特征参数值和冰厚数据进行分析，可能不容易分析两者的相关性。因此对上述数据做平均处理，在0～2 m范围内，每0.1 m为一段，分为20段，并对同一厚度段内的极化特征参数值进行平均，作为这一厚度段的SAR极化特征参数值，用于分析一年平整冰厚度的SAR极化响应特性。
图3 剔除变形冰前后的HEM和SAR影像匹配图对比Fig.3 Comparison of deformed ice unadjusted and eliminated
2.2 相关性分析
为了更好地体现不同海冰SAR极化特性对海冰厚度的响应特性，先将海冰SAR极化特性分为6组，再分别分析它们与海冰厚度的相关性。利用本文数据得到的一年平整冰厚度和极化特性的相关系数如表2所示，其中下标“H、V、R、L”分别表示水平、垂直、右圆和左圆极化。
表2 所用的极化特性参数Table 2 Sea-ice polarization parameters in six groups极化特征参数 相关系数R σHH后向散射系数0.408 σHV 0.569σVH 0.667 σVV 0.033 PhaHH极化相
位-0.157 PhaHV -0.095 PhaVH -0.312 PhaVV -0.413 σHV/σHH极化比0.556 σVH/σVV 0.668 σVV/σHH -0.767 PhaHH—PhaHV极化相位差0.060 PhaHH—PhaVV 0.058 PhaVH—PhaVV -0.155 0.398 ρVV-HH -0.640 ρVV-VH -0.515圆极化ρRR-LL ρVH-HH极化相关系数-0.474 0.718反熵A -0.508 Alpha角α散射熵H H/A/α 分解0.804
从表2可见，一年平整冰的Alpha角与海冰厚度相关性最高，同极化比(σVV/σHH)次之(图4)。Alpha角是从全极化SAR数据中提取出来的，描述了探测目标内部的自由度，说明了目标平均散射的“类型”。以往文献中极少研究这两者的相关性。windriver
图4 海冰厚度和极化特性的相关性Fig.4 Correlations between a-ice thickness and polarization parameters
3 海冰厚度和Alpha角相关性的理论分析
根据本文中数据得到的Alpha角和一年平整冰厚度的相关性较高，由此可以建立两者的经验方程，用以反演海冰厚度。
在海冰生长过程中，海冰表面复介电常数随着海冰厚度的增加而变化，即可以表示为海冰
厚度的函数:ε=ε(h)。其中ε表示海冰表面复介电常数;h表征海冰厚度，m。
海冰的Alpha角主要依赖于地物的表面复介电常数［18］，即Alpha角是海冰表面复介电常数的函数:α=α(ε)。因而，海冰厚度和Alpha角可以借助海冰复介电常数这一中间变量从理论上建立一定的函数关系:α=α(h)。
3.1 海冰厚度和复介电常数的函数关系
在海冰生长增厚的过程中，海冰的盐度不断降低，介电常数也随之减小。为便于分析计算，在不影响研究结果的情况下，用相对介电常数替代介电常数。Vant等［19］发现海冰介电常数ε与海冰厚度满足函数关系:
式中，εi和εb分别为纯冰和卤水的介电常数;Vb为卤水体积分量。
卤水体积分量的表达式可以写为［20］
neutral是什么意思式中，t为海冰温度，0～2 m海冰的温度范围满足此式的要求［21］;Si为海冰盐度。
海冰盐度Si与海冰厚度h的函数关系［22］可写成
devious maids
结合式(2)、(3)，可以推得海冰卤水体积分量随海冰厚度的增加而减小，如图5所示。
综合式(1)～(3)，海冰的介电常数可以表示为海冰厚度的函数，即
其中，ni(i=0，1，2)可通过式(1)～(3)计算得到。由于海冰复介电常数随卤水体积分量的增加而增加，因而会随着海冰厚度的增加而减小，如图6所示。
3.2 海冰复介电常数和α的关系
从海冰后向散射系数中提取的α主要依赖于海冰表面复介电常数。一年平整冰的表面粗糙度较小，因而可根据Bragg散射理论［18］，将α表示为
图5 海冰厚度和卤水体积分量关系曲线Fig.5 Relationship between a-ice thickness and brine volume
图6 海冰厚度和介电常数的关系曲线Fig.6 Relationship between a-ice thickness and dielectric constant
其中
式中，RS和RP分别为平行和垂直于入射面的Bragg散射系数，皆为海冰介电常数ε和入射角θ(rad)的函数［18］。
设 R=RP/RS，则
新年快乐的英语设 y= ε -1，x=cos θ，则
从式(7)可以看出，在入射角θ一定的情况下，随着卤水体积分量的增加而增加，如图7所示，那么也是海冰复介电常数的增函数。
对式(7)进一步化简，并考虑入射角一定的情况下，有
其中为待定系数。由此可得
英文歌曲网
图7 卤水体积分量和的关系曲线Fig.7 Relationship between brine volume and R
从式(9)可以导出cos α在不同入射角下都随着的增加而减小，图8仅以60°入射角为例。因而α是的增函数，也会随着海冰复介电常数的增加而增加。
图8 和cosα的关系曲线Fig.8 Relationship between and cosα
who says3.3 海冰厚度和α的关系
式(9)给出了海冰复介电常数和α的关系。为了便于计算，化简式(9)可得:
其中，bi(i=0，…，6)为待定系数。
将海冰厚度和复介电常数的函数方程式(4)带入式(10)，可得α和海冰厚度的关系式为
spiced式中，h表示海冰厚度;ci(i=0，…，6)为待定系数。
由于式(11)中的高阶项起到的作用可以忽略，因而只保留一阶项:
其中，di(i=0，1)为待定系数。
由于α随着海冰复介电常数的增加而增加，而海冰复介电常数是海冰厚度的减函数，因而α随着海冰厚度的增加而减小。由此可见，α是冰厚的函数，两者具有明显的相关性，这与图4(b)所显示出来的两者强相关性相符。但是，图4(b)中表现出来的是正相关，而理论结果是负相关。这个问题还须做进一步研究。
4 海冰厚度和α反演算法
在海冰厚度和海冰后向散射特征的相关性分析中已将海冰SAR影像中的极化数据和HEM测量的同步一年平整冰厚度数据进行了时空匹配，因而在此将匹配好的数据随机分为两组，一组用于确定待定系数di，另一组用于反演海冰厚度并计算反演误差。
