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

    基于ULS、TLS和超聲測高儀的天然次生林中不同林冠層樹高估測

    2021-07-28 04:07:32趙穎慧楊海城
    關(guān)鍵詞:針葉樹單木闊葉樹

    趙穎慧,楊海城,甄 貞*

    (1.東北林業(yè)大學(xué)林學(xué)院,黑龍江 哈爾濱 150040;2.東北林業(yè)大學(xué)森林生態(tài)系統(tǒng)可持續(xù)經(jīng)營教育部重點實驗室,黑龍江 哈爾濱 150040)

    樹高是森林資源調(diào)查中的重要因子之一,與胸徑(DBH)和樹種等其他因子共同被廣泛用于估測無法以無損方式獲得的生物量等重要森林參數(shù)。在生物量模型構(gòu)建中加入樹高因子作為預(yù)測參數(shù),可以顯著地提高生物量模型的估計精度[1-2],因此,樹高因子測定的準確性至關(guān)重要。目前,樹高的實測方法通常是應(yīng)用測高儀進行實地測量,不僅費時費力,還受林況、測量者的經(jīng)驗及所使用的設(shè)備等多種因素影響[3]。

    激光雷達(light detection and ranging, LiDAR)是一種主動遙感技術(shù),其發(fā)射的激光脈沖在森林中具有很好的穿透性,高采樣密度LiDAR能夠獲取單木三維結(jié)構(gòu)特征,因此被廣泛用于林業(yè)研究[4]。無人機激光雷達數(shù)據(jù)(unmanned aerial vehicle laser scanning, ULS)是近年來新興的遙感數(shù)據(jù),與機載激光雷達數(shù)據(jù)(airborne laser scanning, ALS)相比,具有更低的成本和更高的點云密度[5],在森林冠層結(jié)構(gòu)估測中具有明顯的優(yōu)勢[6-7]。地基激光雷達(terrestrial laser scanning,TLS)生成的3D點云,能夠獲得高質(zhì)量的林下LiDAR影像,可準確地描述林下信息,用于精確提取林下結(jié)構(gòu)[8-10]。有學(xué)者使用TLS提取樹高并進行分析,認為TLS提取的樹高會低于實際值[11-13]。也有學(xué)者認為使用ALS提取的樹高,也會受點密度和森林結(jié)構(gòu)等影響而造成低估[14-15]。另外,Sibona等[16]測量了100株砍伐木樹高,并對砍伐木樹高與野外調(diào)查中使用測距儀獲得的樹高和使用ALS提取的樹高進行了比較,結(jié)果顯示ALS提取的樹高更接近真實值。Wang等[17]比較了ULS和TLS提取的樹高,結(jié)果表明ULS在樹高的估測方面優(yōu)于TLS。Wallace等[18]直接使用ULS點云生成冠層高度模型(canopy height model, CHM),使用標記控制分水嶺算法單木分割提取樹高,結(jié)果表明樹高被低估了0.11 m。由此可見,基于單一數(shù)據(jù)源(ALS、ULS和TLS)提取樹高都有一定的局限性。

    有學(xué)者將林分按照一定的樹木高度為分割點,分別采用不同數(shù)據(jù)源進行估測。例如,Wang等[19]分別用TLS和ALS提取的樹高與樣地實測樹高對比,結(jié)果表明:當(dāng)樹高小于15 m時,TLS提取的樹高與樣地實測樹高更接近;當(dāng)樹高大于15 m時,使用ALS提取的樹高更準確,但該樹高分割點(15 m)的選取比較主觀,僅適合應(yīng)用于芬蘭北部的松樹和樺樹混交林。另外,也有研究者認為樹高測量的誤差主要來自林冠層的遮蔽[20]。林冠層垂直結(jié)構(gòu)信息豐富,體現(xiàn)了森林植被在垂直方向上的層次性和空間配置方式,影響著樹木生長和下層植被分布[21-23]。胡文杰等[24]根據(jù)樣地調(diào)查數(shù)據(jù),采用樹冠光競爭高度法將林分劃分為上林層、中林層和下林層。胡傳偉等[25]通過森林群落學(xué)調(diào)查的方法獲取數(shù)據(jù),應(yīng)用聚類分析和生態(tài)位分析相結(jié)合的方法,研究了天然次生林的樹種組成與垂直結(jié)構(gòu)。對于激光雷達在林冠分層上的應(yīng)用,也有學(xué)者進行了研究[26-30]。有研究將冠層高度分布(camopy height distribution, CHD)用于對林分的冠層劃分,Maltamo等[31]使用改進的直方圖閾值法對CHD直方圖進行計算,實現(xiàn)了對冠層的分層。

    綜上所述,以ALS、ULS和TLS單一數(shù)據(jù)源提取樹高的研究已十分常見,但大多數(shù)只是針對某一林分直接提取樹高,而對林分以CHD進行分層,并根據(jù)不同冠層結(jié)合不同數(shù)據(jù)源提取樹高的對比研究較少。因此,本研究以ULS、TLS和野外樣地調(diào)查實測樹高為數(shù)據(jù)源,以東北林業(yè)大學(xué)帽兒山實驗林場中林施業(yè)區(qū)50 m×50 m(0.25 hm2)樣地為研究對象,基于CHD對林冠層進行分層,針對不同樹木類型(針葉樹和闊葉樹)和不同林冠層(上層和下層)的單木進行樹高估測,分析在不同林冠層中應(yīng)用不同數(shù)據(jù)源探測提取樹高的異同,探索適用于中國北方天然次生林樹高估測的方法,為更準確地估測天然次生林樹高提供依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    研究區(qū)為東北林業(yè)大學(xué)帽兒山實驗林場中林施業(yè)區(qū)50 m×50 m(0.25 hm2)方形樣地(127°34′15″E,45°21′29″N)。林場位于黑龍江省尚志市西北部,地處張廣才嶺西坡,屬于大陸性季風(fēng)氣候,植被是典型的東北天然次生林,原始地帶性頂級群落為紅松闊葉林。由于森林屢遭破壞,原始植被發(fā)生了逆向演替,形成了目前以珍貴闊葉林、楊樺林、柞木林等為主的天然次生林與紅松(Pinuskoraiensis)、落葉松(Larixgmelinii)、樟子松(Pinussylvestrisvar.mongolica)等人工林鑲嵌分布的森林景觀。

    1.2 研究數(shù)據(jù)獲取

    1.2.1 ULS數(shù)據(jù)

    數(shù)據(jù)獲取時間為2019年7月,使用飛馬D200旋翼平臺的一體化高精度航測無人機系統(tǒng),配備高精度LiDAR模塊(D-lidar200)。LiDAR傳感器(RIEGL mini VUX-1UAV)波長為905 nm,最大回波數(shù)為5,最大傳感器頻率為100 Hz,點云密度為104點/m2。

    1.2.2 TLS數(shù)據(jù)

    數(shù)據(jù)獲取時間為2019年10月,由麥格天泓公司開發(fā)的Trimble TX6高速三維激光掃描儀獲得,在樣地(50 m × 50 m)內(nèi)采用蛇形方式布設(shè)多個測點,原則上確保每一株樹都能360°接收到激光掃描儀掃描。高速三維激光掃描儀角度精度80 μrad,掃描速度50萬點/s,最小測程0.6 m,最大測程80 m,測量誤差±2 mm。

    1.2.3 野外調(diào)查樣地數(shù)據(jù)

    2018年7月在帽兒山林場中林施業(yè)區(qū)內(nèi)選擇包含多個優(yōu)勢樹種的50 m×50 m方形樣地,郁閉度為0.8,并于2019年7月復(fù)測。使用RTK(real time kinematic)記錄樣地中所有樹木冠頂(樹冠最高點)在地面上的投影坐標位置(誤差小于0.1 m)。以5 cm為起測徑階進行每木檢尺,記錄包括:胸徑、樹高、冠幅和樹種等因子。使用瑞典Hagolf公司生產(chǎn)的Vertex IV超聲測高儀獲取樹高。樣地中共記錄359株樣木,其中針葉樹144株,主要為紅松;闊葉樹215株,主要為榆樹(Ulmuspumila)、胡桃楸(Juglansmandshurica)、楊樹(PopulusL.)、水曲柳(Fraxinusmandshurica)和白樺(Betulaplatyphylla)等。樣地中不同樹木類型(針葉樹和闊葉樹)的胸徑、樹高和冠幅的統(tǒng)計特征見表1。

    表1 樣地樹木統(tǒng)計特征

    1.3 研究方法

    1.3.1 數(shù)據(jù)預(yù)處理

    ULS點云數(shù)據(jù)預(yù)處理包括:去噪、分離地面點和點云歸一化處理。首先對原始的點云數(shù)據(jù)去噪,之后采用漸進加密三角網(wǎng)算法進行濾波[32],獲得地面點和非地面點。最后,使用地面點生成的數(shù)字高程模型(DEM)對點云進行歸一化處理,得到點云的高程值,若該點為樹頂,則其高程值可視為樹高。

    TLS點云數(shù)據(jù)預(yù)處理包括:去除“飛點”、抽稀、分離地面點和歸一化處理。首先去除“飛點”,即噪點;其次,對點云數(shù)據(jù)采用八叉樹方法[33]進行“抽稀”,通過多次抽稀研究發(fā)現(xiàn),采用13級八叉樹抽稀后獲得的數(shù)據(jù)運算量小又極大地保持了數(shù)據(jù)真實性;最后,采用漸進加密三角網(wǎng)算法對點云進行地面點濾波,然后做歸一化處理。

    1.3.2 基于CHD的林冠分層

    1.3.3 單木樹高的提取

    本研究目的是探討不同林冠層中應(yīng)用不同數(shù)據(jù)源提取樹高的異同,為排除單木樹冠提取算法對結(jié)果的影響,本研究對ULS和TLS數(shù)據(jù)進行人工目視解譯,分離出單株樹并手動量取樹高[19]。當(dāng)使用ULS提取樹高時,有樹高交叉重疊現(xiàn)象:①如果交叉重疊不嚴重,本研究結(jié)合樣地調(diào)查的單木位置,直接提取樹高,計算誤差。②如果交叉重疊嚴重,頂層的點云數(shù)據(jù)是根據(jù)點云位置的最高點與實測樹高的位置匹配,提取樹高;下層的點云則根據(jù)單木冠型來判斷其樹頂點提取樹高,如果冠型難以識別,則記為識別失敗。

    1.3.4 樹高離群值的剔除和比較評價

    在對ULS和TLS數(shù)據(jù)探測提取的樹高與實測樹高進行比較之前,首先要剔除離群值。離群值表明如果3種數(shù)據(jù)源(ULS、TLS及ULS和TLS的交集)提取的樹高之間存在明顯差異,則存在測量錯誤,因此進行比較評價之前要對離群值進行分析[19]。由于樹高交叉重疊現(xiàn)象,在單木樹高提取中被錯誤識別的樹高也會因判定為離群值而被剔除。離群值的判斷方法如下:

    (1)

    (2)

    剔除離群值后,將林冠層分上、下層,分別將ULS和TLS探測提取的單木樹高與樣地實測的單木樹高進行比較,評價指標包括均方根誤差[RMSE,式中記為σ(RMSE)],相對均方根誤差[rRMSE,式中記為σ(rRMSE)]、偏離率[Bias,式中記為σ(Bias)]和相對偏離率[rBias,式中記為σ(rBias)][36],計算公式分別為:

    (3)

    (4)

    (5)

    (6)

    2 結(jié)果與分析

    2.1 基于CHD曲線的林冠分層

    對ULS點云切片后,以點云百分比為橫坐標,點云高度為縱坐標,繪制CHD曲線如圖1A所示。由Lloyd算法計算出的林冠上、下層之間的高度閾值為8.5 m。在CHD中,冠層的相對分布密度被描述為高度的函數(shù)[37]。結(jié)合由圖1B中的點云側(cè)視圖可以看出,在高度為0 m的位置存在少量點云,該區(qū)域可視為地面部分。在8 m以下的高度,點云百分比始終維持在較低的范圍,主要由于樣地中雖然存在大量低矮樹木,但也受到冠層的遮蔽,難以接收到ULS的掃描。隨著高度的增加,點云百分比逐漸增加,當(dāng)高度到達8.5 m時,點云數(shù)量急劇增大,說明此時樹木冠層更加密集,能接收到的點云數(shù)量增加。在9~12 m處為點云百分比較大的位置,結(jié)合ULS的掃描特點,可以認為大部分樹木的枝葉部分都集中在這個高度上,所以能夠獲得大量的點云回波,相應(yīng)地也會對該高度以下的位置有較強的遮蔽作用,在12 m以上,點云百分比逐漸減小,可以認為隨著冠層高度增加樹木逐漸稀疏,很少受到遮擋或不受其他樹木遮擋。

    圖1 CHD曲線(A)和點云視圖(B)

    2.2 離群值的剔除

    按公式(1)、(2)對ULS和TLS探測提取的樹高離群值進行了剔除和比較,結(jié)果如表2所示。

    表2 3種數(shù)據(jù)源的離群值統(tǒng)計

    基于ULS數(shù)據(jù)的離群值共有30個,其中林冠上層有20個,林冠下層有10個。從樹木類型來看,產(chǎn)生離群值的有19株闊葉樹、11株針葉樹。

    基于TLS數(shù)據(jù)的離群值共有11個,主要出現(xiàn)在林冠上層,在林冠下層僅有2株單木被判定為離群值。從樹木類型來看,產(chǎn)生離群值的主要是闊葉樹,針葉樹的離群值共有2個。與實測數(shù)據(jù)對比發(fā)現(xiàn),在11個離群值中,有6個離群值出現(xiàn)在9~14 m的高度上,由圖1的CHD曲線圖可以發(fā)現(xiàn),在這個高度上是點云百分比較大的位置,林冠層枝葉密集,容易產(chǎn)生誤判。此外,TLS數(shù)據(jù)對2株單木的樹高產(chǎn)生了高估,這2株單木都是闊葉樹,高度位于9~12 m的枝葉密度較大的林冠層中。

    ULS數(shù)據(jù)與TLS數(shù)據(jù)產(chǎn)生的離群值交集共有3個,都是9~11 m林冠層的闊葉樹。

    2.3 基于不同數(shù)據(jù)源的樹高比較

    2.3.1 基于ULS的估測樹高與實測樹高比較

    去掉ULS數(shù)據(jù)探測提取樹高的離群值后,共有245株單木,其中林冠上、下層分別有單木230、15株,與樣地實測樹高比較結(jié)果見表3。在ULS探測的單木中,林冠下層僅有1株針葉樹,所以并未對其指標進行計算。

    表3 基于ULS的估測樹高和樣地實測樹高的比較

    總體來看,應(yīng)用ULS數(shù)據(jù)提取的單木樹高和應(yīng)用超聲測高儀實測單木樹高之間的相關(guān)性很強,相關(guān)系數(shù)為0.97。應(yīng)用ULS數(shù)據(jù)提取的單木樹高存在低估現(xiàn)象,特別是闊葉樹被明顯低估了,rBias為-2.45%;針葉樹也存在略微的低估現(xiàn)象,rBias為-0.57%。從估測效果來看,針葉樹要比闊葉樹估測結(jié)果更接近實測樹高,針葉樹RMSE僅為0.41 m,闊葉樹為0.77 m;針葉樹rRMSE為3.44%,闊葉樹為6.23%。

    由表3可知,在不同的林冠層中,不同樹木類型樹高的提取精度有所不同。在林冠上層,針葉樹和闊葉樹的相關(guān)系數(shù)都為0.97,應(yīng)用ULS提取的樹高都被低估了,總體rBias為-1.78%,其中闊葉樹樹高比針葉樹樹高低估得更嚴重一些,針葉樹的rBias為-0.57%,闊葉樹的rBias為-2.73%,在樹高的估測誤差(rRMSE)上,針葉樹要比闊葉樹的估測誤差小2.68%。在林冠下層,估測的闊葉樹樹高存在很大誤差,其RMSE為0.80 m,rRMSE為10.14%,比林冠上層闊葉樹的rRMSE約高出4%。在林冠下層,基于ULS提取的樹高和樣地實測樹高的相關(guān)性比林冠上層的相關(guān)性略低,但仍具有很強的相關(guān)性(R=0.90)。

    應(yīng)用ULS數(shù)據(jù)估測的樹高與實測樹高比較得到的rRMSE隨冠層高度的變化見圖2。對于所有樹木來說,基于ULS提取的樹高誤差(rRMSE)在>9~11 m冠層范圍比≤9 m的范圍有了明顯的減小,rRMSE由10.46%減小到8.49%;在>11~13 m冠層范圍,估測樹高的rRMSE最小,僅為3.78%,這是由于在>11~13 m的冠層范圍,ULS點云和樣地實測均容易發(fā)現(xiàn)樹頂并測量樹高,因此具有較小的誤差;在13 m以上,隨著冠層高度升高,rRMSE略有升高,增加了約1.7%,上升的趨勢要遠小于林冠下層?;赨LS估測的闊葉樹誤差具有相同的變化規(guī)律,說明闊葉樹的誤差變化規(guī)律起主導(dǎo)作用。

    圖2 基于ULS的估測樹高和樣地實測樹高的rRMSE隨冠層高度的變化

    2.3.2 基于TLS的估測樹高與實測樹高比較

    去掉TLS探測樹高的離群值后,共有單木280株,其中林冠上、下層分別有單木227、53株,與樣地實測樹高比較結(jié)果見表4。TLS對林冠下層探測提取的單木樹高的rRMSE要遠遠小于林冠上層樹高,相同林冠層中,針葉樹的rRMSE也要小于闊葉樹。

    表4 基于TLS估測樹高和樣地實測樹高的比較

    總體來看,TLS探測提取的樹高和樣地實測樹高之間的相關(guān)性很強,相關(guān)系數(shù)為0.98,針葉樹和闊葉樹均存在低估現(xiàn)象,特別是闊葉樹,其rBias為-4.07%,被明顯低估了,針葉樹也存在一定的低估現(xiàn)象,rBias為-2.51%。從估測效果來看,針葉樹要比闊葉樹的樹高估測更準確,針葉樹RMSE為0.60 m,闊葉樹為0.95 m;針葉樹rRMSE為5.62%,闊葉樹為8.30%。

    由表4可知,在不同的林冠層中,不同樹木類型的樹高誤差也有所不同。在林冠上層,闊葉樹被低估的程度比針葉樹更明顯,闊葉樹的rBias為-4.70%,針葉樹的rBias為-3.23%。在估測誤差(rRMSE)上,闊葉樹比針葉樹誤差更大,rRMSE比針葉樹高2.69%;在林冠下層,TLS探測提取的針葉樹樹高誤差略小于闊葉樹,rRMSE比闊葉樹低1.16%,針葉樹和闊葉樹的rBias分別為2.00%和2.22%,存在高估現(xiàn)象。

    應(yīng)用TLS數(shù)據(jù)估測的樹高與實測樹高比較得到的rRMSE隨冠層高度的變化見圖3??傮w來看,隨著冠層高度的增加,基于TLS估測樹高誤差(rRMSE)也在不斷增加,在9 m冠層高度以上, rRMSE迅速增大,在>9~17 m的冠層高度上,rRMSE增加了約5.9%?;赥LS估測的闊葉樹和針葉樹樹高的誤差具有同樣的變化規(guī)律,闊葉樹樹高的rRMSE仍然主導(dǎo)著總體變化規(guī)律,而對于針葉樹來說,rRMSE變化更穩(wěn)定一些,rRMSE最小值出現(xiàn)在>7~9 m冠層范圍(3.23%),最大值出現(xiàn)在>13 m冠層(5.52%),變化幅度僅為2.29%,而闊葉樹樹高的誤差變化幅度高達6.71%。

    圖3 基于TLS的估測樹高和樣地實測樹高的rRMSE隨冠層高度的變化

    2.3.3 基于ULS和TLS估測樹高的比較

    基于ULS估測樹高和基于TLS估測樹高的比較如表5所示。由于在林冠下層,應(yīng)用ULS僅探測到1株針葉樹,因此并未對其指標進行計算。

    表5 基于ULS和TLS估測樹高的比較

    由表5可知,分別基于ULS和TLS估測的樹高具有很好的相關(guān)性(總體相關(guān)系數(shù)為0.95)??傮w而言,基于TLS估測的樹高要比基于ULS估測的樹高偏低,rBias為-3.56%,TLS和ULS探測提取的針葉樹和闊葉樹樹高相差不大,rRMSE僅相差0.18%。

    由表5可知,在林冠上層,基于TLS估測的樹高要比基于ULS估測的樹高偏低,rBias為-3.60%,而在冠下層,基于TLS估測的樹高要比基于ULS估測的樹高偏高,rBias為0.69%。這是由于TLS和ULS均受到冠層遮蔽的影響,在林冠上層,TLS不能對樹頂進行完整掃描,因此TLS估測的樹高要比ULS估測的樹高偏低;在林冠下層,ULS對樹木掃描不充分,因此ULS估測的樹高要比TLS估測的樹高偏低。

    分別應(yīng)用ULS和TLS數(shù)據(jù)估測樹高得到的rRMSE隨冠層高度的變化見圖4。由圖4可知,rRMSE變化范圍為5%~8%,在>9~13m冠層范圍的rRMSE最小(約為5%),隨著冠層高度的增加或減小,rRMSE也逐漸增加。闊葉樹的rRMSE具有相同的變化規(guī)律,說明闊葉樹的誤差變化起主導(dǎo)作用。

    圖4 基于ULS和TLS估測樹高的rRMSE隨高度的變化

    3 討 論

    3.1 基于CHD的林冠分層優(yōu)勢

    傳統(tǒng)的林層劃分方法需以實測數(shù)據(jù),并綜合考慮各林層每公頃蓄積量、相鄰林層間林木平均高、各林層平均胸徑以及主林層和其他林層的郁閉度人為地劃分[38],容易受外部環(huán)境及人為主觀因素的影響,分層依據(jù)的數(shù)據(jù)為外業(yè)實測獲得,費時費力,不利于大面積推廣?;赨LS獲取的點云數(shù)據(jù)受外部環(huán)境干擾小,能客觀地描述林分冠層結(jié)構(gòu)。由ULS點云計算得到的CHD可以用來描述冠層元素的垂直分布情況,體現(xiàn)了森林植被在垂直方向上的層次性和空間配置方式[34],本研究使用Lloyd算法計算CHD對林冠層進行分層也獲得了較好的研究結(jié)果,這與Maltamo等[31]的研究結(jié)果一致,因此,基于CHD對林冠層進行劃分能夠較好地反映不同數(shù)據(jù)源(ULS和TLS)估測樹高的適用范圍,而且計算簡便,可以應(yīng)用于大面積的估測,具有一定的推廣意義。

    3.2 離群值產(chǎn)生的原因分析

    ULS產(chǎn)生的離群值主要在林冠下層,是由于冠層的遮蔽對單木冠頂產(chǎn)生了誤判;而在林冠上層,對于部分樹冠相互交叉重疊的樹木,在人工解譯過程中難以正確判斷冠頂?shù)奈恢枚a(chǎn)生離群值。TLS產(chǎn)生的離群值大部分出現(xiàn)在林冠上層,與林冠層枝葉的密集程度密切相關(guān),枝葉越密集,在TLS點云上越容易產(chǎn)生樹頂?shù)恼`判。特別是與周圍樹冠的相互交叉給TLS的樹高提取帶來了困難,由此產(chǎn)生了離群值。實測樹高的誤差也有可能產(chǎn)生離群值,比如具有偏冠、歪斜的單木,有可能在野外調(diào)查中冠頂定位發(fā)生了偏移,從而產(chǎn)生了離群值。與LiDAR系統(tǒng)相比,野外樣地調(diào)查中受多種條件的影響,包括樹種、冠層密集程度和樹冠的形狀等,這些因素都會對單木樹高的精度產(chǎn)生影響[3]。本研究的結(jié)果與Wang等[19]的研究結(jié)果相似:ULS提取樹高的離群值主要是由樹冠部分被覆蓋的單木引起,TLS提取樹高的離群值出現(xiàn)在高林分密度(2 000株/hm2)的優(yōu)勢木中,由于TLS視角的原因,容易產(chǎn)生系統(tǒng)性低估。

    3.3 不同數(shù)據(jù)源對樹高估測精度的影響

    無論在林冠上層還是林冠下層,針葉樹和闊葉樹樹高的所有評價指標(RMSE、rRMSE、Bias、rBias和Pearson相關(guān)系數(shù))均大體相同,說明對于不同樹木種類(針葉樹和闊葉樹),基于ULS估測的樹高與TLS相比,誤差沒有明顯波動;隨著冠層高度的增加,rRMSE均逐漸增加。這主要是由于樹木越高,TLS越不容易掃描到樹頂,與ULS的一致性逐漸減??;而樹高越低,ULS越不容易掃描到樹頂,與TLS的一致性也減小。這種一致性的變化規(guī)律表明:當(dāng)TLS和ULS均能對單木進行準確的掃描時,均具有準確提取樹高的能力[19]。

    在不同的林冠層中,ULS和TLS提取的不同樹木類型的樹高精度有所不同。由表3和表4中對于ULS和TLS在林冠上層提取的樹高精度進行交叉比較發(fā)現(xiàn),在林冠上層,ULS比TLS估測的樹高更接近實測樹高,rRMSE相差2.56%,且ULS和TLS提取的針葉樹樹高要比闊葉樹樹高更接近實測樹高。由于冠層的遮擋,TLS并不能很好地對冠層上部特別是樹頂進行掃描,在本研究中的rBias約為-6.50%,存在低估現(xiàn)象,這與文獻[11-12]的研究結(jié)果一致。在>9~17 m的冠層高度上,rRMSE持續(xù)增大,增加了約5.9%,說明高度越高,TLS點云數(shù)字化能力越弱,提取的樹高誤差越大。這種現(xiàn)象在闊葉樹中的表現(xiàn)非常明顯,而林冠上層的針葉樹情況稍好,但是在點云中也存在樹頂數(shù)字化不完全的情況。本研究中,ULS估測的rRMSE在林冠上層增加了約1.7%,有逐漸增大的趨勢,推測有兩點原因:一方面,林冠上層的樹冠有交叉重疊的現(xiàn)象,影響了ULS對冠頂點的識別;另一方面,有可能是切線法測樹高的缺陷造成的,即超聲測高儀很難正確瞄準高大闊葉樹的頂點而造成估測結(jié)果偏高,且冠層越高趨勢越明顯[3]。在林冠下層,由于林冠層的遮擋,ULS對下層大部分樹木掃描不充分,探測到的數(shù)量比較少,而且難以確定冠頂點的位置,rRMSE達到了10.14%,與實測的樹高相比誤差較大。對表3和表4中ULS和TLS在林冠下層提取的樹高精度進行交叉比較發(fā)現(xiàn),ULS與TLS探測提取的樹高相比,rRMSE高6.31%。對于TLS來說,林冠下層無論針葉樹還是闊葉樹的rRMSE均穩(wěn)定在3%~5%,說明TLS在林冠下層能比較穩(wěn)定地估測樹高,且針葉樹樹高的rRMSE比闊葉樹減少了約1.2%。在林冠下層樹木枝葉密集,對于部分樹冠相互交叉的樹木,在ULS和TLS點云上很難正確發(fā)現(xiàn)樹頂點,這種現(xiàn)象在闊葉樹表現(xiàn)得更加明顯,因為闊葉樹的樹冠更容易相互交叉,而大部分針葉樹有規(guī)則的冠型,樹頂更容易被識別。

    參考文獻(reference):

    [1]ZHAO D H, KANE M, MARKEWITZ D, et al. Additive tree biomass equations for midrotation loblolly pine plantations[J]. For Sci, 2015, 61(4): 613-623. DOI:10.5849/forsci.14-193.

    [2]ZOU W T, ZENG W S, ZHANG L J, et al. Modeling crown biomass for four pine species in China[J]. Forests, 2015, 6(12): 433-449. DOI:10.3390/f6020433.

    [3]LARJAVAARA M, MULLER-LANDAU H C, METCALF J. Measuring tree height: a quantitative comparison of two common field methods in a moist tropical forest[J]. Methods Ecol Evol, 2013, 4(9): 793-801. DOI:10.1111/2041-210x.12071.

    [4]ZHEN Z, QUACKENBUSH L J, ZHANG L J. Trends in automatic individual tree crown detection and delineation: evolution of LiDAR data[J]. Remote Sens, 2016, 8(4): 333. DOI:10.3390/rs8040333.

    [5]ANDERSON K, GASTON K J. Lightweight unmanned aerial vehicles will revolutionize spatial ecology[J]. Front Ecol Environ, 2013, 11(3): 138-146. DOI: 10.1890/120150.

    [7]WALLACE L, LUCIEER A, WATSON C, et al. Development of a UAV-LiDAR system with application to forest inventory[J]. Remote Sens, 2012, 4(6): 1519-1543. DOI: 10.3390/rs4061519.

    [8]HOPKINSON C, CHASMER L, YOUNG-POW C, et al. Assessing forest metrics with a ground-based scanning lidar[J]. Can J For Res, 2004, 34(3): 573-583. DOI: 10.1139/x03-225.

    [9]PFEIFER N, WNTERHALDER D. Modelling of tree cross sections from terrestrial laser scanning data with free-form curves[J]. Pro of Isp Wor Las Sca for For and Lan Ass, 2004, 36(8): 76-81. DOI: 10.1109/TEST.2004.1387399.

    [10]熊妮娜, 王佳. 基于地基激光雷達的活立木材積提取算法[J]. 林業(yè)工程學(xué)報, 2020, 5(6): 143-148. XIONG N N, WANG J. Extratction algorithm for stand volume using ground-based laser scanner[J]. J For Eng, 2020, 5(6): 143-148. DOI:10.13360/j.issn.2096-1359.202001035.

    [11]BREDE B, LAU A, BARTHOLOMEUS H, et al. Comparing Riegl Ricopter UAV LiDAR derived canopy height and DBH with terrestrial LiDAR[J]. Sensors, 2017, 17(10): 2371. DOI: 10.3390/s17102371.

    [12]李丹, 龐勇, 岳彩榮, 等. 基于TLS數(shù)據(jù)的單木胸徑和樹高提取研究[J]. 北京林業(yè)大學(xué)學(xué)報, 2012, 34: 79-86. LI D, PANG Y, YUE C R, et al. Extraction of DBH and height of single tree based on TLS data[J]. Journal of Beijing Forestry University, 2012, 34: 79-86. DOI:10.13332/j.1000-1522.2012.04.027.

    [13]BEYENE S M, HUSSIN Y A, KLOOSTERMAN H E, et al. Fo-rest inventory and aboveground biomass estimation with terrestrial LiDAR in the tropical forest of Malaysia[J]. Can J Remote Sens, 2020, 46(2): 130-145. DOI:10.1080/07038992.2020.1759036.

    [14]GOODWIN N R, COOPS N C, CULVENOR D S. Assessment of forest structure with airborne LiDAR and the effects of platform altitude[J]. Remote Sens Environ, 2006, 103(2): 140-152. DOI:10.1016/j.rse.2006.03.003.

    [15]ANDERSEN H E, REUTEBUCH S E, MCGAUGHEY R J. A rigorous assessment of tree height measurements obtained using airborne lidar and conventional field methods[J]. Can J Remote Sens, 2006, 32(5): 355-366. DOI: 10.5589/m06-030.

    [16]SIBONA E, VITALI A, MELONI F, et al. Direct measurement of tree height provides different results on the assessment of LiDAR accuracy[J]. Forests, 2016, 8(1): 7. DOI: 10.3390/f8010007.

    [18]WALLACE L, MUSK R, LUCIEER A. An assessment of the repeatability of automatic forest inventory metrics derived from UAV-borne laser scanning data[J]. IEEE Trans Geosci Remote Sens, 2014, 52(11): 7160-7169. DOI:10.1109/TGRS.2014.2308208.

    [20]VAGLIO LAURIN G, DING J Q, DISNEY M, et al. Tree height in tropical forest as measured by different ground, proximal, and remote sensing instruments, and impacts on above ground biomass estimates[J]. Int J Appl Earth Obs Geoinformation, 2019, 82: 101899. DOI:10.1016/j.jag.2019.101899.

    [21]TANABE S I, TODA M J, VINOKUROVA A V. Tree shape, fo-rest structure and diversity of drosophilid community: comparison between boreal and temperate birch forests[J]. Ecol Res, 2001, 16(3): 369-385. DOI: 10.1046/j.1440-1703.2001.00402.x.

    [22]MIURA N, JONES S D. Characterizing forest ecological structure using pulse types and heights of airborne laser scanning[J]. Remote Sens Environ, 2010, 114(5): 1069-1076. DOI: 10.1016/j.rse.2009.12.017.

    [23]郎春博, 賈鶴鳴, 邢致愷, 等 . 基于改進粒子群算法的植物冠層圖像分割[J]. 森林工程,2019,35(1):47-52. LANG C B, JIA H M, XING Z K,et al. Multi threshold segmentation of plant canopy image based on improved particle swarm optimization[J]. Forest Engineering,2019,35(1):47-52.

    [24]胡文杰, 崔鴻俠, 王曉榮, 等. 三峽庫區(qū)馬尾松次生林林分結(jié)構(gòu)特征分析[J]. 南京林業(yè)大學(xué)學(xué)報(自然科學(xué)版), 2019, 43(3): 67-76. HU W J, CUI H X, WANG X R, et al. Structure characteristics ofPinusmassonianasecondary forest in the Three Reservoir Area[J]. J Nanjing For Univ (Nat Sci Ed), 2019, 43(3): 67-76. DOI:10.3969/j.issn.1000-2006.201805075.

    [25]胡傳偉, 孫冰, 莊梅梅, 等. 深圳羊臺山近自然風(fēng)景林樹種組成與垂直結(jié)構(gòu)[J]. 南京林業(yè)大學(xué)學(xué)報(自然科學(xué)版), 2010, 34(4): 112-116. HU C W, SUN B, ZHUANG M M, et al. Study on species composition and vertical structure of near-nature scenic forest in Mt.Yangtai, Shenzhen[J]. J Nanjing For Univ (Nat Sci Ed), 2010, 34(4): 112-116. DOI: 10.3969/j.issn.1000-2006.2010.04.025.

    [26]ZIMBLE D A, EVANS D L, CARLSON G C, et al. Characterizing vertical forest structure using small-footprint airborne LiDAR[J]. Remote Sens Environ, 2003, 87(2/3): 171-182. DOI: 10.1016/S0034-4257(03)00139-1.

    [27]趙靜, 李靜, 柳欽火. 森林垂直結(jié)構(gòu)參數(shù)遙感反演綜述[J]. 遙感學(xué)報, 2013, 17(4): 697-716. ZHAO J, LI J, LIU Q H. Review of forest vertical structure parameter inversion based on remote sensing technology[J]. J Remote Sens, 2013, 17(4): 697-716. DOI:10.11834/jrs.20132183.

    [28]ZHAO K G, POPESCU S, NELSON R. Lidar remote sensing of forest biomass: a scale-invariant estimation approach using airborne lasers[J]. Remote Sens Environ, 2009, 113(1): 182-196. DOI:10.1016/j.rse.2008.09.009.

    [29]QIN H M, WANG C, XI X H, et al. Simulating the effects of the airborne lidar scanning angle, flying altitude, and pulse density for forest foliage profile retrieval[J]. Appl Sci, 2017, 7(7): 712. DOI:10.3390/app7070712.

    [30]LEFSKY M A, COHEN W B, ACKER S A, et al. Lidar remote sensing of the canopy structure and biophysical properties of Douglas-fir western hemlock forests[J]. Remote Sens Environ, 1999, 70(3): 339-361. DOI: 10.1016/S0034-4257(99)00052-8.

    [31]MALTAMO M, PACKALéN P, YU X, et al. Identifying and quantifying structural characteristics of heterogeneous boreal forests using laser scanner data[J]. For Ecol Manag, 2005, 216(1/2/3): 41-50. DOI:10.1016/j.foreco.2005.05.034.

    [32]ZHAO X Q, GUO Q H, SU Y J, et al. Improved progressive TIN densification filtering algorithm for airborne LiDAR data in forested areas[J]. ISPRS J Photogramm Remote Sens, 2016, 117: 79-91. DOI:10.1016/j.isprsjprs.2016.03.016.

    [33]SCHNABEL R, KLEIN R. Octree-based point-cloud compression[C]//Eur Sym Point-Based Gra, 2006: 111-120. DOI: 10.2312/SPBG/SPBG06/111-120.

    [34]劉浩, 張崢男, 曹林. 機載激光雷達森林垂直結(jié)構(gòu)剖面參數(shù)的沿海平原人工林林分特征反演[J]. 遙感學(xué)報, 2018, 22(5): 872-888. LIU H, ZHANG Z N, CAO L. Estimating forest stand characteristics in a coastal plain forest plantation based on vertical structure profile parameters derived from ALS data[J]. J Remote Sens, 2018, 22(5): 872-888. DOI:10.11834/jrs.20187465.

    [35]LLOYD S. Least squares quantization in PCM[J]. IEEE Trans Inf Theory, 1982, 28(2): 129-137. DOI:10.1109/TIT.1982.1056489.

    [36]BAZEZEW M N, HUSSIN Y A, KLOOSTERMAN E H. Integrating airborne LiDAR and terrestrial laser scanner forest parameters for accurate above-ground biomass/carbon estimation in Ayer Hitam tropical forest, Malaysia[J]. Int J Appl Earth Obs Geoinformation, 2018, 73: 638-652. DOI:10.1016/j.jag.2018.07.026.

    [37]LEFSKY M A, COHEN W B, ACKER S A, et al. Lidar remote sensing of forest canopy structure and related biophysical parameters at H. J. Andrews Experimental Forest, Oregon, USA[C]//IGARSS ’98. Sensing and Managing the Environment. 1998 IEEE International Geoscience and Remote Sensing. Symposium Proceedings. (Cat. No. 98CH36174). July 6-10, 1998, Seattle, WA, USA. IEEE, 1998: 1252-1254. DOI:10.1109/IGARSS.1998.691367.

    [38]李鳳日. 測樹學(xué)[M]. 4版. 北京: 中國林業(yè)出版社, 2019:64. LI F R.Forest mensuration[M].4th ed. Beijing:China Forestry Publishing House,2019: 64.

    猜你喜歡
    針葉樹單木闊葉樹
    地基與無人機激光雷達結(jié)合提取單木參數(shù)
    杉木與不同闊葉樹種混交生長效應(yīng)分析
    防護林科技(2022年2期)2022-04-25 23:47:32
    融合LiDAR點云與高分影像的單木檢測方法研究
    無人機影像匹配點云單木識別算法
    遙感信息(2019年1期)2019-03-22 01:38:16
    基于雙尺度體元覆蓋密度的TLS點云數(shù)據(jù)單木識別算法
    森林工程(2018年5期)2018-05-14 13:54:30
    針葉樹容器苗木培育技術(shù)
    湖南速生、中生、慢生闊葉樹組單株生長模型構(gòu)建
    針葉樹容器育苗造林技術(shù)研究
    福建省闊葉樹二元材積方程修訂
    閩北山地不同肥料對四個闊葉樹種生長的影響
    一夜夜www| 国产精品乱码一区二三区的特点 | 亚洲人成电影观看| 男人舔女人的私密视频| 50天的宝宝边吃奶边哭怎么回事| 午夜精品在线福利| 99riav亚洲国产免费| 亚洲一区二区三区色噜噜| 18禁黄网站禁片午夜丰满| 黄色视频不卡| 精品人妻1区二区| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品美女久久av网站| 亚洲成人免费电影在线观看| 正在播放国产对白刺激| 真人一进一出gif抽搐免费| av有码第一页| 涩涩av久久男人的天堂| 欧美 亚洲 国产 日韩一| 桃色一区二区三区在线观看| 日本a在线网址| 国产不卡一卡二| 欧美性长视频在线观看| 精品久久久久久成人av| 黄片大片在线免费观看| 日韩三级视频一区二区三区| 欧美日本亚洲视频在线播放| 日韩精品中文字幕看吧| 中文字幕最新亚洲高清| 别揉我奶头~嗯~啊~动态视频| 嫩草影院精品99| 国产精品香港三级国产av潘金莲| 日本黄色视频三级网站网址| 美国免费a级毛片| 九色亚洲精品在线播放| 乱人伦中国视频| 色综合婷婷激情| 日韩欧美在线二视频| 真人一进一出gif抽搐免费| 国产精品av久久久久免费| av欧美777| 一区二区三区高清视频在线| 香蕉久久夜色| 亚洲欧美精品综合久久99| 精品免费久久久久久久清纯| 一区二区三区高清视频在线| 制服人妻中文乱码| 成人亚洲精品一区在线观看| 午夜a级毛片| 亚洲国产欧美日韩在线播放| 韩国av一区二区三区四区| 亚洲国产精品久久男人天堂| 国产极品粉嫩免费观看在线| 亚洲七黄色美女视频| 在线免费观看的www视频| 亚洲伊人色综图| 侵犯人妻中文字幕一二三四区| 亚洲色图综合在线观看| 伊人久久大香线蕉亚洲五| 少妇 在线观看| 亚洲黑人精品在线| 国产亚洲精品第一综合不卡| 黄色成人免费大全| 成人精品一区二区免费| 日本三级黄在线观看| 欧美日韩福利视频一区二区| 男女下面插进去视频免费观看| 免费在线观看完整版高清| 亚洲精品一区av在线观看| 一级黄色大片毛片| 成人亚洲精品av一区二区| 婷婷丁香在线五月| 一夜夜www| АⅤ资源中文在线天堂| 精品国产一区二区三区四区第35| 国产一级毛片七仙女欲春2 | 午夜久久久久精精品| 欧美人与性动交α欧美精品济南到| 日本免费a在线| 国产成年人精品一区二区| 亚洲国产精品999在线| 一区二区三区激情视频| 久久精品亚洲精品国产色婷小说| 精品久久久久久久久久免费视频| 女人被狂操c到高潮| 99在线视频只有这里精品首页| 亚洲aⅴ乱码一区二区在线播放 | 91九色精品人成在线观看| 国产亚洲av高清不卡| 熟妇人妻久久中文字幕3abv| 91老司机精品| www国产在线视频色| 高清毛片免费观看视频网站| 午夜福利免费观看在线| 久久精品国产亚洲av香蕉五月| 此物有八面人人有两片| 国产亚洲av高清不卡| 亚洲va日本ⅴa欧美va伊人久久| 亚洲 欧美一区二区三区| svipshipincom国产片| 欧美丝袜亚洲另类 | 精品国产超薄肉色丝袜足j| 91字幕亚洲| 欧美黑人欧美精品刺激| 一个人观看的视频www高清免费观看 | 欧美黄色淫秽网站| 国产精品久久电影中文字幕| 亚洲精品av麻豆狂野| 精品国产国语对白av| 国产亚洲精品综合一区在线观看 | 村上凉子中文字幕在线| 国产高清视频在线播放一区| 欧美精品亚洲一区二区| 午夜福利影视在线免费观看| 高潮久久久久久久久久久不卡| 国产成人啪精品午夜网站| 操美女的视频在线观看| 99香蕉大伊视频| 久久久久国内视频| 国产高清有码在线观看视频 | bbb黄色大片| 久久性视频一级片| 69av精品久久久久久| 国产成人精品在线电影| 亚洲成av片中文字幕在线观看| av视频免费观看在线观看| 国产真人三级小视频在线观看| 韩国精品一区二区三区| 不卡一级毛片| 黄色成人免费大全| 亚洲国产精品合色在线| 午夜福利影视在线免费观看| 欧美乱色亚洲激情| 黄网站色视频无遮挡免费观看| 岛国在线观看网站| 亚洲国产欧美网| 多毛熟女@视频| 99久久综合精品五月天人人| 欧美黄色淫秽网站| 两性夫妻黄色片| 成人亚洲精品一区在线观看| 国产一卡二卡三卡精品| 久久中文看片网| 成人亚洲精品一区在线观看| 亚洲男人天堂网一区| 国产91精品成人一区二区三区| av在线播放免费不卡| 村上凉子中文字幕在线| av天堂在线播放| 成人国产一区最新在线观看| 悠悠久久av| 黄网站色视频无遮挡免费观看| 一级a爱视频在线免费观看| 国产乱人伦免费视频| 曰老女人黄片| 欧美 亚洲 国产 日韩一| 亚洲专区字幕在线| 精品国产国语对白av| 亚洲中文日韩欧美视频| 中文字幕久久专区| 看免费av毛片| 欧美最黄视频在线播放免费| 日韩大尺度精品在线看网址 | 搡老岳熟女国产| 无遮挡黄片免费观看| 色在线成人网| 亚洲欧美日韩无卡精品| 亚洲av电影在线进入| 国产一区二区三区综合在线观看| www.999成人在线观看| 久久天堂一区二区三区四区| 国产精品久久久人人做人人爽| 一a级毛片在线观看| 精品一品国产午夜福利视频| 成人三级黄色视频| 99re在线观看精品视频| 男人舔女人下体高潮全视频| 男人舔女人下体高潮全视频| 午夜a级毛片| 国产蜜桃级精品一区二区三区| 国产主播在线观看一区二区| 成人三级黄色视频| 国产不卡一卡二| 久久精品国产亚洲av香蕉五月| 少妇被粗大的猛进出69影院| 婷婷六月久久综合丁香| 精品国产美女av久久久久小说| 久久天堂一区二区三区四区| 黄片小视频在线播放| 亚洲男人天堂网一区| 国产一区二区三区综合在线观看| 女生性感内裤真人,穿戴方法视频| 久久久久久久午夜电影| 久久精品影院6| 亚洲中文av在线| 亚洲aⅴ乱码一区二区在线播放 | 999精品在线视频| 亚洲成人免费电影在线观看| 日本 av在线| 国产蜜桃级精品一区二区三区| 国产一卡二卡三卡精品| 99国产精品免费福利视频| 亚洲精品美女久久久久99蜜臀| 19禁男女啪啪无遮挡网站| 91成年电影在线观看| 丝袜美腿诱惑在线| 亚洲九九香蕉| 91av网站免费观看| 精品久久久久久久毛片微露脸| 制服丝袜大香蕉在线| 最近最新中文字幕大全免费视频| 久久久久久久久久久久大奶| 一本久久中文字幕| 一区在线观看完整版| 日本一区二区免费在线视频| 最新在线观看一区二区三区| 91大片在线观看| 中文字幕久久专区| 波多野结衣巨乳人妻| 久久久久久久久中文| 国产成+人综合+亚洲专区| 国产精品国产高清国产av| 一二三四在线观看免费中文在| 亚洲一区二区三区色噜噜| 黄色 视频免费看| 好男人电影高清在线观看| 国产男靠女视频免费网站| 啦啦啦免费观看视频1| 国产欧美日韩一区二区三| 国产高清激情床上av| 亚洲中文字幕一区二区三区有码在线看 | 99精品在免费线老司机午夜| 老鸭窝网址在线观看| 十八禁网站免费在线| 国产av在哪里看| 亚洲一区中文字幕在线| 国产人伦9x9x在线观看| 亚洲色图 男人天堂 中文字幕| 女人高潮潮喷娇喘18禁视频| 悠悠久久av| 亚洲av片天天在线观看| 国产精品野战在线观看| 又黄又粗又硬又大视频| 老熟妇乱子伦视频在线观看| 夜夜看夜夜爽夜夜摸| 成人av一区二区三区在线看| 亚洲av成人不卡在线观看播放网| 人人妻人人澡欧美一区二区 | 亚洲欧美激情综合另类| 欧美日韩亚洲国产一区二区在线观看| 成人18禁高潮啪啪吃奶动态图| 窝窝影院91人妻| 亚洲va日本ⅴa欧美va伊人久久| 久久国产精品影院| 久久久久久免费高清国产稀缺| 日韩免费av在线播放| 国产成人系列免费观看| 国产免费av片在线观看野外av| 国产精品免费一区二区三区在线| 琪琪午夜伦伦电影理论片6080| 狠狠狠狠99中文字幕| 给我免费播放毛片高清在线观看| 人人妻人人澡人人看| 欧美另类亚洲清纯唯美| 九色国产91popny在线| 亚洲av成人一区二区三| 美女国产高潮福利片在线看| 亚洲精品在线美女| 日韩精品青青久久久久久| 日本 av在线| av有码第一页| 人成视频在线观看免费观看| 亚洲五月婷婷丁香| 国产精品秋霞免费鲁丝片| 国产成人啪精品午夜网站| 亚洲 欧美一区二区三区| 757午夜福利合集在线观看| 亚洲精品国产区一区二| 嫩草影视91久久| 久久久国产精品麻豆| 一本久久中文字幕| 日日摸夜夜添夜夜添小说| 身体一侧抽搐| 91字幕亚洲| 少妇被粗大的猛进出69影院| 女性生殖器流出的白浆| 亚洲欧美一区二区三区黑人| 成人手机av| 亚洲第一欧美日韩一区二区三区| 我的亚洲天堂| 老司机深夜福利视频在线观看| 激情在线观看视频在线高清| 日韩免费av在线播放| 天堂√8在线中文| 91av网站免费观看| 两性夫妻黄色片| 欧美一级毛片孕妇| 亚洲欧美日韩高清在线视频| 欧美国产精品va在线观看不卡| 亚洲美女黄片视频| 亚洲精品美女久久久久99蜜臀| 亚洲欧美一区二区三区黑人| tocl精华| 熟女少妇亚洲综合色aaa.| av福利片在线| 色尼玛亚洲综合影院| 国产片内射在线| av电影中文网址| 欧美日韩黄片免| 精品国产一区二区三区四区第35| 一二三四社区在线视频社区8| 丝袜美腿诱惑在线| 男人操女人黄网站| 精品久久久久久成人av| 国产成人精品久久二区二区免费| 国产亚洲av高清不卡| 午夜影院日韩av| 亚洲久久久国产精品| 啪啪无遮挡十八禁网站| 很黄的视频免费| 精品国产乱子伦一区二区三区| 色婷婷久久久亚洲欧美| 久久婷婷成人综合色麻豆| 亚洲 国产 在线| 日韩欧美三级三区| 国产成人精品无人区| 精品一区二区三区四区五区乱码| 日日爽夜夜爽网站| 国产亚洲精品第一综合不卡| 国内毛片毛片毛片毛片毛片| 国产亚洲欧美精品永久| 一区二区三区激情视频| АⅤ资源中文在线天堂| 97超级碰碰碰精品色视频在线观看| 亚洲专区国产一区二区| 亚洲九九香蕉| 国产aⅴ精品一区二区三区波| 日本五十路高清| 久久久久久大精品| 十八禁人妻一区二区| 精品日产1卡2卡| 一区二区三区国产精品乱码| 妹子高潮喷水视频| 一区福利在线观看| 午夜亚洲福利在线播放| 很黄的视频免费| 欧美 亚洲 国产 日韩一| 美女大奶头视频| 色综合站精品国产| 青草久久国产| 国产人伦9x9x在线观看| 一边摸一边抽搐一进一出视频| 午夜福利18| 国产一区二区三区在线臀色熟女| 国产熟女午夜一区二区三区| 国产aⅴ精品一区二区三区波| 老司机在亚洲福利影院| 亚洲欧美日韩无卡精品| 国产精品爽爽va在线观看网站 | 一区二区三区高清视频在线| 91麻豆精品激情在线观看国产| 搡老妇女老女人老熟妇| av片东京热男人的天堂| 日日爽夜夜爽网站| 老汉色∧v一级毛片| 18禁国产床啪视频网站| 母亲3免费完整高清在线观看| aaaaa片日本免费| 日日爽夜夜爽网站| 老司机靠b影院| 91精品三级在线观看| 久久久久九九精品影院| 国产高清videossex| 欧美精品啪啪一区二区三区| 精品久久久久久,| 大陆偷拍与自拍| 人人妻人人澡人人看| 日日摸夜夜添夜夜添小说| 亚洲成人久久性| 国产欧美日韩精品亚洲av| 青草久久国产| 精品一区二区三区av网在线观看| 岛国在线观看网站| 国产又色又爽无遮挡免费看| 国产高清有码在线观看视频 | 国产亚洲精品av在线| 999精品在线视频| 久久中文看片网| 大香蕉久久成人网| 级片在线观看| 亚洲三区欧美一区| 久久精品国产清高在天天线| 麻豆av在线久日| 国产人伦9x9x在线观看| 国产亚洲av高清不卡| 国产一卡二卡三卡精品| 男人的好看免费观看在线视频 | 国产精品香港三级国产av潘金莲| 黄频高清免费视频| 亚洲欧美精品综合久久99| 国产亚洲精品av在线| av天堂在线播放| 欧美黑人精品巨大| 亚洲人成电影观看| 99国产精品99久久久久| 国产精品影院久久| 男女之事视频高清在线观看| 国产成人av激情在线播放| 亚洲美女黄片视频| www.www免费av| 国产精品一区二区在线不卡| 一级a爱片免费观看的视频| 三级毛片av免费| 欧美日韩瑟瑟在线播放| 成人三级黄色视频| cao死你这个sao货| 精品久久久久久成人av| 麻豆成人av在线观看| 亚洲男人天堂网一区| 99香蕉大伊视频| 欧美激情久久久久久爽电影 | 精品国产乱子伦一区二区三区| 欧美一区二区精品小视频在线| 国产一区二区激情短视频| 精品第一国产精品| 日韩免费av在线播放| 99在线人妻在线中文字幕| 波多野结衣一区麻豆| 免费一级毛片在线播放高清视频 | 一夜夜www| 亚洲avbb在线观看| 中文亚洲av片在线观看爽| 亚洲精品粉嫩美女一区| 操美女的视频在线观看| 色老头精品视频在线观看| 91麻豆精品激情在线观看国产| 好男人在线观看高清免费视频 | 九色国产91popny在线| 国产欧美日韩一区二区三区在线| 亚洲性夜色夜夜综合| 欧美日韩亚洲国产一区二区在线观看| 又紧又爽又黄一区二区| 亚洲av成人av| 热99re8久久精品国产| 亚洲va日本ⅴa欧美va伊人久久| 婷婷六月久久综合丁香| 真人一进一出gif抽搐免费| 亚洲av成人一区二区三| 日本欧美视频一区| 国产成人影院久久av| 老司机在亚洲福利影院| 国产黄a三级三级三级人| 欧美久久黑人一区二区| 久久精品国产清高在天天线| 国产亚洲精品第一综合不卡| 十八禁网站免费在线| 精品一区二区三区四区五区乱码| 亚洲午夜理论影院| 美女大奶头视频| 久久久久国内视频| 国产精品电影一区二区三区| 国产成人欧美| 两个人免费观看高清视频| 亚洲人成电影观看| 国产激情欧美一区二区| 老司机靠b影院| 久久婷婷人人爽人人干人人爱 | 国产高清videossex| 色婷婷久久久亚洲欧美| 亚洲视频免费观看视频| 亚洲人成电影免费在线| 亚洲专区字幕在线| 国产精品久久久人人做人人爽| 制服诱惑二区| 午夜日韩欧美国产| 日韩精品免费视频一区二区三区| 男人的好看免费观看在线视频 | 少妇被粗大的猛进出69影院| av福利片在线| 天天一区二区日本电影三级 | 欧美成人免费av一区二区三区| 一边摸一边做爽爽视频免费| 国产亚洲欧美在线一区二区| 中文字幕色久视频| 夜夜看夜夜爽夜夜摸| 日韩视频一区二区在线观看| 国产精品久久电影中文字幕| 999精品在线视频| 国产伦人伦偷精品视频| 久久人妻av系列| 亚洲伊人色综图| 中亚洲国语对白在线视频| 窝窝影院91人妻| av在线播放免费不卡| 成人手机av| 久久精品亚洲精品国产色婷小说| 18禁国产床啪视频网站| 97超级碰碰碰精品色视频在线观看| 国产精品1区2区在线观看.| 亚洲中文av在线| 午夜福利成人在线免费观看| 悠悠久久av| 99国产精品一区二区蜜桃av| 国产精品电影一区二区三区| 午夜福利欧美成人| 99香蕉大伊视频| 免费无遮挡裸体视频| 久久香蕉国产精品| 国产成人啪精品午夜网站| 18禁观看日本| 1024视频免费在线观看| 99久久99久久久精品蜜桃| 午夜免费观看网址| 欧美丝袜亚洲另类 | 亚洲av电影在线进入| 国产片内射在线| 涩涩av久久男人的天堂| 国产亚洲av嫩草精品影院| 午夜日韩欧美国产| 久久精品国产清高在天天线| 国产97色在线日韩免费| 亚洲精品中文字幕在线视频| 亚洲av第一区精品v没综合| 久久久久久久久中文| 久久中文看片网| 久久久久精品国产欧美久久久| 他把我摸到了高潮在线观看| 高清在线国产一区| 久久狼人影院| 亚洲人成电影免费在线| 少妇 在线观看| 极品教师在线免费播放| 人妻丰满熟妇av一区二区三区| 国产又爽黄色视频| 国产国语露脸激情在线看| 久久国产精品影院| 性色av乱码一区二区三区2| 亚洲精品av麻豆狂野| 精品国产美女av久久久久小说| 亚洲一区高清亚洲精品| 不卡一级毛片| 国产精品久久久久久人妻精品电影| 麻豆一二三区av精品| 亚洲第一青青草原| 老司机午夜福利在线观看视频| 日本免费一区二区三区高清不卡 | 亚洲自偷自拍图片 自拍| 国产午夜精品久久久久久| 制服人妻中文乱码| 日本欧美视频一区| aaaaa片日本免费| 亚洲国产欧美日韩在线播放| 色播亚洲综合网| 国产成+人综合+亚洲专区| 国产一区二区在线av高清观看| 欧美国产日韩亚洲一区| 午夜福利18| 精品一区二区三区四区五区乱码| 亚洲一区二区三区色噜噜| 久久精品91无色码中文字幕| 两性午夜刺激爽爽歪歪视频在线观看 | 久久香蕉国产精品| 久久人人爽av亚洲精品天堂| 男女做爰动态图高潮gif福利片 | 国产片内射在线| 少妇粗大呻吟视频| 国产精品日韩av在线免费观看 | 日韩精品中文字幕看吧| 亚洲成人国产一区在线观看| 欧美人与性动交α欧美精品济南到| 午夜免费鲁丝| 男女午夜视频在线观看| 久久亚洲真实| 99精品在免费线老司机午夜| 久久久久亚洲av毛片大全| 69av精品久久久久久| 黄片播放在线免费| 亚洲精品一区av在线观看| 亚洲黑人精品在线| 91老司机精品| 91成人精品电影| 日韩大码丰满熟妇| 亚洲午夜理论影院| 黄色女人牲交| 国产麻豆成人av免费视频| 亚洲av五月六月丁香网| 国产欧美日韩精品亚洲av| 中文字幕色久视频| 亚洲精品在线美女| 久久久久国内视频| 黑人巨大精品欧美一区二区mp4| 久久久久久久久中文| 国产一区二区三区综合在线观看| 精品久久久久久久久久免费视频| 久久中文看片网| 91大片在线观看| 1024香蕉在线观看| 国产aⅴ精品一区二区三区波| 国产精品久久久人人做人人爽| 久久伊人香网站| 国产高清有码在线观看视频 | 91麻豆av在线| 亚洲三区欧美一区| 国产成人一区二区三区免费视频网站| 亚洲欧美日韩无卡精品| 大码成人一级视频| 91麻豆精品激情在线观看国产| 窝窝影院91人妻| 女性生殖器流出的白浆| 亚洲国产欧美网| 久久久久久久午夜电影| 美女国产高潮福利片在线看| 香蕉国产在线看| 日韩欧美国产一区二区入口| 一卡2卡三卡四卡精品乱码亚洲|