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

    中國陸地區(qū)域GF-1 WFV地表反射率產(chǎn)品

    2023-10-19 02:48:16佘文慶張兆明彭燕何國金龍騰飛王桂周
    遙感學(xué)報(bào) 2023年9期
    關(guān)鍵詞:頂角反射率波段

    佘文慶,張兆明,3,彭燕,何國金,3,龍騰飛,3,王桂周

    1.中國科學(xué)院空天信息創(chuàng)新研究院,北京 100094;2.中國科學(xué)院大學(xué),北京 100049;3.海南空天信息研究院海南省地球觀測重點(diǎn)實(shí)驗(yàn)室,三亞 572029

    1 引言

    高分一號(hào)衛(wèi)星(以下簡稱GF-1)是中國高分辨率對(duì)地觀測系統(tǒng)的首顆衛(wèi)星,于2013年4月26日在酒泉衛(wèi)星發(fā)射中心成功發(fā)射。GF-1 衛(wèi)星搭載兩臺(tái)2 m分辨率的全色、8 m分辨率的多光譜相機(jī)和4臺(tái)16 m分辨率的寬幅相機(jī)。其中4臺(tái)寬幅相機(jī)影像拼接合成后,幅寬可達(dá)到800 km(白照廣,2013)。目前GF-1 衛(wèi)星數(shù)據(jù)已廣泛應(yīng)用于林業(yè)、農(nóng)業(yè)和生態(tài)環(huán)境等領(lǐng)域,在一定程度上彌補(bǔ)了中國高分辨率衛(wèi)星數(shù)據(jù)的應(yīng)用短板。

    地表反射率是表征物體表面對(duì)太陽輻射反射的能力,是許多生態(tài)遙感參數(shù)反演的輸入數(shù)據(jù),例如生物量和葉面積指數(shù)等。目前國內(nèi)外公開共享的地表反射率產(chǎn)品主要有250 m、500 m、1 km分辨率的MODIS地表反射率產(chǎn)品(Liang等,2002;Vermote 等,2002)、30 m 分辨率的Landsat 系列地表反射率產(chǎn)品(Roy 等,2010;Vermote 等,2016;彭燕 等,2020)和10 m、20 m、60 m 分辨率的Sentinel-2地表反射率產(chǎn)品(Main-Knorn等,2017)。國產(chǎn)GF-1 WFV數(shù)據(jù)已面向公眾開放共享,但是僅提供L1 級(jí)標(biāo)準(zhǔn)化產(chǎn)品。如果能在GF-1 WFV L1 級(jí)數(shù)據(jù)的基礎(chǔ)上生產(chǎn)16 m 分辨率的地表反射率產(chǎn)品,即可與Landsat系列、Sentinel系列的地表反射率產(chǎn)品協(xié)同形成中國區(qū)域近每日的密集時(shí)間序列中高分辨率地表反射率產(chǎn)品。

    地表反射率反演的關(guān)鍵步驟是大氣校正,大氣校正方法主要分為統(tǒng)計(jì)模型和物理模型(李正強(qiáng) 等,2018)。統(tǒng)計(jì)模型方法是利用地表變量和遙感數(shù)據(jù)之間建立的統(tǒng)計(jì)關(guān)系模型進(jìn)行大氣校正,但由于區(qū)域之間存在差異性,統(tǒng)計(jì)模型通常適用于局部地區(qū)(鄭偉和曾志遠(yuǎn),2004;申茜 等,2021)。物理模型方法的基本原理是利用電磁波在大氣中的輻射傳輸原理建立起來的模型對(duì)遙感影像進(jìn)行大氣校正,其大氣校正精度高且應(yīng)用范圍廣(鄭偉和曾志遠(yuǎn),2004)。基于物理模型的遙感影像大氣校正的一個(gè)核心問題是如何獲取成像時(shí)的氣溶膠光學(xué)厚度AOD(Aerosol Optical Depth)、大氣水蒸汽含量WV(Water Vapor)和臭氧含量TO(Total Ozone)等大氣參數(shù)。通常的做法是采用相關(guān)的大氣產(chǎn)品或利用自身數(shù)據(jù)進(jìn)行反演。如Hu 等(2014)和Peng 等(2016)分別將MOD04(AOD)、MOD05(WV)和MOD07(TO)應(yīng)用于Landsat 5和Landsat 8 影像的大氣校正,地面實(shí)測地表反射率數(shù)據(jù)驗(yàn)證表明了算法的可行性。Martins等(2018)將MCD19A2(AOD 和WV)和MOD08D3(TO)應(yīng)用于CBERS-4 MUX 影像的大氣校正,生產(chǎn)南美洲的CBERS-4 MUX 地表反射率產(chǎn)品,利用Landsat 8 OLI 地表反射率和AERONET 站點(diǎn)數(shù)據(jù)驗(yàn)證了產(chǎn)品的精度,驗(yàn)證結(jié)果表明產(chǎn)品具有較高的可靠性。美國WELD(Web Enabled Landsat Data)團(tuán)隊(duì)利用MODIS 大氣產(chǎn)品(MOD04、MOD05 和MOD07)生產(chǎn)美國區(qū)域的Landsat 地表反射率產(chǎn)品(Roy 等,2010)。這些研究表明了將相關(guān)的MODIS大氣產(chǎn)品應(yīng)用于大區(qū)域遙感影像大氣校正的可行性。在大氣參數(shù)反演方面,大氣水蒸汽含量反演需要數(shù)據(jù)具有水汽吸收通道(Main-Knorn 等,2017),不適用于GF-1 WFV 數(shù)據(jù);氣溶膠光學(xué)厚度反演常用的算法為暗目標(biāo)(Dark Target,DT)法,該算法是利用短波紅外波段受大氣影響較小的特性,根據(jù)暗目標(biāo)區(qū)域短波紅外波段與藍(lán)、紅波段地表反射率之間的經(jīng)驗(yàn)線性關(guān)系來反演AOD(Kaufman 和Sendra,1988),但由于GF-1 WFV 數(shù)據(jù)缺少短波紅外波段,導(dǎo)致DT 方法應(yīng)用于GF-1 WFV 數(shù)據(jù)存在較大困難(Richter等,2006)。

    因此利用MODIS 大氣產(chǎn)品輔助GF-1 WFV 大氣校正成為了一個(gè)重要途徑。然而MODIS 大氣產(chǎn)品在部分區(qū)域存在無效值,尤其是AOD 產(chǎn)品,在一定程度上會(huì)影響地表反射率反演的精度。針對(duì)此問題,通常的做法是直接對(duì)單一AOD 產(chǎn)品采用傳統(tǒng)克里金法插值填充無效值區(qū)域(Hu等,2014;胡勇 等,2018)。該方法雖然能在一定程度上解決無效值問題,但是其計(jì)算復(fù)雜,且對(duì)大量無效值區(qū)域的插值會(huì)出現(xiàn)數(shù)據(jù)失真的現(xiàn)象。目前已發(fā)布和共享了多種不同分辨率的MODIS AOD 產(chǎn)品,若能將這些多源產(chǎn)品在空間上進(jìn)行融合,就可以減少大范圍無效值的情況。

    本文在參考相關(guān)工作的基礎(chǔ)上,開展多源MODIS AOD 產(chǎn)品協(xié)同的GF-1 WFV 大氣校正方法研究。基于多源MODIS AOD 空間融合和動(dòng)態(tài)查找表,設(shè)計(jì)和實(shí)現(xiàn)一種針對(duì)GF-1 WFV 的6S 大氣校正算法,并利用該算法生產(chǎn)和共享中國陸地區(qū)域2020 年地表反射率產(chǎn)品。該產(chǎn)品可為全國植被生長監(jiān)測、地表覆蓋分類和重要生態(tài)遙感參數(shù)反演等研究提供基礎(chǔ)數(shù)據(jù)支撐。

    2 數(shù)據(jù)源介紹

    2.1 GF-1 WFV數(shù)據(jù)

    GF-1 衛(wèi)星傳感器由2 臺(tái)多光譜相機(jī)(PMS)和4臺(tái)寬視場相機(jī)(WFV)組成。本文所用的數(shù)據(jù)為GF-1 WFV,其參數(shù)如下表1 所示。本文選取2020 年覆蓋中國陸地區(qū)域的無云或者少云的GF-1 WFV 正射數(shù)據(jù),總共850 景。所有數(shù)據(jù)來源于中國資源衛(wèi)星應(yīng)用中心(http://www.cresda.com/[2022-10-22])和對(duì)地觀測數(shù)據(jù)共享計(jì)劃(http://ids.ceode.ac.cn/[2022-10-22])。

    表1 GF-1 WFV數(shù)據(jù)參數(shù)Table 1 GF-1 WFV data parameters

    2.2 MODIS大氣產(chǎn)品

    GF-1 WFV 大氣校正所使用的MODIS AOD 產(chǎn)品包括MOD08D3、MOD09CMA 和MCD19A2。其中MOD08D3提供3種全球逐日1° AOD產(chǎn)品(Wei等,2019),本文采用DTB(DT和DB(Deep Blue,深藍(lán))結(jié)合)算法獲取的AOD;MOD09CMA 提供全球逐日0.05° AOD 產(chǎn)品(Vermote 和Kotchenova,2008);MCD19A2提供全球逐日1 km AOD(Lyapustin等,2018)。大氣水蒸汽含量和臭氧含量分別來自MCD19A2 和MOD09CMG 產(chǎn)品。上述產(chǎn)品均可在MODIS產(chǎn)品網(wǎng)站下載和使用(https://ladsweb.modaps.eosdis.nasa.gov/archive/allData/[2022-10-22])。

    2.3 其他數(shù)據(jù)

    其他數(shù)據(jù)包括全球30 m 分辨率 ASTER GDEM(Abrams 等,2020)和重采樣到2.5 nm 步長的GF-1 WFV光譜響應(yīng)函數(shù)。

    3 數(shù)據(jù)處理算法

    中國陸地區(qū)域2020 年GF-1 WFV 地表反射率產(chǎn)品的生產(chǎn)技術(shù)流程如圖1 所示,主要包括GF-1 WFV 預(yù)處理、多源MODIS AOD 產(chǎn)品空間融合、動(dòng)態(tài)查找表構(gòu)建和快速大氣校正等步驟。

    圖1 GF-1 WFV 地表反射率產(chǎn)品生產(chǎn)技術(shù)流程圖Fig.1 Flow chart of GF-1 WFV surface reflectance product generation

    3.1 GF-1 WFV數(shù)據(jù)預(yù)處理

    在對(duì)GF-1 WFV L1 級(jí)數(shù)據(jù)進(jìn)行正射校正的基礎(chǔ)上,利用輻射定標(biāo)系數(shù)和式(1)將衛(wèi)星影像DN值轉(zhuǎn)換為輻射亮度值,然后利用式(2)計(jì)算表觀反射率TOA(Top of Atmosphere reflectance)。

    式中,Le為輻射亮度,Gain和Bias分別為輻射定標(biāo)系數(shù)增量和偏移量。本文中的輻射定標(biāo)系數(shù)來源于包頭輻射定標(biāo)場獲取的數(shù)據(jù)(Ma等,2020)。

    式中,ρTOA為表觀反射率、Le是輻射亮度、d為日地相對(duì)距離,近似取值為1(劉佳 等,2015)、ESUN為大氣外層太陽輻照度,θ為太陽天頂角。

    3.2 多源MODIS AOD產(chǎn)品空間融合

    多源MODIS AOD 產(chǎn)品的空間融合通常是在以確定某一AOD 產(chǎn)品為主數(shù)據(jù),其他AOD 產(chǎn)品為副數(shù)據(jù)的基礎(chǔ)上,可采用直接替換法、線性替換法或泛克里金法UK(Universal Kriging)進(jìn)行替換。直接替換法針對(duì)主數(shù)據(jù)未覆蓋而副數(shù)據(jù)覆蓋的地區(qū),直接用副數(shù)據(jù)的值進(jìn)行替換。線性替換法針對(duì)主數(shù)據(jù)和副數(shù)據(jù)的相同覆蓋區(qū)域,建立線性回歸模型;對(duì)主數(shù)據(jù)未覆蓋而副數(shù)據(jù)覆蓋的區(qū)域,利用線性回歸模型計(jì)算的值進(jìn)行替換。UK 方法是普通克里金法的改進(jìn)版本,它是一種考慮存在漂移的無偏線性估計(jì)的地統(tǒng)計(jì)方法(Armstrong,1984)。UK 法替換針對(duì)主數(shù)據(jù)和副數(shù)據(jù)的相同覆蓋區(qū)域,在泛克里金的漂移項(xiàng)內(nèi)考慮副數(shù)據(jù)的影響因素,構(gòu)建泛克里金模型。對(duì)于主數(shù)據(jù)未覆蓋而副數(shù)據(jù)覆蓋的地區(qū),利用泛克里金模型計(jì)算的值進(jìn)行替換。Chatterjee 等(2010)曾利用UK方法將MODIS AOD 和MISR AOD 數(shù)據(jù)進(jìn)行融合研究,表明了UK方法在AOD融合上的可行性和有效性。

    為了驗(yàn)證不同融合方法的效果,以MCD19A2多個(gè)AOD 產(chǎn)品的融合為例進(jìn)行實(shí)驗(yàn)。MCD19A2 AOD 產(chǎn)品在一天內(nèi)包含多個(gè)時(shí)間的產(chǎn)品。將產(chǎn)品按上午和下午時(shí)間段劃分并取均值,得到上午數(shù)據(jù)和下午數(shù)據(jù)。以上午數(shù)據(jù)為主數(shù)據(jù),下午數(shù)據(jù)為副數(shù)據(jù),分別采用直接法、線性法和泛克里金法進(jìn)行替換,替換結(jié)果如圖2所示。結(jié)果顯示由于氣溶膠隨時(shí)間變化的影響較大,直接替換會(huì)造成融合數(shù)據(jù)出現(xiàn)不連續(xù)和一致性差的問題(圖2(c))。線性替換由于考慮了兩者的線性關(guān)系,一致性得到了一定的提升,但是融合后的數(shù)據(jù)在局部地區(qū)依然會(huì)出現(xiàn)連續(xù)性和一致性較差的問題(圖2(d))。泛克里金替換通過結(jié)合副數(shù)據(jù)漂移影響和主數(shù)據(jù)插值,融合結(jié)果整體上具有較好的連續(xù)性和一致性(圖2(e))。

    圖2 不同替換方法的效果比較Fig.2 Effects of different replacement methods

    因此本文采用泛克里金替換,將MOD08D3、MOD09CMA 和MCD19A2 的AOD 數(shù)據(jù)進(jìn)行空間融合,獲得三者空間覆蓋度的并集產(chǎn)品。首先將MCD19A2 上午和下午的AOD 數(shù)據(jù)空間融合得到MCD19A2(1 km)。然后將MOD08D3 和MOD09 CMA 不同分辨率的AOD 數(shù)據(jù)空間融合得到MOD0809(5 km)。最后將MCD19A2(1 km)和MOD0809(5 km)的AOD 數(shù)據(jù)再進(jìn)行空間融合,得到MOD080919(1km)的AOD 數(shù)據(jù),具體流程如圖3所示。

    3.3 動(dòng)態(tài)查找表構(gòu)建

    6S(Second Simulation of the Satellite Signal in the Solar Spectrum)模型是由法國大氣光學(xué)實(shí)驗(yàn)室Tanre D.和美國馬里蘭大學(xué)地理系Vermote E.在5S模型基礎(chǔ)上發(fā)展而來(Vermote 等,1997)。6S 模型對(duì)大氣輻射傳輸過程進(jìn)行模擬和計(jì)算,得到用于大氣校正的系數(shù)xa,xb和xc,利用式(3)可計(jì)算得到地表反射率。

    式中,ρa(bǔ)c為地表反射率,Le為輻射亮度,xa、xb和xc為大氣校正系數(shù)。

    查找表可近似替代真實(shí)模型,提高大氣校正的效率。傳統(tǒng)查找表的維度和步長是預(yù)先設(shè)定的,不同的維度和步長會(huì)影響查找表的精度和構(gòu)建時(shí)間。傳統(tǒng)查找表在維度上選擇氣溶膠光學(xué)厚度、高程、大氣水蒸汽含量、太陽天頂角和衛(wèi)星天頂角等(Peng 等,2016)。本文在參考傳統(tǒng)查找表基礎(chǔ)上,提出動(dòng)態(tài)查找表方法。

    通過對(duì)2013 年—2019 年4 萬余景GF-1 WFV的星下點(diǎn)太陽和衛(wèi)星天頂角進(jìn)行統(tǒng)計(jì),發(fā)現(xiàn)其星下點(diǎn)太陽天頂角主要在60°以內(nèi),星下點(diǎn)衛(wèi)星天頂角均在50°以內(nèi)(圖4)。一景GF-1 WFV 影像星下點(diǎn)到影像邊緣的太陽天頂角和衛(wèi)星天頂角相差約1°和8°(王中挺 等,2015)。逐像元天頂角與星下點(diǎn)天頂角反演的地表反射率結(jié)果差異較小,因此以星下點(diǎn)天頂角代表整景影像的天頂角。

    圖4 天頂角統(tǒng)計(jì)結(jié)果Fig.4 Statistics of zenith angles

    根據(jù)上述分析結(jié)果,動(dòng)態(tài)查找表在維度上選取氣溶膠光學(xué)厚度、高程和大氣水蒸汽含量。

    動(dòng)態(tài)查找表的步長則是根據(jù)GF-1 WFV 影像的成像時(shí)間和空間范圍提取對(duì)應(yīng)的氣溶膠光學(xué)厚度、高程和大氣水蒸汽影像的數(shù)值范圍,并分別劃分為8 等分、6 等分和4 等分來動(dòng)態(tài)確定。由于每景GF-1 WFV 影像的成像角度不同,根據(jù)每景影像的具體成像條件,均會(huì)構(gòu)建一個(gè)獨(dú)立的查找表。

    3.4 快速大氣校正

    在上述GF-1 WFV 數(shù)據(jù)預(yù)處理、MODIS 大氣產(chǎn)品預(yù)處理和動(dòng)態(tài)查找表構(gòu)建的基礎(chǔ)上,實(shí)現(xiàn)了GF-1 WFV快速大氣校正,具體流程如下:

    (1)根據(jù)GF-1 WFV 影像提取對(duì)應(yīng)成像時(shí)間和空間范圍的MODIS WV 和AOD 產(chǎn)品及DEM 數(shù)據(jù)。利用泛克里金法將多種不同分辨率的MODIS AOD 產(chǎn)品進(jìn)行空間融合,得到1 km 分辨率全覆蓋的MODIS AOD產(chǎn)品。

    (2)提取對(duì)應(yīng)的MODIS WV、泛克里金融合的AOD 和DEM 數(shù)據(jù)的值域范圍以及GF-1 WFV 影像星下點(diǎn)的天頂角,調(diào)用6S模型構(gòu)建動(dòng)態(tài)查找表。

    (3)綜合考慮大氣校正算法的空間成本和時(shí)間效率,將大氣校正輔助數(shù)據(jù)重采樣到64 m。對(duì)所有數(shù)據(jù)分塊,利用動(dòng)態(tài)查找表生成64 m 分辨率的校正系數(shù)圖像。利用線性插值將校正系數(shù)圖像重采樣到16 m,結(jié)合式(3),計(jì)算地表反射率。

    (4)在多進(jìn)程支持下,并行逐塊執(zhí)行步驟(3),最后將每一塊結(jié)果進(jìn)行拼接得到GF-1 WFV 的地表反射率圖像。

    在Python 開發(fā)環(huán)境中按照上述流程編寫代碼實(shí)現(xiàn)GF-1 WFV快速自動(dòng)化地表反射率反演。

    4 結(jié)果與分析

    利用開發(fā)的GF-1 WFV 地表反射率反演算法生產(chǎn)2020 年中國陸地區(qū)域的GF-1 WFV 地表反射率產(chǎn)品,共計(jì)850景。為了驗(yàn)證產(chǎn)品的可靠性,分別采用大氣校正前后影像對(duì)比、與已有地表反射率產(chǎn)品交叉驗(yàn)證以及實(shí)測數(shù)據(jù)驗(yàn)證進(jìn)行產(chǎn)品質(zhì)量分析和精度驗(yàn)證。

    4.1 GF-1 WFV大氣校正前后影像對(duì)比

    將全國陸地區(qū)域GF-1 WFV 的表觀反射率和地表反射率影像采用相同的拉伸比進(jìn)行對(duì)比,結(jié)果如圖5 所示。從圖5 可以看出,大氣校正后影像的清晰度和質(zhì)量都得到了改善,影像具有更好的輻射一致性。

    通過計(jì)算大氣校正前后的影像直方圖,發(fā)現(xiàn)對(duì)于藍(lán)、綠和紅波段,地表反射率的峰值明顯低于表觀反射率峰值(圖6)。這是由于對(duì)于藍(lán)、綠和紅等較短波長,大氣影響主要表現(xiàn)為氣溶膠散射,經(jīng)過大氣校正,消除了氣溶膠散射的影響,地表反射率直方圖峰值整體“左移”。從“左移”的幅度來看,由于藍(lán)光波長最短,受大氣散射影響最大,“左移”的幅度最大,其次是綠光,紅光波段“左移”的幅度最小。

    圖6 GF-1 WFV大氣校正前后的反射率直方圖對(duì)比Fig.6 Comparison of GF-1 WFV reflectance histograms before and after atmospheric correction

    在GF-1 WFV 影像內(nèi)采集林地、農(nóng)田、海水、湖水、裸土和建筑6種典型地物,比較大氣校正前后的反射率,如圖7所示。經(jīng)過大氣校正,各地物的反射率曲線與標(biāo)準(zhǔn)地物反射率曲線更為一致,例如植被在綠波段有小的反射峰,在紅和藍(lán)波段有吸收谷。同時(shí)典型地物可見光波段反射率的差異隨波長的變短而增大,這也符合較短波長波段更容易受到大氣影響的規(guī)律。

    圖7 GF-1 WFV典型地物大氣校正前后反射率差異Fig.7 Reflectance difference of typical land cover with atmospheric correction of GF-1 WFV

    4.2 GF-1 WFV 與Landsat 8 OLI交叉精度驗(yàn)證

    采用USGS 的Landsat 8 OLI Collection1 Level2地表反射率產(chǎn)品(https://earthexplorer.usgs.gov/[2022-10-22])作為參考,對(duì)GF-1 WFV地表反射率產(chǎn)品開展交叉精度驗(yàn)證。

    選取同步過境且覆蓋同一地理區(qū)域的Landsat 8 OLI 和GF-1 WFV 影像,在影像重疊區(qū)域目視解譯得到驗(yàn)證樣本,樣本類型包含植被、裸土和水體(圖8),得到40 組Landsat 8 OLI和GF-1 WFV 驗(yàn)證數(shù)據(jù)(表2)。

    圖8 交叉驗(yàn)證樣區(qū)樣本空間分布Fig.8 Spatial distribution of validation samples in cross-validation

    表2 交叉驗(yàn)證數(shù)據(jù)對(duì)應(yīng)關(guān)系表Table 2 Cross-validation data correspondence table

    在每組驗(yàn)證數(shù)據(jù)的重疊區(qū)域目視解譯得到約30 個(gè)驗(yàn)證樣區(qū),40 組驗(yàn)證數(shù)據(jù)共得到1206 個(gè)驗(yàn)證樣區(qū),其中裸土樣區(qū)740 個(gè)、植被樣區(qū)293 個(gè)、水體樣區(qū)173個(gè)。計(jì)算每個(gè)驗(yàn)證樣區(qū)的均值來代表整個(gè)驗(yàn)證樣區(qū)。選擇均方根誤差(RMSE)和平均絕對(duì)誤差(MAE)作為精度評(píng)估指標(biāo)。

    式中,n為驗(yàn)證樣區(qū)總數(shù),yi為Landsat 8 OLI 地表反射率,為GF-1 WFV地表反射率。

    考慮到GF-1 WFV 數(shù)據(jù)觀測角度與Landsat 8 OLI 觀測角度存在一定差異,利用C 因子BRDF 校正方法(Roy 等,2016)對(duì)GF-1 WFV 數(shù)據(jù)進(jìn)行BRDF 校正,減小兩種衛(wèi)星數(shù)據(jù)觀測角度差異的影響。為減小兩種衛(wèi)星傳感器光譜響應(yīng)函數(shù)差異的影響,本文基于USGS splib07 的光譜測量數(shù)據(jù)(Kokaly 等,2017),利用光譜匹配因子(SBAF)(Chander 等,2013)將GF1 WFV 的地表反射率轉(zhuǎn)換為等效的Landsat 8 OLI 地表反射率。最終交叉驗(yàn)證結(jié)果如圖9所示。

    圖9 交叉驗(yàn)證結(jié)果Fig.9 Results of cross-validation

    從圖9可以看出,藍(lán)、綠、紅和近紅外波段的平均RMSE分別為1.79%、2.33%、2.90%和2.83%,最大RMSE 不超過3%,且GF-1 WFV 和Landsat 8 OLI 各波段地表反射率的相關(guān)性較好,R2均大于0.9。交叉精度驗(yàn)證結(jié)果表明本文生產(chǎn)的GF-1 WFV地表反射率產(chǎn)品具有較好的可靠性。

    4.3 實(shí)測數(shù)據(jù)驗(yàn)證

    為進(jìn)一步驗(yàn)證GF-1 WFV 地表反射率產(chǎn)品精度,采用ASD 光譜儀實(shí)際觀測的地表反射率數(shù)據(jù)開展直接精度驗(yàn)證。實(shí)測數(shù)據(jù)共有6 組,包括2019年9月13日福建南平的水稻、2019年9月24日甘肅張掖的燕麥、2019 年10 月20 日湖南株洲的水體、2019 年9 月15 日福建平潭的沙灘、2019 年9 月20 日甘肅敦煌的戈壁灘和2019 年9 月22 日甘肅嘉峪關(guān)的戈壁灘。

    首先將所有實(shí)測數(shù)據(jù)轉(zhuǎn)化為GF-1 WFV 等效的地表反射率,然后挑選與實(shí)測數(shù)據(jù)同天或者相差3 天內(nèi)的GF-1 WFV 數(shù)據(jù),并提取實(shí)測數(shù)據(jù)對(duì)應(yīng)的驗(yàn)證區(qū)域,計(jì)算驗(yàn)證區(qū)域內(nèi)藍(lán)、綠、紅和近紅外波段的地表反射率均值,最后與實(shí)測數(shù)據(jù)進(jìn)行對(duì)比驗(yàn)證,驗(yàn)證結(jié)果如表3和圖10所示。

    圖10 光譜曲線對(duì)比Fig.10 Comparison of spectral curves

    表3 實(shí)測數(shù)據(jù)驗(yàn)證結(jié)果Table 3 Validation result of ground-measured data

    從表3 和圖10 可以看出,各實(shí)測點(diǎn)各波段的平均RMSE在1.21%—6.14%,各實(shí)測點(diǎn)的測量光譜曲線和GF-1 WFV的光譜曲線趨勢基本一致。對(duì)于植被(水稻和燕麥),實(shí)測反射率和GF-1 WFV反演反射率在藍(lán)、綠和紅波段較為接近,但近紅外波段的反演反射率低于實(shí)測反射率。對(duì)于水體,兩者在藍(lán)、綠和紅波段的差異也較小,但在近紅外波段反演反射率高于實(shí)測反射率。對(duì)于裸地(沙灘和戈壁灘),GF-1 WFV 藍(lán)、綠、紅和近紅外波段的反演反射率與實(shí)測值均非常接近。直接驗(yàn)證結(jié)果表明本產(chǎn)品具有較好的精度和可靠性。但從驗(yàn)證結(jié)果中也可以發(fā)現(xiàn),相較于藍(lán)、綠和紅波段,植被和水體在近紅外波段的反演誤差較大。這可能是由于清潔水體在近紅外波段的反射率非常低(接近于0),而健康植被在近紅外波段具有較高的反射率(超過50%),但是GF-1 WFV 傳感器在近紅外波段的波譜范圍相對(duì)較寬,導(dǎo)致該波段的光譜探測精度和敏感性降低,對(duì)于較高或較低反射特性地物,反演誤差較大。

    5 產(chǎn)品簡介

    中國陸地區(qū)域2020 年GF-1 WFV 地表反射率產(chǎn)品共有850 景,總計(jì)1.9 TB。數(shù)據(jù)格式包括TIF、XML 和json。產(chǎn)品由地表反射率數(shù)據(jù)(*_AIRAC_MOD080919_SR.TIF)、表觀反射率數(shù)據(jù)(*_TOA.TIF)、元數(shù)據(jù)信息(*.XML)、動(dòng)態(tài)查找表(*.json)和區(qū)域DEM、TO、WV、AOD(*.tif)組成。其中地表反射率和表觀反射率數(shù)據(jù)無量綱,除以10000得到真實(shí)的地表反射率和表觀反射率。產(chǎn)品面向國內(nèi)外免費(fèi)共享,可在ftp://bigrs-info.com/gf1wfv_sr/2020內(nèi)下載和使用。

    6 結(jié)論

    GF-1 WFV 是目前應(yīng)用最為廣泛的國產(chǎn)衛(wèi)星數(shù)據(jù)之一。本文利用多源MODIS AOD 融合產(chǎn)品和動(dòng)態(tài)查找表,采用分塊多進(jìn)程的方法設(shè)計(jì)和實(shí)現(xiàn)GF-1 WFV 快速大氣校正算法,并利用該算法生產(chǎn)和共享中國陸地區(qū)域2020 年GF-1 WFV 地表反射率產(chǎn)品。通過定性分析、Landsat 8 OLI地表反射率產(chǎn)品交叉驗(yàn)證和地面實(shí)際測量的地表反射率數(shù)據(jù)對(duì)產(chǎn)品進(jìn)行精度驗(yàn)證和分析。驗(yàn)證結(jié)果表明,本文生產(chǎn)的GF-1 WFV 地表反射率產(chǎn)品的質(zhì)量良好,與Landsat 8 OLI 地表反射率產(chǎn)品具有較高的一致性,基于實(shí)測地表反射率數(shù)據(jù)的直接精度驗(yàn)證結(jié)果表明各波段的均方根誤差(Root Mean Square Error)位于1.21%—6.14%,產(chǎn)品具有較高的精度和可靠性。

    猜你喜歡
    頂角反射率波段
    春日暖陽
    影響Mini LED板油墨層反射率的因素
    近岸水體異源遙感反射率產(chǎn)品的融合方法研究
    一般三棱鏡最大頂角與折射率的關(guān)系
    具有顏色恒常性的光譜反射率重建
    涼亭中的數(shù)學(xué)
    頂角為100°的等腰三角形性質(zhì)的應(yīng)用
    化學(xué)腐蝕硅表面結(jié)構(gòu)反射率影響因素的研究*
    電子器件(2017年2期)2017-04-25 08:58:37
    M87的多波段輻射過程及其能譜擬合
    日常維護(hù)對(duì)L 波段雷達(dá)的重要性
    西藏科技(2015年4期)2015-09-26 12:12:58
    啦啦啦观看免费观看视频高清| 综合色丁香网| 国产美女午夜福利| 中文字幕人妻熟人妻熟丝袜美| 国产午夜精品论理片| 1000部很黄的大片| 日韩三级伦理在线观看| 久久国内精品自在自线图片| 午夜久久久久精精品| 久久鲁丝午夜福利片| 国产女主播在线喷水免费视频网站 | 小说图片视频综合网站| 搡女人真爽免费视频火全软件| 免费看美女性在线毛片视频| 国产视频内射| 色5月婷婷丁香| 大话2 男鬼变身卡| 国产成人精品婷婷| 能在线免费观看的黄片| 中文字幕av在线有码专区| 免费看日本二区| 日本一二三区视频观看| 国产在视频线在精品| 大香蕉久久网| 国产一区二区三区av在线| 男女边吃奶边做爰视频| 简卡轻食公司| 久久这里有精品视频免费| 小说图片视频综合网站| 黑人高潮一二区| 97热精品久久久久久| 精品99又大又爽又粗少妇毛片| 高清视频免费观看一区二区 | 国产一区二区在线av高清观看| 一区二区三区乱码不卡18| 亚洲图色成人| 全区人妻精品视频| 亚洲国产欧洲综合997久久,| av女优亚洲男人天堂| 久久久a久久爽久久v久久| h日本视频在线播放| 欧美日本亚洲视频在线播放| 大话2 男鬼变身卡| 青春草国产在线视频| 狂野欧美白嫩少妇大欣赏| 最近手机中文字幕大全| 一边摸一边抽搐一进一小说| 久久午夜福利片| 99久久无色码亚洲精品果冻| 插阴视频在线观看视频| 好男人在线观看高清免费视频| 亚洲欧美清纯卡通| 亚洲欧美一区二区三区国产| 在线免费观看不下载黄p国产| 美女高潮的动态| 亚洲精品456在线播放app| 欧美高清性xxxxhd video| 人妻系列 视频| 三级毛片av免费| 日韩欧美国产在线观看| 成人高潮视频无遮挡免费网站| 91aial.com中文字幕在线观看| 国产精品精品国产色婷婷| 久久久成人免费电影| 色播亚洲综合网| 欧美另类亚洲清纯唯美| 最近中文字幕2019免费版| 免费观看a级毛片全部| 听说在线观看完整版免费高清| 国产视频内射| 日韩,欧美,国产一区二区三区 | 国产美女午夜福利| 亚洲在久久综合| 国产真实伦视频高清在线观看| 国产精品女同一区二区软件| 婷婷色av中文字幕| 夫妻性生交免费视频一级片| 国产私拍福利视频在线观看| 亚洲国产精品成人久久小说| 久久亚洲国产成人精品v| 国产三级在线视频| 一个人观看的视频www高清免费观看| 欧美极品一区二区三区四区| 高清午夜精品一区二区三区| 久久久久久大精品| 亚洲精品456在线播放app| 中国美白少妇内射xxxbb| 国产午夜精品久久久久久一区二区三区| 美女xxoo啪啪120秒动态图| 建设人人有责人人尽责人人享有的 | 欧美三级亚洲精品| 久久久久久国产a免费观看| 久99久视频精品免费| 久久久久久久久中文| 久久欧美精品欧美久久欧美| 狠狠狠狠99中文字幕| 久久久久久久国产电影| 美女脱内裤让男人舔精品视频| 日日撸夜夜添| 国内精品宾馆在线| 看黄色毛片网站| 国产精品国产三级专区第一集| 99九九线精品视频在线观看视频| 欧美极品一区二区三区四区| 最近中文字幕2019免费版| av又黄又爽大尺度在线免费看 | 亚洲成色77777| 国产私拍福利视频在线观看| 久99久视频精品免费| 在线播放无遮挡| 黄片无遮挡物在线观看| 久久精品久久精品一区二区三区| 亚洲国产精品国产精品| 午夜福利网站1000一区二区三区| 欧美性猛交╳xxx乱大交人| 又爽又黄无遮挡网站| 岛国毛片在线播放| 亚洲av.av天堂| 国产老妇伦熟女老妇高清| 欧美潮喷喷水| 91aial.com中文字幕在线观看| 一二三四中文在线观看免费高清| 99久久成人亚洲精品观看| 久久精品91蜜桃| 又粗又硬又长又爽又黄的视频| 国产人妻一区二区三区在| 国产精品不卡视频一区二区| 久久久成人免费电影| 亚洲av免费高清在线观看| 国产淫语在线视频| 亚洲欧美精品专区久久| 日韩高清综合在线| 欧美日本视频| 婷婷色综合大香蕉| 一本一本综合久久| 欧美激情国产日韩精品一区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品一区二区三区四区久久| 水蜜桃什么品种好| 国产毛片a区久久久久| 精品久久久久久电影网 | 亚洲精品乱码久久久久久按摩| 嫩草影院精品99| 嘟嘟电影网在线观看| 网址你懂的国产日韩在线| 亚洲精品色激情综合| 中文字幕精品亚洲无线码一区| 国产精品99久久久久久久久| 亚洲精品色激情综合| 精品无人区乱码1区二区| 欧美激情国产日韩精品一区| 亚洲av成人av| 丝袜喷水一区| 婷婷色麻豆天堂久久 | 中文亚洲av片在线观看爽| 精品午夜福利在线看| 狂野欧美白嫩少妇大欣赏| 久久鲁丝午夜福利片| 国产精品,欧美在线| 亚洲成人久久爱视频| 国产伦理片在线播放av一区| 国产私拍福利视频在线观看| 久久久国产成人精品二区| 国产精品1区2区在线观看.| 国产欧美另类精品又又久久亚洲欧美| 欧美丝袜亚洲另类| 全区人妻精品视频| 舔av片在线| 亚洲av电影不卡..在线观看| 欧美精品一区二区大全| 日韩成人av中文字幕在线观看| 别揉我奶头 嗯啊视频| 永久免费av网站大全| 在线观看美女被高潮喷水网站| 如何舔出高潮| 永久网站在线| 亚洲aⅴ乱码一区二区在线播放| 久久人妻av系列| 久久精品人妻少妇| 97热精品久久久久久| 精品久久久久久久人妻蜜臀av| 精品人妻视频免费看| 国产精品久久久久久av不卡| 久久久久久久久中文| 午夜a级毛片| 99久国产av精品国产电影| 免费av观看视频| 婷婷六月久久综合丁香| 欧美丝袜亚洲另类| 国产精品久久久久久久电影| 午夜亚洲福利在线播放| 乱码一卡2卡4卡精品| 1000部很黄的大片| 美女内射精品一级片tv| 久久久久免费精品人妻一区二区| 久久精品影院6| 国产乱来视频区| 国产精品1区2区在线观看.| 亚洲欧洲日产国产| 久久草成人影院| 国产熟女欧美一区二区| 美女xxoo啪啪120秒动态图| 午夜久久久久精精品| 国产极品天堂在线| 春色校园在线视频观看| 久久久久免费精品人妻一区二区| 美女国产视频在线观看| 国产单亲对白刺激| 男人舔奶头视频| 国产精品久久久久久av不卡| 欧美+日韩+精品| av免费在线看不卡| 欧美极品一区二区三区四区| 成人毛片a级毛片在线播放| 精品免费久久久久久久清纯| 午夜视频国产福利| 亚洲一级一片aⅴ在线观看| 秋霞在线观看毛片| 国产女主播在线喷水免费视频网站 | 一边摸一边抽搐一进一小说| 国产探花极品一区二区| 久久久久久久午夜电影| 久久久久久久亚洲中文字幕| 在线观看一区二区三区| 人妻系列 视频| 国产av码专区亚洲av| 久久久精品欧美日韩精品| 国产亚洲av片在线观看秒播厂 | 国产成人91sexporn| 热99re8久久精品国产| 亚洲成人av在线免费| 国产视频首页在线观看| 啦啦啦韩国在线观看视频| 少妇被粗大猛烈的视频| 久久99热6这里只有精品| 亚洲电影在线观看av| 美女xxoo啪啪120秒动态图| 国产成人a区在线观看| 亚洲乱码一区二区免费版| 男插女下体视频免费在线播放| 日韩成人伦理影院| 国产精品三级大全| 丝袜喷水一区| 国产成人91sexporn| 欧美3d第一页| 久久久亚洲精品成人影院| 国产精品久久久久久av不卡| 亚洲欧美日韩东京热| 亚洲av成人av| 日本一本二区三区精品| 成人国产麻豆网| 亚洲成人久久爱视频| 成人亚洲欧美一区二区av| 免费看美女性在线毛片视频| 久久精品熟女亚洲av麻豆精品 | 在线免费观看的www视频| 又爽又黄无遮挡网站| 永久免费av网站大全| 嫩草影院入口| 久久鲁丝午夜福利片| 国产精品日韩av在线免费观看| 在线播放国产精品三级| 久久精品夜夜夜夜夜久久蜜豆| 久99久视频精品免费| 最近中文字幕2019免费版| 搡老妇女老女人老熟妇| 中文资源天堂在线| 高清av免费在线| 午夜精品国产一区二区电影 | 日日摸夜夜添夜夜添av毛片| 欧美一区二区国产精品久久精品| 免费观看精品视频网站| 精品午夜福利在线看| 高清毛片免费看| 国产av在哪里看| 好男人在线观看高清免费视频| 久久久久久久久久久免费av| 国产高清有码在线观看视频| 看片在线看免费视频| 国产成人91sexporn| 国产精品伦人一区二区| 亚洲精品久久久久久婷婷小说 | 在线a可以看的网站| 一级毛片久久久久久久久女| 中国美白少妇内射xxxbb| АⅤ资源中文在线天堂| 国产又黄又爽又无遮挡在线| 国产伦精品一区二区三区视频9| 亚洲一区高清亚洲精品| 蜜臀久久99精品久久宅男| 国内精品宾馆在线| 成年女人看的毛片在线观看| 欧美成人一区二区免费高清观看| 黄色一级大片看看| 成人午夜高清在线视频| 久久精品国产亚洲av涩爱| 男女边吃奶边做爰视频| 亚洲av成人精品一区久久| 亚洲av成人精品一二三区| 超碰97精品在线观看| 亚洲精品日韩在线中文字幕| 51国产日韩欧美| 亚洲国产高清在线一区二区三| 91精品伊人久久大香线蕉| 亚洲国产日韩欧美精品在线观看| 中文字幕熟女人妻在线| 成人一区二区视频在线观看| 色综合亚洲欧美另类图片| 国产精品蜜桃在线观看| 精品人妻视频免费看| 免费在线观看成人毛片| 国产成人精品久久久久久| 丰满少妇做爰视频| 99视频精品全部免费 在线| 国产黄片视频在线免费观看| 日韩国内少妇激情av| 最近视频中文字幕2019在线8| 成人综合一区亚洲| 亚洲国产欧洲综合997久久,| 黄色欧美视频在线观看| 国产成年人精品一区二区| 精品久久久久久久久久久久久| 日本色播在线视频| 女的被弄到高潮叫床怎么办| 免费观看a级毛片全部| 国产亚洲一区二区精品| 只有这里有精品99| 三级国产精品片| 男的添女的下面高潮视频| 日韩高清综合在线| 在线a可以看的网站| 免费搜索国产男女视频| 少妇高潮的动态图| 午夜亚洲福利在线播放| 亚洲欧美日韩东京热| 久久综合国产亚洲精品| 亚洲在久久综合| 美女xxoo啪啪120秒动态图| 18+在线观看网站| 中文精品一卡2卡3卡4更新| videos熟女内射| 亚洲怡红院男人天堂| 女的被弄到高潮叫床怎么办| 久久久久精品久久久久真实原创| 一卡2卡三卡四卡精品乱码亚洲| 国产精品福利在线免费观看| 成人二区视频| 亚洲成人中文字幕在线播放| 如何舔出高潮| 免费播放大片免费观看视频在线观看 | 在线播放无遮挡| 色综合站精品国产| 亚洲国产日韩欧美精品在线观看| 青青草视频在线视频观看| 国产又色又爽无遮挡免| 狂野欧美激情性xxxx在线观看| 国产免费又黄又爽又色| 有码 亚洲区| 黄色欧美视频在线观看| 国产一区二区在线av高清观看| 国产探花在线观看一区二区| 国国产精品蜜臀av免费| 国产亚洲精品av在线| or卡值多少钱| 久久亚洲国产成人精品v| 激情 狠狠 欧美| 国产亚洲精品av在线| 久久久精品欧美日韩精品| 97超碰精品成人国产| 久久久久久大精品| 麻豆成人av视频| av在线天堂中文字幕| 天天一区二区日本电影三级| 91在线精品国自产拍蜜月| 精品少妇黑人巨大在线播放 | 九九在线视频观看精品| 免费观看的影片在线观看| 99视频精品全部免费 在线| 一区二区三区乱码不卡18| 国产av码专区亚洲av| 美女cb高潮喷水在线观看| 中文亚洲av片在线观看爽| 午夜免费激情av| 欧美三级亚洲精品| 国产人妻一区二区三区在| 日本免费一区二区三区高清不卡| 狂野欧美白嫩少妇大欣赏| 亚洲激情五月婷婷啪啪| 人妻系列 视频| 汤姆久久久久久久影院中文字幕 | 亚洲精品一区蜜桃| 高清日韩中文字幕在线| 免费观看精品视频网站| 日韩强制内射视频| 麻豆成人av视频| 91在线精品国自产拍蜜月| 色综合色国产| 亚洲va在线va天堂va国产| 国产老妇女一区| 亚洲婷婷狠狠爱综合网| 夜夜看夜夜爽夜夜摸| 在现免费观看毛片| 国产精品综合久久久久久久免费| 日韩一区二区三区影片| 亚洲自拍偷在线| 久久精品国产99精品国产亚洲性色| 亚洲自拍偷在线| 亚洲国产成人一精品久久久| 亚洲av成人精品一二三区| 精品人妻一区二区三区麻豆| 免费不卡的大黄色大毛片视频在线观看 | 少妇的逼水好多| 男女边吃奶边做爰视频| 亚洲国产日韩欧美精品在线观看| 五月玫瑰六月丁香| 亚洲天堂国产精品一区在线| 91av网一区二区| 一区二区三区高清视频在线| 日韩欧美三级三区| 国产真实乱freesex| 少妇熟女aⅴ在线视频| av免费观看日本| 精品无人区乱码1区二区| 国产精品综合久久久久久久免费| 亚洲综合精品二区| 国产色婷婷99| 国产伦精品一区二区三区视频9| 超碰97精品在线观看| 精品国产露脸久久av麻豆 | 人人妻人人澡人人爽人人夜夜 | 高清午夜精品一区二区三区| 成人漫画全彩无遮挡| 国产成人精品一,二区| 免费无遮挡裸体视频| 亚洲欧洲日产国产| 亚洲在线自拍视频| av在线蜜桃| av在线天堂中文字幕| 国产欧美日韩精品一区二区| 久久婷婷人人爽人人干人人爱| 毛片一级片免费看久久久久| 色5月婷婷丁香| 国产精品一及| 国产麻豆成人av免费视频| 大话2 男鬼变身卡| 国模一区二区三区四区视频| 欧美激情久久久久久爽电影| 国产国拍精品亚洲av在线观看| 国产乱人偷精品视频| 免费播放大片免费观看视频在线观看 | 一边摸一边抽搐一进一小说| 日本熟妇午夜| 亚洲国产欧洲综合997久久,| 久久精品国产鲁丝片午夜精品| 久久久久久九九精品二区国产| 国产高清有码在线观看视频| 啦啦啦韩国在线观看视频| 国产成人精品婷婷| 亚洲av成人精品一区久久| 热99在线观看视频| 国产乱来视频区| 91久久精品电影网| 国产亚洲av嫩草精品影院| 一个人看的www免费观看视频| 日日干狠狠操夜夜爽| 午夜视频国产福利| 欧美精品国产亚洲| videos熟女内射| 男女下面进入的视频免费午夜| 精品无人区乱码1区二区| 久久精品人妻少妇| 国产精品一区www在线观看| 天堂中文最新版在线下载 | 亚洲国产日韩欧美精品在线观看| 中文精品一卡2卡3卡4更新| 尤物成人国产欧美一区二区三区| 3wmmmm亚洲av在线观看| 亚洲精品456在线播放app| 亚洲av中文av极速乱| 国产中年淑女户外野战色| 2022亚洲国产成人精品| 一级黄片播放器| 国产免费福利视频在线观看| 亚洲婷婷狠狠爱综合网| 欧美激情久久久久久爽电影| 男女国产视频网站| 秋霞伦理黄片| 长腿黑丝高跟| 国产又黄又爽又无遮挡在线| 日本wwww免费看| 亚洲一区高清亚洲精品| 婷婷色麻豆天堂久久 | 亚洲国产高清在线一区二区三| 日本一二三区视频观看| 亚洲中文字幕一区二区三区有码在线看| 久久久国产成人免费| 22中文网久久字幕| 麻豆精品久久久久久蜜桃| 亚洲成人中文字幕在线播放| 欧美成人午夜免费资源| 在线免费观看的www视频| 久久精品熟女亚洲av麻豆精品 | 亚洲欧洲日产国产| 中文资源天堂在线| 岛国在线免费视频观看| 伦理电影大哥的女人| 一本久久精品| 男的添女的下面高潮视频| 三级国产精品片| 看非洲黑人一级黄片| 老师上课跳d突然被开到最大视频| 免费不卡的大黄色大毛片视频在线观看 | 欧美色视频一区免费| 熟女电影av网| 一级毛片我不卡| 一级黄色大片毛片| av卡一久久| 一级黄色大片毛片| 插阴视频在线观看视频| 国产老妇女一区| 国产成人a∨麻豆精品| 久久亚洲国产成人精品v| 亚洲最大成人av| 国产探花在线观看一区二区| 国产成人精品婷婷| 免费人成在线观看视频色| 久久久久久国产a免费观看| 人人妻人人看人人澡| 欧美性猛交黑人性爽| 中文字幕av在线有码专区| 在线观看美女被高潮喷水网站| 中文字幕久久专区| 国产成人精品婷婷| 天堂网av新在线| 久久久久国产网址| 观看美女的网站| 久久久久国产网址| 天美传媒精品一区二区| 搞女人的毛片| 18禁在线无遮挡免费观看视频| 欧美三级亚洲精品| 亚洲欧美中文字幕日韩二区| 一区二区三区高清视频在线| 国产真实乱freesex| 亚洲精品乱码久久久久久按摩| 男人和女人高潮做爰伦理| 人体艺术视频欧美日本| 午夜福利在线观看吧| 国产av码专区亚洲av| 精品欧美国产一区二区三| 一级黄色大片毛片| 国内精品宾馆在线| 国产亚洲5aaaaa淫片| 99热这里只有精品一区| 村上凉子中文字幕在线| 超碰97精品在线观看| 亚洲欧美精品专区久久| 99久国产av精品| 免费观看a级毛片全部| 99热这里只有精品一区| 国产精品av视频在线免费观看| 国产精品永久免费网站| 欧美精品国产亚洲| 国产成人一区二区在线| 性插视频无遮挡在线免费观看| 日韩,欧美,国产一区二区三区 | 免费看光身美女| 丝袜美腿在线中文| www.av在线官网国产| 六月丁香七月| 永久免费av网站大全| 精品无人区乱码1区二区| 亚洲欧美精品专区久久| 国产精品国产三级国产专区5o | 亚洲精品456在线播放app| 久久精品国产自在天天线| 青春草国产在线视频| 免费观看a级毛片全部| 欧美成人午夜免费资源| 我的女老师完整版在线观看| 国产av一区在线观看免费| 国产精品,欧美在线| 老师上课跳d突然被开到最大视频| 国产一区亚洲一区在线观看| 亚洲欧美成人综合另类久久久 | 国产极品天堂在线| 婷婷色av中文字幕| 午夜a级毛片| 最近最新中文字幕大全电影3| 国产国拍精品亚洲av在线观看| 欧美变态另类bdsm刘玥| 久久精品久久精品一区二区三区| 色网站视频免费| 国产淫片久久久久久久久| 人人妻人人澡欧美一区二区| 蜜桃亚洲精品一区二区三区| 看非洲黑人一级黄片| 国产高清三级在线| a级毛片免费高清观看在线播放| 天天一区二区日本电影三级| 午夜老司机福利剧场| .国产精品久久| 精品欧美国产一区二区三| 蜜臀久久99精品久久宅男| 在线播放国产精品三级| 久久久久久久午夜电影| 欧美不卡视频在线免费观看| 日本一二三区视频观看| 婷婷色麻豆天堂久久 | 亚洲欧洲国产日韩| 丝袜喷水一区| 亚洲av电影在线观看一区二区三区 | 99久久中文字幕三级久久日本| 国产综合懂色| 一级毛片久久久久久久久女|