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

    基于L1中值骨架提取的植物莖干補全研究

    2022-11-29 10:10:16姜金岑李聯(lián)隊王美麗
    關(guān)鍵詞:莖干連接點軸線

    姜金岑,李聯(lián)隊,王美麗,3,4*

    (1.西北農(nóng)林科技大學(xué)信息工程學(xué)院,陜西 楊凌 712100;2.陜西省林業(yè)科學(xué)院,陜西 西安 710082;3.農(nóng)業(yè)農(nóng)村部農(nóng)業(yè)物聯(lián)網(wǎng)重點實驗室,陜西 楊凌 712100;4.陜西省農(nóng)業(yè)信息與智能服務(wù)重點實驗室,陜西 楊凌 712100)

    植物的數(shù)字化與可視化技術(shù)是實現(xiàn)數(shù)字林業(yè)的重要途徑,對植物進(jìn)行可視化建模能夠?qū)χ参锏纳L過程和形態(tài)結(jié)構(gòu)進(jìn)行定量化研究和智能化處理。隨著三維掃描設(shè)備的普及與技術(shù)的發(fā)展,基于植物點云數(shù)據(jù)的三維模型重建越來越受到研究者的重視,且具有精度高、真實感強以及細(xì)節(jié)特征顯著等優(yōu)點[1]。在對植物進(jìn)行三維掃描過程中,由于外界環(huán)境的干擾、視點選擇的不同,掃描得到的三維點云通常存在一定的噪聲,噪聲會使得模型失真,同時也會產(chǎn)生被測物體部分缺失的現(xiàn)象,直接對其進(jìn)行重建將無法得到一個完整的模型。因此植物點云數(shù)據(jù)的處理及補全是植物三維建模中的重點和難點問題。

    通過三維激光掃描技術(shù)對植物進(jìn)行建模的研究一直都是計算機圖形學(xué)的研究重點[2-7]。此外,點云分割是三維點云重建中的關(guān)鍵技術(shù)之一,也是三維重建的前提[8-13]。但重建過程中會發(fā)現(xiàn)掃描得到的三維點云模型常缺失部分細(xì)節(jié)信息,若需進(jìn)行點云修補往往需要提取其骨架[14-19]。其中,Huang等[15]提出了基于 L1 中值的點云骨架提取方法和局部密度權(quán)重來解決點云密度不均勻的問題,但該方法需要調(diào)試大量參數(shù)來獲取精確的骨架,實用性不高;Zhang等[17]使用概率圖模型基于最小距離對骨架節(jié)點間邊的權(quán)重進(jìn)行計算,通過最小生成樹判斷每個節(jié)點間的拓?fù)溥B接情況。但通常這些方法僅適用于人體或其他簡單物體點云數(shù)據(jù),無法對植物點云數(shù)據(jù)的骨架進(jìn)行拓?fù)溥B接。

    近年來,Charles等[20]提出了可以直接處理原始三維點云數(shù)據(jù)的神經(jīng)網(wǎng)絡(luò)PointNet,能完成高層次的點云分類和語義分割任務(wù)。在此基礎(chǔ)上許多研究者通過深度學(xué)習(xí)神經(jīng)網(wǎng)絡(luò)對點云進(jìn)行修補[21-26],但通過深度學(xué)習(xí)的方法往往是數(shù)據(jù)驅(qū)動的,需要對大量的數(shù)據(jù)進(jìn)行訓(xùn)練,無法廣泛應(yīng)用到植物建模等領(lǐng)域中。針對上述研究現(xiàn)狀,本研究擬解決的關(guān)鍵問題是:①構(gòu)造植物點云的拓?fù)溥B接情況,判斷缺失部位連接點;②確定待擬合點的點集,并保證斷連部分銜接平滑且具有真實性;③對缺失部分進(jìn)行擬合并調(diào)整點云密度,使補全結(jié)果與原始點云渾然一體。通過解決這些問題,以期能夠在一定程度上彌補掃描時外界因素所帶來的缺陷,使植物模型更為生動形象且有效地應(yīng)用于教學(xué)、科研等方面,也更為直觀地展現(xiàn)植物的真實形態(tài),為數(shù)字林業(yè)提供依據(jù),助力林業(yè)在信息化時代的發(fā)展。

    1 材料與方法

    1.1 植物點云數(shù)據(jù)選取

    使用三維掃描儀對植物進(jìn)行掃描,受植物形態(tài)結(jié)構(gòu)、遮擋問題和掃描精度的影響,掃描時往往會丟失部分點云數(shù)據(jù),其中植物的葉片與其莖干的連接處通常較細(xì)容易出現(xiàn)缺失情況,如圖1A所示。為了進(jìn)行定量指標(biāo)驗證,本研究使用一些完整的植物CAD模型,通過虛擬掃描技術(shù)獲取其三維點云,去除莖干連接處部分點云后進(jìn)行模擬,如圖1B所示。

    1.2 植物缺失部分定位

    1.2.1 確定缺失部位連接點

    不同植物經(jīng)過掃描后其缺失的點云部分不盡相同,因此需要先定位出缺失的點云部分。對于輸入的植物點云,首先對其坐標(biāo)進(jìn)行歸一化,即將點云坐標(biāo)映射到(-1,1)區(qū)間內(nèi)。然后使用K近鄰算法得到點云中每個點的最近鄰點,計算兩點間的距離,最終可由所有點與其最近鄰點間距離的平均值,得到點云密度。點云密度計算公式:

    (1)

    其中:Dplant表示點云密度,Cplant為點云集合,‖·‖2表示任意兩點間歐氏距離,NCplant表示點云集合中點的數(shù)量。得到點云密度后,可根據(jù)點云密度判斷點云中任意兩點間的距離并對點云進(jìn)行聚類,聚類標(biāo)準(zhǔn)為:

    (2)

    δ=σ·Dplant。

    (3)

    其中:dP,Q為P、Q兩點間距離,δ表示用于聚類分類的距離標(biāo)準(zhǔn),σ表示閾值(默認(rèn)為3.0)。若dP,Q≤δ,則P、Q兩點屬于同一簇;若dP,Q>δ,則P、Q兩點不屬于同一簇。

    使用概率圖模型確定每個點云簇間的拓?fù)溥B接情況,同時確定每個點云簇的連接點。將每個點云簇都視為一個圖結(jié)構(gòu)中的節(jié)點,每兩個節(jié)點之間邊的權(quán)重定義為:

    (4)

    其中:C1、C2是兩個點云簇,WC1,C2為兩個節(jié)點間邊的權(quán)重,P1是集合C1中的點,P2是集合C2中的點,NC1、NC2分別是集合C1、C2中點的數(shù)量。

    如果兩節(jié)點間邊的權(quán)重越小,那么這兩個節(jié)點越有可能相連接;反之,如果兩節(jié)點間邊的權(quán)重越大,那么這兩個節(jié)點越不可能相連接。通常而言,每兩個節(jié)點間邊權(quán)重的定義應(yīng)當(dāng)僅與這兩個節(jié)點對應(yīng)的點云簇中點的最近距離有關(guān)[17],即:

    (5)

    然而植物的葉片、花朵等部位不一定與莖干距離最近,存在葉片與葉片距離最近或葉片與花朵距離最近等情況,故本研究采取最近距離與點的數(shù)量相結(jié)合的權(quán)重定義方式。在概率圖模型中使用Prim算法構(gòu)建最小生成樹,使得生成的樹形結(jié)構(gòu)邊的權(quán)重和最小,即得到概率圖中最大概率的拓?fù)溥B接情況。

    構(gòu)建最小生成樹后,可依據(jù)最小生成樹中節(jié)點的度判定莖干點云。Huang等[15]對點云使用基于L1中值的局部迭代方法,可得到鄰域內(nèi)點的局部中心,并逐漸擴大鄰域范圍,最終獲得點云大致的骨架。然而該方法中需要調(diào)節(jié)的參數(shù)較多,整體流程上不夠魯棒,故本研究僅對莖干點云使用L1中值局部迭代方法提取莖干部分的骨架。對于提取出的莖干骨架點集,整體連線近似于莖干的軸線,點集中的每個點則近似于莖干的中心點,為后續(xù)確定莖干上的連接點及使用Bezier曲線對缺失部分進(jìn)行擬合等步驟奠定基礎(chǔ)。使用該方法提取得到的莖干點云的骨架點集如圖2所示。

    由于點云具有無序性,對后續(xù)步驟的實施造成影響,因此需要對莖干骨架點集進(jìn)行排序。本研究提出一種基于最近鄰點距離的搜索算法,首先迭代最近鄰點判斷輸入點集的端點,進(jìn)而從端點出發(fā)尋找當(dāng)前點的最近鄰點即可對莖干骨架點集進(jìn)行排序,且能滿足旋轉(zhuǎn)不變性。該算法具有較高的魯棒性,能夠應(yīng)對不同順序、不同稀疏程度的骨架點集。通過該算法尋找點集端點并對點集進(jìn)行排序的流程圖如圖3所示。

    對任意輸入點集隨機選取其中一點尋找點集端點的過程如圖4B所示,圖中點顏色的漸變(從紅色到白色)為搜索算法尋找過程,最終對點集排序的結(jié)果如圖4C所示,圖中顏色的漸變?yōu)辄c排序順序的結(jié)果。

    本研究需要對葉片點云及莖干點云間缺失的部分進(jìn)行補全,找到葉片點云及莖干點云中的連接點。莖干的擬合及連接應(yīng)當(dāng)從莖干中心出發(fā),即莖干部分的連接點應(yīng)當(dāng)在莖干骨架點集上。本研究利用葉片點云重心與葉片連接點的連線同莖干軸線的交點作為連接點。對于植物葉片而言,讓連接莖干指向葉片重心在現(xiàn)實中能起到對葉片的支撐作用,因此按該方法定義的莖干骨架連接點在物理層面更具有實際意義。對葉片點云重心的求取是計算葉片點云中所有點坐標(biāo)的平均值。

    但實際上三維空間中兩條直線不一定存在交點,因此對于葉片重心P與葉片連接點Q的連線PQ與莖干軸線不一定存在交點。故在實際算法執(zhí)行中要遍歷莖干骨架點集,尋找點集中距離連線PQ最近的點R,則點R為莖干骨架點集連接點。本研究使用海倫公式計算三維空間中點到直線的距離(d):

    (6)

    (7)

    (8)

    式中,p為半周長。對莖干骨架點集中的每個點R,計算其到連線PQ的距離d,在莖干骨架點集中到連線距離最小的點即為莖干骨架點集連接點。

    1.2.2 確定待擬合點的點集

    確定缺失部位連接點后可得到莖干骨架連接點,但仍需確定莖干骨架點集上參與擬合的待擬合點。根據(jù)植物的生長特征,莖干是由底部向上生長,因此需在莖干連接點下方的點中選取待擬合點。對于已經(jīng)過排序的莖干骨架點集,仍無法確定其順序的方向,在保證旋轉(zhuǎn)不變性的前提下,仍需對莖干骨架的排序順序進(jìn)行調(diào)整,確保其在莖干軸線上從下到上排序的順序。通常對于植物的形態(tài)而言,植物的上半部分存在葉片、花朵等部位較為豐富而下半部分僅存在莖干較為稀疏的現(xiàn)象,因此可通過植物點云中點的密集程度來判斷植物在形態(tài)意義上的上下位置關(guān)系。

    為了近似將植物點云從中線劃分成上下兩個部分,取莖干骨架點集的中點,及當(dāng)前點集順序下一點,以這兩點連線的方向向量(指向莖干骨架點集的中點)表示莖干軸線方向向量,過莖干骨架點集的中點作莖干的橫截面,以該截面為界分為兩個部分的點集,如圖5所示。通過兩部分點集內(nèi)點的數(shù)量來確定上半部分點集與下半部分點集,若上半部分點集內(nèi)的點與莖干的中點連線的方向向量(指向莖干骨架點集的中點)與莖干軸線方向向量相同,則表明當(dāng)前莖干軸線方向與植物形態(tài)方向相同;反之若兩向量方向相反,則表明當(dāng)前莖干軸線方向與植物形態(tài)方向相反,則需將莖干骨架點集順序反轉(zhuǎn),從而將莖干骨架點集的方向與植物形態(tài)意義上的方向調(diào)整一致。

    本研究使用球鄰域搜索選取待擬合點,保證擬合時莖干骨架上參與擬合的點與缺失部分需要擬合生成的點所占比例的權(quán)重近似相等。以莖干骨架點集連接點為球心,以該點與葉片點云連接點間的距離為半徑,該球范圍內(nèi)在莖干骨架點集中的點都將參與擬合。通過這種方法對缺失部分?jǐn)M合得到的莖干半徑與莖干軸線將更為平滑,具有較強的真實感。

    接著,需要確定缺失部分待擬合點,即葉片點云與莖干骨架間的待擬合點。為保證擬合效果更為平滑,葉片點云連接點到莖干骨架點集連接點的距離與莖干骨架點集連接點到上述選擇的待擬合點點集中最下方點的距離應(yīng)相當(dāng)。故依據(jù)莖干骨架點集上的待擬合點的數(shù)量,在葉片點云連接點與莖干骨架點集連接點的連線段上均勻取點,這些點將作為缺失部分的待擬合點參與莖干半徑的擬合。

    1.3 植物缺失部位擬合及補全

    1.3.1 擬合缺失莖干半徑

    缺失莖干往往在現(xiàn)實植物中較細(xì),在掃描結(jié)果中是需要補全的部分。對于該部分莖干,靠近葉片處莖的半徑應(yīng)當(dāng)較小、較貼合葉片的厚度,靠近主軸莖干部分莖半徑應(yīng)當(dāng)較大并且較貼合主軸莖干的粗細(xì)。首先定義葉片點云連接處莖干半徑,從植物主軸莖干連向葉片部分莖干的粗細(xì)應(yīng)當(dāng)是逐漸減小的,對于該部分半徑的擬合應(yīng)由主軸莖干上的半徑(即莖干骨架點集對應(yīng)橫截面擬合圓的半徑)與葉片點云連接處半徑共同決定。

    對于主軸莖干上的半徑,要得到莖干橫截面上近似的輪廓線,而輪廓線則是由莖干點云在莖干橫截面上投影得到的一系列點。在植物原始點云中,以莖干骨架點集中每個點為中心,使用K近鄰(默認(rèn)K為30,即周圍30個點參與擬合)算法計算得到與之相距最近的K個點。以莖干骨架點集中的點P0(x0,y0,z0)及在該點處的方向向量n0(a0,b0,c0)為橫截面的法向量,由三維空間平面點法式得到該橫截面S的方程為:

    a0(x-x0)+b0(y-y0)+c0(z-z0)=0。

    (9)

    由K近鄰算法得到的每個點Pk(xk,yk,zk),過該點且垂直于橫截面直線的方向向量即為橫截面法向量n0,故該直線L的方程為:

    (10)

    由式(9)與式(10)可得K近鄰點集中每個點在莖干橫截面上的投影點P′k的坐標(biāo)為:

    (11)

    由式(11)可以看出,對于投影點P′k(x′k,y′k,z′k)的坐標(biāo)公式具有相似性,且該公式對于三維空間內(nèi)任意一點同樣適用,因此對于點在莖干橫截面上投影的求取算法具有旋轉(zhuǎn)不變性的特征。

    根據(jù)植物的特征,植物的莖干在局部范圍內(nèi)大多呈圓柱體狀,因此以莖干軸線作莖的橫截面,截得植物莖干點云形成的輪廓線近似可看作圓,故需要在橫截面內(nèi)以上述K近鄰點的投影點擬合圓。在三維空間內(nèi)對圓的擬合,需滿足待擬合點點集中所有點在同一平面內(nèi),先擬合球體,再聯(lián)合球與平面的方程得到三維空間內(nèi)圓的方程。使用最小二乘法對點擬合三維空間內(nèi)球體:

    (12)

    其中:Pb(xb,yb,zb)為擬合球體的球心,Q為點P在K近鄰范圍內(nèi)點集中的點,rb為擬合球體的半徑。由式(11)可得點Pb在莖干橫截面上的投影點P′(x′0,y′0,z′0),該點將作為迭代更新優(yōu)化后新的莖干骨架點集中的點,由此可計算擬合球體與莖干橫截面所在平面相交圓,即莖干在點P處擬合圓的半徑rP:

    (13)

    式中,h為擬合球心到莖干橫截面的距離。

    根據(jù)最小二乘法擬合得到的圓如圖6所示。

    對于葉片點云連接處莖干的直徑,應(yīng)當(dāng)近似為葉片點云的厚度。由于點云的無序性且葉片存在彎折情況,因此無法根據(jù)葉片點云的邊界來判定葉片的厚度??紤]到葉片點云與葉片點云連接處周圍莖干部分點云的密度應(yīng)當(dāng)近似,即兩部分在連接處應(yīng)當(dāng)分布均勻且較為平滑,故求取葉片點云中所有點K近鄰(經(jīng)過實驗分析設(shè)定K為3)范圍內(nèi)點間距離平均值的平均值,此距離長度則近似為葉片點云連接處莖干的直徑。葉片點云連接處莖干的半徑(rleaf)為:

    (14)

    其中:Cleaf為葉片點云,P為葉片點云中的點,CKNN(P,K)為點P在葉片點云中K近鄰范圍內(nèi)的點集,Q為該點集中的點,NCleaf為葉片點云中點的數(shù)量。

    由于莖干缺失部分的半徑應(yīng)當(dāng)是逐漸減小的,對于缺失部分的起始部位,即靠近主軸莖干部分的莖干半徑應(yīng)當(dāng)貼合于主軸莖干半徑;對于缺失部分的終止部位,即靠近葉片點云連接處部分的莖干半徑應(yīng)當(dāng)貼合于葉片厚度。故對于待擬合半徑集合:

    R={rt|rt∈Rstem,t∈[0,t0]}∩{rt|rt=rleaf,t=1};

    (15)

    (16)

    其中:Rstem為莖干骨架點集中待擬合點對應(yīng)的莖干半徑集合,即Rstem={rP|P∈Cstem_fit};t為Bezier曲線的參數(shù),t0為莖干骨架點集中最后一個待擬合點對應(yīng)的Bezier曲線參數(shù),NCstem_fit為莖干骨架上待擬合點點集的數(shù)量,NCleaf_fit為缺失部分待擬合點點集的數(shù)量。

    進(jìn)而對待擬合半徑集合R與Bezier曲線參數(shù)t使用最小二乘法擬合二次曲線R′(t),根據(jù)待生成點集中點的數(shù)量計算出待生成的Bezier曲線參數(shù)t′,從而計算出待補全部分的莖干半徑集合R′。

    1.3.2 擬合密度均勻的缺失莖干

    參與Bezier曲線擬合的控制點數(shù)越多,即Bezier曲線的次數(shù)越高,擬合生成的曲線段越貼近控制點,且在空間內(nèi)彎折較多。若將待擬合點點集中所有點均作為控制點參與擬合,得到的缺失部分莖干軸線將過于彎曲,真實感欠佳,如圖7A所示。因此,將擬合二次Bezier曲線作為缺失部分莖干的軸線,該曲線為拋物線,得到的擬合軸線將更為平滑,真實感較強。二次Bezier曲線的3個控制點分別為葉片點云連接點、莖干骨架點集連接點和莖干骨架點集中最下方點。由Bezier曲線鎖定端點控制點性質(zhì),Bezier曲線擬合生成的曲線段必過葉片點云連接點,故莖干在缺失部分與葉片部分的銜接較完整且平滑。最終使用二次Bezier曲線擬合生成的缺失部分莖干軸線段如圖7B所示。

    由二次Bezier曲線擬合生成的曲線段作為缺失部分莖干的軸線,以曲線上的點作關(guān)于該曲線的切線為法向量,在三維空間中做一系列圓,在圓周上均勻取點,所有點構(gòu)成的點集即為補全缺失部分莖干所生成的點集。在三維空間直角坐標(biāo)系中圓的參數(shù)方程為:

    (17)

    其中:圓心為C(x0,y0,z0),圓的半徑為r,向量n1=(a,b1,c1)與向量n2=(a2,b2,c2)為圓面上的兩個單位向量,θ為參數(shù),指圓上某點所對應(yīng)的弧度。

    二次Bezier曲線方程為:

    B(t)=(1-t)2P0+2t(1-t)P1+t2P2。

    (18)

    其導(dǎo)數(shù)為:

    B′(t)=-2(1-t)P0+2(1-2t)P1+2tP2。

    (19)

    擬合缺失部分莖干的生成點集的密度應(yīng)當(dāng)與植物點云密度近似,因此需根據(jù)點云密度確定生成點集中點的密度。生成點集的密度由兩部分組成:

    1)Bezier曲線擬合生成曲線段中點的數(shù)量。計算生成曲線段中點的間隔(d=1/N),其將影響生成點集中點在莖干軸線上的橫向距離。由植物點云密度可得植物點云中最近鄰點的平均距離,而此距離應(yīng)當(dāng)與曲線段中點的間隔相同。由式(19)可得B′(t):

    t∈[0,1]。

    (20)

    其中曲線段的長度為:

    t∈[0,1]。

    (21)

    因此,曲線段中點的間隔為:

    (22)

    由此即可確定Bezier曲線擬合時參數(shù)t′的集合,進(jìn)而確定生成點集中莖干軸線上點的間隔。

    2)圓參數(shù)方程得到圓周上點的數(shù)量。圓周上點的數(shù)量將影響生成點集中在莖干軸線上的縱向間隔,根據(jù)點云密度確定圓周上點間距離:

    (23)

    基于上述點間距離可確定生成點集密度,使生成的缺失部分莖干點云能夠與原始植物點云相匹配。

    本研究采用的方法僅需補全缺失部分點集,而使用Bezier曲線擬合得到的曲線段中既包含缺失部分曲線段也包含莖干內(nèi)部曲線段,因此使用上述方法生成的部分圓周上的點會與原始植物點云中莖干部分的點重合,甚至存在于莖干點云內(nèi)部,故需刪除莖干點云上或其內(nèi)部的生成點。

    通過判斷生成點集中每個點到莖干的距離,若該距離小于該處莖干半徑則將該點刪除;反之,則保留。對于求取某個生成點P到莖干的距離,需先得到莖干骨架點集中到點P距離最近的點Q,及莖干骨架點集中下一點R。根據(jù)海倫公式求得點P到莖干線段QR的距離d。

    遍歷所有生成點集,若d≤Rstem(P)則將當(dāng)前點從生成點集中刪去。

    2 結(jié)果與分析

    2.1 缺失部位連接點確定結(jié)果

    在輸入植物點云并經(jīng)過預(yù)處理后,對點云進(jìn)行聚類,確定點云在植物中所屬的不同部位。對不同點云的分類結(jié)果如圖9所示,其中σ表示閾值,用來調(diào)節(jié)點云聚類的精密程度。σ越小說明兩點間的距離較小才能屬于同一簇;反之,σ越大說明兩點間距離較大也能屬于同一簇。本研究設(shè)定σ的默認(rèn)值為3.0。

    隨后使用最小生成樹判斷點云間的拓?fù)溥B接情況,僅基于最近距離和將最近距離與點的數(shù)量相結(jié)合的結(jié)果如圖10所示。

    將葉片點云與莖干點云進(jìn)行連接,需要先找到相連的葉片點云及莖干點云中的連接點。依據(jù)植物的生長特性,葉片往往是以距離莖干最近的部位與莖干相連接。因此,對于葉片點云連接點而言,該連接點應(yīng)當(dāng)是當(dāng)前集合中與相連集合中所有點的最小距離中最小的點。

    在得到葉片部位的連接點后,根據(jù)本文所述方法確定的莖干骨架點集連接點及其對應(yīng)的最終補全后的莖干效果如圖11所示。

    2.2 待擬合點集確定結(jié)果

    待擬合點點集由莖干骨架上的待擬合點及缺失部分的待擬合點構(gòu)成。兩部分待擬合點如圖12所示,其中:紅點表示缺失部分待擬合點與莖干骨架上待擬合點,這些點都將參與缺失部分莖干半徑的擬合;黃點表示葉片點云連接點、莖干骨架點集連接點和莖干骨架上待擬合點中最下方點,這些點將作為Bezier曲線擬合的控制點參與缺失部分莖干軸線的擬合。

    2.3 擬合結(jié)果密度及細(xì)節(jié)調(diào)整

    獲得待擬合點點集后,通過Bezier曲線對該點集進(jìn)行擬合,并基于曲線方程用在其切線方向上作圓擬合缺失部分莖干點云。根據(jù)本研究所述方法,對擬合的生成點集的密度進(jìn)行調(diào)整,結(jié)果如圖13所示。

    最后,刪去原始莖干點云上或其內(nèi)部的生成點,結(jié)果如圖14所示。

    2.4 植物莖干缺失部分補全結(jié)果

    2.4.1 補全效果分析

    為了能對實驗結(jié)果進(jìn)行可視化比較和定量指標(biāo)驗證,使用一些完整的植物CAD模型,去除其莖干連接處部分點云后進(jìn)行模擬,得到的實驗結(jié)果如圖15所示。其中圖15A為輸入的莖干與葉片連接處缺失的待補全點云,圖15B為通過本文方法擬合的結(jié)果,圖15C為用于參照的真實值點云。

    對于補全后的缺失部分莖干,半徑由莖干至葉片銜接平滑,補全的點云與原始點云密度分布均勻,且補全的莖干整體上彎曲自然,具有一定的物理意義與美觀效用。缺失部分補全的莖干與原始點云渾然一體,具有較強的真實感。

    2.4.2 定量指標(biāo)分析

    采用倒角距離和均方誤差兩種定量指標(biāo)對實驗結(jié)果進(jìn)行分析。

    倒角距離(chamfer distance,CD)用來評估兩個點集S1和S2之間距離(其中S1表示點云真實值數(shù)據(jù),S2表示本實驗的結(jié)果)。eCD表示倒角距離誤差,具體評價公式如下:

    (24)

    式中:|S1|、|S2|分別表示S1、S2點集中元素的數(shù)量。采用均方誤差(mean square error,MSE,式中記為eMSE)來評估實驗結(jié)果,公式如下:

    (25)

    式中,CNN(x)表示S1中的點x在S2中的最近鄰點點集,在此使用K=10的最近鄰搜索,即|CNN(x)|=10。

    在兩種評價指標(biāo)中,本實驗的結(jié)果都取得了較好的結(jié)果,如表1所示(數(shù)值越低表示誤差越小)。

    表1 實驗結(jié)果誤差分析表

    2.4.3 魯棒性分析

    本研究方法適用于大多數(shù)存在莖干與葉片斷連情況的植物點云,提出的算法具有旋轉(zhuǎn)不變性,具有較高的魯棒性。

    為了驗證該方法的魯棒性,對輸入點云進(jìn)行隨機角度的旋轉(zhuǎn),實驗結(jié)果如圖16所示,發(fā)現(xiàn)對于旋轉(zhuǎn)過后的點云進(jìn)行補全仍能達(dá)到預(yù)期效果。

    3 討 論

    本研究以莖葉體植物為研究對象,通過三維掃描設(shè)備得到植物點云數(shù)據(jù)。首先基于點云密度使用并查集對點云進(jìn)行聚類,通過概率圖模型計算點云簇間邊的權(quán)重,使用最小生成樹確定點云的拓?fù)溥B接情況;然后對于已聚類的點云簇,針對莖干部分點云使用基于L1中值的局部迭代方法提取其骨架,提出一種基于最近鄰點距離的搜索算法對莖干骨架點集進(jìn)行排序,進(jìn)而根據(jù)球鄰域搜索確定待擬合點集;最后使用最小二乘法擬合莖干半徑,基于Bezier缺陷擬合缺失部分莖干軸線,以莖干切線方向使用三維參數(shù)方程生成補全點集,對點云進(jìn)行補全,為后續(xù)的三維建模提供數(shù)據(jù)基礎(chǔ)和技術(shù)支撐?;贚1中值骨架提取的對植物莖干缺失部分進(jìn)行擬合補全的方法,具有旋轉(zhuǎn)不變性及較高的魯棒性。實驗結(jié)果在保證莖干骨架點集提取的真實性及莖干半徑計算的準(zhǔn)確性的前提下,完成了對缺失部位的補全,并基于點云密度對生成點的點集進(jìn)行調(diào)整,使得補全莖干與植物整體連接恰當(dāng),渾然一體,可為后續(xù)植物的三維重建與植物的可視化奠定基礎(chǔ)。

    在補全植物點云的過程中,提取植物的骨架是關(guān)鍵步驟之一。由于植物的復(fù)雜性和點云掃描的不完整性,直接使用L1中值局部迭代方法提取的植物骨架會出現(xiàn)骨架不完整、偏移過大等問題。筆者先基于點云密度對點云聚類,得到植物葉片、植物莖干等點云簇,然后僅提取植物莖干點云骨架。但對植物完整骨架的提取仍是一個亟待研究的問題,提取植物點云的骨架將會是對植物點云重建、植物點云分類分割等工作的基礎(chǔ),也是植物可視化領(lǐng)域研究的重難點之一。

    深度學(xué)習(xí)相關(guān)方法具有較高的魯棒性及泛化能力,但多數(shù)方法是數(shù)據(jù)驅(qū)動的,需要大量的相關(guān)數(shù)據(jù)進(jìn)行訓(xùn)練。通過本研究提出的方法可簡單、直接地獲得完整且效果良好的植物點云模型,進(jìn)而構(gòu)建植物模型數(shù)據(jù)集。進(jìn)一步可考慮在植物模型數(shù)據(jù)集中,使用深度學(xué)習(xí)相關(guān)方法進(jìn)行訓(xùn)練,對點云的分類及點云間拓?fù)溥B接情況的判斷。如可使用PointNet等深度學(xué)習(xí)網(wǎng)絡(luò),直接對點云進(jìn)行語義分割,將點云進(jìn)行分類的同時判斷該部分點云的種類,同時基于分割的點云種類及點云的空間集合結(jié)構(gòu)可推導(dǎo)出點云間的拓?fù)滏溄忧闆r,此類方法將具有更高的魯棒性。

    此外,由于植物本身結(jié)構(gòu)較為復(fù)雜,且現(xiàn)實生活中的植物種類繁多,筆者僅針對莖葉體植物完成對其葉片和莖干連接處的補全,在今后將考慮實現(xiàn)對更多類型植物點云的補全工作,增加本算法的泛化能力,實現(xiàn)對多樣性植物點云的處理及建模,從而促進(jìn)林業(yè)信息化的發(fā)展。

    猜你喜歡
    莖干連接點軸線
    基于活立木莖干含水量的楊樹生長狀態(tài)評估模型構(gòu)建*
    曲軸線工件劃傷問題改進(jìn)研究
    基于A3航攝儀的小基高比影像連接點精提取技術(shù)研究
    基于潛熱效應(yīng)的活立木凍融檢測傳感器設(shè)計與實驗
    基于干型差異的云南松子代苗木生物量研究
    種子(2019年5期)2019-07-02 11:37:44
    基于回歸分析的水電機組軸線曲折預(yù)判斷分析
    基于彈性厚粘膠層的結(jié)構(gòu)性連接點響應(yīng)建模和預(yù)測
    汽車文摘(2016年6期)2016-12-07 00:23:38
    巖木瓜莖干化學(xué)成分及PTP1B抑制活性研究
    行書章法淺析(十五)書寫應(yīng)把握行軸線
    基于相關(guān)性篩選原理的公共連接點諧波畸變量的分層量化
    電測與儀表(2015年3期)2015-04-09 11:37:22
    欧美日本亚洲视频在线播放| 亚洲av电影不卡..在线观看| 国产欧美日韩精品亚洲av| www.自偷自拍.com| 中文资源天堂在线| 在线看三级毛片| 视频区欧美日本亚洲| 村上凉子中文字幕在线| 9191精品国产免费久久| 青草久久国产| 午夜亚洲福利在线播放| 香蕉av资源在线| 日韩免费av在线播放| 久久久久久九九精品二区国产 | 久久中文看片网| 亚洲精品久久成人aⅴ小说| 神马国产精品三级电影在线观看 | av在线天堂中文字幕| 日本精品一区二区三区蜜桃| 欧美日本亚洲视频在线播放| 我的老师免费观看完整版| 亚洲精品久久国产高清桃花| 午夜福利欧美成人| 波多野结衣高清无吗| 亚洲成av人片在线播放无| 午夜成年电影在线免费观看| 成人av在线播放网站| 日韩高清综合在线| 国产不卡一卡二| 999精品在线视频| 看免费av毛片| 在线看三级毛片| 久久久久九九精品影院| 波多野结衣高清作品| 欧美日韩黄片免| 国产黄片美女视频| 99久久99久久久精品蜜桃| 精品一区二区三区四区五区乱码| 精品免费久久久久久久清纯| 国产一区二区在线av高清观看| 国产主播在线观看一区二区| 国产精品影院久久| 好看av亚洲va欧美ⅴa在| 国产欧美日韩一区二区精品| 久久精品综合一区二区三区| 黄色女人牲交| 2021天堂中文幕一二区在线观| 日韩国内少妇激情av| 99精品在免费线老司机午夜| 国产麻豆成人av免费视频| 久久精品成人免费网站| 一级a爱片免费观看的视频| 亚洲 欧美 日韩 在线 免费| 欧美日韩国产亚洲二区| 观看免费一级毛片| 桃色一区二区三区在线观看| 中亚洲国语对白在线视频| 久久久水蜜桃国产精品网| av在线播放免费不卡| 亚洲成人免费电影在线观看| 国产亚洲精品综合一区在线观看 | 夜夜看夜夜爽夜夜摸| 国内精品一区二区在线观看| 啦啦啦观看免费观看视频高清| 国产真实乱freesex| 欧美成人性av电影在线观看| 久久久精品欧美日韩精品| av欧美777| 亚洲欧美日韩无卡精品| 身体一侧抽搐| av免费在线观看网站| 女人爽到高潮嗷嗷叫在线视频| 丁香六月欧美| avwww免费| 男人舔奶头视频| 中文亚洲av片在线观看爽| 日韩精品中文字幕看吧| 中文资源天堂在线| 国产亚洲精品av在线| 精品一区二区三区av网在线观看| 桃色一区二区三区在线观看| 男插女下体视频免费在线播放| 十八禁网站免费在线| 99久久精品热视频| 国产在线观看jvid| 午夜福利高清视频| 国产亚洲av高清不卡| 欧美激情久久久久久爽电影| 一区二区三区国产精品乱码| 日韩欧美国产在线观看| 久久草成人影院| 国内精品久久久久久久电影| 欧美成狂野欧美在线观看| 亚洲av电影在线进入| 亚洲,欧美精品.| 在线视频色国产色| 婷婷精品国产亚洲av在线| 成人特级黄色片久久久久久久| 欧美性长视频在线观看| 又黄又爽又免费观看的视频| 亚洲自偷自拍图片 自拍| 色综合婷婷激情| 夜夜躁狠狠躁天天躁| 亚洲人成网站高清观看| 精品国产乱子伦一区二区三区| 最近最新中文字幕大全电影3| 久久久国产成人精品二区| 国产区一区二久久| 久久香蕉激情| 最近最新免费中文字幕在线| 欧美日韩乱码在线| 欧美日韩亚洲综合一区二区三区_| АⅤ资源中文在线天堂| 日韩三级视频一区二区三区| 免费高清视频大片| 床上黄色一级片| 亚洲熟妇熟女久久| 国产成年人精品一区二区| 丰满的人妻完整版| 日本五十路高清| 久久精品人妻少妇| ponron亚洲| 欧美日韩黄片免| 色哟哟哟哟哟哟| 麻豆久久精品国产亚洲av| 一二三四在线观看免费中文在| 在线观看66精品国产| 亚洲一区中文字幕在线| 国语自产精品视频在线第100页| 国产熟女xx| 啦啦啦观看免费观看视频高清| 久久热在线av| 国产精品久久电影中文字幕| 99热这里只有精品一区 | 91国产中文字幕| 亚洲人成电影免费在线| 欧美中文综合在线视频| 国产伦在线观看视频一区| 亚洲自拍偷在线| 国产成年人精品一区二区| 69av精品久久久久久| 久久久久久久精品吃奶| 午夜福利在线观看吧| 男女视频在线观看网站免费 | 亚洲精品一区av在线观看| 长腿黑丝高跟| 久久天躁狠狠躁夜夜2o2o| 日韩大码丰满熟妇| 国产97色在线日韩免费| 国产视频内射| 国产激情久久老熟女| 黑人欧美特级aaaaaa片| www.精华液| 淫妇啪啪啪对白视频| 精品国产乱子伦一区二区三区| 国产精品电影一区二区三区| 99精品在免费线老司机午夜| 丁香六月欧美| 天天一区二区日本电影三级| 欧美日韩精品网址| 亚洲色图av天堂| 亚洲人成网站在线播放欧美日韩| 在线a可以看的网站| 国产精品一区二区精品视频观看| 777久久人妻少妇嫩草av网站| 亚洲欧美日韩东京热| 在线国产一区二区在线| 国产亚洲精品综合一区在线观看 | 久久婷婷人人爽人人干人人爱| 国内精品一区二区在线观看| 在线观看免费视频日本深夜| 啦啦啦免费观看视频1| 狠狠狠狠99中文字幕| 亚洲成a人片在线一区二区| 最近最新中文字幕大全电影3| 在线观看www视频免费| 他把我摸到了高潮在线观看| 嫁个100分男人电影在线观看| 淫妇啪啪啪对白视频| 亚洲成人久久性| 国产av不卡久久| x7x7x7水蜜桃| 国产伦一二天堂av在线观看| 亚洲国产精品999在线| 国产一区二区三区在线臀色熟女| 99精品久久久久人妻精品| 国产久久久一区二区三区| 亚洲第一电影网av| 91老司机精品| 在线国产一区二区在线| 国产亚洲精品综合一区在线观看 | 长腿黑丝高跟| 亚洲无线在线观看| 欧美日韩中文字幕国产精品一区二区三区| 国产黄a三级三级三级人| av在线播放免费不卡| 国产精品影院久久| 亚洲一区二区三区色噜噜| 免费在线观看影片大全网站| 一进一出抽搐gif免费好疼| 亚洲乱码一区二区免费版| 久久精品国产亚洲av香蕉五月| 亚洲,欧美精品.| 久久精品综合一区二区三区| 一个人免费在线观看电影 | 丝袜美腿诱惑在线| av视频在线观看入口| 日韩大码丰满熟妇| x7x7x7水蜜桃| 精品日产1卡2卡| 欧美大码av| 看免费av毛片| 搡老熟女国产l中国老女人| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲七黄色美女视频| 久久久久国内视频| 精品熟女少妇八av免费久了| 国产精品亚洲av一区麻豆| 日本三级黄在线观看| 国产在线精品亚洲第一网站| 看片在线看免费视频| 日日夜夜操网爽| 成人永久免费在线观看视频| 欧美激情久久久久久爽电影| 亚洲中文字幕一区二区三区有码在线看 | 少妇的丰满在线观看| www.熟女人妻精品国产| 亚洲午夜精品一区,二区,三区| 日韩成人在线观看一区二区三区| 亚洲18禁久久av| 1024视频免费在线观看| 少妇的丰满在线观看| 日韩欧美在线二视频| 一级毛片高清免费大全| 久久国产精品影院| 亚洲乱码一区二区免费版| 国产精品一及| 亚洲九九香蕉| 国产av又大| 久久草成人影院| 国产精品一区二区免费欧美| 一本大道久久a久久精品| 午夜两性在线视频| 午夜影院日韩av| 精品久久久久久久毛片微露脸| www.www免费av| 可以免费在线观看a视频的电影网站| 精品第一国产精品| 日韩精品中文字幕看吧| 三级国产精品欧美在线观看 | 国产精品亚洲一级av第二区| 神马国产精品三级电影在线观看 | 老司机午夜福利在线观看视频| 国产久久久一区二区三区| 黄色毛片三级朝国网站| 99久久99久久久精品蜜桃| 成人特级黄色片久久久久久久| 日韩精品免费视频一区二区三区| 久久久久久大精品| a在线观看视频网站| 狂野欧美激情性xxxx| 中文资源天堂在线| 日韩中文字幕欧美一区二区| 日韩欧美精品v在线| 久久久久久久久久黄片| 一夜夜www| 18禁黄网站禁片午夜丰满| 天堂影院成人在线观看| 一边摸一边做爽爽视频免费| 99精品欧美一区二区三区四区| 久久久精品国产亚洲av高清涩受| 啦啦啦观看免费观看视频高清| ponron亚洲| 国产欧美日韩精品亚洲av| 国产亚洲精品久久久久5区| 制服诱惑二区| 成年女人毛片免费观看观看9| 两人在一起打扑克的视频| 国产v大片淫在线免费观看| 欧美最黄视频在线播放免费| 曰老女人黄片| 黄色 视频免费看| 丰满的人妻完整版| 成人午夜高清在线视频| 男插女下体视频免费在线播放| 99久久综合精品五月天人人| 搞女人的毛片| 亚洲av成人不卡在线观看播放网| 88av欧美| 一级片免费观看大全| 免费在线观看亚洲国产| 成人一区二区视频在线观看| 亚洲国产高清在线一区二区三| 99在线人妻在线中文字幕| 99热只有精品国产| 69av精品久久久久久| 国内精品久久久久久久电影| 久久久久国产一级毛片高清牌| svipshipincom国产片| 麻豆一二三区av精品| 久热爱精品视频在线9| 亚洲乱码一区二区免费版| 天堂√8在线中文| 麻豆av在线久日| √禁漫天堂资源中文www| 老司机深夜福利视频在线观看| 麻豆久久精品国产亚洲av| 久久99热这里只有精品18| 亚洲五月天丁香| 狂野欧美激情性xxxx| 成年女人毛片免费观看观看9| 99久久综合精品五月天人人| 免费观看精品视频网站| 久久精品91蜜桃| 亚洲国产精品sss在线观看| 搡老妇女老女人老熟妇| tocl精华| 午夜成年电影在线免费观看| 亚洲成av人片在线播放无| 欧美日韩黄片免| 91麻豆精品激情在线观看国产| 亚洲一码二码三码区别大吗| 免费在线观看成人毛片| 日韩大尺度精品在线看网址| 小说图片视频综合网站| 少妇人妻一区二区三区视频| 亚洲无线在线观看| 午夜精品久久久久久毛片777| 熟女少妇亚洲综合色aaa.| 99久久99久久久精品蜜桃| 国产高清视频在线播放一区| 久久99热这里只有精品18| 亚洲av熟女| 两人在一起打扑克的视频| 亚洲国产日韩欧美精品在线观看 | 国产精品久久久人人做人人爽| 国产精品av久久久久免费| 麻豆成人午夜福利视频| 啦啦啦观看免费观看视频高清| 一级黄色大片毛片| 国产又色又爽无遮挡免费看| 久久人妻福利社区极品人妻图片| 在线观看免费午夜福利视频| 国产成人av激情在线播放| av福利片在线观看| 18禁黄网站禁片免费观看直播| 床上黄色一级片| 午夜精品在线福利| 一本精品99久久精品77| 国产精品1区2区在线观看.| 亚洲国产看品久久| 国产亚洲精品久久久久5区| 在线看三级毛片| 在线国产一区二区在线| 日韩免费av在线播放| 国产黄a三级三级三级人| a在线观看视频网站| 伦理电影免费视频| 在线观看66精品国产| av欧美777| 国产成人精品无人区| 欧美丝袜亚洲另类 | 国产精品香港三级国产av潘金莲| 在线观看美女被高潮喷水网站 | 女生性感内裤真人,穿戴方法视频| 男插女下体视频免费在线播放| 国产成人啪精品午夜网站| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美在线黄色| 18禁美女被吸乳视频| 国产v大片淫在线免费观看| 国产免费男女视频| 国产三级在线视频| 99国产综合亚洲精品| 三级毛片av免费| 真人做人爱边吃奶动态| 国产精品综合久久久久久久免费| 亚洲18禁久久av| 两性夫妻黄色片| 50天的宝宝边吃奶边哭怎么回事| 麻豆久久精品国产亚洲av| 黄片小视频在线播放| 男人的好看免费观看在线视频 | 巨乳人妻的诱惑在线观看| 极品教师在线免费播放| 国产97色在线日韩免费| 国产精品国产高清国产av| 久久这里只有精品中国| 久久久久久久精品吃奶| 久久久国产成人精品二区| 日本 欧美在线| 亚洲黑人精品在线| 国产成人av激情在线播放| 精品电影一区二区在线| 免费看美女性在线毛片视频| 精品久久蜜臀av无| 一本大道久久a久久精品| 亚洲自拍偷在线| 一夜夜www| 国产精品久久电影中文字幕| 亚洲欧美日韩东京热| 中文字幕人成人乱码亚洲影| 日韩av在线大香蕉| 中文在线观看免费www的网站 | 亚洲一码二码三码区别大吗| 国产精品久久电影中文字幕| 亚洲精品在线美女| 天堂动漫精品| 最近在线观看免费完整版| 久久 成人 亚洲| 国产爱豆传媒在线观看 | 亚洲人成网站高清观看| 黄色 视频免费看| a级毛片a级免费在线| 首页视频小说图片口味搜索| 757午夜福利合集在线观看| 色综合亚洲欧美另类图片| 国产真实乱freesex| 麻豆一二三区av精品| 97人妻精品一区二区三区麻豆| 日本 欧美在线| 天堂av国产一区二区熟女人妻 | 深夜精品福利| 午夜精品在线福利| 国产片内射在线| 国产高清视频在线播放一区| 日本a在线网址| 九色国产91popny在线| 色哟哟哟哟哟哟| 日本免费a在线| 欧美国产日韩亚洲一区| 老汉色∧v一级毛片| 天堂√8在线中文| 久久久久久亚洲精品国产蜜桃av| 又紧又爽又黄一区二区| 亚洲国产欧洲综合997久久,| 国产v大片淫在线免费观看| 亚洲专区国产一区二区| 长腿黑丝高跟| 成年女人毛片免费观看观看9| 国内精品一区二区在线观看| 亚洲18禁久久av| 成人国语在线视频| 久久国产精品影院| 国产一区二区在线av高清观看| 精品日产1卡2卡| av福利片在线观看| 国产精品av视频在线免费观看| 亚洲,欧美精品.| 国产成年人精品一区二区| 狂野欧美激情性xxxx| 999精品在线视频| 在线观看美女被高潮喷水网站 | 欧美一区二区国产精品久久精品 | 国产精品电影一区二区三区| 此物有八面人人有两片| 岛国在线免费视频观看| 国产又黄又爽又无遮挡在线| 午夜福利高清视频| 两个人视频免费观看高清| 91大片在线观看| 老汉色∧v一级毛片| 免费观看人在逋| 狂野欧美激情性xxxx| 小说图片视频综合网站| 亚洲av电影在线进入| 国产精品美女特级片免费视频播放器 | 免费看美女性在线毛片视频| 日本一本二区三区精品| 97人妻精品一区二区三区麻豆| 黄片大片在线免费观看| 丁香六月欧美| 午夜福利在线在线| 丝袜美腿诱惑在线| 国产片内射在线| 日韩欧美 国产精品| 国产高清有码在线观看视频 | 国产午夜精品论理片| 国产精品久久久人人做人人爽| 欧美色欧美亚洲另类二区| 亚洲国产日韩欧美精品在线观看 | 少妇的丰满在线观看| 又紧又爽又黄一区二区| 亚洲,欧美精品.| 女警被强在线播放| 老司机深夜福利视频在线观看| 亚洲国产欧洲综合997久久,| 又爽又黄无遮挡网站| 999久久久国产精品视频| 色播亚洲综合网| 99久久99久久久精品蜜桃| 亚洲片人在线观看| 香蕉国产在线看| 欧美日韩国产亚洲二区| 国产精品一区二区三区四区久久| 日本撒尿小便嘘嘘汇集6| 小说图片视频综合网站| 中文字幕熟女人妻在线| 久久午夜综合久久蜜桃| 美女免费视频网站| 久久久久国内视频| 人妻丰满熟妇av一区二区三区| 久久精品夜夜夜夜夜久久蜜豆 | 欧美一级毛片孕妇| 欧美日韩黄片免| 中文字幕av在线有码专区| 国产精品香港三级国产av潘金莲| 亚洲成av人片免费观看| 超碰成人久久| 少妇被粗大的猛进出69影院| 午夜精品久久久久久毛片777| 欧美日韩精品网址| 久久香蕉国产精品| 国内久久婷婷六月综合欲色啪| 亚洲精品在线观看二区| 日本一二三区视频观看| 亚洲va日本ⅴa欧美va伊人久久| АⅤ资源中文在线天堂| 国产av一区在线观看免费| 村上凉子中文字幕在线| 黄片大片在线免费观看| 日韩国内少妇激情av| 变态另类丝袜制服| 亚洲第一欧美日韩一区二区三区| 日韩中文字幕欧美一区二区| 香蕉丝袜av| 国产激情欧美一区二区| 久久99热这里只有精品18| 国产视频内射| 免费在线观看影片大全网站| 啪啪无遮挡十八禁网站| 人人妻人人看人人澡| 欧美午夜高清在线| 人妻久久中文字幕网| 国产三级在线视频| 热99re8久久精品国产| 精品不卡国产一区二区三区| 天天添夜夜摸| 国产av一区在线观看免费| 亚洲美女黄片视频| 禁无遮挡网站| 99国产精品99久久久久| av欧美777| 老司机午夜十八禁免费视频| 一区二区三区国产精品乱码| 久久精品夜夜夜夜夜久久蜜豆 | svipshipincom国产片| 色在线成人网| 非洲黑人性xxxx精品又粗又长| 在线观看66精品国产| 搡老妇女老女人老熟妇| 日韩大尺度精品在线看网址| 国产不卡一卡二| x7x7x7水蜜桃| 日本黄色视频三级网站网址| 亚洲一码二码三码区别大吗| 国产熟女xx| 亚洲熟妇熟女久久| 九色国产91popny在线| 草草在线视频免费看| xxxwww97欧美| 每晚都被弄得嗷嗷叫到高潮| videosex国产| 国产99久久九九免费精品| 1024香蕉在线观看| 亚洲精品国产一区二区精华液| 天天躁夜夜躁狠狠躁躁| 成人手机av| 小说图片视频综合网站| 岛国在线免费视频观看| 欧美绝顶高潮抽搐喷水| 久久久久久国产a免费观看| 久久人妻av系列| 在线观看一区二区三区| 老司机福利观看| 久久亚洲精品不卡| 黄色丝袜av网址大全| 99re在线观看精品视频| 首页视频小说图片口味搜索| 波多野结衣高清无吗| 欧美黄色淫秽网站| 精品人妻1区二区| 亚洲欧美激情综合另类| 伊人久久大香线蕉亚洲五| 99在线视频只有这里精品首页| 国产视频一区二区在线看| 大型av网站在线播放| 久久久久国产精品人妻aⅴ院| 俄罗斯特黄特色一大片| 亚洲成av人片在线播放无| 久久中文字幕人妻熟女| 成人精品一区二区免费| 一个人免费在线观看电影 | 欧美久久黑人一区二区| av在线播放免费不卡| 校园春色视频在线观看| 亚洲欧美日韩无卡精品| 久久久久久亚洲精品国产蜜桃av| 亚洲专区中文字幕在线| 午夜激情福利司机影院| 欧美+亚洲+日韩+国产| 日韩有码中文字幕| 人妻夜夜爽99麻豆av| 真人一进一出gif抽搐免费| 妹子高潮喷水视频| 久久香蕉激情| 午夜福利高清视频| 国产成人av教育| 999久久久国产精品视频| 久9热在线精品视频| 好男人电影高清在线观看| 亚洲人与动物交配视频| 亚洲真实伦在线观看| av片东京热男人的天堂| 国产亚洲精品第一综合不卡|