劉啟康
(1、北京建筑大學(xué) 測繪與城市空間信息學(xué)院,北京100044 2、北京北建大科技有限公司,北京100044)
近年來,目標(biāo)識別在圖像檢索、機器人抓取操作、位置識別與定位等眾多領(lǐng)域得到廣泛的研究。目前,基于圖像的物體識別已經(jīng)有很多研究成果。但是二維圖像的成像過程是從三維空間映射到二維空間,這個過程會丟失大量的信息。而目前三維點云數(shù)據(jù)的獲取已經(jīng)非??旖荩瑫r三維點云數(shù)據(jù)的采集不受光照影響,規(guī)避了二維圖像遇到的光照、姿態(tài)等問題,因此基于點云數(shù)據(jù)的三維物體識別也引起了人們的重視。
目前,已有的三維點云物體識別方法多是通過分析/提取物體的特征點幾何屬性、形狀屬性、結(jié)構(gòu)屬性或者多種屬性的組合等特征進(jìn)行比對、學(xué)習(xí),從而完成物體的識別與分類。根據(jù)所利用特征的不同,已有的三維點云物體識別方法可以分為4類:基于局部特征的物體識別方法、基于全局特征的物體識別方法、基于圖匹配的物體識別方法以及基于機器學(xué)習(xí)的物體識別方法。
基于局部特征的物體識別方法主要是通過局部來識別整體[1]。該方法無需對處理數(shù)據(jù)進(jìn)行分割,往往通過提取物體的關(guān)鍵點,利用局部特征描述符描述特征點,再通過特征匹配、假設(shè)驗證等過程完成物體的識別。
圖1
其中,特征提取與描述是物體識別中非常關(guān)鍵的一步,它將直接影響到物體識別系統(tǒng)的性能。
本文提出的基于局部特征的點云目標(biāo)識別方法是通過對模板點云與場景點云建立全局約束條件,通過八叉樹體素化方式選取特征點,使用SHOT[2]描述符描述模板特征點與場景點的幾何屬性,通過特征匹配代價和幾何匹配代價建立局部匹配模型,求解局部匹配模型最優(yōu)解,最后結(jié)合全局約束條件和局部匹配模型建立最終匹配模型,以獲得最佳匹配目標(biāo)。
本方法提出的識別流程如圖1 所示。
將預(yù)處理后的點云數(shù)據(jù)均勻采樣,以提高PCA 計算軸向和特征值的精確度,減少數(shù)據(jù)采集帶來的誤差影響,然后利用PCA 計算OBB 包圍盒的體積V 以及數(shù)據(jù)特征值Xi,并建立全局約束方程:
基于八叉樹結(jié)構(gòu),將模板點與場景點劃分為體素大小為V的體素結(jié)構(gòu),對于每個體素,統(tǒng)計點云數(shù)量,計算體素的質(zhì)心,以最接近質(zhì)心的點作為特征點。采用此方式選取特征點可以大規(guī)模減少點云數(shù)量,大大提高點云特征描述時的效率,加快目標(biāo)物體的匹配速度。
本文選用局部特征描述子描述點云特征,SHOT 描述子是通過計算特征點局部鄰域曲面的協(xié)方差矩陣特征向量來定義LRF,通過對LRF 重要性的分析,利用加權(quán)線性組合算法來計算協(xié)方差矩陣與符號異性消除。該方法具有旋轉(zhuǎn)和平移不變性,對噪聲和目標(biāo)相互遮擋等問題具有魯棒性,描述性較強。
局部目標(biāo)函數(shù)由兩部分組成:特征匹配代價,幾何匹配代價。
特征匹配代價目前常用的幾種方法有:求交運算、歐氏距離平方、卡方距離、KL 散度。本文將使用卡方距離來表示特征匹配代價,由于其對噪聲的敏感度較低,在高噪聲的情況下性能較高的特性,用以描述直方圖差異,公式如下:
由于局部目標(biāo)函數(shù)是具有整數(shù)約束的非線性目標(biāo)函數(shù),這樣的公式會導(dǎo)致NP 問題,在多項式時間內(nèi)不能有效解決,因此我們采用取絕對值法將引入輔助變量,將非線性目標(biāo)函數(shù)轉(zhuǎn)化成線性目標(biāo)函數(shù):
最終由全局約束方程和局部約束函數(shù)建立總體目標(biāo)函數(shù):
目前比較流行的求解線性規(guī)劃問題的方案有很多,例如:CPLEX、GLPK、Gurobi、LPSolve。本文采用開源、高效、便捷的LPSolve 解決線性規(guī)劃問題,以獲得最終匹配結(jié)果。
本文以地下工井為實驗對象,以A7651、A1790、A7622、A7662、A7752 五組不同類型的工井作為模板數(shù)據(jù),以十組包含模板類型并存在旋轉(zhuǎn)變換的工井作為場景數(shù)據(jù),匹配結(jié)果如圖2 所示。
圖2
實驗數(shù)據(jù)表明:對于A7651、A1790、A7622、A7662、A7752五種不同類型的工井,可以準(zhǔn)確的得到與之相匹配的工井類型A7652、A1791、A7623、A7663、A7753。