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

    適用于Ku波段雷達高度計海冰干舷高度反演的積雪校正方法

    2020-06-12 04:37:38趙全芳孟俊敏劉眉潔
    海洋技術(shù)學(xué)報 2020年2期
    關(guān)鍵詞:干舷高度計海冰

    張 晰,趙全芳,孟俊敏,劉眉潔

    (1.山東科技大學(xué) 測繪科學(xué)與工程學(xué)院,山東 青島 266590;2.自然資源部第一海洋研究所,山東 青島 266061;3.青島大學(xué) 物理學(xué)院,山東 青島 266071)

    海冰直接影響著全球的氣候變化,并通過與大氣的相互作用反饋于全球環(huán)境系統(tǒng)。海冰厚度是重要的氣候環(huán)境變化表現(xiàn)因子。準確估算海冰厚度,對于氣候變化研究、極區(qū)航行保障具有重要意義。利用遙感手段準確估算大范圍海冰厚度始終是當前的國際難題。相比于人工實測、電磁感應(yīng)和仰視聲吶等手段,近年來發(fā)展的雷達高度計測高技術(shù)使得獲取連續(xù)大范圍的海冰厚度成為可能[1]。到目前為止,ERS-1/2、ENVISAT、HY-2、AltiKa、CryoSat-2和Sentinel-3A等雷達高度計相繼發(fā)射,為海冰厚度估算提供了多源的衛(wèi)星觀測手段。在眾多雷達高度計中,工作在Ku波段的CryoSat-2(CS-2)和Sentinel-3A(S3)合成孔徑雷達高度計相較于傳統(tǒng)雷達高度計,具有更高的空間分辨率(沿軌約0.3 km,CS-2/S3交軌約1.5/1.64 km),能夠提供較高分辨率的海冰厚度信息,是當前最先進的衛(wèi)星雷達高度計。

    雷達高度計并不直接計算海冰厚度,是通過計算海冰的冰面出水高度(海冰干舷高度),并結(jié)合浮體定律實現(xiàn)海冰厚度估算[2-3]。因此,準確反演海冰干舷高度對于估算海冰厚度至關(guān)重要。由于衛(wèi)星高度計信號回波的實際時間跟蹤點與預(yù)設(shè)時間跟蹤點存在位置上的偏差,需要對回波波形進行重跟蹤處理,以計算實際時間跟蹤點,進而校正得到衛(wèi)星高度計質(zhì)心到地面點的真實距離。Ricker等[4]使用波形重跟蹤方法指出重跟蹤閾值的選擇對海冰干舷和厚度的估計值有顯著影響。在極區(qū),海冰表面普遍覆蓋著積雪,當雷達高度計探測海冰表面高度時,Ku波段電磁波須穿過覆蓋于海冰之上的積雪層,因此在計算海冰干舷高度時,需考慮電磁波在積雪中穿透速度較慢的距離校正問題。在實際情況下,積雪層受海冰析鹽過程的影響,常會在積雪的上下層形成強鹽度梯度變化[5-7];另外,積雪層內(nèi)粒徑也會與電磁波發(fā)生體散射[8]。這些因素使得微波雷達高度計的電磁波很難完全穿透積雪層直達海冰表面。對于短波長的Ku波段雷達,該現(xiàn)象更為明顯,這將導(dǎo)致雷達信號主散射面向上偏離雪-冰交界面,從而影響海冰干舷高度的估計。但大多數(shù)研究者在使用雷達高度計反演海冰干舷高度時,常假設(shè)Ku波段電磁波可以完全穿透積雪層,直接將雪-冰交界面當作雷達信號的主散射面[2-4,9-12],利用海冰上的整體積雪深度進行距離校正。Armitage和Ridout[13]分析了CS-2雷達高度計電磁波在積雪中的雷達穿透系數(shù),發(fā)現(xiàn)Ku波段電磁波在垂直入射時無法完全穿透積雪層,并指出若不考慮積雪層對電磁波的影響,海冰干舷高度將被高估,從而大大影響海冰厚度的估算。Ricker等[14]發(fā)現(xiàn),未知雷達信號穿透到積雪層中的不確定性可能導(dǎo)致大約0.06~0.12 m的海冰干舷高度偏差。

    因此,選擇合適的冰間水道與海冰的重跟蹤閾值以及確定雷達信號在積雪中的穿透深度,對精確地反演海冰干舷高度起決定性作用?;诖耍疚姆謩e針對CS-2與S3兩種高度計衛(wèi)星開展了確定最優(yōu)波形重跟蹤閾值的分析,并詳細分析了積雪層對Ku波段CS-2和S3合成孔徑雷達高度計主散射面位置的影響,計算了Ku波段電磁波的積雪穿透系數(shù),針對海冰干舷高度精確反演的需求,改進了積雪校正算法。

    1 數(shù)據(jù)

    1.1 CryoSat-2雷達高度計數(shù)據(jù)

    CryoSat-2是歐洲空間局發(fā)射的合成孔徑雷達高度計,其工作在Ku波段,中心頻率為13.575 GHz,帶寬約為320 MHz,空間覆蓋范圍達到南北緯88°,運行重復(fù)子周期為30 d,其沿軌道分辨率可達0.3 km,交軌道分辨率優(yōu)于1.5 km,相比于足跡為10 km的傳統(tǒng)雷達高度計[15],其測量精度顯著提升。該衛(wèi)星有LRM(低分辨率)、SAR和SARin 3種工作模式,本文采用的是SAR模式L1b級數(shù)據(jù)。在該模式下,CS-2的每一個波形采樣窗口大小為60 m,包含有256個距離門[16]。為了進行波形重跟蹤校正,得到海冰/水的表面高程,還需對L1b數(shù)據(jù)進行波形噪聲去除和地球物理校正(地球物理校正包含:干/濕對流層、電離層、逆氣壓、海洋潮汐、長周期平衡潮、海洋負荷潮、固體潮和地心極潮校正等處理)。本文使用的是北極區(qū)域、時間為2017年3月-4月和2018年4月的雷達數(shù)據(jù)。

    1.2 Sentinel-3A雷達高度計數(shù)據(jù)

    Sentinel-3A也是歐空局發(fā)射的合成孔徑雷達高度計,它是第一個100%以SAR模式覆蓋海洋、冰區(qū)和內(nèi)陸水域的雷達高度計。其工作的Ku波段中心頻率與CS-2相同,帶寬約為350 MHz,空間覆蓋范圍為南北緯81.35°,運行重復(fù)子周期為4 d,沿軌道分辨率約為0.3 km,交軌道分辨率大約為1.64 km[17]。本文采用的是SAR模式的L2 WAT Enhanced(Water Enhanced)數(shù)據(jù)。在該模式下波形采樣窗口大小為60 m,包含有128個距離門。與CS-2相同,使用的也是北極區(qū)域2017年3月-4月和2018年4月的數(shù)據(jù),并對這些數(shù)據(jù)進行了與CS-2相同的地球物理校正處理。

    1.3 Operation IceBridge數(shù)據(jù)

    本文使用美國國家航空航天局(NASA)發(fā)布的OIB機載數(shù)據(jù)作為海冰干舷高度的驗證數(shù)據(jù),其每年均會在北極、格陵蘭島和南極上空開展機載探測飛行。該數(shù)據(jù)由機載地形掃描測繪儀(Airborne Topographic Mapper,ATM)、數(shù)字測繪相機(Digital Mapping System Camera,DMS) 和雪雷達(Snow Radar)3種裝備采集得到。ATM用于測量表面高度,空間分辨率為0.4 m,垂直分辨率為0.03 m;DMS用于識別海冰、海水和冰間水道等,空間分辨率為0.1 m;雪雷達用于測量積雪厚度,空間分辨率為40 m,垂直分辨率為0.06 m[18]。綜合上述3種傳感器,機載OIB可提供海冰的激光干舷高度和積雪厚度等數(shù)據(jù)。圖1展示了2017-2018年春季機載OIB北極飛行路線及激光干舷高度分布圖。

    對于NASA發(fā)布的機載OIB數(shù)據(jù),其海冰激光干舷高度、積雪厚度和產(chǎn)品空間分辨率都統(tǒng)一為40 m[18]。在此需要指出的是:機載OIB數(shù)據(jù)攜帶的測高儀器為激光雷達,因此,得到的海冰激光干舷高度并非海冰干舷高度,而是海冰干舷高度與積雪厚度之和。為得到海冰干舷高度,需將OIB得到的激光干舷高度減去雪雷達得到的積雪厚度,才能得到真實的海冰干舷高度。根據(jù)文獻記載,機載OIB的海冰激光干舷高度探測精度為0.015±0.06 m[19],積雪厚度的探測精度為0.01±0.05 m[20]。本文所用的OIB數(shù)據(jù)均與CS-2和S3雷達高度計數(shù)據(jù)時間、地點范圍相同(北極,2017年3月-4月和2018年4月)。由于CS-2、S3雷達高度計與OIB數(shù)據(jù)的空間分辨率不同,為進行逐一比較,將雷達高度計數(shù)據(jù)與OIB數(shù)據(jù)在25 km的網(wǎng)格中取平均,實現(xiàn)數(shù)據(jù)一一對應(yīng)。

    圖1 春季機載OIB北極飛行路線及激光干舷高度分布圖

    2 Ku波段雷達高度計波形重跟蹤閾值確定

    計算海冰干舷高度可通過反演海冰面到雷達質(zhì)心的高程和冰間水道面到雷達質(zhì)心的高程,并計算二者的差值得到。因此,首先需要區(qū)分雷達高度計的回波波形,以識別海冰、開闊水域和冰間水道等不同地物類型,然后再反演各種地物類型表面的高程。

    對于地物類型識別,主要是利用Ocean and Sea Ice SAF(OSI-SAF)發(fā)布的海冰密集度數(shù)據(jù)以及CS-2或S3的波形區(qū)分海冰、開闊水和冰間水道3種類型。常用后向散射系數(shù)(Radar Backscatter Coefficient,Sigma0)、波形前緣寬度(Leading Edge Width,LEW)和脈沖峰值(Pulse Peakiness,PP)3種回波波形特征進行類型識別[21]。Sigma0為雷達高度計接收到的地物的表面后向反射系數(shù);LEW是最大峰值功率5%和95%點位之間的波形前緣寬度;PP是雷達波形最大峰值功率與同一采樣波形里的所有波形總功率的比率[4],計算方法如下:

    式中:NWF代表一個波形內(nèi)的距離門數(shù);WF為雷達波形在第i個距離門處的功率。通常冰間水道的表面較平緩,雷達信號多為鏡面反射,PP值較高;對于開闊水域和海冰則以發(fā)生漫反射為主,PP值相對較低,LEW較寬(開闊水域的LEW相比海冰更大)。常用的分類參數(shù)值設(shè)置[22-23]見表1。

    表1 CryoSat-2和Sentinel-3A的雷達波形分類參數(shù)設(shè)置

    雷達高度計是通過計算信號發(fā)射與返回的時間差和光速來計算雷達到星下地物點之間的距離,因此常需要預(yù)先設(shè)置信號發(fā)射與返回時的時間跟蹤點。但在實際情況中,雷達接收信號的實際時間跟蹤點與預(yù)設(shè)時間跟蹤點之間存在偏差,需要根據(jù)實際情況進行波形重跟蹤校正[24]。本質(zhì)上,波形重跟蹤是將發(fā)生主散射的波形位置定為重跟蹤點,然后計算重跟蹤點和預(yù)設(shè)跟蹤點之間的偏移,以校正雷達高度計質(zhì)心到星下地物點之間的真實距離[25]。

    TFMRA方法(Threshold First Maximum Retracker Algorithm)[26]是海冰干舷高度反演中常用的波形重跟蹤方法。該方法針對冰間水道和海冰等不同的地物類型,經(jīng)驗性地設(shè)置重跟蹤點閾值位置。表2總結(jié)了當前最為主要的一些重跟蹤閾值組合。為了更全面地比較不同重跟蹤閾值組合對干舷反演的影響,并得到最優(yōu)閾值組合,本文在表2的3種閾值組合的基礎(chǔ)上重新構(gòu)建了8種閾值組合,分別對應(yīng)為(40%,40%),(50%,40%),(50%,50%),(60%,40%),(60%,50%),(70%,40%),(70%,50%),(70%,60%)(括號內(nèi)第一個閾值對應(yīng)于冰間水道,第二個閾值對應(yīng)于海冰)。以S3數(shù)據(jù)為例,這些閾值在雷達回波波形中的位置示意圖見圖2。同時,需要指出的是,冰間水道的閾值通常大于海冰的閾值,這是因為對于雷達高度計冰間水道的后向散射常高于海冰[27]。

    表2 文獻中常用的重跟蹤閾值組合

    圖2 Sentinel-3A雷達高度計海冰回波波形中第30~60距離門處歸一化回波功率閾值位置示意圖

    在完成重跟蹤閾值設(shè)置后,利用波形重跟蹤校正后的海冰表面高程減去相應(yīng)的平均海表面高度與海表面高程異常值,得到海冰雷達干舷高度FR(Radar Freeboard)。由于電磁波在穿過積雪層時的速度與真空中不同,所以還需要利用積雪厚度數(shù)據(jù)對雷達海冰干舷高度FR進行距離校正,以得到真正的海冰干舷高度F。因此,當雷達信號完全穿透積雪層時,進行積雪校正的海冰干舷高度計算公式[4]如下:

    式中:cs為積雪層中雷達信號的傳播速度;c為雷達信號在空氣中的傳播速度;hs(1-cs/c)≈0.22hs為積雪校正值。需要注意的是,公式(2)的使用條件為:電磁波能夠完全穿透雪層直達雪-冰交界面。若電磁波無法穿透雪層,應(yīng)用上式計算海冰干舷高度將會不可避免地引入誤差。

    圖3 2017年3月、4月與2018年4月北極地區(qū)不同閾值組合下海冰干舷高度與OIB海冰干舷高度散點圖

    表3 CryoSat-2衛(wèi)星數(shù)據(jù)不同閾值組合計算的海冰干舷高度與OIB海冰干舷高度對比

    將2017年3月、4月與2018年4月CS-2與S3在8種閾值組合方案下得到的海冰干舷高度分別與同期的OIB海冰干舷實測數(shù)據(jù)進行對比,從而確定兩種Ku波段合成孔徑雷達高度計針對海冰干舷高度反演的最優(yōu)波形重跟蹤閾值組合。為保證精度,在本文的處理中,積雪厚度數(shù)據(jù)采用的機載OIB數(shù)據(jù)提供的積雪厚度測量值。圖3和表3給出了不同閾值組合得到的海冰干舷高度與機載OIB得到的海冰干舷高度對比結(jié)果。

    從表3中可以看出,對于CS-2衛(wèi)星數(shù)據(jù),最優(yōu)的閾值組合為70%(冰間水道)、60%(海冰)。在該組合下得到的海冰干舷高度與機載OIB數(shù)據(jù)的海冰干舷高度非常接近。觀測表3可知,其平均絕對值差值和均方根誤差均為最小值,分別為0.060 7 m和0.077 6 m。

    對于S3衛(wèi)星數(shù)據(jù),由表4可知,最優(yōu)的閾值組合為50%(冰間水道)、50%(海冰)。在該組合下的平均絕對值差值和均方根誤差分別為0.075 0 m和0.098 2 m。

    同時,由表3、表4中的數(shù)據(jù)可以發(fā)現(xiàn),CS-2衛(wèi)星數(shù)據(jù)在最優(yōu)閾值組合下得到的海冰干舷高度,無論是反演精度還是相關(guān)性,均高于S3衛(wèi)星數(shù)據(jù)。

    表4 Sentinel-3A衛(wèi)星數(shù)據(jù)不同閾值組合計算的海冰干舷高度與OIB海冰干舷高度對比

    進一步分析圖3,可以發(fā)現(xiàn),無論是哪種閾值組合(即便采用本文得到的最優(yōu)閾值組合),均普遍存在雷達高度計探測的海冰干舷高度大于機載OIB探測的海冰干舷高度的情況。這說明對于Ku波段雷達,積雪層不能視為透明,需要考慮其無法穿透積雪層的情況。當積雪覆蓋在一年冰上時,海冰在析鹽過程中形成的薄鹽水層向上滲入到積雪層中,從而改變了積雪的介電特性,降低了雷達信號的積雪穿透能力[31]。對于覆蓋在多年冰上的積雪,不僅積雪表面的粗糙度更高[32],而且經(jīng)過多年的消融凍結(jié)循環(huán)過程,會形成包含多個不同密度層的積雪層[33]。這些因素綜合在一起,一方面會增加積雪層的后向散射強度,另一方面還會使雷達信號的主散射面由雪-冰交界面上移至積雪層中。這將導(dǎo)致海冰干舷高度被高估,從而引起海冰厚度的高估[14]。

    3 積雪層對海冰干舷高度反演的影響分析

    為詳細評估積雪層對雷達信號的影響,計算了CS-2與S3的雷達穿透系數(shù)f[13]:

    式中:d為積雪厚度hs與式(2)中的積雪校正值0.22hs之和;FR_sat為雷達高度計測的雷達海冰干舷高度;FR_OIB為OIB得到的雷達海冰干舷高度;OIB雷達海冰干舷高度由OIB的海冰干舷高度減去積雪校正值0.22hs得到。式(3)中,當雷達穿透系數(shù)f=0時,表示雷達信號沒有穿入雪層,雷達的主散射面位于空氣-雪交界面;f=1表示雷達信號完全穿透雪層,雷達信號的主散射面位于雪-冰交界面。0<f<1表示雷達信號未完全穿透雪層,雷達信號的主散射面位于積雪層內(nèi)部。

    圖4和表5給出了2017年3月-4月和2018年4月北極區(qū)域CS-2與S3各自的雷達穿透系數(shù)分布與統(tǒng)計結(jié)果。在圖4中還利用高斯分布對雷達穿透系數(shù)進行了擬合,并將雷達穿透系數(shù)的平均值作為擬合曲線的中心值。

    圖4 雷達穿透系數(shù)分布直方圖

    表5 CryoSat-2和Sentinel-3A的雷達穿透系數(shù)平均值

    需要注意的是,圖4存在雷達穿透系數(shù)f<0或者f>1的情況,對于該問題Armitage和Ridout[13]已給出過分析和評估。這主要是由于:(1)在計算海冰干舷時,為了有足夠的覆蓋范圍與機載OIB數(shù)據(jù)進行比較,衛(wèi)星數(shù)據(jù)的平均時間是遠長于機載OIB數(shù)據(jù)的;(2)衛(wèi)星高度計和機載OIB系統(tǒng)的足印大小存在較大差異。由于衛(wèi)星海冰干舷與機載OIB海冰干舷的兩種測量差異,使得衛(wèi)星高度計的海冰干舷分布范圍大于機載OIB的海冰干舷分布范圍,從而出現(xiàn)雷達穿透系數(shù)f不總在[0,1]范圍區(qū)間的情況。

    從圖4和表5可知,Ku波段雷達高度計對一年冰上積雪的穿透能力普遍強于多年冰。例如:CS-2對于一年冰上積雪的雷達穿透系數(shù)f為0.950,對于多年冰為0.889;S3對于一年冰上積雪的雷達穿透系數(shù)f為0.873,對于多年冰為0.856;對于全部海冰類型CS-2的雷達穿透系數(shù)f為0.912,S3的雷達穿透系數(shù)f為0.867。總體而言,CS-2的雷達穿透系數(shù)f大于 S3。

    所以普遍來說,對于Ku波段雷達高度計,其主散射面均位于雪-冰交界面上方的積雪層中,不能再利用公式(2)進行積雪校正,必須發(fā)展新的校正方法。

    4 積雪校正新方法與海冰干舷高度反演

    由雷達穿透系數(shù)f的定義可知,當0<f<1時,雷達高度計測得的海冰干舷高度會比實際的海冰干舷高度高(1-f)hs,所以需要補償這一距離差。同時,還需對雷達信號穿透的那部分積雪厚度進行距離校正0.22fhs。因此,改正的積雪校正模型可表述如下[13]:

    式中:雷達穿透系數(shù)f可取為表5中的統(tǒng)計值,所以式(4)可根據(jù)雷達高度計和海冰類型的不同,分別針對CS-2與S3,得到本文新改進的積雪校正公式。

    對于CS-2一年冰:

    對于CS-2多年冰:

    對于S3一年冰:

    對于S3多年冰:

    為了評估新改進的積雪校正方法的海冰干舷高度反演性能,本文以機載OIB實測海冰干舷高度數(shù)據(jù)為基礎(chǔ),將其與改進的積雪校正法的海冰干舷反演結(jié)果、通用積雪校正方法(公式(2))的反演結(jié)果和ESA發(fā)布的海冰干舷高度產(chǎn)品進行了對比分析。ESA海冰干舷高度產(chǎn)品是歐空局利用CS-2和S3數(shù)據(jù)制作的海冰干舷高度產(chǎn)品(L2級),該數(shù)據(jù)可公開下載(ftp://science-pds.cryosat.esa.int)。

    表6 通用積雪校正方法、提出的積雪校正方法、ESA海冰干舷高度產(chǎn)品與機載OIB數(shù)據(jù)的海冰干舷高度對比

    表6和圖5給出的是3種海冰干舷反演方法的結(jié)果與機載OIB海冰干舷高度的差異統(tǒng)計表和分布圖。根據(jù)圖與表的統(tǒng)計結(jié)果可知,對于ESA發(fā)布的海冰干舷高度產(chǎn)品,CS-2發(fā)布的海冰干舷高度產(chǎn)品值明顯偏大,而S3發(fā)布的海冰干舷高度產(chǎn)品值明顯偏小,但總體上S3海冰干舷高度產(chǎn)品的誤差低于CS-2。對于通用積雪校正方法,CS-2和S3反演的海冰干舷高度均偏高。而本文改進的積雪校正方法顯著校正了通用積雪校正方法存在的高估問題,CS-2和S3海冰干舷高度的反演精度無論是平均絕對值差值還是均方根誤差均小于前兩種方法。另外需要指出的是,本文所改進的積雪校正方法對于一年冰和多年冰,其精度均有所提高,呈現(xiàn)了良好的穩(wěn)定性。

    圖5 海冰干舷高度差異分布圖

    5 結(jié)論與討論

    工作在Ku波段的CryoSat-2和Sentinel-3A合成孔徑雷達高度計是當前最先進的高度計。本文針對雷達回波重跟蹤閾值的位置,以及積雪對雷達信號穿透能力的干擾影響,發(fā)展了改進的積雪校正算法,提高了海冰干舷高度反演精度。主要結(jié)論有:

    (1)對于CS-2雷達高度計最優(yōu)的波形重跟蹤閾值組合為70%(冰間水道)、60%(海冰);對于S3雷達高度計,其最優(yōu)的閾值組合為50%(冰間水道),50%(海冰)。

    (2)Ku波段雷達信號未完全穿透一年冰與多年冰上的積雪層,并且CS-2的雷達穿透系數(shù)大于S3。CS-2對于一年冰上積雪層的雷達穿透系數(shù)為0.95,對于多年冰為0.889;S3對于一年冰上積雪層的雷達穿透系數(shù)為0.873,對于多年冰為0.856;對于全部海冰類型而言,CS-2的雷達穿透系數(shù)為0.912,S3的雷達穿透系數(shù)為0.867。

    (3)相比于通用積雪校正法和ESA產(chǎn)品法,在本文的改進積雪校正法中得到的海冰干舷高度與OIB數(shù)據(jù)之間的誤差最小,精度最高。

    本文的研究僅使用了北極春季的數(shù)據(jù),為提高算法的普適性,未來將會利用北極其他季節(jié)的數(shù)據(jù)開展雷達穿透系數(shù)分析,從而進一步優(yōu)化算法。在下一步的研究工作中,準備利用本文改進的方法開展南極地區(qū)的海冰干舷高度反演與校正。

    致謝:感謝歐洲空間局提供的CryoSat-2和Sentinel-3A雷達高度計數(shù)據(jù)及海冰干舷高度數(shù)據(jù);感謝美國冰雪中心提供Operation IceBridge數(shù)據(jù)。

    猜你喜歡
    干舷高度計海冰
    基于漂流浮標的南大洋衛(wèi)星高度計有效波高研究
    海洋通報(2021年3期)2021-08-14 02:20:46
    末次盛冰期以來巴倫支海-喀拉海古海洋環(huán)境及海冰研究進展
    海洋通報(2021年3期)2021-08-14 02:20:38
    MIMU/GNSS/ODO/高度計/航姿儀組合導(dǎo)航微系統(tǒng)硬件設(shè)計
    航天控制(2020年5期)2020-03-29 02:10:36
    關(guān)于冬季和北大西洋冬季干舷的探討
    中國船檢(2019年11期)2019-12-26 03:18:40
    型深建造誤差對勘劃載重線標志的影響
    中國船檢(2019年10期)2019-10-31 07:07:40
    同化衛(wèi)星高度計觀測對CAS-ESM-C上層海洋溫度模擬的改進
    基于SIFT-SVM的北冰洋海冰識別研究
    合成孔徑雷達高度計與傳統(tǒng)高度計精度比對分析與機載試驗驗證
    載重線公約的規(guī)定、釋義及應(yīng)用(三)
    中國船檢(2015年4期)2015-05-31 09:27:14
    應(yīng)用MODIS數(shù)據(jù)監(jiān)測河北省近海海域海冰
    河北遙感(2014年4期)2014-07-10 13:54:59
    很黄的视频免费| 男人舔女人下体高潮全视频| 亚洲成人中文字幕在线播放| 岛国视频午夜一区免费看| 久久6这里有精品| 国产亚洲精品久久久com| 日本黄色片子视频| 一级毛片高清免费大全| 亚洲av成人不卡在线观看播放网| 午夜免费观看网址| 欧美成人一区二区免费高清观看| 欧美国产日韩亚洲一区| 亚洲国产精品合色在线| 亚洲欧美日韩无卡精品| 国产精品爽爽va在线观看网站| 欧美乱妇无乱码| 丰满人妻熟妇乱又伦精品不卡| 国产亚洲精品综合一区在线观看| 亚洲一区二区三区色噜噜| 久久婷婷人人爽人人干人人爱| 丰满乱子伦码专区| 男女之事视频高清在线观看| 久久香蕉国产精品| 成人亚洲精品av一区二区| 麻豆国产av国片精品| 老鸭窝网址在线观看| 啪啪无遮挡十八禁网站| 欧美成狂野欧美在线观看| 久久久久性生活片| 欧美日韩一级在线毛片| 最新在线观看一区二区三区| 成人鲁丝片一二三区免费| 中文字幕人妻熟人妻熟丝袜美 | 成人欧美大片| 真人一进一出gif抽搐免费| 19禁男女啪啪无遮挡网站| 欧美中文综合在线视频| 久久欧美精品欧美久久欧美| www.熟女人妻精品国产| 免费在线观看影片大全网站| 亚洲欧美日韩高清在线视频| 亚洲人成电影免费在线| 69av精品久久久久久| 在线观看66精品国产| 日韩精品青青久久久久久| 成人一区二区视频在线观看| 毛片女人毛片| 日本 av在线| 波多野结衣高清无吗| 亚洲国产欧洲综合997久久,| 亚洲人成电影免费在线| 亚洲在线观看片| 麻豆成人av在线观看| 男女做爰动态图高潮gif福利片| 免费大片18禁| 国产单亲对白刺激| 亚洲va日本ⅴa欧美va伊人久久| 国产精品久久久久久亚洲av鲁大| 欧美性感艳星| 色视频www国产| 国产激情偷乱视频一区二区| 亚洲中文字幕日韩| 三级男女做爰猛烈吃奶摸视频| 九九热线精品视视频播放| 精品一区二区三区人妻视频| netflix在线观看网站| aaaaa片日本免费| 久99久视频精品免费| 两人在一起打扑克的视频| xxxwww97欧美| 国产亚洲av嫩草精品影院| 国产精品女同一区二区软件 | 亚洲av电影不卡..在线观看| 久久久久久大精品| 久久久久性生活片| 深爱激情五月婷婷| 天美传媒精品一区二区| 天堂动漫精品| 99国产极品粉嫩在线观看| 国产老妇女一区| 内射极品少妇av片p| 可以在线观看毛片的网站| 麻豆一二三区av精品| 欧美区成人在线视频| 色噜噜av男人的天堂激情| 美女大奶头视频| 国产精品永久免费网站| 亚洲欧美激情综合另类| 俺也久久电影网| 午夜福利在线观看吧| 国产野战对白在线观看| 色老头精品视频在线观看| 国产精品99久久久久久久久| 一个人看的www免费观看视频| 国产熟女xx| 成熟少妇高潮喷水视频| 亚洲人与动物交配视频| 国产成人影院久久av| 啪啪无遮挡十八禁网站| 欧美日韩亚洲国产一区二区在线观看| 成熟少妇高潮喷水视频| av天堂中文字幕网| 脱女人内裤的视频| 观看免费一级毛片| 亚洲国产中文字幕在线视频| 国产97色在线日韩免费| 中文字幕高清在线视频| 在线视频色国产色| 国产三级黄色录像| 国产中年淑女户外野战色| 人人妻人人看人人澡| 人人妻人人澡欧美一区二区| 国产精品一区二区三区四区久久| 亚洲av免费在线观看| 国产成人a区在线观看| tocl精华| 精品一区二区三区视频在线 | 久久久国产成人精品二区| 亚洲精品久久国产高清桃花| 亚洲熟妇中文字幕五十中出| 国产黄a三级三级三级人| 久久中文看片网| 岛国视频午夜一区免费看| 少妇丰满av| 最好的美女福利视频网| 日韩人妻高清精品专区| 欧美最新免费一区二区三区 | 欧美日韩乱码在线| 免费看十八禁软件| 免费在线观看成人毛片| 精品久久久久久成人av| 男人和女人高潮做爰伦理| 岛国在线免费视频观看| 国产三级在线视频| 99国产精品一区二区三区| 日本黄色视频三级网站网址| 久久国产乱子伦精品免费另类| 国产精品一区二区三区四区久久| 黄片大片在线免费观看| 最新美女视频免费是黄的| 在线国产一区二区在线| 国产乱人视频| 99热只有精品国产| 听说在线观看完整版免费高清| 欧美乱妇无乱码| 国产色爽女视频免费观看| 久久香蕉精品热| 亚洲自拍偷在线| 噜噜噜噜噜久久久久久91| 午夜精品久久久久久毛片777| 成熟少妇高潮喷水视频| 国产激情欧美一区二区| 熟女电影av网| 在线观看一区二区三区| 免费大片18禁| 色视频www国产| 真人做人爱边吃奶动态| 熟女电影av网| 男人舔女人下体高潮全视频| 99国产综合亚洲精品| 桃色一区二区三区在线观看| 亚洲第一电影网av| 午夜日韩欧美国产| 亚洲电影在线观看av| 五月伊人婷婷丁香| 亚洲精品美女久久久久99蜜臀| 99久久无色码亚洲精品果冻| 国产单亲对白刺激| 成人高潮视频无遮挡免费网站| 老汉色∧v一级毛片| 黄片小视频在线播放| 免费搜索国产男女视频| 蜜桃亚洲精品一区二区三区| 国产综合懂色| 国产精品乱码一区二三区的特点| www日本在线高清视频| 亚洲电影在线观看av| 久久中文看片网| 一个人免费在线观看电影| 欧美日韩乱码在线| 麻豆国产av国片精品| 亚洲av电影不卡..在线观看| 蜜桃亚洲精品一区二区三区| 一级毛片高清免费大全| 国产老妇女一区| 一区福利在线观看| 手机成人av网站| 国产高清三级在线| 国产精华一区二区三区| 亚洲国产欧洲综合997久久,| 悠悠久久av| 可以在线观看毛片的网站| 三级毛片av免费| 欧美性猛交黑人性爽| 国产av在哪里看| 欧美乱色亚洲激情| 日本在线视频免费播放| 一本一本综合久久| 国产黄a三级三级三级人| 日韩欧美国产一区二区入口| 在线播放无遮挡| 中出人妻视频一区二区| 最近最新免费中文字幕在线| 午夜免费观看网址| 国产高清视频在线播放一区| 亚洲成av人片在线播放无| av片东京热男人的天堂| 亚洲av第一区精品v没综合| 国产伦精品一区二区三区视频9 | 国产精品99久久久久久久久| 国产单亲对白刺激| 日韩高清综合在线| 久久精品影院6| 免费看a级黄色片| 成人国产综合亚洲| 成人18禁在线播放| 男人的好看免费观看在线视频| 亚洲av美国av| 91字幕亚洲| 18禁黄网站禁片免费观看直播| 3wmmmm亚洲av在线观看| 99在线人妻在线中文字幕| 国产一区二区激情短视频| 国产成人影院久久av| 黄色片一级片一级黄色片| 亚洲第一欧美日韩一区二区三区| 天堂影院成人在线观看| 免费观看精品视频网站| 国产精品99久久久久久久久| 老汉色∧v一级毛片| 欧美午夜高清在线| 99国产综合亚洲精品| 精品熟女少妇八av免费久了| 亚洲第一电影网av| av黄色大香蕉| 欧美日本视频| 日本成人三级电影网站| 亚洲 欧美 日韩 在线 免费| 亚洲精品在线观看二区| 中文字幕熟女人妻在线| e午夜精品久久久久久久| 国产亚洲精品一区二区www| 国产一区在线观看成人免费| 女人十人毛片免费观看3o分钟| 琪琪午夜伦伦电影理论片6080| 亚洲真实伦在线观看| 精品福利观看| 少妇的逼水好多| 人妻久久中文字幕网| 操出白浆在线播放| 天天躁日日操中文字幕| 女生性感内裤真人,穿戴方法视频| АⅤ资源中文在线天堂| 国产真人三级小视频在线观看| 久久精品91蜜桃| 国产成+人综合+亚洲专区| 最新美女视频免费是黄的| 又爽又黄无遮挡网站| 男人和女人高潮做爰伦理| 两个人视频免费观看高清| 成年女人永久免费观看视频| av国产免费在线观看| 在线观看av片永久免费下载| 一级黄色大片毛片| 国产伦在线观看视频一区| 精品人妻一区二区三区麻豆 | 99视频精品全部免费 在线| 小蜜桃在线观看免费完整版高清| 欧美中文日本在线观看视频| 久久久国产成人免费| 欧美成人a在线观看| 偷拍熟女少妇极品色| 中文字幕人妻丝袜一区二区| 午夜福利欧美成人| 97碰自拍视频| 身体一侧抽搐| 亚洲成人精品中文字幕电影| 男女做爰动态图高潮gif福利片| 国产精品 国内视频| av黄色大香蕉| 亚洲最大成人中文| 精华霜和精华液先用哪个| 变态另类丝袜制服| 少妇人妻一区二区三区视频| 国产精品久久久久久久电影 | 日韩欧美一区二区三区在线观看| 此物有八面人人有两片| h日本视频在线播放| 熟妇人妻久久中文字幕3abv| 日韩欧美在线乱码| 国产激情偷乱视频一区二区| av黄色大香蕉| www.熟女人妻精品国产| 欧美成人免费av一区二区三区| bbb黄色大片| 午夜精品在线福利| 午夜老司机福利剧场| 91久久精品国产一区二区成人 | 男女视频在线观看网站免费| 深爱激情五月婷婷| 日日干狠狠操夜夜爽| 女人被狂操c到高潮| 国产伦精品一区二区三区四那| 亚洲专区国产一区二区| 又粗又爽又猛毛片免费看| 九色成人免费人妻av| 少妇裸体淫交视频免费看高清| 国产一区二区在线av高清观看| 一进一出好大好爽视频| 啪啪无遮挡十八禁网站| 成人18禁在线播放| 国产亚洲精品久久久久久毛片| 亚洲av免费在线观看| 一区二区三区免费毛片| 两个人看的免费小视频| 亚洲一区二区三区色噜噜| 长腿黑丝高跟| 免费人成视频x8x8入口观看| 99热这里只有是精品50| 久久国产乱子伦精品免费另类| 国产一区二区亚洲精品在线观看| 国产成人av教育| 精品一区二区三区av网在线观看| 亚洲av二区三区四区| 男人舔女人下体高潮全视频| 久9热在线精品视频| 制服人妻中文乱码| 欧美日韩黄片免| 欧美一区二区国产精品久久精品| 夜夜夜夜夜久久久久| 在线a可以看的网站| 他把我摸到了高潮在线观看| 偷拍熟女少妇极品色| 欧美乱妇无乱码| 99热这里只有是精品50| 成年女人看的毛片在线观看| 国产私拍福利视频在线观看| 欧美一级a爱片免费观看看| 少妇熟女aⅴ在线视频| 日本五十路高清| 久久精品影院6| 一区二区三区免费毛片| 中出人妻视频一区二区| 好看av亚洲va欧美ⅴa在| 99视频精品全部免费 在线| 五月玫瑰六月丁香| 精品久久久久久久末码| 黄色丝袜av网址大全| 免费无遮挡裸体视频| 国产高清激情床上av| 欧美黑人巨大hd| 亚洲天堂国产精品一区在线| 波野结衣二区三区在线 | 亚洲专区中文字幕在线| 老汉色av国产亚洲站长工具| 国产成人av教育| 少妇高潮的动态图| 亚洲18禁久久av| 午夜免费成人在线视频| 麻豆国产av国片精品| 中文字幕高清在线视频| 久久精品国产自在天天线| 天美传媒精品一区二区| 国产综合懂色| 国产高潮美女av| 亚洲18禁久久av| 国产成人aa在线观看| 亚洲精品一卡2卡三卡4卡5卡| 国产欧美日韩精品亚洲av| 99久久综合精品五月天人人| 久久国产精品影院| 男女做爰动态图高潮gif福利片| 精品欧美国产一区二区三| 国产探花极品一区二区| 日韩欧美三级三区| 99久久综合精品五月天人人| 亚洲 欧美 日韩 在线 免费| 午夜亚洲福利在线播放| 国产精品永久免费网站| 日本黄大片高清| 亚洲欧美一区二区三区黑人| 两个人看的免费小视频| 精品久久久久久,| 最近最新中文字幕大全免费视频| 国产 一区 欧美 日韩| 亚洲美女黄片视频| 精品一区二区三区视频在线 | 在线十欧美十亚洲十日本专区| 一区二区三区激情视频| 久久中文看片网| 欧美黄色淫秽网站| 成人高潮视频无遮挡免费网站| 日本五十路高清| 一边摸一边抽搐一进一小说| 人妻丰满熟妇av一区二区三区| 久久欧美精品欧美久久欧美| 99riav亚洲国产免费| 欧美最新免费一区二区三区 | 色播亚洲综合网| 99国产精品一区二区三区| 51午夜福利影视在线观看| 欧美精品啪啪一区二区三区| 色精品久久人妻99蜜桃| 亚洲人成网站在线播放欧美日韩| tocl精华| 此物有八面人人有两片| 女人被狂操c到高潮| av国产免费在线观看| 校园春色视频在线观看| 香蕉av资源在线| 一本精品99久久精品77| 色吧在线观看| 亚洲人成伊人成综合网2020| 91在线观看av| 成年版毛片免费区| 高清日韩中文字幕在线| 日本一二三区视频观看| 国产在线精品亚洲第一网站| 午夜影院日韩av| 国产精品久久视频播放| 级片在线观看| 男人舔奶头视频| 欧美性猛交黑人性爽| 很黄的视频免费| 一级毛片高清免费大全| 中文字幕人成人乱码亚洲影| 日本a在线网址| 欧美中文日本在线观看视频| 亚洲精品一卡2卡三卡4卡5卡| avwww免费| 欧美成人一区二区免费高清观看| 日本在线视频免费播放| 久久久久亚洲av毛片大全| 母亲3免费完整高清在线观看| xxxwww97欧美| netflix在线观看网站| 最好的美女福利视频网| 日韩欧美免费精品| 国内久久婷婷六月综合欲色啪| 日本五十路高清| 欧美色视频一区免费| 国产精品99久久99久久久不卡| 淫秽高清视频在线观看| 亚洲人成网站高清观看| 国产免费av片在线观看野外av| 蜜桃亚洲精品一区二区三区| 日本免费一区二区三区高清不卡| 国产黄a三级三级三级人| 久久久国产成人精品二区| 亚洲专区国产一区二区| 亚洲国产精品sss在线观看| 国产亚洲精品久久久com| 国产精品亚洲美女久久久| 在线视频色国产色| 成人av在线播放网站| 亚洲国产欧洲综合997久久,| 蜜桃亚洲精品一区二区三区| avwww免费| 亚洲国产精品久久男人天堂| 日韩欧美三级三区| 国产精品野战在线观看| 亚洲狠狠婷婷综合久久图片| 无限看片的www在线观看| 岛国在线免费视频观看| 国产精品亚洲一级av第二区| 午夜精品在线福利| 此物有八面人人有两片| 首页视频小说图片口味搜索| 国产主播在线观看一区二区| 在线播放国产精品三级| 岛国视频午夜一区免费看| 高潮久久久久久久久久久不卡| 亚洲18禁久久av| 成人一区二区视频在线观看| 露出奶头的视频| 色综合婷婷激情| 乱人视频在线观看| 又黄又粗又硬又大视频| 亚洲黑人精品在线| 欧美日韩精品网址| 18禁裸乳无遮挡免费网站照片| 动漫黄色视频在线观看| 亚洲人成电影免费在线| 日韩成人在线观看一区二区三区| 黄色丝袜av网址大全| 九色成人免费人妻av| 一本综合久久免费| 精品久久久久久久末码| 国内少妇人妻偷人精品xxx网站| 亚洲一区高清亚洲精品| 男女做爰动态图高潮gif福利片| 夜夜躁狠狠躁天天躁| 亚洲性夜色夜夜综合| 欧美激情久久久久久爽电影| 国产亚洲av嫩草精品影院| 91字幕亚洲| 精品午夜福利视频在线观看一区| 宅男免费午夜| 久久精品国产综合久久久| 精品久久久久久久毛片微露脸| 欧美一级a爱片免费观看看| 亚洲在线观看片| 伊人久久大香线蕉亚洲五| 欧美成人a在线观看| 欧美绝顶高潮抽搐喷水| 99国产综合亚洲精品| 成人欧美大片| 国产精品香港三级国产av潘金莲| 草草在线视频免费看| 美女高潮的动态| 亚洲真实伦在线观看| 欧美丝袜亚洲另类 | 黄色女人牲交| 两人在一起打扑克的视频| 欧美大码av| 亚洲乱码一区二区免费版| 热99re8久久精品国产| 12—13女人毛片做爰片一| 精品熟女少妇八av免费久了| 精品久久久久久,| 最新美女视频免费是黄的| 久久久精品欧美日韩精品| 亚洲色图av天堂| 亚洲人成电影免费在线| 男女下面进入的视频免费午夜| 高潮久久久久久久久久久不卡| 国产精品女同一区二区软件 | 日日摸夜夜添夜夜添小说| 欧美日韩乱码在线| 欧美日韩黄片免| 少妇的丰满在线观看| 宅男免费午夜| 成人永久免费在线观看视频| 99久久无色码亚洲精品果冻| 波野结衣二区三区在线 | 国产精品久久久久久久电影 | av黄色大香蕉| 国产不卡一卡二| 亚洲不卡免费看| 精品国产美女av久久久久小说| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品成人久久久久久| 欧美3d第一页| 日本五十路高清| 两人在一起打扑克的视频| 亚洲美女黄片视频| 久久精品国产99精品国产亚洲性色| 中文字幕人成人乱码亚洲影| 国产精品 欧美亚洲| 精品国产三级普通话版| 久久草成人影院| av欧美777| 精品久久久久久久久久久久久| 久久精品国产自在天天线| 99久久综合精品五月天人人| av国产免费在线观看| 国产精品久久久久久久久免 | 老熟妇仑乱视频hdxx| 特大巨黑吊av在线直播| 免费高清视频大片| 99国产精品一区二区三区| 久久九九热精品免费| 无遮挡黄片免费观看| 天堂影院成人在线观看| 免费看美女性在线毛片视频| 亚洲中文字幕一区二区三区有码在线看| 欧美黄色淫秽网站| 波多野结衣高清作品| 国产伦在线观看视频一区| 黄片大片在线免费观看| 色在线成人网| 看免费av毛片| 欧美zozozo另类| 大型黄色视频在线免费观看| 精品久久久久久成人av| 午夜福利免费观看在线| 婷婷精品国产亚洲av| 热99re8久久精品国产| 婷婷六月久久综合丁香| 看黄色毛片网站| 桃红色精品国产亚洲av| 日韩成人在线观看一区二区三区| 一级毛片女人18水好多| 欧美激情在线99| 日韩免费av在线播放| 99久久成人亚洲精品观看| 九九热线精品视视频播放| 中文亚洲av片在线观看爽| 夜夜夜夜夜久久久久| 亚洲熟妇中文字幕五十中出| 757午夜福利合集在线观看| 国产成人aa在线观看| 九九久久精品国产亚洲av麻豆| 在线观看舔阴道视频| 日本一二三区视频观看| 亚洲精品日韩av片在线观看 | 亚洲国产精品sss在线观看| 在线播放国产精品三级| 亚洲熟妇中文字幕五十中出| 亚洲国产精品成人综合色| 精品人妻一区二区三区麻豆 | 国产免费一级a男人的天堂| 国产淫片久久久久久久久 | av专区在线播放| 久久久国产精品麻豆| 亚洲午夜理论影院| 一夜夜www| 91久久精品电影网| 精品久久久久久成人av| 黄色视频,在线免费观看| 全区人妻精品视频| 色av中文字幕| 嫩草影院精品99| 一区二区三区激情视频| 精品99又大又爽又粗少妇毛片 |