• <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)面擬合
    天天影视国产精品| 日本-黄色视频高清免费观看| 亚洲精品日韩在线中文字幕| 多毛熟女@视频| av电影中文网址| 亚洲国产成人一精品久久久| 伊人久久精品亚洲午夜| av在线观看视频网站免费| 亚洲精品色激情综合| 午夜福利视频精品| 久久久久久久久久久丰满| 国产淫语在线视频| 男的添女的下面高潮视频| 18禁裸乳无遮挡动漫免费视频| 精品少妇久久久久久888优播| 少妇熟女欧美另类| 日韩精品免费视频一区二区三区 | 男女免费视频国产| 18禁裸乳无遮挡动漫免费视频| 丝袜美足系列| 内地一区二区视频在线| 欧美最新免费一区二区三区| 99热网站在线观看| 99久久综合免费| 另类精品久久| 亚洲丝袜综合中文字幕| 欧美成人精品欧美一级黄| freevideosex欧美| 日韩,欧美,国产一区二区三区| 亚洲不卡免费看| 人妻系列 视频| 一本一本综合久久| 2021少妇久久久久久久久久久| a级毛片黄视频| 国产精品国产三级国产专区5o| 亚洲av综合色区一区| 精品亚洲乱码少妇综合久久| 免费少妇av软件| 国产精品国产av在线观看| 精品人妻熟女毛片av久久网站| 你懂的网址亚洲精品在线观看| 伊人亚洲综合成人网| 欧美精品一区二区免费开放| 在线观看三级黄色| 亚洲精品国产av成人精品| 高清在线视频一区二区三区| 国产精品.久久久| 精品国产乱码久久久久久小说| 黄色怎么调成土黄色| 韩国高清视频一区二区三区| 亚洲性久久影院| 最新的欧美精品一区二区| 有码 亚洲区| 国产亚洲精品久久久com| 91精品国产九色| 亚洲性久久影院| av网站免费在线观看视频| 久久人妻熟女aⅴ| 黄色视频在线播放观看不卡| 插阴视频在线观看视频| 国产亚洲欧美精品永久| 啦啦啦在线观看免费高清www| 久久久亚洲精品成人影院| 如日韩欧美国产精品一区二区三区 | 欧美最新免费一区二区三区| 日韩一本色道免费dvd| 久久久久精品久久久久真实原创| 一级片'在线观看视频| 亚洲,一卡二卡三卡| 免费黄色在线免费观看| 国产伦理片在线播放av一区| 欧美精品一区二区免费开放| 欧美性感艳星| 嘟嘟电影网在线观看| 精品少妇黑人巨大在线播放| 国产极品粉嫩免费观看在线 | 国产精品麻豆人妻色哟哟久久| videosex国产| 国产精品秋霞免费鲁丝片| 少妇猛男粗大的猛烈进出视频| 99热这里只有精品一区| 久久av网站| a级毛片在线看网站| 中文字幕最新亚洲高清| 中文字幕免费在线视频6| 性高湖久久久久久久久免费观看| 国产伦精品一区二区三区视频9| 精品久久久精品久久久| 热99久久久久精品小说推荐| 最黄视频免费看| av卡一久久| 亚洲av在线观看美女高潮| 亚洲成人av在线免费| 成人国产av品久久久| 在线看a的网站| 国产国拍精品亚洲av在线观看| 人妻少妇偷人精品九色| 亚洲无线观看免费| 国产精品秋霞免费鲁丝片| 午夜老司机福利剧场| 亚洲国产毛片av蜜桃av| 最新的欧美精品一区二区| 国产成人av激情在线播放 | 亚洲精品,欧美精品| 久久人人爽av亚洲精品天堂| 久久精品国产鲁丝片午夜精品| 99国产综合亚洲精品| 成人国语在线视频| 一边摸一边做爽爽视频免费| 亚洲av二区三区四区| 99热国产这里只有精品6| 热re99久久国产66热| 97在线人人人人妻| 国产永久视频网站| 免费观看性生交大片5| 女性生殖器流出的白浆| 精品久久蜜臀av无| 午夜福利影视在线免费观看| 人妻 亚洲 视频| 久久免费观看电影| 精品少妇内射三级| 一区二区日韩欧美中文字幕 | 国产伦精品一区二区三区视频9| 日韩免费高清中文字幕av| 国产免费又黄又爽又色| freevideosex欧美| 一区二区三区乱码不卡18| 国产高清三级在线| 大香蕉久久成人网| 国产高清有码在线观看视频| tube8黄色片| 天天躁夜夜躁狠狠久久av| 久久影院123| 如何舔出高潮| 熟女电影av网| 综合色丁香网| 国产成人精品一,二区| 久久久久久久久久久久大奶| 国产成人a∨麻豆精品| 一级爰片在线观看| 国产精品无大码| 在线观看免费高清a一片| 精品亚洲成国产av| 街头女战士在线观看网站| 久久青草综合色| 久久精品久久久久久久性| 国产精品久久久久久精品古装| 日本vs欧美在线观看视频| 十分钟在线观看高清视频www| 精品久久久久久久久av| 国产免费视频播放在线视频| 最近的中文字幕免费完整| av黄色大香蕉| 在线精品无人区一区二区三| 寂寞人妻少妇视频99o| 国产极品粉嫩免费观看在线 | 亚洲精品国产av成人精品| 18禁动态无遮挡网站| 久久毛片免费看一区二区三区| 永久网站在线| 极品人妻少妇av视频| 久久97久久精品| 极品少妇高潮喷水抽搐| 视频在线观看一区二区三区| 91精品伊人久久大香线蕉| 久久久久久久亚洲中文字幕| 汤姆久久久久久久影院中文字幕| 午夜老司机福利剧场| av视频免费观看在线观看| 欧美日韩视频精品一区| 久久久久网色| 国产成人freesex在线| 十分钟在线观看高清视频www| 久久精品夜色国产| 久久久久久久久久人人人人人人| 男女国产视频网站| 日日啪夜夜爽| 欧美日韩成人在线一区二区| 热99国产精品久久久久久7| 国产不卡av网站在线观看| 少妇被粗大猛烈的视频| 日韩亚洲欧美综合| 国产午夜精品一二区理论片| 亚洲av日韩在线播放| 亚洲精品中文字幕在线视频| 在线免费观看不下载黄p国产| 美女视频免费永久观看网站| 日韩电影二区| 国产成人精品在线电影| 国产午夜精品久久久久久一区二区三区| 一级毛片电影观看| 美女国产高潮福利片在线看| 少妇精品久久久久久久| 国产在线免费精品| av播播在线观看一区| 中国国产av一级| 黄色视频在线播放观看不卡| 成年女人在线观看亚洲视频| 欧美三级亚洲精品| 91在线精品国自产拍蜜月| 又黄又爽又刺激的免费视频.| 国产黄色免费在线视频| 日韩免费高清中文字幕av| a级片在线免费高清观看视频| 精品国产一区二区久久| 麻豆精品久久久久久蜜桃| 亚洲成色77777| 一级a做视频免费观看| 国产精品蜜桃在线观看| 亚洲精品乱码久久久久久按摩| 黄色毛片三级朝国网站| 亚洲精华国产精华液的使用体验| 亚洲精品色激情综合| 欧美亚洲日本最大视频资源| 22中文网久久字幕| 美女内射精品一级片tv| 国产精品99久久久久久久久| av又黄又爽大尺度在线免费看| 国产亚洲一区二区精品| 成年美女黄网站色视频大全免费 | 考比视频在线观看| av又黄又爽大尺度在线免费看| 国产一级毛片在线| a级毛片在线看网站| 美女福利国产在线| 久久久欧美国产精品| 亚洲av.av天堂| 国产亚洲一区二区精品| 久久久亚洲精品成人影院| 18禁在线无遮挡免费观看视频| 午夜视频国产福利| 自拍欧美九色日韩亚洲蝌蚪91| 婷婷成人精品国产| 亚洲激情五月婷婷啪啪| 国产精品久久久久成人av| 久久精品国产亚洲av涩爱| 高清av免费在线| 欧美成人午夜免费资源| 哪个播放器可以免费观看大片| 老司机亚洲免费影院| 看免费成人av毛片| 日韩av免费高清视频| 纯流量卡能插随身wifi吗| 91久久精品电影网| 婷婷色综合www| 国产精品99久久久久久久久| 色哟哟·www| 99久久人妻综合| 2022亚洲国产成人精品| 中文字幕精品免费在线观看视频 | 成人黄色视频免费在线看| av女优亚洲男人天堂| 成人免费观看视频高清| 天堂8中文在线网| 日本欧美视频一区| 中国三级夫妇交换| 亚洲国产精品国产精品| 美女主播在线视频| 十八禁网站网址无遮挡| 搡老乐熟女国产| 91成人精品电影| 欧美亚洲 丝袜 人妻 在线| 国产精品久久久久成人av| 亚洲精品美女久久av网站| 久久久精品94久久精品| 午夜影院在线不卡| 久久久久久伊人网av| 国产69精品久久久久777片| 成年美女黄网站色视频大全免费 | 大话2 男鬼变身卡| 国产精品久久久久久av不卡| 欧美激情极品国产一区二区三区 | 麻豆乱淫一区二区| 亚洲激情五月婷婷啪啪| 水蜜桃什么品种好| 中文天堂在线官网| 免费观看的影片在线观看| 亚洲国产日韩一区二区| 日本av手机在线免费观看| 精品视频人人做人人爽| 国产国语露脸激情在线看| 99热全是精品| 精品人妻偷拍中文字幕| 人妻少妇偷人精品九色| 你懂的网址亚洲精品在线观看| 一二三四中文在线观看免费高清| 熟女av电影| av又黄又爽大尺度在线免费看| 亚洲图色成人| 一区二区av电影网| 蜜桃久久精品国产亚洲av| 在线观看免费日韩欧美大片 | 亚洲av国产av综合av卡| 黄片播放在线免费| 国内精品宾馆在线| av一本久久久久| 一本大道久久a久久精品| 国产深夜福利视频在线观看| 亚洲成人av在线免费| 国产国拍精品亚洲av在线观看| 天堂俺去俺来也www色官网| 91精品三级在线观看| 日本wwww免费看| 国产老妇伦熟女老妇高清| kizo精华| 搡女人真爽免费视频火全软件| 亚洲av免费高清在线观看| 成人午夜精彩视频在线观看| 中文欧美无线码| 嫩草影院入口| 日韩一区二区视频免费看| 免费人妻精品一区二区三区视频| 啦啦啦视频在线资源免费观看| tube8黄色片| 成人手机av| 丝袜美足系列| 久久久久国产网址| 9色porny在线观看| 午夜免费男女啪啪视频观看| 国国产精品蜜臀av免费| av专区在线播放| 亚洲一级一片aⅴ在线观看| 乱人伦中国视频| 最黄视频免费看| 亚洲精品中文字幕在线视频| 亚洲熟女精品中文字幕| 亚洲成人一二三区av| 只有这里有精品99| 性色av一级| 一级,二级,三级黄色视频| 我要看黄色一级片免费的| 亚洲精品一二三| 成年女人在线观看亚洲视频| 国产精品国产av在线观看| 亚洲精品视频女| 国产精品女同一区二区软件| 乱码一卡2卡4卡精品| 亚洲精品aⅴ在线观看| 亚洲欧洲日产国产| 老司机亚洲免费影院| 亚洲欧美色中文字幕在线| 三级国产精品欧美在线观看| 夫妻午夜视频| 国产乱来视频区| 人妻制服诱惑在线中文字幕| 精品久久久久久久久亚洲| 亚洲精品国产av成人精品| 国产精品不卡视频一区二区| 99热国产这里只有精品6| 中文乱码字字幕精品一区二区三区| 欧美+日韩+精品| 七月丁香在线播放| 韩国高清视频一区二区三区| 精品国产露脸久久av麻豆| 99九九线精品视频在线观看视频| 极品少妇高潮喷水抽搐| 18在线观看网站| 亚洲高清免费不卡视频| 赤兔流量卡办理| 国产成人免费观看mmmm| 丰满饥渴人妻一区二区三| 插逼视频在线观看| 国产亚洲最大av| 十八禁高潮呻吟视频| 男女啪啪激烈高潮av片| 成人18禁高潮啪啪吃奶动态图 | 日韩制服骚丝袜av| 热99国产精品久久久久久7| www.av在线官网国产| 国产高清国产精品国产三级| 晚上一个人看的免费电影| av在线观看视频网站免费| 国产精品一区www在线观看| 国产成人freesex在线| 久久人人爽人人爽人人片va| a级毛色黄片| 国产女主播在线喷水免费视频网站| 国产精品久久久久久av不卡| 日本黄色日本黄色录像| 一边亲一边摸免费视频| 狠狠精品人妻久久久久久综合| 国产永久视频网站| 美女福利国产在线| 麻豆成人av视频| 精品午夜福利在线看| 秋霞在线观看毛片| 一级毛片黄色毛片免费观看视频| 成人漫画全彩无遮挡| 亚洲精品一区蜜桃| 少妇人妻 视频| 内地一区二区视频在线| 色吧在线观看| 夜夜爽夜夜爽视频| 国产精品三级大全| 亚洲欧美成人精品一区二区| 国产综合精华液| 99re6热这里在线精品视频| 最近中文字幕高清免费大全6| 一个人看视频在线观看www免费| 校园人妻丝袜中文字幕| 人人妻人人澡人人爽人人夜夜| 免费观看无遮挡的男女| 国产av码专区亚洲av| 亚洲精品国产av成人精品| 日本黄大片高清| 五月天丁香电影| 国产永久视频网站| 另类精品久久| 亚洲国产av影院在线观看| 乱码一卡2卡4卡精品| 亚洲av福利一区| 久久这里有精品视频免费| 国产高清国产精品国产三级| 久久久久国产网址| 插阴视频在线观看视频| 美女大奶头黄色视频| 51国产日韩欧美| 最新的欧美精品一区二区| 热re99久久精品国产66热6| 制服诱惑二区| 伊人久久国产一区二区| 国产亚洲欧美精品永久| 日日摸夜夜添夜夜添av毛片| 免费高清在线观看视频在线观看| 亚洲国产欧美日韩在线播放| 国产极品粉嫩免费观看在线 | 一区在线观看完整版| 18禁观看日本| 日韩中文字幕视频在线看片| 热99久久久久精品小说推荐| 久久青草综合色| 欧美最新免费一区二区三区| 国产成人a∨麻豆精品| 91精品三级在线观看| 国产精品国产av在线观看| 人人妻人人澡人人看| 嘟嘟电影网在线观看| 超色免费av| 午夜老司机福利剧场| 久久这里有精品视频免费| 这个男人来自地球电影免费观看 | 亚洲四区av| 欧美+日韩+精品| 久久久久人妻精品一区果冻| 少妇精品久久久久久久| 国产精品秋霞免费鲁丝片| 国产精品久久久久久精品古装| 午夜日本视频在线| 日韩亚洲欧美综合| 欧美最新免费一区二区三区| 亚洲色图综合在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 丝袜脚勾引网站| 插逼视频在线观看| 高清不卡的av网站| 高清视频免费观看一区二区| 亚洲精品色激情综合| 日日啪夜夜爽| 乱码一卡2卡4卡精品| 高清视频免费观看一区二区| 各种免费的搞黄视频| 久久精品国产鲁丝片午夜精品| 亚洲国产欧美在线一区| 人妻一区二区av| 黄色毛片三级朝国网站| 国产免费福利视频在线观看| 久久99热6这里只有精品| 又大又黄又爽视频免费| 日本91视频免费播放| 多毛熟女@视频| 精品少妇黑人巨大在线播放| 国产在视频线精品| 最新的欧美精品一区二区| 久热久热在线精品观看| 久久久久久伊人网av| 亚洲av成人精品一区久久| 2018国产大陆天天弄谢| 中文字幕久久专区| 又大又黄又爽视频免费| 777米奇影视久久| 精品99又大又爽又粗少妇毛片| 久久久亚洲精品成人影院| 18禁在线无遮挡免费观看视频| 亚洲欧美成人精品一区二区| 又粗又硬又长又爽又黄的视频| 晚上一个人看的免费电影| 日韩av不卡免费在线播放| 一区二区三区免费毛片| 中国美白少妇内射xxxbb| 成人国产av品久久久| 99国产综合亚洲精品| 人妻少妇偷人精品九色| 少妇被粗大猛烈的视频| 男人添女人高潮全过程视频| 久久综合国产亚洲精品| 欧美日本中文国产一区发布| 免费日韩欧美在线观看| 亚洲精品日韩av片在线观看| 国产成人精品一,二区| 9色porny在线观看| 熟女电影av网| 国产精品一二三区在线看| 亚洲欧美成人综合另类久久久| 丁香六月天网| xxx大片免费视频| 午夜福利视频在线观看免费| av一本久久久久| 久久鲁丝午夜福利片| 精品少妇久久久久久888优播| 不卡视频在线观看欧美| 夜夜骑夜夜射夜夜干| av在线观看视频网站免费| 亚洲国产精品国产精品| 亚洲欧洲精品一区二区精品久久久 | 超色免费av| 亚洲av在线观看美女高潮| 多毛熟女@视频| 性高湖久久久久久久久免费观看| 人妻少妇偷人精品九色| 免费播放大片免费观看视频在线观看| 男女无遮挡免费网站观看| 国产综合精华液| 国产精品久久久久久久电影| 人妻系列 视频| 久久久久久久久大av| 午夜视频国产福利| 大片免费播放器 马上看| 日韩一本色道免费dvd| 亚洲美女搞黄在线观看| av一本久久久久| 波野结衣二区三区在线| 韩国av在线不卡| 久久久久久伊人网av| 欧美老熟妇乱子伦牲交| 夜夜骑夜夜射夜夜干| 亚洲美女黄色视频免费看| 欧美精品高潮呻吟av久久| 国产成人av激情在线播放 | 精品少妇内射三级| 亚洲精品亚洲一区二区| 久久精品国产亚洲av天美| 一级片'在线观看视频| 免费大片黄手机在线观看| 亚洲欧美中文字幕日韩二区| 人人妻人人澡人人看| 午夜视频国产福利| av国产久精品久网站免费入址| 国产极品天堂在线| 精品人妻一区二区三区麻豆| 9色porny在线观看| 亚洲欧美一区二区三区黑人 | 少妇 在线观看| 亚洲怡红院男人天堂| kizo精华| 99久久精品国产国产毛片| 亚洲熟女精品中文字幕| 成年女人在线观看亚洲视频| 毛片一级片免费看久久久久| 男女边摸边吃奶| 97在线人人人人妻| 黄片播放在线免费| 久久99一区二区三区| 欧美 日韩 精品 国产| 精品久久国产蜜桃| 免费高清在线观看日韩| 久久99热6这里只有精品| 亚洲精品色激情综合| 久久久久国产网址| 狂野欧美白嫩少妇大欣赏| 久久久久久久久大av| 欧美人与善性xxx| 久久久久网色| 精品国产一区二区久久| 日韩成人伦理影院| 美女cb高潮喷水在线观看| 欧美人与性动交α欧美精品济南到 | 久久人人爽av亚洲精品天堂| 在线看a的网站| 蜜桃久久精品国产亚洲av| 久久久久人妻精品一区果冻| 国产精品一区二区三区四区免费观看| 18禁在线无遮挡免费观看视频| 欧美三级亚洲精品| 免费观看的影片在线观看| 亚洲av日韩在线播放| 亚洲欧美一区二区三区国产| 国产永久视频网站| 国产日韩欧美视频二区| 免费观看a级毛片全部| 国产精品久久久久久久久免| 欧美日韩视频精品一区| 午夜精品国产一区二区电影| 婷婷色综合大香蕉| 久久精品国产亚洲av天美| 精品午夜福利在线看| www.av在线官网国产| 少妇熟女欧美另类| 日韩中文字幕视频在线看片| 女性生殖器流出的白浆| 建设人人有责人人尽责人人享有的| 亚洲综合色网址| 波野结衣二区三区在线| 超碰97精品在线观看| 夫妻午夜视频| 你懂的网址亚洲精品在线观看| 女的被弄到高潮叫床怎么办| 亚洲欧洲精品一区二区精品久久久 | 精品国产一区二区三区久久久樱花| 国产日韩欧美视频二区| 国产男人的电影天堂91| 亚洲成色77777| 另类精品久久| 97在线视频观看| 欧美三级亚洲精品| 蜜桃久久精品国产亚洲av|