一种基于航空激光测高的海冰干舷估测方法及系统

1.本发明属于大地测量与测绘工程领域，特别涉及一种基于航空激光测高的海冰干舷估测方法及系统。


背景技术：

2.海冰是地球系统的重要组成部分，在调解极区大气和海洋相互作用方面发挥着重要作用，获取准确的海冰参数对于研究气候变化相关科学问题至关重要。海冰厚度是海冰众参数中最重要的参数之一，也是最难以观测的参数。海冰干舷高度是指海冰表面相对于局部海平面的高度，是获取海冰厚度和体积的基础数据。大范围的海冰干舷获取手段主要有卫星测高和航空测高等。
3.相较于卫星测高能够获取半球尺度的海冰干舷数据，航空测高方法受飞机航程等因素限制，其观测范围有限，但是其观测数据具有更高的分辨率和更高的精度，其脚点远小于卫星脚点大小。考虑到冬季海冰冻结期内，冰间水道较少且大部分冰间水道宽度在米级到几十米之间，因此卫星测高难以准确识别冰间水道，而航空激光测高数据脚点大小为1m，条带扫描脚点间隔3-4米，能够精确识别冰间水道，得到瞬时海面高程和海冰干舷。因此，航空激光测高为海域尺度海冰变化研究提供了高分辨率、高精度的观测数据。
4.利用航空激光测高估测海冰干舷时，首先提取位于冰间水道内的激光脚点，以此计算局部海面高程，然后根据海冰表面高程和局部海面高程之差，得到海冰干舷高度。一般利用脚点表面反射率进行冰间水道内激光脚点的提取。但是，在冰间水道内，通常存在一定的尼罗冰、灰白冰。此类冰型上的激光脚点表面反射率与冰间水道上的脚点反射率呈现类似的分布特征。因此，利用表面反射率难以有效区分冰间水道和尼罗冰、灰白冰。这将会造成局部海面高值的高估，进而导致海冰干舷的低估。为此，本发明的目的在于提出一种基于航空激光测高的海冰干舷估测方法，能够有效去除尼罗冰、灰白冰对局部海面高计算的影响，以提高海冰干舷估测的精度。


技术实现要素：

5.本发明的目的在于提出一种基于航空激光测高的海冰干舷估测技术方案。
6.本发明所提供的一种基于航空激光测高数据的北极海冰干舷估测方法，包括以下步骤，
7.步骤1，航空激光测高数据处理，包括以下子步骤，
8.步骤1.1，提取航线上激光脚点参数，包括经度、纬度、高程、轨道号和表面反射率；
9.步骤1.2，利用潮汐模型，对激光测高数据进行潮汐校正；
10.步骤1.3，利用平均海平面模型和高程数据，计算激光脚点的剩余高程；
11.步骤2，根据表面反射率，进行冰间水道识别，提取位于冰间水道内的激光脚点；
12.步骤3，局部海面高计算，包括以下子步骤，
13.步骤3.1，根据预设的大步长，对沿轨激光脚点数据进行分段，计算各分段冰间水
道内激光脚点剩余高程平均值和标准差，根据数据分布情况设置标准差倍数作为阈值，剔除剩余高程异常值相应的冰间水道内激光脚点；
14.步骤3.2，根据预设的小步长，对沿轨激光脚点数据进行分段，提取各分段冰间水道内激光脚点剩余高程值最小的若干个点，并以这若干个点的平均剩余高程值作为该分段的局部海面高；步骤3.3，对于不含冰间水道的分段，利用相邻分段局部海面高插值得到该段的局部海面高；
15.步骤4，利用剩余高程和局部海面高，计算并输出海冰干舷预测结果。
16.而且，步骤3中计算局部海面高时，利用数理统计方法对冰间水道内脚点异常值进行了剔除，并通过选取剩余高程值最低的激光脚点作为计算局部海面高的系点，从而避免了冰间水道内薄冰对局部海面高值的高估，提高了海冰干舷估测的可靠性。
17.而且，预设的大步长为50km，预设的小步长为20km。
18.而且，步骤2中，当表面反射率r《0.25，判定相应激光脚点位于冰间水道。
19.而且，步骤3.1中，设置2倍标准差作为阈值，剔除剩余高程在2倍标准差之外的冰间水道内激光脚点。
20.而且，步骤3.2中，根据预设的小步长，对沿轨激光脚点数据进行分段，提取各分段冰间水道内激光脚点剩余高程值最小的10个点，并以这10个点的平均剩余高程值作为该分段的局部海面高。
21.另一方面，本发明还提供一种基于航空激光测高数据的北极海冰干舷估测系统，用于实现如上所述的一种基于航空激光测高数据的北极海冰干舷估测方法。
22.而且，包括以下模块，
23.第一模块，用于航空激光测高数据处理，包括以下操作，
24.提取航线上激光脚点参数，包括经度、纬度、高程、轨道号和表面反射率；
25.利用潮汐模型，对激光测高数据进行潮汐校正；
26.利用平均海平面模型和高程数据，计算激光脚点的剩余高程；
27.第二模块，用于根据表面反射率，进行冰间水道识别，提取位于冰间水道内的激光脚点；
28.第三模块，用于局部海面高计算，包括以下操作，
29.根据预设的大步长，对沿轨激光脚点数据进行分段，计算各分段冰间水道内激光脚点剩余高程平均值和标准差，根据数据分布情况设置标准差倍数作为阈值，剔除剩余高程异常值相应的冰间水道内激光脚点；
30.根据预设的小步长，对沿轨激光脚点数据进行分段，提取各分段冰间水道内激光脚点剩余高程值最小的若干个点，并以这若干个点的平均剩余高程值作为该分段的局部海面高；
31.对于不含冰间水道的分段，利用相邻分段局部海面高插值得到该段的局部海面高；
32.第四模块，用于利用剩余高程和局部海面高，计算并输出海冰干舷预测结果。
33.或者，包括处理器和存储器，存储器用于存储程序指令，处理器用于调用存储器中的存储指令执行如上所述的一种基于航空激光测高数据的北极海冰干舷估测方法。
34.或者，包括可读存储介质，所述可读存储介质上存储有计算机程序，所述计算机程
序执行时，实现如上所述的一种基于航空激光测高数据的北极海冰干舷估测方法。
35.本发明提出的一种基于航空激光测高的海冰干舷估测方案，能够有效去除尼罗冰、灰白冰对局部海面高计算的影响，以提高海冰干舷估测结果的可靠性。
附图说明
36.图1为本发明实施例的流程图。
具体实施方式
37.下面结合附图和实施例对本发明做进一步详细描述。
38.参见图1，本发明实施例提供的一种基于航空激光测高的海冰干舷估测方法，包括以下步骤：
39.步骤1，航空激光测高数据处理，包括以下子步骤：
40.步骤1.1，提取包括经度、纬度、高程、轨道号、表面反射率等参数：实施例基于icebridge航空激光测高数据，以icebridge航线20140428_1为例，提取航线上激光脚点的经度(lon)、纬度(lat)、高程(h)、表面反射率(r)、轨道号(orb)等；
41.步骤1.2，利用潮汐模型，对激光测高数据进行潮汐校正：实施例利用tpx09潮汐模型，对激光脚点进行潮汐校正，得到潮汐改正后的脚点高程h
corr
：
42.h
corr
＝h-h
tide
43.其中，h为激光脚点高程，h
tide
为潮汐改正，
44.步骤1.3，利用平均海平面模型和潮汐改正后的脚点高程数据，计算激光脚点的剩余高程hr：
45.hr＝h
corr-h
mss
46.其中，h
mss
为平均海平面高改正，h
corr
为潮汐改正后的脚点高程，
47.步骤2，根据表面反射率，进行冰间水道识别，提取冰间水道内的激光脚点：在不同反射面(海冰、冰间水道)上激光脚点表面反射率呈现不同的特征，因此可根据表面反射率r，进行冰间水道识别，提取冰间水道内激光脚点。在实施例中，优选设置当r《0.25，判定该脚点位于冰间水道。
48.步骤3，局部海面高计算，实施例采用的优先实施方式，包括以下子步骤，
49.步骤3.1，以50km为步长，对沿轨激光脚点数据进行分段，计算各分段冰间水道内激光脚点剩余高程平均值h
50km
和标准差σ，剔除剩余高程在2倍标准差之外的冰间水道内激光脚点hr＜h
50km-2σ或hr＞h
50km
+2σ；
50.步骤3.2，以20km为步长，对沿轨激光脚点数据进行分段，提取各分段冰间水道内激光脚点剩余高程值最小的10个点，并以这10个点的平均剩余高程值作为该分段的局部海面高h
ssh
；
51.步骤3.3，对于不含冰间水道的分段(实施例为20km分段)，利用相邻分段局部海面高插值得到该段的局部海面高；
52.步骤3中计算局部海面高时，利用数理统计方法对冰间水道内脚点异常值进行了剔除，并通过选取剩余高程值最低的激光脚点作为计算局部海面高的系点，从而避免了冰间水道内薄冰对局部海面高值的高估，提高了海冰干舷估测的可靠性。
53.步骤4，利用沿轨剩余高程hr和局部海面高h
ssh
，计算沿轨海冰干舷hf：
54.hf＝h
r-h
ssh
55.具体实施时，用于可以根据情况设置阈值、步长等取值，或者采用经验值。
56.具体实施时，本发明技术方案提出的方法可由本领域技术人员采用计算机软件技术实现自动运行流程，实现方法的系统装置例如存储本发明技术方案相应计算机程序的计算机可读存储介质以及包括运行相应计算机程序的计算机设备，也应当在本发明的保护范围内。
57.在一些可能的实施例中，提供一种基于航空激光测高数据的北极海冰干舷估测系统，包括以下模块，
58.第一模块，用于航空激光测高数据处理，包括以下操作，
59.提取航线上激光脚点参数，包括经度、纬度、高程、轨道号和表面反射率；
60.利用潮汐模型，对激光测高数据进行潮汐校正；
61.利用平均海平面模型和高程数据，计算激光脚点的剩余高程；
62.第二模块，用于根据表面反射率，进行冰间水道识别，提取位于冰间水道内的激光脚点；
63.第三模块，用于局部海面高计算，包括以下操作，
64.根据预设的大步长，对沿轨激光脚点数据进行分段，计算各分段冰间水道内激光脚点剩余高程平均值和标准差，根据数据分布情况设置标准差倍数作为阈值，剔除剩余高程异常值相应的冰间水道内激光脚点；
65.根据预设的小步长，对沿轨激光脚点数据进行分段，提取各分段冰间水道内激光脚点剩余高程值最小的若干个点，并以这若干个点的平均剩余高程值作为该分段的局部海面高；
66.对于不含冰间水道的分段，利用相邻分段局部海面高插值得到该段的局部海面高；
67.第四模块，用于利用剩余高程和局部海面高，计算并输出海冰干舷预测结果。
68.在一些可能的实施例中，提供一种基于航空激光测高数据的北极海冰干舷估测系统，包括处理器和存储器，存储器用于存储程序指令，处理器用于调用存储器中的存储指令执行如上所述的一种基于航空激光测高数据的北极海冰干舷估测方法。
69.在一些可能的实施例中，提供一种基于航空激光测高数据的北极海冰干舷估测系统，包括可读存储介质，所述可读存储介质上存储有计算机程序，所述计算机程序执行时，实现如上所述的一种基于航空激光测高数据的北极海冰干舷估测方法。
70.尽管上面结合附图对本发明的较佳实施例进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，并不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可以作出很多形式，这些均属于本发明的保护范围之内。