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

    基于航空點云的落葉松林分蓄積量建模和應(yīng)用測試初探

    2021-10-08 08:19:02高金萍孫忠秋于慧娜高顯連吳發(fā)云劉迎春
    林業(yè)資源管理 2021年4期
    關(guān)鍵詞:蓄積量蓄積落葉松

    高金萍,孫忠秋,于慧娜,高顯連,吳發(fā)云,劉迎春

    (國家林業(yè)和草原局調(diào)查規(guī)劃設(shè)計院,北京 100714)

    1 研究背景

    森林蓄積量是森林資源監(jiān)測的核心調(diào)查內(nèi)容之一,它反映了一個國家或地區(qū)森林資源總規(guī)?;蛩?也是衡量森林生態(tài)狀況的重要指標(biāo)之一。我國從“十二五”開始,將其作為約束性指標(biāo)納入國民經(jīng)濟與社會發(fā)展規(guī)劃。從20世紀(jì)70年代開始,我國建立的全國森林資源清查體系和各省(地區(qū))森林資源規(guī)劃設(shè)計調(diào)查制度,均是以野外地面調(diào)查為主體,通過專業(yè)的地面調(diào)查隊伍,獲取權(quán)威、可靠的胸徑、樹高等重要指標(biāo),再結(jié)合本地樹種蓄積模型,推算得到森林蓄積量指標(biāo)。由于我國森林分布范圍廣、面積大、大多地處偏僻遙遠的林區(qū)或山區(qū),傳統(tǒng)的地面調(diào)查需要獲取每株樣木胸徑和樹高,調(diào)查任務(wù)繁重、條件艱苦、工作強度極大,需要耗費大量的人力、物力和財力[1]。當(dāng)前,森林資源監(jiān)測工作面臨調(diào)查人員匱乏、調(diào)查隊伍結(jié)構(gòu)老化、調(diào)查成本日趨高漲的困境,迫切需要創(chuàng)新突破和升級換代監(jiān)測技術(shù)和裝備設(shè)施,以適應(yīng)新時期森林資源監(jiān)測工作的新形勢和新要求[2]。

    對地觀測與定位技術(shù)的興起,為森林資源動態(tài)監(jiān)測提供了先進感知手段。近些年來,光學(xué)遙感技術(shù)的發(fā)展解決了傳統(tǒng)地面調(diào)查存在的時間和人力耗費過高的不足,基于光學(xué)遙感影像中提取的植被指數(shù)、紋理特征等信息,通過構(gòu)建地面樣地實測蓄積與遙感變量之間的回歸模型,可實現(xiàn)對森林蓄積量的定量估計[3-4]。但由于光學(xué)遙感難以獲取森林三維結(jié)構(gòu)信息,且在森林郁閉度高的區(qū)域會存在光譜信號飽和等問題[5],在實際工作中還是難以滿足森林資源調(diào)查精度的要求。

    激光雷達技術(shù)可以直接、精確地獲取植被的三維結(jié)構(gòu)信息,有效彌補光學(xué)遙感林業(yè)應(yīng)用能力的不足[6]。國外自20世紀(jì)80年代開始開展森林參數(shù)估測研究[7],Keiko等[8]利用激光雷達點云估測日本西部闊葉林的樣地蓄積量,Vincenzo[9]估測了意大利米蘭都市區(qū)周邊闊葉森林的蓄積量及地上生物量。國內(nèi)近10年來,利用激光雷達點云高度、密度等信息構(gòu)建森林結(jié)構(gòu)、估測森林蓄積的研究逐步擴大化。2014年,曹林[10]借助機載激光雷達點云,估測了江蘇常熟虞山林場北亞熱帶天然次生林森林參數(shù);2018年,劉浩等[11]借助機載激光雷達點云,預(yù)測了中國東南沿海平原人工林林分特征;2020年,曾偉生等[12]基于東北林區(qū)紅松林樣地的機載激光雷達數(shù)據(jù)和地面實測數(shù)據(jù),建立了紅松林分蓄積量模型,預(yù)估精度達到森林資源調(diào)查相關(guān)技術(shù)規(guī)定要求,并開展了航空蓄積表理論研究。隨著近年來激光雷達數(shù)據(jù)獲取成本日趨降低,激光雷達技術(shù)由科學(xué)研究轉(zhuǎn)為業(yè)務(wù)探索時機已趨成熟,但盡管建模理論研究成果較多,大多都局限于模型和精度理論探討,結(jié)合實際生產(chǎn)實踐業(yè)務(wù)數(shù)據(jù)對模型進行應(yīng)用測試和分析評價的較少。

    2017年,原國家林業(yè)局調(diào)查規(guī)劃設(shè)計院(以下簡稱國家林草局規(guī)劃院)啟動了“星機地”綜合試驗計劃,擬在我國東北、華北、西北、南方、西南、西藏6個區(qū)域(12個省左右)開展有人機、無人機和地面調(diào)查相結(jié)合的“星機地”綜合試驗,以獲取全國典型森林類型(包括占全國85%蓄積量的主要樹種)的激光雷達及地面樣地數(shù)據(jù)等信息(約4 500個樣地)。該計劃的一個重要目標(biāo)就是開展激光點云林分蓄積量建模和業(yè)務(wù)化應(yīng)用研究,最終形成覆蓋全國主要森林分區(qū)和主要樹種的、面向林業(yè)生產(chǎn)實踐的航空蓄積表模型體系。2017—2020年,先后在湖南張家界周邊、東北虎豹國家公園(以下簡稱東北虎豹公園)、福建武夷山國家公園和海南熱帶雨林國家公園開展了空地一體化調(diào)查。本文將以東北虎豹公園中落葉松為例,探索基于航空點云的林分蓄積量建模和業(yè)務(wù)化應(yīng)用測試方案,為其他樹種林分蓄積量建模和小班蓄積量估測提供技術(shù)路線和方法,為逐步構(gòu)建我國面向生產(chǎn)的分林區(qū)、分樹種的航空蓄積表模型體系提供數(shù)據(jù)和模型支撐。

    2 研究區(qū)與數(shù)據(jù)概況

    2.1 研究區(qū)簡介

    東北虎豹公園地跨黑龍江和吉林兩省,其中吉林省片區(qū)占70%面積以上,東起吉林省琿春林業(yè)局青龍臺林場,西至吉林省汪清縣林業(yè)局南溝林場,南自吉林省琿春林業(yè)局敬信林場,北到黑龍江省東京城林業(yè)局奮斗林場,總面積近1.5萬km2。公園處于亞洲溫帶針闊葉混交林生態(tài)系統(tǒng)的中心地帶,區(qū)域內(nèi)有眾多的溫帶森林植被類型,是北半球溫帶區(qū)生物多樣性最豐富的地區(qū)之一。2018年9—10月,國家林草局規(guī)劃院在東北虎豹公園開展了有人機、無人機飛行作業(yè)及地面樣地調(diào)查工作,獲取了研究區(qū)大量有人機和無人機激光點云數(shù)據(jù)及對應(yīng)的地面樣地、樣木調(diào)查數(shù)據(jù),為基于航空點云林分蓄積模型和小班林分蓄積量估測研究提供了珍貴的空地一體化調(diào)查基礎(chǔ)數(shù)據(jù)。

    2.2 研究區(qū)數(shù)據(jù)分析

    根據(jù)東北內(nèi)蒙古重點國有林區(qū)二類調(diào)查主要樹種統(tǒng)計結(jié)果,東北虎豹公園空地調(diào)查主要樹種為落葉松(Larixspp.)、云杉(Piceaspp.)、冷杉(Abiesspp.)、樺木樹(Betulaspp.)、楊樹(Populusspp.)、椴樹樹(Tiliaspp.)、櫟類(Quercusspp.)7種純林樹種和針葉混、針闊混、闊葉混3類混交林分類型。其中落葉松是研究區(qū)內(nèi)分布最為廣泛的樹種之一,且根據(jù)第9次全國森林資源清查結(jié)果,在針葉林中,落葉松分布面積居第一位,蓄積量居第四位,是我國重要的針葉樹種[12]。

    2.2.1有無人機激光雷達調(diào)查數(shù)據(jù)

    激光雷達調(diào)查數(shù)據(jù)以有人機激光雷達點云為主,對有人機未覆蓋的區(qū)域采用無人機補飛。有人機飛行區(qū)域主要分布在大興溝、天橋嶺、穆棱、東京城等4個林業(yè)局及部分林場(總面積約5 600km2);無人機飛行區(qū)域包含17個1 km2的航空調(diào)查大樣地,主要分布在大興溝、天橋嶺、汪清、琿春等4個林業(yè)局。激光雷達傳感器主要參數(shù)及調(diào)查數(shù)據(jù)分布范圍分別如表1和圖1所示。

    表1 激光雷達傳感器參數(shù)

    圖1 激光雷達點云數(shù)據(jù)覆蓋范圍

    2.2.2地面調(diào)查樣地數(shù)據(jù)

    在東北虎豹公園的航空(包括有人機和無人機)飛行覆蓋區(qū)域內(nèi),對需要建模的7個純林樹種和3個混交林布設(shè)地面調(diào)查樣地??紤]到不同樹高、郁閉度等級對林分蓄積量的影響,將每個樹種劃分5個樹高級(高、較高、中等、較低、低)及3個郁閉度等級(0.2~0.5,0.5~0.8,0.8以上),樣地布設(shè)時盡可能保證建模樣本涵蓋所有樹高級和郁閉度等級,以提高模型的適用性。本研究共布設(shè)并調(diào)查完成488個圓形樣地(半徑15m),其中落葉松樣地57個,其主要調(diào)查因子如表2所示。表3展示了部分落葉松樣地的關(guān)鍵調(diào)查因子數(shù)據(jù)。

    表2 地面調(diào)查樣地數(shù)據(jù)表結(jié)構(gòu)(部分因子)

    表3 樣地主要調(diào)查因子示例(部分樣地)

    2.2.3落葉松小班二類本底調(diào)查成果

    2016—2018年,東北內(nèi)蒙古重點國有林區(qū)完成了二類調(diào)查,東北虎豹公園建立了覆蓋全境的、時間節(jié)點統(tǒng)一至2018年底的二類調(diào)查本底。以大興溝林業(yè)局為例,基于二類調(diào)查本底成果,提取大興溝林業(yè)局內(nèi)優(yōu)勢樹種為落葉松的小班,全部落葉松小班分布圖如圖2所示。

    圖2 大興溝林業(yè)局落葉松小班空間分布圖

    3 研究方法和結(jié)果

    根據(jù)落葉松地面樣地邊界,提取落葉松樣地激光點云,建立地面樣地實測蓄積與激光點云中提取的高度、密度、覆蓋度等變量之間的回歸關(guān)系,構(gòu)建基于航空點云調(diào)查數(shù)據(jù)的落葉松林分蓄積量回歸模型并分析其理論精度。

    為推進理論模型業(yè)務(wù)化研究,嘗試?yán)媒⒌哪P统晒麑ρ芯繀^(qū)落葉松小班單位面積蓄積量進行估算;通過與二類調(diào)查獲取的蓄積量數(shù)據(jù)對比分析,評估基于航空點云的林分蓄積量模型相對精度,初步評估基于航空點云的林分蓄積量模型進入業(yè)務(wù)應(yīng)用的潛力和能力。本研究先在較小的范圍(21個落葉松小班)開展模型應(yīng)用測試,取得較為理想效果;然后將測試范圍擴大到2018年開展過二類實地調(diào)查的大興溝林業(yè)局(約260多個落葉松小班)??傮w技術(shù)路線如圖3所示。

    圖3 總體技術(shù)路線

    3.1 建模變量提取

    基于精確分類后的激光雷達調(diào)查數(shù)據(jù)和落葉松地面樣地范圍,裁切出與地面樣地位置和大小一致的落葉松激光點云數(shù)據(jù)。利用點云中分類得到的地面點通過TIN插值生成地面高程模型(DEM),基于DEM的地面高程信息對點云進行歸一化處理,得到歸一化后的落葉松點云數(shù)據(jù)。從落葉松點云數(shù)據(jù)中提取建模所需變量,包括系統(tǒng)提供的46個高度變量、10個密度變量、1個覆蓋度變量,共計57個統(tǒng)計變量,各類變量內(nèi)涵及分析如表4所示。

    表4 激光點云數(shù)據(jù)初始變量分析

    3.2 蓄積估測模型構(gòu)建與評價

    建立蓄積量與點云變量之間的關(guān)系,一般采用以下非線性回歸模型:

    V=β0x1β1x2β2…xnβn+ε

    (1)

    式中:V為蓄積量;β0,β1,β2,…,βn為待定回歸系數(shù);x1,x2,…,xn為建模變量;ε為誤差項。

    由于上述非線性模型求解較為復(fù)雜,通常會采用對數(shù)回歸的方法,也就是將所有變量進行對數(shù)變換,將多元非線性問題轉(zhuǎn)化為多元線性回歸來解決:

    lnV=α0+α1lnx1+α2lnx2+…+αnlnxn+δ

    (2)

    式中:V為蓄積量;α0,α1,α2,…,αn為待定回歸系數(shù);x1,x2,…,xn為建模變量;δ為誤差項。

    利用40個落葉松樣地進行建模,余下17個樣地用來驗證模型估測精度。采用多元線性回歸方法時,運用逐步回歸法合適變量。由于提取的變量較多,一些與蓄積相關(guān)性不高的變量可以在逐步回歸之前先舍去。因此,將提取的變量與需要預(yù)測的蓄積量進行Pearson相關(guān)分析,根據(jù)Pearson相關(guān)系數(shù)大小,初步剔除Pearson相關(guān)系數(shù)較低(<0.3)的變量。

    利用逐步回歸得到的回歸系數(shù)構(gòu)建林分蓄積量模型,采用決定系數(shù)(R2)、均方根誤差(RMSE)、相對均方根誤差(rRMSE)、平均絕對誤差(MAE)、平均絕對百分比誤差(MAPE)、總體相對誤差(TRE)來衡量模型估測精度。

    利用逐步回歸篩選進入模型的變量,分別是50%高度百分位數(shù)(h50th)和覆蓋度(CC),建模結(jié)果為:

    lnV=0.918lnh50th+1.303lnCC+3.091

    (3)

    利用建立的模型估測17個落葉松驗證樣地的單位面積蓄積,繪制基于模型的蓄積估測值和地面樣地蓄積實測值的散點圖(圖4),得到的R2為0.84,RMSE為23.77m3/hm2,rRMSE為0.18,MAE為20.08 m3/hm2,MAPE為16.21%,TRE為4.5%。從驗證結(jié)果來看,基于激光點云的落葉松林分蓄積量模型具有較高的精度。按照《森林資源規(guī)劃設(shè)計調(diào)查技術(shù)規(guī)程》[14]中小班調(diào)查蓄積量的精度等級標(biāo)準(zhǔn),A,B,C 三級要求相對誤差分別不超過15%,20%和25%,基于激光點云的落葉松林分蓄積量模型接近A級、可以滿足B 級的精度要求。

    圖4 驗證樣地蓄積估測與實測值散點圖

    3.3 應(yīng)用模型估算樣地所在落葉松小班單位面積蓄積

    基于東北虎豹公園二類小班調(diào)查本底成果,提取落葉松地面樣地所在的二類小班21個(小班平均面積約15hm2);疊加空地一體化調(diào)查獲取的激光雷達點云數(shù)據(jù)后,獲取這21個二類小班的航空點云數(shù)據(jù);然后利用已建立的落葉松林分蓄積量模型估測小班單位面積蓄積量。經(jīng)分析,小班蓄積估測值與小班實際調(diào)查值的散點圖如圖5所示,落葉松小班二類蓄積量與基于激光點云的模型估算得到的蓄積量比較情況如表5所示。

    圖5 21個落葉松小班蓄積估測與二類本底調(diào)查蓄積散點圖

    表5 落葉松小班估測蓄積與二類調(diào)查結(jié)果對比分析表

    (續(xù)表)

    由此可以看出,二者的R2為0.78,呈現(xiàn)較好的相關(guān)性;小班蓄積總量估測相對精度為89.4%。

    3.3 應(yīng)用模型估算大興溝林業(yè)局二類調(diào)查落葉松小班單位面積蓄積

    基于東北虎豹公園二類小班調(diào)查本底成果,提取2018年開展實地調(diào)查的全部落葉松小班266個(絕大部分落在大興溝林業(yè)局,小班平均面積約7.22hm2);疊加空地一體化調(diào)查獲取的激光雷達點云數(shù)據(jù)后,獲取這266個二類小班的航空點云數(shù)據(jù);然后利用已建立的落葉松林分蓄積量模型估測小班單位面積蓄積量。經(jīng)分析,小班蓄積估測值與小班實際調(diào)查值的散點圖如圖6所示。

    圖6 266個落葉松小班蓄積估測與二類本底調(diào)查蓄積散點圖

    266個小班蓄積模型估測總量 298 627.58hm2,二類調(diào)查得到小班蓄積總量 301 914.00hm2,相差約 3 286.42 hm2,小班蓄積總量估測相對精度為98.9%。東北虎豹公園二類調(diào)查采用角規(guī)樣地抽樣推算小班蓄積、布設(shè)控制樣地控制總體精度的調(diào)查方式,從總體相對精度來看,可以認(rèn)為模型應(yīng)用有效性較好。

    4 結(jié)論與討論

    本文以東北林區(qū)廣泛分布的落葉松為重點案例,探索了激光雷達調(diào)查方式下的落葉松林分蓄積量模型研建方法,并結(jié)合東北內(nèi)蒙古重點國有林區(qū)二類調(diào)查成果,運用建模成果估算了該地區(qū)落葉松二類小班單位面積蓄積,是激光雷達科研模型成果在業(yè)務(wù)化應(yīng)用方向上的一次有效嘗試。從研究結(jié)果和研究區(qū)該樹種的應(yīng)用情況來看,激光點云的落葉松林分蓄積量模型預(yù)估精度達到《森林資源規(guī)劃設(shè)計調(diào)查技術(shù)規(guī)程》的要求,可以在實踐中推廣應(yīng)用。

    鑒于我國地理區(qū)域差異大、森林類型多、森林狀況復(fù)雜,要想建立起覆蓋全國主要樹種的、面向生產(chǎn)應(yīng)用的航空蓄積表模型體系,需要十年如一日花大氣力研究每個樹種模型的普適性、精度以及參數(shù)的簡單易得和穩(wěn)定性。就本研究來說,由于經(jīng)費和飛行區(qū)域受限,本次調(diào)查主要覆蓋東北虎豹公園西部,東部地區(qū)較少,樣本代表性有一定影響,且每個建模樹種樣本數(shù)均在50個上下。后續(xù)將采用“抽十法”或“留一法”等交叉驗證方法,以提高模型穩(wěn)定度。此外,本次模型應(yīng)用測試范圍為266個小班,數(shù)據(jù)量中等。后續(xù)將對整個東北虎豹公園研究區(qū)的落葉松小班(約2 000多個)進行模型應(yīng)用和測試,分析大規(guī)模業(yè)務(wù)應(yīng)用中模型精度、運算效率等方面可能出現(xiàn)的問題,在模型迭代應(yīng)用中實現(xiàn)模型優(yōu)化和效率提升。

    猜你喜歡
    蓄積量蓄積落葉松
    山西落葉松雜交良種逾10萬畝
    長白落葉松離體再生體系的建立
    一元立木材積表計算蓄積量誤差的探討
    樟子松人工林林分蓄積量計算方法的對比分析
    神奇的落葉松提取物
    藏藥佐太中汞在小鼠體內(nèi)的蓄積
    中成藥(2017年7期)2017-11-22 07:32:46
    林分蓄積量估算方法研究
    淺談藏醫(yī)“疾病蓄積”與健康養(yǎng)生
    2015年湖南省活立木蓄積量、森林覆蓋率排名前10位的縣市區(qū)
    堅持創(chuàng)新 蓄積行業(yè)永續(xù)發(fā)展動能
    国产精品二区激情视频| 大话2 男鬼变身卡| 中文字幕亚洲精品专区| 国产99久久九九免费精品| 亚洲精品国产一区二区精华液| 成年美女黄网站色视频大全免费| 亚洲成av片中文字幕在线观看| 宅男免费午夜| 热99久久久久精品小说推荐| 国产一卡二卡三卡精品| 熟女av电影| 久久人人97超碰香蕉20202| 真人做人爱边吃奶动态| 国产亚洲午夜精品一区二区久久| 亚洲欧洲日产国产| 精品一区二区三区av网在线观看 | 国产男女内射视频| 超色免费av| 国产精品 国内视频| 亚洲,欧美,日韩| 最近手机中文字幕大全| av网站在线播放免费| 久久99精品国语久久久| 亚洲美女黄色视频免费看| 久久久久久久久久久久大奶| 高清视频免费观看一区二区| 亚洲av在线观看美女高潮| 老司机深夜福利视频在线观看 | 国产精品一区二区在线不卡| 国产视频首页在线观看| 亚洲国产中文字幕在线视频| 亚洲黑人精品在线| 国产男女超爽视频在线观看| 岛国毛片在线播放| 国产精品二区激情视频| 十八禁网站网址无遮挡| 一二三四社区在线视频社区8| 99精国产麻豆久久婷婷| 日本a在线网址| 视频区欧美日本亚洲| 99国产精品免费福利视频| 91精品三级在线观看| 深夜精品福利| 在线观看人妻少妇| 精品一区二区三卡| 80岁老熟妇乱子伦牲交| 欧美黑人欧美精品刺激| 久久精品人人爽人人爽视色| 色网站视频免费| 国产一卡二卡三卡精品| 亚洲欧美一区二区三区国产| 一边摸一边抽搐一进一出视频| 午夜激情av网站| 亚洲精品一卡2卡三卡4卡5卡 | bbb黄色大片| 大型av网站在线播放| 黄片小视频在线播放| 国产xxxxx性猛交| 久久久久久久国产电影| 亚洲成人免费电影在线观看 | 日韩,欧美,国产一区二区三区| 国产淫语在线视频| 精品少妇一区二区三区视频日本电影| 亚洲精品中文字幕在线视频| 亚洲av综合色区一区| 欧美性长视频在线观看| 久久精品久久久久久噜噜老黄| 久久人妻福利社区极品人妻图片 | 首页视频小说图片口味搜索 | 国产黄色视频一区二区在线观看| 亚洲成人手机| 国产视频首页在线观看| 首页视频小说图片口味搜索 | 麻豆乱淫一区二区| 免费在线观看视频国产中文字幕亚洲 | 蜜桃国产av成人99| 免费看av在线观看网站| 成人国产一区最新在线观看 | 久久性视频一级片| 菩萨蛮人人尽说江南好唐韦庄| 黑丝袜美女国产一区| h视频一区二区三区| 精品视频人人做人人爽| 黄片小视频在线播放| 大香蕉久久网| 大型av网站在线播放| 国产极品粉嫩免费观看在线| 日本wwww免费看| 男女边摸边吃奶| 国产精品成人在线| 亚洲国产av影院在线观看| 亚洲国产欧美在线一区| 亚洲欧美日韩高清在线视频 | 久久久久网色| 婷婷色综合www| 精品福利观看| 国产在线一区二区三区精| av欧美777| 精品国产乱码久久久久久男人| cao死你这个sao货| 久久精品亚洲av国产电影网| 精品免费久久久久久久清纯 | 国产在线观看jvid| 亚洲精品一区蜜桃| 精品人妻在线不人妻| 精品人妻在线不人妻| 国产在线视频一区二区| 成人亚洲精品一区在线观看| 大话2 男鬼变身卡| 日韩制服骚丝袜av| 国产精品欧美亚洲77777| 亚洲国产中文字幕在线视频| 久久综合国产亚洲精品| 国产精品熟女久久久久浪| 亚洲精品成人av观看孕妇| 午夜福利视频在线观看免费| 亚洲国产精品999| 肉色欧美久久久久久久蜜桃| 麻豆国产av国片精品| 色综合欧美亚洲国产小说| 天堂8中文在线网| 亚洲国产欧美网| 色播在线永久视频| 精品少妇久久久久久888优播| 久久热在线av| 免费日韩欧美在线观看| 每晚都被弄得嗷嗷叫到高潮| 国产野战对白在线观看| 黄色 视频免费看| 男人爽女人下面视频在线观看| 亚洲三区欧美一区| 亚洲欧美色中文字幕在线| 国产伦理片在线播放av一区| 天天添夜夜摸| 如日韩欧美国产精品一区二区三区| 性少妇av在线| 中文精品一卡2卡3卡4更新| 久久这里只有精品19| 黄色片一级片一级黄色片| 精品欧美一区二区三区在线| 校园人妻丝袜中文字幕| 中文字幕精品免费在线观看视频| 国产亚洲一区二区精品| 国产亚洲精品第一综合不卡| 女人高潮潮喷娇喘18禁视频| 男人舔女人的私密视频| 天天躁夜夜躁狠狠躁躁| 欧美+亚洲+日韩+国产| 日本五十路高清| 无遮挡黄片免费观看| 免费少妇av软件| avwww免费| 精品视频人人做人人爽| 国产免费又黄又爽又色| 91精品三级在线观看| 国产成人精品无人区| 国产av精品麻豆| 女性生殖器流出的白浆| 成人免费观看视频高清| 视频区欧美日本亚洲| 精品久久久久久久毛片微露脸 | 视频区欧美日本亚洲| 久久女婷五月综合色啪小说| 成在线人永久免费视频| 香蕉丝袜av| 只有这里有精品99| 最新在线观看一区二区三区 | 不卡av一区二区三区| 丰满少妇做爰视频| 精品卡一卡二卡四卡免费| 天堂8中文在线网| 亚洲自偷自拍图片 自拍| 蜜桃国产av成人99| h视频一区二区三区| 亚洲国产精品一区二区三区在线| 亚洲五月婷婷丁香| 亚洲熟女精品中文字幕| 国产高清videossex| 免费在线观看视频国产中文字幕亚洲 | 18禁裸乳无遮挡动漫免费视频| 18禁黄网站禁片午夜丰满| 日韩制服骚丝袜av| 无限看片的www在线观看| 啦啦啦啦在线视频资源| 国产不卡av网站在线观看| 国产免费现黄频在线看| 免费观看人在逋| 伊人亚洲综合成人网| 一区二区三区精品91| 在线观看免费午夜福利视频| 亚洲欧洲精品一区二区精品久久久| 超碰成人久久| 国产免费又黄又爽又色| 国产免费福利视频在线观看| 免费不卡黄色视频| 精品少妇久久久久久888优播| 一级黄色大片毛片| 日本五十路高清| 一个人免费看片子| 国产真人三级小视频在线观看| 熟女av电影| 久久人妻熟女aⅴ| 婷婷色麻豆天堂久久| 亚洲欧洲国产日韩| 国产成人a∨麻豆精品| 欧美日韩黄片免| 少妇被粗大的猛进出69影院| 日韩av不卡免费在线播放| 久久这里只有精品19| 亚洲国产欧美在线一区| 亚洲成人手机| 亚洲专区中文字幕在线| 捣出白浆h1v1| 中国美女看黄片| 性少妇av在线| 国产成人av激情在线播放| 9191精品国产免费久久| 超碰97精品在线观看| 999久久久国产精品视频| 两个人免费观看高清视频| 精品国产一区二区三区四区第35| 精品久久久久久久毛片微露脸 | 日韩一卡2卡3卡4卡2021年| 成人18禁高潮啪啪吃奶动态图| 男人爽女人下面视频在线观看| 亚洲av成人不卡在线观看播放网 | 晚上一个人看的免费电影| 一级黄色大片毛片| 亚洲国产看品久久| 国产欧美亚洲国产| 日本wwww免费看| 99久久99久久久精品蜜桃| 中文字幕亚洲精品专区| 在线天堂中文资源库| 这个男人来自地球电影免费观看| 别揉我奶头~嗯~啊~动态视频 | 久久天躁狠狠躁夜夜2o2o | 国产日韩欧美视频二区| 亚洲av在线观看美女高潮| 免费看不卡的av| 丝袜美腿诱惑在线| 电影成人av| 久久久久视频综合| 高清欧美精品videossex| 韩国精品一区二区三区| 国产精品久久久av美女十八| 人妻 亚洲 视频| 欧美日本中文国产一区发布| 成人亚洲欧美一区二区av| 免费看十八禁软件| 久久久精品国产亚洲av高清涩受| 欧美中文综合在线视频| 亚洲 欧美一区二区三区| 女性被躁到高潮视频| 18禁观看日本| 欧美av亚洲av综合av国产av| 国产成人精品在线电影| 欧美精品一区二区免费开放| 亚洲 欧美一区二区三区| 中文字幕制服av| 黄色视频不卡| 日韩制服骚丝袜av| 91九色精品人成在线观看| 黄色视频在线播放观看不卡| 日韩一区二区三区影片| 秋霞在线观看毛片| 丰满饥渴人妻一区二区三| 另类亚洲欧美激情| 青草久久国产| 亚洲欧美日韩高清在线视频 | 国产精品.久久久| 国产国语露脸激情在线看| 午夜日韩欧美国产| 亚洲av成人不卡在线观看播放网 | 成人亚洲精品一区在线观看| 国产成人影院久久av| 99香蕉大伊视频| 中文精品一卡2卡3卡4更新| 免费在线观看日本一区| 久热爱精品视频在线9| 成人午夜精彩视频在线观看| 亚洲五月婷婷丁香| 真人做人爱边吃奶动态| 精品人妻1区二区| cao死你这个sao货| 搡老乐熟女国产| 国产日韩欧美视频二区| 久久毛片免费看一区二区三区| 国产亚洲欧美在线一区二区| 国产精品偷伦视频观看了| av片东京热男人的天堂| 少妇猛男粗大的猛烈进出视频| 一区二区三区激情视频| 国产精品香港三级国产av潘金莲 | 国语对白做爰xxxⅹ性视频网站| av在线老鸭窝| 久久99一区二区三区| svipshipincom国产片| 九色亚洲精品在线播放| 女人高潮潮喷娇喘18禁视频| 日韩,欧美,国产一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 中国国产av一级| 水蜜桃什么品种好| 久久狼人影院| 免费看十八禁软件| 黑人欧美特级aaaaaa片| 天天躁夜夜躁狠狠躁躁| 性色av乱码一区二区三区2| av一本久久久久| 麻豆乱淫一区二区| 久久久精品免费免费高清| 岛国毛片在线播放| 久久九九热精品免费| 欧美另类一区| 亚洲熟女毛片儿| 精品人妻熟女毛片av久久网站| 女警被强在线播放| 午夜久久久在线观看| 老司机影院成人| videos熟女内射| 在线看a的网站| www.熟女人妻精品国产| 曰老女人黄片| 热99久久久久精品小说推荐| av片东京热男人的天堂| 亚洲中文字幕日韩| 国产真人三级小视频在线观看| 如日韩欧美国产精品一区二区三区| 少妇猛男粗大的猛烈进出视频| 另类亚洲欧美激情| 国产一区二区在线观看av| 在线 av 中文字幕| 国产精品.久久久| 亚洲综合色网址| 中文精品一卡2卡3卡4更新| 男女高潮啪啪啪动态图| 少妇粗大呻吟视频| 国产成人精品久久二区二区91| xxxhd国产人妻xxx| 午夜福利视频在线观看免费| 最新的欧美精品一区二区| 久久国产亚洲av麻豆专区| 少妇被粗大的猛进出69影院| 波多野结衣一区麻豆| 亚洲自偷自拍图片 自拍| 欧美激情 高清一区二区三区| 亚洲伊人久久精品综合| 欧美日韩av久久| 久久精品国产亚洲av涩爱| 桃花免费在线播放| 亚洲专区国产一区二区| 人人妻人人爽人人添夜夜欢视频| 久久这里只有精品19| 777久久人妻少妇嫩草av网站| 成年人免费黄色播放视频| 成人国语在线视频| 99热全是精品| 90打野战视频偷拍视频| 一级片'在线观看视频| 一级黄片播放器| 黄色视频不卡| 日韩中文字幕视频在线看片| 侵犯人妻中文字幕一二三四区| 国产在线视频一区二区| 日韩大片免费观看网站| 亚洲国产精品成人久久小说| 热99久久久久精品小说推荐| 亚洲av成人精品一二三区| 热99国产精品久久久久久7| 黄片小视频在线播放| 极品人妻少妇av视频| a 毛片基地| a级片在线免费高清观看视频| 精品国产乱码久久久久久小说| 日本欧美国产在线视频| 亚洲男人天堂网一区| 丝袜美足系列| 亚洲精品久久成人aⅴ小说| 91国产中文字幕| 亚洲精品国产av成人精品| 亚洲成人国产一区在线观看 | av国产久精品久网站免费入址| 考比视频在线观看| 又大又黄又爽视频免费| 欧美日韩精品网址| 亚洲国产日韩一区二区| 亚洲成国产人片在线观看| 麻豆乱淫一区二区| 午夜福利在线免费观看网站| 国产精品免费大片| 亚洲 国产 在线| 美女大奶头黄色视频| 久久 成人 亚洲| 男人添女人高潮全过程视频| 99香蕉大伊视频| 一级a爱视频在线免费观看| av有码第一页| 欧美少妇被猛烈插入视频| 狠狠精品人妻久久久久久综合| 大型av网站在线播放| 国产成人精品久久久久久| 国产激情久久老熟女| 99久久人妻综合| 1024视频免费在线观看| 在线观看免费高清a一片| 国产不卡av网站在线观看| 日韩一卡2卡3卡4卡2021年| 亚洲熟女毛片儿| 天堂中文最新版在线下载| 啦啦啦中文免费视频观看日本| 国产xxxxx性猛交| 日本av免费视频播放| 欧美 亚洲 国产 日韩一| 亚洲第一青青草原| 日韩大片免费观看网站| 天堂俺去俺来也www色官网| 亚洲av电影在线进入| 成人影院久久| 午夜免费成人在线视频| 后天国语完整版免费观看| 欧美国产精品va在线观看不卡| 99re6热这里在线精品视频| 99热网站在线观看| 秋霞在线观看毛片| 成人亚洲精品一区在线观看| www.精华液| 韩国精品一区二区三区| 1024视频免费在线观看| 91九色精品人成在线观看| 深夜精品福利| 爱豆传媒免费全集在线观看| 中文字幕制服av| 男女边摸边吃奶| 夜夜骑夜夜射夜夜干| www.999成人在线观看| 国产亚洲欧美精品永久| 国产欧美日韩一区二区三区在线| 99国产精品一区二区蜜桃av | 亚洲精品av麻豆狂野| 午夜免费观看性视频| 又大又黄又爽视频免费| 久久国产精品影院| 国产男女超爽视频在线观看| 精品第一国产精品| 日韩视频在线欧美| av福利片在线| 亚洲精品第二区| 精品国产乱码久久久久久男人| 国产成人av激情在线播放| 欧美在线黄色| 亚洲国产av影院在线观看| 亚洲三区欧美一区| cao死你这个sao货| 欧美日韩国产mv在线观看视频| 最近中文字幕2019免费版| 香蕉国产在线看| 国产老妇伦熟女老妇高清| 日本黄色日本黄色录像| 国产男女内射视频| 亚洲 国产 在线| 欧美日韩成人在线一区二区| 老司机深夜福利视频在线观看 | 亚洲国产精品999| 国产男人的电影天堂91| 自线自在国产av| 久久99精品国语久久久| 一边摸一边做爽爽视频免费| 婷婷色麻豆天堂久久| 国产精品一区二区精品视频观看| 黑丝袜美女国产一区| 精品卡一卡二卡四卡免费| 日韩熟女老妇一区二区性免费视频| 久久国产亚洲av麻豆专区| 国产三级黄色录像| 日日摸夜夜添夜夜爱| 一二三四社区在线视频社区8| 久久影院123| 热99久久久久精品小说推荐| √禁漫天堂资源中文www| 人妻 亚洲 视频| 国产成人啪精品午夜网站| 国产一区有黄有色的免费视频| 中文字幕人妻熟女乱码| 国产亚洲一区二区精品| 日韩人妻精品一区2区三区| 欧美 亚洲 国产 日韩一| 亚洲九九香蕉| 午夜福利影视在线免费观看| 黄色一级大片看看| 久久国产精品人妻蜜桃| 日韩精品免费视频一区二区三区| 99久久精品国产亚洲精品| 丝袜在线中文字幕| 91国产中文字幕| 久久人人爽人人片av| 天天躁夜夜躁狠狠久久av| 久久精品亚洲熟妇少妇任你| 美女脱内裤让男人舔精品视频| 黄色毛片三级朝国网站| 国产成人精品在线电影| 国产1区2区3区精品| 亚洲国产日韩一区二区| 久久热在线av| 少妇被粗大的猛进出69影院| 777米奇影视久久| 国产又爽黄色视频| 国产精品人妻久久久影院| 欧美在线黄色| 亚洲欧洲日产国产| 另类亚洲欧美激情| 少妇的丰满在线观看| 国产野战对白在线观看| 成年女人毛片免费观看观看9 | 狂野欧美激情性bbbbbb| 中文字幕人妻丝袜制服| 我的亚洲天堂| 热99久久久久精品小说推荐| 自拍欧美九色日韩亚洲蝌蚪91| av天堂在线播放| 国产麻豆69| 欧美+亚洲+日韩+国产| 99精品久久久久人妻精品| 一区二区三区激情视频| 天天躁夜夜躁狠狠久久av| 国产精品久久久久成人av| 看十八女毛片水多多多| 日本av手机在线免费观看| 亚洲精品一二三| 国产精品 欧美亚洲| 日韩 欧美 亚洲 中文字幕| 日本vs欧美在线观看视频| 飞空精品影院首页| 午夜av观看不卡| a 毛片基地| 最近最新中文字幕大全免费视频 | 老汉色av国产亚洲站长工具| 搡老乐熟女国产| 久热这里只有精品99| 自拍欧美九色日韩亚洲蝌蚪91| 精品人妻1区二区| e午夜精品久久久久久久| 国产在线观看jvid| 丰满人妻熟妇乱又伦精品不卡| av欧美777| www.精华液| 精品人妻在线不人妻| www.999成人在线观看| 成人三级做爰电影| 日韩电影二区| 熟女少妇亚洲综合色aaa.| 日本欧美视频一区| 在线av久久热| 国产在线观看jvid| 爱豆传媒免费全集在线观看| 欧美精品高潮呻吟av久久| 女人精品久久久久毛片| 精品亚洲乱码少妇综合久久| 国产极品粉嫩免费观看在线| 国产高清视频在线播放一区 | 亚洲av欧美aⅴ国产| 久久久久网色| 免费观看av网站的网址| 夜夜骑夜夜射夜夜干| 国产成人精品久久久久久| 人妻人人澡人人爽人人| 人体艺术视频欧美日本| av国产久精品久网站免费入址| 色94色欧美一区二区| 9191精品国产免费久久| 亚洲av电影在线进入| 久久精品国产综合久久久| 欧美日韩亚洲综合一区二区三区_| 夜夜骑夜夜射夜夜干| 男女免费视频国产| 一级a爱视频在线免费观看| 亚洲精品日韩在线中文字幕| 亚洲专区国产一区二区| 久久精品成人免费网站| 一区二区三区激情视频| 最新的欧美精品一区二区| 秋霞在线观看毛片| 18在线观看网站| 热re99久久精品国产66热6| 精品熟女少妇八av免费久了| 精品人妻1区二区| 狂野欧美激情性xxxx| 一区福利在线观看| 国产精品一区二区在线不卡| 深夜精品福利| 国产免费一区二区三区四区乱码| 香蕉丝袜av| 男人爽女人下面视频在线观看| 免费高清在线观看视频在线观看| 欧美亚洲日本最大视频资源| 国产一区二区在线观看av| 国产精品麻豆人妻色哟哟久久| 又粗又硬又长又爽又黄的视频| 午夜福利乱码中文字幕| 亚洲精品久久午夜乱码| 免费高清在线观看日韩| 成年人午夜在线观看视频| 免费高清在线观看视频在线观看| 亚洲专区中文字幕在线| 成人国语在线视频| 久久狼人影院| 制服人妻中文乱码| 岛国毛片在线播放| 女人爽到高潮嗷嗷叫在线视频| 国产成人精品久久久久久| 手机成人av网站| 免费高清在线观看视频在线观看| 美女视频免费永久观看网站| 久久精品久久久久久噜噜老黄| 国产亚洲精品第一综合不卡|