• <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ā)展動能
    国产久久久一区二区三区| 色综合色国产| av免费观看日本| 啦啦啦在线观看免费高清www| 久久99热这里只频精品6学生| 国产精品99久久99久久久不卡 | 三级国产精品欧美在线观看| 亚洲欧美成人综合另类久久久| 亚洲精品国产成人久久av| 在现免费观看毛片| 国产一区有黄有色的免费视频| 成年版毛片免费区| 精品人妻视频免费看| 欧美少妇被猛烈插入视频| 日本欧美国产在线视频| 国产日韩欧美亚洲二区| 亚洲经典国产精华液单| 高清日韩中文字幕在线| 日本一二三区视频观看| 在线观看一区二区三区激情| 久久久久久久久久人人人人人人| 国产伦精品一区二区三区视频9| 少妇裸体淫交视频免费看高清| 在线免费观看不下载黄p国产| 大又大粗又爽又黄少妇毛片口| 人体艺术视频欧美日本| 最近最新中文字幕免费大全7| 久久韩国三级中文字幕| 中文字幕免费在线视频6| 伦精品一区二区三区| 亚洲精品亚洲一区二区| 欧美日韩精品成人综合77777| 久久久成人免费电影| 日本wwww免费看| 日日撸夜夜添| 日韩强制内射视频| 一级毛片久久久久久久久女| 亚洲欧洲日产国产| 日韩大片免费观看网站| 亚洲自偷自拍三级| 噜噜噜噜噜久久久久久91| 国产精品久久久久久av不卡| 一级毛片电影观看| 国产精品三级大全| 99热6这里只有精品| 老女人水多毛片| 日韩一区二区视频免费看| 麻豆成人av视频| 赤兔流量卡办理| 日韩欧美一区视频在线观看 | 在线观看三级黄色| 久久久久久国产a免费观看| 2018国产大陆天天弄谢| 久久久久久国产a免费观看| 日韩强制内射视频| 噜噜噜噜噜久久久久久91| 欧美日韩视频高清一区二区三区二| 国产精品爽爽va在线观看网站| 毛片一级片免费看久久久久| 国产 一区 欧美 日韩| 99精国产麻豆久久婷婷| 免费电影在线观看免费观看| 五月天丁香电影| 国产一区有黄有色的免费视频| 夜夜看夜夜爽夜夜摸| 一级毛片 在线播放| 久久午夜福利片| 丝瓜视频免费看黄片| 成人免费观看视频高清| 色吧在线观看| 国产精品久久久久久精品电影小说 | 国产精品99久久99久久久不卡 | 国产黄片视频在线免费观看| 97在线人人人人妻| 麻豆久久精品国产亚洲av| 一本色道久久久久久精品综合| 18禁裸乳无遮挡免费网站照片| 久久久久久久久久久丰满| 日本午夜av视频| 日韩免费高清中文字幕av| 亚洲精品日韩在线中文字幕| 免费人成在线观看视频色| 国产黄色视频一区二区在线观看| 男人爽女人下面视频在线观看| 51国产日韩欧美| 国产精品一区www在线观看| 观看美女的网站| 又大又黄又爽视频免费| 搡女人真爽免费视频火全软件| 成人漫画全彩无遮挡| 亚洲国产精品成人久久小说| 国产一区有黄有色的免费视频| 人人妻人人澡人人爽人人夜夜| 精品国产乱码久久久久久小说| 午夜免费鲁丝| 内地一区二区视频在线| 国产精品久久久久久久电影| 国产成年人精品一区二区| 18禁裸乳无遮挡免费网站照片| av黄色大香蕉| 在线看a的网站| 少妇人妻久久综合中文| 国产精品熟女久久久久浪| 国产高清三级在线| 精品人妻偷拍中文字幕| 一级毛片久久久久久久久女| 高清毛片免费看| 亚洲av二区三区四区| 亚洲av成人精品一二三区| 一二三四中文在线观看免费高清| 日日啪夜夜爽| 两个人的视频大全免费| 嫩草影院精品99| 麻豆精品久久久久久蜜桃| 免费av观看视频| 午夜福利在线在线| 丝袜美腿在线中文| av在线老鸭窝| 国产乱人视频| 日韩成人av中文字幕在线观看| 成人二区视频| 亚洲成色77777| 亚洲av成人精品一二三区| 亚洲国产最新在线播放| 日韩电影二区| 亚洲自拍偷在线| 久久久久久久精品精品| 亚洲国产精品国产精品| 男的添女的下面高潮视频| 精品久久久久久久久av| 欧美3d第一页| av天堂中文字幕网| 国产国拍精品亚洲av在线观看| 亚洲欧美中文字幕日韩二区| 久久久久久久久久成人| 亚洲伊人久久精品综合| 国产 精品1| 国产成人午夜福利电影在线观看| 亚洲精品日韩av片在线观看| 一级毛片我不卡| 免费黄频网站在线观看国产| 亚洲精品自拍成人| 亚洲精品久久午夜乱码| 观看免费一级毛片| 亚洲精品亚洲一区二区| 六月丁香七月| 美女被艹到高潮喷水动态| 日韩一本色道免费dvd| 久久久精品欧美日韩精品| 啦啦啦啦在线视频资源| 最近中文字幕2019免费版| 啦啦啦在线观看免费高清www| 五月天丁香电影| 国产又色又爽无遮挡免| 看非洲黑人一级黄片| 大片免费播放器 马上看| videos熟女内射| 国产免费福利视频在线观看| 偷拍熟女少妇极品色| 18禁动态无遮挡网站| 国产精品.久久久| 99热网站在线观看| 51国产日韩欧美| 一个人看视频在线观看www免费| 日韩伦理黄色片| 日韩欧美精品免费久久| 欧美高清成人免费视频www| 国产一区二区三区av在线| 亚洲高清免费不卡视频| 80岁老熟妇乱子伦牲交| 黄片wwwwww| videossex国产| 美女被艹到高潮喷水动态| 亚洲欧美精品专区久久| 在线观看免费高清a一片| 1000部很黄的大片| 国产高清有码在线观看视频| 午夜精品一区二区三区免费看| 日韩av不卡免费在线播放| 欧美一区二区亚洲| 老师上课跳d突然被开到最大视频| 边亲边吃奶的免费视频| 一个人观看的视频www高清免费观看| 内地一区二区视频在线| 男的添女的下面高潮视频| 午夜亚洲福利在线播放| 国产日韩欧美在线精品| 亚洲天堂av无毛| 91精品伊人久久大香线蕉| 久久精品国产鲁丝片午夜精品| 国产午夜福利久久久久久| 一级毛片电影观看| 精品人妻一区二区三区麻豆| 汤姆久久久久久久影院中文字幕| 午夜福利视频1000在线观看| 特级一级黄色大片| 国产成人精品婷婷| 成人亚洲精品av一区二区| 高清午夜精品一区二区三区| 校园人妻丝袜中文字幕| 亚洲国产精品国产精品| 国产黄频视频在线观看| 毛片女人毛片| 国产探花极品一区二区| 亚洲性久久影院| 看非洲黑人一级黄片| 在现免费观看毛片| 五月开心婷婷网| 国产国拍精品亚洲av在线观看| 亚洲精品中文字幕在线视频 | 午夜免费男女啪啪视频观看| www.色视频.com| 熟女人妻精品中文字幕| 波野结衣二区三区在线| 偷拍熟女少妇极品色| videos熟女内射| 午夜日本视频在线| 99久久精品国产国产毛片| 日日摸夜夜添夜夜添av毛片| 免费不卡的大黄色大毛片视频在线观看| av女优亚洲男人天堂| 欧美人与善性xxx| 日日撸夜夜添| 少妇熟女欧美另类| 国产成人午夜福利电影在线观看| 亚洲av.av天堂| 欧美成人一区二区免费高清观看| 精品人妻一区二区三区麻豆| 精华霜和精华液先用哪个| 中文乱码字字幕精品一区二区三区| 国产毛片在线视频| 一级片'在线观看视频| 中文乱码字字幕精品一区二区三区| 日韩国内少妇激情av| 亚洲国产欧美在线一区| 免费黄网站久久成人精品| 人妻夜夜爽99麻豆av| 欧美日韩视频高清一区二区三区二| 观看免费一级毛片| 亚洲av成人精品一二三区| 秋霞在线观看毛片| 亚洲电影在线观看av| 我要看日韩黄色一级片| 亚洲人成网站高清观看| 亚洲欧洲国产日韩| 蜜桃亚洲精品一区二区三区| 亚洲一区二区三区欧美精品 | 国内少妇人妻偷人精品xxx网站| 欧美zozozo另类| 美女被艹到高潮喷水动态| 国产精品秋霞免费鲁丝片| 亚洲人与动物交配视频| 我的老师免费观看完整版| 亚洲欧美日韩卡通动漫| 午夜免费观看性视频| 亚洲欧美成人综合另类久久久| 久久人人爽人人片av| .国产精品久久| 国产男人的电影天堂91| 国产乱人偷精品视频| 欧美亚洲 丝袜 人妻 在线| 一本久久精品| 岛国毛片在线播放| 日本爱情动作片www.在线观看| 99久久精品一区二区三区| 亚洲国产精品专区欧美| 亚洲av不卡在线观看| 国内少妇人妻偷人精品xxx网站| 国产免费一区二区三区四区乱码| 亚洲国产精品999| 国产一区亚洲一区在线观看| 亚洲av在线观看美女高潮| 国产极品天堂在线| 国产精品一区www在线观看| 人体艺术视频欧美日本| 色播亚洲综合网| 人人妻人人看人人澡| 亚洲自偷自拍三级| 真实男女啪啪啪动态图| 日韩中字成人| 成人鲁丝片一二三区免费| 国产女主播在线喷水免费视频网站| 亚洲精品影视一区二区三区av| 伦精品一区二区三区| 国产一区二区亚洲精品在线观看| 日韩视频在线欧美| 亚洲自偷自拍三级| 国产精品麻豆人妻色哟哟久久| 日韩中字成人| 国产v大片淫在线免费观看| 国产黄色免费在线视频| 亚洲国产色片| 亚洲欧美精品自产自拍| 国产亚洲精品久久久com| 亚洲成色77777| 黄色怎么调成土黄色| 80岁老熟妇乱子伦牲交| 色哟哟·www| 国产日韩欧美亚洲二区| 久久99热这里只有精品18| 蜜桃亚洲精品一区二区三区| 国产 一区精品| 日韩 亚洲 欧美在线| 亚洲欧美一区二区三区黑人 | 嫩草影院精品99| 日韩电影二区| 国产成人福利小说| 日韩视频在线欧美| 久久久国产一区二区| 国产乱人偷精品视频| 久久精品国产自在天天线| 啦啦啦啦在线视频资源| 欧美日韩在线观看h| 国产 精品1| 午夜福利在线观看免费完整高清在| 毛片一级片免费看久久久久| 白带黄色成豆腐渣| 亚洲第一区二区三区不卡| 国产亚洲av嫩草精品影院| 国产黄片视频在线免费观看| 嫩草影院新地址| 国产日韩欧美亚洲二区| 亚洲一区二区三区欧美精品 | 国产伦在线观看视频一区| 男女下面进入的视频免费午夜| 一级黄片播放器| 九草在线视频观看| 日本午夜av视频| 丰满人妻一区二区三区视频av| 国产成年人精品一区二区| www.色视频.com| 亚洲欧美日韩东京热| 欧美+日韩+精品| 日韩视频在线欧美| 51国产日韩欧美| 91精品伊人久久大香线蕉| 亚洲最大成人手机在线| 99热这里只有精品一区| 午夜视频国产福利| 国产精品一区二区在线观看99| 只有这里有精品99| 18禁在线播放成人免费| 晚上一个人看的免费电影| 91精品一卡2卡3卡4卡| 亚洲精品久久久久久婷婷小说| 激情 狠狠 欧美| 麻豆成人午夜福利视频| 国产在线男女| 国产极品天堂在线| 狂野欧美激情性xxxx在线观看| 99热全是精品| 国产成人精品福利久久| 亚洲国产欧美在线一区| 免费黄色在线免费观看| 日本熟妇午夜| 免费电影在线观看免费观看| 毛片女人毛片| 18禁动态无遮挡网站| 99久久人妻综合| 99热全是精品| 一级av片app| 新久久久久国产一级毛片| 久久久久性生活片| 久久久色成人| 亚洲四区av| 久久久久久久国产电影| 亚洲欧美成人精品一区二区| 纵有疾风起免费观看全集完整版| 一二三四中文在线观看免费高清| 亚洲经典国产精华液单| 午夜福利视频精品| 精品国产一区二区三区久久久樱花 | 国产精品蜜桃在线观看| 黄片wwwwww| 国产精品久久久久久久久免| 亚洲欧美日韩无卡精品| 国产精品久久久久久久久免| 男女下面进入的视频免费午夜| 久久久久久久国产电影| 免费观看a级毛片全部| 69人妻影院| 亚洲av不卡在线观看| 欧美日本视频| 国产精品一区二区性色av| 在线a可以看的网站| 黄色一级大片看看| 日韩欧美 国产精品| 日本黄大片高清| 免费大片黄手机在线观看| 国产又色又爽无遮挡免| 国产亚洲5aaaaa淫片| 国产黄片视频在线免费观看| 性色av一级| 成人高潮视频无遮挡免费网站| 精品一区在线观看国产| 亚洲无线观看免费| 国产亚洲最大av| 亚洲人成网站在线观看播放| 舔av片在线| 高清视频免费观看一区二区| 亚洲aⅴ乱码一区二区在线播放| 亚洲色图av天堂| 少妇猛男粗大的猛烈进出视频 | 好男人视频免费观看在线| 国产老妇女一区| 高清欧美精品videossex| 伊人久久精品亚洲午夜| 精品国产乱码久久久久久小说| 91久久精品国产一区二区成人| 亚洲精品456在线播放app| 精品久久久久久久人妻蜜臀av| 日韩一区二区视频免费看| 国产精品人妻久久久久久| 真实男女啪啪啪动态图| 国产久久久一区二区三区| 久久精品国产亚洲av涩爱| 免费播放大片免费观看视频在线观看| 99久国产av精品国产电影| 日韩电影二区| 欧美日本视频| 老司机影院毛片| 国产男女内射视频| 国产中年淑女户外野战色| 国产探花极品一区二区| 美女主播在线视频| 中文乱码字字幕精品一区二区三区| 国产伦理片在线播放av一区| 精品一区二区免费观看| 秋霞在线观看毛片| 成年女人看的毛片在线观看| 观看美女的网站| 国产在视频线精品| 伦精品一区二区三区| 777米奇影视久久| 午夜精品国产一区二区电影 | 亚洲av.av天堂| 国产v大片淫在线免费观看| 国产永久视频网站| 成人无遮挡网站| 亚洲国产精品999| 国产精品人妻久久久影院| 啦啦啦啦在线视频资源| 波野结衣二区三区在线| 精品久久久久久电影网| 91久久精品国产一区二区成人| 91aial.com中文字幕在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品久久久久久久电影| 成年人午夜在线观看视频| 国内少妇人妻偷人精品xxx网站| 日韩一区二区三区影片| 身体一侧抽搐| 六月丁香七月| 国产伦在线观看视频一区| 五月玫瑰六月丁香| 久久人人爽人人片av| 国产成人一区二区在线| 三级国产精品欧美在线观看| 亚洲精品久久久久久婷婷小说| 一级片'在线观看视频| 一本久久精品| 97精品久久久久久久久久精品| av线在线观看网站| 大又大粗又爽又黄少妇毛片口| 交换朋友夫妻互换小说| 女人十人毛片免费观看3o分钟| 在线观看人妻少妇| 男女边摸边吃奶| 午夜老司机福利剧场| 制服丝袜香蕉在线| 国产视频内射| 搡女人真爽免费视频火全软件| 国产综合懂色| 噜噜噜噜噜久久久久久91| 国产欧美亚洲国产| 日韩欧美精品免费久久| 国产永久视频网站| 中文字幕亚洲精品专区| 欧美精品国产亚洲| 久久亚洲国产成人精品v| 在线观看免费高清a一片| 天堂中文最新版在线下载 | 免费黄频网站在线观看国产| 草草在线视频免费看| 国产精品久久久久久精品电影小说 | 日本wwww免费看| 亚洲四区av| 亚洲在久久综合| 国产精品国产三级国产av玫瑰| 久久久久网色| 人人妻人人爽人人添夜夜欢视频 | 成人漫画全彩无遮挡| 成人美女网站在线观看视频| 看非洲黑人一级黄片| 丝袜脚勾引网站| 国内精品美女久久久久久| 亚洲人与动物交配视频| 男女下面进入的视频免费午夜| 国产色婷婷99| 久久综合国产亚洲精品| 在线a可以看的网站| 国产黄色视频一区二区在线观看| 亚洲成人中文字幕在线播放| 观看免费一级毛片| 18禁在线播放成人免费| 在线a可以看的网站| 狂野欧美白嫩少妇大欣赏| 日本午夜av视频| 国产日韩欧美在线精品| 午夜亚洲福利在线播放| 日本免费在线观看一区| 日韩免费高清中文字幕av| 国产av码专区亚洲av| 好男人视频免费观看在线| 日韩国内少妇激情av| 欧美国产精品一级二级三级 | 日韩av不卡免费在线播放| 男男h啪啪无遮挡| 久久久欧美国产精品| 亚洲国产欧美在线一区| 中文在线观看免费www的网站| 免费高清在线观看视频在线观看| 免费观看的影片在线观看| 色综合色国产| 男男h啪啪无遮挡| h日本视频在线播放| 色播亚洲综合网| 亚洲人成网站在线观看播放| 国产成人精品婷婷| 老女人水多毛片| 亚洲熟女精品中文字幕| 黄色视频在线播放观看不卡| 精品酒店卫生间| 久久综合国产亚洲精品| 亚洲av中文字字幕乱码综合| 国产精品久久久久久精品电影| 人妻系列 视频| 在线a可以看的网站| 国精品久久久久久国模美| 欧美老熟妇乱子伦牲交| 校园人妻丝袜中文字幕| 亚洲,欧美,日韩| 听说在线观看完整版免费高清| 尾随美女入室| 一级二级三级毛片免费看| 午夜精品一区二区三区免费看| 国产精品成人在线| 国产在视频线精品| 国产黄片美女视频| 欧美性感艳星| 精品一区二区三卡| 草草在线视频免费看| 欧美 日韩 精品 国产| 哪个播放器可以免费观看大片| 免费看日本二区| 欧美日韩在线观看h| 男人狂女人下面高潮的视频| av网站免费在线观看视频| 亚洲国产最新在线播放| 亚洲成色77777| 18禁裸乳无遮挡免费网站照片| 欧美激情国产日韩精品一区| 一区二区av电影网| 丝瓜视频免费看黄片| av国产免费在线观看| 免费黄网站久久成人精品| 深夜a级毛片| 日韩一区二区三区影片| 国产黄频视频在线观看| 久久久久久久国产电影| 免费播放大片免费观看视频在线观看| 中文在线观看免费www的网站| 97超视频在线观看视频| 99热这里只有是精品50| 久久久久久国产a免费观看| 免费av观看视频| 一本久久精品| 国产精品久久久久久久电影| 免费观看的影片在线观看| 一级毛片黄色毛片免费观看视频| 欧美xxxx黑人xx丫x性爽| 国产精品.久久久| 91久久精品国产一区二区成人| 黄色配什么色好看| 国产黄片视频在线免费观看| 亚洲在久久综合| 免费观看性生交大片5| 国产精品99久久久久久久久| 激情 狠狠 欧美| 尾随美女入室| 国产综合精华液| 欧美日本视频| 久久人人爽av亚洲精品天堂 | 超碰97精品在线观看| 神马国产精品三级电影在线观看| 国产美女午夜福利| 亚洲精品乱码久久久久久按摩| 简卡轻食公司| 亚洲精品视频女| 久久99热6这里只有精品| 亚洲成人精品中文字幕电影| 嘟嘟电影网在线观看| 3wmmmm亚洲av在线观看| 女的被弄到高潮叫床怎么办| 少妇被粗大猛烈的视频| 成年人午夜在线观看视频| 三级男女做爰猛烈吃奶摸视频| 午夜老司机福利剧场| av.在线天堂| 国产老妇伦熟女老妇高清| 午夜福利高清视频| 免费在线观看成人毛片| 欧美亚洲 丝袜 人妻 在线| 欧美丝袜亚洲另类| 亚洲av中文av极速乱|