[發明專利]一種利用GNSS大氣掩星彎曲角數據反演對流層頂參數的方法

在審

• 申請號: 202010185847.7
• 申請日: 2020-03-17
• 公開（公告）號: CN111323797A
• 公開（公告）日: 2020-06-23
• 發明（設計）人: 白偉華;孫越強;劉梓琰;杜起飛;劉黎軍;李偉;王先毅;蔡躍榮;曹光偉;夏俊明;孟祥廣;柳聰亮;趙丹陽;尹聰;胡鵬;王冬偉;劉成;吳春俊;李福;喬顥;程雙雙;朱光武
• 申請（專利權）人: 中國科學院國家空間科學中心
• 主分類號: G01S19/39
• 分類號: G01S19/39
• 代理公司: 北京方安思達知識產權代理有限公司 11472
• 代理人: 陳琳琳;楊青

1.一種利用GNSS大氣掩星彎曲角數據反演大氣對流層頂參數的方法，所述方法包括：

獲取GNSS掩星彎曲角數據并進行預處理；

對預處理后的GNSS掩星彎曲角數據進行截斷；

對截斷后的GNSS掩星彎曲角廓線取自然對數，使之成為近似線性的隨碰撞高度減小而增大的曲線；

將對數化GNSS掩星彎曲角進行協方差變換，得到協方差(CT)廓線；

尋找CT廓線的最大極大值，若該值大于第一閾值；則該點的碰撞高度為對流層頂高度，否則，尋找CT曲線的次大極大值，若該值大于第二閾值，則該點的碰撞高度為對流層頂高度，否則，放棄使用該CT廓線。

2.根據權利要求1所述的利用GNSS大氣掩星彎曲角數據反演對流層頂參數的方法，其特征在于，所述獲取GNSS掩星彎曲角數據并進行預處理；具體為：

獲取GNSS掩星彎曲角數據，從中選取GNSS掩星彎曲角廓線，要求其彎曲角與碰撞高度在5km-35km高度范圍內不存在無效值。

3.根據權利要求2所述的利用GNSS大氣掩星彎曲角數據反演對流層頂參數的方法，其特征在于，所述對GNSS掩星彎曲角數據進行截斷；具體為：截取35km以下的彎曲角廓線。

4.根據權利要求3所述的利用GNSS大氣掩星彎曲角數據反演對流層頂參數的方法，其特征在于，所述將對數化GNSS掩星彎曲角進行協方差變換，具體包括：

根據下述公式：

對自然對數彎曲角廓線做協方差變換，得到CT廓線；

其中，2a為協方差變換的寬度，z_b與z_t分別為自然對數彎曲角的下限與上限，f(z)為自然對數彎曲角廓線，CT(b)為CT廓線，b為該點對應的碰撞高度，h()是突出擾動的數學變換函數，z為變量。

5.根據權利要求4所述的利用GNSS大氣掩星彎曲角數據反演對流層頂參數的方法，其特征在于，所述第一閾值為1.05×最大極大值點上下5km內CT的平均值；所述第二閾值為1.05×次大極大值點上下5km內CT的平均值。

6.根據權利要求5所述的利用GNSS大氣掩星彎曲角數據反演對流層頂參數的方法，其特征在于，所述方法還包括：若CT廓線的最大極大值點滿足要求，則QC值標為1；若CT廓線的次大極大值點滿足要求，則QC值標為0。

7.根據權利要求1-6之一所述的利用GNSS大氣掩星彎曲角數據反演對流層頂參數的方法，其特征在于，所述方法還包括：當獲取到對流層頂高度，通過該值獲取對流層頂溫度。