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

    基于地基激光雷達的欒樹分形特征分析

    2022-07-26 06:21:16李輝林沂孟祥爽史振偉蔡萬園
    關鍵詞:欒樹分枝分形

    李輝,林沂,孟祥爽,史振偉,蔡萬園

    基于地基激光雷達的欒樹分形特征分析

    李輝,林沂*,孟祥爽,史振偉,蔡萬園

    北京大學地球與空間科學學院遙感與地理信息系統(tǒng)研究所空間信息集成與3S工程應用北京市重點實驗室, 北京 100871

    分枝結構的幾何特征是樹木分枝按照一定規(guī)律組合而成的形態(tài),體現(xiàn)了個體樹木在生理、生態(tài)角度對于環(huán)境的適應能力,以及對光、水、養(yǎng)分、空間等生長影響因子的競爭能力,是生態(tài)領域研究的一個重要內(nèi)容。地基激光雷達(Terrestrial Laser Scanning, TLS)技術的發(fā)展為獲取單木結構的參數(shù)提供了一種便捷、無損的綠色方式。本文使用基于TLS點云的方法來代替耗費大量人力的方式來獲取樹木分枝的結構參數(shù),并進行分形特征的分析,發(fā)現(xiàn)在欒樹的生長發(fā)育階段,為了充分利用光照條件和生長空間,隨著DBH的增加,分枝率和枝徑比同步上升,但分枝角度呈先上升后下降的趨勢。

    地基激光雷達; 欒樹; 特征分析

    樹木可以被看作各種組件組合而成的集合體[1],不同分枝特征及樹葉、芽的動態(tài)變化反映了樹木對光、水、養(yǎng)分與空間等資源的利用[2],表達了樹木在不同生長階段為了適應環(huán)境的最終選擇的應對策略[3]。對于樹木分形特征的研究,是為了了解樹木對于環(huán)境空間的利用能力,探討生長環(huán)境與個體生長之間的相互關系[4,5]。以往的研究表明,樹木的分形特征和冠層構型反映了其對生存環(huán)境的妥協(xié)策略[6,7]。而經(jīng)過幾十年的發(fā)展,人們對不同樹木構型和分枝特征的認識和研究已逐步深入[8]。孫書存等通過實測數(shù)據(jù)分析了遼東栓皮櫟的分形特征[9]。張文輝等通過分析沙參的分枝特征,發(fā)現(xiàn)該理論可用于研究瀕危種群的分布[10]。馬克明等利用分形理論研究了東北草原格局,并提出了計算分支格局的模型[11]。趙相健等研究發(fā)現(xiàn),競爭地位不同的太白紅杉的分枝格局差異性與地形具備較強的關聯(lián)性[12]。林勇等通過分析不同離海距離下的木麻黃分支格局,得出持續(xù)性風力對其影響較大,冠層形狀和分枝長度都體現(xiàn)了對環(huán)境的響應特征[13]。但由于測量方式的限制,該領域的研究也多基于草本植物,雖在之后擴大到高大喬木,但依然存在樣本數(shù)據(jù)較少,難以大規(guī)模進行研究的問題。

    地基激光雷達系統(tǒng)又稱地面激光掃描(Terrestrial Laser Scanning, TLS)系統(tǒng)[14,15],可通過發(fā)射激光脈沖獲取目標表面的三維坐標[16]。地基激光雷達掃描并結合相應算法是獲取樹木的三維結構的便捷方式,為實現(xiàn)可重復,無創(chuàng)且客觀的獲取樹木分枝結構的提供了新的方式[15,17]。利用TLS掃描能夠便捷的估算樹木的結構參數(shù)[18,19],使用一些算法和定量結構模型,樹木胸徑[20,21]、樹高[22-24]、樹冠直徑和面積[25,26]可以從單木點云中估算出來。Dassot M等獲取了42個不同樹種、不同大小的樣品的地基激光雷達數(shù)據(jù),從點云數(shù)據(jù)中選取基徑大于7 cm的樹枝,通過柱面擬合方法進行基徑預測建模[27]。R?der A等對13 m以下的31棵針葉樹每輪的枝高和枝數(shù)進行了自動識別和估計[28]。Raumonen P等提出了一種模型可以近似地得到整棵樹或部分樹的分枝結構和大小特性[29]。Lau A等首次利用定量結構模型提取樹木分枝參數(shù)[30],Chen S等利用QSM(Quantitative Structure Model,QSM)模型進行生物量估算[31]。倪文健等基于樹干點云垂直連續(xù)分布的特征提出了一種樹干識別方法[32]。王寧寧等由點云數(shù)據(jù)得到楊樹的樹冠特征參數(shù)并進行分析[33]。劉金鵬等基于Skeltre算法提出一種利用柱形擬合法修正分支點處分支延伸方向的新方法[34]。截止目前,利用TLS點云對分枝分形特征的研究還比較少,本文旨在使用TLS點云的來代替耗費人力的方式來獲取樹木分枝的結構參數(shù),并進行分形特征的研究。

    欒樹俗稱“燈籠樹”,屬于落葉喬木,耐寒、耐旱[35]。在我國多作為園林綠化樹木,尤其在北方地區(qū)分布較廣[36],同時其根系發(fā)達,抗污染性強,具備修復工業(yè)污染土地的能力[37],葉子可作為染料的原材料[38],欒樹花不僅可以藥用還能做油料[39],具有較高的生態(tài)價值和經(jīng)濟價值。因此對于欒樹的相關研究不僅具備經(jīng)濟意義,在目前追求低碳環(huán)保、建立綠色中國的大背景下,還具備較強的生態(tài)學意義。本文以北京市昌平區(qū)的40棵10年生綠化欒樹為研究對象,同時獲取點云數(shù)據(jù)和人工測量數(shù)據(jù),通過QSM定量結構模型提取每棵樹木的分枝長度、直徑等結構參數(shù),驗證QSM算法可靠性的同時,針對欒樹研究其分枝格局和特征。

    圖 1 技術路線圖

    1 材料與方法

    1.1 人工測量數(shù)據(jù)

    研究區(qū)域位于北京市昌平區(qū),將研究對象40棵欒樹按順序依次編號1至40號,測量的參數(shù)主要有枝下木的長度、直徑,1級和2級分枝的長度、直徑,對于每一級分枝按照逆時針進行排序,子分枝與母分枝對應記錄。為了便于統(tǒng)計,將枝下木定義為0級分枝,由0級分枝生長出來的為1級分枝,源自1級分枝的定義為2級分枝,同樣,源自于2級分枝和3級分枝的定義為3級分枝和4級分枝。分枝長度定義為基部到分枝末端之間的距離。直徑定義為基部直徑和末端直徑的平均值,通過分別測量分枝基部和末端的周長后折算得出。由于隨著分枝層級的增加,分枝直徑逐漸遞減,人工測量誤差會相應增大,因此只測量到2級分枝。即便如此,40棵樹木的人工測量在4人協(xié)作的情況下,共計耗費2 d時間,且隨著測量時間增長,測量誤差增大的概率也相應增加。

    表 1 實測數(shù)據(jù)統(tǒng)計

    人工實測樹木的分枝結構參數(shù)統(tǒng)計見表1。從實測數(shù)據(jù)里看到,樹干(0級分枝)的直徑多在15 cm左右且分布相對均勻,1級分枝和2級分枝的長度和直徑波動均較大。從0級分枝到2級分枝,直徑和長度的均值均呈現(xiàn)逐級遞減的特征。

    1.2 TLS數(shù)據(jù)獲取

    TLS數(shù)據(jù)使用徠卡BLK360地面激光掃描儀獲取。該儀器采用的激光波長為830 nm,掃描速率360 000/s個點,掃描速度3 min內(nèi)可完成一站全景掃描,掃描距離0.6~60 m,點云精度距離10 m處的點間距為6 mm,測距精度10 m范圍內(nèi)誤差為4 mm。具體結構參數(shù)見表2所示。為了獲取最佳效果的點云數(shù)據(jù),本實驗中每個測站距離樹干約2 m~3 m。采樣工作由1人獨立完成,共計耗費6 h,相比人工測量大幅節(jié)約了人力和時間,且數(shù)據(jù)更加完善。

    1.3 TLS數(shù)據(jù)預處理

    點云數(shù)據(jù)預處理主要包括多站點數(shù)據(jù)的配準和噪點的去除。掃描獲得的點云數(shù)據(jù)可通過對應軟件讀取,使用CloudCompare點云處理軟件對樹木周邊地物進行人工交互去除后,再用CSF濾波去除地面點、SOR濾波清理點云噪聲,并對每棵樹的點云進行目視檢查,確保去除噪聲時不丟失樹木的任何部分數(shù)據(jù)。最后進行單木分割。

    本文分離原始點云中地面點與非地面點時采用了布料模擬濾波算法(Cloth Simulation Filtering, CSF),算法的概述如圖2所示。該算法是由張吳明等于2016年提出,典型的優(yōu)點是適用性強和形象化,不需要根據(jù)地形來調整太多參數(shù),簡化地面濾波的處理流程[40]。Yilmaz CS等[41]、張昌賽等[42]通過將CSF算法[43]與Photoscan、Global Mapper和Lastools等多種濾波算法相比較,得出CSF算法在精度、穩(wěn)定性和適用性上均表現(xiàn)較強的結論。

    1.4 樹木分枝結構參數(shù)提取的QSM算法

    通過定量結構模型(Quantitative Structure Model, QSM)完成對欒樹分枝的重建,從而提取1、2級分枝的長度、半徑等結構參數(shù),并與實測數(shù)值進行匹配,驗證其精度與可靠度的同時,為進一步的分枝結構分析提供數(shù)據(jù)基礎。

    基于單木點云,通過QSM算法對分枝結構進行三維重建。QSM算法通過一系列圓柱來擬合分枝,將屬于同一分枝的圓柱長度之和作為分枝長度,將屬于同一分枝的第一個和最后一個圓柱的直徑平均值作為分枝直徑。完成對樹木三維結構的重建之后,基于重建樹木的分枝長度與實測樹木分枝長度來進行分枝匹配,對于每個分枝,將長度與其最接近的重建分枝作為匹配對象,若重建分枝與實測分枝的長度相對差異超過50%,則不對該分枝進行配對。

    研究采用Raumonen P等提出的用于重建樹木精細結構的定量模型QSM[29],是一種利用地面激光掃描得到的樹干和樹枝的點云,快速、自動地建立高精度樹木模型的新方法。該方法通過輸入多站點掃描的一棵樹的點云,然后建立樹的柔性圓柱模型,對樹的可見部分進行魯棒重構,并能完整的獲取樹木拓撲分支結構。主要步驟如圖3所示,其單樹點云及其重建結果如圖4所示,圖4(a)為其中1號欒樹預處理后點云,圖4(b)為初步重建結果。

    圖3 定量模型QSM主要流程圖

    圖4 QSM重建結果

    1.5 基于TLS的欒樹分形特征分析

    為進行欒樹的分枝特征分析,本文使用改進后QSM算法提取前四個級別的分枝,針對這些分枝分別提取長度、直徑、枝徑比(Ratio of Branch Diameter,RBD)、分枝率(Bifurcation Ratio,BR)和分枝角度等參數(shù),進行分形特征分析。枝徑比的計算公式如下:RBD=D+1/D(1)

    其中RBD為第級分枝的枝徑比,D+1和D分別為第+1級和第級分枝的直徑。

    分枝率包括總體分枝率(Overall Bifurcation Ratio,OBR)和逐步分枝率(Stepwise Bifurcation Ratio,SBR),總體分枝率計算公式如下:=(N-N)/(N-1) (2)

    其中N為所有級別分枝總數(shù),1為1級分枝個數(shù),N為最后1級分枝個數(shù)。

    逐步分枝率(SBR)則可根據(jù)下式計算:SBR=N/N+1(3)

    其中SBR為第級分枝的分枝率,N為第級分枝的個數(shù),N+1為第+1級分枝的個數(shù)。

    2 結果與分析

    2.1 分枝重建精度評價

    表3 分枝重建精度評價

    表3顯示運用QSM算法對不同直徑級別的分枝重建精度,可以發(fā)現(xiàn)對于分枝直徑小于5 cm、大于5 cm且小于10 cm和大于10 cm三個級別,分枝長度提取的誤差均較小,且提取結果均小于實際分枝長度,相對誤差分別為-2.88%、-4.05%和-7.15%。QSM算法對分枝直徑提取的誤差較大,對于三個直徑級別而言,提取結果均小于實際的分枝直徑,且相對誤差均在50%左右。該結論也與前人研究相一致,即QSM對于直徑小于20 cm的分枝直徑估算誤差相對較大[44]。研究通過QSM算法完成了欒樹分枝結構重建,驗證了QSM算法在直徑偏小(直徑小于20 cm)的分枝重建時,對直徑產(chǎn)生明顯的低估,且重建直徑均小于實測直徑。

    2.2 不同分枝級別的分枝個數(shù)和長度統(tǒng)計分析

    表4 不同DBH組別欒樹各級分枝個數(shù)與長度

    圖5 不同DBH級別欒樹各級分枝個數(shù)

    表4展示了不同DBH組別欒樹各級分枝的個數(shù)與長度,可以發(fā)現(xiàn)對于不同組別,隨著DBH的增加,1級分枝的長度有所增加,而其他級別的分枝沒有隨著DBH的變化產(chǎn)生明顯變化。隨著分枝級別的加深,不同DBH組別的分枝長度均呈現(xiàn)下降趨勢。就分枝個數(shù)而言,隨著DBH的增加,2級及以后級別的分枝個數(shù)均呈現(xiàn)上升趨勢;而隨著分枝級別的增加,分枝個數(shù)呈現(xiàn)先上升后下降的趨勢,如圖5所示。由此可以得出,研究區(qū)欒樹的DBH越大,1級分枝的長度越長,同時也有助于各級別分枝個數(shù)的增長。研究區(qū)欒樹正處于生長階段,因此隨著DBH的增長,1級分枝的長度增加,同時各級分枝的個數(shù)也有所上升,以此促進上部枝條對光照條件的競爭,從而促進欒樹的生長發(fā)育。

    2.3 不同分枝級別的枝徑比

    表5 不同DBH級別欒樹各級分枝枝徑比

    圖6 不同DBH級別欒樹各級分枝枝徑比

    表5和圖6展示了不同DBH級別欒樹各級分枝枝徑比的變化,可以發(fā)現(xiàn)對于不同的DBH組別,在前四級分枝,每一級分枝的枝徑比隨著DBH的增加均呈現(xiàn)下降趨勢,由此可以得出,隨著欒樹發(fā)育階段的深入,后一級分枝與前一級分枝的直徑差異在增大。在同一DBH組別內(nèi),前三級分枝直徑比逐漸增大,在DBH小于15 cm的范圍內(nèi),前四級分枝直徑比均隨著級別增加而增大,但DBH在15 cm以上的的組別里,4級分枝的枝徑比均開始下降,即3級分枝的枝徑比達峰值。通過對比不同分枝級別的枝徑比可以發(fā)現(xiàn),1級分枝的枝徑比均較低,后三級分枝的枝徑比相差不大,隨著分枝級別的加深,枝徑比呈現(xiàn)先上升后平緩的變化趨勢。由此可以得出,1級分枝通常具有較大的直徑,2級分枝的直徑與1級分枝直徑間的差異較大,隨著分枝級別的上升,后一級別與前一級別分枝的直徑差距在逐漸減小。

    2.4 不同分枝級別的分枝率

    表 6 不同DBH級別欒樹各級分枝率

    圖 7 不同DBH級別欒樹各級分枝率

    表6和圖7展示了不同DBH級別欒樹各級分枝率變化,可以發(fā)現(xiàn)對于不同級別,1級分枝與其他幾個級別的分枝相比具有較高的分枝率,2、3、4級分枝的分枝率相互差別比較小,但隨著分枝級別的上升,分枝率總體呈現(xiàn)下降的趨勢。對于不同的DBH組別,隨著發(fā)育階段的加深,1、2、3、4級分枝的分枝率隨著DBH的增加整體呈現(xiàn)上升趨勢,尤其是1級分枝的分枝率增長幅度較大。由此可見,在欒樹的生長發(fā)育階段,隨著DBH的增長,有助于上級分枝分離出更多的下級分枝,從而增強對光照的競爭,促進樹木整體的發(fā)育??傮w上,隨著DBH的增加,研究區(qū)欒樹的分枝率呈下降趨勢。

    2.5 不同分枝級別的分枝角度

    表7 不同DBH級別欒樹各級分枝角度

    圖8 不同DBH級別欒樹各級分枝角度

    分枝角度對于樹冠的形態(tài)具有重要作用,樹木分枝在樹冠中的分布方向影響著樹冠對光照和CO2等條件的利用。圖8和表7展示了不同DBH級別欒樹各級分枝角度的變化,可以發(fā)現(xiàn)對于不同的分枝級別,隨著級別的增加,分枝角度呈現(xiàn)上升趨勢,到高級別分枝的變化則較平緩,以便更好地利用生長空間并吸收光照。對于不同的DBH組別,能夠發(fā)現(xiàn)隨著DBH的增加,分枝角度呈現(xiàn)先上升后下降的趨勢。對于不同的分枝界別而言,DBH在15~17 cm的組別具有最大的分枝角度。由此可以得出,在樹木生長的前期分枝角度逐漸增大,隨著樹木生長過程的深入,分枝角度又逐漸變小,這有助于分枝更好地利用空間并吸收光照。

    3 討論

    本文基于欒樹的TLS點云使用QSM算法提取前四級分枝的結構參數(shù),進行了欒樹分枝的分形特征分析。結果表明,研究區(qū)的欒樹正處于生長階段,隨著發(fā)育階段的深入,欒樹的DBH增長,1級分枝的長度也隨之增加,同時各級分枝的個數(shù)也有所上升。后一級分枝與前一級分枝的直徑差異在增大,枝徑比上升。欒樹前四級分枝的分枝率隨著DBH的增加表顯示出上升趨勢,其中1級分枝的分枝率具有較大的增長幅度,分枝角度則隨著DBH的增加表現(xiàn)出先上升后下降的趨勢。以上分枝格局有助于欒樹分離出更多的下級分枝,增強對光照的競爭和對發(fā)育空間的利用,從而促進對樹木整體的發(fā)育。探討欒樹的分枝格局將有助于增進其對光照和生長空間的利用情況,并揭示其應對環(huán)境變化的響應,從而有助于進行對欒樹的精確管理和科學經(jīng)營。應當注意的是,本文選擇的欒樹為同齡10年樹,均處于生長發(fā)育階段,不同樹木之間結構的差異性較小,因此對比度更高,但是在更大范圍進行數(shù)據(jù)采集時,可能就會面臨對比難度增大的問題。同時,盡管QSM算法在分枝參數(shù)提取上有很強的能力,但是實驗證明其在分枝直徑方面的誤差較大,因此后續(xù)研究需要對如何提升直徑的提取能力以及對長度等參數(shù)的精度提高上作出努力。再者,如果可以針對同采樣區(qū)的欒樹進行年度的跟蹤研究,則可以更加系統(tǒng)科學的了解不同年齡的欒樹分形特征以及隨樹齡增長所產(chǎn)生的動態(tài)變化,后從而在更廣泛的角度討論TLS在進行樹木分枝結構研究方面的應用潛力。

    4 結論

    本實驗在開展過程中,發(fā)現(xiàn)人工測量需要人員多、耗時長、測量級別少的缺點,且如果對較大樹木進行測量,非破壞難以完成,而地基激光雷達相對節(jié)約大量人力和時間。本文基于QSM算法的重建結果,進行了欒樹冠層分形特征的分析,通過對比不同組別欒樹分枝特征的差異,發(fā)現(xiàn)在欒樹的生長發(fā)育階段,隨著DBH的增加,欒樹1級分枝的長度上升,同時具有更多的次級分枝數(shù)目,分枝率和枝徑比也同步上升。而隨著DBH的增加,各級分枝角度均呈現(xiàn)先上升后下降的趨勢,以更充分地利用光照條件和生長空間,促進欒樹的發(fā)育。因此,本研究證明基于TLS點云能夠實現(xiàn)對欒樹分枝結構參數(shù)的精確獲取,從而開展多級別分枝的分形特征研究,大幅節(jié)省人工測量所需要的大量人力物力,同時增進對欒樹分枝結構與生態(tài)屬性聯(lián)系的理解,為樹木的精細化管理與科學經(jīng)營提供關鍵的決策依據(jù)。

    [1] Ohsawa M, Nitta I. Patterning of subtropical warm-temperate evergreen broad-leaved forests in east Asian mountains with special reference to shoot phenology [J]. Tropics, 1997,6(4):317-334

    [2] Wilkerson GG. Plant growth modeling and applications [J]. Agricultural Systems, 2004,82(1):93-94

    [3] 孫棟元,趙成義,王麗娟,等.荒漠植物構型研究進展[J].水土保持研究,2011(5):281-287

    [4] 徐程揚.不同光環(huán)境下紫椴幼樹樹冠結構的可塑性響應[J].應用生態(tài)學報,2001,12(3):339-343

    [5] 李俊清,臧潤國,蔣有緒.歐洲水青岡(L.)構筑型與形態(tài)多樣性研究[J].生態(tài)學報,2001(1):151-155

    [6] 陳波,宋永昌,達良俊.木本植物的構型及其在植物生態(tài)學研究的進展[J].生態(tài)學雜志,2002(3):53-57,29

    [7] 鐘章成,曾波.植物種群生態(tài)研究進展[J].西南師范大學學報:自然科學版,2001,26(2):230-236

    [8] 黃文娟,李志軍,楊趙平,等,胡楊異形葉結構型性狀及其與胸徑關系[J].生態(tài)學雜志,2010,29(12):2347-2352

    [9] 孫書存,陳靈芝.遼東櫟植冠的構型分析[J].植物生態(tài)學報,1999,23(5):433-440

    [10] 張文輝,祖元剛,馬克明.裂葉沙參與泡沙參種群分布格局分形特征的分析[J].植物生態(tài)學報,1999,23(1):31-39

    [11] 馬克明,祖元剛.興安落葉松分枝格局的分形特征[J].植物研究,2000,20(2):235-241

    [12] 趙相健,王孝安.太白紅杉分枝格局的可塑性研究[J].西北植物學報,2005,25(1):113-117

    [13] 林勇明,俞偉,劉奕,等.不同距海處木麻黃分枝格局特征及冠形分析[J].應用與環(huán)境生物學報,2013,19(4):587-592

    [14] Grau E. Estimation of 3D vegetation density with Terrestrial Laser Scanning data using voxels. A sensitivity analysis of influencing parameters [J]. Remote Sensing of Environment, 2017,191:373-388

    [15] Malhi , Jackson, Bentley,. New perspectives on the ecology of tree structure and tree communities through terrestrial laser scanning [J]. Interface Focus, 2018,8(2):20170052

    [16] Bazezew MN, Hussin YA, Kloosterman EH. Integrating Airborne LiDAR and Terrestrial Laser Scanner forest parameters for accurate above-ground biomass/carbon estimation in Ayer Hitam tropical forest, Malaysia [J]. International Journal of Applied Earth Observation & Geoinformation, 2018,73:638-652

    [17] Wilkes P, Lau A, Disney M,. Data acquisition considerations for terrestrial laser scanning of forest plots [J]. Remote Sensing of Environment, 2017,196:140-153

    [18] Takoudjou SM, Ploton P, Sonk? B,. Using terrestrial laser scanning data to estimate large tropical trees biomass and calibrate allometric models: A comparison with traditional destructive approach [J]. Methods in Ecology & Evolution, 2017,9(4):905-916

    [19] Stovall, Vorster, Anderson,. Non-destructive aboveground biomass estimation of coniferous trees using terrestrial LiDAR [J]. Remote Sensing of Environment, 2017,200:31-42

    [20] Tansey K, Semles N, Anstee A,. Estimating tree and stand variables in a Corsican Pine woodland from terrestrial. laser scanner data [J]. International Journal of Remote Sensing, 2009,30(19):5195-5209

    [21] Huang H, Li Z, Gong P,. Automated methods for measuring DBH and tree heights with a commercial scanning。Lidar [J]. Photogrammetric Engineering & Remote Sensing, 2011,77(3):219-227

    [22] Burt A, Disney MI, Armston P,. Rapid characterisation of forest structure from TLS and 3D modelling. in. Geoscience and Remote Sensing Symposium (IGARSS) [J]. IEEE International. 2013, PP: 3387-3390

    [23] Krooks A , Kaasalainen S , Kankare V, e. Predicting tree structure from tree height using terrestrial laser scanning. and quantitative structure models [J]. Silva Fenn, 2014, 48: 1125

    [24] Brede B, Lau A, Bartholomeus HM,. Comparing RIEGL RiCOPTER UAV LiDAR derived canopy height and DBH with terrestrial LiDAR [J]. Sensors, 2017,17:2371

    [25] Zhao F, Guo Q, Kelly M. Allometric equation choice impacts lidar-based forest biomass estimates: A case study from the Sierra National Forest, CA [J]. Agricultural & Forest Meteorology, 2012,165:64-72

    [26] Srinivasan S, Popescu SC, Eriksson M, et al. Terrestrial laser scanning as an effective tool to retrieve tree level height, crown width, and stem diameter [J].Remote Sensing, 2015,7(2): 1877-1896

    [27] Dassot M, Colin A, Santenoise P,. Terrestrial laser scanning for measuring the solid wood volume, including branches, of adult standing trees in the forest environment [J]. Computers & Electronics in Agriculture, 2012,89:86-93

    [28] R?der A,Klemmt HJ., Seifert S,Integration of the forest growth simulator SILVA in practical forest management. EnviroInfo (1). 2007:239-246

    [29] Raumonen P, Kaasalainen M, ?kerblom M,. Fast automatic precision tree models from terrestrial laser scanner data [J]. Remote Sensing, 2013,5(2):491-520

    [30] Lau A, Martius C., Bartholomeus H, et. Estimating architecture-based metabolic scaling exponents of tropical trees using terrestrial LiDAR and 3D modelling [J]. Forest Ecology and Management, 2019,439:132-145

    [31] Chen S, Feng Z., Chen P,. Nondestructive estimation of the above-ground biomass of multiple tree species in boreal forests of China using Terrestrial Laser Scanning [J]. Forests, 2019,10(11):936-961

    [32] 倪文儉,過志峰,孫國清,等.基于地基激光雷達數(shù)據(jù)的單木結構參數(shù)提取研究[J].高技術通訊,2010,20(2):191-198

    [33] 王寧寧,尹文廣,黃秦軍,等.三維掃描技術在獲取楊樹樹冠結構特征參數(shù)上的應用[J].林業(yè)科學,2015,51(5):108-116

    [34] 劉金鵬,張懷清,劉閩,等.基于點云骨架的單木分枝參數(shù)提取研究[J].北京林業(yè)大學學報,2016,38(10):15-20

    [35] 楊吉華,張永濤,王貴霞,等.欒樹、黃連木、黃櫨水分生理生態(tài)特性的研究[J].水土保持學報,2002(4):152-154

    [36] 張勝,湯少展,彭兆豐.欒樹中有效成分的提取、應用及其生物修復綜述[J].江蘇林業(yè)科技,2020,47(2):48-51

    [37] 馮獻賓,董倩,李旭新,等.黃連木和黃山欒樹的抗寒性[J].應用生態(tài)學報,2011,22(5):1141-1146

    [38] 馬廣恩,申雅維,魯學照,等.欒樹抗菌有效成分的研究[J].中草藥,1998(2):84-85

    [39] 翟梅枝,郭景麗,王磊,等.欒樹花黃酮類化合物的提取工藝研究[J].西北林學院學報,2010,25(2):136-139

    [40] Zhao JQ, Ledoux H, Stoter J,. HSW: Heuristic Shrink-wrapping for automatically repairing solid-based City GML LOD2 building models [J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2018,146:289-304

    [41] Yilmaz CS, Yilmaz VO, Güng R. Investigating the performances of commercial and non-commercial software for ground filtering of UAV-based point clouds [J]. International Journal of Remote Sensing, 2018:5016-5042

    [42] 張昌賽,劉正軍,楊樹文,等.基于LiDAR數(shù)據(jù)的布料模擬濾波算法的適用性分析[J].激光技術,2018,42(3):410-416

    [43] Zhang WM, Qi JB, Wan P,. An easy-to-use airborne LiDAR data filtering method based on cloth simulation [J]. Remote Sensing, 2016,8(6):501

    [44] Lau A, Bentley LP, Martius C,. Quantifying branch architecture of tropical trees using terrestrial LiDAR and 3D modelling [J]. Trees, 2018,32(5):1219-1231

    Fractal Characteristics ofon Ground-based LiDAR

    LI Hui, LIN Yi*, MENG Xiang-shuang, SHI Zhen-wei, CAI Wan-yuan

    100871,

    The geometric characteristics of tree meristems describes the structures that branches form according to certain laws. It reflects the adaptability of individual tree to the environment, and the competitive power to factors that affect growth such as light, water and growth space.The advancement of Terrestrial Laser Scanning (TLS) technology has made it possible to obtain parameters of single wood constructions in a convenient, non-destructive, and environmentally friendly manner. In this study, we used TLS point cloud based method instead of original labor-intensive methods to estimate structural parameters of tree branches, and perform further analysis of fractal characteristics. The results showed that to make better use of light and growth space during the growth stage, branching ratio and branch diameter ratio rised synchronously as DBH icreasd,while branching angles increased first, and then decreased。

    LiDAR; Koelreuteria paniculata; characteristics analysis

    P237

    A

    1000-2324(2022)03-0475-09

    10.3969/j.issn.1000-2324.2022.03.021

    2021-12-14

    2022-02-07

    國家自然科學基金面上項目(32171782,31870531)

    李輝(1983-),男,碩士研究生,講師,主要從事激光雷達遙感研究. E-mail:Hui_li@pku.edu.cn

    Author for correspondence. E-mail:yi.lin@pku.edu.cn

    猜你喜歡
    欒樹分枝分形
    一株吊蘭
    感受分形
    隨季節(jié)而變奏的欒樹
    分形之美
    我和欒樹有個約會
    我和欒樹有個約會
    分形空間上廣義凸函數(shù)的新Simpson型不等式及應用
    帶移民和拯救的二次加權分枝過程的有關性質
    受控兩性分枝過程
    欒樹
    長江文藝(2017年2期)2017-02-13 19:17:24
    大片免费播放器 马上看| 狂野欧美激情性xxxx在线观看| 中文资源天堂在线| 精品一品国产午夜福利视频| 丁香六月天网| 女的被弄到高潮叫床怎么办| 欧美日韩亚洲高清精品| 久久久久久久国产电影| 成人亚洲精品一区在线观看| 国产一区有黄有色的免费视频| 黑人巨大精品欧美一区二区蜜桃 | 国产成人精品一,二区| 亚洲不卡免费看| 涩涩av久久男人的天堂| 另类精品久久| 黑人巨大精品欧美一区二区蜜桃 | 国产综合精华液| 9色porny在线观看| 亚洲欧洲精品一区二区精品久久久 | 成人国产麻豆网| 91久久精品电影网| av播播在线观看一区| 精华霜和精华液先用哪个| 日韩精品有码人妻一区| 亚洲av中文av极速乱| 在线天堂最新版资源| 中文字幕人妻丝袜制服| 人体艺术视频欧美日本| 精品亚洲成a人片在线观看| 久久国产亚洲av麻豆专区| 亚洲一区二区三区欧美精品| 精品国产露脸久久av麻豆| 日韩一本色道免费dvd| 亚洲在久久综合| 国产一区二区三区综合在线观看 | 9色porny在线观看| 男的添女的下面高潮视频| 国产在线视频一区二区| 中文字幕免费在线视频6| 国产在视频线精品| 最近最新中文字幕免费大全7| 国产成人精品无人区| 99re6热这里在线精品视频| 男女无遮挡免费网站观看| 国产精品久久久久成人av| 久久 成人 亚洲| 久久人人爽av亚洲精品天堂| av免费在线看不卡| 免费播放大片免费观看视频在线观看| 国产熟女午夜一区二区三区 | 男人狂女人下面高潮的视频| 日韩亚洲欧美综合| 校园人妻丝袜中文字幕| 亚洲激情五月婷婷啪啪| 伊人久久精品亚洲午夜| 国产伦在线观看视频一区| 国产日韩欧美亚洲二区| 51国产日韩欧美| 一级毛片我不卡| 亚洲精品日韩在线中文字幕| h日本视频在线播放| 97在线人人人人妻| 啦啦啦啦在线视频资源| 亚洲综合精品二区| 欧美日韩视频高清一区二区三区二| 黑人高潮一二区| 国产在线视频一区二区| 欧美激情国产日韩精品一区| 男女无遮挡免费网站观看| 乱人伦中国视频| 久久精品熟女亚洲av麻豆精品| 精品99又大又爽又粗少妇毛片| 狂野欧美激情性bbbbbb| 久久久久国产精品人妻一区二区| 三上悠亚av全集在线观看 | 国产乱人偷精品视频| 简卡轻食公司| 性色av一级| 日日摸夜夜添夜夜爱| 欧美国产精品一级二级三级 | 欧美精品亚洲一区二区| 美女xxoo啪啪120秒动态图| av免费观看日本| av播播在线观看一区| 久久国产精品大桥未久av | 欧美最新免费一区二区三区| 成人18禁高潮啪啪吃奶动态图 | 欧美区成人在线视频| 婷婷色麻豆天堂久久| 99久久精品热视频| 亚洲精品乱久久久久久| 黄色怎么调成土黄色| 亚洲av欧美aⅴ国产| 精华霜和精华液先用哪个| 噜噜噜噜噜久久久久久91| 久久午夜福利片| 三级国产精品片| 亚洲不卡免费看| 日韩成人伦理影院| 亚洲电影在线观看av| 美女主播在线视频| 国产精品免费大片| 国产视频首页在线观看| 国产在视频线精品| 亚洲精品456在线播放app| 亚洲av免费高清在线观看| 晚上一个人看的免费电影| 亚洲三级黄色毛片| 精品人妻熟女毛片av久久网站| 赤兔流量卡办理| 久久午夜福利片| 春色校园在线视频观看| 一本大道久久a久久精品| 国产视频内射| 高清欧美精品videossex| 简卡轻食公司| 熟妇人妻不卡中文字幕| 在线观看www视频免费| 性色avwww在线观看| 久久久国产一区二区| 久久6这里有精品| 丰满人妻一区二区三区视频av| 女人精品久久久久毛片| 国产极品天堂在线| 99热国产这里只有精品6| 国产精品偷伦视频观看了| 永久免费av网站大全| 日日啪夜夜撸| 日本免费在线观看一区| 街头女战士在线观看网站| 亚洲电影在线观看av| 久久亚洲国产成人精品v| 国产一区有黄有色的免费视频| 午夜免费男女啪啪视频观看| 人人妻人人看人人澡| 免费人成在线观看视频色| 中文字幕免费在线视频6| 欧美xxxx性猛交bbbb| 91精品国产九色| 99热全是精品| 三级国产精品片| 亚洲成人一二三区av| 夜夜爽夜夜爽视频| 亚洲三级黄色毛片| 国产淫语在线视频| 久久精品久久久久久久性| 91成人精品电影| 亚洲第一av免费看| 国产日韩欧美亚洲二区| 乱码一卡2卡4卡精品| 久久精品久久久久久久性| 精品99又大又爽又粗少妇毛片| 国产乱来视频区| 狂野欧美白嫩少妇大欣赏| 一级av片app| 人妻系列 视频| 亚洲精品乱码久久久v下载方式| 日日啪夜夜爽| 老熟女久久久| tube8黄色片| 一个人免费看片子| 内地一区二区视频在线| 中文字幕制服av| 秋霞伦理黄片| 嘟嘟电影网在线观看| 国产成人aa在线观看| 美女xxoo啪啪120秒动态图| a级毛片免费高清观看在线播放| 亚洲怡红院男人天堂| 成人特级av手机在线观看| 97在线视频观看| 久久久久人妻精品一区果冻| 免费久久久久久久精品成人欧美视频 | 午夜福利在线观看免费完整高清在| 观看美女的网站| 欧美性感艳星| 日本欧美国产在线视频| 麻豆成人午夜福利视频| 黄色毛片三级朝国网站 | 亚洲真实伦在线观看| 永久免费av网站大全| 在线天堂最新版资源| 日本免费在线观看一区| 伊人亚洲综合成人网| 日产精品乱码卡一卡2卡三| 美女国产视频在线观看| 亚洲欧美成人综合另类久久久| 日韩大片免费观看网站| 大片电影免费在线观看免费| 欧美国产精品一级二级三级 | 国产日韩欧美亚洲二区| 日本91视频免费播放| 国产成人精品一,二区| 乱人伦中国视频| 国产精品不卡视频一区二区| www.av在线官网国产| 国产在线视频一区二区| 亚洲经典国产精华液单| 国产av精品麻豆| av福利片在线观看| freevideosex欧美| 亚洲精品久久久久久婷婷小说| 免费人妻精品一区二区三区视频| 3wmmmm亚洲av在线观看| 色婷婷av一区二区三区视频| 嫩草影院入口| 大香蕉97超碰在线| 热99国产精品久久久久久7| 在现免费观看毛片| 少妇熟女欧美另类| 九九爱精品视频在线观看| 国产日韩欧美亚洲二区| 大又大粗又爽又黄少妇毛片口| 久久精品熟女亚洲av麻豆精品| 免费少妇av软件| 日韩在线高清观看一区二区三区| 三级经典国产精品| 欧美老熟妇乱子伦牲交| 一个人看视频在线观看www免费| 亚洲中文av在线| 99re6热这里在线精品视频| 欧美97在线视频| 日本欧美视频一区| 午夜福利乱码中文字幕| 精品一品国产午夜福利视频| 欧美xxⅹ黑人| 永久免费av网站大全| 精品福利永久在线观看| 美女脱内裤让男人舔精品视频| 黄色 视频免费看| 精品少妇内射三级| 久热爱精品视频在线9| 啦啦啦在线免费观看视频4| 永久免费av网站大全| 亚洲成国产人片在线观看| 91九色精品人成在线观看| 美女中出高潮动态图| 午夜福利在线观看吧| av天堂久久9| 久久久久国产精品人妻一区二区| 一级a爱视频在线免费观看| 国产主播在线观看一区二区| 日韩大码丰满熟妇| 后天国语完整版免费观看| 欧美黄色淫秽网站| 在线观看舔阴道视频| 国产精品免费视频内射| 日日夜夜操网爽| 精品少妇久久久久久888优播| 欧美乱码精品一区二区三区| 亚洲av成人不卡在线观看播放网 | 天天躁狠狠躁夜夜躁狠狠躁| 一本大道久久a久久精品| 女人爽到高潮嗷嗷叫在线视频| 老司机深夜福利视频在线观看 | 三上悠亚av全集在线观看| 久久久久久久大尺度免费视频| 高清视频免费观看一区二区| 99国产综合亚洲精品| 少妇裸体淫交视频免费看高清 | 精品久久久精品久久久| 又黄又粗又硬又大视频| 亚洲精品一卡2卡三卡4卡5卡 | 成人免费观看视频高清| 深夜精品福利| 狠狠狠狠99中文字幕| 俄罗斯特黄特色一大片| e午夜精品久久久久久久| 黑人猛操日本美女一级片| 精品第一国产精品| 国产免费一区二区三区四区乱码| 下体分泌物呈黄色| 一边摸一边做爽爽视频免费| 欧美变态另类bdsm刘玥| 国产一区二区三区综合在线观看| 欧美激情高清一区二区三区| 69精品国产乱码久久久| 自拍欧美九色日韩亚洲蝌蚪91| 桃花免费在线播放| 婷婷丁香在线五月| 黄片播放在线免费| 亚洲成人手机| 一本大道久久a久久精品| 少妇粗大呻吟视频| 久久亚洲国产成人精品v| 日韩欧美免费精品| 男女之事视频高清在线观看| 成年人黄色毛片网站| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲av电影在线观看一区二区三区| 久久久久精品国产欧美久久久 | 亚洲七黄色美女视频| 久久久久国内视频| 亚洲色图 男人天堂 中文字幕| 亚洲av男天堂| 黄色怎么调成土黄色| 成人三级做爰电影| 美女扒开内裤让男人捅视频| 老司机福利观看| 欧美日韩亚洲国产一区二区在线观看 | 久久女婷五月综合色啪小说| 国内毛片毛片毛片毛片毛片| 精品国产乱码久久久久久小说| 成人国产av品久久久| 午夜久久久在线观看| 视频区图区小说| 桃花免费在线播放| 女性被躁到高潮视频| 国产野战对白在线观看| 久久精品熟女亚洲av麻豆精品| 岛国毛片在线播放| 国产成人精品久久二区二区91| 久久九九热精品免费| 黄色a级毛片大全视频| 搡老熟女国产l中国老女人| 九色亚洲精品在线播放| 亚洲欧美日韩另类电影网站| 91成人精品电影| 免费观看av网站的网址| 亚洲av国产av综合av卡| 菩萨蛮人人尽说江南好唐韦庄| av在线播放精品| 久久人人爽人人片av| 国产xxxxx性猛交| 欧美日本中文国产一区发布| 久久av网站| 精品免费久久久久久久清纯 | 狂野欧美激情性xxxx| 中文字幕最新亚洲高清| 精品高清国产在线一区| 亚洲精品一区蜜桃| 每晚都被弄得嗷嗷叫到高潮| 精品一区二区三卡| 99精品欧美一区二区三区四区| 国产淫语在线视频| 捣出白浆h1v1| 久久精品国产亚洲av高清一级| 青草久久国产| 看免费av毛片| 亚洲精品国产色婷婷电影| 又黄又粗又硬又大视频| 精品一区二区三区av网在线观看 | 97在线人人人人妻| 人妻一区二区av| 中文字幕人妻丝袜一区二区| 久久狼人影院| 18禁国产床啪视频网站| 每晚都被弄得嗷嗷叫到高潮| 天天躁夜夜躁狠狠躁躁| 亚洲,欧美精品.| 国产成人精品无人区| 菩萨蛮人人尽说江南好唐韦庄| 丝袜人妻中文字幕| 欧美日韩福利视频一区二区| 丁香六月欧美| 欧美人与性动交α欧美软件| 女人被躁到高潮嗷嗷叫费观| 嫩草影视91久久| 69av精品久久久久久 | 亚洲中文av在线| 午夜福利在线免费观看网站| 精品熟女少妇八av免费久了| 亚洲av电影在线观看一区二区三区| 欧美 亚洲 国产 日韩一| 国产精品影院久久| 他把我摸到了高潮在线观看 | 国产一区二区 视频在线| 国产成人免费无遮挡视频| 国产亚洲av片在线观看秒播厂| a级片在线免费高清观看视频| 久久人妻熟女aⅴ| 美女脱内裤让男人舔精品视频| 一本—道久久a久久精品蜜桃钙片| 中文字幕高清在线视频| 亚洲av成人一区二区三| 国产精品久久久久久精品古装| 久久国产精品大桥未久av| 美女主播在线视频| 黄色怎么调成土黄色| 久久久精品区二区三区| 久久精品熟女亚洲av麻豆精品| 男女边摸边吃奶| 国产精品成人在线| 女人精品久久久久毛片| 一区二区三区激情视频| 免费黄频网站在线观看国产| 成人免费观看视频高清| 亚洲av日韩在线播放| 女警被强在线播放| 国产精品熟女久久久久浪| 国产亚洲欧美在线一区二区| www.精华液| 韩国精品一区二区三区| 亚洲国产欧美一区二区综合| 亚洲 欧美一区二区三区| 精品一区二区三卡| 高清黄色对白视频在线免费看| 日本猛色少妇xxxxx猛交久久| 三上悠亚av全集在线观看| 成人亚洲精品一区在线观看| 欧美精品一区二区大全| 精品卡一卡二卡四卡免费| 涩涩av久久男人的天堂| 精品国产乱码久久久久久男人| 欧美日韩视频精品一区| 国产一区二区三区综合在线观看| 成人手机av| 亚洲成国产人片在线观看| 亚洲精品成人av观看孕妇| 考比视频在线观看| 亚洲激情五月婷婷啪啪| 中文字幕最新亚洲高清| 亚洲精品中文字幕在线视频| 国产精品二区激情视频| 国产av精品麻豆| 黄片播放在线免费| av不卡在线播放| 久久人人爽人人片av| 亚洲国产精品999| 又大又爽又粗| 日韩制服骚丝袜av| 一区二区三区乱码不卡18| videos熟女内射| 俄罗斯特黄特色一大片| 色综合欧美亚洲国产小说| 纯流量卡能插随身wifi吗| 国产成人精品在线电影| 国产成人一区二区三区免费视频网站| 午夜激情久久久久久久| 91精品伊人久久大香线蕉| 婷婷色av中文字幕| 亚洲欧美精品自产自拍| 中文字幕色久视频| 亚洲国产精品成人久久小说| 日韩电影二区| 精品国产一区二区三区四区第35| 狠狠婷婷综合久久久久久88av| 久9热在线精品视频| 亚洲午夜精品一区,二区,三区| 国产精品久久久人人做人人爽| 亚洲精品中文字幕一二三四区 | 午夜福利在线免费观看网站| 久久免费观看电影| 伦理电影免费视频| 国产高清视频在线播放一区 | 亚洲专区中文字幕在线| 午夜福利在线免费观看网站| 一进一出抽搐动态| 中国国产av一级| 又黄又粗又硬又大视频| 老司机靠b影院| 免费看十八禁软件| 两个人看的免费小视频| 女人被躁到高潮嗷嗷叫费观| 欧美激情高清一区二区三区| 久久久久国产一级毛片高清牌| 精品人妻在线不人妻| www.av在线官网国产| 十八禁人妻一区二区| 精品视频人人做人人爽| 满18在线观看网站| 国产精品 欧美亚洲| 免费日韩欧美在线观看| 国产日韩欧美在线精品| 韩国精品一区二区三区| 精品久久久久久久毛片微露脸 | 久久国产精品大桥未久av| 最近最新中文字幕大全免费视频| 欧美精品一区二区免费开放| 欧美黄色淫秽网站| 欧美精品啪啪一区二区三区 | 亚洲五月婷婷丁香| 亚洲国产中文字幕在线视频| 男人舔女人的私密视频| av欧美777| 中亚洲国语对白在线视频| 国产成人欧美在线观看 | 桃花免费在线播放| 久久久久久久大尺度免费视频| 久久国产精品大桥未久av| 久久久久视频综合| 老司机影院成人| 精品视频人人做人人爽| 亚洲伊人久久精品综合| 午夜福利影视在线免费观看| 成年动漫av网址| 美女大奶头黄色视频| 脱女人内裤的视频| 99国产极品粉嫩在线观看| 啪啪无遮挡十八禁网站| 欧美一级毛片孕妇| 蜜桃国产av成人99| 国内毛片毛片毛片毛片毛片| 色综合欧美亚洲国产小说| 五月天丁香电影| 久久女婷五月综合色啪小说| 成人国产av品久久久| 天堂中文最新版在线下载| 精品久久久久久电影网| 两性午夜刺激爽爽歪歪视频在线观看 | 视频区欧美日本亚洲| 亚洲国产欧美在线一区| 女人高潮潮喷娇喘18禁视频| 一级,二级,三级黄色视频| 91大片在线观看| 亚洲专区国产一区二区| 热99久久久久精品小说推荐| www.精华液| 韩国精品一区二区三区| 丝袜美腿诱惑在线| 中文欧美无线码| 国产成人精品久久二区二区免费| 青春草视频在线免费观看| 高清黄色对白视频在线免费看| 国产男人的电影天堂91| 自线自在国产av| 亚洲精品中文字幕一二三四区 | 国产97色在线日韩免费| 青春草亚洲视频在线观看| av在线app专区| 免费女性裸体啪啪无遮挡网站| e午夜精品久久久久久久| 黑人猛操日本美女一级片| 国产主播在线观看一区二区| www.精华液| 热99久久久久精品小说推荐| 精品一区在线观看国产| 日韩精品免费视频一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 动漫黄色视频在线观看| 国产一区二区三区av在线| 水蜜桃什么品种好| 欧美日本中文国产一区发布| 久久久久久久久免费视频了| 美女视频免费永久观看网站| 视频在线观看一区二区三区| netflix在线观看网站| 亚洲欧美色中文字幕在线| 国产精品亚洲av一区麻豆| 女人爽到高潮嗷嗷叫在线视频| 久久久国产精品麻豆| 免费久久久久久久精品成人欧美视频| avwww免费| 99热网站在线观看| 一个人免费看片子| 国产精品成人在线| 婷婷成人精品国产| 人人妻人人澡人人爽人人夜夜| 各种免费的搞黄视频| 美女福利国产在线| 亚洲少妇的诱惑av| 一区二区av电影网| 热99国产精品久久久久久7| 国产欧美日韩一区二区三区在线| 人成视频在线观看免费观看| 啦啦啦在线免费观看视频4| 国产精品99久久99久久久不卡| 欧美激情极品国产一区二区三区| tocl精华| 不卡av一区二区三区| 日韩欧美一区视频在线观看| 精品久久久久久电影网| 男人爽女人下面视频在线观看| 中文字幕人妻丝袜制服| 国产亚洲欧美精品永久| 91字幕亚洲| 操出白浆在线播放| 国产精品一区二区精品视频观看| 国产精品久久久久久人妻精品电影 | 欧美精品一区二区大全| 黄片播放在线免费| 超碰97精品在线观看| 国产不卡av网站在线观看| 国产高清国产精品国产三级| 天天躁夜夜躁狠狠躁躁| av超薄肉色丝袜交足视频| 狂野欧美激情性bbbbbb| 五月开心婷婷网| 中亚洲国语对白在线视频| 日韩中文字幕视频在线看片| 777久久人妻少妇嫩草av网站| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品一二三区在线看| 日本av免费视频播放| 欧美大码av| 无限看片的www在线观看| 精品免费久久久久久久清纯 | 国产精品麻豆人妻色哟哟久久| 涩涩av久久男人的天堂| 亚洲成人国产一区在线观看| 纵有疾风起免费观看全集完整版| 欧美 亚洲 国产 日韩一| 人妻 亚洲 视频| 日本vs欧美在线观看视频| 国产精品久久久久久精品电影小说| 青春草视频在线免费观看| 国产有黄有色有爽视频| 欧美性长视频在线观看| 日韩,欧美,国产一区二区三区| 69精品国产乱码久久久| 侵犯人妻中文字幕一二三四区| 99re6热这里在线精品视频| tocl精华| 老熟妇仑乱视频hdxx| 99re6热这里在线精品视频| 日韩一卡2卡3卡4卡2021年| 9热在线视频观看99| 国产人伦9x9x在线观看| 国产日韩欧美亚洲二区| 一级黄色大片毛片| 男人添女人高潮全过程视频| www.熟女人妻精品国产| av线在线观看网站| 亚洲色图 男人天堂 中文字幕| 丰满饥渴人妻一区二区三|