[發明專利]一種高程異常內插推估方法在審
| 申請號： | 202010139946.1 | 申請日： | 2020-03-03 |
| 公開（公告）號： | CN111366130A | 公開（公告）日： | 2020-07-03 |
| 發明（設計）人： | 宋雷;趙碩;陳旭 | 申請（專利權）人： | 山東交通學院 |
| 主分類號： | G01C5/00 | 分類號： | G01C5/00 |
| 代理公司： | 北京匯捷知識產權代理事務所(普通合伙) 11531 | 代理人： | 馬金華 |
| 地址： | 250357 山東*** | 國省代碼： | 山東;37 |
| 權利要求書： | 查看更多 | 說明書： | 查看更多 |
| 摘要： | |||
| 搜索關鍵詞： | 一種 高程 異常 內插 方法 | ||
1.一種高程異常內插推估方法，其特征在于，包括以下步驟：
(1)擬合控制點的選取，以離散高程異常待定點為中心，在似大地水準面格網數值模型中選取與該點的大地經度差和大地緯度差分別小于預定數目格網間距的格網結點的大地坐標和高程異常；
所述預定數目格網間距的格網結點以兩個格網間距為條件選取獲得；
(2)點位投影，將選取的擬合控制點和高程異常待定點的大地坐標進行高斯投影，獲得點位的投影平面坐標；
(3)計算擬合參數，根據擬合控制點的投影平面坐標計算擬合區域的高程異常二次多項式擬合參數；
(4)計算待定點的高程異常，根據區域擬合參數和高程異常待定點的平面坐標計算該點的高程異常。
2.根據權利要求1所述的高程異常內插推估方法，其特征在于，所述步驟(1)包括：
設定離散高程異常待定點D的大地坐標為(BD,LD)，選取的擬合控制點大地坐標為(Bi,Li)，似大地水準面數值模型經、緯度格網間距為Δβ，則點位選取依據如下算法：
按照以上擬合控制點選取原則，在似大地水準面格網中取出16個格網節點高程異常作為擬合控制點。
3.根據權利要求1所述的高程異常內插推估方法，其特征在于，所述步驟(2)包括：
選用高程異常待定點實際經度作為投影的中央子午線，并且以中央子午線和赤道的交點o作為坐標原點，以中央子午線的投影為縱坐標x軸，以赤道的投影為橫坐標y軸,形成高斯平面直角坐標系。
4.根據權利要求3所述的高程異常內插推估方法，其特征在于，在點位投影過程中，地球橢球的長半軸：a＝6378137m；地球扁率：1/298.257 222 101，地球自轉角速度ω＝7.292115×10-5rad/s；地心引力常數GM＝3.986004418×1014m3/s2；
點位投影就是根據橢球面上A點的大地坐標(B,L)，求該點在高斯投影平面上的直角坐標(x,y)，即的坐標變換，精度精確至0.5mm時的公式為:
式中:X為子午線弧長；l為橢球面上P點的經度與中央子午線L0的經度差：l＝L-L0,P點在中央子午線之東,l為正，在西則為負；N為地球橢球的卯酉圈曲率半徑,t和η為簡記符號,式(2)中其它量的計算公式為：
t＝tanB
η2＝e′2cos2B
在式(2)中，設有子午線上兩點O和H，O在赤道上，H的緯度為B，a為橢球長半軸，b為橢球短半軸,O和H間的子午線弧長X計算公式
式(3)中，參數的計算公式為：
5.根據權利要求1所述的高程異常內插推估方法，其特征在于，所述步驟(3)包括：
假定在高程異常待定點鄰近區域的高程異常ζ與坐標(x,y)存在如下函數關系：
式中a0,a1,a2,a3,a4,a5為多項式系數,
利用擬合控制點的投影平面坐標和高程異常組成多項式擬合方程，16個擬合控制點組成16個方程的方程組，根據最小二乘原理求解6個多項式系數的值。
6.根據權利要求1所述的高程異常內插推估方法，其特征在于，所述步驟(4)包括：
D點高程異常根據經緯度由多項式函數可求得
該專利技術資料僅供研究查看技術是否侵權等信息，商用須獲得專利權人授權。該專利全部權利屬于山東交通學院，未經山東交通學院許可，擅自商用是侵權行為。如果您想購買此專利、獲得商業授權和技術合作，請聯系【客服】
本文鏈接：http://www.szxzyx.cn/pat/books/202010139946.1/1.html，轉載請聲明來源鉆瓜專利網。
