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

    基于多源遙感數(shù)據(jù)的森林生物量定量評價(jià)研究

    2016-06-04 03:14:10徐天蜀岳彩榮章皖秋張王菲西南林業(yè)大學(xué)林學(xué)院云南昆明650224
    關(guān)鍵詞:云南松

    徐天蜀 岳彩榮 章皖秋 張王菲 袁 華(西南林業(yè)大學(xué)林學(xué)院,云南昆明650224)

    ?

    基于多源遙感數(shù)據(jù)的森林生物量定量評價(jià)研究

    徐天蜀 岳彩榮 章皖秋 張王菲 袁 華
    (西南林業(yè)大學(xué)林學(xué)院,云南昆明650224)

    摘要:以ALOS PALSAR L波段雙極化FBD微波遙感數(shù)據(jù)及多光譜光學(xué)遙感數(shù)據(jù)AVNIR-2為基礎(chǔ),對數(shù)據(jù)進(jìn)行預(yù)處理,利用地面云南松林樣地坐標(biāo),提取HH、HV雙極化后向散射系數(shù)及極化比值3個(gè)因子,結(jié)合光學(xué)遙感數(shù)據(jù)提取4個(gè)波段值及NDVI、RVI 2個(gè)植被指數(shù),作為云南松林生物量估測因子。分別以微波數(shù)據(jù)、光學(xué)數(shù)據(jù)、微波及光學(xué)數(shù)據(jù)結(jié)合的多源遙感數(shù)據(jù),建立3個(gè)云南松林生物量估測模型。結(jié)果表明:所建模型經(jīng)方差分析均達(dá)到顯著相關(guān)或極顯著相關(guān)水平;PALSAR L波段雙極化后向散射系數(shù),可以反映森林生物量的變化,但反演精度有待進(jìn)一步提高;AVNIR-2數(shù)據(jù)模型優(yōu)于PALSAR L波段雙極化數(shù)據(jù)模型;多源數(shù)據(jù)模型與光學(xué)數(shù)據(jù)模型的估測精度相近。

    關(guān)鍵詞:森林生物量;微波數(shù)據(jù);光學(xué)數(shù)據(jù);多源遙感數(shù)據(jù);定量評價(jià);云南松

    遙感技術(shù)是較大尺度森林生物量定量評價(jià)及動(dòng)態(tài)監(jiān)測的重要工具,目前已有大量的研究和應(yīng)用,但所應(yīng)用的數(shù)據(jù)以光學(xué)遙感數(shù)據(jù)居多。由于可見光和近紅外光譜只與綠葉生物量產(chǎn)生反應(yīng),理論上只能獲取森林葉生物量的信息,無法反映森林冠層以下的樹干信息,在利用光譜信息估測森林生物量時(shí)實(shí)際是利用了森林生物量與綠葉生物量的相對生長關(guān)系達(dá)到間接估計(jì)森林生物量,然而這種對應(yīng)關(guān)系在森林結(jié)構(gòu)復(fù)雜或者估測對象包含多個(gè)樹種時(shí)會(huì)出現(xiàn)較大的誤差或錯(cuò)誤[1-2]。微波具有穿透樹冠的能力,主要與森林生物量的主體—枝和樹干發(fā)生作用,也能與樹葉發(fā)生作用,因而微波遙感技術(shù)為森林生物量全面和精確估測提供了可行的工具[2-3]。因此,本研究以ALOS PALSAR L波段雙極化微波遙感數(shù)據(jù)以及ALOS AVNIR多光譜數(shù)據(jù)為信息源,以云南松(Pinus yunnanensis)林為研究對象,探索多源數(shù)據(jù)結(jié)合應(yīng)用于森林生物量定量估測的方法和技術(shù)。

    1 研究區(qū)概況

    研究區(qū)位于云南省昆明市宜良縣境內(nèi),選在宜良縣花園林場。屬山地地貌,平均海拔1 700~1 900 m,中亞熱帶季風(fēng)氣候,平均氣溫16.3℃,夏無酷暑,冬無嚴(yán)寒,年平均氣溫16.3℃,年平均日照2 177.3 h,每年平均降雨量912.2 mm,年平均相對濕度75%,全年無霜期260 d左右。主要植被有云南松、華山松(Pinus armandi)、杉木(Cunninghamia lanceolata)、云南油杉(Keteleeria evelyniana)、桉樹(Eucalyptus robusta)等。

    2 數(shù)據(jù)來源

    2.1微波數(shù)據(jù)

    本研究的數(shù)據(jù)為2010年7月獲取ALOS PALSAR L波段的FBD Level 1.1雙極化數(shù)據(jù),包括水平極化HH、交叉極化HV。圖像中心點(diǎn)坐標(biāo):北緯24.76°,東經(jīng)103.26°;入射角34.3°,其中方位向分辨率為3.18 m,距離向分辨率為9.37 m,像幅18 432×4 640[4]。

    2.2多光譜數(shù)據(jù)

    采用2010年7月獲取的ALOS衛(wèi)星的新型可見光和近紅外輻射計(jì)AVNIR-2數(shù)據(jù),包括藍(lán)、綠、紅、近紅外4個(gè)波段,分辨率為10 m。

    2.3云南松林生物量數(shù)據(jù)

    分別在2010年6月及2011年7月進(jìn)行地面樣地調(diào)查,針對云南松純林及混交林,采用典型選樣法布設(shè)樣地。利用角規(guī)控制檢尺的方法調(diào)查樣地單位面積林木蓄積量,分別對不同樹種采用蓄積量-生物量轉(zhuǎn)換方程[4]計(jì)算單位面積生物量,同時(shí)記錄樣地中心點(diǎn)的GPS坐標(biāo),用于與遙感數(shù)據(jù)匹配。本研究共調(diào)查樣地100個(gè)。

    3 基礎(chǔ)數(shù)據(jù)處理

    3.1微波數(shù)據(jù)預(yù)處理

    微波遙感數(shù)據(jù)的成像方式為斜距成像,地物后向散射受地形影響較大,產(chǎn)生透視收縮、疊掩、陰影、頂部位移等幾何畸變,因此在應(yīng)用時(shí)均需要進(jìn)行數(shù)據(jù)預(yù)處理,本文以Gamma軟件為處理平臺(tái)進(jìn)行數(shù)據(jù)處理。

    1)輻射定標(biāo)。微波圖像的每個(gè)像元值都能表示成唯一確定的后向散射信號功率值,這一過程稱為輻射定標(biāo)。為了充分利用多通道多極化SAR數(shù)據(jù)信息,必須先對這些數(shù)據(jù)進(jìn)行輻射定標(biāo),才能應(yīng)用于地表散射特性的分析[5-6]。根據(jù)日本ALOS衛(wèi)星公司提供的定標(biāo)參數(shù),Level 1.1數(shù)據(jù)采用下列定標(biāo)公式:

    式中:X為像元讀數(shù)值;Y為后向散射系數(shù)(dB)。

    2)正射校正。微波圖像是根據(jù)多普勒原理進(jìn)行斜距成像的,SAR圖像的正射校正就是要將SAR圖像的距離-多普勒坐標(biāo)轉(zhuǎn)換為帶有地理編碼并糾正地形影響后的正射微波圖像。本研究利用帶有地圖坐標(biāo)的DEM,產(chǎn)生模擬SAR圖像,對真實(shí)SAR圖像進(jìn)行地理編碼及正射校正,獲取δ0圖像[7]。

    3)斑點(diǎn)噪聲抑制/濾波處理。相干斑點(diǎn)噪聲是微波影像固有的特征。當(dāng)微波照射在粗糙地表,返回的信號是由一個(gè)圖像分辨單元內(nèi)的多個(gè)散射體反射回波組成的,形成“斑點(diǎn)噪聲”,從而影響SAR圖像質(zhì)量[5-6]。常用的去噪方法有多視處理、濾波處理。本研究在正射校正時(shí)多視處理的基礎(chǔ)上,采用中值濾波進(jìn)行處理,濾波窗口為5×5。

    4)地形輻射校正。在不平坦地區(qū),地形對后向散射系數(shù)會(huì)產(chǎn)生影響,稱為地形引起的輻射畸變。利用局部入射角及校正后的像元面積進(jìn)行地形輻射校正,獲取γ0圖像。預(yù)處理后的圖像見圖1。

    3.2多光譜數(shù)據(jù)預(yù)處理

    以已經(jīng)做過幾何精校正的Quick-Bird影像為基準(zhǔn),均勻地選取了26個(gè)GCP點(diǎn),選擇二次多項(xiàng)式作為糾正方程,對AVNIR-2影像進(jìn)行幾何精校正,校正后的誤差控制在1個(gè)像元內(nèi)。

    圖1 預(yù)處理后的γ0L-hv圖像與地面樣地位置Fig.1 Preprocessed γ0L-hvimage and the location of sample plots

    4 森林生物量定量估測

    4.1微波遙感數(shù)據(jù)估測[8-17]

    1)遙感估測因子的提取。對經(jīng)過預(yù)處理的HH、HV 2種極化的γ0圖像,利用地面樣地GPS坐標(biāo)提取后向散射系數(shù),得、及極化比值。

    2)模型建立。以HH、HV后向散射系數(shù)及極化比值HV/ HH 3個(gè)因子為自變量,以樣地生物量的自然對數(shù)值為因變量,建立多元線性回歸模型,建模樣本為80個(gè),模型如下:

    3)方差分析及線性回歸關(guān)系顯著性檢驗(yàn)。根據(jù)方差分析(表1),遙感因子與森林生物量具有顯著相關(guān)關(guān)系,相關(guān)系數(shù)R=0.362 6,可用于森林生物量估測。

    表1 方差分析表Table 1 Variance analysis

    4.2多光譜光學(xué)遙感數(shù)據(jù)估測

    1)遙感估測因子的選擇。ALOS衛(wèi)星的AVNIR-2數(shù)據(jù)包括藍(lán)光B1、綠光B2、紅光B3、近紅外光B44個(gè)波段,在幾何校正的基礎(chǔ)上,以樣地GPS坐標(biāo)為準(zhǔn),提取4個(gè)波段的光譜值作為生物量估測因子。

    不同的光譜通道所獲得的植被信息有著各自的特點(diǎn),僅用個(gè)別波段或多個(gè)單波段數(shù)據(jù)來提取植被信息是相當(dāng)局限的,因而常常選用多光譜遙感數(shù)據(jù)經(jīng)線性或非線性組合,產(chǎn)生對植被長勢、生物量等有一定指示意義的植被指數(shù)加以分析。本研究利用4個(gè)光譜通道提取了歸一化植被指數(shù)NDVI及比值植被指數(shù)RVI 2個(gè)植被指數(shù)作為估測因子。

    2)模型的建立。以上述光譜值及植被指數(shù)共6個(gè)因子作為自變量,以樣地生物量作為因變量,建立多元線性回歸模型,模型如下:

    3)方差分析及線性回歸關(guān)系顯著性檢驗(yàn)。根據(jù)方差分析(表1),遙感因子與森林生物量具有極顯著相關(guān)關(guān)系,相關(guān)系數(shù)R=0.522 1,可用于森林生物量估測。

    4.3多源遙感數(shù)據(jù)結(jié)合的估測

    1)遙感因子的選擇。將微波數(shù)據(jù)與多光譜數(shù)據(jù)結(jié)合,作為森林生物量的估測因子,即:微波數(shù)據(jù)的、及極化比值3個(gè)因子;多光譜數(shù)據(jù)的B1、B2、B3、B44個(gè)光譜值;NDVI、RVI 2個(gè)植被指數(shù)共9個(gè)因子作為估測因子。

    2)模型的建立。以上述9個(gè)因子作為模型的自變量,以樣地生物量的對數(shù)值作為模型因變量,建立多元線性回歸模型。模型如下:

    3)方差分析及線性回歸關(guān)系顯著性檢驗(yàn)。根據(jù)方差分析(表1),遙感因子與森林生物量具有極顯著相關(guān)關(guān)系,相關(guān)系數(shù)R=0.536 6,可用于森林生物量估測。

    4.4精度檢驗(yàn)與評價(jià)

    1)精度檢驗(yàn)。利用隨機(jī)選取的20個(gè)檢驗(yàn)樣本及已建立的3個(gè)模型,分別計(jì)算樣地生物量的預(yù)測值。

    精度評價(jià)的指標(biāo)采用平均相對誤差(MRE)和均方根誤差(RMSE)。計(jì)算公式如下:

    3個(gè)模型預(yù)測值精度檢驗(yàn)見表2。

    表2 模型預(yù)測精度檢驗(yàn)表Table 2 Models prediction accuracy test

    2)精度評價(jià)。從相關(guān)系數(shù)分析,多源數(shù)據(jù)模型略優(yōu)于多光譜模型,多源數(shù)據(jù)模型及多光譜模型優(yōu)于微波數(shù)據(jù)模型。

    從平均相對誤差MRE分析,微波數(shù)據(jù)模型的平均相對誤差為36.8%,多光譜數(shù)據(jù)模型的平均相對誤差為29.1%,多源數(shù)據(jù)模型的平均相對誤差為29.9%??梢钥闯龆嘣磾?shù)據(jù)模型與多光譜數(shù)據(jù)模型優(yōu)于單純的微波數(shù)據(jù)模型,多源數(shù)據(jù)模型略優(yōu)于多光譜數(shù)據(jù)模型。

    從均方根誤差RMSE分析,多光譜數(shù)據(jù)模型為11.86,多源數(shù)據(jù)模型為13.21,微波數(shù)據(jù)模型為14.17,單純的微波數(shù)據(jù)模型的RMSE最大。

    5 結(jié)論與討論

    1)利用微波數(shù)據(jù)模型、多光譜數(shù)據(jù)模型、多源數(shù)據(jù)模型3種方法建立的森林生物量估測模型經(jīng)方差分析檢驗(yàn)均達(dá)到了顯著相關(guān)或極顯著相關(guān)水平,可用于森林生物量估測。

    2)微波雙極化數(shù)據(jù)的總后向散射系數(shù),可以反映森林生物量的變化,但直接利用總后向散射系數(shù)進(jìn)行反演的精度有待進(jìn)一步提高。

    3)AVNIR-2光學(xué)數(shù)據(jù)模型優(yōu)于PALSAR L波段雙極化微波數(shù)據(jù)模型。AVNIR-2光學(xué)模型與多源數(shù)據(jù)模型對生物量估測精度差異不大。

    4)微波后向散射系數(shù)不僅包括森林冠層的反射能量,還包括樹干、枝的多次體散射以及林下地面的表面散射。單純的利用總后向散射系數(shù)構(gòu)建生物量模型,不能確定各個(gè)散射成分對生物量的影響,會(huì)產(chǎn)生較大的誤差??蛇M(jìn)一步利用全極化數(shù)據(jù)進(jìn)行目標(biāo)分解,提取不同的散射分量,尋找與生物量相關(guān)性高的因子,挖掘微波數(shù)據(jù)應(yīng)用的潛力,提高估測精度。本研究采用多元線性回歸模型用于建模,而微波數(shù)據(jù)與多光譜數(shù)據(jù)與生物量的關(guān)系可能為非線性關(guān)系,可進(jìn)一步采用非線性模型擬合,提高估測精度。

    [參考文獻(xiàn)]

    [1]黎夏,葉嘉安,王樹功,等.紅樹林濕地植被生物量的雷達(dá)遙感估算[J].遙感學(xué)報(bào),2006,10(3):337-336.

    [2]陳爾學(xué).合成孔徑雷達(dá)森林生物量估測研究進(jìn)展[J].世界林業(yè)研究,1999,12(6):18-23.

    [3]Ranson K J,Sun G.Mapping biomass of a northern forest using multifrequency SAR data[J].Geoscience &Remote Sensing IEEE Transactions on,1994,32(2):388-396.

    [4]黃從德,張健,楊萬勤,等.四川省及重慶地區(qū)森林植被碳儲(chǔ)量動(dòng)態(tài)[J].生態(tài)學(xué)報(bào),2008,28(3):966-975.

    [5]尤紅建,付琨.合成孔徑雷達(dá)圖像精準(zhǔn)處理[M].北京:科學(xué)出版社,2011.

    [6]王超,張紅.全極化合成孔徑雷達(dá)圖像處理[M].北京:科學(xué)出版社,2008.

    [7]徐天蜀,岳彩榮,張王菲.山區(qū)星載ALOS PALSAR微波遙感圖像的正射校正[J].農(nóng)業(yè)工程學(xué)報(bào),2012,28(1):173-178.

    [8]楊永恬,李增元,陳爾學(xué),等.基于ALOS PALSAR數(shù)據(jù)的森林蓄積量估測技術(shù)研究[J].林業(yè)資源管理,2010(2):113-117.

    [9]廖靜娟,沈國狀.基于多極化SAR圖像的鄱陽湖濕地地表淹沒狀況動(dòng)態(tài)變化分析[J].遙感技術(shù)與應(yīng)用,2008,23(4):373-377.

    [10]張遠(yuǎn),張中浩,蘇世亮,等.基于微波冠層散射模型的水稻生物量遙感估算[J].農(nóng)業(yè)工程學(xué)報(bào),2011,27(9):100-105.

    [11]Gamma Remote Sensing AG.Land Application Tools-LAT.Documentation User′s Guide(Version 1.0)[R].Gümligen,2009.

    [12]Santoro M,F(xiàn)ransson J E S,Eriksson L E B,et al.Signatures of ALOS PALSAR L-band backscatter in Swedish forest[J].IEEE Transactions on Geoscience &Remote Sensing,2009,47(12):4001-4019.

    [13]Zhang Y,Huang H,Chen X,et al.Mapping paddy rice biomass using ALOS/ PALSAR imagery[C]/ / Education Technology and Training,2008.and 2008 International Workshop on Geoscience and Remote Sensing.ETT and GRS 2008.International Workshop on IEEE, 2008:207-210.

    [14]Mitchard E T A,Saatchi S S,Woodhouse I H,et al.Using satellite radar backscatter to predict above-ground woody biomass:a consistent relationship across four different African landscapes[J].Geophysical Research Letters,2009,36(23):195-215.

    [15]Kasischke E S,Christensen N L,Bourgeau-Chavez L L.Correlating radar backscatter with components of biomass in loblolly pine forests[J].IEEE Transactions on Geoscience&Remote Sensing,1995,33(3):643-659.

    [16]Rauste Y.Multi-temporal JERS SAR data in boreal forest biomass mapping[J].Remote Sensing of Environment,2005,97(2):263-275.

    [17]Kurvonen L,Pulliainen J,Hallikainen M.Retrieval of biomass in boreal forests from multitemporal ERS-1 and JERS-1 SAR images[J].Geoscience&Remote Sensing IEEE Transactions on,1999,37(1):198-205.

    (責(zé)任編輯曹龍)

    第1作者:徐天蜀(1964-),女,博士,教授。研究方向:森林資源管理及林業(yè)“3S”技術(shù)應(yīng)用研究。Email:tsxue64@163.com。

    Research on Quantitative Evaluation of Forest Biomass Based on Multi-Source Remote Sensing Data

    Xu Tianshu,Yue Cairong,Zhang Wanqiu,Zhang Wangfei,Yuan Hua
    (College of Forestry,Southwest Forestry University,Kumming Yunnan 650224,China)

    Abstract:Based on ALOS PALSAR L-band dual-polarization FBD microwave remote sensing data and multispectral optical remote sensing data AVNIR-2,preprocessing was conducted,then three microwave data factors including HH and HV dual-polarization backscattering coefficients and polarization ratio,and four original bands and two vegetation indices of NDVI,RVI which combined with optical remote sensing data were extracted with the coordinates of Pinus yunnanensis forest sample plots.After that,the three estimation models of Pinus yunnanensis forest biomass were built by using microwave data,optical data,and the multi-source remote sensing data combined with microwave and optical data.Research suggests that and the models reached the significant or extremely significant correlation level by variance analysis.PALSAR L-band dual-polarization backscattering coefficient could be reflected the changes of forest biomass,but inversion accuracy needs to be improved.AVNIR-2 data model was better than PALSAR L-band dual-polarization data model.The estimation accuracy of multi-source data model was similar to optical data model.

    Key words:forest biomass,microwave data,optical data,multi-source remote sensing data,quantitative evaluation,Pinus yunnanensis

    基金項(xiàng)目:國家自然科學(xué)基金資助項(xiàng)目(31260156、30960302)資助;西南林業(yè)大學(xué)林學(xué)一級學(xué)科資助;西南地區(qū)生物多樣性保育國家林業(yè)局重點(diǎn)實(shí)驗(yàn)室開放基金資助。

    收稿日期:2015-12-29

    doi:10.11929/ j.issn.2095-1914.2016.03.022

    中圖分類號:S771.8

    文獻(xiàn)標(biāo)志碼:A

    文章編號:2095-1914(2016)03-0126-05

    猜你喜歡
    云南松
    不同生長狀況對云南松產(chǎn)脂的影響
    云南松正交膠合木膠合工藝探索?
    昆明市西山林場云南松縱坑切梢小蠹危害調(diào)查及綜合防治
    基于干型差異的云南松子代苗木生物量研究
    種子(2019年5期)2019-07-02 11:37:44
    云南松幼苗對水分脅迫與干旱鍛煉的研究
    云南松天然更新幼苗和幼樹期的生長動(dòng)態(tài)分析
    云南松及其變種栽培與適應(yīng)性研究進(jìn)展
    云南云龍?zhí)斐刈匀槐Wo(hù)區(qū)云南松干形研究
    綠色科技(2018年9期)2018-01-29 06:01:52
    一年生云南松不同家系苗生物量研究
    一年生云南松不同家系苗生物量研究
    国产免费又黄又爽又色| 国产精品一二三区在线看| 视频区图区小说| 人妻 亚洲 视频| 亚洲欧美中文字幕日韩二区| 最黄视频免费看| 中国美白少妇内射xxxbb| 精品久久久久久久末码| 在线看a的网站| 99热这里只有精品一区| 我的女老师完整版在线观看| 亚洲欧美日韩东京热| 狂野欧美白嫩少妇大欣赏| 视频中文字幕在线观看| 人人妻人人爽人人添夜夜欢视频 | 成年av动漫网址| 婷婷色麻豆天堂久久| 久久婷婷青草| 国产美女午夜福利| 97在线视频观看| 国产亚洲5aaaaa淫片| 小蜜桃在线观看免费完整版高清| av福利片在线观看| a 毛片基地| 一级a做视频免费观看| 久久久精品94久久精品| 2021少妇久久久久久久久久久| 久久国产乱子免费精品| 免费大片黄手机在线观看| 久久精品国产鲁丝片午夜精品| av在线播放精品| 亚洲av.av天堂| 亚洲欧美一区二区三区黑人 | 赤兔流量卡办理| 日本av手机在线免费观看| 国内揄拍国产精品人妻在线| 精品酒店卫生间| 最近最新中文字幕大全电影3| 在现免费观看毛片| 亚州av有码| 一级毛片黄色毛片免费观看视频| 成人美女网站在线观看视频| 97精品久久久久久久久久精品| 欧美日韩视频高清一区二区三区二| 精品人妻视频免费看| 久久久久国产网址| 三级国产精品片| 久久午夜福利片| 亚洲av成人精品一二三区| 人妻一区二区av| 国产一级毛片在线| 成人无遮挡网站| 有码 亚洲区| 成年美女黄网站色视频大全免费 | 亚洲精品亚洲一区二区| 蜜桃在线观看..| 国产在线一区二区三区精| av不卡在线播放| 久久影院123| 高清日韩中文字幕在线| 国产精品久久久久久久电影| 人人妻人人看人人澡| 亚洲av在线观看美女高潮| 伦理电影免费视频| 免费看不卡的av| 男人狂女人下面高潮的视频| 日韩欧美精品免费久久| 国产精品一区二区在线不卡| 久久久久久久久大av| 亚洲激情五月婷婷啪啪| 国产欧美亚洲国产| 欧美97在线视频| a 毛片基地| 日本av手机在线免费观看| 高清黄色对白视频在线免费看 | 亚洲欧美日韩无卡精品| 涩涩av久久男人的天堂| 高清日韩中文字幕在线| 欧美区成人在线视频| 久久影院123| 交换朋友夫妻互换小说| 中文字幕制服av| 久久99热这里只有精品18| 91午夜精品亚洲一区二区三区| 久久久久久久大尺度免费视频| 各种免费的搞黄视频| 国产白丝娇喘喷水9色精品| 99久久人妻综合| 一级毛片我不卡| 精品久久国产蜜桃| 免费观看a级毛片全部| 精品国产露脸久久av麻豆| 亚洲国产毛片av蜜桃av| 亚洲国产欧美在线一区| 国产免费一级a男人的天堂| 国产免费又黄又爽又色| 五月天丁香电影| 日韩中文字幕视频在线看片 | 青春草亚洲视频在线观看| 亚洲色图av天堂| 成年av动漫网址| 精品少妇久久久久久888优播| 97超碰精品成人国产| 热99国产精品久久久久久7| 欧美日韩精品成人综合77777| 国产又色又爽无遮挡免| 欧美精品亚洲一区二区| 青春草亚洲视频在线观看| 亚洲美女视频黄频| 高清午夜精品一区二区三区| 三级国产精品片| 国产精品欧美亚洲77777| 伦理电影大哥的女人| 99久久人妻综合| 毛片女人毛片| 亚洲国产精品一区三区| 一区在线观看完整版| 亚洲欧美日韩东京热| 久久久久性生活片| 午夜福利视频精品| 自拍欧美九色日韩亚洲蝌蚪91 | 秋霞在线观看毛片| 中国三级夫妇交换| 婷婷色麻豆天堂久久| 日韩国内少妇激情av| 成年免费大片在线观看| 国产精品欧美亚洲77777| 1000部很黄的大片| 日韩一区二区视频免费看| 亚洲精品国产av成人精品| 在线观看一区二区三区| 九草在线视频观看| 99热这里只有精品一区| 男人和女人高潮做爰伦理| 久久久久视频综合| 国产男女超爽视频在线观看| 精品99又大又爽又粗少妇毛片| 国产免费一区二区三区四区乱码| 观看美女的网站| 国产精品一区二区三区四区免费观看| 在线免费十八禁| 夫妻午夜视频| 欧美97在线视频| 亚洲人与动物交配视频| av一本久久久久| 欧美日韩视频精品一区| 午夜视频国产福利| 水蜜桃什么品种好| 多毛熟女@视频| 成人黄色视频免费在线看| 免费观看在线日韩| 久久毛片免费看一区二区三区| 简卡轻食公司| 久久ye,这里只有精品| 婷婷色综合大香蕉| 久久亚洲国产成人精品v| 国产色爽女视频免费观看| 女性被躁到高潮视频| 国国产精品蜜臀av免费| 亚洲av男天堂| 中文欧美无线码| av在线蜜桃| 五月开心婷婷网| 涩涩av久久男人的天堂| 男人狂女人下面高潮的视频| 久久国内精品自在自线图片| 国产精品女同一区二区软件| 黑人猛操日本美女一级片| 各种免费的搞黄视频| 亚洲怡红院男人天堂| 国产午夜精品久久久久久一区二区三区| 日本黄大片高清| 国产成人精品福利久久| 亚洲性久久影院| av专区在线播放| 久久精品国产a三级三级三级| 免费看av在线观看网站| 狂野欧美激情性bbbbbb| 婷婷色综合www| 成人漫画全彩无遮挡| 激情五月婷婷亚洲| av一本久久久久| 九色成人免费人妻av| 中文精品一卡2卡3卡4更新| 久久久久久久久大av| 日韩大片免费观看网站| 国产精品福利在线免费观看| 天天躁夜夜躁狠狠久久av| 91狼人影院| 国产欧美日韩精品一区二区| 伊人久久国产一区二区| 久热久热在线精品观看| 男人狂女人下面高潮的视频| 女性被躁到高潮视频| 午夜老司机福利剧场| 在线观看国产h片| av免费在线看不卡| 我的老师免费观看完整版| 一级毛片黄色毛片免费观看视频| 国产伦理片在线播放av一区| 成人美女网站在线观看视频| 晚上一个人看的免费电影| 2021少妇久久久久久久久久久| 亚洲精品久久午夜乱码| 黑人猛操日本美女一级片| 日本猛色少妇xxxxx猛交久久| 日韩伦理黄色片| 91久久精品电影网| 国产大屁股一区二区在线视频| 免费看av在线观看网站| 大香蕉97超碰在线| 久热这里只有精品99| 青春草亚洲视频在线观看| 久久99热这里只有精品18| 亚洲av在线观看美女高潮| 18禁在线播放成人免费| 我的女老师完整版在线观看| 国产黄色免费在线视频| 亚洲欧美精品自产自拍| 岛国毛片在线播放| 成年女人在线观看亚洲视频| 热re99久久精品国产66热6| 大片电影免费在线观看免费| 久久久久精品久久久久真实原创| 久久97久久精品| 亚洲怡红院男人天堂| 男女啪啪激烈高潮av片| 大陆偷拍与自拍| 涩涩av久久男人的天堂| 国产深夜福利视频在线观看| 久久人人爽人人片av| 丝袜脚勾引网站| 久久午夜福利片| 久久韩国三级中文字幕| 舔av片在线| 免费av中文字幕在线| 一区二区av电影网| 国产精品秋霞免费鲁丝片| 十分钟在线观看高清视频www | 国产在线免费精品| 精品一区二区三卡| 一级毛片aaaaaa免费看小| 国产日韩欧美在线精品| 蜜桃亚洲精品一区二区三区| 国产精品一二三区在线看| 青春草国产在线视频| 亚洲国产精品成人久久小说| 亚洲国产精品成人久久小说| 亚洲综合精品二区| 黄色视频在线播放观看不卡| 男人爽女人下面视频在线观看| 大香蕉97超碰在线| 久久久久视频综合| 久久国产乱子免费精品| 国产 一区精品| 国产一区二区三区av在线| 91aial.com中文字幕在线观看| 亚洲精品久久久久久婷婷小说| 久久久亚洲精品成人影院| 日本-黄色视频高清免费观看| 身体一侧抽搐| 成人亚洲精品一区在线观看 | 日本一二三区视频观看| 久久6这里有精品| 少妇裸体淫交视频免费看高清| 国产精品女同一区二区软件| 自拍欧美九色日韩亚洲蝌蚪91 | 校园人妻丝袜中文字幕| 亚洲欧美日韩无卡精品| 黄色日韩在线| 免费大片黄手机在线观看| 精品酒店卫生间| 国产日韩欧美亚洲二区| 韩国av在线不卡| 日本欧美国产在线视频| 一区二区三区免费毛片| 久久精品久久精品一区二区三区| 亚洲av成人精品一二三区| 国产 一区精品| 五月伊人婷婷丁香| 成人免费观看视频高清| 欧美+日韩+精品| 成人一区二区视频在线观看| 亚洲不卡免费看| av国产精品久久久久影院| 久久久久视频综合| 亚洲自偷自拍三级| 久久久久精品久久久久真实原创| 国产精品成人在线| 亚洲天堂av无毛| 日本与韩国留学比较| 菩萨蛮人人尽说江南好唐韦庄| 日本午夜av视频| 免费观看在线日韩| 永久网站在线| 日韩人妻高清精品专区| 一区二区三区精品91| av在线观看视频网站免费| 免费人妻精品一区二区三区视频| 一区在线观看完整版| 精品久久久久久久久亚洲| 久久99精品国语久久久| 亚洲国产精品专区欧美| 亚洲精品亚洲一区二区| 1000部很黄的大片| 亚洲美女黄色视频免费看| 欧美bdsm另类| 欧美bdsm另类| 午夜福利在线观看免费完整高清在| 色综合色国产| 人妻夜夜爽99麻豆av| 亚洲av国产av综合av卡| 少妇被粗大猛烈的视频| av又黄又爽大尺度在线免费看| 爱豆传媒免费全集在线观看| 亚洲成人中文字幕在线播放| 大片电影免费在线观看免费| 一级毛片我不卡| 韩国av在线不卡| 国产精品99久久99久久久不卡 | 十八禁网站网址无遮挡 | 亚洲精品国产av成人精品| 久久精品久久久久久久性| 男女无遮挡免费网站观看| 日本与韩国留学比较| 精品久久久精品久久久| 国产精品不卡视频一区二区| 亚洲美女视频黄频| 亚洲欧美一区二区三区国产| 伦理电影大哥的女人| av一本久久久久| 伊人久久国产一区二区| 汤姆久久久久久久影院中文字幕| 男人和女人高潮做爰伦理| 十八禁网站网址无遮挡 | 国产一区二区三区av在线| 欧美日韩国产mv在线观看视频 | 大码成人一级视频| 人妻系列 视频| 国产有黄有色有爽视频| 国产精品女同一区二区软件| 99热这里只有是精品50| 久久国产乱子免费精品| 麻豆成人午夜福利视频| 岛国毛片在线播放| 一区二区av电影网| a级毛色黄片| 欧美日本视频| 成人国产av品久久久| 欧美日韩国产mv在线观看视频 | 久久精品国产自在天天线| 性色avwww在线观看| 欧美成人午夜免费资源| 99视频精品全部免费 在线| 久久人人爽av亚洲精品天堂 | 在线观看一区二区三区| 久久精品夜色国产| 色综合色国产| 18禁在线无遮挡免费观看视频| 国产免费又黄又爽又色| 草草在线视频免费看| 秋霞在线观看毛片| av网站免费在线观看视频| 国产av国产精品国产| 久久毛片免费看一区二区三区| 亚洲欧美日韩无卡精品| 国产v大片淫在线免费观看| 黑丝袜美女国产一区| 日韩av免费高清视频| 在线观看人妻少妇| 日本av免费视频播放| 国产精品久久久久久av不卡| 99国产精品免费福利视频| 亚洲国产欧美人成| 久久久久久久大尺度免费视频| 欧美日韩视频精品一区| 久久久久国产精品人妻一区二区| 日本-黄色视频高清免费观看| 国产精品女同一区二区软件| 80岁老熟妇乱子伦牲交| 五月玫瑰六月丁香| 精品99又大又爽又粗少妇毛片| 欧美精品人与动牲交sv欧美| 卡戴珊不雅视频在线播放| 久久这里有精品视频免费| 交换朋友夫妻互换小说| 最近的中文字幕免费完整| 熟女电影av网| av国产精品久久久久影院| 欧美97在线视频| 少妇人妻 视频| 成年免费大片在线观看| 欧美日韩综合久久久久久| 一级毛片电影观看| 欧美变态另类bdsm刘玥| 七月丁香在线播放| 国产精品久久久久久久久免| 只有这里有精品99| 国产美女午夜福利| 欧美xxⅹ黑人| 欧美xxxx黑人xx丫x性爽| 99国产精品免费福利视频| 黑人猛操日本美女一级片| 99热这里只有是精品50| www.色视频.com| 观看免费一级毛片| 国内揄拍国产精品人妻在线| 国产亚洲一区二区精品| 国产女主播在线喷水免费视频网站| 亚洲国产最新在线播放| av福利片在线观看| 美女内射精品一级片tv| 国产欧美亚洲国产| 久久久国产一区二区| 少妇 在线观看| 国产老妇伦熟女老妇高清| 国产中年淑女户外野战色| 成人午夜精彩视频在线观看| 你懂的网址亚洲精品在线观看| 成人国产av品久久久| 91在线精品国自产拍蜜月| 国产精品国产av在线观看| 国产精品国产三级国产av玫瑰| 亚洲人成网站在线播| 亚洲美女黄色视频免费看| 少妇人妻精品综合一区二区| 欧美区成人在线视频| 毛片一级片免费看久久久久| 91精品一卡2卡3卡4卡| 看非洲黑人一级黄片| 人妻系列 视频| 永久免费av网站大全| 国产精品久久久久久av不卡| 五月开心婷婷网| 91精品国产国语对白视频| 亚洲在久久综合| 国产精品无大码| 丰满少妇做爰视频| 中文在线观看免费www的网站| 亚洲真实伦在线观看| 多毛熟女@视频| 日韩免费高清中文字幕av| 一级黄片播放器| 色5月婷婷丁香| 久久久久精品久久久久真实原创| 亚洲综合精品二区| 赤兔流量卡办理| 成年av动漫网址| 亚洲中文av在线| 又粗又硬又长又爽又黄的视频| 啦啦啦啦在线视频资源| 国产毛片在线视频| 蜜臀久久99精品久久宅男| 夫妻午夜视频| 日韩一区二区三区影片| 亚洲第一av免费看| av.在线天堂| 男女边摸边吃奶| 在线观看国产h片| 91精品国产九色| 国产成人午夜福利电影在线观看| 人体艺术视频欧美日本| 精品午夜福利在线看| 99热这里只有是精品在线观看| 欧美高清成人免费视频www| 免费高清在线观看视频在线观看| 国内揄拍国产精品人妻在线| 在线观看免费高清a一片| 亚洲一级一片aⅴ在线观看| 亚洲va在线va天堂va国产| 韩国高清视频一区二区三区| 你懂的网址亚洲精品在线观看| 有码 亚洲区| 欧美 日韩 精品 国产| 国产白丝娇喘喷水9色精品| 国产精品国产三级国产av玫瑰| 成年人午夜在线观看视频| 亚洲欧美成人精品一区二区| 亚洲欧美一区二区三区黑人 | 一级爰片在线观看| 久久av网站| 精华霜和精华液先用哪个| 青青草视频在线视频观看| 亚洲精品国产av蜜桃| 亚洲,欧美,日韩| 亚洲欧美成人综合另类久久久| 少妇裸体淫交视频免费看高清| 性色av一级| 一个人看视频在线观看www免费| 久久久成人免费电影| 精品99又大又爽又粗少妇毛片| 亚洲在久久综合| 各种免费的搞黄视频| 亚洲人与动物交配视频| 欧美老熟妇乱子伦牲交| 久久久午夜欧美精品| 一级毛片电影观看| 下体分泌物呈黄色| 久久久欧美国产精品| 亚洲av欧美aⅴ国产| 成人18禁高潮啪啪吃奶动态图 | 美女中出高潮动态图| 交换朋友夫妻互换小说| 18禁裸乳无遮挡动漫免费视频| 天天躁日日操中文字幕| 观看免费一级毛片| 99久久精品一区二区三区| 国产永久视频网站| 丰满乱子伦码专区| 99久国产av精品国产电影| tube8黄色片| 美女脱内裤让男人舔精品视频| 各种免费的搞黄视频| 久久青草综合色| 国产欧美日韩一区二区三区在线 | 菩萨蛮人人尽说江南好唐韦庄| 美女福利国产在线 | 日本av免费视频播放| 欧美日韩国产mv在线观看视频 | 十分钟在线观看高清视频www | 国产成人精品婷婷| 纯流量卡能插随身wifi吗| 亚洲精品aⅴ在线观看| 天天躁日日操中文字幕| 深爱激情五月婷婷| 国产亚洲精品久久久com| 亚洲国产日韩一区二区| 国产国拍精品亚洲av在线观看| 国产欧美日韩一区二区三区在线 | 岛国毛片在线播放| 插逼视频在线观看| 久久久成人免费电影| av女优亚洲男人天堂| 亚洲人成网站在线观看播放| 国产黄频视频在线观看| 午夜免费鲁丝| 简卡轻食公司| 人人妻人人爽人人添夜夜欢视频 | 免费观看av网站的网址| 又粗又硬又长又爽又黄的视频| 欧美精品亚洲一区二区| 亚洲精品一二三| 自拍偷自拍亚洲精品老妇| 久久6这里有精品| 熟妇人妻不卡中文字幕| 欧美区成人在线视频| 亚洲精品色激情综合| 美女cb高潮喷水在线观看| 国产日韩欧美在线精品| 亚洲精品乱码久久久久久按摩| 最近手机中文字幕大全| 一个人看的www免费观看视频| 内地一区二区视频在线| 精品熟女少妇av免费看| 久久6这里有精品| 国产伦理片在线播放av一区| 国产男女内射视频| 国产色爽女视频免费观看| 国产伦精品一区二区三区视频9| 色视频www国产| 日韩av在线免费看完整版不卡| 一本一本综合久久| 久久国产精品大桥未久av | 国产亚洲91精品色在线| 新久久久久国产一级毛片| 精品国产乱码久久久久久小说| 亚洲精品第二区| 日韩人妻高清精品专区| 国产精品熟女久久久久浪| 在线天堂最新版资源| 天天躁日日操中文字幕| 国产探花极品一区二区| 久久综合国产亚洲精品| 国产精品人妻久久久久久| 色哟哟·www| 欧美97在线视频| 国产精品人妻久久久影院| 精品亚洲成a人片在线观看 | 性色av一级| 国产亚洲5aaaaa淫片| 久久久久国产精品人妻一区二区| 纵有疾风起免费观看全集完整版| 中文字幕精品免费在线观看视频 | 精品酒店卫生间| 色视频在线一区二区三区| 欧美变态另类bdsm刘玥| 国产亚洲91精品色在线| 免费少妇av软件| av又黄又爽大尺度在线免费看| 免费黄频网站在线观看国产| 国产成人91sexporn| 日日啪夜夜爽| 亚洲av男天堂| 精品人妻一区二区三区麻豆| 亚洲怡红院男人天堂| av线在线观看网站| 国产伦精品一区二区三区视频9| 亚洲人成网站在线播| 26uuu在线亚洲综合色| 精品99又大又爽又粗少妇毛片| 日韩中字成人| 美女视频免费永久观看网站| 精华霜和精华液先用哪个| 欧美 日韩 精品 国产| 国产一区二区在线观看日韩| 亚洲av不卡在线观看| 亚洲成人中文字幕在线播放| 如何舔出高潮| 又大又黄又爽视频免费| 国产深夜福利视频在线观看| 九九久久精品国产亚洲av麻豆| 直男gayav资源| 我的女老师完整版在线观看| 老女人水多毛片| 亚洲国产欧美在线一区|