劉信偉 ，靖常峰，羅德利，杜明義，蔡國印

(1.北京建筑工程學(xué)院測(cè)繪學(xué)院，北京 100044;2.北京市市政專業(yè)設(shè)計(jì)院股份公司，北京 100037)

1 引言

特征點(diǎn)是幾何形狀的特征基元，它不因坐標(biāo)系的改變而變化。傳統(tǒng)特征點(diǎn)提取方法，如目視判讀或相似度匹配等，由于人為因素和相似度函數(shù)的誤差，所提取特征點(diǎn)精度受限。Woo［1］認(rèn)為測(cè)量點(diǎn)的法矢或曲率的突變是區(qū)域的邊界，提出將法矢或曲率的突變點(diǎn)作為特征點(diǎn)。馬驪溟［2］采用高斯曲率的方法，在散亂點(diǎn)云數(shù)據(jù)中提取特征點(diǎn)。Huang［3］在完成數(shù)據(jù)點(diǎn)三角網(wǎng)格化的基礎(chǔ)上，估算各測(cè)點(diǎn)的法矢和曲率，把曲率極值點(diǎn)作為邊界特征點(diǎn)。這些方法均是直接提取方法，直接采用掃描測(cè)量點(diǎn)作為特征點(diǎn)，精度受限于測(cè)量誤差，因此提取的特征點(diǎn)未必是曲面真正的特征點(diǎn)。本文從這個(gè)問題出發(fā)，首先對(duì)點(diǎn)云數(shù)據(jù)進(jìn)行預(yù)處理，去除噪聲數(shù)據(jù)和非感興趣數(shù)據(jù);然后確定適當(dāng)大小的鄰域，通過計(jì)算兩個(gè)相鄰點(diǎn)之間的高斯曲率K和平均曲率H，判斷出工作區(qū)內(nèi)特征點(diǎn)的大致位置，最后擬合目標(biāo)物體的局部曲面，解算出該曲面的極值特征點(diǎn)。

2 點(diǎn)云的預(yù)處理

三維激光掃描儀獲取的點(diǎn)云數(shù)據(jù)，通常包含測(cè)量過程中產(chǎn)生的噪聲數(shù)據(jù)，影響表面重建［4～5］。此外，在三維點(diǎn)云數(shù)據(jù)中，還存在掃描物體之外其他物體的點(diǎn)云數(shù)據(jù)，這些數(shù)據(jù)雖然不屬于噪聲數(shù)據(jù)，但我們不感興趣，需要將其去除［6～8］。因此，我們?cè)诰植奎c(diǎn)云數(shù)據(jù)的擬合之前，首先進(jìn)行點(diǎn)云數(shù)據(jù)的預(yù)處理［9～12］。點(diǎn)云的預(yù)處理的步驟如圖1所示。

圖1 點(diǎn)云預(yù)處理步驟

點(diǎn)云預(yù)處理工作是特征點(diǎn)提取的先決條件。預(yù)處理效果好，則特征點(diǎn)提取和后處理精度高、效率高，擬合曲面更接近被掃描物體的真實(shí)情況，所提取的曲面極值特征點(diǎn)更接近于物體的真實(shí)特征點(diǎn)，否則反之。

3 曲面極值特征點(diǎn)的獲取

選用三角格網(wǎng)模型模擬表達(dá)被掃描物體表面是一種常用方法，模型簡潔并且可以很好的表達(dá)高度不規(guī)則物體表面的拓?fù)潢P(guān)系［13］。三角格網(wǎng)模型如圖2所示。

圖2 河南洛陽一佛像的貼圖后三角格網(wǎng)模型

每個(gè)頂點(diǎn)Di，在其周圍有N個(gè)網(wǎng)格頂點(diǎn)，基于它們的局部關(guān)系確定頂點(diǎn)Di(xi，yi，zi)。用函數(shù)式z=f(x，y)表示Di的局部曲面。二階多項(xiàng)式:

采用最小二乘法，求解多項(xiàng)式函數(shù)式(1)的系數(shù)。多項(xiàng)式(1)中的系數(shù)確定，準(zhǔn)確的頂點(diǎn)D'i(Di在擬合曲面中的位置)的坐標(biāo)位置也將確定?？梢杂萌缦露囗?xiàng)式表示。

頂點(diǎn)Di附近的N個(gè)鄰近點(diǎn)的精度和合適的擬合多項(xiàng)式函數(shù)是擬合最佳曲面的關(guān)鍵。根據(jù)需要，選擇合適的控制因素以獲得期望的結(jié)果。多次平滑也可以應(yīng)用到三角格網(wǎng)建立的過程中，以獲得更加光滑的表面模型。平滑區(qū)域的大小通常是直徑2個(gè)～3個(gè)三角格網(wǎng)的圓曲面。這樣的曲面可以大大降低局部噪聲對(duì)明顯局部特征的影響。

三角格網(wǎng)曲面經(jīng)過平滑后，削弱了測(cè)量噪聲，消除了小的幾何特征，較好地保存了主要的表面特征數(shù)據(jù)。這樣較明顯地減小特征匹配搜索過程的復(fù)雜性，獲得可靠的初始剛體變換值。通過平滑確實(shí)改變了測(cè)量數(shù)據(jù)點(diǎn)的位置，但是它的改變是十分微小的(我們可以通過擬合曲面函數(shù)式解求這個(gè)微小變化)，因此，平滑而引起的局部幾何變形不會(huì)對(duì)最終的配準(zhǔn)產(chǎn)生影響。

得到較好的擬合曲面后，根據(jù)高等數(shù)學(xué)中解求極值點(diǎn)的方法，進(jìn)行計(jì)算。根據(jù)式(2)，可得:

由式(3)和式(4)得到極值點(diǎn)的x值和y值，把x值和y值代入式(2)解算出z值，這樣就得到了曲面的極值特征點(diǎn)坐標(biāo)(x，y，z)。

4 擬合曲面的參數(shù)求解

對(duì)取出的三維點(diǎn)云數(shù)據(jù)進(jìn)行分析處理，然后根據(jù)最小二乘原理，解求函數(shù)式式(2)中的系數(shù)，確定曲面的顯示表達(dá)式。

式中:n為點(diǎn)數(shù);N為系數(shù)個(gè)數(shù):n－N為多余觀測(cè)。

設(shè)定一個(gè)限差ε作為評(píng)定精度的標(biāo)準(zhǔn)。本文在做實(shí)驗(yàn)時(shí)，限差ε的取值是點(diǎn)云掃描精度的1/2。若δ＞ε，則說明存在粗差，精度不可取，應(yīng)對(duì)每個(gè)測(cè)量點(diǎn)的平差殘余誤差vz進(jìn)行比較檢查，最大者為粗差，將其剔除或重新選點(diǎn)后再進(jìn)行平差，直至滿足δ＜ε為止。

這樣求解出的參數(shù)，可以與式(2)聯(lián)合解算，求解擬合空間曲面的極值特征點(diǎn)。

5 實(shí)驗(yàn)結(jié)果及精度分析

5.1 實(shí)驗(yàn)

利用Geomagic Studio 9.0從佛像的三維點(diǎn)云數(shù)據(jù)中取出部分的點(diǎn)云數(shù)據(jù)。保存為.obj格式的文件，然后再另存為.txt格式的文本文件。

通過上面介紹的參數(shù)求解的方法，求解擬合曲面的函數(shù)表達(dá)式，進(jìn)而求得擬合曲面的極值點(diǎn)坐標(biāo)。同時(shí)分析擬合曲面及其極值特征點(diǎn)的精度。

編程求解出的擬合曲面的函數(shù)式為:

5.2 精度分析

將取出的點(diǎn)云數(shù)據(jù)中未參與擬合計(jì)算的部分點(diǎn)作為檢核點(diǎn)，進(jìn)行檢核，求解其擬合精度。表1為部分點(diǎn)的檢核情況。

擬合曲面的精度分析 表1

根據(jù)式(8)，求解擬合曲面的擬合精度為:0.0009374。擬合曲面極值點(diǎn)坐標(biāo)為:(－0.080337，－2.220626，2.749954)。

6 討論及結(jié)論

預(yù)處理的必要性。一般來說，三角形網(wǎng)格是通過含有噪聲數(shù)據(jù)的點(diǎn)云數(shù)據(jù)構(gòu)建的。由于頂點(diǎn)位置處可能存有噪聲，導(dǎo)致產(chǎn)生大量的小特征碎片，每個(gè)特征碎片含幾個(gè)三角形切面，具有相同的局部表面類型。這些特征碎片不是局部特征的真實(shí)描述，并且在其他網(wǎng)格上也沒有與之對(duì)應(yīng)的特征，增加了特征匹配的虛假率。為了提高特征匹配的效率和可靠性，需要采取一些算法去除這些小的特征碎片。

曲面擬合前初始判斷的必要性。要想進(jìn)一步提高結(jié)果的可靠性，還應(yīng)對(duì)曲面的頂點(diǎn)或是谷點(diǎn)進(jìn)行初始判斷。求取每兩個(gè)相鄰點(diǎn)之間的高斯曲率K和平均曲率H，高斯曲率K和平均曲率H的不同值的組合代表8種不同的曲面類型。局部曲面類型可以分為8類，可以提供一個(gè)離散的搜索區(qū)間，判斷曲面類型。然后，在判斷后的頂點(diǎn)附近取10個(gè)左右的點(diǎn)云數(shù)據(jù)，用這些數(shù)據(jù)進(jìn)行曲面擬合，求得曲面極值特征點(diǎn)，這樣所求結(jié)果可靠性會(huì)更高。

在做實(shí)驗(yàn)時(shí)，點(diǎn)云的預(yù)處理，坐標(biāo)系的轉(zhuǎn)換，曲面特征的初始判斷以及曲面的擬合等，每一步都十分關(guān)鍵。每一步處理的好壞，都直接影響提取特征點(diǎn)的精度。

［1］Woo H，Kang E，Wang Sem － yung，et al.A New Segmentation Method for Point Cloud Data［J］.International Journal of Machine Tools andManufacture(S0890 －6955)，2002，42(2):167～178.

［2］馬驪溟，徐毅，李澤湘.基于高斯曲率極值點(diǎn)的散亂點(diǎn)云數(shù)據(jù)特征點(diǎn)提取.系統(tǒng)仿真學(xué)報(bào)，2008，20(9):2341～2344.

［3］Huang J，Menq C H.Automatic Data Segmentation for Geometric Feature Extraction from Unorganized 3－D Coordinate Points［J］.IEEE Transactions on Robotics and Automation(S1042 －296X)，2001，17(3):268 ～279.

［4］張毅，劉旭敏，隋穎等.基于K－近鄰點(diǎn)云去噪算法的研究與改進(jìn)［J］.計(jì)算機(jī)應(yīng)用，2009，29(4):1011～1012.

［5］T.Tasdizen，R.Whitaker，P.Burchard，S.Osher.Geometric Surface Smoothing via Anisotropic Diffusion of Normals.Proceedings of IEEE Conference on Visualization，2002，125 ～132.

［6］劉含波.基于散亂點(diǎn)云數(shù)據(jù)的隱式曲面重建研究［D］.哈爾濱:哈爾濱工業(yè)大學(xué)，2009.

［7］Hanbo Liu，Xin Wang and Wenyi Qiang.Implicit Surface Reconstruction from 3D Scattered Points Based on Variational Level Set Method.ISSCAA，2008，641 ～644.

［8］Wang Xin，Yang Jian，Liu Hanbo and Ma Yan.A Fast Stereo Matching Algorithm for Real－ time Robot Application.2007 IEEE InternationalConference on Robotics and Biomimetics，2007，908 ～913.

［9］張鴻飛，程效軍，賈東峰.多視點(diǎn)散亂點(diǎn)云配準(zhǔn)及壓縮改進(jìn)算法研究［J］.測(cè)繪通報(bào)，2012(2):43～47.

［10］吳世雄，王文，陳子辰等.大規(guī)模掃描測(cè)點(diǎn)的自適應(yīng)數(shù)據(jù)壓縮［J］.浙江大學(xué)學(xué)報(bào)·工學(xué)版，2004，38(9):1200～1204.

［11］邵正偉，席平.基于八叉樹編碼的點(diǎn)云數(shù)據(jù)精簡方法［J］.工程圖學(xué)學(xué)報(bào)，2010(4):73～76.

［12］師振中，王秀英，劉錫國.逆向工程中點(diǎn)云數(shù)據(jù)壓縮算法的研究與改進(jìn)［J］.江蘇大學(xué)學(xué)報(bào)·自然科學(xué)版，2006，27(B09):35 ～39.

［13］N.Li，P.cheng，M.A.Sutton，S.R.McNeill.Three － dimensional Point Cloud Registration by Matching Surface Features with Relaxation Labling Method.Society for Experimental Mechanics，2005，2:71 ～82.