基于被动微波多系点反演海冰厚度的方法及系统

1.本发明属于摄影测量和遥感领域，特别涉及基于被动微波多系点反演海冰厚度的方法及系统。


背景技术：

2.在过去的100年中，世界各地气象站记录的地表长时序温度数据表明，全球气温上升了0.6℃,而北极地区的局部变暖速度在1961年至1990年期间甚至超过了每十年0.5℃。两极地区覆盖有全球最大范围的冰雪地表，这使得在全球变暖的大背景下，两极地区对气候变化表现最为敏感。尤其是极地海冰作为气候系统的重要因子，对全球气候变化有着指示器作用。由于海洋-大气系统的相互作用，海冰的变化特征相比于陆地冰雪更加明显，如海冰范围变小、厚度变薄、冰型变化显著以及漂移和融化再冻结等。
3.北极海冰作为海洋和大气之间的天然屏障，调节海洋和大气的热量和辐射平衡，使区域保持低能量状态，成为地球的热汇和冷源，以一种复杂的方式影响全球气候系统，进而影响其循环和温度。通过雪、冰的高反照率等因素相关的反馈系统对外部驱动因素的快速响应，北极海冰变化提供了全球气候变化的早期指标。因此，获取准确的海冰参数(厚度、密集度、范围、反照率等)至关重要。
4.海冰厚度作为海冰的重要参数之一，对气候变化和实际应用均有十分重要的意义。海冰厚度对大气-海冰-海洋耦合作用过程的影响尤为显著，决定了海-气能量与水汽的交换过程和速率，主导着海冰的物理力学特征，影响海冰的运动、形变以及海冰的冻结与消融过程。近年来北极海冰变化十分活跃，海冰覆盖范围在不断退缩，冰层厚度在不断变薄，而厚度小于半米的薄冰更是在海洋-大气热交换中占主导地位。因此，能够准确反演薄冰厚度十分重要。若需获取大范围的海冰厚度，对于厚冰，主要利用卫星测高数据反演；对于薄冰，主要利用遥感数据反演。
5.被动微波遥感技术因为不受天气影响逐渐成为众多学者反演大尺度海冰厚度的方法，依据被动微波亮温值反演海冰厚度，原理在于目标物体的微波辐射与自身的物理性质相关，由于海冰内部结构复杂，其微波辐射能量与海水不同，因此微波辐射计探测到的亮温值可以轻易的分辨出海冰与海水，并且亮温值随海冰厚度的增长而增长，并逐渐趋于稳定。传统的方法是利用单系点对反演研究区域甚至是全北极海冰厚度初值，然而，不同海域的海水、海冰由于温度、盐度的不同会影响亮温值的不同，进而导致反演海冰厚度的不确定性较大。因此本发明提出构建多系点对，代表不同海域，共同反演海冰厚度，降低不确定性。


技术实现要素：

6.本发明的目的在于针对现有薄冰厚度估算方法中存在的不确定性，提供基于被动微波的多系点反演海冰厚度的方法及系统。
7.本发明所提供一种基于被动微波多系点反演海冰厚度的方法，包括以下步骤，
8.步骤1，卫星数据处理，提取被动微波数据包括经度、纬度、多角度多极化下的亮温
观测值；提取海冰密集度数据包括经度、纬度和海冰密集度；
9.步骤2，海冰亮温值整合，包括以下子步骤，
10.步骤2.1，根据阈值法，进行射频干扰识别；
11.步骤2.2，利用剔除射频干扰后的多角度多极化亮温观测值，计算得出海冰亮温值强度；
12.步骤2.3，根据预设的分辨率，计算每个格网单元范围内的亮温观测值做平均；
13.步骤3，海冰厚度反演，包括以下子步骤，
14.步骤3.1，以海冰密集度产品格网为基础，进行数据处理；
15.步骤3.2，将步骤2中提取的海冰亮温值，根据最近邻法插值到对应格网单元内；
16.步骤3.3，在海冰冻结期内，构建格网单元内的海冰亮温值与海冰密集度的长时间序列，得出系点对，即稳定的海水亮温值和海冰亮温值；
17.步骤3.4，得出系点对，引入半经验趋近方程，与亮温观测值方程联立，带入格网单元内的亮温观测值时间序列，根据最小二乘优化，得出衰减因子；
18.步骤4，在海冰-海水交界区域挑选重复步骤3，获取不同海域的不同系点对，最后，应用全部的系点对反演海冰厚度。
19.而且，步骤2.2中，计算得出海冰亮温值强度t
b,i
如下，
[0020][0021]
其中，t
b,h
、t
b,v
分别是地球参考框架中的水平和垂直极化的亮度温度。
[0022]
而且，步骤3.4中，得出系点对，引入半经验趋近方程：
[0023]
t
obs
＝t
1-(t
1-t0)exp(-γd)
[0024][0025]
式中，t0、t1为海水、海冰系点，t
obs
为亮温观测值，γ为衰减因子，d为海冰厚度；
[0026]
与亮温观测值方程联立：
[0027]
t
abs
＝ce
ice
t
ice
(1c)e
water
t
water
[0028]
式中，c为海冰密集度，亮温观测值t
obs
取决于海水亮温值t
water
、海冰亮温值t
ice
及海水发射率e
water
、海冰发射率e
ice
。
[0029]
另一方面，本发明还提供一种基于多系点的被动微波海冰厚度反演系统，用于实现如上所述的一种基于多系点的被动微波海冰厚度反演方法。
[0030]
而且，包括处理器和存储器，存储器用于存储程序指令，处理器用于调用存储器中的存储指令执行如上所述的一种基于多系点的被动微波海冰厚度反演方法。
[0031]
或者，包括可读存储介质，所述可读存储介质上存储有计算机程序，所述计算机程序执行时，实现如上所述的一种基于多系点的被动微波海冰厚度反演方法。
[0032]
本发明提出的基于被动微波多系点对反演海冰厚度的方法，能够避免海冰厚度反演中由单系点反演带来的不确定性，提高海冰厚度反演的可靠性。
附图说明
[0033]
图1为本发明实施例的流程图；
[0034]
图2为本发明实施例的不同海冰厚度下的相干和非相干发射率示意图。
具体实施方式
[0035]
下面结合附图和实施例对本发明做进一步详细描述。
[0036]
本发明实施例提供的一种基于被动微波多系点反演海冰厚度的方法，实现基于多系点的被动微波海冰厚度反演，如图1所示：
[0037]
步骤1，卫星被动微波数据处理，包括提取海冰区域内各格网单元的相关数据；
[0038]
实施例基于土壤湿度与海洋盐度(smos)卫星被动微波亮温数据，以北极区域海冰区为例，提取了40
°
～225
°
e，66
°
～85
°
n范围内各格网单元的经度(lon)、纬度(lat)、亮温(tb)、极化方式(p)、入射角度(inc)、法拉第旋转角(fa)、快照编号(n)等。
[0039]
步骤2，海冰亮温值整合
[0040]
实施例进一步提出的实现过程包括子步骤如下：
[0041]
步骤2.1，射频干扰识别
[0042]
实施例提出根据阈值法，进行射频干扰识别。利用被动微波亮温数据反演海冰厚度时，射频干扰识别是关键。在接收微波信号过程中，由于雷达和无线电传输的影响，部分smos观测值会受射频干扰(rfi)的影响。同时入射角度、极化模式和卫星的上升、下降轨道也会对rfi造成影响。为了简化，应用阈值法在水平和垂直极化亮温观测值中对rfi进行判断。优选地，如果其中任何极化亮温观测值在快照内超出300k或者小于0k，则此快照被视为已被rfi污染。
[0043]
步骤2.2，利用剔除射频干扰后的多角度多极化亮温观测值，计算得出海冰亮温值强度tb,i：
[0044][0045]
其中，t
b,h
、t
b,v
分别是地球参考框架中的水平和垂直极化的亮度温度。
[0046]
步骤2.3，根据预设的海冰亮温值格网单位分辨率，计算每个格网内的亮温均值：实施例中优选以15km
×
15km格网单元为分辨率，计算每个格网内的亮温均值。
[0047]
步骤3，海冰厚度反演
[0048]
依据被动微波亮温值反演海冰厚度，原理在于目标物体的微波辐射与自身的物理性质相关，由于海冰内部结构复杂，其微波辐射能量与海水不同，因此微波辐射计探测到的亮温值可以分辨出海冰与海水，并且亮温值随海冰厚度的增长而增长，并逐渐趋于稳定。图2为利用典型北极环境模拟不同的海冰厚度下相干/非相干发射率方程，可以看出相干发射率具有较强的周期性，单一发射率对应不同的海冰厚度，非相干发射率随海冰厚度增长并趋于稳定，适用于反演薄冰厚度。本发明实施例进一步提出的基于多系点被动微波亮温反演海冰厚度方法如下：
[0049]
步骤3.1，根据海冰密集度的分辨率进行数据处理：实施例中优选以海冰密集度12.5km
×
12.5km格网单元为基础，进行数据处理；
[0050]
现有被动微波反演海冰厚度分辨率为12.5km。在选取系点对的过程中，需要构建同一格网单元内的海冰亮温和海冰密集度时间序列；根据时间序列来确定海冰冻结初期和末期，依次确定海水、海冰系点值；为了同步海冰亮温值和海冰密集度，所以将12.5km格网
单元为基础，进行数据处理。
[0051]
步骤3.2，将步骤2中提取的海冰亮温值，根据最近邻法插值到对应格网单元内。
[0052]
步骤3.3，在海冰冻结期内，构建格网单元内的海冰亮温值与海冰密集度的长时间序列，得出系点对，即稳定的海水亮温值和海冰亮温值。
[0053]
步骤3.4，得出系点对，引入半经验趋近方程：
[0054]
t
obs
＝t
1-(t
1-t0)exp(-γd)
[0055][0056]
式中t0、t1为海水、海冰系点，t
obs
为亮温观测值，γ为衰减因子，d为海冰厚度。
[0057]
与亮温观测值方程联立：
[0058]
t
abs
＝ce
ice
t
ice
(1-c)e
water
t
water
[0059]
式中，c为海冰密集度，亮温观测值t
obs
取决于海水亮温值t
water
、海冰亮温值t
ice
及海水发射率e
water
、海冰发射率e
ice
。
[0060]
带入此格网单元内的亮温观测值时间序列，根据最小二乘优化，得出衰减因子。
[0061]
步骤4，空白数据填补
[0062]
实施例沿海冰-海水交界区域，在不同海域重复步骤3.3，获取不同海域的系点对。最后，应用全部的系点对反演海冰厚度。最后得到的2010年北极部分区域的海冰厚度。
[0063]
步骤1中的北极区域海冰区是实施例的实验区域，实验区域包含不同的海域，选取单系点对反演海冰厚度，不足以反应实验区域内的全部海冰、海水性质，因此重复步骤3.3,构建不同海域各格网单元的海冰亮温和海冰密集度的时间序列，确定各海域内的代表系点对。根据全部的系点对，共同反演实验区域内的海冰厚度。
[0064]
具体实施时，本发明技术方案提出的方法可由本领域技术人员采用计算机软件技术实现自动运行流程，实现方法的系统装置例如存储本发明技术方案相应计算机程序的计算机可读存储介质以及包括运行相应计算机程序的计算机设备，也应当在本发明的保护范围内。
[0065]
在一些可能的实施例中，提供一种基于多系点的被动微波海冰厚度反演系统，包括处理器和存储器，存储器用于存储程序指令，处理器用于调用存储器中的存储指令执行如上所述的一种基于多系点的被动微波海冰厚度反演方法。
[0066]
在一些可能的实施例中，提供一种基于多系点的被动微波海冰厚度反演系统，包括可读存储介质，所述可读存储介质上存储有计算机程序，所述计算机程序执行时，实现如上所述的一种基于多系点的被动微波海冰厚度反演方法。
[0067]
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代，但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。