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

    距離約束的網格曲面曲線設計方法?

    2020-01-02 03:45:58俞成海馬文娟何利力
    軟件學報 2020年10期
    關鍵詞:光順流形插值

    金 耀,宋 丹 ,俞成海,馬文娟 ,宋 瀅 ,何利力

    1(浙江理工大學 信息學院,浙江 杭州 310018)

    2(天津大學 電氣自動化與信息工程學院,天津 300072)

    在網格曲面上進行曲線設計是進行幾何處理與編輯的一項重要內容,且有著廣泛的應用背景,如網格處理中的部件分割[1,2]、虛擬手術中的模擬切割[3,4]、制造業(yè)中的圖案設計[5-7]以及流形上的計算幾何[8-11]等.這種曲線有多種設計方式,其中常見的有插值與光順兩種:用戶在網格曲面上輸入幾個控制點或一條初始曲線,系統便生成一條附著于曲面的光滑曲線,使之經過這些輸入點(插值)或逼近于初始曲線(光順).換言之,在網格曲面上設計曲線,需同時滿足流形約束(曲線位于曲面上)、光滑約束、插值點約束(若有)等,而其中的流形約束使其較之于歐式空間的曲線設計問題更為復雜,難度更大.

    現有針對網格曲面的曲線設計方法大致可以分為3 類:投影法[12-14]、光順法[1,2,4]以及參數化法[6,15].投影法通常將光滑的空間曲線迭代地投影到網格流形上,方法簡單、高效,但其迭代過程緩慢且投影步驟往往不太魯棒.光順法即對給定曲線在保持流形約束下進行光順,方法通常較為魯棒,但其效率較低.參數化法即將局部曲面參數化到歐式平面,在平面上設計曲線并將結果映射回曲面,方法往往魯棒、高效且靈活,但其受限于局部區(qū)域,難以進行大范圍的曲線設計.

    針對上述不足,本文提出了一種高效且實用的曲線設計方法.該方法的思想是:將流形約束轉化為距離約束,用切平面逼近局部曲面,將距離約束松弛為點到切平面的距離,并與光滑約束、插值點約束共同作為優(yōu)化目標進行求解.求解時采用“整體-局部”交替迭代的思想,并借助Gauss-Newton 法控制其收斂行為.該方法不僅提高了魯棒性與計算效率,而且具有形狀可控、適用范圍更廣等特點,能夠應用于多種場合.

    1 相關工作

    曲線設計是計算機圖形學與計算機輔助設計的一個充滿活力且具有較長歷史的研究話題,現有的大部分研究工作著眼于歐氏空間的曲線設計,其相關的理論與算法被深入研究并已趨于成熟[16,17].隨著數字幾何的廣泛應用以及工業(yè)界與學術界需求的增加,在流形空間,尤其是在曲面上進行曲線設計,越來越多地引起人們的關注.但是,針對這類方法的理論與算法未及前者成熟.根據曲面表示形式的不同,這些方法可分為兩大類:針對參數曲面的方法與針對網格曲面的方法.前者通常采用參數化法:由于參數曲面具有一個天然的參數域,因此通??稍趨涤蛏显O計曲線并將結果映射至原曲面空間[18-20].但是絕大多數方法僅考慮具有單個參數域的參數曲面,對于形狀與拓撲復雜的曲面則不適用.由于在實際應用中網格曲面更為廣泛,因此本文僅討論與研究內容相關的針對網格曲面的工作.

    網格曲面上的曲線由于其所在的曲面具有分段線性特點,通常表示為一條分段線性的折線段.其中,具有廣泛應用背景的離散測地線便是其中一種,且一直是圖形學相關領域的研究熱點[9-11].但測地線旨在計算“局部最短”的曲線,而本文所研究的是離散曲面上的“光滑”曲線,因此本文僅考慮網格曲面上自由曲線的設計.現有方法大致分為以下3 類.

    · 投影法

    這類方法將流形約束進行松弛,首先生成經過插值點的光滑曲線,然后將其迭代地投影到流形上.Hofer 等人[12,13]提出了一種基于能量極小化的樣條曲線插值方法,他們將曲線離散成折線,把問題描述成一個帶流形約束的優(yōu)化問題,其中,目標函數是樣條能量與測地能量的組合;并采用投影梯度法進行求解,即將松弛流形約束得到的空間曲線迭代地投影到曲面上.該方法與流形的表示形式無關,可用于諸如隱式曲面、參數曲面與網格曲面等,但其投影步驟并不能保證一定能找到理想的投影點,魯棒性較差.同時,由于投影采用了移動最小二乘法(MLS)局部擬合法,效率較低.Pottmann[14]隨后提出了高階樣條能量,并將該方法進行推廣在曲面上擬合曲線.Wallner 和Pottmann[21]分別提出兩種方法,將歐氏空間中設計細分曲線的線性細分法推廣至流形曲面:第1種提出用“測地平均法”替代傳統的線性平均法,但需計算大量的測地線,較為耗時;第2 種是在外圍空間設計細分曲線,然后將結果投影到曲面上,但是投影算法對于復雜曲面并不魯棒.Morera 等人[22]與Sarlabous 等人[23]均采用測地細分策略,分別提出測地Bézier 曲線與測地圓錐Bézier 曲線的細分曲線生成算法.該方法雖然魯棒,但涉及到測地線的計算,時間復雜度較高.劉斌等人[5]提出了一種測地B 樣條插值方法,通過插值點反求控制點,并運用最近點投影法將其映射到網格曲面上,然后采用插值誤差反向補償法使曲線逐漸逼近曲面.但該方法同樣耗時且投影算法可能存在多解,從而無法處理復雜曲面.Bock 等人[24]為生成曲邊網格,在線性網格邊上采樣點,將其投影到曲面上,并對投影點運用高階樣條擬合,其中,投影使用了光滑的Phong 法向場,但生成曲線并不要求嚴格地位于曲面.Panozzo 等人[25]借助Frechet 均值定義,將歐式空間的加權平均法推廣到流形空間,并用于設計樣條曲線.他們把網格曲面嵌入到高維空間并在該空間進行加權平均,然后采用Phong 法向場將其投影到嵌入網格.該方法無需迭代且所得曲線光滑性好,但其高維嵌入步驟極為耗時.綜上所述,投影法直觀且易于實現,但在計算效率或魯棒性方面表現較差.

    · 光順法

    這類方法通常被用于優(yōu)化特征邊,往往首先松弛光滑約束,然后在光滑能量驅動下對初始曲線進行光順,并保持曲線始終在曲面上.Jung 等人[26]將圖像中的主動輪廓模型用于網格曲面,對分割線進行光順,通過優(yōu)化表示曲線形狀的內部能量和網格特征的外部能量迭代地演化曲線.Benhabiles 等人[27]改進了該方法中的外部能量,對其進行平均和歸一化處理,所得到的曲線能夠更好地與特征邊對齊.Kaplansky 等人[2]采用測地流方程光順封閉曲線,通過修改擴散準則(方程)以適應不同需求.Wu 等人[28]與Zhang 等人[29]基于測地曲率流方程演化曲線形狀,得到了較好的結果,但其效率難以達到交互要求.由于這些方法所生成的曲線均沿著網格邊,雖然對于提取特征邊有益,但難以設計一般的光滑曲線.Ji 等人[1]也采用了主動輪廓模型,通過構造新的蛇形能量,并結合曲線節(jié)點分裂、移動和移除等拓撲操作使曲線能夠穿越網格面片,所構造的曲線更為光滑.Lawonn 等人[4]提出了一種網格域曲線的拉普拉斯光順算子,通過降低曲線的測地曲率以平滑曲線.該方法能夠控制曲線的光滑度以及與初始曲線的貼近度,但是由于采用了局部光順策略,需要較多的迭代次數且無法插值控制點.總而言之,這類方法魯棒性較好,所設計的曲線通常具有良好的光滑性,但是由于在嚴格的流形約束下進行光順,導致其效率低下,且難以用于設計插值曲線.

    · 參數化法

    這類方法運用平面局部參數化方法,在參數域進行曲線設計,并將結果映射到原曲面.Lee 等人[15]提出了“幾何蛇形”法演變曲線形狀:他們對初始曲線所在的局部區(qū)域進行參數化,并在參數域上運用蛇形能量演化曲線,并將結果映射回網格空間.進而基于極小化原則與局部顯著性提出“智能剪刀”法,并用于網格切割[30].朱文明等人[31]采用保角映射構造參數域,并在參數域上生成細分曲線,最終將結果映射回網格曲面.劉斌等人提出了基于指數映射的測地B 樣條曲線設計方法[6,32],他們將源曲線控制頂點映射到切空間,根據曲線遷移前后控制頂點法保持坐標不變的原則,建立曲線遷移前后控制頂點的對應關系,實現曲線的幾何變換與樣條曲線陣列的構造.但該方法僅限于局部區(qū)域,且需頻繁更新參數域.這類方法可照搬歐式空間的曲線設計方法,簡單且靈活,然而現有方法大多考慮局部參數化方法,難以設計大范圍的曲線,且未曾考慮參數化形變誤差對結果的影響.

    本文提出的方法屬于投影法,但它將流形約束松弛為點到切平面的距離約束進行求解,克服了上述3 類方法存在的缺點:相比于傳統的投影法,該方法同時優(yōu)化光滑與流形能量,使曲線能夠逐步光順的同時,穩(wěn)步地貼近曲面,從而使投影步驟的魯棒性和可靠性大為提高;相比于光順法,該方法由于松弛了流形約束,優(yōu)化自由度更大,且計算效率更高;相比于參數化法,該方法能夠處理局部參數化法難以適用的具有復雜拓撲的網格曲面.

    2 曲線設計方法

    2.1 問題描述

    給定一個嵌入于三維歐氏空間R3的光滑曲面S以及一堆有序點陣{pj}?S,要求一條位于該曲面上且插值于{pj}的光滑曲線c∈S.若采用參數方程表示曲線:c(t)→R3,并運用3 次樣條能量度量曲線的光滑度[12,15],則求解曲面上的光滑插值曲線問題可轉化為計算如下位置與流形約束下的優(yōu)化問題:

    其中,x(t)∈S表示位于曲面S的曲線,x″(t)為曲線關于位置參數t的二階導,其幾何意義與曲線的曲率正相關.

    當S為網格曲面時,其上的曲線可表示為分段線性的離散曲線c=〈V,E〉,其中,V={q1,q2,…,qn}?S為曲線上的頂點集,E為形成曲線的邊集.

    ·當曲線為開曲線時,E==1,2,...,n-1}?S;

    ·當其為閉曲線時,E==1,2,...,n}?S(%為取模運算).

    由于曲線是離散的,公式(1)中曲線的二階導數x″(t)并不存在,因此需首先對該公式進行離散化,并采用數值方法進行求解.

    2.2 流形約束轉化

    對于離散網格曲面S,由于其沒有解析表達式,導致公式(1)中的流形約束難以顯式表示,從而給求解該帶約束的優(yōu)化問題帶來較大的困難.為此,本文引入一個距離函數D:R3×S→R,用于度量空間中任意點到流形S的距離,從而將流形約束轉化為距離約束,即等價地表示為曲線上任意點q∈x(t)到曲面S的距離為0:

    而距離函數D則可定義為點到局部曲面的最短距離:

    其中,d為點到點的歐式距離函數,B(q,r)為q的半徑為r的球形鄰域.

    2.3 離散化

    公式(1)中帶約束的優(yōu)化問題可用數值方法求解,但首先需進行離散化.為方便對二階微分表示的3 次樣條能量進行離散化,類似文獻[12],本文采用弦長參數化方法計算采樣密度,即對相鄰兩個插值點所形成的線段pi pi+1,確定其采樣點個數Ni為

    其中,D為離散步長,亦即每條線段的長度(默認情況下設置為網格平均邊長的五分之一),用于控制曲線的采樣密度.由此,該曲線可離散成折線段,其有序采樣點為{qj},其中,=pi(ki為插值點在采樣點中的序號).由于該離散方式使得曲線上的采樣點能近似均勻分布,因此公式(1)的目標函數可近似地離散化表示為均權拉普拉斯能量[33]:

    其中,對于開曲線,M表示折線段中除卻首尾端點的離散點個數;而對于閉曲線,M表示所有離散點個數.于是,公式(1)可改寫為如下離散表達式:

    為求解該問題,本文運用罰函數法,結合距離公式(3),并將公式(5)中硬約束表示的插值點約束與流形約束分別作成軟約束,則公式(5)進一步轉化為如下優(yōu)化問題:

    其中,

    · 參數ω默認取作大的數值(108),使之近似滿足插值約束;

    ·λ由用戶控制,用于平衡光滑性與流形約束的滿足程度.

    2.4 數值求解

    由公式(7.1)、公式(7.2)所定義的優(yōu)化問題可知:每一對qj與j′q都是彼此相關的,其非線性程度較高,因此直接求解難度較大.為此,本文計算相鄰插值點之間的Dijsktra 路徑作為初始曲線,根據公式(4)所得的采樣點數量在曲線上均勻采樣離散點{qj},并設置{q′j}={qj};然后借鑒局部/整體交替迭代的優(yōu)化策略[34]進行求解,即將該兩組優(yōu)化變量進行分離,并把每次迭代分解為整體階段和局部階段,直至前后兩次計算得到的曲線差異很小為止,其計算過程如下.

    整體階段.固定{q′j},求解{qj}.此時即計算能量函數(7.1)的極小值.由于該式是關于{qj}的二次函數,因此可轉化為矩陣方程組AQ=B,其中,A為常系數矩陣.

    局部階段.固定{qj},求解{q′j}.此時對每一個qj,求解問題(7.2),即在qj的球形鄰域與網格曲面的交集中搜索最近點.由于計算球體與網格的交集以及在該交集處的最近點較為耗時,因此可采用拓撲鄰域替代球體幾何鄰域,即將前一次迭代得到的投影點(n為迭代次數)所在三角形面片f0的r-環(huán)拓撲鄰域的面片集合F作為交集.在該交集中,搜索qj到該局部區(qū)域的最近點距離:

    其中,點到三角形的距離d(qj,fi)可采用文獻[35]中介紹的方法.具體地,以f0作為種子點,以廣度優(yōu)先的方式遍歷F,對每個面片計算該點在面片上的法向投影:若投影點在其內部,則將其作為投影點并結束搜索;否則,計算點到面片的距離并更新最近點,直至遍歷完所有面片.當局部區(qū)域非凸時,有可能存在多個最近點,此時可選取使局部拉普拉斯能量,即最小的最近點作為候選點.

    易證,上述交替迭代的優(yōu)化策略可使能量函數(7.1)單調下降,從而收斂.具體地,在每一次迭代中,其整體階段求解能量函數極小值能獲得比當前能量更低的值;而在局部階段,由于固定{qj},能量函數(7.1)的前2 項為常數,第3 項由于D(qj,S)≤||qj-q′j||,因此也能得到更低的值.從而對于每一次迭代,該策略均將使能量值下降.

    但由于q′j由前一階段曲線采樣點計算得到,而非當前qj的投影點,直接用歐式距離||qj-q′j||度量點到曲面的距離并不精確.因此,上述算法雖能收斂,但收斂速度較為緩慢.本文采用該方法在駱駝頭模型上設計一條經過3 個插值點的閉合曲線,如圖1 所示.隨著迭代的進行,曲線變得越來越光滑,但其變化過程緩慢,迭代次數過多.雖然該算法在理論上是收斂的,且其整體階段僅需高效地回代計算更新曲線坐標,但由于其局部階段的投影計算耗時相對較多,若迭代次數過多將嚴重影響計算效率.在本實驗中,達到最大迭代次數(1 000)依然未能收斂.

    Fig.1 Iteration results and its energy graph with point-to-point distance constraints (λ=0.1)圖1 基于點到點距離約束的迭代結果與能量圖(λ=0.1)

    為此,本文借鑒用于曲面配準的迭代最近點(ICP)[36]的算法思想,用對應點處的切平面逼近局部曲面,并采用點到切平面的距離近似度量點到曲面的距離.設q′j處的切平面為Tj,其對應的法向為N(q′j),則qj到曲面S的距離可近似表示如下:

    其示意圖如圖2 所示.

    Fig.2 Illustration of the distance from a point to its tangent plane圖2 點到切平面距離的示意圖

    從而,公式(7.1)可近似轉化為如下優(yōu)化問題:

    相對于點到點的距離約束,點到切平面的距離約束由于約束曲線采樣點在近切平面上移動,使得一方面給予其足夠的自由度進行光順,另一方面又能貼近網格曲面,因此能夠較快地實現曲線的光順.

    類似基于點到點距離約束的數值方法,本文同樣采用“整體-局部”交替迭代的策略求解公式(10).然而,若直接采用該策略,將可能導致曲線在整體階段“過度”光順而偏離曲面,且不能保證算法收斂——雖然能量函數值在整體階段得到下降,但在局部投影階段可能上升.為此,本文借助Gauss-Newton 法[37]的思想,采用局部線搜索控制其收斂行為.利用公式(10)的Hessian 矩陣H(E(q))計算當前曲線頂點的迭代方向,并結合線性回溯搜索法控制迭代步長.即:當能量函數值大于當前值時,采用對步長折半的方式保證能量函數在迭代過程中單調下降,從而使算法能夠穩(wěn)定收斂.由于公式(10)確定的能量函數為平方和形式,其對應的Hessian 矩陣是對稱正定的,因此其每次迭代的“整體”階段均可轉化為凸問題進行求解,無需進行正則化.其迭代求解的每一步通過如下公式更新曲線頂點位置:

    其中,ω∈(0,1]為迭代步長,通過折半線性回溯法確定大小.

    為了提高計算效率,本文分別從整體和局部階段兩方面進行考慮.整體階段需求解一個線性系統,利用公式(11)中Hessian 矩陣的結構在迭代中始終保持不變的性質,通過對矩陣結構進行預分解以加速.局部階段需搜索最近點,本文用切平面法向量與局部鄰域面片的交點進行近似,一旦找到交點則結束,不必遍歷剩余面片.由于交點不一定存在,此時可采用最近點替代.

    由上述算法得到的曲線,其每一段離散的折線將非常接近曲面但往往不嚴格位于曲面上(對位于不同面片的相鄰兩個點形成的線段).為此,本文采用文獻[38]中介紹的投影法將其映射到曲面上:對于每一條未處于網格表面的線段,由一個端點(起始點)向另一個端點(終止點)逐步構造切割平面,與相關的網格邊進行求交,并將交點更新為起始點,逐步往前傳播.如此循環(huán),直至起始點與終止點重疊.具體細節(jié)可參考文獻[38].

    綜上所述,本文提出的基于距離約束的曲線優(yōu)化算法可描述如下.

    算法1.基于距離約束的曲線優(yōu)化算法.

    3 實驗結果與討論

    本文用C++語言實現了上述算法,并在8 核Intel i7-7700U 處理器、8G 內存的Linux 虛擬機環(huán)境下運行.其中,算法所涉及的線性方程組的求解采用Eigen 庫[39]提供的稀疏Cholesky 分解法——LLT 算法.下面,本文通過一系列實驗展示該算法的特點.

    3.1 算法的收斂性實驗

    本文在人臉模型上設計了一條經過3 個插值點的開曲線.該曲線在迭代20 次后收斂,圖3 給出了迭代的中間結果.由圖可見:隨著迭代次數的增加,曲線在曲面上滑動而趨于光滑,最終收斂并達到穩(wěn)定狀態(tài).大量實驗結果表明:本算法在步長控制策略的作用下均能穩(wěn)定地收斂,且在絕大多數情況下均在梯度條件下收斂.

    Fig.3 Sampled curves on the face model generated during the iterations (λ=0.1)圖3 人臉模型上插值曲線的迭代序列(λ=0.1)

    圖4 所示為3 個迭代曲線圖,分別為能量函數值(公式(10))、對應的梯度模長的平方(施加了log 函數)以及平均距離的變化.隨著迭代的進行,能量函數值逐步下降;其梯度值盡管在局部跳動,但總體亦呈現下降趨勢;曲線到網格的平均距離值也較為穩(wěn)定地下降,且當算法收斂時,其平均距離也達到較小值.此時利用投影法將其映射到曲面上,對曲線光滑度的影響較小.

    本文算法較為耗時部分為迭代中局部階段的投影點計算.該算法在前次迭代的投影點所在面片的r-環(huán)拓撲鄰域搜索投影點.為加速起見,本文采用法向與鄰域的交點替代最近點以過濾部分面片方式,并比較了兩種策略(r=2),實驗結果如圖5(a)所示.由圖可見:兩者所設計的曲線形狀較為相似,幾乎重合,但法向交點策略迭代88 次,耗時104ms,而最近點策略迭代101 次,耗時431ms.因此,前者能在不影響結果的前提下有效地提升效率.本文同時比較了不同鄰域大小對結果的影響,如圖5(b)所示(其中的紅點為插值點,下同).由圖可見:1-環(huán)鄰域結果與其他結果有一定差異,而2-環(huán)與3-環(huán)鄰域的結果幾無區(qū)別,三者的迭代次數分別為126、95 與88 次.由此說明,r對算法收斂快慢有一定影響.這是由于r較小時,曲線移動的步長相應減少,因此需要更多的迭代次數到達極小值點;而當r≥2 時,其結果差異較小,雖然能在一定程度上加快收斂,但同時也意味著將增加局部搜索次數.綜合權衡算法效率,本文將r默認設置為2.

    Fig.4 Generated curves during the iteration on the face model圖4 人臉模型上曲線的迭代變化過程

    Fig.5 Comparison results with different strategies of local neighboring search圖5 采用不同局部鄰域搜索策略的結果比較

    3.2 可控參數λ

    在公式(10)中,參數λ用于平衡曲線的光滑度與流形約束的滿足程度.當λ過大時,流形約束能量將占主導.此時,算法行為將與光滑法(例如文獻[4])類似.曲線在迭代時易陷入局部極小值,即曲線得不到充分光順,同時也將增大方程條件數,使得數值性質變差.當λ過小時,光滑約束能量將占主導,此時,算法行為將與傳統的投影法(如文獻[12])相似.雖然曲線得到了充分的光順,但是距離曲面較遠,投影將破壞曲線的光滑度.當λ取值適宜時,算法將兼具兩類方法的優(yōu)勢.圖6 展示了插值點相同但λ取值不同情況下,人臉模型上的插值曲線.

    Fig.6 Interpolatory curves on the face model with different values of λ圖6 λ取不同值時,人臉模型上的插值曲線

    由圖可見:當λ取值為1 和10 時,曲線偏離初始曲線較近,未能得到充分光順;當其取值為0.1 和0.01 時,曲線偏離初始曲線較遠,但光滑度得到了改善.為了定量比較各曲線的光滑度,本文采用表示曲線曲率的離散拉普拉斯能量,即公式(5)進行度量.λ從0.01 變化到10,各曲線的光滑度分別為[0.0200,0.0133,0.0237,0.0247].通過大量實驗發(fā)現,λ=0.1 時獲得的曲線能夠取得較好的光滑度.因此,本文實驗中將參數λ默認取作0.1.

    3.3 比較實驗

    相比于經典的基于能量極小化的投影法[12],本文算法具有眾多優(yōu)點.本文實現了文獻[12]中描述的算法(部分參數原文未給出,其結果可能與原文有差異),并與本文方法進行了比較.

    · 首先,本文方法效率較高.文獻[12]中的方法首先松弛流形約束,然后逐步將其投影到流形上,其中,投影步驟采用了耗時的移動最小二乘法(MLS)估算局部曲率,使得效率大為降低.圖7(a)所示為兩種方法在雕像上的曲線設計結果,其中,本文方法(右圖)耗時117.2ms,文獻[12]中的方法(左圖)耗時281.9ms(迭代22 次收斂);

    · 其次,兩者采用同樣的離散拉普拉斯能量作為優(yōu)化目標,但本文方法得到的結果通常比文獻[12]的結果更為光滑.文獻[12]在逐步投影時并未考慮光滑約束,且最終還需將隱式曲面上的曲線投影到網格上,因此在一定程度上影響了其光滑性;而本文方法中光滑步驟與投影步驟交替進行,且由于曲面上大部分采樣點落在網格面上,很多時候實際投影點即為該點本身,因此能夠有效地保持光滑性.圖7(a)所示結果的光滑性在視覺上的差異并不明顯,為此本文采用公式(5)來度量光滑度并對其進行比較:左圖的光滑度值為2.49-6,右圖的光滑度值為2.76-6;

    · 再次,本文方法更為魯棒.文獻[12]在具有尖銳邊特征邊的模型上設計曲線時,可能因法曲率估算不夠準確使其步長估算不正確,從而導致投影錯誤,使曲線出現較明顯的鋸齒狀(如圖7(b)左所示),而本文方法不受特征邊影響(如圖7(b)右所示).

    Fig.7 Comparison results of the algorithm[12] and our method圖7 文獻[12]算法與本文算法的結果比較

    本文同時比較了文獻[24]中提出的基于Phong 投影的B 樣條設計方法.

    該方法首先運用多維尺度變換法(multidimensional scaling)將網格曲面嵌入到八維歐式空間,然后在該空間設計B 樣條曲線,并運用Phong 投影法將其映射到嵌入網格上,最后將所得曲線映射回原網格.文獻[24]僅給出了B 樣條的正向設計方法,即根據控制點設計曲線.為方便比較,本文實現了B 樣條的逆向算法設計插值曲線(采用夾持端點條件)[40],并與本文方法進行比較,其結果如圖8 所示.由圖8 可見:兩者結果雖然在曲線外形上存在一定差異,但視覺上均較為光滑.就曲線設計步驟而言,文獻[24]由于無需迭代,效率較高,在該實例中(曲線的離散點數目相同)僅耗時1ms;而本文算法耗時6ms.但是該方法需計算點點測地距離與網格的高維嵌入,計算復雜性很高,在實際應用中可能會帶來不便.

    相比于最新的基于局部拉普拉斯光順的曲線設計算法[4],本文算法亦具有較大優(yōu)勢.文獻[4]的算法有一個參數t∈(0,1],用于控制與初始曲線的逼近程度.為使其具有近似插值效果,將t統一取為1,與本文算法進行比較.首先,本文算法可設計插值曲線與閉合曲線,而文獻[4]算法則不能,其比較結果如圖9 與圖10 所示.圖9 比較了兩種方法在瓶子與梨模型上所生成的開曲線.前者算法得到的曲線沒有經過所有控制點,而本文算法則能插值于控制點;圖10 比較了兩種方法在腳與狗頭模型上所生成的閉曲線.前者算法由于采用了局部拉普拉斯光順算法,無法處理閉合曲線,因此所生成的曲線在首尾連接處出現尖角,而本文算法得到的曲線則在視覺上處處光滑.同時,前者算法所生成曲線的采樣點均在網格邊上,無法做到自適應細分,因此其光滑度受到網格分辨率的影響,如圖9(b)、圖10(b)所示,曲線上的“折角”較為明顯.此外,在多數情況下,本文算法效率較高.對于圖9、圖10 中的4 個例子,本文算法的運行時間分別為27.8、25.4、28.4 與44.0(ms),而文獻[4]算法則分別需要123、27、68、14(ms).

    Fig.8 Comparison results of the algorithm[24] and our method on the wings of the gargoyle model圖8 文獻[24]中算法與本文算法在怪獸翅膀模型的結果比較

    Fig.9 Comarison results for open curves of the algorithm[4] and ours圖9 文獻[4]中算法與本文算法所設計的開曲線的比較結果

    Fig.10 Comarison results for close curves of the algorithm[4] and ours圖10 文獻[4]中算法與本文算法所設計的閉曲線的比較結果

    相比于基于參數化的曲線設計方法[6,15],本文算法適用范圍廣泛.本文實現了文獻[15]中提出的基于參數化的幾何蛇形算法,該方法通過優(yōu)化內部能量和外部能量的組合演化給定曲線.其中,內部能量表示為曲線一階導與二階導能量的線性組合,外部能量與網格特征相關;而其二階導能量即為本文所采用的能量函數(見公式(5)),區(qū)別在于該算法在參數域,而本文算法直接在曲面上計算曲線.為在同等條件下進行比較,本文對該算法進行簡化,僅運用二階導能量優(yōu)化曲線,采用最小二乘共形映射(LSCM)[41]在初始曲線的5-環(huán)鄰域計算局部參數域.由于該優(yōu)化方程沒有外部能量,因此可簡化為線性方程組進行高效求解.圖11 給出了具有2 個虧格的“8 字”模型上所設計的曲線的比較結果.由圖11(a)與圖11(b)的結果可以看出,兩者均能在曲面上設計光滑曲線.文獻[15]中算法由于無需迭代,在效率上比本文算法更有優(yōu)勢,其僅耗時9ms,而本文算法耗時85.3ms.但是文獻[15]無法處理圖11(c)所示曲線——其局部鄰域需經過切割才能進行參數化,而在切割線處難以保證曲線的連續(xù)性.而本文算法可在任意虧格的曲面上自由地設計曲線,與模型的拓撲復雜性無關.

    Fig.11 Interpolatory curve on the eight model圖11 “8 字”模型上的插值曲線

    3.4 應 用

    本文方法不僅能夠設計插值曲線,且能光順給定的曲線.該問題即等價于設計逼近于初始曲線c0∈S的光滑曲線.因此僅需對插值曲線模型稍加改造,即將公式(10)中的插值點約束修改成數據擬合項,并以軟約束的形式加入到目標函數,使得曲線在得到光順的同時,盡可能逼近初始曲線,即.其中,為采樣自c0與qj對應的點;α可根據實際需求(逼近程度)確定,默認設置為0.01.圖12(a)所示為鴨子脖子部位的初始分割線,由于其分割邊界由網格邊依次相連而成,其曲線形狀的鋸齒狀較為明顯;而經過本文算法光順后,可得到圖12(b)所示的曲線,視覺效果上變得更為光滑.

    Fig.12 Curve smoothing result on the duck model圖12 鴨子模型分割線的光順結果

    本文曲線設計算法可用于多種場合,圖13 給出了幾個應用實例.其中,圖13(a)展示了在魚身模型上交互勾勒五角星圖案的例子;圖13(b)給出了在牛模型的身體部位切割孔洞的例子,可用于虛擬手術,觀察器官內部環(huán)境;圖13(c)所示為部件切割的例子,可用于交互分割模型部件.

    Fig.13 Applications of our curve design method圖13 曲線設計的應用實例

    本文算法最為耗時的部分為迭代中所涉及的線性方程組的求解與局部投影計算.由于該算法在整體與局部階段均采用了加速策略,因此求解效率相對較高.表1列出了在各個模型上設計曲線的運行時間.由表1可見:算法效率與網格模型的規(guī)模相關度較小,而主要取決于曲線的離散采樣點的數目,且總體效率可滿足交互應用的要求.

    Table 1 Cost time of generating curve on various models表1 各個模型上曲線的生成時間

    4 結 語

    本文提出了一種針對網格曲面的高效且實用的曲線設計方法,該方法將流形約束松弛為點到切平面距離的約束,并將所有約束共同描述成一個優(yōu)化問題進行求解.該方法魯棒、高效,適用范圍廣,可用于任意復雜拓撲(如多虧格)的網格曲面.此外,相比于前沿算法,本文方法不僅能夠設計閉合與開曲線、插值與光順曲線,而且在效率上占有明顯優(yōu)勢,能夠滿足交互響應要求.

    本文方法通過采用局部線搜索策略控制迭代步長,能夠保證算法收斂,但是該策略由于對步長進行了限制,有時未能充分松弛曲線形狀,從而可能陷入局部最優(yōu)解.在今后,我們將探索更好的收斂策略,在保證收斂性的同時,更好地控制曲線形狀.此外,本文算法的復雜度與曲線采樣點的個數相關,對于設計采樣密度較大的曲線,其效率將成為一個瓶頸.因此,這也是未來需要研究的工作.

    猜你喜歡
    光順流形插值
    緊流形上的Schr?dinger算子的譜間隙估計
    迷向表示分為6個不可約直和的旗流形上不變愛因斯坦度量
    平面網格銑削加工光順刀軌快速生成方法
    Nearly Kaehler流形S3×S3上的切觸拉格朗日子流形
    基于Sinc插值與相關譜的縱橫波速度比掃描方法
    HDSHM系統船體型線光順應用經驗
    一種改進FFT多譜線插值諧波分析方法
    基于四項最低旁瓣Nuttall窗的插值FFT諧波分析
    基于面法向量譜變換的網格光順算法
    計算機工程(2015年4期)2015-07-05 08:27:42
    基于多故障流形的旋轉機械故障診斷
    国产亚洲av片在线观看秒播厂| 18禁动态无遮挡网站| av在线app专区| 伊人久久精品亚洲午夜| 如日韩欧美国产精品一区二区三区 | 久热这里只有精品99| 国产在线免费精品| 99精国产麻豆久久婷婷| 亚洲一级一片aⅴ在线观看| 日本爱情动作片www.在线观看| 最近的中文字幕免费完整| 日本91视频免费播放| 人人妻人人澡人人爽人人夜夜| 日产精品乱码卡一卡2卡三| 久久人人爽人人爽人人片va| 最近手机中文字幕大全| 国产亚洲精品第一综合不卡 | 欧美人与性动交α欧美精品济南到 | 亚洲在久久综合| 国产黄色免费在线视频| 亚洲国产精品一区三区| 久热久热在线精品观看| 欧美日本中文国产一区发布| 视频区图区小说| 日本爱情动作片www.在线观看| 啦啦啦啦在线视频资源| www.av在线官网国产| 亚洲精品日韩在线中文字幕| 欧美另类一区| 国产日韩欧美视频二区| 热99久久久久精品小说推荐| 亚洲精品aⅴ在线观看| 18禁裸乳无遮挡动漫免费视频| 99久久综合免费| 中文字幕亚洲精品专区| 亚洲人成77777在线视频| 久久久久精品久久久久真实原创| 啦啦啦视频在线资源免费观看| 另类精品久久| 日韩,欧美,国产一区二区三区| 国产一区二区三区av在线| 99九九线精品视频在线观看视频| 如何舔出高潮| 国产欧美亚洲国产| 日本wwww免费看| 91久久精品国产一区二区成人| 特大巨黑吊av在线直播| 久久久久久久久久人人人人人人| 国产精品人妻久久久影院| 国语对白做爰xxxⅹ性视频网站| 亚洲不卡免费看| 亚洲精品美女久久av网站| 国产精品一国产av| 精品少妇久久久久久888优播| 韩国高清视频一区二区三区| 久久精品国产亚洲av涩爱| 中文欧美无线码| 青青草视频在线视频观看| 麻豆乱淫一区二区| 日韩熟女老妇一区二区性免费视频| 精品久久久久久久久亚洲| 校园人妻丝袜中文字幕| 国产在线一区二区三区精| 狠狠婷婷综合久久久久久88av| 99热国产这里只有精品6| 免费观看性生交大片5| 飞空精品影院首页| 国产精品99久久久久久久久| 亚洲欧洲国产日韩| 美女内射精品一级片tv| 97在线人人人人妻| 大香蕉久久网| 一本色道久久久久久精品综合| 一本大道久久a久久精品| 亚洲av成人精品一区久久| 免费观看在线日韩| 国产片内射在线| 少妇 在线观看| 成年女人在线观看亚洲视频| 永久免费av网站大全| 啦啦啦在线观看免费高清www| 久久久久网色| av专区在线播放| 亚洲精品视频女| 制服人妻中文乱码| 久久久久久人妻| 国产成人精品在线电影| 久久久久视频综合| 成年人午夜在线观看视频| 大又大粗又爽又黄少妇毛片口| 母亲3免费完整高清在线观看 | 婷婷色麻豆天堂久久| 亚洲精品国产色婷婷电影| 丰满饥渴人妻一区二区三| 日韩,欧美,国产一区二区三区| 久久精品久久久久久噜噜老黄| 91久久精品国产一区二区成人| 亚洲av在线观看美女高潮| 国产熟女欧美一区二区| 久久久国产一区二区| 在线 av 中文字幕| 久热这里只有精品99| 欧美亚洲 丝袜 人妻 在线| 美女xxoo啪啪120秒动态图| 欧美亚洲日本最大视频资源| 人体艺术视频欧美日本| 91久久精品国产一区二区三区| 少妇人妻精品综合一区二区| www.av在线官网国产| 菩萨蛮人人尽说江南好唐韦庄| 一级毛片我不卡| 国产成人精品一,二区| 国产免费福利视频在线观看| 人人澡人人妻人| 中文字幕免费在线视频6| 久久久久久久久久人人人人人人| 欧美国产精品一级二级三级| 久久久久久人妻| av.在线天堂| 男女边吃奶边做爰视频| 如何舔出高潮| 国产女主播在线喷水免费视频网站| 男女高潮啪啪啪动态图| 九九爱精品视频在线观看| 亚洲av二区三区四区| 女人精品久久久久毛片| 性色av一级| www.av在线官网国产| 美女国产视频在线观看| 97精品久久久久久久久久精品| a级片在线免费高清观看视频| 亚洲天堂av无毛| 一二三四中文在线观看免费高清| 成人综合一区亚洲| 制服诱惑二区| 亚洲av免费高清在线观看| 考比视频在线观看| 又大又黄又爽视频免费| 丝袜脚勾引网站| 在线播放无遮挡| 3wmmmm亚洲av在线观看| 国语对白做爰xxxⅹ性视频网站| 日韩精品免费视频一区二区三区 | 亚洲国产精品专区欧美| 亚洲精品第二区| 亚洲精品美女久久av网站| 一级毛片 在线播放| 下体分泌物呈黄色| 国产精品.久久久| 色吧在线观看| 丰满饥渴人妻一区二区三| 国产白丝娇喘喷水9色精品| 爱豆传媒免费全集在线观看| 久久精品久久久久久久性| videosex国产| 伦理电影免费视频| 乱码一卡2卡4卡精品| 天天操日日干夜夜撸| 18禁在线播放成人免费| tube8黄色片| 国产黄色免费在线视频| 一本久久精品| 亚洲性久久影院| 精品一区二区免费观看| 性色av一级| 久久久久久久久久人人人人人人| 一级二级三级毛片免费看| 内地一区二区视频在线| 日韩一本色道免费dvd| av国产精品久久久久影院| 精品人妻一区二区三区麻豆| 国产精品偷伦视频观看了| 91精品国产国语对白视频| 国产日韩欧美在线精品| 日本午夜av视频| 亚洲欧美清纯卡通| 国产探花极品一区二区| 2018国产大陆天天弄谢| 免费黄网站久久成人精品| 欧美xxⅹ黑人| 久久97久久精品| 国产高清三级在线| 26uuu在线亚洲综合色| av不卡在线播放| 久久亚洲国产成人精品v| 国产国语露脸激情在线看| 亚洲av免费高清在线观看| 尾随美女入室| 大码成人一级视频| 亚洲av不卡在线观看| 大话2 男鬼变身卡| 一二三四中文在线观看免费高清| 亚洲精品一二三| 99热网站在线观看| 妹子高潮喷水视频| 免费观看在线日韩| 91精品伊人久久大香线蕉| 成人综合一区亚洲| 少妇人妻久久综合中文| 三级国产精品片| 欧美最新免费一区二区三区| 成人漫画全彩无遮挡| 亚洲国产精品国产精品| 欧美日韩在线观看h| av电影中文网址| 国产伦理片在线播放av一区| 五月玫瑰六月丁香| 国产亚洲精品久久久com| 九色亚洲精品在线播放| 亚洲在久久综合| 纵有疾风起免费观看全集完整版| 免费不卡的大黄色大毛片视频在线观看| 国产亚洲欧美精品永久| 天美传媒精品一区二区| a 毛片基地| 亚洲av成人精品一区久久| 成人亚洲精品一区在线观看| 18在线观看网站| 国国产精品蜜臀av免费| 午夜激情av网站| 午夜免费观看性视频| 日韩av不卡免费在线播放| 中文精品一卡2卡3卡4更新| 性色avwww在线观看| av在线app专区| 美女脱内裤让男人舔精品视频| 日本vs欧美在线观看视频| 国模一区二区三区四区视频| 久久久久久久久久久丰满| 国产免费现黄频在线看| 草草在线视频免费看| 久久青草综合色| 97在线人人人人妻| 99热国产这里只有精品6| 狂野欧美激情性bbbbbb| 高清午夜精品一区二区三区| 日日爽夜夜爽网站| 欧美变态另类bdsm刘玥| 精品午夜福利在线看| 亚洲精品中文字幕在线视频| 亚洲av中文av极速乱| 午夜老司机福利剧场| 欧美激情 高清一区二区三区| 桃花免费在线播放| 在线亚洲精品国产二区图片欧美 | 一本一本综合久久| 插逼视频在线观看| 久久韩国三级中文字幕| 国产精品99久久久久久久久| 久久ye,这里只有精品| av黄色大香蕉| 少妇丰满av| 麻豆精品久久久久久蜜桃| 一级毛片黄色毛片免费观看视频| 卡戴珊不雅视频在线播放| 另类亚洲欧美激情| 免费不卡的大黄色大毛片视频在线观看| 亚洲国产欧美日韩在线播放| 最新的欧美精品一区二区| 国产精品久久久久成人av| 日本黄色日本黄色录像| 国产精品三级大全| 一级爰片在线观看| 中国三级夫妇交换| 老熟女久久久| 最黄视频免费看| 国产 精品1| 色网站视频免费| 色视频在线一区二区三区| 99久久精品一区二区三区| 日本wwww免费看| 另类精品久久| 亚洲av中文av极速乱| 一级毛片我不卡| 成年女人在线观看亚洲视频| 国产色婷婷99| 国产亚洲精品久久久com| 国内精品宾馆在线| 日本欧美视频一区| 天天躁夜夜躁狠狠久久av| 日韩强制内射视频| 中文字幕亚洲精品专区| 色94色欧美一区二区| 日韩不卡一区二区三区视频在线| 夜夜骑夜夜射夜夜干| 在线观看人妻少妇| 国产永久视频网站| 亚洲欧美日韩卡通动漫| 久久久国产欧美日韩av| 18在线观看网站| 美女福利国产在线| 新久久久久国产一级毛片| 午夜福利网站1000一区二区三区| 日韩大片免费观看网站| 国产精品嫩草影院av在线观看| 在线观看免费视频网站a站| 在线看a的网站| 久久久国产一区二区| 99久国产av精品国产电影| 成人综合一区亚洲| 美女大奶头黄色视频| 天美传媒精品一区二区| 成年av动漫网址| 美女脱内裤让男人舔精品视频| 精品久久久久久电影网| 久久综合国产亚洲精品| 欧美丝袜亚洲另类| 一个人看视频在线观看www免费| 又大又黄又爽视频免费| 2021少妇久久久久久久久久久| 久久久精品免费免费高清| 大香蕉久久成人网| 亚洲婷婷狠狠爱综合网| 精品一区二区三卡| 亚洲精品自拍成人| 久久青草综合色| 国产亚洲av片在线观看秒播厂| 国产精品 国内视频| 国产成人av激情在线播放 | 成年女人在线观看亚洲视频| 欧美精品一区二区免费开放| 国产成人a∨麻豆精品| 国产成人精品一,二区| 成人手机av| 亚洲精品国产av成人精品| 久久亚洲国产成人精品v| 国产一区亚洲一区在线观看| 97精品久久久久久久久久精品| 国产免费现黄频在线看| 九色成人免费人妻av| 国产精品嫩草影院av在线观看| 久久精品熟女亚洲av麻豆精品| 男的添女的下面高潮视频| 国产精品女同一区二区软件| 免费黄网站久久成人精品| 欧美日韩一区二区视频在线观看视频在线| 97超视频在线观看视频| 午夜福利影视在线免费观看| h视频一区二区三区| av网站免费在线观看视频| 视频中文字幕在线观看| 老女人水多毛片| 久久久午夜欧美精品| av天堂久久9| 国产乱人偷精品视频| 久久热精品热| 精品卡一卡二卡四卡免费| 久久久久国产网址| 亚洲精品av麻豆狂野| 国产一区二区三区综合在线观看 | 啦啦啦啦在线视频资源| 欧美人与性动交α欧美精品济南到 | 久久女婷五月综合色啪小说| 免费观看无遮挡的男女| 波野结衣二区三区在线| 日韩免费高清中文字幕av| 国产精品99久久久久久久久| 亚洲国产精品一区三区| 日韩成人伦理影院| 亚洲情色 制服丝袜| 欧美老熟妇乱子伦牲交| 天堂8中文在线网| 欧美老熟妇乱子伦牲交| 日韩成人av中文字幕在线观看| 欧美+日韩+精品| 国产精品一区二区在线不卡| 国产亚洲精品第一综合不卡 | 国产精品偷伦视频观看了| 热re99久久国产66热| 欧美精品高潮呻吟av久久| 亚洲av中文av极速乱| 大片免费播放器 马上看| 久久精品熟女亚洲av麻豆精品| 欧美日韩精品成人综合77777| 午夜福利影视在线免费观看| 蜜桃在线观看..| 中国三级夫妇交换| 亚洲精品中文字幕在线视频| 色网站视频免费| 成人18禁高潮啪啪吃奶动态图 | a级片在线免费高清观看视频| 免费播放大片免费观看视频在线观看| 亚洲精品日本国产第一区| 九草在线视频观看| 七月丁香在线播放| 免费看av在线观看网站| 亚洲一区二区三区欧美精品| 青春草国产在线视频| av播播在线观看一区| 婷婷色av中文字幕| 国产一区二区在线观看日韩| 亚洲美女搞黄在线观看| 国产成人91sexporn| 国产亚洲av片在线观看秒播厂| 成人国语在线视频| 一级毛片aaaaaa免费看小| 一级毛片电影观看| 人人妻人人爽人人添夜夜欢视频| 亚洲精品久久午夜乱码| 大香蕉久久网| 最近手机中文字幕大全| 亚洲av欧美aⅴ国产| 午夜激情av网站| 婷婷色麻豆天堂久久| 久久久精品区二区三区| 26uuu在线亚洲综合色| 亚洲图色成人| 大码成人一级视频| 精品人妻在线不人妻| av天堂久久9| 久久久久久久久久成人| av卡一久久| 大陆偷拍与自拍| 亚洲欧美清纯卡通| 天天躁夜夜躁狠狠久久av| 一级黄片播放器| 制服诱惑二区| 亚洲国产精品国产精品| 亚洲av福利一区| 国产淫语在线视频| 18禁在线播放成人免费| 女性生殖器流出的白浆| 精品久久久精品久久久| 中文天堂在线官网| 妹子高潮喷水视频| 日韩制服骚丝袜av| 黄色欧美视频在线观看| 在线观看国产h片| 日本黄色日本黄色录像| 国模一区二区三区四区视频| 亚洲精品aⅴ在线观看| 九色亚洲精品在线播放| 欧美日韩视频精品一区| 97在线人人人人妻| 国产精品免费大片| 欧美人与善性xxx| 午夜福利视频精品| 寂寞人妻少妇视频99o| 精品人妻熟女av久视频| 亚洲色图 男人天堂 中文字幕 | 一级毛片aaaaaa免费看小| 91久久精品国产一区二区三区| 久久久久久久久久成人| 99九九在线精品视频| 久久人妻熟女aⅴ| 人妻少妇偷人精品九色| 欧美日韩视频精品一区| 精品亚洲成a人片在线观看| 亚洲精品色激情综合| 一级毛片aaaaaa免费看小| 成人漫画全彩无遮挡| 国产成人精品福利久久| 亚洲精品国产av成人精品| 久久久久久久久久久久大奶| 欧美亚洲日本最大视频资源| 精品久久久久久电影网| 精品一品国产午夜福利视频| 美女xxoo啪啪120秒动态图| 女性生殖器流出的白浆| 五月天丁香电影| 欧美日韩成人在线一区二区| 亚洲精品乱码久久久v下载方式| 久久久久视频综合| 日韩成人伦理影院| 精品卡一卡二卡四卡免费| 精品少妇久久久久久888优播| 伦理电影大哥的女人| 久久国内精品自在自线图片| 久久久午夜欧美精品| 中文字幕亚洲精品专区| 亚洲一级一片aⅴ在线观看| 久久久精品免费免费高清| 久久久久久久精品精品| 国产av码专区亚洲av| 日韩 亚洲 欧美在线| 一区在线观看完整版| 91成人精品电影| 久久久精品94久久精品| 女性生殖器流出的白浆| 男人爽女人下面视频在线观看| 国产精品偷伦视频观看了| 丝瓜视频免费看黄片| 亚洲av福利一区| av视频免费观看在线观看| 中文精品一卡2卡3卡4更新| 国产免费又黄又爽又色| 美女主播在线视频| 中文字幕最新亚洲高清| 欧美 亚洲 国产 日韩一| 纵有疾风起免费观看全集完整版| 女性被躁到高潮视频| 精品熟女少妇av免费看| 日韩亚洲欧美综合| 美女脱内裤让男人舔精品视频| 久久国内精品自在自线图片| 国产免费一级a男人的天堂| 欧美精品亚洲一区二区| 亚洲精品日韩在线中文字幕| 麻豆精品久久久久久蜜桃| 欧美亚洲 丝袜 人妻 在线| 51国产日韩欧美| 国产永久视频网站| 国产老妇伦熟女老妇高清| 国产免费福利视频在线观看| 热99国产精品久久久久久7| 欧美三级亚洲精品| 一区二区三区乱码不卡18| .国产精品久久| videossex国产| 色94色欧美一区二区| 成人18禁高潮啪啪吃奶动态图 | 岛国毛片在线播放| 人人妻人人添人人爽欧美一区卜| 97超碰精品成人国产| 丰满乱子伦码专区| 久久久久国产网址| 曰老女人黄片| 国产一区二区在线观看日韩| 国产淫语在线视频| 少妇猛男粗大的猛烈进出视频| 两个人的视频大全免费| 亚洲国产成人一精品久久久| 日韩av免费高清视频| 黄色配什么色好看| 在线 av 中文字幕| 欧美bdsm另类| 在线观看免费视频网站a站| 中文字幕av电影在线播放| 永久网站在线| 香蕉精品网在线| 又粗又硬又长又爽又黄的视频| 99re6热这里在线精品视频| 国产在视频线精品| 日本av免费视频播放| 在线观看一区二区三区激情| 免费观看性生交大片5| 成人毛片60女人毛片免费| 成人二区视频| 啦啦啦在线观看免费高清www| 国产一区二区三区综合在线观看 | 一区二区日韩欧美中文字幕 | 久久亚洲国产成人精品v| 午夜老司机福利剧场| 只有这里有精品99| 久久久a久久爽久久v久久| 欧美日韩视频高清一区二区三区二| 久久人妻熟女aⅴ| 婷婷色av中文字幕| 啦啦啦在线观看免费高清www| 人妻制服诱惑在线中文字幕| 午夜激情av网站| 亚洲怡红院男人天堂| 久久狼人影院| 啦啦啦视频在线资源免费观看| 母亲3免费完整高清在线观看 | av电影中文网址| 国产午夜精品久久久久久一区二区三区| 国产高清三级在线| 久久精品久久精品一区二区三区| 男的添女的下面高潮视频| 亚洲精品乱码久久久v下载方式| 一级片'在线观看视频| 九草在线视频观看| 最近的中文字幕免费完整| 亚洲综合色惰| 欧美人与善性xxx| 精品一区二区三区视频在线| 亚洲欧美成人综合另类久久久| 日韩欧美一区视频在线观看| 精品卡一卡二卡四卡免费| 只有这里有精品99| 国产免费又黄又爽又色| 国产免费视频播放在线视频| 久久久午夜欧美精品| 又黄又爽又刺激的免费视频.| 日本av手机在线免费观看| 在线看a的网站| 亚洲美女搞黄在线观看| 中国国产av一级| 亚洲欧美一区二区三区黑人 | 欧美丝袜亚洲另类| 啦啦啦在线观看免费高清www| 成人免费观看视频高清| 青春草国产在线视频| 国产男人的电影天堂91| 成人亚洲精品一区在线观看| 国产精品一二三区在线看| 啦啦啦在线观看免费高清www| 日韩熟女老妇一区二区性免费视频| 久久人人爽人人片av| 国产精品国产三级国产av玫瑰| 不卡视频在线观看欧美| 最近2019中文字幕mv第一页| 国产成人91sexporn| 精品人妻熟女毛片av久久网站| 欧美精品亚洲一区二区| 99久久中文字幕三级久久日本| 久久精品国产亚洲av涩爱| 国产成人aa在线观看| 日韩大片免费观看网站| 九草在线视频观看| 久久精品国产亚洲网站| 国产成人精品婷婷| 亚洲人成77777在线视频| 超碰97精品在线观看| √禁漫天堂资源中文www| 国产av精品麻豆| 久久久久久久大尺度免费视频| 日韩强制内射视频| 寂寞人妻少妇视频99o| 99久久人妻综合| 黄片播放在线免费| 天天躁夜夜躁狠狠久久av| 又黄又爽又刺激的免费视频.| 久久国产亚洲av麻豆专区|