• <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
    国产一区二区三区综合在线观看 | 男人操女人黄网站| 精品一品国产午夜福利视频| 亚洲美女搞黄在线观看| 51国产日韩欧美| 美女内射精品一级片tv| 人人澡人人妻人| 亚洲国产欧美在线一区| 亚洲精品色激情综合| 日本与韩国留学比较| 在线天堂中文资源库| 考比视频在线观看| 一个人免费看片子| 乱码一卡2卡4卡精品| 欧美亚洲日本最大视频资源| 丁香六月天网| 精品一区二区三区视频在线| 日韩免费高清中文字幕av| 中文字幕av电影在线播放| 久久久久人妻精品一区果冻| 亚洲av电影在线观看一区二区三区| 国产精品秋霞免费鲁丝片| 日韩av在线免费看完整版不卡| 满18在线观看网站| videosex国产| 久久久久精品人妻al黑| 午夜福利乱码中文字幕| 国产色婷婷99| 又黄又爽又刺激的免费视频.| 中文字幕亚洲精品专区| 丝袜脚勾引网站| 精品国产国语对白av| 免费观看性生交大片5| 人人澡人人妻人| 久久影院123| 97在线视频观看| 国产成人精品婷婷| 久久久久精品久久久久真实原创| 免费高清在线观看视频在线观看| 在线看a的网站| 一级毛片黄色毛片免费观看视频| 国产精品.久久久| 免费大片黄手机在线观看| 免费看不卡的av| 午夜福利乱码中文字幕| 亚洲三级黄色毛片| 男女午夜视频在线观看 | 精品人妻熟女毛片av久久网站| 亚洲丝袜综合中文字幕| 涩涩av久久男人的天堂| 自线自在国产av| 中文字幕制服av| a级毛色黄片| 一二三四在线观看免费中文在 | 日日啪夜夜爽| 国产无遮挡羞羞视频在线观看| 亚洲欧美中文字幕日韩二区| 亚洲欧美清纯卡通| 啦啦啦视频在线资源免费观看| 天堂8中文在线网| 国产一区亚洲一区在线观看| 国产成人免费观看mmmm| 22中文网久久字幕| 精品福利永久在线观看| 国产av精品麻豆| 国产白丝娇喘喷水9色精品| 国产精品不卡视频一区二区| 久久狼人影院| 久久这里有精品视频免费| 丝瓜视频免费看黄片| 啦啦啦啦在线视频资源| 97人妻天天添夜夜摸| 女性被躁到高潮视频| 亚洲精品国产av成人精品| 欧美日韩国产mv在线观看视频| 亚洲国产精品一区三区| 欧美性感艳星| 久久影院123| 国产精品女同一区二区软件| 激情五月婷婷亚洲| 自拍欧美九色日韩亚洲蝌蚪91| 成人国产av品久久久| 亚洲精品自拍成人| 香蕉丝袜av| 全区人妻精品视频| 你懂的网址亚洲精品在线观看| 99视频精品全部免费 在线| 国产精品一区二区在线观看99| 亚洲成av片中文字幕在线观看 | 亚洲四区av| 大香蕉97超碰在线| 下体分泌物呈黄色| 欧美丝袜亚洲另类| 在线亚洲精品国产二区图片欧美| 日韩制服丝袜自拍偷拍| 女的被弄到高潮叫床怎么办| 美女视频免费永久观看网站| 晚上一个人看的免费电影| 欧美 日韩 精品 国产| 如何舔出高潮| 日韩制服丝袜自拍偷拍| 新久久久久国产一级毛片| 亚洲国产欧美日韩在线播放| 99久久人妻综合| 乱码一卡2卡4卡精品| 十八禁高潮呻吟视频| 伊人久久国产一区二区| 亚洲精品自拍成人| 亚洲av中文av极速乱| 亚洲美女黄色视频免费看| 看非洲黑人一级黄片| 一级毛片 在线播放| 亚洲av电影在线观看一区二区三区| 18禁裸乳无遮挡动漫免费视频| 极品人妻少妇av视频| 亚洲国产精品专区欧美| 女人精品久久久久毛片| 大香蕉久久成人网| 国产毛片在线视频| 少妇被粗大猛烈的视频| 亚洲av电影在线观看一区二区三区| 色视频在线一区二区三区| 中文天堂在线官网| 亚洲一码二码三码区别大吗| 男的添女的下面高潮视频| 18在线观看网站| 在线观看免费日韩欧美大片| 精品国产乱码久久久久久小说| 一区二区av电影网| 99久久中文字幕三级久久日本| 桃花免费在线播放| 国产国拍精品亚洲av在线观看| 欧美国产精品一级二级三级| 91在线精品国自产拍蜜月| 在线观看免费高清a一片| 亚洲国产最新在线播放| 色婷婷久久久亚洲欧美| 51国产日韩欧美| 亚洲成人一二三区av| 久久精品国产a三级三级三级| 久久久a久久爽久久v久久| 日韩电影二区| 亚洲性久久影院| 免费在线观看完整版高清| 免费大片黄手机在线观看| 欧美日韩一区二区视频在线观看视频在线| 国产有黄有色有爽视频| 免费观看在线日韩| 日本黄大片高清| 只有这里有精品99| 草草在线视频免费看| 黑人高潮一二区| 亚洲国产精品一区三区| 中文字幕av电影在线播放| 一本大道久久a久久精品| 五月天丁香电影| 最近的中文字幕免费完整| 男人舔女人的私密视频| 精品视频人人做人人爽| 亚洲,欧美,日韩| 国产成人精品福利久久| 国产免费一级a男人的天堂| 日本av免费视频播放| 日韩三级伦理在线观看| 青青草视频在线视频观看| 男人爽女人下面视频在线观看| 九九爱精品视频在线观看| 巨乳人妻的诱惑在线观看| 久久久久久久久久久免费av| 久久99一区二区三区| 黄片无遮挡物在线观看| 久久精品国产自在天天线| 精品久久蜜臀av无| 久久青草综合色| 一二三四在线观看免费中文在 | 热99久久久久精品小说推荐| 亚洲中文av在线| 国产永久视频网站| 欧美日韩视频精品一区| 亚洲综合色网址| 一二三四中文在线观看免费高清| av不卡在线播放| 天堂俺去俺来也www色官网| 乱人伦中国视频| 国产又色又爽无遮挡免| 性高湖久久久久久久久免费观看| 交换朋友夫妻互换小说| 国产熟女午夜一区二区三区| 欧美少妇被猛烈插入视频| 久久韩国三级中文字幕| 久久精品夜色国产| 亚洲少妇的诱惑av| 人妻一区二区av| 最近最新中文字幕大全免费视频 | 国产一区二区三区av在线| 国产激情久久老熟女| 国产免费一区二区三区四区乱码| 国产乱来视频区| 涩涩av久久男人的天堂| 欧美+日韩+精品| 国产成人精品久久久久久| 最黄视频免费看| 亚洲精品国产色婷婷电影| 波多野结衣一区麻豆| 成人手机av| 亚洲av.av天堂| av天堂久久9| 99久久精品国产国产毛片| 久久99热6这里只有精品| 天天影视国产精品| 捣出白浆h1v1| 精品福利永久在线观看| 亚洲av.av天堂| 人人澡人人妻人| av女优亚洲男人天堂| 欧美成人精品欧美一级黄| 亚洲国产毛片av蜜桃av| 国产精品人妻久久久久久| 99国产综合亚洲精品| 国产熟女午夜一区二区三区| 精品一区在线观看国产| av不卡在线播放| 丝袜在线中文字幕| 精品视频人人做人人爽| 国产精品一区二区在线观看99| 国产免费一级a男人的天堂| 一区在线观看完整版| 大片电影免费在线观看免费| 亚洲精品av麻豆狂野| 五月天丁香电影| av在线app专区| 亚洲性久久影院| 久久99一区二区三区| 亚洲天堂av无毛| 夫妻午夜视频| 久久韩国三级中文字幕| 国产精品国产三级国产专区5o| 国产高清三级在线| 午夜福利,免费看| av网站免费在线观看视频| 亚洲综合色惰| 亚洲精品av麻豆狂野| 日韩 亚洲 欧美在线| 两个人免费观看高清视频| 亚洲欧美精品自产自拍| 日韩中字成人| 少妇的逼好多水| 在线观看免费高清a一片| 国产精品一区www在线观看| 乱人伦中国视频| 午夜福利,免费看| 欧美精品一区二区大全| 亚洲国产精品专区欧美| 亚洲欧洲精品一区二区精品久久久 | 亚洲一码二码三码区别大吗| 五月伊人婷婷丁香| 免费久久久久久久精品成人欧美视频 | 亚洲成国产人片在线观看| 黄片播放在线免费| 熟女人妻精品中文字幕| 伊人久久国产一区二区| 日韩,欧美,国产一区二区三区| 国产精品 国内视频| 99久久中文字幕三级久久日本| 国产成人免费无遮挡视频| 天堂俺去俺来也www色官网| 成年人免费黄色播放视频| 青春草视频在线免费观看| 伦理电影大哥的女人| 一本—道久久a久久精品蜜桃钙片| 久久人人爽人人片av| 久久精品国产亚洲av涩爱| 国产黄频视频在线观看| 久久韩国三级中文字幕| 亚洲精品久久成人aⅴ小说| 国产 一区精品| 毛片一级片免费看久久久久| videos熟女内射| 亚洲精品视频女| 久久久久人妻精品一区果冻| 久热久热在线精品观看| 99国产综合亚洲精品| 成人国产麻豆网| 人人妻人人爽人人添夜夜欢视频| 亚洲国产精品一区二区三区在线| 国产麻豆69| av国产精品久久久久影院| 少妇人妻精品综合一区二区| 久久精品熟女亚洲av麻豆精品| 国产在线免费精品| 成人黄色视频免费在线看| 精品第一国产精品| 国产老妇伦熟女老妇高清| av在线app专区| 亚洲国产日韩一区二区| 免费人妻精品一区二区三区视频| 国产亚洲午夜精品一区二区久久| 一级毛片我不卡| 韩国精品一区二区三区 | 日本午夜av视频| 国产成人av激情在线播放| 九色亚洲精品在线播放| 青春草视频在线免费观看| 天堂中文最新版在线下载| 亚洲成国产人片在线观看| 一区二区三区乱码不卡18| 男女午夜视频在线观看 | av在线播放精品| 久久精品久久久久久久性| 国产伦理片在线播放av一区| 国产精品国产三级专区第一集| 成人影院久久| 最近的中文字幕免费完整| 80岁老熟妇乱子伦牲交| 免费女性裸体啪啪无遮挡网站| 久久这里只有精品19| 美女内射精品一级片tv| 纯流量卡能插随身wifi吗| 欧美人与善性xxx| 日产精品乱码卡一卡2卡三| 国产片内射在线| 日日啪夜夜爽| 国产白丝娇喘喷水9色精品| 美女主播在线视频| 日本爱情动作片www.在线观看| 欧美3d第一页| 人人妻人人澡人人爽人人夜夜| 人人妻人人澡人人看| 性色avwww在线观看| 乱人伦中国视频| 色94色欧美一区二区| 少妇人妻 视频| 国产成人免费无遮挡视频| 高清毛片免费看| 少妇高潮的动态图| 亚洲欧美成人综合另类久久久| 亚洲 欧美一区二区三区| 午夜精品国产一区二区电影| 黄色一级大片看看| 日韩精品有码人妻一区| 久久97久久精品| 欧美日韩亚洲高清精品| 久久精品国产亚洲av涩爱| 女人久久www免费人成看片| 国产精品无大码| 一本色道久久久久久精品综合| 9191精品国产免费久久| 国产高清三级在线| 男女下面插进去视频免费观看 | 女的被弄到高潮叫床怎么办| 九草在线视频观看| 亚洲国产av新网站| 丰满少妇做爰视频| 男女无遮挡免费网站观看| 日产精品乱码卡一卡2卡三| 校园人妻丝袜中文字幕| 日韩制服丝袜自拍偷拍| 99精国产麻豆久久婷婷| av又黄又爽大尺度在线免费看| 97人妻天天添夜夜摸| 免费看不卡的av| 91国产中文字幕| 成人漫画全彩无遮挡| 晚上一个人看的免费电影| 色网站视频免费| 久久人人爽人人爽人人片va| 熟女人妻精品中文字幕| a级片在线免费高清观看视频| 亚洲色图 男人天堂 中文字幕 | 啦啦啦啦在线视频资源| 女性生殖器流出的白浆| 日韩不卡一区二区三区视频在线| 最后的刺客免费高清国语| 成年女人在线观看亚洲视频| 免费高清在线观看视频在线观看| 综合色丁香网| 久久精品国产鲁丝片午夜精品| 满18在线观看网站| av一本久久久久| 如何舔出高潮| 日本免费在线观看一区| 又粗又硬又长又爽又黄的视频| 黄色视频在线播放观看不卡| 亚洲国产精品999| 久久精品久久久久久噜噜老黄| 日韩欧美精品免费久久| 久久婷婷青草| 69精品国产乱码久久久| 亚洲成人手机| 国产色婷婷99| 日韩av在线免费看完整版不卡| 久久韩国三级中文字幕| 男人添女人高潮全过程视频| 国产极品天堂在线| 中文字幕av电影在线播放| 亚洲精品一区蜜桃| 最近最新中文字幕大全免费视频 | 色哟哟·www| 亚洲av综合色区一区| 亚洲经典国产精华液单| 免费久久久久久久精品成人欧美视频 | xxxhd国产人妻xxx| 男人舔女人的私密视频| 精品人妻偷拍中文字幕| 欧美性感艳星| 亚洲少妇的诱惑av| 亚洲欧美一区二区三区黑人 | 女人精品久久久久毛片| 一本久久精品| av视频免费观看在线观看| 免费人妻精品一区二区三区视频| 欧美性感艳星| 精品亚洲乱码少妇综合久久| 在线看a的网站| 久久热在线av| 国产成人欧美| 插逼视频在线观看| 亚洲内射少妇av| 天堂俺去俺来也www色官网| 大香蕉久久网| 免费不卡的大黄色大毛片视频在线观看| 熟女人妻精品中文字幕| 亚洲,欧美精品.| 菩萨蛮人人尽说江南好唐韦庄| 日本欧美国产在线视频| 亚洲国产精品999| 国产欧美日韩综合在线一区二区| 99久国产av精品国产电影| 欧美人与善性xxx| 99久久中文字幕三级久久日本| 水蜜桃什么品种好| 性高湖久久久久久久久免费观看| 草草在线视频免费看| 9热在线视频观看99| 免费av不卡在线播放| 亚洲综合色网址| 超碰97精品在线观看| 国产有黄有色有爽视频| 国产亚洲最大av| 欧美日韩亚洲高清精品| 永久网站在线| 国产综合精华液| av线在线观看网站| 日本黄大片高清| 观看美女的网站| 国产精品久久久av美女十八| 国产精品 国内视频| 91精品伊人久久大香线蕉| 久久久国产欧美日韩av| 香蕉国产在线看| av女优亚洲男人天堂| 久久久久国产网址| 国产男女内射视频| 黄网站色视频无遮挡免费观看| 国产av码专区亚洲av| 在线观看三级黄色| 制服诱惑二区| 精品久久久精品久久久| 欧美激情 高清一区二区三区| av免费在线看不卡| 韩国精品一区二区三区 | 草草在线视频免费看| 熟女人妻精品中文字幕| 国产精品一区二区在线观看99| 国产激情久久老熟女| 在线亚洲精品国产二区图片欧美| 边亲边吃奶的免费视频| 国产精品久久久久久av不卡| 少妇的逼水好多| 国产成人免费观看mmmm| 在线观看www视频免费| 九九在线视频观看精品| 久久99热这里只频精品6学生| 精品视频人人做人人爽| 国产极品粉嫩免费观看在线| 51国产日韩欧美| 成人亚洲欧美一区二区av| 看免费av毛片| 成人午夜精彩视频在线观看| 一区二区av电影网| 少妇熟女欧美另类| 制服丝袜香蕉在线| xxxhd国产人妻xxx| 午夜福利视频在线观看免费| 国产免费一级a男人的天堂| 高清av免费在线| 欧美日韩av久久| 国产不卡av网站在线观看| 国产精品秋霞免费鲁丝片| 欧美人与性动交α欧美软件 | 香蕉国产在线看| 亚洲国产欧美在线一区| 欧美国产精品va在线观看不卡| 久久精品aⅴ一区二区三区四区 | 亚洲精品美女久久久久99蜜臀 | 日本av手机在线免费观看| 久久久精品免费免费高清| 日日撸夜夜添| 久久久久久久国产电影| 欧美精品一区二区免费开放| 欧美 亚洲 国产 日韩一| 国产男女超爽视频在线观看| av国产精品久久久久影院| 色5月婷婷丁香| 免费av中文字幕在线| 久久久久久久亚洲中文字幕| 草草在线视频免费看| 婷婷成人精品国产| 亚洲av国产av综合av卡| 99久久精品国产国产毛片| 看十八女毛片水多多多| 久久97久久精品| 国产毛片在线视频| 亚洲精品国产色婷婷电影| 青春草视频在线免费观看| 一级毛片 在线播放| 久久人人97超碰香蕉20202| 国产 一区精品| 人妻系列 视频| 亚洲一区二区三区欧美精品| 国产成人aa在线观看| 满18在线观看网站| 午夜福利影视在线免费观看| 波野结衣二区三区在线| 日日爽夜夜爽网站| 亚洲欧美成人精品一区二区| 多毛熟女@视频| 久久婷婷青草| 亚洲av福利一区| 精品人妻一区二区三区麻豆| 日韩中字成人| 另类精品久久| 婷婷色麻豆天堂久久| 国产精品99久久99久久久不卡 | 国内精品宾馆在线| 欧美日韩综合久久久久久| 在线观看免费高清a一片| 亚洲中文av在线| 免费观看无遮挡的男女| 99国产精品免费福利视频| 人人澡人人妻人| 18+在线观看网站| av一本久久久久| 日韩一区二区视频免费看| 日本猛色少妇xxxxx猛交久久| 欧美97在线视频| 爱豆传媒免费全集在线观看| 精品人妻一区二区三区麻豆| 久久这里只有精品19| 国产深夜福利视频在线观看| 免费人妻精品一区二区三区视频| 老司机影院成人| 90打野战视频偷拍视频| 国产av国产精品国产| 99国产精品免费福利视频| 搡老乐熟女国产| 热99久久久久精品小说推荐| 久久久精品94久久精品| 亚洲国产看品久久| 大片电影免费在线观看免费| 丝袜喷水一区| 久久久久久伊人网av| 日日啪夜夜爽| av片东京热男人的天堂| av线在线观看网站| 久久这里只有精品19| 黑人巨大精品欧美一区二区蜜桃 | 亚洲欧美一区二区三区国产| 久久人妻熟女aⅴ| 国产精品久久久久成人av| 久久久久久久久久人人人人人人| 大片免费播放器 马上看| 成人亚洲精品一区在线观看| 丝袜美足系列| 亚洲人成77777在线视频| 99久久综合免费| 美女脱内裤让男人舔精品视频| 午夜91福利影院| 夫妻性生交免费视频一级片| 91精品国产国语对白视频| 性高湖久久久久久久久免费观看| 久久国产亚洲av麻豆专区| 九草在线视频观看| 亚洲欧美日韩卡通动漫| 国产一区二区激情短视频 | 中文字幕人妻熟女乱码| 高清毛片免费看| 极品人妻少妇av视频| 内地一区二区视频在线| 日韩精品有码人妻一区| 亚洲精品国产av蜜桃| 免费观看a级毛片全部| 晚上一个人看的免费电影| 亚洲天堂av无毛| 天天躁夜夜躁狠狠躁躁| 日本wwww免费看| 男女午夜视频在线观看 | 国产爽快片一区二区三区| 精品一区二区三区四区五区乱码 | 女性被躁到高潮视频| 亚洲一码二码三码区别大吗| 亚洲精品乱码久久久久久按摩| 国产精品久久久久久精品电影小说| 大话2 男鬼变身卡| 丝袜脚勾引网站| 18在线观看网站| 亚洲av综合色区一区| 亚洲色图 男人天堂 中文字幕 | 欧美人与善性xxx| 国产免费视频播放在线视频| 免费观看av网站的网址| 在线免费观看不下载黄p国产| 国产 一区精品| 午夜精品国产一区二区电影| 国产成人91sexporn|