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

    基于無人機(jī)影像的天山云杉林單木樹冠信息提取

    2020-09-21 09:30:02楊勇強(qiáng)王振錫師玉霞高亞利
    新疆農(nóng)業(yè)科學(xué) 2020年8期
    關(guān)鍵詞:單木分水嶺郁閉度

    楊勇強(qiáng),王振錫,師玉霞,連 玲,高亞利

    (新疆農(nóng)業(yè)大學(xué)林學(xué)與園藝學(xué)院,/新疆教育廳干旱區(qū)林業(yè)生態(tài)與產(chǎn)業(yè)技術(shù)重點(diǎn)實(shí)驗(yàn)室, 烏魯木齊 830052)

    0 引 言

    【研究意義】樹冠是單株林木主干以上部位的總稱,是樹木與外界進(jìn)行營(yíng)養(yǎng)交換的主要場(chǎng)所[1],其生長(zhǎng)的優(yōu)劣程度直接影響森林生態(tài)系統(tǒng)對(duì)地上部分資源的利用能力,與林木生長(zhǎng)有著密切的關(guān)系[2]。天山云杉(PiceaSchrenkianavartianshanica)是新疆山地森林的重要組成部分,主要分布在天山北麓,對(duì)新疆的生態(tài)平衡、物種多樣性和水源涵養(yǎng)有著至關(guān)重要的作用[3-4]。而對(duì)采伐后天山云杉單株樹冠的提取可以估測(cè)林分密度并推測(cè)地上生物量,樹冠信息的精確提取是構(gòu)建林業(yè)遙感定量監(jiān)測(cè)平臺(tái)的基礎(chǔ)[5]。【前人研究進(jìn)展】早期的遙感影像分辨率低,單個(gè)像素的大小比樹冠大,無法準(zhǔn)確識(shí)別單木樹冠信息[6]。但隨著遙感技術(shù)的發(fā)展及影像分辨率的提高,使得提取單木樹冠信息成為可能。目前,獲取單株樹冠信息的方法很多,如谷地跟蹤法[7-9]、區(qū)域增長(zhǎng)法[10-11]、分水嶺分割法[12-13]和局部射線法[14]等。隨著科學(xué)技術(shù)的發(fā)展使數(shù)據(jù)獲取變得多樣化,研究者進(jìn)而將激光雷達(dá)、高分辨率雷達(dá)[15-16]以及點(diǎn)云數(shù)據(jù)[17]融入單木樹冠提取。其中分水嶺分割方法應(yīng)用范圍較廣,然而傳統(tǒng)的分水嶺分割提取樹冠易造成過分割現(xiàn)象。很多研究在傳統(tǒng)分水嶺分割方法上進(jìn)行改進(jìn),郭昱杉等[18]將樹冠頂點(diǎn)作為標(biāo)記,應(yīng)用標(biāo)記控制分水嶺分割進(jìn)行樹冠信息提取,提取精度有所提升;滕文秀等[19]提出迭代H-minima改進(jìn)分水嶺分割方法提取樹冠信息,比標(biāo)記控制分水嶺分割算法F測(cè)度提高31.99%;于旭宅等[20]基于NDVI植被指數(shù)對(duì)樹冠信息提取并有效剔除背景及輸電線路等的影響,樣本精度達(dá)到88.3%。【本研究切入點(diǎn)】對(duì)具有復(fù)雜立地條件的山區(qū)森林樹冠提取的研究較少。標(biāo)記控制分水嶺分割方法的關(guān)鍵在于樹冠頂點(diǎn)和樹冠輪廓的正確識(shí)別。天山云杉樹冠通常是中心高、周圍低的近似圓球形,在無人機(jī)遙感影像上表現(xiàn)為樹冠頂點(diǎn)出較亮(光譜值較高),四周較暗(光譜之較低)的特征,這對(duì)利用局部最大值法探測(cè)樹冠頂點(diǎn)有利。天山云杉林的年齡和林分密度的不同,會(huì)造成樹冠大小不一、樹冠相互連接及遮擋的現(xiàn)象,對(duì)樹冠輪廓識(shí)別不利,導(dǎo)致過分割現(xiàn)象。研究傳統(tǒng)分水嶺分割中過分割現(xiàn)象,實(shí)現(xiàn)天山云杉單木樹冠自動(dòng)提取?!緮M解決的關(guān)鍵問題】在傳統(tǒng)分水嶺分割方法的基礎(chǔ)上進(jìn)行優(yōu)化:基于無人機(jī)遙感影像,利用高斯-拉普拉斯算子(Laplacian of Gaussian,LOG)結(jié)合最大類間方差尋找最優(yōu)閾值對(duì)影像進(jìn)行處理,并使用標(biāo)記控制分水嶺方法對(duì)樹冠信息進(jìn)行提取,為天山云杉林的伐后更新及恢復(fù)評(píng)價(jià)提供參考。

    1 材料與方法

    1.1 材 料

    研究區(qū)位于天山北麓中段的新疆農(nóng)業(yè)大學(xué)實(shí)習(xí)林場(chǎng)(N43°16′~44°07′,E86°46′~87°56′),地處頭屯河上游,北臨準(zhǔn)噶爾盆地,平均海拔約2 200 m,林場(chǎng)內(nèi)年平均溫度3℃,年降水量500~600 mm,年降水分布不均,主要體現(xiàn)在季節(jié)不同上,春夏偏多,2009年調(diào)查表明,總面積10 046.91 hm2,林地面積6 269.06 hm2,樣地主要分布在東北和西北坡,以天山云杉為主,林下土壤為普通灰褐色森林土。圖1

    2017年8月,在新疆農(nóng)業(yè)大學(xué)實(shí)習(xí)林場(chǎng)進(jìn)行野外典型樣地調(diào)查,根據(jù)新疆農(nóng)業(yè)大學(xué)實(shí)習(xí)林場(chǎng)二類調(diào)查結(jié)果,按照低(0.20 ~ 0.39)中(0.40 ~ 0.69)高(0.7以上)郁閉度對(duì)主伐跡地進(jìn)行挑選,每組挑選3塊典型樣地進(jìn)行研究,測(cè)量并記錄樣方內(nèi)所有天山云杉的樹高、胸徑和郁閉度等信息。其中林分郁閉度采用Win SCANOPY 2016冠層分析儀獲取。并用亞米級(jí)GPS記錄樣地中心點(diǎn)和四個(gè)角的位置信息,以便將樣地位置信息準(zhǔn)確的對(duì)應(yīng)至無人機(jī)遙感影像中。

    研究區(qū)影像數(shù)據(jù)獲取時(shí)間為2017年冬季,無人機(jī)進(jìn)行航拍,航拍地點(diǎn)為新疆農(nóng)業(yè)大學(xué)實(shí)習(xí)林場(chǎng)主伐跡地區(qū)域,覆蓋面積大約36.2 km2,所使用的無人機(jī)型號(hào)為UNIT-AIR-FW3,航拍高度為923 m,搭載相機(jī)型號(hào)為DSC-RX1RM2(35 mm),像素大小是4.53 μm×4.53 μm,拍攝影像尺寸為7 952×5 304,等效焦距為35 mm,空間分辨率為0.118 m,共拍攝1 200幅影像。對(duì)無人機(jī)影像的預(yù)處理包括:影像質(zhì)量檢測(cè),特征點(diǎn)提取,生成數(shù)字高程模型(digital elevation model,DEM),圖像拼接等操作。

    圖1 研究區(qū)影像Fig.1 Image map of the study area

    1.2 方 法

    1.2.1 試驗(yàn)設(shè)計(jì)

    研究通過對(duì)天山云杉的單木樹冠信息提取,對(duì)實(shí)施天然林保護(hù)工程后天山云杉林伐后更新恢復(fù)進(jìn)行評(píng)價(jià),設(shè)計(jì)步驟:對(duì)無人機(jī)原始影像進(jìn)行圖像預(yù)處理;對(duì)預(yù)處理后的影像的前景(樹冠)進(jìn)行增強(qiáng);采用LOG與形態(tài)學(xué)濾波處理原始影像,得到梯度圖像并獲取區(qū)域最大值,進(jìn)而標(biāo)記樹冠頂點(diǎn)與樹冠邊緣;進(jìn)行標(biāo)記控制分水嶺分割提取單木樹冠信息;將分割后的矢量數(shù)據(jù)導(dǎo)入ArcGIS 10.0中進(jìn)行面積統(tǒng)計(jì),以目視解譯勾繪的樹冠面積為參考計(jì)算其精度。

    1.2.2 構(gòu)建梯度幅值圖像

    傳統(tǒng)的分水嶺分割算法與圖像本身沒有很大的相關(guān)性,與圖像的梯度有很強(qiáng)的關(guān)系[21],研究提取樹冠的關(guān)鍵是準(zhǔn)確識(shí)別樹冠的邊緣位置,對(duì)于邊緣檢測(cè)算子現(xiàn)有Sobel算子、Laplacian算子等,對(duì)Laplacian算子和Sobel算子進(jìn)行對(duì)比,經(jīng)過觀察發(fā)現(xiàn)Laplacian算子在邊緣的細(xì)節(jié)識(shí)別與連續(xù)性上要優(yōu)于Sobel算子,但缺點(diǎn)是對(duì)噪聲很敏感。在利用Laplacian算子提取邊緣之前,先對(duì)原始影像進(jìn)行高斯平滑操作。圖2

    圖2 2種算法的邊緣檢測(cè)Fig.2 Edge detection graph of the two algorithms

    1.2.3 形態(tài)學(xué)濾波

    標(biāo)記控制分水嶺分割的關(guān)鍵是選取正確的標(biāo)記。由于研究區(qū)的天山云杉林是天然林,其中天山云杉的年齡、立地條件的不同,會(huì)造成樹冠形狀及大小不一致,以及樹冠內(nèi)部復(fù)雜多變,在圖像中會(huì)產(chǎn)生大量噪聲,這不利于樹冠頂點(diǎn)的探測(cè)。因此,在進(jìn)行樹冠頂點(diǎn)探測(cè)之前,需對(duì)原始圖像進(jìn)行平滑處理。傳統(tǒng)的線性濾波器對(duì)圖像中的低頻分量保留效果好,但不能有效的抑制紋理細(xì)節(jié),造成后續(xù)分水嶺分割時(shí)過分割的現(xiàn)象比較嚴(yán)重。數(shù)學(xué)形態(tài)學(xué)(Mathematical Morphology)是非線性濾波器,使用這種方法對(duì)原始影像進(jìn)行降噪處理時(shí),不會(huì)產(chǎn)生邊緣偏移,因此,使用形態(tài)學(xué)重建的開閉操作對(duì)圖像進(jìn)行降噪處理。

    在形態(tài)學(xué)重建的開閉運(yùn)算的過程中,結(jié)構(gòu)元素的選擇很重要,結(jié)構(gòu)元素選擇過大,會(huì)分離出一些較小的樹冠,造成樹冠提取值偏小,結(jié)構(gòu)元素過小,在較大的樹冠中會(huì)有樹枝被作作樹冠頂點(diǎn),造成過分割現(xiàn)象,選擇合適的結(jié)構(gòu)元素對(duì)樹冠提取的精度有重要的影響。天山云杉在無人機(jī)遙感影像中類似圓形,經(jīng)過多次試驗(yàn),最終確定結(jié)構(gòu)元素為圓盤形,尺寸大小為11個(gè)像元。經(jīng)過上述步驟,原始影像中的噪聲已基本消除,有利于局部最大值的探測(cè),其中a,c,e分別為疏林區(qū)、中郁閉度林區(qū)和密林區(qū)經(jīng)過形態(tài)學(xué)濾波后探測(cè)出的局部最大值圖像,b,d,f分別為疏林區(qū)、中郁閉度林區(qū)和密林區(qū)未經(jīng)形態(tài)學(xué)處理探測(cè)出的局部最大值圖像。圖3

    圖3 局部最大值探測(cè)Fig.3 local maximum detection results

    1.2.4 樹冠標(biāo)記提取

    經(jīng)過基于形態(tài)學(xué)重建的開、閉運(yùn)算后,非樹冠區(qū)域的噪聲和局部極值點(diǎn)明顯減少,但在高郁閉度林分中會(huì)出現(xiàn)樹冠之間相互重疊現(xiàn)象,直接提取會(huì)使一個(gè)樹冠被分成2個(gè)或多個(gè)樹冠,造成過分割。提取局域最大值作為樹冠頂點(diǎn),形成標(biāo)記,作為分水嶺變換的起始點(diǎn)。為防止過分割現(xiàn)象,需要將樹冠與非樹冠區(qū)域區(qū)分開,首先在無人機(jī)原始影像的灰度圖像中給定一個(gè)閾值,將影像中的樹冠與非樹冠區(qū)域分割開。對(duì)比多種閾值分割方法,采用最大類間方差閾值法(OTSU)進(jìn)行最優(yōu)閾值的選擇,該方法使樹冠部分區(qū)域的灰度值、非樹冠區(qū)域的灰度值以及整幅影像的平均灰度值之間的差異最大,能很好的區(qū)別圖像的前景與背景[22]。

    1.2.5 單木樹冠分割

    實(shí)驗(yàn)以拉普拉斯邊緣算子獲取的梯度幅值圖像為基礎(chǔ),采用形態(tài)學(xué)濾波對(duì)樹冠頂點(diǎn)進(jìn)行標(biāo)記,并用最大類間方差閾值法對(duì)樹冠輪廓進(jìn)行提取約束分割。最終使用優(yōu)化后的標(biāo)記控制分水嶺分割方法對(duì)單株樹冠進(jìn)行提取,得到分割后的圖像。

    1.2.6 精度評(píng)價(jià)

    提取的單木樹冠面積與在ArcGIS 10.0軟件中目視解譯勾繪的單木樹冠的面積對(duì)比分析,通過兩者的差異來評(píng)價(jià)自動(dòng)勾繪樹冠的精度,將兩者的關(guān)系分為5類[23],即匹配、接近匹配、丟失、合并和分割[18],將參考樹冠與分割樹冠重合達(dá)50%以上視為匹配,低于50%視為漏分。匹配及接近匹配為正確分割,丟失和合并視為漏分誤差,分割對(duì)應(yīng)錯(cuò)分誤差,用以下公式來表示:

    (4)

    (5)

    (6)

    式中:Ad代表準(zhǔn)確率;Ar代表召回率;F代表F測(cè)度;Nc表示正確分割樹冠個(gè)數(shù);Nr表示參考樹冠個(gè)數(shù);Nd代表自動(dòng)勾繪樹冠個(gè)數(shù)總數(shù)。

    選取的不同郁閉度下9塊樣地的無人機(jī)影像,在Matlab軟件首先進(jìn)行圖像前、背景的互換,再用拉普拉斯算子計(jì)算圖像的灰度梯度圖,并且對(duì)原圖像進(jìn)行形態(tài)學(xué)平滑濾波操作,進(jìn)行降噪處理對(duì)偽樹冠進(jìn)行清除,提取區(qū)域最大值作為樹冠頂點(diǎn),形成標(biāo)記圖像,最后利用標(biāo)記控制分水嶺分割算法進(jìn)行樹冠信息提取。

    在不同郁閉度下均勻隨機(jī)選取30棵樹,通過目視解譯勾繪單木樹冠,之后將自動(dòng)分割的樹冠與目視解譯勾繪的樹冠矢量圖在Arc GIS10.0軟件中進(jìn)行樹冠東西,南北冠幅計(jì)算,計(jì)算出每株樹的樹冠面積。

    2 結(jié)果與分析

    2.1 自動(dòng)勾繪

    研究表明,在疏林區(qū)和中郁閉度林區(qū)由于樹冠邊緣清晰,樹冠重疊現(xiàn)象較少,樹冠信息提取良好,密林區(qū)因?yàn)槠淞址置芏却?,林木年齡跨度大,造成樹冠之間相互重疊,遮擋小樹等現(xiàn)象發(fā)生,在密林區(qū)樹冠分割精度較低。圖4~6

    圖4 疏林區(qū)樹冠提取Fig.4 Canopy extraction results in sparse forest areas

    圖5 中郁閉度林區(qū)樹冠提取Fig.5 Canopy extraction results in medium canopy density forest area

    圖6 密林區(qū)樹冠提取Fig.6 Tree canopy extraction results in dense forest area

    2.2 單株樹冠信息提取

    研究表明,自動(dòng)勾繪的樹冠面積與目視解譯勾繪的樹冠面積相對(duì)誤差最大為0.10,平均相對(duì)誤差為0.04,自動(dòng)勾繪的樹冠面積略低或等于目視解譯后的樹冠面積;自動(dòng)勾繪的樹冠面積與目視解譯勾繪的樹冠面積相對(duì)誤差最大為0.28,平均相對(duì)誤差為0.10,自動(dòng)勾繪的樹冠面積略低或等于目視解譯后的樹冠面積;自動(dòng)勾繪的樹冠面積與目視解譯勾繪的樹冠面積相對(duì)誤差最大為0.43,平均相對(duì)誤差為0.15,在內(nèi)部結(jié)構(gòu)較為復(fù)雜的密林區(qū)存在些許過分割現(xiàn)象;優(yōu)化后的分割算法對(duì)疏林區(qū)、中林區(qū)這種內(nèi)部結(jié)構(gòu)簡(jiǎn)單的林分提取較好,對(duì)于密林區(qū)結(jié)構(gòu)復(fù)雜的林分提取不是很理想。表1~3

    表1 樹冠面積計(jì)算Table 1 Comparison of every tree crown

    表2 樹冠面積計(jì)算Table 2 Comparison of every tree crown

    表3 樹冠面積計(jì)算Table 3 Comparison of every tree crown

    2.3 精度評(píng)價(jià)

    研究表明,疏林區(qū)總共分割出404棵樹冠,目視解譯勾繪出404棵樹冠,其中匹配樹冠個(gè)數(shù)達(dá)到397棵,準(zhǔn)確率為98.26%,召回率為98.26%,F(xiàn)測(cè)度為98.26%,中郁閉度林分共分割出520棵樹冠,目視解譯勾繪出524棵樹冠,其準(zhǔn)確率為93.26%,召回率為92.56%,F(xiàn)測(cè)度為92.91%,密林區(qū)分割出610棵樹冠,目視解譯勾繪出621棵樹冠,由于樹冠相互遮擋、重疊等原因,其準(zhǔn)確率為88.36%,召回率為86.79%,F(xiàn)測(cè)度為87.57%。表4

    表4 精度統(tǒng)計(jì)Table 4 The accuracy statistics of different forest areas(%)

    3 討 論

    單株樹冠的提取,對(duì)評(píng)估樹木生長(zhǎng)狀況以及林分密度等有重要的意義,在森林調(diào)查中必不可少,對(duì)樹冠信息的獲取是建立森林定量化管理平臺(tái)的關(guān)鍵[5],通過分析前人利用標(biāo)記控制分水嶺分割提取單木樹冠的方法,得出其存在易受噪聲影響和過分割的缺點(diǎn),研究對(duì)標(biāo)記控制分水嶺算法進(jìn)行優(yōu)化,使提取樹冠面積的精度有了很大提高,過分割現(xiàn)象在疏、中郁閉度林分中基本消除,為快速動(dòng)態(tài)監(jiān)測(cè)伐后天山云杉林的生長(zhǎng)狀況提供有效技術(shù)手段。

    分水嶺分割方法是以梯度圖像為基礎(chǔ)進(jìn)行分割,在樹冠邊緣位置圖像的灰度值會(huì)發(fā)生驟變,很多學(xué)者[9、24、19]均是利用Sobel算子進(jìn)行梯度圖像提取,Sobel算子對(duì)水平直線檢測(cè)較精準(zhǔn),而天山云杉在遙感影像中結(jié)構(gòu)類似圓狀,對(duì)比各種邊緣檢測(cè)算子后,發(fā)現(xiàn)Laplacian算子在樹冠邊緣檢測(cè)時(shí)表現(xiàn)良好,連續(xù)性好,又因其對(duì)微弱邊緣有較高的反映,并在高郁閉度林分中對(duì)于樹冠重疊區(qū)域的邊緣檢測(cè)精準(zhǔn)度較高。因此,運(yùn)用Laplacian算子對(duì)原始影像進(jìn)行邊緣提取。但其也有一定的缺點(diǎn),Laplacian算子是二階導(dǎo)數(shù)算子對(duì)噪聲非常的敏感,在影像拍攝過程中,由于傳感器的偏移以及電磁干擾等會(huì)使圖像產(chǎn)生較多噪聲,會(huì)對(duì)樹冠邊緣檢測(cè)有所干擾,研究使用高斯濾波先對(duì)影像進(jìn)行濾波操作,有效降低了噪聲的干擾。

    在標(biāo)記前景與背景時(shí)。張海濤[25]使用最小變換(H-minima)技術(shù)對(duì)圖像進(jìn)行分割,該技術(shù)需要人為設(shè)定限制條件,設(shè)定閾值過小,會(huì)產(chǎn)生很多噪聲,設(shè)定閾值過大,對(duì)象邊緣信息會(huì)減弱,由于該文所研究的天山云杉林為天然林,情況較為復(fù)雜,使用最小變換技術(shù)對(duì)圖像分割效果不理想。研究選最大類間方差法來尋找圖像的最優(yōu)閾值,這種方法的優(yōu)勢(shì)是不需要人為的設(shè)定參數(shù),能夠客觀的選擇多區(qū)域的最優(yōu)閾值,并且使各個(gè)類之間分離程度最大,在研究區(qū)分樹冠與非樹冠區(qū)域具有良好的表現(xiàn)。

    使用研究?jī)?yōu)化后的標(biāo)記控制分水嶺分割方法對(duì)樹冠信息進(jìn)行提取,在疏、中郁閉度林分中過分割現(xiàn)象基本消除,在密林區(qū)由于樹冠相互重疊,存在部分過分割現(xiàn)象。但在疏、中、密林區(qū)提取單株樹冠的F測(cè)度均在85%以上,提取精度較高,對(duì)林業(yè)調(diào)查具有一定的參考價(jià)值。

    4 結(jié) 論

    通過改進(jìn)標(biāo)記控制分水嶺分割算法,提取了伐后更新的天山云杉單木樹冠信息,提取單木樹冠信息F測(cè)度在疏林區(qū)為98.26%,中郁閉度林區(qū)為92.91%,密林區(qū)為87.57%,在一定程度上有效解決了標(biāo)記控制分水嶺分割在樹冠提取中的過分割現(xiàn)象,可較好提取新疆山區(qū)天山云杉林的樹冠信息。

    猜你喜歡
    單木分水嶺郁閉度
    地基與無人機(jī)激光雷達(dá)結(jié)合提取單木參數(shù)
    融合LiDAR點(diǎn)云與高分影像的單木檢測(cè)方法研究
    不同郁閉度馬尾松林下種植射干的生長(zhǎng)效果分析
    無人機(jī)影像匹配點(diǎn)云單木識(shí)別算法
    遙感信息(2019年1期)2019-03-22 01:38:16
    郁閉度與七指毛桃生長(zhǎng)的相關(guān)性分析
    2019,一定是個(gè)分水嶺!
    寬葉藍(lán)靛果葉片不同海拔和郁閉度的遺傳變異研究
    森林工程(2018年3期)2018-06-26 03:40:46
    基于雙尺度體元覆蓋密度的TLS點(diǎn)云數(shù)據(jù)單木識(shí)別算法
    森林工程(2018年5期)2018-05-14 13:54:30
    淺談?dòng)糸]度在森林經(jīng)營(yíng)中的科學(xué)應(yīng)用
    “華北第一隧”——張涿高速分水嶺隧道貫通
    舔av片在线| 国产精品综合久久久久久久免费| 欧美最新免费一区二区三区| 日韩成人伦理影院| 小蜜桃在线观看免费完整版高清| 看黄色毛片网站| 国内精品一区二区在线观看| 国产精品国产高清国产av| 国产白丝娇喘喷水9色精品| av天堂中文字幕网| 老熟妇仑乱视频hdxx| 欧美又色又爽又黄视频| 国产精品无大码| 精品乱码久久久久久99久播| 国产精品久久视频播放| 国产又黄又爽又无遮挡在线| 在线a可以看的网站| 久久亚洲精品不卡| 午夜老司机福利剧场| 免费黄网站久久成人精品| 特级一级黄色大片| 亚洲成人精品中文字幕电影| 99国产精品一区二区蜜桃av| 别揉我奶头 嗯啊视频| 久久鲁丝午夜福利片| 嫩草影视91久久| 国产一区二区在线观看日韩| 日本色播在线视频| 乱码一卡2卡4卡精品| 精品国内亚洲2022精品成人| 麻豆一二三区av精品| av专区在线播放| 内射极品少妇av片p| 亚洲经典国产精华液单| 色哟哟哟哟哟哟| 亚洲激情五月婷婷啪啪| 男女下面进入的视频免费午夜| 久久人人爽人人爽人人片va| 国产成人a∨麻豆精品| 亚洲欧美日韩高清专用| 欧美丝袜亚洲另类| 狂野欧美白嫩少妇大欣赏| 又爽又黄无遮挡网站| 九九爱精品视频在线观看| 亚洲三级黄色毛片| 欧美日本视频| 免费无遮挡裸体视频| 变态另类成人亚洲欧美熟女| 我的老师免费观看完整版| 国产真实伦视频高清在线观看| 精品乱码久久久久久99久播| av视频在线观看入口| 热99re8久久精品国产| av在线老鸭窝| 国产探花极品一区二区| 蜜臀久久99精品久久宅男| 久久久久久大精品| aaaaa片日本免费| 日韩在线高清观看一区二区三区| 美女大奶头视频| 精品免费久久久久久久清纯| 亚洲一区高清亚洲精品| 精品少妇黑人巨大在线播放 | 国产精品日韩av在线免费观看| 给我免费播放毛片高清在线观看| 精品久久久噜噜| 欧美三级亚洲精品| 男女边吃奶边做爰视频| 免费看av在线观看网站| 精品福利观看| 丰满的人妻完整版| 成人av在线播放网站| 欧美最黄视频在线播放免费| av中文乱码字幕在线| 中文字幕熟女人妻在线| 久久久久性生活片| 久久久久免费精品人妻一区二区| 一a级毛片在线观看| 卡戴珊不雅视频在线播放| 99热这里只有精品一区| 你懂的网址亚洲精品在线观看 | avwww免费| 卡戴珊不雅视频在线播放| 亚洲图色成人| 国产女主播在线喷水免费视频网站 | 色噜噜av男人的天堂激情| av在线播放精品| 成人亚洲精品av一区二区| 国产毛片a区久久久久| 色噜噜av男人的天堂激情| 亚洲高清免费不卡视频| 级片在线观看| 日本黄色视频三级网站网址| 亚洲精品国产成人久久av| 日本黄色片子视频| 国产老妇女一区| 国产一区亚洲一区在线观看| 日韩 亚洲 欧美在线| 小蜜桃在线观看免费完整版高清| 少妇丰满av| 18禁在线无遮挡免费观看视频 | 九色成人免费人妻av| 最近2019中文字幕mv第一页| 午夜免费男女啪啪视频观看 | 在线观看美女被高潮喷水网站| 国产精品亚洲美女久久久| 人妻制服诱惑在线中文字幕| 国产高潮美女av| 亚洲美女视频黄频| 亚洲无线在线观看| 久久久久国内视频| 99久国产av精品| 赤兔流量卡办理| 欧美三级亚洲精品| 亚洲精品一区av在线观看| 国产私拍福利视频在线观看| 亚洲最大成人av| 精品久久久久久久末码| 非洲黑人性xxxx精品又粗又长| 午夜福利成人在线免费观看| 精品国内亚洲2022精品成人| 插逼视频在线观看| 一级a爱片免费观看的视频| 亚洲内射少妇av| 午夜福利在线观看吧| 午夜激情福利司机影院| 精品久久久噜噜| 九九在线视频观看精品| 亚洲av免费在线观看| 久久久色成人| 免费看光身美女| 久久久久久伊人网av| 亚洲欧美日韩高清专用| 91精品国产九色| .国产精品久久| 国产高清视频在线观看网站| 美女 人体艺术 gogo| 搞女人的毛片| 99在线人妻在线中文字幕| 免费在线观看成人毛片| 亚洲专区国产一区二区| 久久久久久大精品| 精品久久久久久久末码| 99久国产av精品| 亚洲性夜色夜夜综合| 一夜夜www| 日韩欧美精品免费久久| 国产精品爽爽va在线观看网站| 最近手机中文字幕大全| 人妻久久中文字幕网| 中文字幕av在线有码专区| 亚洲国产精品合色在线| 亚洲精品一区av在线观看| 91久久精品国产一区二区三区| 能在线免费观看的黄片| 国产精品电影一区二区三区| 国产高清不卡午夜福利| eeuss影院久久| 免费看美女性在线毛片视频| 国产精品美女特级片免费视频播放器| 亚洲天堂国产精品一区在线| 黄色日韩在线| 国产伦精品一区二区三区视频9| 午夜亚洲福利在线播放| 在线免费观看不下载黄p国产| 中文字幕久久专区| 在线观看免费视频日本深夜| 亚洲精品日韩av片在线观看| 最近视频中文字幕2019在线8| 午夜爱爱视频在线播放| 亚洲欧美日韩高清在线视频| 亚洲人成网站在线播放欧美日韩| 九九久久精品国产亚洲av麻豆| 一区二区三区免费毛片| 国内久久婷婷六月综合欲色啪| 人妻夜夜爽99麻豆av| 我要搜黄色片| 久久久精品大字幕| 人人妻,人人澡人人爽秒播| 日本a在线网址| 精品久久久久久久久久免费视频| 97碰自拍视频| 午夜免费激情av| 内地一区二区视频在线| 亚洲av中文av极速乱| 激情 狠狠 欧美| 久久久久久久久大av| 在线免费十八禁| 熟妇人妻久久中文字幕3abv| 亚州av有码| 欧美色欧美亚洲另类二区| 在线a可以看的网站| 日日撸夜夜添| 97超级碰碰碰精品色视频在线观看| 中文资源天堂在线| 久久久成人免费电影| 最近的中文字幕免费完整| 成人高潮视频无遮挡免费网站| 久久久久久伊人网av| 久久精品国产亚洲网站| 精品熟女少妇av免费看| 国产精品女同一区二区软件| 免费看av在线观看网站| 日本黄色视频三级网站网址| 一个人观看的视频www高清免费观看| 精品久久久噜噜| 日本色播在线视频| 免费观看的影片在线观看| 国产乱人偷精品视频| 欧美+亚洲+日韩+国产| 久久天躁狠狠躁夜夜2o2o| 国产激情偷乱视频一区二区| 亚洲av美国av| 中文字幕免费在线视频6| 国产高潮美女av| 91久久精品国产一区二区成人| 国产伦精品一区二区三区四那| 成人国产麻豆网| 国产片特级美女逼逼视频| 精品少妇黑人巨大在线播放 | 国产老妇女一区| 舔av片在线| 久久人妻av系列| 日韩欧美精品v在线| 禁无遮挡网站| 中国美白少妇内射xxxbb| 日韩欧美三级三区| 国产精品永久免费网站| 2021天堂中文幕一二区在线观| 一级毛片我不卡| 日韩成人伦理影院| 国国产精品蜜臀av免费| 亚洲天堂国产精品一区在线| 男女下面进入的视频免费午夜| 欧美xxxx性猛交bbbb| 成年版毛片免费区| 在线观看美女被高潮喷水网站| 成人美女网站在线观看视频| 久久精品国产清高在天天线| 亚洲av第一区精品v没综合| 成年女人毛片免费观看观看9| 久久精品国产亚洲av天美| 亚洲av一区综合| a级毛色黄片| 美女被艹到高潮喷水动态| 在线观看av片永久免费下载| 国产精品三级大全| 久久人人爽人人片av| 赤兔流量卡办理| 大型黄色视频在线免费观看| 亚洲内射少妇av| 亚洲国产精品成人综合色| 久久精品夜色国产| 中文字幕免费在线视频6| 赤兔流量卡办理| 亚洲激情五月婷婷啪啪| 噜噜噜噜噜久久久久久91| 亚洲高清免费不卡视频| 中文在线观看免费www的网站| 欧美日韩综合久久久久久| 欧美精品国产亚洲| 国产麻豆成人av免费视频| 国产激情偷乱视频一区二区| 无遮挡黄片免费观看| av在线天堂中文字幕| 99热全是精品| 亚洲欧美日韩东京热| 亚洲精品一区av在线观看| 丝袜美腿在线中文| 国产色爽女视频免费观看| 精品人妻视频免费看| 亚洲欧美清纯卡通| 亚洲精品一卡2卡三卡4卡5卡| 久久久久久久午夜电影| 十八禁国产超污无遮挡网站| 在线观看一区二区三区| 久久人人爽人人片av| av在线蜜桃| 六月丁香七月| 亚洲av二区三区四区| 舔av片在线| 国产aⅴ精品一区二区三区波| 精华霜和精华液先用哪个| 亚洲四区av| 国产色婷婷99| 看片在线看免费视频| 又黄又爽又刺激的免费视频.| 99久久久亚洲精品蜜臀av| 最近视频中文字幕2019在线8| 99热6这里只有精品| 亚洲av中文字字幕乱码综合| 国产三级中文精品| а√天堂www在线а√下载| 成人亚洲精品av一区二区| 成人亚洲欧美一区二区av| 十八禁网站免费在线| 欧美丝袜亚洲另类| 免费看光身美女| 高清日韩中文字幕在线| 国产亚洲精品久久久com| 亚洲av中文字字幕乱码综合| 亚洲成人精品中文字幕电影| 国产精品久久视频播放| 啦啦啦观看免费观看视频高清| 女同久久另类99精品国产91| 国产精品一及| 国产成人aa在线观看| 国产伦一二天堂av在线观看| 午夜a级毛片| 亚洲精品影视一区二区三区av| 天堂av国产一区二区熟女人妻| 久久欧美精品欧美久久欧美| 免费观看人在逋| 亚洲一区高清亚洲精品| 波多野结衣高清无吗| 久久久国产成人免费| 亚洲专区国产一区二区| 在线观看午夜福利视频| 欧美人与善性xxx| 有码 亚洲区| 亚洲av成人精品一区久久| 九九久久精品国产亚洲av麻豆| 女的被弄到高潮叫床怎么办| 蜜臀久久99精品久久宅男| 国产麻豆成人av免费视频| 永久网站在线| 亚洲四区av| 日本成人三级电影网站| 欧美成人一区二区免费高清观看| 亚洲经典国产精华液单| 欧美xxxx黑人xx丫x性爽| av在线亚洲专区| 国产一区二区在线av高清观看| 精品久久久噜噜| 91麻豆精品激情在线观看国产| 免费看日本二区| 久久精品人妻少妇| 欧美丝袜亚洲另类| 日本与韩国留学比较| 久久精品夜色国产| 国产乱人视频| 一级黄片播放器| 又粗又爽又猛毛片免费看| 最近在线观看免费完整版| 亚洲国产精品合色在线| 国产成人a∨麻豆精品| 人妻制服诱惑在线中文字幕| 91久久精品国产一区二区成人| 真人做人爱边吃奶动态| 国产一区二区三区在线臀色熟女| 精品久久久久久久人妻蜜臀av| 免费高清视频大片| 国产激情偷乱视频一区二区| 色综合站精品国产| 亚洲人成网站在线观看播放| 久久久久久大精品| 国产午夜精品论理片| 舔av片在线| 日韩成人伦理影院| 亚洲精品国产成人久久av| 少妇的逼水好多| 亚洲精品粉嫩美女一区| 少妇的逼水好多| 久久久久久久午夜电影| 中文资源天堂在线| 在线观看美女被高潮喷水网站| 久久久久免费精品人妻一区二区| 我要搜黄色片| 99久久九九国产精品国产免费| 亚洲av第一区精品v没综合| 国产精品99久久久久久久久| 亚洲七黄色美女视频| 伊人久久精品亚洲午夜| 国产精品久久久久久av不卡| 国产精品一区二区三区四区免费观看 | 精品无人区乱码1区二区| 国产熟女欧美一区二区| 久久久久九九精品影院| 国产精品亚洲一级av第二区| 成人无遮挡网站| 亚洲人成网站在线播放欧美日韩| 久久国内精品自在自线图片| ponron亚洲| 可以在线观看的亚洲视频| 国产爱豆传媒在线观看| 国产成人aa在线观看| 国产三级中文精品| 国产精品一二三区在线看| 亚洲真实伦在线观看| 国产中年淑女户外野战色| 午夜福利在线在线| 色视频www国产| 午夜影院日韩av| 亚洲成a人片在线一区二区| 免费看a级黄色片| 夜夜夜夜夜久久久久| 欧美日本亚洲视频在线播放| 91精品国产九色| www日本黄色视频网| 内地一区二区视频在线| 国产精品一及| 成人欧美大片| 国产一区二区三区在线臀色熟女| 日本熟妇午夜| 俄罗斯特黄特色一大片| 亚洲人成网站高清观看| 久久国产乱子免费精品| 91久久精品国产一区二区成人| 我要看日韩黄色一级片| 免费黄网站久久成人精品| 成年版毛片免费区| 三级经典国产精品| 精品午夜福利在线看| 精品一区二区三区视频在线| 久久久久国产精品人妻aⅴ院| 亚洲内射少妇av| 别揉我奶头 嗯啊视频| 少妇丰满av| 日本撒尿小便嘘嘘汇集6| 九九在线视频观看精品| 久久6这里有精品| 日日摸夜夜添夜夜添小说| 热99在线观看视频| 欧美色欧美亚洲另类二区| 精品熟女少妇av免费看| 夜夜爽天天搞| 国产午夜精品论理片| 国产男人的电影天堂91| 成年免费大片在线观看| 麻豆久久精品国产亚洲av| 性色avwww在线观看| 一个人看视频在线观看www免费| 麻豆国产97在线/欧美| 日本一本二区三区精品| 午夜影院日韩av| 秋霞在线观看毛片| 亚洲中文字幕一区二区三区有码在线看| 精品免费久久久久久久清纯| 97超碰精品成人国产| 五月玫瑰六月丁香| 国产91av在线免费观看| 国产真实伦视频高清在线观看| 精品人妻熟女av久视频| 日本-黄色视频高清免费观看| 欧美在线一区亚洲| 尾随美女入室| 变态另类丝袜制服| 成人美女网站在线观看视频| 国内揄拍国产精品人妻在线| 国产淫片久久久久久久久| 国产日本99.免费观看| 少妇熟女aⅴ在线视频| 亚洲国产欧美人成| 老女人水多毛片| 最新中文字幕久久久久| 亚洲欧美中文字幕日韩二区| 一区二区三区高清视频在线| 在现免费观看毛片| 成人亚洲欧美一区二区av| 美女内射精品一级片tv| 黄色一级大片看看| 免费搜索国产男女视频| 麻豆精品久久久久久蜜桃| 国产大屁股一区二区在线视频| 亚洲欧美成人综合另类久久久 | 女同久久另类99精品国产91| 免费观看人在逋| 日韩一区二区视频免费看| 97超视频在线观看视频| 亚洲精品456在线播放app| 国产男靠女视频免费网站| 国产成人aa在线观看| 国产极品精品免费视频能看的| 亚洲五月天丁香| 国产成人freesex在线 | 在线播放无遮挡| 亚洲熟妇熟女久久| 欧美一级a爱片免费观看看| av免费在线看不卡| 亚洲不卡免费看| 日本 av在线| 午夜精品国产一区二区电影 | 国产在视频线在精品| 日韩精品有码人妻一区| 12—13女人毛片做爰片一| 最近手机中文字幕大全| 欧美区成人在线视频| 免费搜索国产男女视频| 免费看光身美女| 激情 狠狠 欧美| 久久久欧美国产精品| 18禁在线播放成人免费| 国产视频一区二区在线看| 最近在线观看免费完整版| 亚洲自拍偷在线| 成人亚洲精品av一区二区| 免费不卡的大黄色大毛片视频在线观看 | 亚洲国产精品sss在线观看| 久久精品国产亚洲av天美| 国产一区二区在线av高清观看| 欧美激情国产日韩精品一区| 禁无遮挡网站| 免费看a级黄色片| 国产乱人偷精品视频| 国产成人福利小说| 国产精品一区二区三区四区免费观看 | 少妇熟女aⅴ在线视频| 高清毛片免费观看视频网站| 波多野结衣高清无吗| 欧美成人免费av一区二区三区| 黑人高潮一二区| 国产精品电影一区二区三区| 欧美激情在线99| 午夜福利在线观看免费完整高清在 | 亚洲欧美精品综合久久99| 国产蜜桃级精品一区二区三区| 天堂av国产一区二区熟女人妻| 亚洲av第一区精品v没综合| 全区人妻精品视频| 久久久国产成人免费| 1000部很黄的大片| 午夜免费激情av| 搡老熟女国产l中国老女人| 亚洲av不卡在线观看| 我的老师免费观看完整版| 国产精华一区二区三区| 精品一区二区三区人妻视频| 一边摸一边抽搐一进一小说| 亚洲人成网站在线观看播放| 成年女人永久免费观看视频| 如何舔出高潮| 蜜臀久久99精品久久宅男| 亚洲一区高清亚洲精品| 午夜福利在线在线| 午夜精品在线福利| 成年版毛片免费区| 亚洲aⅴ乱码一区二区在线播放| 熟妇人妻久久中文字幕3abv| 午夜亚洲福利在线播放| 日韩成人伦理影院| 亚洲国产欧美人成| 中文在线观看免费www的网站| 人人妻人人澡欧美一区二区| 亚洲国产色片| 日本黄大片高清| 18禁在线无遮挡免费观看视频 | 久久这里只有精品中国| 亚洲不卡免费看| 日本与韩国留学比较| 性欧美人与动物交配| 欧美一区二区精品小视频在线| 激情 狠狠 欧美| 特级一级黄色大片| 97超视频在线观看视频| 欧美另类亚洲清纯唯美| 午夜福利高清视频| 国产人妻一区二区三区在| 国产高清视频在线播放一区| 超碰av人人做人人爽久久| 伊人久久精品亚洲午夜| 亚洲av成人av| 欧美成人精品欧美一级黄| 国产午夜精品久久久久久一区二区三区 | 国产男人的电影天堂91| 午夜精品一区二区三区免费看| 无遮挡黄片免费观看| 中文字幕免费在线视频6| 男女下面进入的视频免费午夜| 色尼玛亚洲综合影院| 1000部很黄的大片| 一区二区三区四区激情视频 | 成年av动漫网址| 色噜噜av男人的天堂激情| 亚洲欧美中文字幕日韩二区| 精品一区二区三区视频在线观看免费| 春色校园在线视频观看| 内射极品少妇av片p| 国产精品精品国产色婷婷| 欧美绝顶高潮抽搐喷水| 精品免费久久久久久久清纯| 搞女人的毛片| 黑人高潮一二区| 亚洲在线观看片| 亚洲精品国产成人久久av| 欧美人与善性xxx| 六月丁香七月| 国内精品美女久久久久久| 美女xxoo啪啪120秒动态图| 一个人看的www免费观看视频| 六月丁香七月| 日本黄色片子视频| 久久久欧美国产精品| 亚洲美女视频黄频| 亚洲成a人片在线一区二区| 男人舔奶头视频| 欧美成人精品欧美一级黄| 亚洲av.av天堂| 黄色欧美视频在线观看| 丰满人妻一区二区三区视频av| 97超碰精品成人国产| 色哟哟·www| 成人精品一区二区免费| 一个人观看的视频www高清免费观看| 国产午夜福利久久久久久| 美女高潮的动态| 午夜福利高清视频| 欧美3d第一页| 香蕉av资源在线| 精品久久久久久久久久免费视频| 色综合站精品国产| 欧美高清成人免费视频www| 国产国拍精品亚洲av在线观看| 久久精品国产亚洲av涩爱 | 亚洲av一区综合| 精品福利观看|