• <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
    一年生云南松不同家系苗生物量研究
    一年生云南松不同家系苗生物量研究
    少妇人妻 视频| 成人影院久久| 成人毛片a级毛片在线播放| 伊人久久精品亚洲午夜| 一个人看视频在线观看www免费| 在线播放无遮挡| 久久久久久人妻| 国产v大片淫在线免费观看| 国产精品久久久久久精品古装| 久久人人爽人人片av| videossex国产| 日本欧美国产在线视频| 精品久久久久久久末码| 在线亚洲精品国产二区图片欧美 | 亚洲精品456在线播放app| 久热这里只有精品99| 欧美三级亚洲精品| 亚洲精品乱久久久久久| 97热精品久久久久久| 久久久久久久大尺度免费视频| 久久久精品免费免费高清| 成年美女黄网站色视频大全免费 | 22中文网久久字幕| 日本猛色少妇xxxxx猛交久久| 性色avwww在线观看| 女的被弄到高潮叫床怎么办| 亚洲成色77777| 一区在线观看完整版| 在线观看国产h片| 精品人妻熟女av久视频| 少妇的逼水好多| 亚洲精品aⅴ在线观看| 99久久中文字幕三级久久日本| 国产免费视频播放在线视频| videos熟女内射| 成人免费观看视频高清| 日本wwww免费看| 久久6这里有精品| 婷婷色综合www| 青青草视频在线视频观看| 精品亚洲成a人片在线观看 | 亚洲国产精品成人久久小说| 日韩av免费高清视频| 2018国产大陆天天弄谢| 2018国产大陆天天弄谢| 激情 狠狠 欧美| 精品少妇久久久久久888优播| 久久久久人妻精品一区果冻| 麻豆精品久久久久久蜜桃| 天堂俺去俺来也www色官网| 特大巨黑吊av在线直播| 国产精品久久久久久精品电影小说 | 男女下面进入的视频免费午夜| 一本—道久久a久久精品蜜桃钙片| 永久网站在线| 免费av不卡在线播放| 亚洲欧美精品自产自拍| 国产av码专区亚洲av| 免费观看av网站的网址| 日日啪夜夜撸| 亚洲欧美成人综合另类久久久| 亚洲精品色激情综合| 男女边吃奶边做爰视频| 啦啦啦中文免费视频观看日本| 97热精品久久久久久| 国产精品久久久久久精品电影小说 | 蜜桃久久精品国产亚洲av| 亚洲人成网站高清观看| 国产精品国产av在线观看| 免费大片黄手机在线观看| 男人舔奶头视频| 99久久中文字幕三级久久日本| 亚洲va在线va天堂va国产| 精品亚洲乱码少妇综合久久| 亚洲,欧美,日韩| h视频一区二区三区| 亚洲国产欧美在线一区| av免费在线看不卡| 日本av手机在线免费观看| 男人和女人高潮做爰伦理| 卡戴珊不雅视频在线播放| 久久久久网色| 久热久热在线精品观看| h视频一区二区三区| 亚洲av二区三区四区| av.在线天堂| 亚洲国产高清在线一区二区三| 中文字幕av成人在线电影| 精华霜和精华液先用哪个| 亚洲av福利一区| 一级片'在线观看视频| 亚洲精品自拍成人| 一区二区av电影网| 免费不卡的大黄色大毛片视频在线观看| 国产一区亚洲一区在线观看| 久久综合国产亚洲精品| 国产片特级美女逼逼视频| 熟妇人妻不卡中文字幕| 乱码一卡2卡4卡精品| 日韩在线高清观看一区二区三区| 欧美高清成人免费视频www| 伦理电影大哥的女人| 国产免费一级a男人的天堂| 精品人妻视频免费看| 国产精品国产三级国产av玫瑰| 精品久久国产蜜桃| 卡戴珊不雅视频在线播放| 99热6这里只有精品| 乱码一卡2卡4卡精品| 人妻 亚洲 视频| 大香蕉久久网| 两个人的视频大全免费| 亚洲欧美精品专区久久| 国产精品av视频在线免费观看| 高清日韩中文字幕在线| 99热网站在线观看| av免费观看日本| 精品人妻视频免费看| 男女下面进入的视频免费午夜| 亚洲精品日韩av片在线观看| 成人毛片a级毛片在线播放| 边亲边吃奶的免费视频| 永久免费av网站大全| av.在线天堂| 国产伦在线观看视频一区| 久久精品国产亚洲av涩爱| 99re6热这里在线精品视频| 精品国产露脸久久av麻豆| 亚洲美女黄色视频免费看| 爱豆传媒免费全集在线观看| 嫩草影院入口| 免费在线观看成人毛片| a级毛色黄片| 免费不卡的大黄色大毛片视频在线观看| 亚洲av男天堂| 多毛熟女@视频| 少妇丰满av| 亚洲综合色惰| a级毛片免费高清观看在线播放| 成人漫画全彩无遮挡| 国产男女内射视频| 99九九线精品视频在线观看视频| 97在线视频观看| 美女cb高潮喷水在线观看| 亚洲精品456在线播放app| 赤兔流量卡办理| 一级毛片黄色毛片免费观看视频| 国产一区二区在线观看日韩| 91久久精品国产一区二区三区| 99久久综合免费| 九色成人免费人妻av| 日本-黄色视频高清免费观看| 亚洲精品自拍成人| 亚洲伊人久久精品综合| 国产精品成人在线| 亚洲国产高清在线一区二区三| 蜜臀久久99精品久久宅男| 欧美精品一区二区免费开放| 少妇人妻一区二区三区视频| 一级a做视频免费观看| 春色校园在线视频观看| 女的被弄到高潮叫床怎么办| 97超碰精品成人国产| 国产免费一级a男人的天堂| 国产成人91sexporn| 国精品久久久久久国模美| 亚洲精品第二区| 一级毛片aaaaaa免费看小| 女人久久www免费人成看片| 只有这里有精品99| 欧美国产精品一级二级三级 | 国产亚洲91精品色在线| 免费黄色在线免费观看| 大陆偷拍与自拍| 在线观看三级黄色| 嘟嘟电影网在线观看| 国产av码专区亚洲av| 国产亚洲午夜精品一区二区久久| 久久久久网色| 最近的中文字幕免费完整| 只有这里有精品99| 国产一区二区在线观看日韩| 国产男人的电影天堂91| 卡戴珊不雅视频在线播放| 最近的中文字幕免费完整| 久久久欧美国产精品| 又大又黄又爽视频免费| 国产片特级美女逼逼视频| 日韩欧美一区视频在线观看 | a级毛色黄片| 亚洲av国产av综合av卡| 欧美亚洲 丝袜 人妻 在线| www.色视频.com| 91在线精品国自产拍蜜月| 如何舔出高潮| 亚洲性久久影院| 日韩一本色道免费dvd| 免费av中文字幕在线| 国产精品爽爽va在线观看网站| 久久人人爽人人爽人人片va| 日韩中文字幕视频在线看片 | 亚洲人与动物交配视频| 超碰av人人做人人爽久久| 成人影院久久| 青春草亚洲视频在线观看| 97超视频在线观看视频| 久久久久久伊人网av| 日日撸夜夜添| 亚州av有码| 啦啦啦在线观看免费高清www| 有码 亚洲区| 国产高清国产精品国产三级 | 成人亚洲精品一区在线观看 | 观看美女的网站| 久久精品夜色国产| 国产成人午夜福利电影在线观看| 亚洲精品视频女| 国产精品人妻久久久影院| 日日摸夜夜添夜夜爱| 国产精品蜜桃在线观看| 午夜福利高清视频| 一级毛片我不卡| 久久久欧美国产精品| 看非洲黑人一级黄片| 日本一二三区视频观看| 人人妻人人看人人澡| 午夜福利网站1000一区二区三区| 搡女人真爽免费视频火全软件| 国精品久久久久久国模美| tube8黄色片| 1000部很黄的大片| 亚洲内射少妇av| 男人狂女人下面高潮的视频| 日韩欧美 国产精品| 啦啦啦中文免费视频观看日本| 国产成人精品一,二区| 亚洲真实伦在线观看| 国语对白做爰xxxⅹ性视频网站| 80岁老熟妇乱子伦牲交| 日日撸夜夜添| 纯流量卡能插随身wifi吗| 男的添女的下面高潮视频| 日本av手机在线免费观看| 亚洲第一av免费看| 欧美日韩精品成人综合77777| 老司机影院成人| 国产成人一区二区在线| 久久久久视频综合| 极品教师在线视频| 国产极品天堂在线| 最近手机中文字幕大全| 欧美老熟妇乱子伦牲交| 在线观看三级黄色| 18禁动态无遮挡网站| 女性生殖器流出的白浆| 日韩欧美一区视频在线观看 | 亚洲欧美精品自产自拍| 一级毛片电影观看| 中文字幕免费在线视频6| 国产av码专区亚洲av| 免费av中文字幕在线| 在线亚洲精品国产二区图片欧美 | 内地一区二区视频在线| 两个人的视频大全免费| 日韩成人av中文字幕在线观看| 精品亚洲乱码少妇综合久久| 91久久精品电影网| 精品一区在线观看国产| 亚洲精品成人av观看孕妇| 一级毛片aaaaaa免费看小| 超碰97精品在线观看| 尾随美女入室| 熟女电影av网| 日韩精品有码人妻一区| 亚洲aⅴ乱码一区二区在线播放| 国产欧美日韩一区二区三区在线 | 嫩草影院入口| 国产大屁股一区二区在线视频| 日本-黄色视频高清免费观看| 高清av免费在线| 欧美精品亚洲一区二区| 91精品一卡2卡3卡4卡| 国产午夜精品久久久久久一区二区三区| 激情 狠狠 欧美| 欧美3d第一页| 看免费成人av毛片| 免费播放大片免费观看视频在线观看| 国精品久久久久久国模美| 亚洲精品aⅴ在线观看| 久久99精品国语久久久| 国产真实伦视频高清在线观看| 大码成人一级视频| 啦啦啦视频在线资源免费观看| 超碰97精品在线观看| 性高湖久久久久久久久免费观看| 精品久久久久久电影网| 一区二区三区四区激情视频| 精品一区二区三区视频在线| 久久婷婷青草| 亚洲av男天堂| 久久久久久久久久久丰满| 日日啪夜夜撸| 久久久久国产精品人妻一区二区| 国产精品久久久久久久电影| 亚洲aⅴ乱码一区二区在线播放| 日韩成人伦理影院| 2018国产大陆天天弄谢| 日韩视频在线欧美| 啦啦啦在线观看免费高清www| 亚洲精品乱码久久久v下载方式| 一区二区三区精品91| 国产精品一区www在线观看| 国产有黄有色有爽视频| 国产精品久久久久久av不卡| 成人高潮视频无遮挡免费网站| 舔av片在线| av专区在线播放| 久久精品久久久久久噜噜老黄| 性色av一级| 欧美精品人与动牲交sv欧美| 国内精品宾馆在线| 一级毛片黄色毛片免费观看视频| 插阴视频在线观看视频| 国产av精品麻豆| 久久97久久精品| 尾随美女入室| 偷拍熟女少妇极品色| 国产又色又爽无遮挡免| 麻豆国产97在线/欧美| 国产一区有黄有色的免费视频| 国产成人91sexporn| 午夜精品国产一区二区电影| 少妇被粗大猛烈的视频| 熟妇人妻不卡中文字幕| 亚洲欧洲国产日韩| 欧美极品一区二区三区四区| 国内揄拍国产精品人妻在线| 国产v大片淫在线免费观看| 久久久午夜欧美精品| 免费在线观看成人毛片| 免费高清在线观看视频在线观看| 99热这里只有是精品在线观看| 中文字幕av成人在线电影| 亚洲精品日本国产第一区| 精品国产一区二区三区久久久樱花 | 校园人妻丝袜中文字幕| 国产成人freesex在线| 天堂8中文在线网| 久久亚洲国产成人精品v| 色哟哟·www| 国产淫语在线视频| 免费av中文字幕在线| 国产欧美日韩精品一区二区| 激情五月婷婷亚洲| 日日摸夜夜添夜夜添av毛片| 777米奇影视久久| 极品少妇高潮喷水抽搐| 亚洲精品乱久久久久久| 成人综合一区亚洲| 夫妻午夜视频| 97超碰精品成人国产| 国产视频首页在线观看| 久久久欧美国产精品| 精品一区在线观看国产| 尾随美女入室| 久久青草综合色| 免费看光身美女| 日本爱情动作片www.在线观看| 亚洲久久久国产精品| 国产综合精华液| 欧美日本视频| 免费大片18禁| 九草在线视频观看| 久久午夜福利片| 成人亚洲欧美一区二区av| 国模一区二区三区四区视频| 欧美最新免费一区二区三区| 精品久久久久久久末码| 亚洲成色77777| 久久午夜福利片| av天堂中文字幕网| 毛片女人毛片| 国产永久视频网站| 国产精品一二三区在线看| 国产精品福利在线免费观看| 特大巨黑吊av在线直播| 草草在线视频免费看| 老熟女久久久| 十分钟在线观看高清视频www | 国产色婷婷99| 国产精品免费大片| 国产又色又爽无遮挡免| 一级毛片久久久久久久久女| 黄色一级大片看看| 97热精品久久久久久| 美女脱内裤让男人舔精品视频| 丰满少妇做爰视频| 国产成人精品久久久久久| 在现免费观看毛片| 欧美三级亚洲精品| 2021少妇久久久久久久久久久| 国产免费福利视频在线观看| 国产精品99久久99久久久不卡 | kizo精华| 久久99热这里只频精品6学生| 日韩,欧美,国产一区二区三区| 一区二区三区四区激情视频| 色综合色国产| 能在线免费看毛片的网站| 久久久久久久久久人人人人人人| 一区二区三区乱码不卡18| 色网站视频免费| 日韩免费高清中文字幕av| 精品熟女少妇av免费看| 99久国产av精品国产电影| av.在线天堂| 亚洲综合精品二区| av在线蜜桃| 日韩 亚洲 欧美在线| freevideosex欧美| 欧美成人精品欧美一级黄| 亚洲欧美成人精品一区二区| 国产亚洲91精品色在线| tube8黄色片| 十分钟在线观看高清视频www | 99久久精品一区二区三区| 狂野欧美激情性xxxx在线观看| 国产又色又爽无遮挡免| 嫩草影院新地址| 九草在线视频观看| 最后的刺客免费高清国语| 国产亚洲5aaaaa淫片| 蜜桃亚洲精品一区二区三区| 少妇被粗大猛烈的视频| 老师上课跳d突然被开到最大视频| 国产欧美另类精品又又久久亚洲欧美| 精品一品国产午夜福利视频| 18禁在线播放成人免费| 国产欧美亚洲国产| 久久午夜福利片| 人人妻人人澡人人爽人人夜夜| 精品久久久精品久久久| 身体一侧抽搐| 在线免费观看不下载黄p国产| 亚洲欧美精品专区久久| 亚洲不卡免费看| 国产一区亚洲一区在线观看| 国产精品国产av在线观看| 中国国产av一级| 人人妻人人添人人爽欧美一区卜 | 国产精品一区二区性色av| 另类亚洲欧美激情| 搡女人真爽免费视频火全软件| 男人狂女人下面高潮的视频| 亚洲精品视频女| 老司机影院成人| 91午夜精品亚洲一区二区三区| 少妇熟女欧美另类| 久久精品国产自在天天线| 涩涩av久久男人的天堂| 多毛熟女@视频| 人妻一区二区av| 一二三四中文在线观看免费高清| 日本免费在线观看一区| 亚洲人与动物交配视频| 五月伊人婷婷丁香| 老师上课跳d突然被开到最大视频| 一级片'在线观看视频| 日韩在线高清观看一区二区三区| 久久6这里有精品| 插逼视频在线观看| 狂野欧美激情性bbbbbb| 国产真实伦视频高清在线观看| 国产精品久久久久久久久免| 亚洲av.av天堂| 精品一品国产午夜福利视频| 在线观看美女被高潮喷水网站| 亚洲中文av在线| av不卡在线播放| 亚洲一区二区三区欧美精品| 18禁动态无遮挡网站| 日韩大片免费观看网站| 国产又色又爽无遮挡免| 国产成人免费无遮挡视频| 久久女婷五月综合色啪小说| 一二三四中文在线观看免费高清| 精品99又大又爽又粗少妇毛片| 中文字幕人妻熟人妻熟丝袜美| 边亲边吃奶的免费视频| 国产成人a∨麻豆精品| 日日摸夜夜添夜夜添av毛片| av又黄又爽大尺度在线免费看| 欧美成人午夜免费资源| 国产精品福利在线免费观看| 国产精品久久久久成人av| 国产黄片美女视频| 免费观看性生交大片5| 免费看不卡的av| 国产国拍精品亚洲av在线观看| 久久久a久久爽久久v久久| 自拍偷自拍亚洲精品老妇| 久久久久网色| 三级经典国产精品| 91在线精品国自产拍蜜月| 91精品一卡2卡3卡4卡| 高清av免费在线| 久久久久久久亚洲中文字幕| 日韩制服骚丝袜av| 亚洲av不卡在线观看| 欧美日韩在线观看h| 久久久久精品性色| 大话2 男鬼变身卡| 国产亚洲精品久久久com| av天堂中文字幕网| 国产成人a∨麻豆精品| www.av在线官网国产| 亚洲aⅴ乱码一区二区在线播放| 亚洲最大成人中文| 久久6这里有精品| 伦理电影大哥的女人| 免费观看性生交大片5| 内射极品少妇av片p| 成人国产av品久久久| 汤姆久久久久久久影院中文字幕| 美女内射精品一级片tv| 国产欧美亚洲国产| 精品少妇久久久久久888优播| 内射极品少妇av片p| 亚洲精品国产成人久久av| 国产欧美日韩精品一区二区| 欧美另类一区| 久久这里有精品视频免费| 亚洲av福利一区| 亚洲av成人精品一区久久| 久久久久久伊人网av| 欧美高清性xxxxhd video| 黄片wwwwww| 男人狂女人下面高潮的视频| 亚洲精品日韩av片在线观看| 欧美成人精品欧美一级黄| 街头女战士在线观看网站| 少妇 在线观看| 久久久午夜欧美精品| 国产精品熟女久久久久浪| 久久久久久久久久久免费av| 91精品伊人久久大香线蕉| 日本黄色日本黄色录像| 毛片一级片免费看久久久久| 欧美区成人在线视频| 热99国产精品久久久久久7| 大话2 男鬼变身卡| 99久久精品国产国产毛片| 全区人妻精品视频| 男男h啪啪无遮挡| 黄色配什么色好看| 黄色怎么调成土黄色| 十八禁网站网址无遮挡 | 看非洲黑人一级黄片| 欧美极品一区二区三区四区| 久久久久精品久久久久真实原创| 蜜桃亚洲精品一区二区三区| 日本与韩国留学比较| 久久精品夜色国产| 22中文网久久字幕| 国产人妻一区二区三区在| 亚洲av日韩在线播放| 日本一二三区视频观看| 日本免费在线观看一区| 国产精品嫩草影院av在线观看| 中文天堂在线官网| 丝袜脚勾引网站| 国产乱人视频| 国产精品一区www在线观看| 午夜日本视频在线| 久久精品夜色国产| 免费高清在线观看视频在线观看| 国产乱来视频区| 哪个播放器可以免费观看大片| 99热这里只有精品一区| 黑人高潮一二区| 国产伦精品一区二区三区四那| 97在线视频观看| 国产午夜精品一二区理论片| 舔av片在线| 久久久久国产网址| 嫩草影院入口| av天堂中文字幕网| 老司机影院成人| 日日啪夜夜撸| 香蕉精品网在线| 亚洲伊人久久精品综合| 97在线人人人人妻| 国产精品人妻久久久影院| 午夜福利网站1000一区二区三区| 国产美女午夜福利| 亚洲精品,欧美精品| 超碰97精品在线观看| 亚洲欧洲国产日韩| 熟妇人妻不卡中文字幕| 丰满人妻一区二区三区视频av| 99久久人妻综合| 亚洲av福利一区| 欧美老熟妇乱子伦牲交| 自拍偷自拍亚洲精品老妇| 一区二区三区免费毛片| 熟女电影av网| 亚洲精品日本国产第一区| 亚洲精品亚洲一区二区| 亚洲精品乱码久久久v下载方式| 少妇的逼好多水| 97超碰精品成人国产| 好男人视频免费观看在线| 身体一侧抽搐| 国产精品伦人一区二区| 大码成人一级视频|