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

    海洋白帽的高空間分辨率光學(xué)遙感估算分析

    2023-02-26 07:43:08趙碧丁靜劉建強(qiáng)焦俊男唐君陸應(yīng)誠(chéng)
    遙感學(xué)報(bào) 2023年1期
    關(guān)鍵詞:白帽覆蓋率海面

    趙碧,丁靜,劉建強(qiáng),焦俊男,唐君,陸應(yīng)誠(chéng)

    1.南京大學(xué) 國(guó)際地球系統(tǒng)科學(xué)研究所, 南京 210023;

    2.國(guó)家衛(wèi)星海洋應(yīng)用中心, 北京 100081;

    3.自然資源部空間海洋遙感與應(yīng)用重點(diǎn)實(shí)驗(yàn)室, 北京 100081

    1 引 言

    海面在風(fēng)的作用下產(chǎn)生海浪,隨著海表風(fēng)速增加到某一臨界值,海浪破碎并在波峰處產(chǎn)生大量的氣泡,從而形成目視特征顯著的白色氣泡水體,被稱為海洋白帽或白冠(Oceanic whitecaps,本文中統(tǒng)稱為白帽)(Raizer 等,1994;Melville,1996;Pivaev等,2021)。海洋白帽是海面波浪破碎的直觀表現(xiàn)(Callaghan 等,2008a;Salisbury 等,2014;Paget等,2015),是海氣交換的重要媒介和通道(Frouin 等,1996 和2001;Callaghan 等,2008a),在海氣通量交換(Asher 和Wanninkhof,1998;Goddijn-Murphy等,2016;Monahan等,1983)、海洋能量傳遞(Phillips,1985;Melville,1996)、海浪發(fā)展及波形破碎(Hanson 和Phillips,1999;Melville,1996;Myrhaug 和Holmedal,2008)等方面具有指示作用,是海洋環(huán)境監(jiān)測(cè)研究的重要對(duì)象。

    描述海洋白帽的屬性參數(shù)主要有白帽反射率、白帽覆蓋率等,海洋白帽在光學(xué)到微波譜段內(nèi)的反射特性基本得到闡明(Koepke,1984;Kokhanovsky,2004;Moore 等,2000;Nicolas 等,2001),其在可見(jiàn)光、近紅外波段的反射率高于背景水體,因此海洋白帽在現(xiàn)場(chǎng)觀測(cè)資料中易于目視識(shí)別和提取(Koepke,1984;Frouin 等,1996,2001)。海洋白帽的現(xiàn)場(chǎng)觀測(cè)分析基于原位觀測(cè)的照片或視頻資料,通過(guò)統(tǒng)計(jì)一定范圍內(nèi)白帽面積和發(fā)生頻率等特征,從而揭示該區(qū)域環(huán)境動(dòng)力的現(xiàn)象與影響因素(Monahan等,1983;Goddijn-Murphy等,2011;Callaghan 等,2008a;Jia 和Zhao,2019;Zhang 等,2015)?;诂F(xiàn)場(chǎng)觀測(cè)資料能進(jìn)一步給出白帽覆蓋率,即白帽面積占總觀測(cè)面積的比例(Whitecap coverage 或Whitecap fraction)從而可建立起白帽覆蓋率與風(fēng)速之間的統(tǒng)計(jì)關(guān)系模型(Monahan 和Muircheartaigh,1980;Callaghan 等,2008b;Lafon等,2004)。不同海域的“白帽覆蓋率—風(fēng)速”模型研究,進(jìn)一步揭示了風(fēng)對(duì)海面(波浪破碎、海面粗糙度、海氣交換等)的作用機(jī)制、影響程度及區(qū)域分異特點(diǎn)(Monahan等,1983;Goddijn-Murphy等,2011;Callaghan等,2008a;Jia和Zhao,2019;Zhang等,2015)。

    衛(wèi)星光學(xué)遙感技術(shù)開(kāi)展海洋白帽的觀測(cè)分析,是厘清不同區(qū)域海洋白帽的分布變化、揭示海洋環(huán)境動(dòng)力特征的有效技術(shù)支撐(Koepke,1984;Kokhanovsky,2004;Moore 等,2000;Nicolas 等,2001)。在海洋白帽遙感研究中,除了云霧的影響,圖像空間分辨率、海洋水色差異、海面耀光反射、海洋環(huán)境動(dòng)力等均會(huì)產(chǎn)生影響(Monahan等,1983;Goddijn-Murphy等,2011;Callaghan等,2008a;Myrhaug 和Holmedal,2008)。如在粗空間分辨率的MODIS 數(shù)據(jù)中,難以直接從圖像中鑒別海洋白帽,目前多利用統(tǒng)計(jì)模型來(lái)估算白帽反射信號(hào)貢獻(xiàn),提高M(jìn)ODIS 數(shù)據(jù)的大氣校正精度(Gordon和Wang,1994;Hu等,2000),因此,要從海洋遙感的高空間分辨率數(shù)據(jù)中(海洋遙感與陸地遙感在空間分辨率的需求與認(rèn)知上存在差異,不展開(kāi)討論)進(jìn)行白帽的識(shí)別提取、厘清海洋白帽的光學(xué)遙感特征。海洋白帽能在美國(guó)陸地衛(wèi)星(Landsat 8)30 m 空間分辨率的陸地成像儀OLI(Operational Land Imager)圖像中清晰識(shí)別(Kubryakov 等,2021),也能在歐空局哨兵2 號(hào)衛(wèi)星(Sentinel-2 A/B)10 m 空間分辨率的多光譜成像儀MSI(Multi Spectral Instrument)圖像中清晰展現(xiàn)(如圖1所示)。

    圖1 海洋白帽的照片與遙感影像Fig.1 Oceanic whitecaps in photograph and remotely sensed images

    開(kāi)展高空間分辨率海洋光學(xué)遙感數(shù)據(jù)中海洋白帽的識(shí)別提取與定量反演研究,挖掘光學(xué)遙感數(shù)據(jù)在海洋環(huán)境動(dòng)力(如風(fēng)場(chǎng)、流場(chǎng)等)監(jiān)測(cè)中的應(yīng)用潛力,能拓展海洋光學(xué)遙感的研究范疇。本研究基于Sentinel-2 MSI 高分辨率數(shù)據(jù),進(jìn)行海面白帽覆蓋率的光學(xué)遙感反演,探索白帽覆蓋率對(duì)海面環(huán)境動(dòng)力背景的響應(yīng),分析白帽覆蓋率反演算法的適用條件、不確定性和干擾因素;面向中國(guó)首個(gè)海洋水色業(yè)務(wù)衛(wèi)星星座(海洋一號(hào)C/D衛(wèi)星,HY-1C/D)的高空間分辨率光學(xué)載荷——海岸帶成像儀CZI(Coastal Zone Imager),開(kāi)展不同風(fēng)速背景下,CZI 圖像中白帽的可識(shí)別性研究,分析白帽反射信號(hào)在CZI各波段中的信號(hào)貢獻(xiàn)差異,討論海洋白帽在CZI數(shù)據(jù)中的存在特點(diǎn);從而為國(guó)產(chǎn)自主海洋光學(xué)衛(wèi)星的精細(xì)應(yīng)用、海洋環(huán)境動(dòng)力與海氣交互監(jiān)測(cè)提供參考。

    2 數(shù)據(jù)與方法

    2.1 數(shù)據(jù)及預(yù)處理

    選用Sentinel-2 衛(wèi)星的MSI 數(shù)據(jù),其10 m 空間分辨率的影像上能夠清晰識(shí)別出海洋白帽。相對(duì)應(yīng)的背景風(fēng)速數(shù)據(jù),選擇歐洲中期天氣預(yù)報(bào)中心ECMWF(European Centre for Medium-Range Weather Forecasts)的ERA5 再分析風(fēng)速產(chǎn)品,其空間分辨率為0.25°、時(shí)間分辨率為1 h。數(shù)據(jù)選擇及預(yù)處理需基于如下幾點(diǎn)考慮:(1)先根據(jù)背景風(fēng)速優(yōu)選MSI 數(shù)據(jù)。白帽的產(chǎn)生需要一定的風(fēng)速條件,據(jù)現(xiàn)場(chǎng)統(tǒng)計(jì),在6 m/s 以上的風(fēng)速會(huì)產(chǎn)生海洋白帽,風(fēng)速大于7—8 m/s 時(shí)白帽迅速增多(Monahan 和Zietlow,1969),為了保證在MSI 影像中檢測(cè)到白帽,本研究選擇8 m/s 以上風(fēng)速條件的影像。(2)對(duì)于高空間分辨率數(shù)據(jù),太陽(yáng)耀光反射的影響較粗空間分辨率數(shù)據(jù)更為復(fù)雜(Wu 等,2022)。為了盡量避免耀光反射的影響,選擇MSI數(shù)據(jù)的θm大于40°的影像。θm為海面太陽(yáng)鏡面反射角與傳感器觀測(cè)角度之間的夾角,由傳感器與太陽(yáng)的幾何角度決定,是海面耀光反射強(qiáng)弱的直接指針,其大小可表征平均耀光反射的強(qiáng)弱,θm越大表示耀光反射越弱(Hu 等,2009;Lu 等,2016;Wen等,2018),θm的計(jì)算公式如下:

    式中,θ0、θ和φ分別為太陽(yáng)天頂角、傳感器天頂角和相對(duì)方位角。(3)數(shù)據(jù)預(yù)處理。利用ACOLITE(版本20190326)軟件對(duì)MSI L1C 級(jí)別數(shù)據(jù)進(jìn)行瑞利校正獲得10 m 分辨率的數(shù)據(jù)(Rrc)。鑒于云會(huì)明顯干擾白帽覆蓋率的估算,因此選擇無(wú)云或少云的影像,對(duì)圖像中少量的云進(jìn)行掩膜,并對(duì)掩膜結(jié)果進(jìn)行5個(gè)像元的膨脹處理。

    對(duì)中國(guó)東黃海2021年MSI 數(shù)據(jù)進(jìn)行篩選,研究主要使用MSI 842 nm 波段生成Rrc反射率產(chǎn)品,可忽略背景水體差異,具體區(qū)域和影像時(shí)間如圖2所示。圖2(a)、2(b)、2(c)分別由2景、5景、2 景MSI 影像拼接而成。由于影像數(shù)據(jù)的成像時(shí)間在2:16—2:41 UTC,因此將風(fēng)速產(chǎn)品2:00 UTC 及3:00 UTC的數(shù)據(jù)取平均。

    2.2 海面白帽覆蓋率估算

    光學(xué)衛(wèi)星數(shù)據(jù)中,大氣頂層總輻亮度Lt可表示為

    式中,Lr、La和Lra分別表示瑞利散射、氣溶膠散射、瑞利—?dú)馊苣z散射輻亮度,LWC示白帽反射輻亮度,Lw表示離水輻亮度,Lg表示耀光反射輻亮度,T、t分別為大氣直射、散射透過(guò)率(Gordon和Wang,1994)。根據(jù)公式:

    可將輻亮度轉(zhuǎn)換為反射率,式(3)中R為反射率,L為輻亮度,F(xiàn)0為太陽(yáng)輻照度,θ0為太陽(yáng)高度角。海面反射率又分為受白帽影響及不受影響部分(Koepk,1984),結(jié)合式(2),海面反射率可表達(dá)為

    式中,Rr、Ra和Rra分別表示瑞利散射、氣溶膠散射、瑞利—?dú)馊苣z散射反射率,Rwc示白帽反射率,Rw表示離水反射率,Rg表示耀光反射率,T、t分別為大氣直射、散射透過(guò)率,W為白帽覆蓋率(Whitecap coverage),表示海面被白帽覆蓋的部分占總面積的比例,由式(4)可得

    式中,Rb=Ra+Rra+tRw,代表了來(lái)自背景海水與大氣的反射率貢獻(xiàn)。由于近紅外波段(842 nm)中離水輻亮度可以忽略(即tRw≈0),且忽略耀光反射的影響(圖2 中數(shù)據(jù)的θm均大于40°,則TRg= 0)。因此,針對(duì)圖2 中白帽覆蓋率可根據(jù)下式(6)進(jìn)行估算:

    圖2 Sentinel-2衛(wèi)星MSI近紅外波段Rrc(842 nm)影像圖Fig.2 Sentinel-2 MSI near infrared Rrc(842 nm)images

    近紅外波段白帽的反射率(Rwc)約為0.55(Koepke,1984;Frouin 等,1996 和2001),海洋白帽的光學(xué)特征相對(duì)穩(wěn)定,具有較為明確的光譜反射率特征(Kokhanovsky,2004;Garaba 和Dierssen,2018),而圖像上白帽反射率變化則是與圖像中白帽覆蓋率密切相關(guān)。近紅外波段大氣散射透過(guò)率(t)為0.7—0.8(Kneizys 等,1981),上述值僅可用于背景海水與大氣均勻穩(wěn)定且不受海面耀光反射干擾的區(qū)域。此外,利用式(6)進(jìn)行白帽覆蓋率(W)反演的關(guān)鍵在于準(zhǔn)確的背景反射信號(hào)Rb的獲取。本研究中,假定在4 km×4 km的區(qū)域內(nèi),海水和大氣的光學(xué)屬性是均勻的,對(duì)圖像進(jìn)行400×400(4 km×4 km)的窗口濾波,獲得窗口內(nèi)的最小值作為背景反射率,可進(jìn)一步計(jì)算得出圖像的白帽覆蓋率(W)。特別需要說(shuō)明的是,上述公式和參數(shù)有著嚴(yán)格的適用條件,本文原理分析中的數(shù)據(jù)邊界清楚,可以適用,但是在未來(lái)的實(shí)際業(yè)務(wù)應(yīng)用中,還需要發(fā)展業(yè)務(wù)化應(yīng)用算法。

    3 基于MSI影像的白帽覆蓋率估算與分析

    3.1 白帽覆蓋率

    基于MSI 842 nm 波段Rrc數(shù)據(jù)(圖2)反演的白帽覆蓋率如圖3所示,其中云、陸地、輻射沙脊等已經(jīng)被掩膜處理。白帽覆蓋率作為海面波浪破碎狀況的直觀體現(xiàn),一方面反映了風(fēng)速這一主要?jiǎng)恿σ蛩?,與風(fēng)速大小呈正相關(guān);另一方面,又受到海洋渦旋、內(nèi)波、鋒面等的調(diào)制(Monahan等,1983;Goddijn-Murphy 等,2011;Callaghan等,2008a)。如在圖3(a)中,2021年1月13日的白帽覆蓋率的變化主要受到海面風(fēng)速的影響,而在圖3(b)和3(c)中,在2021年1月27日和2021年2月21日的白帽覆蓋率在局部區(qū)域又受到海洋鋒面、渦旋等的影響(具體分析見(jiàn)3.3 節(jié))??傮w而言,圖3 中平均白帽覆蓋率在1%上下(對(duì)應(yīng)風(fēng)速8—12 m/s),這與白帽的現(xiàn)場(chǎng)觀測(cè)(基于海面實(shí)拍照片或視頻等資料)分析結(jié)果較為接近(Blanchard,1983)。

    圖3 基于MSI數(shù)據(jù)的白帽覆蓋率估算結(jié)果Fig.3 MSI images derived whitecap coverage

    為了驗(yàn)證白帽覆蓋率算法的估算精度,本文針對(duì)圖2(a)中紅框框選出的12 個(gè)區(qū)域開(kāi)展精度評(píng)價(jià)。12 個(gè)500×400 像元驗(yàn)證區(qū)域的842 nm 波段Rrc圖像見(jiàn)圖4 第一和第三列,對(duì)這12 個(gè)區(qū)域圖像采用人工目視解譯的方法,通過(guò)選取合適閾值分割出圖像中的海洋白帽,并計(jì)算出該區(qū)域中的白帽覆蓋率(圖4 第二、第四列)。將12 個(gè)驗(yàn)證區(qū)域的算法反演結(jié)果與目視解譯結(jié)果進(jìn)行對(duì)比分析(如圖5所示),算法反演的白帽覆蓋率與目視解譯的白帽覆蓋率具有較好的一致性(R2=0.85),說(shuō)明提取結(jié)果的精度可靠驗(yàn)證擬合線低于1∶1 線,說(shuō)明本研究算法在弱耀光條件下存在一定的低估。

    3.2 基于白帽覆蓋率的海面風(fēng)速反演

    海面白帽是風(fēng)生波浪破碎的直觀表現(xiàn),在其他海洋環(huán)境動(dòng)力影響較為均勻的區(qū)域,白帽覆蓋率與海面風(fēng)速(海面10 m 高處的風(fēng)速,U10)具有密切的統(tǒng)計(jì)關(guān)系(Monahan 和Muircheartaigh,1980;Anguelova 和Webster,2006;Anguelova 和Bettenhausen,2019)。基于海面實(shí)測(cè)數(shù)據(jù),構(gòu)建“白帽覆蓋率—風(fēng)速”統(tǒng)計(jì)模型,有利于對(duì)區(qū)域海洋環(huán)境動(dòng)力等的深入研究。需要指出的是,本研究的重點(diǎn)在于能否反演出海面風(fēng)速,并與已有風(fēng)速產(chǎn)品進(jìn)行對(duì)比分析,但基于本區(qū)域海面實(shí)測(cè)結(jié)果的“白帽覆蓋率—風(fēng)速”統(tǒng)計(jì)模型缺乏,因此選擇Monahan 和Muircheartaigh(1980)的白帽覆蓋率—風(fēng)速模型來(lái)反演海面風(fēng)速。模型公式如下:

    式中,W代表白帽覆蓋率,U10為風(fēng)速,對(duì)其進(jìn)行變換可以給出風(fēng)速反演模型如下:

    利用式(8)結(jié)合圖3(a)的海面白帽覆蓋率提取結(jié)果,可以得到該區(qū)域風(fēng)速反演結(jié)果(如圖6(a),空間分辨率為4 km),圖6(b)為該區(qū)域同時(shí)段的ERA5 海面再分析風(fēng)速產(chǎn)品(時(shí)間差<15 min,空間分辨率為0.25°)。對(duì)比圖6(a)和6(b),整體風(fēng)速為東北方向風(fēng)速高,向西南方向風(fēng)速遞減,MSI推導(dǎo)的風(fēng)速結(jié)果因?yàn)榻K沿海輻射沙脊的影響會(huì)出現(xiàn)一些反演異常值(見(jiàn)圖6(a)左下角)。總體上,基于MSI 白帽覆蓋率反演的海面風(fēng)速與ERA5海面風(fēng)速產(chǎn)品的空間變化趨勢(shì)與風(fēng)速變化區(qū)間基本一致。給出二者的頻率累計(jì)曲線對(duì)比如圖6(c),用于進(jìn)一步評(píng)估風(fēng)速反演的精度。海面風(fēng)速的頻率累計(jì)曲線圖可以說(shuō)明,MSI白帽覆蓋率反演的海面風(fēng)速與ERA5 海面風(fēng)速產(chǎn)品變化基本一致,但MSI 反演的海面風(fēng)速略低于ERA5 風(fēng)速產(chǎn)品的值,二者之間的差異在風(fēng)速為10 m/s左右達(dá)到最大,為1.6 m/s,大部分差異在1 m/s 以內(nèi),這主要是由于算法對(duì)MSI圖像中白帽覆蓋率的低估造成(圖5)。

    圖5 目視解譯與算法估算白帽覆蓋率對(duì)比圖Fig.5 Comparison of visual interpretation and algorithmic estimation of whitecap coverage from MSI images

    圖6 MSI白帽反演風(fēng)速及其與ERA5風(fēng)速產(chǎn)品的對(duì)比Fig.6 Comparison between whitecaps derived wind speed based on MSI images and wind speed product of ERA5

    3.3 其他環(huán)境動(dòng)力對(duì)白帽覆蓋率的調(diào)制

    海洋白帽主要反映的是風(fēng)生波浪的破碎過(guò)程,雖然海面風(fēng)場(chǎng)是白帽產(chǎn)生的首要因素,但風(fēng)生波浪也會(huì)被渦旋、內(nèi)波等其他海洋環(huán)境動(dòng)力所調(diào)制,從而改變白帽的發(fā)生概率(Thorpe 等,1987;Romero,2019)。圖7分別反映了海洋鋒面、渦旋、內(nèi)波和船只尾跡對(duì)白帽覆蓋率的影響:圖7左列為MSIRrc(842 nm)影像(圖7(a)與圖7(b)左列分別為圖2(b)和2(c)中紅框區(qū)域,圖7(c)左列為圖8(a)中紅框區(qū)域);圖7右列為白帽覆蓋率。

    MSI 842 nm 的開(kāi)闊海洋觀測(cè)數(shù)據(jù)中,雖然該數(shù)據(jù)的觀測(cè)角度避開(kāi)了強(qiáng)耀光反射(圖2中數(shù)據(jù)的θm均大于40°),但在高空間分辨率光學(xué)影像中,海面耀光反射差異實(shí)際上仍存在,并對(duì)影像總信號(hào)有貢獻(xiàn)(Wu 等,2022)。如圖7(a)所示,海洋鋒線(圖中紅色虛線)兩側(cè)呈現(xiàn)明顯的水團(tuán)圖像特征差異,圖7(a)中白帽覆蓋率在不同水團(tuán)之間也表現(xiàn)出類(lèi)似的差異,雖然海水離水輻亮度差異可以忽略,但是白帽覆蓋率差異和耀光反射差異同時(shí)存在;圖7(b)和(c)反映了受海洋渦旋與內(nèi)波的影響,白帽覆蓋率會(huì)存在一定程度的升高。易于理解的是,相同風(fēng)速條件下,內(nèi)波的波峰處較其他位置,波浪更易于破碎而產(chǎn)生白帽(Thorpe等,1987;Kubryakov等,2021)。

    圖7 不同海洋環(huán)境動(dòng)力對(duì)白帽的影響Fig.7 Influence of different marine environmental dynamics on whitecaps

    海洋白帽是風(fēng)生波浪破碎的產(chǎn)物,這種波浪破碎也受到海面鋒面、渦旋、內(nèi)波等各種其他海洋動(dòng)力因素的調(diào)制,在高空間分辨率的MSI 影像上,海洋白帽及其覆蓋率所反映的海洋與大氣背景動(dòng)力特征能清晰展現(xiàn)。

    3.4 耀光反射對(duì)白帽覆蓋率估算的影響

    海面耀光反射差異一直是海洋光學(xué)遙感研究關(guān)注的重點(diǎn)問(wèn)題(Cox 和Munk,1954;Jackson 和Alpers,2010;Kay 等,2009;Harmel 等,2018)?;诮?jīng)典統(tǒng)計(jì)模型(Cox 和Munk,1954)估算的耀光反射率,雖可以準(zhǔn)確描述粗空間分辨率衛(wèi)星光學(xué)遙感數(shù)據(jù)(如250 m 的MODIS 數(shù)據(jù))中耀光反射的強(qiáng)弱;但在高空間分辨率的海洋光學(xué)遙感數(shù)據(jù)中(如10 m 的MSI數(shù)據(jù)),耀光反射規(guī)律會(huì)有所不同(Wu 等,2022),表現(xiàn)出較大的離散性,無(wú)法通過(guò)統(tǒng)計(jì)模型或簡(jiǎn)單的空間濾波進(jìn)行剔除(Jackson 和Alpers,2010)。本研究白帽覆蓋率估算算法中的4 km × 4 km 窗口濾波方法,本質(zhì)是一種平滑處理(主要用于消除海水和大氣的影響),不能直接從總體反射率中有效消除耀光反射的影響,因此無(wú)法處理強(qiáng)耀光反射的圖像。

    為了闡明耀光反射對(duì)海面白帽覆蓋率估算的影響,選擇2021年6月2日的6 景MSI 數(shù)據(jù)(數(shù)據(jù)預(yù)處理與圖2 和圖3 一致)。如圖8(a)中所示,θm從左到右逐漸減小(由32°減小到16°),耀光反射不斷增強(qiáng);圖8(b)為算法估算的白帽覆蓋率,與耀光強(qiáng)度空間變化相似;圖8(c)為該區(qū)域同時(shí)段的ERA5風(fēng)速產(chǎn)品,在無(wú)明顯水團(tuán)、渦旋等海洋動(dòng)力環(huán)境的調(diào)制下,風(fēng)速與白帽覆蓋率的空間分異顯然不同(圖8(b))。此外,在相似風(fēng)速下(8—12 m/s),圖8(b)中強(qiáng)耀光反射下的白帽覆蓋率較圖3(a)的結(jié)果有一個(gè)數(shù)量級(jí)上的差異;對(duì)圖8(a)與圖8(b)做采樣分析,以θm為橫軸(圖8(d)),不難發(fā)現(xiàn)樣點(diǎn)的Rrc(842 nm)反射率與白帽覆蓋率(W)的變化趨勢(shì)一致。

    圖8 Sentinel-2 MSI影像中白帽覆蓋率與耀光變化分析Fig.8 Whitecap coverage and sunglint variation in Sentinel-2 MSI images

    海面耀光反射是白帽識(shí)別估算的主要干擾,本研究算法只適用于弱耀光反射下(θm大于40°)的MSI圖像,如何有效消除耀光反射的影響,也是今后亟待關(guān)注的重點(diǎn)問(wèn)題之一。

    4 CZI數(shù)據(jù)中海洋白帽的響應(yīng)

    HY-1C/D 是中國(guó)首個(gè)海洋水色業(yè)務(wù)衛(wèi)星星座,CZI載荷能提供星下點(diǎn)空間分辨率為50 m、中心波長(zhǎng)位于460 nm、560 nm、650 nm 和825 nm 的4 波段海洋光學(xué)遙感數(shù)據(jù)(沈亞峰 等,2020)。本文所用的CZI 數(shù)據(jù)均為L(zhǎng)2A 級(jí)的瑞利校正反射率產(chǎn)品(Rrc),由于像元的混合效應(yīng),相對(duì)較低風(fēng)速條件下的白帽難以被CZI圖像有效識(shí)別,隨著風(fēng)速的增加,海洋白帽逐漸被有效檢測(cè)出,即CZI圖像中檢測(cè)出海洋白帽的最低風(fēng)速要大于MSI圖像。

    圖9 為選取的開(kāi)闊海洋CZIRrc(825 nm)圖像,圖像都為海面弱耀光反射條件,海洋環(huán)境動(dòng)力影響不明顯,且涵蓋不同風(fēng)速條件下的影像,該區(qū)域同時(shí)段的ERA5風(fēng)速均已經(jīng)標(biāo)注在Rrc(825 nm)圖像上。對(duì)CZI圖像上的海洋白帽進(jìn)行識(shí)別并計(jì)算其覆蓋率,圖9(a)和圖9(b)中海洋白帽不能被有效識(shí)別;在較高風(fēng)速下(圖9(c)、(d)、(e)、(f))海洋白帽能夠從CZI 圖像上檢測(cè)出。50 m 空間分辨率CZI圖像中可識(shí)別的海面白帽,其所對(duì)應(yīng)的海面風(fēng)速約9 m/s 以上,隨著風(fēng)速的增加,其白帽覆蓋率也隨之增加(圖9(c)、(d)、(e)、(f)),這也反映了CZI數(shù)據(jù)具有對(duì)高風(fēng)速下的海面白帽覆蓋率進(jìn)行遙感監(jiān)測(cè)與估算的應(yīng)用潛力。

    圖9 不同風(fēng)速條件下HY-1C/D CZI影像中的白帽及其目視解譯結(jié)果Fig.9 Whitecaps in HY-1C/D CZI images and corresponding visual interpretation results at different wind speed(U10)

    為了進(jìn)一步分析闡明CZI圖像中的海洋白帽反射率與背景海水反射率的差異,本研究對(duì)圖9 中4 幅CZI 圖像(圖9(c)、(d)、(e)、(f))的4 個(gè)波段(460 nm、560 nm、650 nm 和825 nm),對(duì)其海洋白帽反射率與背景海水反射率分別進(jìn)行統(tǒng)計(jì)分析,結(jié)果如圖10所示。CZI 的4 個(gè)波段圖像中,白帽反射率均明顯高于背景海水反射率,反射率差異分別為5.8%、7.8%、8.3%和7.3%。

    圖10 HY-1C/D CZI不同波段的海洋白帽與背景海水反射率差異Fig.10 Reflection difference between whitecaps and background seawater in different HY-1C/D CZI bands

    5 結(jié) 論

    海洋白帽是海面風(fēng)生波浪破碎的直觀表現(xiàn),其區(qū)域性特征參數(shù)——白帽覆蓋率,不僅與海面風(fēng)速有著密切的統(tǒng)計(jì)關(guān)系,也會(huì)受到其他海洋環(huán)境動(dòng)力的調(diào)制。海洋白帽對(duì)入射光具有較高的反射散射特征,是海洋光學(xué)遙感數(shù)據(jù)的大氣校正一個(gè)重要環(huán)節(jié)。隨著海洋光學(xué)遙感數(shù)據(jù)空間分辨率的提高,海洋白帽能被有效的探測(cè),從而開(kāi)展更為豐富的研究與應(yīng)用探索。本文針對(duì)Sentinel-2 MSI數(shù)據(jù),分析海面白帽反射率的信號(hào)貢獻(xiàn),利用一種用于弱耀光反射條件下的白帽覆蓋率估算方法,開(kāi)展MSI圖像上白帽覆蓋率反演?;贛SI圖像估算的不同區(qū)域海面白帽覆蓋率,與海面風(fēng)速的變化具有良好的空間一致性和統(tǒng)計(jì)關(guān)系,基于白帽覆蓋率估算的海面風(fēng)速與ERA5再分析海洋風(fēng)速產(chǎn)品對(duì)比驗(yàn)證表明,其相對(duì)偏差不大于1.6 m/s。此外,不同水團(tuán)、鋒面、渦旋、內(nèi)波等其他海洋環(huán)境動(dòng)力對(duì)海面白帽的調(diào)制作用,能顯著改變遙感數(shù)據(jù)反演的白帽覆蓋率;與此同時(shí),MSI圖像中耀光反射對(duì)海面白帽遙感的干擾與影響也不容忽視。最后,基于HY-1C/D 衛(wèi)星CZI 載荷數(shù)據(jù)開(kāi)展海面白帽的響應(yīng)分析,結(jié)果表明:海面風(fēng)速大于9 m/s 時(shí),在無(wú)其他海洋環(huán)境動(dòng)力的調(diào)制下,海面白帽可被50 m 空間分辨率的CZI 載荷所探測(cè),且白帽覆蓋率隨著風(fēng)速的增加而增加;不同CZI波段圖像中,白帽反射率與背景海水反射率有5.8%—8.3%的差異。這表明,HY-1C/D 衛(wèi)星CZI 載荷可成為海洋白帽深入研究的可用數(shù)據(jù)源,能用于較高風(fēng)速下的海面風(fēng)速反演??傊S著海洋光學(xué)遙感載荷空間分辨率的提升,通過(guò)海洋白帽的有效探測(cè),能為海面風(fēng)場(chǎng)及其他海洋環(huán)境動(dòng)力監(jiān)測(cè)提供新的技術(shù)方法參考,也必將在海面風(fēng)場(chǎng)反演、海氣交換、海洋動(dòng)力環(huán)境研究等方面發(fā)揮越來(lái)越重要的作用。

    志 謝感謝美國(guó)南佛羅里達(dá)大學(xué)海洋科學(xué)學(xué)院胡傳民教授、武漢大學(xué)測(cè)繪遙感信息工程國(guó)家重點(diǎn)實(shí)驗(yàn)室田禮喬教授、中國(guó)科學(xué)院煙臺(tái)海岸帶研究所邢前國(guó)研究員對(duì)本研究的指導(dǎo);感謝歐空局和歐洲中期天氣預(yù)報(bào)中心提供的MSI和ERA5數(shù)據(jù)(https://scihub.copernicus.eu/dhus/#/home,https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-single-levels?tab=form);感謝國(guó)家衛(wèi)星海洋應(yīng)用中心提供的CZI 數(shù)據(jù)(http://www.nsoas.org.cn/)。

    猜你喜歡
    白帽覆蓋率海面
    民政部等16部門(mén):到2025年村級(jí)綜合服務(wù)設(shè)施覆蓋率超80%
    鳥(niǎo)
    我國(guó)全面實(shí)施種業(yè)振興行動(dòng) 農(nóng)作物良種覆蓋率超過(guò)96%
    海面床,輕輕搖
    第六章 邂逅“胖胖號(hào)”
    孤 獨(dú)(外一首)
    繞轉(zhuǎn)小舌頭
    海面上的“一千座埃菲爾鐵塔”
    基于噴丸隨機(jī)模型的表面覆蓋率計(jì)算方法
    涉嫌者的謊話
    亚洲自偷自拍三级| 天堂俺去俺来也www色官网| 日本黄色日本黄色录像| 日日摸夜夜添夜夜爱| 精品熟女少妇av免费看| 国产欧美日韩一区二区三区在线 | 国产乱人偷精品视频| 人妻夜夜爽99麻豆av| av卡一久久| 亚洲国产精品999| 五月玫瑰六月丁香| 亚洲国产精品国产精品| 亚洲美女黄色视频免费看| 嫩草影院入口| 日本wwww免费看| 亚洲最大成人中文| 中文乱码字字幕精品一区二区三区| 视频区图区小说| 蜜桃在线观看..| 男人和女人高潮做爰伦理| 大陆偷拍与自拍| 美女福利国产在线 | 国产 一区 欧美 日韩| 亚洲精品日韩av片在线观看| 国产片特级美女逼逼视频| 国产黄片美女视频| 国产黄片视频在线免费观看| 精品视频人人做人人爽| 各种免费的搞黄视频| 老熟女久久久| 国产女主播在线喷水免费视频网站| 亚洲国产最新在线播放| 十八禁网站网址无遮挡 | 三级经典国产精品| 久久综合国产亚洲精品| 久久这里有精品视频免费| 国产国拍精品亚洲av在线观看| 精品亚洲成a人片在线观看 | 噜噜噜噜噜久久久久久91| 欧美一级a爱片免费观看看| 精品久久久久久电影网| 国产精品久久久久久久久免| 国产亚洲5aaaaa淫片| 狂野欧美激情性bbbbbb| 美女脱内裤让男人舔精品视频| 老司机影院成人| 久久精品国产鲁丝片午夜精品| 一级毛片 在线播放| 女人久久www免费人成看片| 啦啦啦中文免费视频观看日本| 免费看光身美女| 黄色配什么色好看| 高清在线视频一区二区三区| 日日啪夜夜爽| 亚洲自偷自拍三级| 精品一区二区三卡| 大片电影免费在线观看免费| 精品一区二区三卡| 国产精品欧美亚洲77777| 欧美97在线视频| 高清av免费在线| 97超视频在线观看视频| 成人毛片a级毛片在线播放| 国产探花极品一区二区| 久久精品国产鲁丝片午夜精品| 国产精品99久久99久久久不卡 | 久久久久视频综合| 国产伦理片在线播放av一区| 建设人人有责人人尽责人人享有的 | 亚洲精品久久久久久婷婷小说| 少妇 在线观看| 久久99精品国语久久久| 国产亚洲一区二区精品| av播播在线观看一区| 国产精品久久久久久精品古装| 五月伊人婷婷丁香| 国产91av在线免费观看| 亚洲欧洲国产日韩| 亚洲国产精品国产精品| 国产一区亚洲一区在线观看| 国产亚洲av片在线观看秒播厂| 国产精品国产三级国产av玫瑰| 亚洲第一av免费看| 黑人高潮一二区| 秋霞伦理黄片| 99九九线精品视频在线观看视频| 一个人看视频在线观看www免费| 国产av国产精品国产| 女性生殖器流出的白浆| 青春草亚洲视频在线观看| 亚洲无线观看免费| 国产免费视频播放在线视频| 99精国产麻豆久久婷婷| 大香蕉97超碰在线| 亚洲精品乱码久久久久久按摩| 国产乱来视频区| 亚洲精品日韩av片在线观看| 欧美区成人在线视频| 亚洲第一av免费看| 国产精品国产三级国产专区5o| av国产精品久久久久影院| 日本av免费视频播放| 在线观看人妻少妇| 亚洲精品国产成人久久av| 国产成人午夜福利电影在线观看| 偷拍熟女少妇极品色| 亚洲中文av在线| 亚洲四区av| 午夜福利高清视频| 欧美丝袜亚洲另类| 尾随美女入室| 亚洲欧美中文字幕日韩二区| 亚洲av电影在线观看一区二区三区| 精品人妻视频免费看| 久久久精品免费免费高清| 国产亚洲一区二区精品| 18禁裸乳无遮挡免费网站照片| 黄色怎么调成土黄色| 免费观看a级毛片全部| 欧美xxⅹ黑人| 特大巨黑吊av在线直播| 啦啦啦啦在线视频资源| 国产成人91sexporn| 91精品一卡2卡3卡4卡| 成人漫画全彩无遮挡| 成人亚洲欧美一区二区av| 欧美xxxx性猛交bbbb| 亚洲欧美成人综合另类久久久| 免费黄频网站在线观看国产| 免费久久久久久久精品成人欧美视频 | 少妇丰满av| 观看免费一级毛片| 免费av中文字幕在线| 国产精品一区二区三区四区免费观看| 六月丁香七月| 天天躁夜夜躁狠狠久久av| 成年av动漫网址| 日本猛色少妇xxxxx猛交久久| 99re6热这里在线精品视频| 一级毛片我不卡| 一区二区三区精品91| 内地一区二区视频在线| 国产亚洲欧美精品永久| 色综合色国产| 国产精品嫩草影院av在线观看| 亚洲真实伦在线观看| 日日啪夜夜爽| 成人18禁高潮啪啪吃奶动态图 | 亚洲高清免费不卡视频| av又黄又爽大尺度在线免费看| 男人舔奶头视频| 啦啦啦啦在线视频资源| 菩萨蛮人人尽说江南好唐韦庄| 中文资源天堂在线| 久久久久网色| av播播在线观看一区| 联通29元200g的流量卡| 美女内射精品一级片tv| 一区二区三区精品91| 男人添女人高潮全过程视频| 久久国产亚洲av麻豆专区| 777米奇影视久久| 精品午夜福利在线看| a级毛片免费高清观看在线播放| 最近中文字幕高清免费大全6| 国产高清有码在线观看视频| 婷婷色综合www| 麻豆国产97在线/欧美| 最近最新中文字幕免费大全7| 又粗又硬又长又爽又黄的视频| 在线观看三级黄色| 午夜激情久久久久久久| 国产亚洲91精品色在线| 精品久久久久久久久亚洲| 男的添女的下面高潮视频| 不卡视频在线观看欧美| 插逼视频在线观看| 久久精品人妻少妇| 只有这里有精品99| 国产精品一及| 日韩中文字幕视频在线看片 | 亚洲精品国产色婷婷电影| 亚洲精品国产色婷婷电影| 亚洲国产精品一区三区| 亚洲,一卡二卡三卡| a级毛片免费高清观看在线播放| 有码 亚洲区| 日本猛色少妇xxxxx猛交久久| 亚洲欧美精品自产自拍| 三级经典国产精品| 黄色配什么色好看| 亚洲色图av天堂| 亚洲精品一二三| 精品酒店卫生间| 亚洲性久久影院| 亚洲精品一区蜜桃| 色视频www国产| 最近手机中文字幕大全| 有码 亚洲区| 少妇的逼好多水| 一级毛片久久久久久久久女| 草草在线视频免费看| 99热国产这里只有精品6| 日本-黄色视频高清免费观看| 最近中文字幕2019免费版| 成人毛片60女人毛片免费| 看非洲黑人一级黄片| 亚洲av男天堂| 亚洲精品国产成人久久av| 好男人视频免费观看在线| 水蜜桃什么品种好| 国产视频内射| 狠狠精品人妻久久久久久综合| 久久久a久久爽久久v久久| 久久久欧美国产精品| 黄色欧美视频在线观看| 久久这里有精品视频免费| 99精国产麻豆久久婷婷| 亚洲av欧美aⅴ国产| 一本色道久久久久久精品综合| 亚洲精品日本国产第一区| 午夜视频国产福利| 日韩伦理黄色片| 国语对白做爰xxxⅹ性视频网站| 中文在线观看免费www的网站| 亚洲一区二区三区欧美精品| 日日摸夜夜添夜夜爱| 久久久久久久久久久丰满| 天堂中文最新版在线下载| 97在线视频观看| 免费在线观看成人毛片| 永久网站在线| 丰满迷人的少妇在线观看| 如何舔出高潮| 男女免费视频国产| 精品一区在线观看国产| 内地一区二区视频在线| av在线蜜桃| 国产亚洲最大av| 免费看日本二区| 大香蕉久久网| 久久久久久久精品精品| 国产日韩欧美亚洲二区| 成人毛片60女人毛片免费| 久久综合国产亚洲精品| 免费黄频网站在线观看国产| 久久99精品国语久久久| 国产成人a∨麻豆精品| 国产精品麻豆人妻色哟哟久久| 91午夜精品亚洲一区二区三区| 视频中文字幕在线观看| 中文在线观看免费www的网站| 男女啪啪激烈高潮av片| 日本wwww免费看| 日韩制服骚丝袜av| 91久久精品国产一区二区三区| 国内精品宾馆在线| 久久久久精品性色| 国产精品久久久久久精品古装| 亚洲国产精品国产精品| av专区在线播放| 夫妻性生交免费视频一级片| 一级毛片黄色毛片免费观看视频| 精品国产一区二区三区久久久樱花 | 一边亲一边摸免费视频| 日韩大片免费观看网站| 日韩制服骚丝袜av| 国产男女内射视频| 我要看黄色一级片免费的| 欧美xxxx黑人xx丫x性爽| 国产在线一区二区三区精| 亚洲在久久综合| 一区二区三区精品91| 青青草视频在线视频观看| 在线免费观看不下载黄p国产| 国产免费又黄又爽又色| 亚洲国产高清在线一区二区三| 18禁裸乳无遮挡动漫免费视频| 亚洲av免费高清在线观看| 久热这里只有精品99| 久久99精品国语久久久| 在线免费十八禁| 欧美人与善性xxx| 99久久中文字幕三级久久日本| 51国产日韩欧美| 亚洲精品一二三| 18禁裸乳无遮挡动漫免费视频| 久久精品国产亚洲网站| 欧美日韩一区二区视频在线观看视频在线| 一级a做视频免费观看| 国产乱人偷精品视频| 男人舔奶头视频| 亚洲成人手机| 亚洲无线观看免费| 亚洲国产av新网站| 久久精品国产a三级三级三级| 边亲边吃奶的免费视频| 久久精品久久精品一区二区三区| 国产 一区 欧美 日韩| 免费人妻精品一区二区三区视频| 亚洲精品国产av成人精品| 免费人成在线观看视频色| 精品人妻视频免费看| 午夜老司机福利剧场| 欧美变态另类bdsm刘玥| 久久久久久伊人网av| 亚洲综合精品二区| 日韩国内少妇激情av| 高清av免费在线| 中文字幕制服av| av在线蜜桃| 久久97久久精品| 久久精品熟女亚洲av麻豆精品| 美女主播在线视频| 亚洲精品第二区| 三级国产精品欧美在线观看| 天美传媒精品一区二区| 黄色欧美视频在线观看| 久久精品国产自在天天线| 尤物成人国产欧美一区二区三区| 男人舔奶头视频| 国产精品久久久久成人av| 国产精品99久久久久久久久| 国产男人的电影天堂91| 国产欧美日韩一区二区三区在线 | 亚洲欧美一区二区三区国产| 国产成人freesex在线| 在线观看三级黄色| 久久精品国产a三级三级三级| 日韩人妻高清精品专区| 国模一区二区三区四区视频| 欧美变态另类bdsm刘玥| a级毛片免费高清观看在线播放| 国产乱人偷精品视频| 一区二区三区精品91| 午夜老司机福利剧场| 欧美变态另类bdsm刘玥| 色吧在线观看| 国产精品伦人一区二区| 日韩一本色道免费dvd| 一区二区三区免费毛片| 国产成人freesex在线| 又大又黄又爽视频免费| 国产久久久一区二区三区| 日产精品乱码卡一卡2卡三| 王馨瑶露胸无遮挡在线观看| 欧美老熟妇乱子伦牲交| 啦啦啦在线观看免费高清www| 欧美 日韩 精品 国产| 午夜福利高清视频| 国产在视频线精品| 看十八女毛片水多多多| av一本久久久久| 亚洲精品一区蜜桃| 免费高清在线观看视频在线观看| 亚洲欧美清纯卡通| 国产欧美亚洲国产| 久久ye,这里只有精品| 又黄又爽又刺激的免费视频.| av国产免费在线观看| 亚洲人与动物交配视频| 男女免费视频国产| 能在线免费看毛片的网站| 人妻一区二区av| 麻豆国产97在线/欧美| av在线app专区| 欧美日韩视频精品一区| 少妇 在线观看| 肉色欧美久久久久久久蜜桃| 天天躁日日操中文字幕| 国产精品99久久久久久久久| 久久ye,这里只有精品| 中文在线观看免费www的网站| 自拍偷自拍亚洲精品老妇| 高清欧美精品videossex| 国产欧美亚洲国产| 在线观看国产h片| tube8黄色片| 日韩中字成人| 啦啦啦啦在线视频资源| 日韩av不卡免费在线播放| 日本av手机在线免费观看| 在线观看免费视频网站a站| 日本猛色少妇xxxxx猛交久久| 乱码一卡2卡4卡精品| 一区在线观看完整版| 国产 精品1| 国产熟女欧美一区二区| 日日摸夜夜添夜夜爱| 人人妻人人看人人澡| 日韩大片免费观看网站| 纯流量卡能插随身wifi吗| 亚洲精品自拍成人| 毛片女人毛片| 国产黄片美女视频| 亚洲欧美日韩卡通动漫| 欧美人与善性xxx| 久久久久久久久大av| 超碰97精品在线观看| 男的添女的下面高潮视频| 亚洲国产最新在线播放| 黄色配什么色好看| 欧美日韩视频高清一区二区三区二| 青春草视频在线免费观看| 视频区图区小说| 干丝袜人妻中文字幕| 欧美精品人与动牲交sv欧美| 99久久精品国产国产毛片| 国产精品一区二区三区四区免费观看| 欧美性感艳星| 欧美变态另类bdsm刘玥| xxx大片免费视频| 免费黄网站久久成人精品| 天天躁夜夜躁狠狠久久av| 夫妻午夜视频| videossex国产| 中文字幕亚洲精品专区| 午夜福利在线在线| 最后的刺客免费高清国语| 岛国毛片在线播放| 在线观看美女被高潮喷水网站| 日韩电影二区| 亚洲国产欧美在线一区| 国产亚洲av片在线观看秒播厂| 婷婷色麻豆天堂久久| 久久99热这里只频精品6学生| av黄色大香蕉| 亚洲精品自拍成人| 啦啦啦视频在线资源免费观看| 熟妇人妻不卡中文字幕| 亚洲国产高清在线一区二区三| 美女国产视频在线观看| 一本一本综合久久| 中文字幕人妻熟人妻熟丝袜美| 男女免费视频国产| 国产视频首页在线观看| 又黄又爽又刺激的免费视频.| 日本wwww免费看| 国产精品一二三区在线看| 边亲边吃奶的免费视频| 简卡轻食公司| 欧美精品一区二区大全| 最近手机中文字幕大全| 色网站视频免费| 亚洲色图综合在线观看| 久久av网站| 一级毛片电影观看| 在线观看一区二区三区激情| 亚洲精品乱久久久久久| 国产v大片淫在线免费观看| 国产精品国产三级专区第一集| 狂野欧美白嫩少妇大欣赏| 亚洲欧美精品自产自拍| 亚洲四区av| 日韩av在线免费看完整版不卡| 性色av一级| 男女啪啪激烈高潮av片| 十分钟在线观看高清视频www | 欧美一级a爱片免费观看看| 精品一区在线观看国产| 少妇高潮的动态图| 免费大片黄手机在线观看| 亚洲国产欧美在线一区| 国产高清三级在线| 有码 亚洲区| 午夜视频国产福利| 午夜激情久久久久久久| 成人综合一区亚洲| 黑人高潮一二区| 高清在线视频一区二区三区| 国产成人免费观看mmmm| 亚洲成人一二三区av| 欧美日韩在线观看h| 午夜激情福利司机影院| 亚洲第一区二区三区不卡| 午夜老司机福利剧场| 亚洲成色77777| 啦啦啦视频在线资源免费观看| 十八禁网站网址无遮挡 | 久久韩国三级中文字幕| 大香蕉久久网| 亚洲欧美日韩卡通动漫| 国产在线视频一区二区| 天堂8中文在线网| 亚洲无线观看免费| 亚洲国产精品一区三区| kizo精华| 舔av片在线| 国产成人午夜福利电影在线观看| 色哟哟·www| 成人一区二区视频在线观看| 亚洲av福利一区| 一级爰片在线观看| 国产一区亚洲一区在线观看| 国产精品麻豆人妻色哟哟久久| 97超视频在线观看视频| 亚洲第一av免费看| 在线观看美女被高潮喷水网站| 大陆偷拍与自拍| 97精品久久久久久久久久精品| 少妇丰满av| 男人添女人高潮全过程视频| 最近中文字幕高清免费大全6| 欧美日韩一区二区视频在线观看视频在线| a级毛色黄片| www.色视频.com| h视频一区二区三区| 色视频在线一区二区三区| 国产有黄有色有爽视频| 久久久久精品性色| 日韩人妻高清精品专区| 国产精品国产三级国产专区5o| 国产一区二区三区av在线| 特大巨黑吊av在线直播| 26uuu在线亚洲综合色| 亚洲美女搞黄在线观看| 视频区图区小说| 国产欧美日韩一区二区三区在线 | 精品一品国产午夜福利视频| 亚洲国产欧美在线一区| 夜夜爽夜夜爽视频| 777米奇影视久久| 中文字幕久久专区| 2018国产大陆天天弄谢| 黄色欧美视频在线观看| 精品久久久久久久久av| 亚洲av中文字字幕乱码综合| 少妇的逼好多水| 午夜视频国产福利| 久久ye,这里只有精品| 欧美老熟妇乱子伦牲交| 我要看日韩黄色一级片| 亚洲,欧美,日韩| 欧美zozozo另类| 黑丝袜美女国产一区| 国产在视频线精品| 国产精品偷伦视频观看了| 新久久久久国产一级毛片| 亚洲av电影在线观看一区二区三区| 欧美区成人在线视频| 免费在线观看成人毛片| 成人国产av品久久久| 精品人妻一区二区三区麻豆| 久久久久人妻精品一区果冻| 女性生殖器流出的白浆| 亚洲精品日本国产第一区| 国产一区二区在线观看日韩| 日韩电影二区| av卡一久久| 91午夜精品亚洲一区二区三区| 直男gayav资源| 国国产精品蜜臀av免费| 网址你懂的国产日韩在线| 97精品久久久久久久久久精品| 国产精品一区二区在线观看99| www.av在线官网国产| 亚洲高清免费不卡视频| 国产极品天堂在线| 国产色爽女视频免费观看| 色视频在线一区二区三区| 国产成人精品久久久久久| 欧美 日韩 精品 国产| 国产男女内射视频| 欧美最新免费一区二区三区| 亚洲精品成人av观看孕妇| 国产精品一及| 高清视频免费观看一区二区| 久久国内精品自在自线图片| 亚洲美女黄色视频免费看| 中文字幕久久专区| 91精品伊人久久大香线蕉| 婷婷色麻豆天堂久久| 久久久a久久爽久久v久久| 日本爱情动作片www.在线观看| 日韩欧美 国产精品| 狠狠精品人妻久久久久久综合| 国产精品三级大全| 成人免费观看视频高清| 亚洲精品成人av观看孕妇| 一级毛片 在线播放| 亚洲在久久综合| 丝袜脚勾引网站| 日韩制服骚丝袜av| 人人妻人人爽人人添夜夜欢视频 | 国产精品一区二区三区四区免费观看| 欧美精品亚洲一区二区| 精品一区二区免费观看| 午夜视频国产福利| 欧美成人一区二区免费高清观看| 中文欧美无线码| 成人一区二区视频在线观看| 美女福利国产在线 | 在线 av 中文字幕| 一区二区av电影网| 国产成人午夜福利电影在线观看| 免费观看的影片在线观看| 欧美日韩精品成人综合77777| 久久精品熟女亚洲av麻豆精品| 91精品伊人久久大香线蕉| 免费人妻精品一区二区三区视频| 成人综合一区亚洲| 久久青草综合色| 国内少妇人妻偷人精品xxx网站| 伊人久久国产一区二区| 天美传媒精品一区二区| 日韩一本色道免费dvd| 观看免费一级毛片| 国产免费视频播放在线视频| 只有这里有精品99| 尤物成人国产欧美一区二区三区| 人人妻人人添人人爽欧美一区卜 | 婷婷色综合www| 99热国产这里只有精品6| 日韩成人伦理影院| 国产淫片久久久久久久久|