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

    利用MODIS多通道反照率產(chǎn)品估算OCO-2氧氣A吸收帶陸表反照率方法

    2023-05-17 05:03:30楊潔李四維王慶鑫
    遙感學(xué)報 2023年4期
    關(guān)鍵詞:反照率香農(nóng)方根

    楊潔, 李四維, 王慶鑫

    1. 武漢大學(xué) 遙感信息工程學(xué)院, 武漢 430079;2. 武漢大學(xué) 測繪遙感信息工程國家重點實驗室, 武漢 430079

    1 引 言

    衛(wèi)星接收到的輻亮度由大氣各組分與下墊面的反射共同組成,其中下墊面反射率作為大氣組分反演的重要干擾項,在痕量氣體、氣溶膠、云層等衛(wèi)星遙感中影響著定量反演的精度,通常需要提前確定,例如Veefkind等(2006)反演臭氧總量、Wang 等(2008)反演云量、云高的工作。因此,準(zhǔn)確估算下墊面反射率對大氣組分的反演至關(guān)重要。

    基于衛(wèi)星晴空觀測估算陸表反照率或反射率的方法已有較多研究。Herman 和Celarier (1997)首先基于TOMS(Total Ozone Mapping Spectrometer)數(shù)據(jù)構(gòu)建了340 nm 和380 nm 下墊面反射率數(shù)據(jù)庫,它能提供幾乎覆蓋全球范圍的1.25°×1.0°分辨率的數(shù)據(jù)。類似的,Koelemeijer 等(2003)、Kleipool 等(2008)、Tilstra 等(2017)基于GOME(Global Ozone Monitoring Experiment)、OMI (Ozone Monitoring Instrument)、GOME-2 和SCIAMACHY(Scanning Imaging Absorption Spectrometer for Atmospheric Chartography)的數(shù)據(jù)分別構(gòu)建了4 個分辨率大約為1.0°×1.0°的下墊面反射率數(shù)據(jù)庫。除數(shù)據(jù)源差異(不同的中心波長、空間分辨率,以及產(chǎn)品年份)外,它們都基于長時間衛(wèi)星觀測數(shù)據(jù),并假設(shè)陸表反射在年際間同比不變,采用最小值或眾數(shù)值等統(tǒng)計方法合成的朗伯等效反射率月平均產(chǎn)品??偟膩碚f,上述產(chǎn)品的時空分辨率較低,并且朗伯等效反射率在大氣遙感中可能導(dǎo)致不可忽視的反演偏差(Lorente等,2018)。

    云遙感與晴空下大氣痕量氣體或氣溶膠遙感的應(yīng)用場景不同,衛(wèi)星視野中的云層會阻礙對陸表的觀測,不能直接提供準(zhǔn)確的陸表反照率,但基于時序的估算方法能盡可能地填補數(shù)據(jù)空缺,例如基于MODIS(Moderate-resolution Imaging Spectroradiometer)多日晴空觀測構(gòu)建陸表二向反射率函數(shù)的MOD43 和MAIAC(Multi-Angle Implementation of Atmospheric Correction)算法(Lyapustin等,2018;Strahler等,2021)。

    除上述基于衛(wèi)星觀測構(gòu)建自身通道陸表反射數(shù)據(jù)庫的方法外,根據(jù)其他載荷或其他通道的陸表反射產(chǎn)品,也能轉(zhuǎn)換得到目標(biāo)載荷通道的陸表反射率或反照率數(shù)據(jù)。例如,在Kaufman 等(1997) 基 于MODIS 觀 測 和Ge 等(2019) 基 于AHI(Advanced Himawari Imager)觀測的氣溶膠反演工作中,都基于短波紅外波段的表觀反射率,通過線性變換得到了暗像元的藍光和紅光波段的地表反射率;在OMI 的NO2柱濃度反演和AATSR(Advanced Along-Track Scanning Radiometer)的氣溶膠反演算法中,都直接使用了MODIS 鄰近通道的陸表反射產(chǎn)品進行反演(Russell 等,2011;Shi等,2017);此外,在不同載荷的相似通道之間,也能在顧及儀器響應(yīng)的差異后,通過轉(zhuǎn)換得到目標(biāo)載荷的陸表反射率數(shù)據(jù)(Claverie 等,2018;Zhang等,2018;Wang等,2020)。

    本文的目標(biāo)是在考慮陸表各向異性反射的前提下,討論用于OCO-2(Orbiting Carbon Observatory 2)云屬性反演的氧氣A 吸收帶陸表反照率的估算方法。雖然梁順林等(2020)提出一種新的陸表反射率和反照率的估算方法,但目前仍缺少可用的OCO-2 氧氣A 吸收帶陸表反照率產(chǎn)品,一方面是因為OCO-2 衛(wèi)星的分辨率為1.29 km×2.25 km,遠高于上述基于GOME、GOME-2、SCIAMACHY等衛(wèi)星數(shù)據(jù)建立的陸表反射先驗數(shù)據(jù)庫(包含氧氣A 吸收帶)的空間分辨率,而使用欠采樣的陸表反照率會損害反演的質(zhì)量(Heckel等,2011);另一方面是因為OCO-2觀測的刈幅極窄(不足10 km)、重返周期很長(16日),無法采用基于時序的估算方法來填補有云天氣下的數(shù)據(jù)空缺。因而,當(dāng)下基于OCO-2 氧氣A 吸收帶的云反演或者氣溶膠反演的研究被局限于海洋區(qū)域(Zeng 等,2020;Yang等,2021),即使在陸地區(qū)域進行反演,其反演結(jié)果也有較大的偏差(Nanda 等,2020)。此外,Worden 等(2017)也指出當(dāng)前基于OCO-2 反演XCO2的精度也受到陸表反照率的干擾。

    獲取準(zhǔn)確的陸表反照率的難點源于地表的復(fù)雜多樣性。地表覆蓋類型在空間上的差異、在一個月內(nèi)的變化、陸表的二向反射特性,以及觀測通道的中心波長等都是確定陸表反照率所需要考慮的因素(Bacour等,2020)。MODIS的MCD43C3產(chǎn)品能提供與OCO-2 空間分辨率相近的、顧及了陸表反射各向異性的多通道陸表反照率,雖然缺少針對氧氣A 吸收帶的觀測通道,但MODIS 多通道反照率和OCO-2 氧氣A 吸收帶陸表反照率都是地表覆蓋類型在光學(xué)性質(zhì)上的表征,它們受相同的時空因素作用,也服從相似的二向反射規(guī)律。本文著重討論由MODIS 陸表反照率估計OCO-2 氧氣A 吸收帶陸表反照率的可行性,實現(xiàn)并測試轉(zhuǎn)換精度,分析轉(zhuǎn)換的誤差源。基于MODIS 陸表反射產(chǎn)品的轉(zhuǎn)換方法不僅能為OCO-2 的云遙感提供必需的陸表反照率,也可為基于OCO-2 觀測的氣溶膠或氣體反演提供必要的陸表反射率。

    2 數(shù) 據(jù)

    2.1 MCD43C3、MCD12C1產(chǎn)品

    搭載于Terra 和Aqua 兩顆衛(wèi)星上的MODIS 傳感器,裝備了0.4—14 μm 波長范圍內(nèi)的36 個觀測通道,刈幅寬達2330 km,每日的觀測幾乎覆蓋全球(Justice等,2002)。

    MCD43C3 產(chǎn)品由MODIS 觀測反演得到,記錄了地方時正午時刻多個通道的陸表反照率。反演算法針對可能存在的云層干擾,采用與有云日鄰近的無云日觀測來生產(chǎn)其陸表反照率數(shù)據(jù),以削弱云層等干擾(Wang 等,2018)。單個MCD43C3文件的數(shù)據(jù)覆蓋全球,空間分辨率為0.05°,時間分辨率為1 d。除MODIS觀測通道#1—7外,MCD43C3產(chǎn)品還提供了合成的寬波段黑空/白空反照率產(chǎn)品,包括可見光通道、近紅外通道和短波通道,各通道的波長范圍見表1。此外,產(chǎn)品中還有數(shù)據(jù)質(zhì)量標(biāo)識、正午時刻的太陽天頂角等輔助信息。

    表1 MCD43C3產(chǎn)品的光譜通道Table 1 Spectral bands of the MCD43C3 product

    MCD12C1 產(chǎn)品是基于MODIS 數(shù)據(jù)生產(chǎn)的地表覆蓋類型產(chǎn)品(Friedl等,2010),包含3種分類方案:國際地圈生物圈計劃(IGBP)全球植被分類方案(Loveland 和Belward,1997)、馬里蘭大學(xué)(UMD)植被分類方案(Hansen等,2000)、MODIS葉面積指數(shù)(LAI)方案(Lotsch等,2003)。IGBP采用的陸表覆蓋類型參見表2。單個MCD12C1 文件的數(shù)據(jù)覆蓋全球,空間分辨率為0.05°,時間分辨率為1 a。

    表2 國際地圈生物圈計劃(IGBP)全球植被分類方案Table 2 IGBP global vegetation classification scheme

    2.2 OCO2_L2_Lite_FP產(chǎn)品

    OCO-2 衛(wèi)星提供了氧氣A 吸收帶、二氧化碳強吸收帶和弱吸收帶的輻亮度觀測,其晴空觀測可用于反演二氧化碳產(chǎn)品,例如OCO2_L2_Lite_FP產(chǎn)品。這3 個通道的陸表反射特性(O’Dell 等,2018) 也與二氧化碳參數(shù)一同反演。然而,OCO2_L2_Lite_FP產(chǎn)品只提供了晴空或云量極少時的陸表反射數(shù)據(jù),因為只有篩選出的晴空觀測才能被用于反演二氧化碳參數(shù)。缺少有云時的陸表反照率數(shù)據(jù)是本文估計OCO-2 氧氣A 吸收帶陸表反照率的原因。

    OCO2_L2_Lite_FP 產(chǎn)品的反演算法采用形狀固定的反射率函數(shù)(Rahman 等,1993;Dubovik 等,2011),以及反射率函數(shù)的幅度來表征陸表的二向反射特性(O’Dell 等,2018)。除了反射率函數(shù)的幅度,OCO2_L2_Lite_FP產(chǎn)品還提供了觀測時刻的太陽天頂角等輔助信息。

    為分析MODIS 通道與OCO-2 氧氣A 吸收帶陸表反照率的變換規(guī)律,可基于時空最鄰近原則從MCD43C3 產(chǎn)品中挑選出與OCO2_L2_Lite_FP 產(chǎn)品時空匹配的MODIS 各通道的陸表反照率數(shù)據(jù)。需要注意的是,MCD43C3 產(chǎn)品提供的是正午時刻的陸表反照率,而OCO2_L2_Lite_FP 產(chǎn)品提供的是陸表反射率函數(shù),所以需要依據(jù)反射率函數(shù)計算正午時刻的陸表反照率,與MCD43C3 產(chǎn)品的黑空反照率相對應(yīng)。

    3 方 法

    本節(jié)首先論證由MODIS 向OCO-2 轉(zhuǎn)換陸表反照率的可行性,然后基于香農(nóng)熵的分析證明MODIS 的多通道反照率數(shù)據(jù)能提供地表覆蓋類型信息,且不影響上述轉(zhuǎn)換所能提供的最優(yōu)精度,最后介紹能實現(xiàn)上述精度的人工神經(jīng)網(wǎng)絡(luò)多通道模型。

    3.1 轉(zhuǎn)換的可行性

    為分析由MODIS 通道向OCO-2 氧氣A 吸收帶陸表反照率轉(zhuǎn)換的可行性,本文由2016 年內(nèi)時空匹配的OCO2_L2_Lite_FP 和MCD43C3 產(chǎn)品統(tǒng)計了兩者之間的相關(guān)系數(shù),如圖1所示。

    圖1 OCO-2衛(wèi)星氧氣A吸收帶的陸表反照率與MODIS多通道反照率間的相關(guān)系數(shù)Fig. 1 Correlation coefficients for the land surface albedos of the OCO-2 oxygen A-band and of the MODIS channels

    圖1 中氧氣A 吸收帶的陸表反照率與通道#2的黑空反照率的相關(guān)性最高,因為其中心波長(858 nm)與氧氣A 吸收帶的中心波長(770 nm)最為接近。兩者之間的高相關(guān)性(高達0.91)也說明可以由MODIS 通道#2 的反照率通過線性變換估計OCO-2 氧氣A 吸收帶的陸表反照率(基于樣本測試的均方根誤差為0.034)。此外,MODIS 的觀測通道(#1—#7)與合成的寬波段通道之間的相關(guān)系數(shù)也很高,說明可以只選擇觀測通道的反照率作為轉(zhuǎn)換模型的輸入。

    3.2 地表覆蓋類型的干擾

    將2016 年的陸表反照率數(shù)據(jù)按MCD12C1 產(chǎn)品提供的IGBP 地表覆蓋類型分組,并對分組后的數(shù)據(jù)集分別建立基于通道#2 的線性轉(zhuǎn)換模型,可得到16組不同的線性變換參數(shù),如圖2所示。圖2中的百分比代表各分組中數(shù)據(jù)的占比,其中分組7、8、9、10、12、16 的占比超過了5%,是主要的地表覆蓋類型。圖中橢圓代表各分組中反照率散點所聚集的區(qū)域,橢圓的長軸是散點的主方向,其長半軸的長度是散點在主方向上標(biāo)準(zhǔn)差的3倍,其短半軸的長度是散點在垂直于主方向上標(biāo)準(zhǔn)差的3倍。

    圖2 按地表覆蓋類型分組建立的線性轉(zhuǎn)換模型Fig. 2 The linear models based on land cover types

    橢圓長軸的方向代表對應(yīng)線性模型的斜率,右側(cè)的箱型圖說明了各分組間斜率的分散程度。除雪和冰類型(編號16)外,其他地表覆蓋類型的斜率在0.95—1.25,6 種主要地表覆蓋類型間的最大斜率差接近0.2。只有當(dāng)所有分組的線性變換參數(shù)(斜率、截距)盡可能地一致時,基于單通道的線性模型才能忽略地表覆蓋類型之間的差異。

    不一致的斜率說明按地表覆蓋類型分別建立的轉(zhuǎn)換模型比統(tǒng)一的轉(zhuǎn)換模型更有利于提升轉(zhuǎn)換的準(zhǔn)確度。然而,提供地表覆蓋類型的MCD12C1產(chǎn)品的時間間隔是1 a,且即使在年末也不能提供當(dāng)年的產(chǎn)品,所以不能將MCD12C1 產(chǎn)品中的地表覆蓋類型作為轉(zhuǎn)換模型的輸入。此外,較長的橢圓短半軸也說明基于單通道的簡單線性變換會造成較大的誤差,例如雪和冰類型。

    若在轉(zhuǎn)換模型中增加與通道#2 相關(guān)性較弱的通道#3 和通道#7,三通道的組合可能提供更多的信息。圖3展示了由三通道的黑空反照率繪制的散點圖,圖中橫軸是MODIS 通道#3 與通道#2 的比值,縱軸是通道#7與通道#2的比值。

    圖3 MODIS反照率比值的分布圖Fig. 3 Distribution of MODIS albedo ratios

    各地表覆蓋類型的散點,尤其是6種主要地表覆蓋類型的散點分布在圖3中不同的位置,且呈現(xiàn)聚集趨勢,說明MODIS 通道#2、通道#3 和通道#7的組合中蘊含著地表覆蓋類型的信息,那么建立基于MODIS 多通道反照率的轉(zhuǎn)換模型可能提供更加準(zhǔn)確的氧氣A 吸收帶陸表反照率數(shù)據(jù)。圖3 中6 種主要地表覆蓋類型的集中區(qū)域雖然相互重疊,但若引入更多通道的反照率數(shù)據(jù),則可能提升區(qū)分地表覆蓋類型的能力,也意味著更高的轉(zhuǎn)換精度。

    3.3 轉(zhuǎn)換的最優(yōu)精度

    本小節(jié)依據(jù)香農(nóng)熵量化不同輸入通道轉(zhuǎn)換模型的不確定度,并確定可實現(xiàn)的最優(yōu)轉(zhuǎn)換精度,以判定多通道反照率能否替代地表覆蓋類型,且提供相似的轉(zhuǎn)換精度。

    香農(nóng)熵是對信息所具有的不確定性的度量,其數(shù)值越大,說明這條信息的不確定度也越大,那么需要更多的信息才能進一步確定(Shannon,1997)。香農(nóng)熵(H)被定義為

    式中,函數(shù)P(xi)指事件xi發(fā)生的概率。例如,若已知所有N項事件的發(fā)生概率相同,那么由式(1)知,事件x1將發(fā)生這一表述的香農(nóng)熵為log2(N);同理,若已知其中事件x1的發(fā)生概率為100%,那么由式(1)知,事件x1將發(fā)生這一表述的香農(nóng)熵為-log2(1) = 0,即無需更多信息就能確認(rèn)。

    與之類似,香農(nóng)熵的概念也可用于量化氧氣A吸收帶陸表反照率估值的不確定性。圖4展示了以MODIS 通道#2 的黑空反照率間隔0.01 分組的氧氣A吸收帶陸表反照率的離散概率函數(shù)(太陽天頂角為30°)。相比于分組前,分組后各概率函數(shù)的跨度明顯減小,峰值位置也隨通道#2 的黑空反照率逐漸右移。圖5(a)中a.1 與之對應(yīng),展示了分組計算的香農(nóng)熵分布,集中在2.6—4.0,峰值位于3.0。除將太陽天頂角作為分組條件外,圖例BW#2指代將MODIS 通道#2 的黑/白空反照率也作為分組條件。

    圖4 OCO-2氧氣A吸收帶陸表反照率的離散概率函數(shù)Fig. 4 Discrete probability function of surface albedo on OCO-2 oxygen A-band

    為量化地表覆蓋類型信息對減小氧氣A 吸收帶陸表反照率估值的不確定度的貢獻,本節(jié)使用MCD12C1產(chǎn)品中的3種分類方案的結(jié)果(LC),以及年積日(Day)和經(jīng)緯度信息(Pos)來綜合表征地表覆蓋類型的光學(xué)特性。引入年積日是為了解釋地表覆蓋類型在一年內(nèi)隨季節(jié)的可能變化,例如落葉林的四季變化等。引入經(jīng)緯度是為了解釋某一類地表覆蓋類型在不同地理人文環(huán)境中的可能差異,例如不同氣候環(huán)境中種植的不同農(nóng)作物等。添加地表覆蓋類型后的香農(nóng)熵分布如圖5所示,集中在0.8—3.2,峰值位于2.4。減小的香農(nóng)熵說明引入地表覆蓋類型信息確實能減小氧氣A吸收帶陸表反照率估值的不確定度,這對在通道間準(zhǔn)確轉(zhuǎn)換陸表反照率是有利的。

    然而及時獲取地表覆蓋類型產(chǎn)品是有困難的,所以3.2小節(jié)提出使用MODIS的通道#2、通道#3和通道#7 的反照率數(shù)據(jù)替代地表覆蓋類型的信息,效果見圖5。使用三通道反照率的香農(nóng)熵分布與使用地表覆蓋類型的香農(nóng)熵分布(圖5)非常相似,也集中在0.8—3.2,峰值也位于2.4 附近。相似的香農(nóng)熵分布說明使用多通道反照率可以抑制多樣的地表覆蓋類型對估計氧氣A 吸收帶陸表反照率的干擾。

    本文也將MCD32C3 產(chǎn)品提供的7 個觀測通道的黑空反照率都作為分組條件也進行了測試,得到的香農(nóng)熵更小,主要分布在0.6—3.0,峰值位于1.4 附近,見圖5。左移的峰值說明除通道#2、通道#3 和通道#7 外,其他通道中仍蘊含有未被利用的信息。值為1.4 的香農(nóng)熵等價于跨度為0.026 的均勻分布概率函數(shù)的香農(nóng)熵,這可能是基于MODIS 多通道反照率數(shù)據(jù)轉(zhuǎn)換OCO-2 氧氣A 吸收帶陸表反照率所能達到的最高精度。

    圖5 也展示了同時使用MODIS 七通道黑空反照率和地表覆蓋類型信息進行氧氣A 吸收帶陸表反照率估計的香農(nóng)熵分布,其分布與不添加地表覆蓋類型信息的香農(nóng)熵分布相似,說明只使用MCD43C3 產(chǎn)品估計OCO-2 氧氣A 吸收帶陸表反照率是可行的,無需等待MCD12C1產(chǎn)品的發(fā)布。

    3.4 基于神經(jīng)網(wǎng)絡(luò)的多通道模型

    前文指出使用MODIS 七通道反照率來估計OCO-2 氧氣A 吸收帶陸表反照率能獲得最好的估計精度,并且由基于樣本統(tǒng)計的香農(nóng)熵分布可以判斷多數(shù)估值的絕對誤差應(yīng)在0.026以內(nèi)。

    人工神經(jīng)網(wǎng)絡(luò)非常適合將由大樣本統(tǒng)計的非線性規(guī)律轉(zhuǎn)換為計算模型(梁順林 等,2016)。本文選擇BP神經(jīng)網(wǎng)絡(luò)作為訓(xùn)練模型,以MCD43C3產(chǎn)品中的MODIS 七通道黑空/白空反照率以及正午時刻的太陽天頂角作為輸入,將由OCO2_L2_Lite_FP產(chǎn)品得到的正午時刻OCO-2 氧氣A 吸收帶陸表反照率作為輸出,學(xué)習(xí)它們之間的映射關(guān)系。圖6展示了測試結(jié)果,OCO-2 氧氣A 吸收帶陸表反照率估值的均方根誤差(RMSE)為0.026,相關(guān)系數(shù)(R)為0.94。該測試得到的均方根誤差(0.026)與3.3小節(jié)中得到的香農(nóng)熵(1.4)吻合,說明訓(xùn)練模型已充分利用輸入信息,達到了其所能提供的最佳估值精度。

    圖6 基于多通道模型的轉(zhuǎn)換結(jié)果與參考值的對比Fig. 6 Comparison of the reference albedo and the estimation from the multi-channel model

    陸表的性質(zhì)變化會影響多個波段上的輻射特性,因此多通道模型能更準(zhǔn)確地反映陸表的物化特性,并更準(zhǔn)確地估算陸表氧氣A 帶反照率以及其的變化規(guī)律。以反照率年變化較大的落葉闊葉林為例,圖7(a)展示了位于北溫帶(23°26′N—66°34′N)內(nèi)落葉闊葉林的MODIS通道#2黑空反照率(黑點),以及經(jīng)多項式擬合得到的年變化曲線(藍線)。此外,還展示了經(jīng)多項式擬合的氧氣A吸收帶估算值的年變化曲線(紅線)。通道#2 與氧氣A 吸收帶的反照率在12、1、2 月因葉片完全脫落而位于低值,在3—5 月間因生長出新葉而迅速上升,于6—8 月達到頂峰,在9—11 月因落葉而下降,周而復(fù)始。4—9月間,擬合的通道#2與氧氣A 吸收帶的陸表反照率的差異較小,這與圖7(b)中4—9 月紅邊斜率較大是一致的,因為期間闊葉林的葉面積大,對近紅外輻射的反射強。

    圖7 2017年落葉闊葉林譜反照率的變化Fig. 7 The spectral albedo of deciduous broad-leaved forest in 2017

    4 結(jié)果與分析

    本文提出的多通道模型本質(zhì)上是一個統(tǒng)計模型,統(tǒng)計模型的訓(xùn)練集和測試集應(yīng)避免在時間和空間上重疊,否則不利于評價統(tǒng)計模型的可用性。

    4.1 基于時空的驗證

    為說明3.4 小節(jié)中提出的多通道模型并不依賴地理因素,即多通道模型不會因為測試數(shù)據(jù)與訓(xùn)練數(shù)據(jù)在空間上不重疊而失效,設(shè)計如下實驗:采用2016 年經(jīng)度范圍180°—90°W、0°—90°E 的數(shù)據(jù)作為訓(xùn)練集以訓(xùn)練多通道模型,采用2016 年經(jīng)度范圍90°W—0°、90°E—180°的數(shù)據(jù)作為測試集以評估其準(zhǔn)確性,測試結(jié)果如圖8所示。訓(xùn)練集的均方根誤差為0.024,相關(guān)系數(shù)為0.96;測試集的結(jié)果與之非常接近,均方根誤差為0.026,相關(guān)系數(shù)為0.93。

    圖8 對多通道模型基于空間的驗證Fig. 8 Validation of the multi-channel model on new space

    為說明多通道模型并不依賴時間因素,即多通道模型不會因為測試數(shù)據(jù)與訓(xùn)練數(shù)據(jù)在時間上不重疊而失效,設(shè)計如下實驗:采用2016 年的數(shù)據(jù)作為訓(xùn)練集訓(xùn)練多通道模型,采用2017 年的數(shù)據(jù)作為測試集評估其準(zhǔn)確性,測試如圖9所示。訓(xùn)練集的均方根誤差為0.025,相關(guān)系數(shù)為0.95;測試集的結(jié)果與之非常接近,均方根誤差為0.026,相關(guān)系數(shù)為0.94。若采用2018 年或2019 年的數(shù)據(jù)測試,氧氣A 吸收帶陸表反照率估值的均方根誤差均為0.025,相關(guān)系數(shù)分別為0.94和0.95。

    圖9 對多通道模型基于時間的驗證Fig. 9 Validation of the multi-channel model on new time

    上述基于時空的驗證充分說明多通道模型在實踐中的可用性,即使訓(xùn)練數(shù)據(jù)與測試數(shù)據(jù)在空間和時間上都不重疊,多通道模型也表現(xiàn)較好。例如,使用2016 年180°—90°W,0°—90°E 范圍內(nèi)的數(shù)據(jù)訓(xùn)練多通道模型,在2019 年90°W—0°,90°E—180°范圍內(nèi)的數(shù)據(jù)上進行測試,測試結(jié)果的均方根誤差為0.025,相關(guān)系數(shù)為0.94。

    在地表覆蓋類型有變化的區(qū)域,時空無關(guān)的模型能基于最新的MODIS 觀測資料,依靠其掌握的各地表覆蓋類型在多通道反照率上的轉(zhuǎn)換規(guī)律,而不是反照率的時空變化規(guī)律,對OCO-2 氧氣A吸收帶陸表反照率進行估計。時空無關(guān)的模型能發(fā)現(xiàn)地物在時空上的變化,而基于歷史、基于地理位置的模型并不適合捕捉這樣的變化。

    4.2 基于地表類型的驗證

    圖10 展示了多通道模型在各地表覆蓋類型上的表現(xiàn)。多數(shù)地表覆蓋類型的均方根誤差小于0.03,均方根誤差最小的地表覆蓋類型是密灌叢,值為0.015,只有雪和冰類型的均方根誤差遠大于其他類型,值為0.070。多數(shù)類型的相關(guān)系數(shù)大于0.85,相關(guān)系數(shù)最大的地表覆蓋類型是裸地或低植被覆蓋地,值為0.96,雪和冰類型的相關(guān)系數(shù)最小,值為0.84。雪和冰類型可能與其他地表覆蓋類型的二向反射特性有較大的差異,但在OCO2_L2_Lite_FP的反演算法中并未進行區(qū)分,使用的是同一個陸表反射率形狀函數(shù),這可能是導(dǎo)致雪和冰類型所對應(yīng)的氧氣A 吸收帶陸表反照率估值不準(zhǔn)確的原因之一。

    圖10 按IGBP地表覆蓋類型分組展示的多通道模型的測試結(jié)果Fig. 10 Test results of multi-channel models grouped by IGBP land cover type

    改用多通道模型后,OCO-2 氧氣A 吸收帶陸表反照率估值的準(zhǔn)確度在各分組的測試中都有提升(改進前后的對比見表3)。相比于單通道線性模型,多通道模型在所有分組上的均方根誤差都有所降低,相關(guān)系數(shù)也都有所升高。其中,裸地或低植被覆蓋地分組的相關(guān)系數(shù)在兩個模型中都是最高的,但多通道模型的相關(guān)系數(shù)和均方根誤差仍有改善。即使在表現(xiàn)不佳的雪和冰分組中,采用多通道模型也能改善估值的相關(guān)系數(shù)和均方根誤差。

    表3 單通道線性模型與多通道模型的性能在各地表覆蓋類型上的對比Table 3 Comparison between single-channel linear model and multi-channel model in each land cover type

    4.3 輸入數(shù)據(jù)質(zhì)量與轉(zhuǎn)換誤差

    轉(zhuǎn)換結(jié)果的誤差來源包括多通道模型的誤差和輸入數(shù)據(jù)的誤差。作為參考值的氧氣A 吸收帶陸表反照率來自O(shè)CO-2 衛(wèi)星產(chǎn)品,也可能是不夠準(zhǔn)確的。

    評估轉(zhuǎn)換精度的參考值是由OCO2_L2_Lite_FP產(chǎn)品中假設(shè)的反射率函數(shù)(ρ)的固定形狀和反演的幅度參數(shù)(s)決定的。OCO2_L2_Lite_FP 反演不區(qū)分陸地上的地表覆蓋類型,使用同一個形狀函數(shù),真實陸表上形狀不同的反射率函數(shù)可能是氧氣A 吸收帶陸表反照率參考值不準(zhǔn)確的原因之一。此外,反射率函數(shù)的幅度參數(shù)是由晴空條件下單一角度的觀測反演得到,而反照率是對所有反射方向的積分,缺少多角度觀測的約束也是導(dǎo)致氧氣A 吸收帶陸表反照率參考值不準(zhǔn)確的原因之一。

    轉(zhuǎn)換結(jié)果的誤差也可能來自MCD43C3 產(chǎn)品的七通道反照率,它們是多通道模型的輸入。在生成MCD43C3 的反照率格網(wǎng)化產(chǎn)品時,按網(wǎng)格內(nèi)反照率的反演值和空缺值的比例定義了其產(chǎn)品質(zhì)量。標(biāo)識0代表最佳的質(zhì)量,其判定條件是網(wǎng)格內(nèi)的反照率數(shù)據(jù)都是由觀測反演得到。標(biāo)識1—4 對應(yīng)的產(chǎn)品質(zhì)量逐步下降,其中反演數(shù)據(jù)的比例逐漸降低。按MCD43C3 產(chǎn)品的質(zhì)量標(biāo)識分組統(tǒng)計均方根誤差,結(jié)果如圖11 所示。在輸入數(shù)據(jù)質(zhì)量最好的分組中,OCO-2 氧氣A 吸收帶估值的均方根誤差大約是0.02,質(zhì)量標(biāo)識為3 和4 的輸入數(shù)據(jù)的均方根誤差都超過0.05。

    圖11 基于不同質(zhì)量的MCD43C3產(chǎn)品估算氧氣A帶反照率的均方根誤差Fig. 11 The RMSE of oxygen A band albedo estimated based on MCD43C3 products of each quality

    將上述估值應(yīng)用于需要最優(yōu)估計的遙感方法中時,參考圖11可由MCD43C3的產(chǎn)品質(zhì)量標(biāo)識確定OCO-2 氧氣A 吸收帶陸表反照率的先驗誤差。此外,當(dāng)前也缺少對OCO2_L2_Lite_FP 產(chǎn)品中陸表反照率數(shù)據(jù)質(zhì)量的驗證研究,而圖11 至少可以說明產(chǎn)品中氧氣A 吸收帶的陸表反照率數(shù)據(jù)在可能受到云干擾的前提下,均方根誤差小于0.06;在晴空天氣下,數(shù)據(jù)的均方根誤差小于0.02。

    圖10的各子圖中都存在遠離聚集區(qū)域的散點,它們對應(yīng)于MCD43C3 產(chǎn)品中質(zhì)量較差的數(shù)據(jù)。按MCD43C3 產(chǎn)品質(zhì)量標(biāo)識與地表覆蓋類型分別統(tǒng)計各分組中的均方根誤差與相關(guān)系數(shù),結(jié)果如表4所示。當(dāng)輸入最佳質(zhì)量的數(shù)據(jù)時,超過半數(shù)的地表覆蓋類型所對應(yīng)的均方根誤差小于0.02,雪和冰類型的均方根誤差也大幅下降致0.029,與其他類型的均方根誤差也相差不大。對于產(chǎn)品質(zhì)量為1—4 的數(shù)據(jù),半數(shù)地表覆蓋類型的均方根誤差大于0.035,雪和冰類型的均方根誤差為0.076,遠大于其他類型的均方根誤差。圖11 和表4 說明輸入的MODIS 七通道反照率數(shù)據(jù)的質(zhì)量是決定OCO-2 氧氣A吸收帶陸表反照率估值精度最重要的因素。

    表4 在各地表覆蓋類型上不同質(zhì)量的輸入導(dǎo)致的誤差Table 4 Error due to different quality of input data

    5 結(jié) 論

    基于相關(guān)性的分析可知,根據(jù)MODIS 通道#2的反照率數(shù)據(jù)可以估計OCO-2 氧氣A 吸收帶的陸表反照率,但地表覆蓋類型是干擾通道間陸表反照率轉(zhuǎn)換的重要因素。在不同輸入的轉(zhuǎn)換模型間進一步比較香農(nóng)熵后可知:(1)基于MOIDS通道#2 反照率和地表覆蓋類型的OCO-2 反照率估值的不確定度與基于MODIS 三通道(#2,#3,#7)的不確定度相似;(2)基于MODIS 七通道(#1—#7)的不確定度更小,并且再加入地表覆蓋類型也不能減小估值的不確定度;(3)基于七通道統(tǒng)計的香農(nóng)熵峰值為1.4,所以多數(shù)OCO-2 氧氣A 吸收帶反照率估值的絕對誤差應(yīng)在0.026 以內(nèi)?;谏鲜龇治?,本文提出了基于BP 神經(jīng)網(wǎng)絡(luò)的多通道模型,該模型以MCD43C3 產(chǎn)品中的MODIS 七通道黑空/白空反照率以及正午時刻的太陽天頂角作為輸入,使用由OCO2_L2_Lite_FP 產(chǎn)品得到的正午時刻陸表反照率作為輸出進行訓(xùn)練。

    為驗證基于人工神經(jīng)網(wǎng)絡(luò)的多通道模型在不同時間和不同空間上的可用性。本文在覆蓋半個地球的數(shù)據(jù)集訓(xùn)練模型,在另半個地球上的數(shù)據(jù)集上進行測試;另外在2016年的全球數(shù)據(jù)集訓(xùn)練,在2017 年全球數(shù)據(jù)集上測試;兩次測試的相關(guān)系數(shù)都超過0.93,均方根誤差都是0.026。此外,本文也驗證了多通道模型在不同地表覆蓋類型上的有效性。在分組統(tǒng)計的相關(guān)系數(shù)和均方根誤差中,不論是表現(xiàn)最好的裸地或低植被覆蓋地類型,還是較差的雪和冰類型,多通道模型都優(yōu)于單通道線性模型。

    多通道模型的轉(zhuǎn)換誤差與輸入數(shù)據(jù)的產(chǎn)品質(zhì)量密切相關(guān)?;谧罴奄|(zhì)量的MCD43C3 產(chǎn)品的轉(zhuǎn)換結(jié)果的均方根誤差優(yōu)于0.02,均方根誤差隨輸入數(shù)據(jù)產(chǎn)品質(zhì)量的下降而增大,直至超過0.05。此外,上述統(tǒng)計得到的均方根誤差可能與氧氣A吸收帶陸表反照率參考值的誤差有關(guān)。

    本文提供了一種估算OCO-2 氧氣A 吸收帶陸表反照率的新思路,該思路不僅可應(yīng)用于OCO-2對云的遙感中地表反照率影響的去除,還可有助于OCO-2 對氣體、氣溶膠的反演中地表反照率的確定。

    猜你喜歡
    反照率香農(nóng)方根
    基于藍天空反照率的氣溶膠輻射強迫模擬
    方根拓展探究
    薩吾爾山木斯島冰川反照率時空變化特征研究
    冰川凍土(2022年6期)2022-02-12 08:31:06
    大衛(wèi),不可以
    長江三角洲地區(qū)大氣氣溶膠柱單次散射反照率特性研究
    校園恩仇錄:小混混和易拉罐女王的故事
    均方根嵌入式容積粒子PHD 多目標(biāo)跟蹤方法
    艾米麗的呼嚕
    揭開心算方根之謎
    基于香農(nóng)熵的超細(xì)粉體填料混合均勻度的評價研究
    中國塑料(2015年9期)2015-10-14 01:12:18
    a级毛片黄视频| 亚洲精品国产区一区二| 久久久精品94久久精品| 视频在线观看一区二区三区| h视频一区二区三区| 亚洲一区中文字幕在线| 国产精品 欧美亚洲| 黄片大片在线免费观看| 亚洲精品美女久久久久99蜜臀| 久久久久久久大尺度免费视频| 亚洲国产精品一区三区| a在线观看视频网站| 国产91精品成人一区二区三区 | 高清欧美精品videossex| 女人精品久久久久毛片| 高清黄色对白视频在线免费看| 日韩欧美一区视频在线观看| 国产亚洲av高清不卡| 亚洲精品久久午夜乱码| 最黄视频免费看| 国产精品自产拍在线观看55亚洲 | 久久精品亚洲av国产电影网| 高清视频免费观看一区二区| 纯流量卡能插随身wifi吗| 人妻人人澡人人爽人人| 亚洲国产中文字幕在线视频| 十八禁人妻一区二区| 久久久水蜜桃国产精品网| 韩国高清视频一区二区三区| 久久 成人 亚洲| 国产一区二区在线观看av| 国产不卡av网站在线观看| 岛国在线观看网站| 欧美日韩成人在线一区二区| 亚洲中文av在线| 美女午夜性视频免费| 国产精品99久久99久久久不卡| 黑人欧美特级aaaaaa片| 精品一区二区三区四区五区乱码| 久久99一区二区三区| 日本精品一区二区三区蜜桃| 嫩草影视91久久| 亚洲一区中文字幕在线| 我的亚洲天堂| 免费不卡黄色视频| 母亲3免费完整高清在线观看| a级片在线免费高清观看视频| 男人舔女人的私密视频| 两性午夜刺激爽爽歪歪视频在线观看 | 可以免费在线观看a视频的电影网站| 亚洲国产精品一区二区三区在线| 亚洲精品久久成人aⅴ小说| 亚洲精品中文字幕在线视频| 精品国产乱码久久久久久小说| 免费在线观看完整版高清| 日韩欧美国产一区二区入口| 免费日韩欧美在线观看| 青青草视频在线视频观看| 女人精品久久久久毛片| 超色免费av| 国产成人精品久久二区二区91| 亚洲,欧美精品.| 97在线人人人人妻| 99热全是精品| 午夜91福利影院| 中文字幕最新亚洲高清| 亚洲国产av影院在线观看| a级毛片在线看网站| 精品少妇黑人巨大在线播放| 搡老岳熟女国产| 亚洲国产欧美日韩在线播放| 我要看黄色一级片免费的| 搡老乐熟女国产| 国产一区二区在线观看av| 亚洲精华国产精华精| 女人精品久久久久毛片| 亚洲av电影在线观看一区二区三区| 操美女的视频在线观看| 久久久久视频综合| 中国国产av一级| 亚洲精品国产一区二区精华液| 法律面前人人平等表现在哪些方面 | 丝袜美腿诱惑在线| 欧美老熟妇乱子伦牲交| 午夜福利,免费看| 极品少妇高潮喷水抽搐| 日韩一区二区三区影片| 亚洲自偷自拍图片 自拍| 久久综合国产亚洲精品| 国产精品亚洲av一区麻豆| 国产一区二区三区av在线| 日本a在线网址| 欧美激情极品国产一区二区三区| 久久精品aⅴ一区二区三区四区| 成年av动漫网址| 国产日韩欧美视频二区| 老司机深夜福利视频在线观看 | 他把我摸到了高潮在线观看 | 人妻 亚洲 视频| 午夜精品国产一区二区电影| 午夜91福利影院| 一本大道久久a久久精品| 91精品伊人久久大香线蕉| 国精品久久久久久国模美| 久久午夜综合久久蜜桃| 久久久精品94久久精品| www.精华液| 熟女少妇亚洲综合色aaa.| 少妇的丰满在线观看| 18禁裸乳无遮挡动漫免费视频| 日韩制服骚丝袜av| 日本一区二区免费在线视频| 国产精品偷伦视频观看了| 美女大奶头黄色视频| 纵有疾风起免费观看全集完整版| 少妇裸体淫交视频免费看高清 | 欧美大码av| 欧美大码av| 国产高清videossex| 免费久久久久久久精品成人欧美视频| 婷婷丁香在线五月| 老司机午夜福利在线观看视频 | 欧美日韩亚洲综合一区二区三区_| 亚洲av日韩精品久久久久久密| www.av在线官网国产| 一二三四社区在线视频社区8| 亚洲精品美女久久久久99蜜臀| 一个人免费在线观看的高清视频 | 老司机午夜十八禁免费视频| 91精品国产国语对白视频| 中文字幕制服av| 久久久精品免费免费高清| 国产欧美亚洲国产| 操出白浆在线播放| 精品第一国产精品| 9热在线视频观看99| 亚洲伊人色综图| 黄色视频不卡| 成年人黄色毛片网站| 亚洲欧洲精品一区二区精品久久久| av电影中文网址| 老熟女久久久| 一个人免费在线观看的高清视频 | 中文字幕色久视频| 9色porny在线观看| 9热在线视频观看99| 欧美日韩亚洲高清精品| 99re6热这里在线精品视频| 国产欧美亚洲国产| 欧美激情高清一区二区三区| 国产欧美亚洲国产| 男女之事视频高清在线观看| 久久久久久亚洲精品国产蜜桃av| av在线播放精品| 久久精品国产综合久久久| 久热爱精品视频在线9| 精品亚洲成a人片在线观看| 久久精品人人爽人人爽视色| 亚洲人成电影免费在线| 久久 成人 亚洲| 亚洲精品在线美女| 亚洲国产精品成人久久小说| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美日本中文国产一区发布| 欧美在线一区亚洲| 亚洲国产av新网站| 婷婷丁香在线五月| 在线 av 中文字幕| 亚洲精品中文字幕在线视频| 免费黄频网站在线观看国产| 亚洲专区中文字幕在线| 一边摸一边做爽爽视频免费| 久久久久久人人人人人| 不卡一级毛片| 18禁国产床啪视频网站| 欧美在线黄色| 成人亚洲精品一区在线观看| 午夜影院在线不卡| 婷婷成人精品国产| 亚洲三区欧美一区| www.999成人在线观看| 国产欧美日韩一区二区精品| 一本色道久久久久久精品综合| 大片电影免费在线观看免费| 亚洲中文日韩欧美视频| 亚洲成人免费电影在线观看| 男女午夜视频在线观看| 老熟妇乱子伦视频在线观看 | 欧美精品亚洲一区二区| 欧美日韩亚洲高清精品| 啦啦啦免费观看视频1| avwww免费| 丝袜美腿诱惑在线| 亚洲欧美一区二区三区黑人| 精品国产一区二区久久| 精品一区二区三区av网在线观看 | 亚洲天堂av无毛| 美女扒开内裤让男人捅视频| 亚洲欧美日韩高清在线视频 | 动漫黄色视频在线观看| 好男人电影高清在线观看| 老汉色av国产亚洲站长工具| 99精品欧美一区二区三区四区| 精品人妻熟女毛片av久久网站| 成人亚洲精品一区在线观看| 亚洲国产av新网站| 熟女少妇亚洲综合色aaa.| 91老司机精品| 精品国产国语对白av| 在线十欧美十亚洲十日本专区| 免费黄频网站在线观看国产| 最近最新中文字幕大全免费视频| 最新的欧美精品一区二区| 久9热在线精品视频| 夜夜夜夜夜久久久久| 久久久精品94久久精品| 日日爽夜夜爽网站| 女性被躁到高潮视频| 亚洲精品国产区一区二| 一个人免费看片子| 高清视频免费观看一区二区| 国产男女超爽视频在线观看| 久久免费观看电影| 亚洲欧美成人综合另类久久久| 美女扒开内裤让男人捅视频| 看免费av毛片| 亚洲va日本ⅴa欧美va伊人久久 | 一区二区日韩欧美中文字幕| 国产亚洲欧美在线一区二区| 精品福利永久在线观看| 欧美日韩一级在线毛片| 中文精品一卡2卡3卡4更新| 高潮久久久久久久久久久不卡| 我要看黄色一级片免费的| 午夜两性在线视频| av线在线观看网站| 80岁老熟妇乱子伦牲交| 人人妻人人澡人人爽人人夜夜| 欧美精品人与动牲交sv欧美| 亚洲人成电影免费在线| 国产成人啪精品午夜网站| 国产成人免费无遮挡视频| 曰老女人黄片| 欧美+亚洲+日韩+国产| 久久精品国产亚洲av香蕉五月 | 91精品伊人久久大香线蕉| 国产精品久久久久久精品电影小说| 悠悠久久av| 久久亚洲精品不卡| 色94色欧美一区二区| 免费高清在线观看视频在线观看| 免费不卡黄色视频| 国产成人av激情在线播放| 五月天丁香电影| 少妇粗大呻吟视频| 久久精品aⅴ一区二区三区四区| 91成年电影在线观看| 1024视频免费在线观看| 亚洲伊人久久精品综合| 成人18禁高潮啪啪吃奶动态图| 欧美人与性动交α欧美软件| 亚洲国产成人一精品久久久| 久久精品aⅴ一区二区三区四区| 精品久久久久久久毛片微露脸 | 国产精品欧美亚洲77777| 亚洲av成人一区二区三| 日韩一区二区三区影片| 波多野结衣av一区二区av| 国产欧美日韩精品亚洲av| 伊人亚洲综合成人网| 两人在一起打扑克的视频| 黄色a级毛片大全视频| 国产精品久久久久久精品电影小说| 丝瓜视频免费看黄片| 久久中文看片网| 精品免费久久久久久久清纯 | 99久久人妻综合| 黄片小视频在线播放| 香蕉丝袜av| 久9热在线精品视频| svipshipincom国产片| 丰满少妇做爰视频| 飞空精品影院首页| 国产av一区二区精品久久| 免费女性裸体啪啪无遮挡网站| 日韩有码中文字幕| 搡老乐熟女国产| 国产老妇伦熟女老妇高清| 丝瓜视频免费看黄片| 十八禁网站网址无遮挡| 免费看十八禁软件| 国产男人的电影天堂91| 色播在线永久视频| 热re99久久精品国产66热6| 国产主播在线观看一区二区| 19禁男女啪啪无遮挡网站| 美女高潮到喷水免费观看| 国产91精品成人一区二区三区 | 狂野欧美激情性xxxx| 亚洲av成人一区二区三| 久久天堂一区二区三区四区| 麻豆国产av国片精品| 搡老乐熟女国产| 国产高清视频在线播放一区 | 丰满少妇做爰视频| www.999成人在线观看| 夫妻午夜视频| 日韩视频在线欧美| 老司机午夜十八禁免费视频| 国产一级毛片在线| 永久免费av网站大全| 免费av中文字幕在线| 欧美一级毛片孕妇| 精品亚洲成国产av| 三级毛片av免费| 丰满饥渴人妻一区二区三| 欧美在线一区亚洲| 午夜福利影视在线免费观看| 岛国毛片在线播放| 精品国内亚洲2022精品成人 | 女性被躁到高潮视频| 亚洲国产日韩一区二区| 亚洲av电影在线观看一区二区三区| 国产精品久久久av美女十八| 欧美黄色淫秽网站| 新久久久久国产一级毛片| 别揉我奶头~嗯~啊~动态视频 | 成年人黄色毛片网站| a级毛片黄视频| 国产一区二区三区在线臀色熟女 | 男女之事视频高清在线观看| 国产在视频线精品| 在线观看免费视频网站a站| 亚洲国产欧美一区二区综合| 国产精品久久久久久精品电影小说| 国产伦人伦偷精品视频| 高清欧美精品videossex| 五月开心婷婷网| 9色porny在线观看| 人妻一区二区av| 欧美少妇被猛烈插入视频| 亚洲国产毛片av蜜桃av| 爱豆传媒免费全集在线观看| 精品人妻1区二区| av福利片在线| 三上悠亚av全集在线观看| 中文欧美无线码| 亚洲欧美日韩另类电影网站| 久久精品成人免费网站| 黑人巨大精品欧美一区二区mp4| 俄罗斯特黄特色一大片| 搡老熟女国产l中国老女人| 亚洲中文av在线| 淫妇啪啪啪对白视频 | 99久久精品国产亚洲精品| 午夜精品国产一区二区电影| 99热国产这里只有精品6| 亚洲精品中文字幕在线视频| 欧美日韩一级在线毛片| 国产欧美日韩一区二区三区在线| 欧美日韩一级在线毛片| www.自偷自拍.com| 国产一区二区在线观看av| 狠狠精品人妻久久久久久综合| 亚洲美女黄色视频免费看| 国产精品国产三级国产专区5o| 一级,二级,三级黄色视频| 9热在线视频观看99| 国产免费福利视频在线观看| 亚洲一区中文字幕在线| 国产精品久久久久久人妻精品电影 | 精品一区二区三区四区五区乱码| 女人高潮潮喷娇喘18禁视频| 久久久久久亚洲精品国产蜜桃av| 免费黄频网站在线观看国产| av网站在线播放免费| 国产日韩欧美在线精品| 久久久国产一区二区| 巨乳人妻的诱惑在线观看| 狂野欧美激情性xxxx| 午夜影院在线不卡| 2018国产大陆天天弄谢| 成年动漫av网址| 午夜激情久久久久久久| 一区福利在线观看| 国产成人啪精品午夜网站| 成年动漫av网址| svipshipincom国产片| 国产主播在线观看一区二区| av一本久久久久| 80岁老熟妇乱子伦牲交| av天堂久久9| 国产不卡av网站在线观看| 国产高清videossex| 他把我摸到了高潮在线观看 | 一边摸一边做爽爽视频免费| 久久精品亚洲av国产电影网| 国产亚洲一区二区精品| 男人操女人黄网站| 亚洲专区中文字幕在线| 日本vs欧美在线观看视频| 久久人人97超碰香蕉20202| 99久久国产精品久久久| 久久人人97超碰香蕉20202| 菩萨蛮人人尽说江南好唐韦庄| 99久久人妻综合| 极品人妻少妇av视频| 国产野战对白在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 日韩视频在线欧美| 成人影院久久| 中文字幕人妻丝袜一区二区| 日日摸夜夜添夜夜添小说| 岛国毛片在线播放| av电影中文网址| 在线看a的网站| 一边摸一边抽搐一进一出视频| 人成视频在线观看免费观看| 秋霞在线观看毛片| 他把我摸到了高潮在线观看 | 亚洲成人免费电影在线观看| 亚洲天堂av无毛| 亚洲国产成人一精品久久久| 日本wwww免费看| 丝袜喷水一区| 日本猛色少妇xxxxx猛交久久| 人妻 亚洲 视频| 国产人伦9x9x在线观看| 亚洲久久久国产精品| 免费一级毛片在线播放高清视频 | 精品久久蜜臀av无| 午夜免费观看性视频| av一本久久久久| 高清视频免费观看一区二区| 久久国产亚洲av麻豆专区| 久久中文字幕一级| 国产精品一区二区精品视频观看| 大香蕉久久成人网| 999久久久国产精品视频| 夜夜夜夜夜久久久久| 国产日韩欧美在线精品| 欧美大码av| 精品少妇黑人巨大在线播放| 亚洲欧美一区二区三区久久| av网站免费在线观看视频| 亚洲精品成人av观看孕妇| 久久精品亚洲av国产电影网| 嫩草影视91久久| 两人在一起打扑克的视频| 在线观看一区二区三区激情| 黑人猛操日本美女一级片| 日本精品一区二区三区蜜桃| 中文字幕高清在线视频| 99国产综合亚洲精品| 老司机亚洲免费影院| 最近最新免费中文字幕在线| 国产精品久久久久久精品古装| 美女主播在线视频| 国产亚洲精品第一综合不卡| 黑人猛操日本美女一级片| 在线精品无人区一区二区三| 女人爽到高潮嗷嗷叫在线视频| 日本a在线网址| 九色亚洲精品在线播放| 制服诱惑二区| 国产深夜福利视频在线观看| 久久久久国产一级毛片高清牌| 久热这里只有精品99| 国产一区二区 视频在线| 亚洲精品粉嫩美女一区| 中文字幕精品免费在线观看视频| 国产成人影院久久av| 悠悠久久av| 日本撒尿小便嘘嘘汇集6| avwww免费| 亚洲国产精品一区三区| 欧美日韩中文字幕国产精品一区二区三区 | 久久综合国产亚洲精品| 午夜福利免费观看在线| 女性被躁到高潮视频| 国产精品成人在线| 久久精品国产a三级三级三级| 久久久久国产一级毛片高清牌| 国产91精品成人一区二区三区 | 黑人巨大精品欧美一区二区蜜桃| 婷婷丁香在线五月| 搡老岳熟女国产| 国产人伦9x9x在线观看| 91成年电影在线观看| 欧美性长视频在线观看| 欧美亚洲 丝袜 人妻 在线| 好男人电影高清在线观看| 狂野欧美激情性xxxx| 国产在线观看jvid| 中文字幕最新亚洲高清| 国产精品.久久久| 欧美人与性动交α欧美软件| 最新的欧美精品一区二区| 老熟妇仑乱视频hdxx| 首页视频小说图片口味搜索| 久久天堂一区二区三区四区| 亚洲av成人一区二区三| 男女边摸边吃奶| 女人被躁到高潮嗷嗷叫费观| 亚洲综合色网址| 久久av网站| 大片免费播放器 马上看| 国产成人精品久久二区二区91| 国产精品一二三区在线看| 日本撒尿小便嘘嘘汇集6| 午夜福利一区二区在线看| 一级毛片女人18水好多| 91av网站免费观看| a 毛片基地| tocl精华| tube8黄色片| 国产精品一区二区在线不卡| 99精国产麻豆久久婷婷| 国产精品免费视频内射| 国产欧美日韩综合在线一区二区| 国产日韩欧美在线精品| 日本av手机在线免费观看| 免费观看av网站的网址| 国产精品一二三区在线看| 夜夜夜夜夜久久久久| a级片在线免费高清观看视频| 热re99久久国产66热| 欧美在线一区亚洲| 亚洲专区中文字幕在线| svipshipincom国产片| 女性生殖器流出的白浆| 国产精品1区2区在线观看. | 免费高清在线观看日韩| 国产视频一区二区在线看| 久久久精品国产亚洲av高清涩受| 午夜免费观看性视频| 91精品伊人久久大香线蕉| 欧美在线一区亚洲| 91字幕亚洲| 国产激情久久老熟女| 国产亚洲欧美在线一区二区| 国产精品久久久久久精品电影小说| av福利片在线| 80岁老熟妇乱子伦牲交| 精品国产国语对白av| 麻豆av在线久日| 嫩草影视91久久| 亚洲精品国产色婷婷电影| 国产免费现黄频在线看| 国产亚洲av高清不卡| 亚洲免费av在线视频| 亚洲成人国产一区在线观看| 夜夜骑夜夜射夜夜干| 桃红色精品国产亚洲av| 国产亚洲精品第一综合不卡| 青青草视频在线视频观看| 亚洲欧美一区二区三区黑人| 婷婷色av中文字幕| av又黄又爽大尺度在线免费看| 亚洲精品美女久久av网站| 色婷婷久久久亚洲欧美| 久久久国产成人免费| 啦啦啦免费观看视频1| 久久99热这里只频精品6学生| 免费女性裸体啪啪无遮挡网站| 啦啦啦啦在线视频资源| 日韩中文字幕视频在线看片| 一二三四社区在线视频社区8| 五月天丁香电影| 亚洲一卡2卡3卡4卡5卡精品中文| 男人添女人高潮全过程视频| 中文字幕av电影在线播放| 欧美精品一区二区大全| 好男人电影高清在线观看| 亚洲欧洲精品一区二区精品久久久| 俄罗斯特黄特色一大片| 999精品在线视频| 视频区图区小说| 考比视频在线观看| 亚洲av成人一区二区三| 大香蕉久久成人网| 国产真人三级小视频在线观看| 国产99久久九九免费精品| 一级片免费观看大全| 亚洲少妇的诱惑av| 两个人看的免费小视频| 天堂中文最新版在线下载| 欧美亚洲日本最大视频资源| 午夜激情久久久久久久| 亚洲av电影在线进入| 亚洲精品第二区| 国产区一区二久久| 男女床上黄色一级片免费看| 女性被躁到高潮视频| 精品国产一区二区久久| 亚洲精品一二三| 免费日韩欧美在线观看| 在线观看免费午夜福利视频| 精品一品国产午夜福利视频| 精品少妇内射三级| 亚洲人成电影观看| 国产成人系列免费观看| 一区二区三区四区激情视频| 亚洲精品日韩在线中文字幕| 在线观看人妻少妇| 亚洲精品久久午夜乱码| 美女大奶头黄色视频| 精品亚洲乱码少妇综合久久| 高潮久久久久久久久久久不卡| 永久免费av网站大全| 极品少妇高潮喷水抽搐| 我的亚洲天堂| 亚洲av欧美aⅴ国产| 丁香六月天网| 久久青草综合色|