• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    向下延拓航空重力數(shù)據(jù)的Tikhonov雙參數(shù)正則化法

    2011-01-04 07:58:40鄧凱亮黃謨濤暴景陽歐陽永忠陸秀平劉傳勇
    測繪學(xué)報 2011年6期
    關(guān)鍵詞:水準(zhǔn)面正則重力

    鄧凱亮,黃謨濤,暴景陽,歐陽永忠,陸秀平,劉傳勇

    1.大連艦艇學(xué)院 海測工程系,遼寧大連116018;2.海軍海洋測繪研究所,天津300061

    向下延拓航空重力數(shù)據(jù)的Tikhonov雙參數(shù)正則化法

    鄧凱亮1,2,黃謨濤2,暴景陽1,歐陽永忠2,陸秀平2,劉傳勇2

    1.大連艦艇學(xué)院 海測工程系,遼寧大連116018;2.海軍海洋測繪研究所,天津300061

    為避免正則化參數(shù)對向下延拓過程可靠成分的修正影響,提出Tikhonov雙參數(shù)正則化法。引進(jìn)截斷參數(shù),將法矩陣的奇異值分為相對較大的奇異值(可靠部分)和相對較小的奇異值(不可靠部分);引進(jìn)正則化參數(shù),只對法矩陣的小奇異值進(jìn)行修正,以抑制高頻誤差對向下延拓解的影響。采用廣義交互確認(rèn)法(GCV)確定截斷參數(shù)和正則化參數(shù)?;贓GM2008重力場模型仿真一組航空重力數(shù)據(jù),驗證該方法對航空重力數(shù)據(jù)向下延拓過程的有效性。

    向下延拓;航空重力數(shù)據(jù);Tikhonov雙參數(shù)正則化法;GCV法

    1 引 言

    航空重力測量不僅能快速經(jīng)濟(jì)地獲得中高頻重力場信息,而且能遠(yuǎn)距離測量,有效地突破測量范圍的局限性,是衛(wèi)星重力測量和海面重力測量的重要補充[1-2]。

    在大地水準(zhǔn)面等應(yīng)用上,需要將航空重力數(shù)據(jù)向下延拓到大地水準(zhǔn)面上,并采用球近似。向下延拓的基本方法是逆Poisson方法,這是一個信號放大的非平穩(wěn)過程,很小的觀測噪聲往往引起待估參數(shù)較大的誤差,屬于不適定問題[2]。對此,國內(nèi)外學(xué)者作了許多有益的研究[3-10],提出迭代法[3]、配置法[4]、倒錐法[5]、直接代表法[6]、Tikhonov正則化法[7-10]等;文獻(xiàn)[7—8]對上述方法在航空重力測量數(shù)據(jù)向下延拓中的應(yīng)用進(jìn)行了比較分析,表明造成向下延拓結(jié)果不適定的主要原因是觀測誤差的高頻部分被法矩陣的小奇異值嚴(yán)重放大。Tikhonov正則化方法通過正則化參數(shù)對法矩陣小奇異值的修正,能有效抑制高頻誤差對參數(shù)估值的影響,同時驗證了該方法的優(yōu)越性。但是正則化參數(shù)不僅修正了較小的奇異值,而且對相對較大的奇異值也進(jìn)行了修正,這樣使模型中的可靠成分也受到了畸變的影響[11-12]。

    為了避免正則化參數(shù)對向下延拓過程可靠成分的修正影響,筆者提出Tikhonov雙參數(shù)正則化法?;舅悸肥牵涸谙蛳卵油刈V分解的基礎(chǔ)上,引進(jìn)截斷參數(shù),將法矩陣的奇異值分為相對較大的奇異值(可靠部分)和相對較小的奇異值(不可靠部分);引進(jìn)正則化參數(shù),只對法矩陣的小奇異值進(jìn)行修正,以抑制高頻誤差對向下延拓解的影響。參數(shù)的選擇采用廣義交互確認(rèn)法(generalized cross-validation,GCV)?;贓GM2008重力場模型仿真了一組航空重力數(shù)據(jù),以驗證該方法對航空重力數(shù)據(jù)向下延拓過程的有效性。

    2 向下延拓原理

    根據(jù)位理論和Dirichlet邊值理論,重力異常的向上延拓可由Poisson積分公式[3]計算,公式為

    式中,R為地球半徑;h是延拓的高度;r為延拓高度處的球半徑(r=R+h);ω為以R為半徑的球面;K(r,ψ,R)是積分核函數(shù),其公式[3]為

    式中,ψ是計算點與球面流動點之間的角距,其表達(dá)式為

    一般情況下,航空重力異常和大地水準(zhǔn)面上的重力異常均是離散值,故將式(1)離散化,表達(dá)式可化為

    式中,ωc是有效積分范圍;N為積分范圍ωc內(nèi)的測量點數(shù);Δgω-ωc(r,φ,λ)為遠(yuǎn)區(qū)影響,引起的誤差為截斷誤差,可由重力場模型近似計算。

    則式(3)可寫為

    當(dāng)航空重力異常Δgh(r)向下延拓到大地水準(zhǔn)面時,Δgh(r)是已知觀測量,大地水準(zhǔn)面上的重力異常Δg(R,θj,λj)是待求量,這是Poisson積分的逆問題。

    積分區(qū)域內(nèi)大地水準(zhǔn)面上的重力異常的最小二乘解為

    3 Tikhonov雙參數(shù)正則化法

    航空重力測量數(shù)據(jù)的向下延拓是信號放大的不適定過程,很小的觀測噪聲往往能引起待估參數(shù)較大的誤差,這正是正則化法所要解決的問題。在眾多的正則化方法中,以Tikhonov正則化法應(yīng)用最為廣泛。

    3.1 向下延拓模型的譜分解

    在航空測量數(shù)據(jù)的向下延拓中,觀測方程系數(shù)矩陣的奇異值單調(diào)遞減,最小二乘解被放大的觀測噪聲污染。對系數(shù)矩陣A作奇異值分解(SVD)

    式中,U、V分別是A的左右奇異向量,U=[u1u2… un],V=[v1v2… vn],且滿足是對角矩陣,其對角線上的元素為A的遞減的奇異值λi。將式(6)代入式(5)中得最小二乘解的譜分解形式

    3.2 Tikhonov雙參數(shù)正則化法

    Tikhonov正則化法的實質(zhì)是尋找適當(dāng)?shù)臑V波因子,用相鄰的適定解去逼近原問題的解[14]。對于觀測方程式(4)而言,其正則化函數(shù)為

    式中,α>0是正則化參數(shù);‖x‖表示x的范數(shù)。根據(jù)式(9)的約束條件,可得Tikhonov正則化解

    Tikhonov正則化解的濾波因子為

    將式(10)進(jìn)行譜分解,得到

    由式(12)、式(13)可以看出,正則化參數(shù)不僅對小的奇異值進(jìn)行了修正,也對相對較大的奇異值即模型的可靠部分進(jìn)行了修正,使得最小二乘解發(fā)生了畸變。

    對此,設(shè)計了Tikhonov雙參數(shù)正則化法。主要思路是:引進(jìn)Tikhonov截斷參數(shù),將奇異值分為相對較大的奇異值和相對較小的奇異值,即可靠部分和不可靠部分;沿用正則化參數(shù),使其只對模型的不可靠部分進(jìn)行修正。這樣既抑制了小奇異值對觀測噪聲高頻段的放大影響,又避免模型的可靠部分受到修正影響。

    Tikhonov雙參數(shù)正則化法的濾波因子為

    式中,k為截斷參數(shù)。

    Tikhonov雙參數(shù)正則化法的譜分解形式為

    Tikhonov雙參數(shù)正則化解的均方誤差為

    由式(16)可以看出,Tikhonov雙參數(shù)正則化法的估值誤差包括兩部分:前者是測量誤差引起的估值誤差,隨著截斷參數(shù)k的增大而減小,隨著正則化參數(shù)α的增大而減??;后者是正則化引起的估值誤差,隨著截斷參數(shù)k的增大而減小,隨著正則化參數(shù)α的增大而增大。因此,存在最佳的截斷參數(shù)k和正則化參數(shù)α,使得由測量誤差與正則化誤差引起的估值誤差達(dá)到最佳的平衡。

    3.3 正則化參數(shù)的選取

    近年來,統(tǒng)計學(xué)界和大地測量學(xué)界提出了多種方法選擇正則化參數(shù)。廣義交互確認(rèn)法(generalized cross-validation,GCV)得到了廣泛的應(yīng)用[15-16],該法是基于“遺漏一個”(leave-out-one)的思想。該思想轉(zhuǎn)化為一個便于操作的最小化問題時,就是選擇合適正則化參數(shù)GCV函數(shù)最小。該函數(shù)定義為

    Qα是所謂的影響矩陣,由Axα=QαL 定義。Wahba給出了求得GCV函數(shù)的另一種形式

    式中,n為觀測值個數(shù);tr為矩陣的跡。最佳的截斷參數(shù)k和正則化參數(shù)α對應(yīng)于GCV函數(shù)的最小值。

    對于Tikhonov雙參數(shù)正則化法,基于Qα的定義,容易得到

    式中

    式中,VN、UN由Poisson方程的法矩陣的譜分解式N=UNdiag(λN)VTN得到;λN(k,α)由式(20)得到

    在數(shù)據(jù)量大的情況下執(zhí)行GCV,應(yīng)考慮正則化參數(shù)的收斂速度[17]:在得到法矩陣的奇異值λN后,假設(shè)正則化參數(shù)為q個,則取αq=λn,ratio=,修正參數(shù)為

    在本文的仿真試驗中,法矩陣的奇異值λN分布圖見圖1,所考慮的正則化參數(shù)α的分布圖見圖2。

    圖1 法矩陣的奇異值分布圖Fig.1 Singular values characteristic of the Poisson normal matrix

    圖2 正則化參數(shù)分布圖Fig.2 Characteristic regularization parameters

    4 仿真試驗及分析

    為了驗證基于Tikhonov雙參數(shù)正則化法對航空重力數(shù)據(jù)進(jìn)行向下延拓結(jié)果的有效性,設(shè)計了航行高度為5km的仿真試驗。

    4.1 數(shù)據(jù)準(zhǔn)備

    基于EGM2008重力場模型計算的重力異常作為仿真試驗的數(shù)據(jù)基礎(chǔ)。EGM2008重力場模型[18]是由NGA(national geospatial-intelligence agency)釋放的全球超高階地球重力場模型,由衛(wèi)星重力測量、衛(wèi)星測高和大地水準(zhǔn)面重力觀測等資料聯(lián)合解算得到,模型階數(shù)達(dá)到2 160。EGM2008重力場模型導(dǎo)出的重力異常在我國大陸的總體精度為10.5mGal[19](1Gal=0.01m/s2)。

    基于EGM2008重力場模型計算某區(qū)域2 160階的大地水準(zhǔn)面上的重力異常Δg0(圖3)和5km高空面的重力異常Δg5(圖4)。它的范圍為1°× 1°,格網(wǎng)間距為2′×2′。仿效移去-恢復(fù)技術(shù)的應(yīng)用,取EGM2008重力場模型的360階作為參考模型,得到參考重力異常Δg′0和參考航空重力異常Δg′5。各項重力異常的統(tǒng)計特性見表1。

    圖3 仿真的大地水準(zhǔn)面重力異常Δg0Fig.3 Simulative ground gravity dataΔg0

    圖4 仿真的航空重力異常Δg5Fig.4 Simulative airborne gravity dataΔg5

    表1 仿真區(qū)域各項重力異常特性的統(tǒng)計Tab.1 Statistics of gravity anomalies at simulation area mGal

    4.2 試驗步驟

    為了模擬航空重力測量數(shù)據(jù)的處理結(jié)果的誤差,對仿真的航空重力異常Δg5引入三種觀測誤差分別是零均值的白噪聲:e1(σ=±1mGal),e2(σ=±3mGal),e3(σ=±5mGal)。試驗步驟如下:

    (1)模擬產(chǎn)生3個量級的誤差并加入Δg5中,得到含誤差的仿真航空重力異常Δg5e1、Δg5e2、Δg5e3,并從中移去參考航空重力異常Δg′5,得到殘差航空重力異常δΔg5e1、δΔg5e2、δΔg5e3。

    (2)基于Tikhonov雙參數(shù)正則化法,利用GCV法,求得各殘差航空重力異常δΔg5e1、δΔg5e2、δΔg5e3向下延拓時的截斷參數(shù)k和正則化參數(shù)α(圖5)。為了驗證雙參數(shù)法的效果,設(shè)計了Tikhonov法的兩種方案,方案1是只有正則化參數(shù)的Tikhonov法,方案2是Tikhonov雙參數(shù)正則化法。

    (3)按式(12)對各殘差航空重力異常δΔg5e1、 δΔg5e2、δΔg5e3向下延拓,得到延拓后的殘差大地水準(zhǔn)面重力異常δΔg0e1、δΔg0e2、δΔg0e3,并恢復(fù)參考重力異常Δg′0,得到大地水準(zhǔn)面重力異常Δg0e1、Δg0e2、Δg0e3。

    圖5 各誤差條件下兩種方案確定的正則化參數(shù)示意圖Fig.5 Regularization parameters given by two projects each under the conditions of three error levels

    4.3 結(jié)果比較與分析

    為了分析Tikhonov雙參數(shù)正則化法對航空重力數(shù)據(jù)基于逆Poisson積分方程的向下延拓模型的有效性,將向下延拓得到的大Δg0e1、Δg0e2、Δg0e3與仿真的大地水準(zhǔn)面重力異常Δg0進(jìn)行比較,以檢驗向下延拓后數(shù)據(jù)的精度。

    為了驗證本方法的效果,設(shè)計了三種計算方法:方法1,直接逆Poisson積分方程的向下延拓法;方法2,基于GCV法的只有正則化參數(shù)的Tikhonov法的向下延拓法;方法3,基于GCV法的Tikhonov雙參數(shù)正則化法的向下延拓法。

    基于三種計算方法獲得的大地水準(zhǔn)面重力異常與原始海面重力異常比較的差值統(tǒng)計結(jié)果見表2。

    表2 與模擬的重力異常真值Δg0比較的差值統(tǒng)計Tab.2 Statistics of the differences between calculated values(Δg0e1,Δg0e2,Δg0e3)and the simulative true valuesΔg0in the comparison

    基于三種計算方法獲得的大地水準(zhǔn)面上的重力異常見圖6。

    由表2可以看出:

    (1)方法1的法矩陣的條件數(shù)達(dá)到了14 886.62,遠(yuǎn)遠(yuǎn)大于1 000的病態(tài)界點,表明此時的法矩陣是嚴(yán)重病態(tài)的。

    (2)在設(shè)計的三種誤差條件下,方法2使得法矩陣的條件數(shù)控制在70以內(nèi),表明此時的向下延拓結(jié)果具有良好的穩(wěn)定性;方法3將法矩陣的條件數(shù)控制在30以內(nèi),表明該方法能有效抑制微小誤差被小奇異值放大程度,在保證向下延拓結(jié)果的穩(wěn)定性方面具有優(yōu)越性。

    (3)在設(shè)計的三種誤差條件下,采用方法2,標(biāo)準(zhǔn)差分別由原來的20.28mGal、54.53mGal、83.94mGal降低到4.36mGal、6.45mGal、 8.18mGal,而方法3,能進(jìn)一步提高成果的精度,分別降到了3.27mGal、4.29mGal和6.08mGal。表明方法2和方法3都可有效地抑制觀測誤差對向下延拓結(jié)果的影響,提高成果的精度,后者的改善效果更加明顯。

    圖6 各誤差條件下三種方法的向下延拓成果Fig.6 Downward continuation results given by three methods each under the conditions of three error levels

    5 結(jié)束語

    利用直接逆Poisson積分方程進(jìn)行航空重力數(shù)據(jù)向下延拓,延拓結(jié)果極不穩(wěn)定(用條件數(shù)表示)。在Tikhonov正則化處理過程中,基于GCV法得到與誤差環(huán)境相匹配的正則化參數(shù),能抑制小奇異值對觀測誤差的放大影響,改進(jìn)延拓結(jié)果的穩(wěn)定性和精度。但是此時正則化參數(shù)對法矩陣的可靠部分即相對大的奇異值也進(jìn)行了改正,這樣將導(dǎo)致模型的可靠部分發(fā)生畸變。

    在對向下延拓過程法矩陣進(jìn)行譜分解的基礎(chǔ)上,為避免正則化參數(shù)對法矩陣奇異值可靠成分的修正影響,提出Tikhonov雙參數(shù)正則化法。即在正則化參數(shù)的基礎(chǔ)上,引進(jìn)截斷參數(shù),將法矩陣的奇異值分為相對較大的奇異值(可靠部分)和相對較小的奇異值(不可靠部分),而正則化參數(shù)只對法矩陣的不可靠部分進(jìn)行修正,以抑制高頻誤差對向下延拓解的影響。

    Tikhonov雙參數(shù)正則化法的關(guān)鍵是獲取最優(yōu)的截斷參數(shù)和正則化參數(shù)。采用GCV函數(shù)作為判斷標(biāo)準(zhǔn),GCV函數(shù)隨著截斷參數(shù)和正則化參數(shù)的同時引進(jìn)而作了相應(yīng)的改進(jìn)。通過基于EGM2008重力場模型仿真的航空重力數(shù)據(jù)進(jìn)行向下延拓試驗,驗證了Tikhonov雙參數(shù)正則化法要優(yōu)于僅有正則化參數(shù)的Tikhonov正則化法。

    [1] SUN Zhongmiao.Theory,Methods and Applications of Airborne Gravimetry[D].Zhengzhou:Information Engineering University,2004.(孫中苗.航空重力測量理論、方法及應(yīng)用研究[D].鄭州:信息工程大學(xué),2004.)

    [2] KERN M.An Analysis of the Combination and Downward Continuation of Satellite,Airborne and Terrestrial Gravity Data[D].Calgary:University of Calgary,2003.

    [3] HEISKANEN W A,MORITZ H.Physical Geodesy[M].Beijing:Surveying and Mapping Press,1979.(海斯卡涅W A,莫里斯H.物理大地測量[M].北京:測繪出版社,1979.)

    [4] KELLER W,HIRSCH M.Downward Continuation Versus Free-air-reduction in Airborne Gravimetry[C]∥Geodesy and Physics of the Earth:Geodetic Contributions to Geodynamics.Berlin:Springer-Verlag,1992:266-272.

    [5] YANG Qiangwen,WU Xiaoping.“Qusai-cone Method”,a New Way of Downward Continuation of Airborne Gravimetry Data[J].Journal of the PLA Institute of Surveying and Mapping,1998,15(3):169-172.(楊強文,吳曉平.航空重力數(shù)據(jù)向下延拓的倒錐法[J].解放軍測繪學(xué)院學(xué)報,1998,15(3):169-172.)

    [6] SHI Pan,WANG Xingtao.The Frequency Domain Analysis of the Determination of Terrestrial Mean Gravity Anomaly Using Airborne Gravimetry[J].Acta Geodaetica et Cartographica Sinica,1995,24(4):301-308.(石磐,王興濤.空中測量地面平均重力異常的頻域分析[J].測繪學(xué)報,1995,24(4):301-308.)

    [7] WANG Xingtao,XIA Zheren,SHI Pan,et al,A Comparison of Different Downward Continuation Methods for Airborne Gravity Data[J].Chinese Journal of Geophysics,2004,47(6):1017-1021.(王興濤,夏哲仁,石磐,等.航空重力測量數(shù)據(jù)向下延拓方法比較[J].地球物理學(xué)報,2004,47(6):1017-1021.)

    [8] WANG Xingtao,SHI Pan,ZHU Feizhou.Regularization Methods and Spectral Decomposition for the Downward Continuation of Airborne Gravity Data[J].Acta Geodaetica et Cartographica Sinica,2004,33(1):33-38.(王興濤,石磐,朱非洲.航空重力測量數(shù)據(jù)向下延拓的正則化算法及其譜分解[J].測繪學(xué)報,2004,33(1):33-38.)

    [9] HAO Yanling,CHENG Yi,SUN Feng,et al.Simulation Research on Tikhonov Regularization Algorithm in Downward Continuation[J].Chinese Journal of Scientific Instrument,2008,29(8):2190-2194.(郝燕玲,成怡,孫楓,等.Tikhonov正則化向下延拓算法仿真實驗研究[J].儀器儀表學(xué)報,2008,29(8):2190-2194.)

    [10] CHENG Yi,HAO Yanling,LIU Fanming.Analysis of Effect and Downward Continuation of Airborne Gravity Data[J].Journal of System Simulation,2008,20(3):605-609.(成怡,郝燕玲,劉繁明.航空重力測量數(shù)據(jù)向下延拓及其影響因素分析[J].系統(tǒng)仿真學(xué)報,2008,20(3):605-609.)

    [11] WANG Zhenjie.Research on the Regularization Solutions of Ill-Posed Problems in Geodesy[D].Wuhan:Institute of Geodesy and Geophysics,Chinese Academy of Sciences,2003.(王振杰.大地測量中不適定問題的正則化解法研究[D].武漢:中國科學(xué)院測量與地球物理研究所,2003.)

    [12] DENG Kailiang.Research on the Procession,Combination and Application of the Muti-source Gravity Data on the Sea[D].DaLian:Dalian Naval Academy,2011.(鄧凱亮.海域多源重力數(shù)據(jù)的處理、融合及應(yīng)用研究[D].大連:海軍大連艦艇學(xué)院,2011.)

    [13] SHEN Yunzhong,XU Houze.Spectral Decomposition Formula of Regularization Solution for Ill-posed Equation[J].Journal of Geodesy and Geodynamics,2002,33(3):11-14.(沈云中,許厚澤.不適定方程正則化算法的譜分解式[J].大地測量學(xué)與地球動力學(xué),2002,33(3):11-14.)

    [14] LIU Jijun.Regularization Method and Application of the Ill-posed Equation[M].Beijing:Science Press,2005.(劉繼軍.不適定問題的正則化方法及應(yīng)用[M].北京:科學(xué)出版社,2005.)

    [15] KUSCHE J,KLEES R.Regularization of Gravity Field Estimation from Satellite Gravity Gradients[J].Journal of Geodesy,2002,76(6-7):359-368.

    [16] LIU Xianglin,DITMAR P.Smoothing a Satellite Orbit on the Basis of B-spline and Regularization[J].Chinese Journal of Geophysics,2006,49(1):99-105.(柳響林,DITMAR P.基于B-spline和正則化算法的低軌衛(wèi)星軌道平滑[J].地球物理學(xué)報,2006,49(1):99-105.)

    [17] GOLUB G H,VON MATT U.Generalized Cross-validation for Large-scale Problems[J].Journal of Computational and Graphical Statistics,1997,6(1):1-34.

    [18] PAVLIS N K,HOLMES S A,KENYON S C,et al.An Earth Gravitational Model to Degree 2160:EGM2008[J].Geophysical Research Abstracts,2008,10:13-18.

    [19] ZHANG Chuanyin,GUO Chunxi,CHEN Junyong,et al.EGM2008and Its Application Analysis in Chinese Mainland[J].Acta Geodaetica et Cartographica Sinica,2009,38(4):283-289.(章傳銀,郭春喜,陳俊勇,等.EGM2008地球重力場模型在中國大陸適用性分析[J].測繪學(xué)報,2009, 38(4):283-289.)

    Tikhonov Two-parameter Regulation Algorithm in Downward Continuation of Airborne Gravity Data

    DENG Kailiang1,2,HUANG Motao2,BAO Jingyang1,OUYANG Yongzhong2,LU Xiuping2,LIU Chuanyong2
    1.Department of Hydrography and Cartography,Dalian Naval Academy,Dalian116018,China;2.Naval Institute of Hydrographic Surveying and Charting,Tianjin 300061,China

    In order to avoid the correct effect provided by regulation parameter on the downward continuation credible portion,Tikhonov two-parameter regulation algorithm is presented.Using the truncation parameter,the singular values of the normal matrix are divided into the relative large singular values(the credible portion)and the relative small singular values(the incredible portion).Using regulation parameter,only the relative small singular values are corrected,and the effect induced by the high frequency error is suppressed.The truncation parameter and the regulation parameter are determined by GCV method.Using the simulative airborne gravity data based on the EGM2008,as true values of the gravity field,the results of numerical calculation examples have clearly demonstrated that this method for downward continuation of airborne gravity data is fairly efficient,and shown significant advantages.

    downward continuation;airborne gravity data;Tikhonov two-parameter regulation algorithm;GCV method

    DENG Kailiang(1983—),male,PhD candidate,majors in marine geodetic survey.

    1001-1595(2011)06-0690-07

    P223

    A

    國家863計劃(2006AA06A202;2009AA121405);國家自然科學(xué)基金(41074002;61071006);國家重大海洋勘測專項資助(4200301)

    叢樹平)

    2010-10-11

    2011-03-17

    鄧凱亮(1983—),男,博士生,研究方向為海洋大地測量。

    E-mail:dengkailiang036@163.com

    猜你喜歡
    水準(zhǔn)面正則重力
    瘋狂過山車——重力是什么
    剩余有限Minimax可解群的4階正則自同構(gòu)
    類似于VNL環(huán)的環(huán)
    仰斜式重力擋土墻穩(wěn)定計算復(fù)核
    一張紙的承重力有多大?
    GPS似大地水準(zhǔn)面精化及精度分析
    有限秩的可解群的正則自同構(gòu)
    重力異常向上延拓中Poisson積分離散化方法比較
    顧及完全球面布格異常梯度項改正的我國似大地水準(zhǔn)面精化
    基于自適應(yīng)最小二乘配置的區(qū)域似大地水準(zhǔn)面擬合
    欧美精品一区二区免费开放| 国产成人免费观看mmmm| 性色av乱码一区二区三区2| 国产伦理片在线播放av一区| 亚洲精品国产av蜜桃| 一边摸一边做爽爽视频免费| 欧美黄色淫秽网站| 天天躁夜夜躁狠狠久久av| 五月开心婷婷网| 国产伦理片在线播放av一区| 亚洲国产看品久久| 嫩草影视91久久| 在线观看www视频免费| 下体分泌物呈黄色| 一本—道久久a久久精品蜜桃钙片| 蜜桃在线观看..| 不卡av一区二区三区| 亚洲人成电影免费在线| 国产又色又爽无遮挡免| 精品卡一卡二卡四卡免费| 香蕉国产在线看| 精品国产一区二区三区四区第35| 黄色怎么调成土黄色| 十八禁人妻一区二区| av线在线观看网站| 青青草视频在线视频观看| 欧美精品人与动牲交sv欧美| 国产成人系列免费观看| 男人爽女人下面视频在线观看| 亚洲欧美一区二区三区久久| 久久久久精品国产欧美久久久 | av在线app专区| 欧美日韩av久久| 不卡av一区二区三区| 最近手机中文字幕大全| 99久久99久久久精品蜜桃| 午夜免费鲁丝| 人人澡人人妻人| 中国国产av一级| 美女扒开内裤让男人捅视频| xxx大片免费视频| 少妇人妻 视频| 亚洲av成人精品一二三区| 美女主播在线视频| 久9热在线精品视频| 1024视频免费在线观看| 99国产精品一区二区蜜桃av | 亚洲精品一二三| 精品人妻在线不人妻| 赤兔流量卡办理| 国产野战对白在线观看| 国产真人三级小视频在线观看| 国产伦理片在线播放av一区| 日本欧美视频一区| 久久久国产欧美日韩av| 人妻 亚洲 视频| 亚洲国产看品久久| 国产午夜精品一二区理论片| 日韩 欧美 亚洲 中文字幕| av欧美777| 亚洲国产毛片av蜜桃av| 久久久国产欧美日韩av| 一级a爱视频在线免费观看| 一本大道久久a久久精品| 五月天丁香电影| 中文乱码字字幕精品一区二区三区| 2021少妇久久久久久久久久久| 国产黄色视频一区二区在线观看| 热99国产精品久久久久久7| 日韩 亚洲 欧美在线| 国产黄色视频一区二区在线观看| 一级黄色大片毛片| 天堂俺去俺来也www色官网| 精品少妇一区二区三区视频日本电影| 一级黄色大片毛片| 2021少妇久久久久久久久久久| 黄色视频在线播放观看不卡| 曰老女人黄片| 97人妻天天添夜夜摸| 91九色精品人成在线观看| 欧美日韩国产mv在线观看视频| 高清不卡的av网站| 脱女人内裤的视频| 一本久久精品| 国产视频一区二区在线看| 韩国高清视频一区二区三区| 亚洲精品在线美女| 999精品在线视频| 三上悠亚av全集在线观看| 亚洲精品中文字幕在线视频| 999久久久国产精品视频| 老熟女久久久| 国产精品麻豆人妻色哟哟久久| 亚洲伊人色综图| 无遮挡黄片免费观看| 国产在线观看jvid| 两个人免费观看高清视频| 男女床上黄色一级片免费看| www.av在线官网国产| 91精品三级在线观看| 91精品三级在线观看| 亚洲九九香蕉| 大型av网站在线播放| 1024视频免费在线观看| 欧美日韩福利视频一区二区| 日本wwww免费看| 国产男人的电影天堂91| 伦理电影免费视频| 久久久久久久精品精品| 国产免费视频播放在线视频| 女人精品久久久久毛片| 99热网站在线观看| 亚洲免费av在线视频| 只有这里有精品99| 国产成人影院久久av| 亚洲精品美女久久久久99蜜臀 | 一区二区av电影网| 国产精品三级大全| 亚洲精品久久午夜乱码| 丝袜美足系列| 久久人人爽人人片av| 久久久久久久国产电影| 精品国产一区二区三区四区第35| 欧美黑人精品巨大| 黄色一级大片看看| 成人亚洲欧美一区二区av| 久久久亚洲精品成人影院| 18禁黄网站禁片午夜丰满| 欧美日韩亚洲国产一区二区在线观看 | 久久久久久亚洲精品国产蜜桃av| 国产高清videossex| 色播在线永久视频| 一本久久精品| 国产精品一二三区在线看| 日韩中文字幕欧美一区二区 | 99国产精品一区二区蜜桃av | 母亲3免费完整高清在线观看| 欧美日韩视频精品一区| 亚洲免费av在线视频| 国产精品 欧美亚洲| 青青草视频在线视频观看| 亚洲成人手机| 日本欧美国产在线视频| 日韩中文字幕欧美一区二区 | 欧美精品一区二区免费开放| 香蕉国产在线看| 一级毛片 在线播放| 亚洲欧洲日产国产| 十分钟在线观看高清视频www| 欧美国产精品一级二级三级| 我的亚洲天堂| 高清不卡的av网站| 热99国产精品久久久久久7| 建设人人有责人人尽责人人享有的| 少妇粗大呻吟视频| 免费少妇av软件| 91精品三级在线观看| 香蕉丝袜av| 婷婷色综合大香蕉| 两人在一起打扑克的视频| 中国国产av一级| 国产免费视频播放在线视频| 欧美精品亚洲一区二区| 狠狠精品人妻久久久久久综合| 极品少妇高潮喷水抽搐| 亚洲国产精品999| 亚洲久久久国产精品| 欧美xxⅹ黑人| 免费人妻精品一区二区三区视频| 亚洲国产av影院在线观看| 激情五月婷婷亚洲| av国产精品久久久久影院| 精品国产国语对白av| 亚洲精品第二区| 91字幕亚洲| 国产精品麻豆人妻色哟哟久久| 男女床上黄色一级片免费看| 久久久久久久久久久久大奶| 少妇人妻 视频| 首页视频小说图片口味搜索 | 18禁观看日本| 欧美av亚洲av综合av国产av| 国产伦理片在线播放av一区| 美女脱内裤让男人舔精品视频| 久久九九热精品免费| 18禁裸乳无遮挡动漫免费视频| 搡老乐熟女国产| 不卡av一区二区三区| 久久久久久人人人人人| 日韩视频在线欧美| 久久国产精品男人的天堂亚洲| 深夜精品福利| 久久人妻熟女aⅴ| 亚洲色图 男人天堂 中文字幕| 99久久精品国产亚洲精品| 日本av手机在线免费观看| 亚洲中文日韩欧美视频| 妹子高潮喷水视频| 成年动漫av网址| 老汉色av国产亚洲站长工具| 国产成人精品无人区| 成年人黄色毛片网站| www.精华液| 9191精品国产免费久久| 国产视频一区二区在线看| 亚洲成人免费av在线播放| 99国产精品免费福利视频| 色网站视频免费| 美女主播在线视频| 国产精品人妻久久久影院| 日日摸夜夜添夜夜爱| 欧美 亚洲 国产 日韩一| 午夜视频精品福利| 两个人免费观看高清视频| 99国产综合亚洲精品| 女警被强在线播放| √禁漫天堂资源中文www| 中文字幕av电影在线播放| 亚洲精品一卡2卡三卡4卡5卡 | svipshipincom国产片| 男男h啪啪无遮挡| 在线观看www视频免费| 亚洲第一av免费看| 日韩熟女老妇一区二区性免费视频| 亚洲精品自拍成人| 国产成人欧美在线观看 | 青青草视频在线视频观看| 国产av国产精品国产| 亚洲熟女毛片儿| 亚洲精品国产av蜜桃| 另类亚洲欧美激情| 国产女主播在线喷水免费视频网站| av片东京热男人的天堂| av国产精品久久久久影院| 亚洲精品一卡2卡三卡4卡5卡 | 国产免费一区二区三区四区乱码| 午夜视频精品福利| 黄色一级大片看看| 在线观看免费日韩欧美大片| 国产视频首页在线观看| 水蜜桃什么品种好| 日韩欧美一区视频在线观看| 国产熟女欧美一区二区| 女人精品久久久久毛片| 91九色精品人成在线观看| 91字幕亚洲| 老汉色av国产亚洲站长工具| 波多野结衣一区麻豆| 人人妻人人添人人爽欧美一区卜| 亚洲熟女毛片儿| 尾随美女入室| 亚洲一码二码三码区别大吗| 大码成人一级视频| 国产精品熟女久久久久浪| 丝瓜视频免费看黄片| 亚洲av日韩精品久久久久久密 | 亚洲精品国产一区二区精华液| 一级毛片我不卡| 乱人伦中国视频| 久久久精品区二区三区| 两个人免费观看高清视频| 高清黄色对白视频在线免费看| 狂野欧美激情性xxxx| 无限看片的www在线观看| 97在线人人人人妻| 欧美激情高清一区二区三区| 18禁裸乳无遮挡动漫免费视频| 一本综合久久免费| 免费少妇av软件| 国产女主播在线喷水免费视频网站| 美女扒开内裤让男人捅视频| 亚洲国产精品成人久久小说| 两个人免费观看高清视频| 久久天躁狠狠躁夜夜2o2o | 久久精品aⅴ一区二区三区四区| 亚洲黑人精品在线| 男女边吃奶边做爰视频| 秋霞在线观看毛片| 午夜福利,免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 高潮久久久久久久久久久不卡| 免费在线观看影片大全网站 | 亚洲国产av影院在线观看| 99九九在线精品视频| 亚洲成人免费av在线播放| 精品免费久久久久久久清纯 | av电影中文网址| 无遮挡黄片免费观看| 国产福利在线免费观看视频| 久久这里只有精品19| 精品久久久精品久久久| 欧美精品高潮呻吟av久久| 国产精品久久久久久人妻精品电影 | av不卡在线播放| 日韩精品免费视频一区二区三区| 欧美精品一区二区大全| 免费av中文字幕在线| 97精品久久久久久久久久精品| 欧美黑人欧美精品刺激| 波多野结衣一区麻豆| 丝袜美足系列| 午夜影院在线不卡| 1024视频免费在线观看| 久久综合国产亚洲精品| 中文字幕制服av| 最近中文字幕2019免费版| 熟女少妇亚洲综合色aaa.| 妹子高潮喷水视频| 亚洲成色77777| 99久久精品国产亚洲精品| 亚洲国产最新在线播放| 女人高潮潮喷娇喘18禁视频| 狠狠婷婷综合久久久久久88av| 777米奇影视久久| 国产精品免费视频内射| 人妻一区二区av| 国产成人av教育| 高清欧美精品videossex| 国产真人三级小视频在线观看| 中国美女看黄片| av在线播放精品| videos熟女内射| 精品福利永久在线观看| 黄片播放在线免费| 新久久久久国产一级毛片| 一级片免费观看大全| 亚洲精品久久午夜乱码| 极品少妇高潮喷水抽搐| 亚洲国产欧美一区二区综合| 国产精品一区二区免费欧美 | 久久精品国产a三级三级三级| 丝袜美腿诱惑在线| 黑人猛操日本美女一级片| 亚洲欧美激情在线| 久久女婷五月综合色啪小说| 成年女人毛片免费观看观看9 | 精品欧美一区二区三区在线| 久久久久国产精品人妻一区二区| 亚洲精品美女久久av网站| 亚洲免费av在线视频| 99热全是精品| 久久女婷五月综合色啪小说| 男女边吃奶边做爰视频| 美女视频免费永久观看网站| 欧美日韩亚洲综合一区二区三区_| 国产精品欧美亚洲77777| 另类精品久久| 精品少妇黑人巨大在线播放| 亚洲欧洲国产日韩| 成年人午夜在线观看视频| 另类亚洲欧美激情| 午夜老司机福利片| 97精品久久久久久久久久精品| 国产在线视频一区二区| 亚洲成国产人片在线观看| 精品福利观看| 一本一本久久a久久精品综合妖精| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧美精品综合一区二区三区| 国产成人精品在线电影| 熟女少妇亚洲综合色aaa.| 国产高清视频在线播放一区 | 99国产精品一区二区蜜桃av | 国产精品久久久久久精品电影小说| 18禁国产床啪视频网站| 中文字幕人妻熟女乱码| 女性被躁到高潮视频| 亚洲成av片中文字幕在线观看| 老司机午夜十八禁免费视频| 美女国产高潮福利片在线看| 欧美黄色片欧美黄色片| 少妇人妻 视频| 国产精品久久久久久精品电影小说| 91精品三级在线观看| 国产一区亚洲一区在线观看| 妹子高潮喷水视频| 中文字幕另类日韩欧美亚洲嫩草| 精品一区在线观看国产| 黑人巨大精品欧美一区二区蜜桃| 国产精品久久久久久人妻精品电影 | 激情视频va一区二区三区| 欧美日韩亚洲综合一区二区三区_| 黑人欧美特级aaaaaa片| 欧美av亚洲av综合av国产av| 一区在线观看完整版| 丰满少妇做爰视频| 国产精品久久久av美女十八| 成人黄色视频免费在线看| 亚洲一区中文字幕在线| 亚洲中文av在线| 激情五月婷婷亚洲| 涩涩av久久男人的天堂| 免费观看a级毛片全部| 国产精品熟女久久久久浪| 亚洲国产精品999| 精品久久久久久电影网| 亚洲精品成人av观看孕妇| 国产高清不卡午夜福利| 国产视频首页在线观看| 亚洲欧美色中文字幕在线| 免费在线观看黄色视频的| 777久久人妻少妇嫩草av网站| 啦啦啦 在线观看视频| 亚洲精品国产一区二区精华液| 国产高清视频在线播放一区 | 久久鲁丝午夜福利片| 国产精品三级大全| 亚洲精品国产区一区二| av网站免费在线观看视频| 成人18禁高潮啪啪吃奶动态图| 涩涩av久久男人的天堂| 日本色播在线视频| 在现免费观看毛片| 久久精品人人爽人人爽视色| 美女扒开内裤让男人捅视频| 91字幕亚洲| 日韩一本色道免费dvd| 亚洲欧美一区二区三区久久| 午夜免费男女啪啪视频观看| 女警被强在线播放| 亚洲欧美精品综合一区二区三区| 视频区欧美日本亚洲| 一级毛片电影观看| 自线自在国产av| 免费看十八禁软件| 亚洲国产最新在线播放| 国产欧美日韩精品亚洲av| √禁漫天堂资源中文www| 涩涩av久久男人的天堂| 欧美人与性动交α欧美精品济南到| 欧美国产精品va在线观看不卡| 国产在线免费精品| 成人18禁高潮啪啪吃奶动态图| 在线观看一区二区三区激情| 少妇人妻 视频| 亚洲精品国产区一区二| 这个男人来自地球电影免费观看| 又黄又粗又硬又大视频| 亚洲av国产av综合av卡| 首页视频小说图片口味搜索 | 最新的欧美精品一区二区| 99九九在线精品视频| 精品少妇黑人巨大在线播放| 男女边吃奶边做爰视频| 久久久国产欧美日韩av| 人人妻人人澡人人看| 婷婷色av中文字幕| 嫁个100分男人电影在线观看 | 国产成人精品久久久久久| netflix在线观看网站| 亚洲国产日韩一区二区| 久久久国产一区二区| 中文字幕制服av| 欧美日韩成人在线一区二区| a级毛片黄视频| 天天添夜夜摸| 日韩av免费高清视频| 天天躁狠狠躁夜夜躁狠狠躁| 日韩伦理黄色片| 99热全是精品| 亚洲av欧美aⅴ国产| 热99国产精品久久久久久7| 国产成人精品在线电影| 午夜免费鲁丝| 亚洲成人国产一区在线观看 | 久久国产精品男人的天堂亚洲| 精品国产一区二区久久| 亚洲国产欧美日韩在线播放| 婷婷色麻豆天堂久久| 午夜激情av网站| 日本a在线网址| 国产成人精品久久二区二区免费| 亚洲久久久国产精品| 成人影院久久| 午夜福利一区二区在线看| 日韩中文字幕视频在线看片| videos熟女内射| 免费不卡黄色视频| 国产在线一区二区三区精| 嫁个100分男人电影在线观看 | 亚洲精品自拍成人| 午夜激情av网站| 日韩人妻精品一区2区三区| 手机成人av网站| 51午夜福利影视在线观看| 一本久久精品| 老司机影院毛片| 精品少妇内射三级| 欧美激情极品国产一区二区三区| xxxhd国产人妻xxx| 亚洲,欧美,日韩| 青春草亚洲视频在线观看| 亚洲视频免费观看视频| www.熟女人妻精品国产| 日韩 欧美 亚洲 中文字幕| 国产成人一区二区三区免费视频网站 | 国产日韩欧美视频二区| 久久天躁狠狠躁夜夜2o2o | 国产女主播在线喷水免费视频网站| 无限看片的www在线观看| 国产精品av久久久久免费| 99九九在线精品视频| 国产亚洲欧美在线一区二区| 国精品久久久久久国模美| 国产一区二区激情短视频 | 大片免费播放器 马上看| 久久久国产精品麻豆| 色播在线永久视频| 90打野战视频偷拍视频| 黄色毛片三级朝国网站| 久久鲁丝午夜福利片| 日本av手机在线免费观看| 久久av网站| 青青草视频在线视频观看| 亚洲男人天堂网一区| 黄色a级毛片大全视频| 高清视频免费观看一区二区| 久久人人爽av亚洲精品天堂| 人人妻人人爽人人添夜夜欢视频| 久久人人爽av亚洲精品天堂| 亚洲av美国av| 精品第一国产精品| 久久av网站| 青青草视频在线视频观看| 亚洲av美国av| 50天的宝宝边吃奶边哭怎么回事| 欧美激情 高清一区二区三区| 久久影院123| 欧美黄色片欧美黄色片| 人人澡人人妻人| 捣出白浆h1v1| 啦啦啦视频在线资源免费观看| 亚洲中文av在线| 国产免费一区二区三区四区乱码| 午夜福利视频在线观看免费| av片东京热男人的天堂| videos熟女内射| 日韩大片免费观看网站| 好男人电影高清在线观看| 又大又黄又爽视频免费| 成人国产一区最新在线观看 | 欧美 日韩 精品 国产| 最黄视频免费看| 日韩制服丝袜自拍偷拍| 免费不卡黄色视频| tube8黄色片| 午夜激情久久久久久久| 精品人妻在线不人妻| 成人三级做爰电影| 天天影视国产精品| 婷婷丁香在线五月| 国产精品久久久av美女十八| 在线av久久热| 国产日韩欧美在线精品| 高清av免费在线| 国产午夜精品一二区理论片| 在线观看国产h片| 老汉色av国产亚洲站长工具| 国产1区2区3区精品| 黄色一级大片看看| 国产成人免费无遮挡视频| 天天添夜夜摸| 老司机深夜福利视频在线观看 | 久久精品成人免费网站| 激情视频va一区二区三区| 99香蕉大伊视频| 韩国精品一区二区三区| 女人精品久久久久毛片| 日韩视频在线欧美| 欧美精品av麻豆av| 老汉色av国产亚洲站长工具| 亚洲成人手机| 国产免费福利视频在线观看| 国产精品一区二区精品视频观看| 亚洲五月婷婷丁香| 国产欧美日韩一区二区三区在线| 欧美精品一区二区免费开放| 九色亚洲精品在线播放| av网站免费在线观看视频| 国产精品久久久久久精品电影小说| 国产精品久久久久久人妻精品电影 | 韩国高清视频一区二区三区| 国产精品一国产av| 狂野欧美激情性bbbbbb| 天堂中文最新版在线下载| 一边摸一边抽搐一进一出视频| 视频在线观看一区二区三区| 国产无遮挡羞羞视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 老熟女久久久| 久久久久国产一级毛片高清牌| 国产av一区二区精品久久| 成人18禁高潮啪啪吃奶动态图| 亚洲精品美女久久久久99蜜臀 | 国产精品.久久久| 最近手机中文字幕大全| 人成视频在线观看免费观看| 国产精品国产三级国产专区5o| 中文字幕人妻熟女乱码| 国产野战对白在线观看| 一边亲一边摸免费视频| 黄片小视频在线播放| 国产1区2区3区精品| 国产亚洲精品久久久久5区| 欧美精品亚洲一区二区| 一区二区三区精品91| 999久久久国产精品视频| 亚洲成国产人片在线观看| 在线观看免费视频网站a站| 亚洲情色 制服丝袜| 国产欧美亚洲国产| 午夜福利影视在线免费观看| 欧美精品人与动牲交sv欧美| 超碰97精品在线观看| 亚洲熟女精品中文字幕|