基于GPS双频信号的雪水当量估计方法及装置

基于gps双频信号的雪水当量估计方法及装置
技术领域
1.本发明属于气象监测技术领域，具体涉及一种基于gps双频信号的雪水当量估计方法及装置。


背景技术：

2.积雪是许多国家与地区重要的淡水来源，也是气候变化与水文研究的重点内容。雪水当量作为重要的积雪参数，能够直接表征区域内水循环变化及能量平衡状况。因此，对于雪水当量的稳定、精准、可持续性估计对进行气候、生态变化研究、水资源科学管理等研究意义重大。
3.目前，雪水当量主要通过人工或站点观测的方式获取，相比雪水当量，积雪深度在野外更容易测量，因此，利用雪深估计雪水当量成为了获取雪水当量估计值的一种重要方式。相关技术中，gps-ir技术利用易于获取的gps数据实现积雪深度变化的监测，但是，在利用信号进行雪深反演的过程中，其采用抗差估计等方法进行粗差的剔除，这样无法削弱信号间的系统差异；尤其是对常用于雪深监测的gps l1与l2c信号来说，信号强度及信噪比噪声的影响会导致两种信号在反演精度、数据数量上有一定差异，进而导致雪水当量的估计结果存在较大偏差。


技术实现要素：

4.为了解决现有技术中存在的上述问题，本发明提供了一种基于gps双频信号的雪水当量估计方法及装置。本发明要解决的技术问题通过以下技术方案实现：
5.第一方面，本发明提供一种基于gps双频信号的雪水当量估计方法，应用于gps接收机；
6.所述雪水当量估计方法包括：
7.接收卫星发射的l1信号，并获得所述l1信号的第一信噪比观测值；
8.根据所述第一信噪比观测值和所述l1信号的波长，确定第一反射高度；所述第一反射高度包括水文年每个年积日中l1信号经积雪表面反射后的反射高度；
9.接收卫星发射的l2c信号，并获得所述l2c信号的第二信噪比观测值；
10.根据所述第二信噪比观测值和所述l2c信号的波长，确定第二反射高度；所述第二反射高度包括至少部分水文年年积日中l2c信号经积雪表面反射后的反射高度；
11.对将l2c信号缺测天数对应的第一反射高度进行转换，得到l2c信号缺测天数对应的第三反射高度；
12.根据第二反射高度和第三反射高度，确定水文年中每个年积日的积雪深度；
13.根据所述积雪深度，估计水文年中每个年积日的雪水当量。
14.在本发明的一个实施例中，所述根据所述第一信噪比观测值和所述l1信号的波长，确定第一反射高度的步骤，包括：
15.根据第一信噪比观测值，获得第一信噪比观测值的残差序列；
16.对第一信噪比观测值的残差序列进行频谱分析，确定所述第一信噪比观测值的振荡频率；
17.根据所述l1信号的波长以及所述振荡频率，按照如下公式确定第一反射高度：
[0018][0019]
式中，f1表示所述第一信噪比观测值的振荡频率，λ1表示l1信号的波长，h
r1
表示第一反射高度；其中，所述第一反射高度包括水文年每个年积日中l1信号经积雪表面反射后的反射高度。
[0020]
在本发明的一个实施例中，所述对将l2c信号缺测天数对应的第一反射高度进行转换，得到l2c信号缺测天数对应的第三反射高度的步骤，包括：
[0021]
获取l2c信号缺测天数对应的第一反射高度；
[0022]
确定第一线性拟合系数及第二线性拟合系数，利用所述第一线性拟合系数和所述第二线性拟合系数，将l2c信号缺测天数对应的第一反射高度转换得到l2c信号缺测天数对应的第三反射高度。
[0023]
在本发明的一个实施例中，按照如下公式对l2c信号缺测天数对应的第一反射高度进行转换：
[0024]hr3
＝a
·h′
r1
b
[0025]
其中，a、b分别表示第一线性拟合系数和第二线性拟合系数，h
′
r1
表示l2c信号缺测天数对应的第一反射高度，h
r3
表示l2c信号缺测天数对应的第三反射高度。
[0026]
在本发明的一个实施例中，所述根据第二反射高度和第三反射高度，确定水文年中每个年积日的积雪深度的步骤，包括：
[0027]
计算水文年中每个年积日对应的第二反射高度或第三反射高度与gps接收机自身高度的差值，得到水文年中每个年积日的积雪深度。
[0028]
在本发明的一个实施例中，所述gps接收机自身高度为：无雪期内l2c信号每日反射高度的平均值。
[0029]
在本发明的一个实施例中，根据所述积雪深度，按照如下公式估计水文年中每个年积日的雪水当量：
[0030][0031]
式中，ρ0为雪密度常数，ρ
max
为所述水文年中雪密度的最大值，k1、k2为预设的经验系数，doyi为所述水文年中的第i个年积日，hi为所述水文年中的第i个年积日的积雪深度，为所述水文年中的第i个年积日的雪水当量。
[0032]
第二方面，本发明还提供一种基于gps双频信号的雪水当量估计装置，应用于gps接收机；
[0033]
所述雪水当量估计方法包括：
[0034]
第一接收单元，用于接收卫星发射的l1信号，并获得所述l1信号的第一信噪比观测值；
[0035]
第一确定单元，用于根据所述第一信噪比观测值和所述l1信号的波长，确定第一反射高度；所述第一反射高度包括水文年每个年积日中l1信号经积雪表面反射后的反射高
度；
[0036]
第二接收单元，用于接收卫星发射的l2c信号，并获得所述l2c信号的第二信噪比观测值；
[0037]
第二确定单元，用于根据所述第二信噪比观测值和所述l2c信号的波长，确定第二反射高度；所述第二反射高度包括至少部分水文年年积日中l2c信号经积雪表面反射后的反射高度；
[0038]
转换单元，用于对将l2c信号缺测天数对应的第一反射高度进行转换，得到l2c信号缺测天数对应的第三反射高度；
[0039]
深度确定单元，用于根据第二反射高度和第三反射高度，确定水文年中每个年积日的积雪深度；
[0040]
估计单元，用于根据所述积雪深度，估计水文年中每个年积日的雪水当量。
[0041]
与现有技术相比，本发明的有益效果在于：
[0042]
本发明提供了一种基于gps双频信号的雪水当量估计方法及装置，通过对l1信号的第一反射高度和l2c信号的第二反射高度建立线性模型，将l2c反演精度与l1反演数据量进行有效结合，克服了现有技术中l1信号雪水当量监测精度低、而l2c信号数据量少的问题，不仅能够实现gps-ir积雪深度反演，也提高了雪水当量估计结果的准确性。
[0043]
以下将结合附图及实施例对本发明做进一步详细说明。
附图说明
[0044]
图1是本发明实施例提供的基于gps双频信号的雪水当量估计方法的一种流程图；
[0045]
图2是本发明实施例提供的基于gps双频信号的雪水当量估计方法的一种示意图；
[0046]
图3是本发明实施例提供的雪水当量估计结果的一种示意图；
[0047]
图4是本发明实施例提供的基于gps双频信号的雪水当量估计装置的一种结构示意图。
具体实施方式
[0048]
下面结合具体实施例对本发明做进一步详细的描述，但本发明的实施方式不限于此。
[0049]
图1是本发明实施例提供的基于gps双频信号的雪水当量估计方法的一种流程图，图2是本发明实施例提供的基于gps双频信号的雪水当量估计方法的一种示意图。请参见图1-2，一种基于gps双频信号的雪水当量估计方法，其特征在于，应用于gps接收机；
[0050]
所述雪水当量估计方法包括：
[0051]
s1、接收卫星发射的l1信号，并获得l1信号的第一信噪比观测值；
[0052]
s2、根据第一信噪比观测值和l1信号的波长，确定第一反射高度；第一反射高度包括水文年每个年积日中l1信号经积雪表面反射后的反射高度；
[0053]
s3、接收卫星发射的l2c信号，并获得l2c信号的第二信噪比观测值；
[0054]
s4、根据第二信噪比观测值和l2c信号的波长，确定第二反射高度；第二反射高度包括至少部分水文年年积日中l2c信号经积雪表面反射后的反射高度；
[0055]
s5、对将l2c信号缺测天数对应的第一反射高度进行转换，得到l2c信号缺测天数
对应的第三反射高度；
[0056]
s6、根据第二反射高度和第三反射高度，确定水文年中每个年积日的积雪深度；
[0057]
s7、根据积雪深度，估计水文年中每个年积日的雪水当量。
[0058]
具体而言，在本发明提供的基于gps双频信号的雪水当量估计方法中，gps接收机首先接收gps卫星或gnss卫星发射的l1信号，部分l1信号直接被gps接收机接收、另一部分l1信号经积雪表面反射后被gps接收机接收，gps接收机可获得l1信号的信噪比观测值，即第一信噪比观测值。步骤s2中，gps接收机根据第一信噪比观测值及l1信号的波长，计算第一反射高度；示例性地，l1信号的波长为19.1cm。
[0059]
进一步地，gps接收机接收gps卫星或gnss卫星发射的l2c信号，得到l2c信号的第二信噪比观测值，其中，l2c信号的波长可以为24.4cm，根据第二信噪比观测值及l2c信号的波长，即可计算得到第二反射高度。
[0060]
需要说明的是，第一反射高度包括水文年每个年积日中l1信号经积雪表面反射后的反射高度，而由于l2c信号存在数据量少的问题，因此第二反射高度可能只包括部分水文年年积日中l2c信号经积雪表面反射后的反射高度。针对l2c信号缺测的这部分天数，本实施例对缺测天数对应的第一反射高度进行转换，得到l2c信号缺测天数对应的第三反射高度，这样第二反射高度和第三反射高度就包括水文年每个年积日中l2c信号经积雪表面反射后的反射高度。
[0061]
上述步骤s6～s7中，gps接收机根据第二反射高度和第三反射高度，确定水文年中每个年积日的积雪深度，并根据积雪深度估计水文年中每个年积日的雪水当量。
[0062]
可见，本发明通过对l1信号的第一反射高度和l2c信号的第二反射高度建立线性模型，将l2c反演精度与l1反演数据量进行有效结合，克服了现有技术中l1信号雪水当量监测精度低、而l2c信号数据量少的问题，不仅能够实现gps-ir积雪深度反演，也提高了雪水当量估计结果的准确性。
[0063]
可选地，上述步骤s2中，根据第一信噪比观测值和l1信号的波长，确定第一反射高度的步骤，包括：
[0064]
s201、根据第一信噪比观测值，获得第一信噪比观测值的残差序列；
[0065]
s202、对第一信噪比观测值的残差序列进行频谱分析，确定第一信噪比观测值的振荡频率；
[0066]
s203、根据l1信号的波长以及振荡频率，按照如下公式确定第一反射高度：
[0067][0068]
式中，f1表示第一信噪比观测值的振荡频率，λ1表示l1信号的波长，h
r1
表示第一反射高度；其中，第一反射高度包括水文年每个年积日中l1信号经积雪表面反射后的反射高度。
[0069]
本实施例中，首先根据第一信噪比观测值获得第一信噪比观测值的残差序列，接着对非等间隔采样的第一信噪比观测值的残差序列进行频谱分析，以确定第一信噪比观测值的振荡频率f1，而第一反射高度因此根据l1信号的波长和振荡频率f1可计算得到第一反射高度。
[0070]
同样地，对于l2c信号，需要获得第二信噪比观测值的残差序列，进而通过频谱分析确定第二信噪比观测值的振荡频率f2，并根据l2c信号的波长和振荡频率f2计算第二反射高度。由于上述过程与第一反射高度的计算过程相同，故此处不再赘述。
[0071]
可选地，上述步骤s5中，对将l2c信号缺测天数对应的第一反射高度进行转换，得到l2c信号缺测天数对应的第三反射高度的步骤，包括：
[0072]
获取l2c信号缺测天数对应的第一反射高度；
[0073]
确定第一线性拟合系数及第二线性拟合系数，利用第一线性拟合系数和第二线性拟合系数，将l2c信号缺测天数对应的第一反射高度转换得到l2c信号缺测天数对应的第三反射高度。
[0074]
上述步骤s5中，对于l2c信号的缺测天数，可以按照如下公式对l2c信号缺测天数对应的第一反射高度进行转换：
[0075]hr3
＝a
·h′
r1
b
[0076]
其中，a、b分别表示第一线性拟合系数和第二线性拟合系数，h
′
r1
表示l2c信号缺测天数对应的第一反射高度，h
r3
表示l2c信号缺测天数对应的第三反射高度。
[0077]
应当理解，通过对l1信号的第一反射高度和l2c信号的第二反射高度建立线性模型，可以利用缺测天数对应的第一反射高度转换得到l2c信号缺测天数对应的第三反射高度，第二反射高度与第三反射高度包括水文年每个年积日中l2c信号经积雪表面反射后的反射高度，因此本实施例将l2c反演精度与l1反演数据量进行有效结合，克服了现有技术中l1信号雪水当量监测精度低、而l2c信号数据量少的问题，
[0078]
可选地，根据第二反射高度和第三反射高度，确定水文年中每个年积日的积雪深度的步骤，包括：
[0079]
计算水文年中每个年积日对应的第二反射高度或第三反射高度与gps接收机自身高度的差值，得到水文年中每个年积日的积雪深度。
[0080]
如图2所示，若要计算水文年中某个年积日的积雪深度，则需要将该年积日对应的第二反射高度或第三反射高度与gps接收机的高度作差。由于本实施例是利用l2c信号的反射高度确定积雪深度，因此gps接收机自身高度可以为：无雪期内l2c信号每日反射高度的平均值h’。
[0081]
进一步地，上述步骤s7中，根据积雪深度，按照如下公式估计水文年中每个年积日的雪水当量：
[0082][0083]
式中，ρ0为雪密度常数，指雪季初期积雪深度极小时对应的雪密度，ρ
max
为水文年中雪密度的最大值，k1、k2为预设的经验系数，doyi为水文年中的第i个年积日，hi为水文年中的第i个年积日的积雪深度，为水文年中的第i个年积日的雪水当量。
[0084]
图3是本发明实施例提供的雪水当量估计结果的一种示意图。如图3所示，利用上述雪水当量估计方法基于l1信号及l2c信号进行雪水当量的估计，由图3(a)、(c)、(e)中雪水当量时间序列可以看出l1、l2c以及得到的雪水当量与snotel雪水当量有近似的变化趋势。显然，本发明提供的雪水当量估计方法相对反演效果较好的l2c雪深的增加了近25.8％的反演值，有效结合了l1的反演数据量以及l2c的反演精度，能够满足对于大部分应用对于
雪水当量估计的需求。
[0085]
图4是本发明实施例提供的基于gps双频信号的雪水当量估计装置的一种结构示意图。如图4所示，本发明实施例提供一种基于gps双频信号的雪水当量估计装置，应用于gps接收机；
[0086]
所述雪水当量估计方法包括：
[0087]
第一接收单元410，用于接收卫星发射的l1信号，并获得所述l1信号的第一信噪比观测值；
[0088]
第一确定单元420，用于根据所述第一信噪比观测值和所述l1信号的波长，确定第一反射高度；所述第一反射高度包括水文年每个年积日中l1信号经积雪表面反射后的反射高度；
[0089]
第二接收单元430，用于接收卫星发射的l2c信号，并获得所述l2c信号的第二信噪比观测值；
[0090]
第二确定单元440，用于根据所述第二信噪比观测值和所述l2c信号的波长，确定第二反射高度；所述第二反射高度包括至少部分水文年年积日中l2c信号经积雪表面反射后的反射高度；
[0091]
转换单元450，用于对将l2c信号缺测天数对应的第一反射高度进行转换，得到l2c信号缺测天数对应的第三反射高度；
[0092]
深度确定单元460，用于根据第二反射高度和第三反射高度，确定水文年中每个年积日的积雪深度；
[0093]
估计单元470，用于根据所述积雪深度，估计水文年中每个年积日的雪水当量。
[0094]
通过上述各实施例可知，本发明的有益效果在于：
[0095]
本发明提供了一种基于gps双频信号的雪水当量估计方法及装置，通过对l1信号的第一反射高度和l2c信号的第二反射高度建立线性模型，将l2c反演精度与l1反演数据量进行有效结合，克服了现有技术中l1信号雪水当量监测精度低、而l2c信号数据量少的问题，不仅能够实现gps-ir积雪深度反演，也提高了雪水当量估计结果的准确性。
[0096]
在本发明的描述中，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本发明的描述中，“多个”的含义是两个或两个以上，除非另有明确具体的限定。
[0097]
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。此外，本领域的技术人员可以将本说明书中描述的不同实施例或示例进行接合和组合。
[0098]
尽管在此结合各实施例对本技术进行了描述，然而，在实施所要求保护的本技术过程中，本领域技术人员通过查看所述附图、公开内容、以及所附权利要求书，可理解并实现所述公开实施例的其他变化。在权利要求中，“包括”(comprising)一词不排除其他组成部分或步骤，“一”或“一个”不排除多个的情况。单个处理器或其他单元可以实现权利要求
中列举的若干项功能。相互不同的从属权利要求中记载了某些措施，但这并不表示这些措施不能组合起来产生良好的效果。
[0099]
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。