陳 濤 劉桂昌 張 倩 魏 朗
(長安大學(xué)汽車學(xué)院交通部汽車運(yùn)輸安全保障技術(shù)重點(diǎn)實(shí)驗(yàn)室,陜西 西安 710064)
目前,便攜式三維掃描儀的應(yīng)用非常廣泛,主要用于逆向工程、實(shí)物拍攝、創(chuàng)建產(chǎn)品的三維模型和有限元分析的數(shù)據(jù)捕捉等。這些領(lǐng)域?qū)θS掃描儀的數(shù)據(jù)精度和分辨率精度都有很高的要求,因此,在選用掃描儀時必須確認(rèn)其精度和分辨率,才能獲得滿意的三維數(shù)據(jù)模型[1]。
隨著掃描儀技術(shù)的日趨成熟,尚缺少與便攜式激光掃描儀分辨率指標(biāo)相關(guān)的國際或國家標(biāo)準(zhǔn)出臺。國內(nèi)外對分辨率的研究主要傾向于對分辨率影響因素的研究。如通過對比不同種類掃描儀測試結(jié)果來分辨掃描儀分辨率優(yōu)劣[2],分析掃描儀掃描距離和角度對掃描儀分辨率的影響[3-4],以及利用計(jì)算公式來計(jì)算一定面積上點(diǎn)云的密度與掃描儀分辨率間的關(guān)系[5]。但通過研究掃描儀分辨率的影響因素,仍不能給出一個簡單、直觀、可操作的儀器使用判別方法。本文在進(jìn)行大量試驗(yàn)的基礎(chǔ)上,借助科學(xué)的統(tǒng)計(jì)方法,給出了一種可行的方案。
便攜式的激光掃描儀掃描系統(tǒng)主要由掃描儀、控制器和電源供應(yīng)系統(tǒng)三部分組成。激光掃描儀本身主要包括激光測距系統(tǒng)和激光掃描系統(tǒng),同時,也集成了CCD以及儀器內(nèi)部控制和校正等系統(tǒng)[6]。
三維激光掃描儀所采用的最基本的方法有飛行時間法(time of flying,ToF)和三角測量法。飛行時間法利用激光發(fā)射到接收之間的時間延遲來計(jì)算距離。但由于光波的飛行速度達(dá)3×108m/s,為達(dá)到較高的距離測量精度,對于定時系統(tǒng)的時間分辨率有特別高的要求,較難在技術(shù)上得到實(shí)現(xiàn)。三角測量法則以傳統(tǒng)的三角測量為基礎(chǔ)。由于三維面型對結(jié)構(gòu)照明光束產(chǎn)生的空間調(diào)制改變了成像光束的角度,即改變了成像光點(diǎn)在接收裝置陣列上的位置,因此通過對成像光點(diǎn)位置以及系統(tǒng)光路的幾何參數(shù)的確定,可最終計(jì)算出掃描對象與激光發(fā)射器的距離[7]。因此,近年來激光三維技術(shù)的發(fā)展主要在于三角測量法的更新。
單獨(dú)的三角測量方式要求投影光軸和成像光軸之間保持恒定的夾角,且必須附加一維或二維的相對掃描來獲取目標(biāo)點(diǎn)的三維數(shù)據(jù)。Handyscan EXAscan等掃描儀大多采用雙三角測量方式來解決單三角測量系統(tǒng)帶來的弊端。
三維掃描儀的原理圖如圖1所示[8-10]。利用三角形邊與角的特定關(guān)系,可以計(jì)算出目標(biāo)點(diǎn)的三維坐標(biāo)。
圖1 三維掃描儀的原理圖Fig.1 Principle of the 3D scanner
分辨率是一個籠統(tǒng)的概念,泛指測量或顯示系統(tǒng)對目標(biāo)物細(xì)節(jié)的分辨能力,較高的分辨率可以更好地體現(xiàn)目標(biāo)物的細(xì)節(jié)。
在不同的領(lǐng)域,分辨率的標(biāo)準(zhǔn)也不同。通常情況下,分辨率指每英寸(1英寸=25.4mm)上可產(chǎn)生的點(diǎn)數(shù)或像素點(diǎn)數(shù),即每英寸點(diǎn)(dots per inch,DPI)。對于光學(xué)設(shè)備,特別是三維掃描儀設(shè)備,分辨率就是對細(xì)節(jié)分辨能力的一種度量,即掃描儀能夠區(qū)分的最小特征參數(shù),也是描述點(diǎn)云中目標(biāo)細(xì)微程度的指標(biāo)[5]。
與平面輸入輸出設(shè)備的數(shù)據(jù)不同,三維掃描的數(shù)據(jù)屬于矢量化的數(shù)據(jù),盡管文獻(xiàn)[5]中提出使用點(diǎn)云密度的方法來衡量分辨率,但實(shí)際操作起來還是具有一定的難度。在此,取相鄰點(diǎn)與點(diǎn)之間的空間距離作為掃描儀的分辨率的指標(biāo)。
鑒定掃描儀分辨率的最重要的工作就是數(shù)據(jù)的處理工作。當(dāng)三維掃描儀工作時,CCD攝像頭拾取掃描對象的過程即拍照的過程。依據(jù)其工作原理,在選取樣本時,通常要考慮選取點(diǎn)云分布近似經(jīng)緯方向的最近的兩點(diǎn)之間的空間距離,而對于對角方向的距離則不予考慮。由于點(diǎn)云的數(shù)據(jù)量龐大,無法逐一檢查,因此通常采用隨機(jī)抽樣的方式。
手持式三維激光掃描儀的分辨率的鑒定流程如下。
①連接掃描儀至計(jì)算機(jī),做好掃描的前期工作。
②獲取測試數(shù)據(jù)。分辨率等級受制于容積框的大小,一定的容積框大小又有三個分辨率級別可選,由低到高依次是L、M、H三個等級。在此基礎(chǔ)上,還有高倍鏡模式(如L+、M+、H+),就是在原有的掃描模式下將分辨率水平提高4倍。將容積框的大小設(shè)定為最小(最小只能到100 mm),在分辨率級別最高的情況下,名義最高分辨率可達(dá)0.20 mm,開高倍鏡則最高分辨率可達(dá)0.05 mm。
③ 導(dǎo)入掃描數(shù)據(jù)至Geomagic Qualify,獲取樣本。對于大數(shù)量的母體,要估計(jì)母體置信區(qū)間,通常選取大子樣作為樣本。根據(jù)實(shí)際經(jīng)驗(yàn),一般認(rèn)為n≥45的子樣是大子樣。
④建立分析模型,處理樣本數(shù)據(jù)。
⑤根據(jù)分析結(jié)果生成分辨率數(shù)據(jù)模型,判定分辨率水平。
加拿大Creaform公司生產(chǎn)的Handyscan EXAscan高精度系列掃描儀是一種高精度的手持式自定位三維激光掃描儀。它在兩個普通的CCD鏡頭的基礎(chǔ)上又增加了一個CCD鏡頭,使掃描儀的最大分辨率在原來的基礎(chǔ)上增加了4倍,這極大地豐富了被掃描對象的細(xì)節(jié)。
Handyscan 3D掃描儀采用的是自定位模式。在掃描過程中,掃描儀會實(shí)時捕捉目標(biāo)點(diǎn),進(jìn)而計(jì)算和記錄其各自相對于掃描儀的位置。這些目標(biāo)點(diǎn)會對掃描物體進(jìn)行定位,定位后即可通過Handyscan 3D掃描儀采集物體表面的三維尺寸數(shù)據(jù)。
Handyscan 3D掃描儀的定位原理是利用不在同一直線上的四個點(diǎn)來確定一個三維坐標(biāo)系,這樣在掃描前或掃描過程中,掃描儀只要能同時捕捉到合適的四個點(diǎn),就能確定這四個點(diǎn)所組成的坐標(biāo)系區(qū)域采樣點(diǎn)的坐標(biāo)值。而所有的定位點(diǎn)又互有聯(lián)系,故所有的小坐標(biāo)系就可以統(tǒng)一成一個擁有共同坐標(biāo)系的空間,采樣點(diǎn)的坐標(biāo)也會隨之轉(zhuǎn)化為公共坐標(biāo)的坐標(biāo)值。因此,采集信息之前需要對目標(biāo)進(jìn)行貼點(diǎn)標(biāo)記。這種掃描特性的優(yōu)點(diǎn)在于被掃描物體可以是運(yùn)動的、空間位置是自由的。
測試Handyscan掃描儀的分辨率等級,必要時需要開啟第三個CCD鏡頭,這對于計(jì)算機(jī)硬件配置有很高的要求。本次試驗(yàn)使用的是HP Z800圖形工作站。工作站的主要配置為:2個英特爾至強(qiáng)四核處理器X5570、12 GB內(nèi)存(DDR3-1333 MHE ECC)和3D NVIDIA Quadro FX5800專業(yè)顯卡。
掃描儀的掃描對象是標(biāo)準(zhǔn)塊。打開VXSCAN掃描軟件(儀器自帶驅(qū)動程序)設(shè)定面細(xì)節(jié)。設(shè)定內(nèi)容包括設(shè)定容積框的大小和分辨率的級別。掃描試驗(yàn)獲取的樣本數(shù)據(jù)涵蓋了容積框邊長分別為100 mm、300 mm、500 mm、700 mm、900 mm、1 000 mm 在 L、M、H和H+四種等級模式下的掃描數(shù)據(jù)。
Geomagic Qualify通過德國標(biāo)準(zhǔn)計(jì)量機(jī)構(gòu)PTB(Physikalisch-Technische Bundesanstalt)認(rèn)證,可以迅速檢測產(chǎn)品的計(jì)算機(jī)輔助設(shè)計(jì)(CAD)模型和制造的產(chǎn)品之間的差異。憑借Geomagic Qualify提供的高級分析功能,可快速判斷出掃描儀的分辨率狀況,為處理掃描數(shù)據(jù)提供理想的工作平臺。
影響掃描精度的因素有很多,影響分辨率的因素也同樣很多,包括被掃描物的反光、掃描距離、溫度及對象顏色等[2,11-12]。
為減少外界干擾因素對測試結(jié)果的影響,應(yīng)盡量在點(diǎn)云分布較均勻的部位選取樣本。同時,為了區(qū)分對平面和曲面細(xì)節(jié)的提取能力,分別對樣本進(jìn)行選取;樣本容量維持45。
樣本數(shù)據(jù)均值匯總表如表1所示,分辨率等級包含低(L)、中(M)、高(H)和高等級開啟高倍鏡(H+)這四個等級。
表1 樣本數(shù)據(jù)均值匯總表Tab.1 Summary of the means of sample data mm
建立如表2所示的雙因素重復(fù)試驗(yàn)方差分析模型,分析不同分辨率等級、容積框邊長及平面和曲面對掃描數(shù)據(jù)分辨率是否有顯著影響。
表2 雙因素重復(fù)試驗(yàn)方差分析表Tab.2 Two-factor repeated tests ANOVA
根據(jù)雙因素重復(fù)試驗(yàn)方差分析模型,可得出以下假設(shè)。
① 假設(shè) H01:α1=α2=… =αr=0(αi稱為因子 A在水平Ai的效應(yīng))。
② 假設(shè)H02:β1=β2=… =βs=0(βj稱為因子 B 在水平Bj的效應(yīng))。
③ 假設(shè) H03:γij=0,i=1,2,…,r;j=1,2,…,s(γij稱為因子A、B在組合水平Ai×Bj的交互影響)。
若假設(shè)成立,則表明因子A和B或A和B的交互作用對試驗(yàn)結(jié)果無顯著影響;否則,因子對試驗(yàn)結(jié)果有顯著影響。
重復(fù)試驗(yàn)二元方差分析的F分布公式如下:
式中:QA為因子A引起的離差平方和;QB為因子B引起的離差平方和;QI為因子A、B交互作用引起的離差平方和;QE為誤差;r、s、c為常數(shù)。式(1)服從自由度為的F分布;式(2)服從自由度為的F分布;式(3)服從自由度為的F分布。一旦給定置信度水平α,查表即可得F分布的值。若,即拒絕H01,認(rèn)為因子A對試驗(yàn)結(jié)果有顯著影響;否則,接受H01,認(rèn)為因子A對試驗(yàn)結(jié)果無顯著影響。同理,若,即拒絕H02,否則接受 H02;若,即拒絕H03,否則接受 H03。
根據(jù)所建立的方差分析模型,其所涉及的三個因子包括掃描部位、容積框大小、分辨率等級。理論上,要分析這三個因子對分辨率的影響(包括單個或交互影響),應(yīng)該進(jìn)行三次組合,即控制一個因子,研究另外兩個因子。
設(shè)定顯著性水平α為5%,控制分辨率等級,并先將其設(shè)定為L等級,因子A表示不同的掃描部位,即平面和曲面兩種;因子B設(shè)定為在L等級下容積框的邊長,為100 mm、300 mm、500 mm、700 mm、900 mm、1 000 mm。根據(jù)建立的模型,可以計(jì)算出因子A的F分布FA=40.174 3和因子B的F分布FB=8 669.787 7,兩因子交互作用的F 分布FI=112.803 6,查 F 分布表得到 FA>F0.05(1 528)=3.84,F(xiàn)B>F0.05(5 528)=2.21,F(xiàn)I> F0.05(5 528)。故可認(rèn)定在L等級下,不同掃描部位、設(shè)定不同容積框邊長對分辨率的影響是顯著的。同理,可得其他三個分辨率等級F值,結(jié)果如表3所示。只有在M等級時掃描部位對分辨率影響不顯著,除此之外影響皆顯著。
表3 各分辨率等級下F值匯總表Tab.3 Summary of F-value at various resolution levels
對于類似可再設(shè)定容積框邊長討論掃描部位和分辨率等級對分辨率的影響,以及設(shè)定分辨率等級討論容積框邊長和掃描部位對分辨率的影響。為方便統(tǒng)計(jì),建議使用SPSS等統(tǒng)計(jì)軟件。
由以上分析可知,影響分辨率有三個內(nèi)部因素:掃描部位、設(shè)定的容積框邊長和設(shè)定的分辨率等級,且都對分辨率影響顯著。
分辨率線性模型表達(dá)式為:
式中:D為實(shí)際分辨率大小,mm;L為設(shè)定的容積框邊長,mm。
試驗(yàn)數(shù)據(jù)的分析結(jié)果表明,掃描不同部位對分辨率的影響是顯著的。
為了對分辨率有一個更直觀的了解,在這里采用加權(quán)平均法,即根據(jù)被掃描對象的特征估算曲面部位與平面部位所占的比例給予不同的權(quán)值。以加權(quán)后的結(jié)果擬合,既大大簡化模型、提高效率,又具有一定的可信度。標(biāo)準(zhǔn)塊表面曲面處約占總面積的20%、平面處約占80%,故可取加權(quán)值分別為0.2和0.8。擬合結(jié)果如表4所示。
表4 擬合參數(shù)匯總表Tab.4 Summary of the fitting parameters
由表4可知,該數(shù)據(jù)模型的擬合度均超過0.979,F(xiàn)值的顯著水平也都接近于0,回歸方程具有統(tǒng)計(jì)學(xué)意義。
試驗(yàn)統(tǒng)計(jì)獲得三維激光掃描儀Handyscan EXAscan的分辨率模型與廠家標(biāo)定的名義分辨率基本吻合:最高測試分辨率為0.049 7 mm,名義最高分辨率為0.05 mm;分辨率與掃描范圍間呈線性關(guān)系。此外,從實(shí)用角度看,H+級分辨率(最高)掃描時對計(jì)算機(jī)配置要求較高,難以實(shí)現(xiàn)便攜的要求。
通過掃描標(biāo)準(zhǔn)塊、獲取點(diǎn)云模型數(shù)據(jù)后,應(yīng)用Geomagic Qualify軟件處理、分析并建立模型,以判別三維掃描儀的分辨率水平,是一種切實(shí)可行的分辨率測試方法。該方法對于了解、認(rèn)證掃描儀的分辨率具有重要的指導(dǎo)意義。
[1]林家明,楊隆榮.CCD攝像機(jī)技術(shù)的發(fā)展趨勢及應(yīng)用前景[J].光學(xué)技術(shù),1999,11(6):43 -47.
[2]Boehler W,Marbs A.Investigating laser scanner accuracy[R].CIPA Symposium,Turkey,October 2003.
[3]Ahokas E,Yu Xiaowei,Oksanen J,et al.Optimization of the scanning angle for countrywide laser scanning[C]//ISPRS WGⅢ/3,Ⅲ/4,Ⅴ/3 Workshop“Laser scanning 2005”,Enschede,the Netherlands,2005.
[4]楊偉,劉春,劉大杰.激光掃描數(shù)據(jù)三維坐標(biāo)轉(zhuǎn)換的精度分析[J].工程勘察,2004(3):61 -63.
[5]朱凌,石若明.地面三維激光掃描點(diǎn)云分辨率研究[J].遙感學(xué)報(bào),2008,5(3):405 -410.
[6]徐進(jìn)軍,余明輝,鄭炎兵.地面三維激光掃描儀應(yīng)用綜述[J].工程勘察,2008(12):31 -34.
[7]Lichti D,Harvey B.The effects of reflecting surface material properties on time-of-flight laser scanner measurements[C]∥Symposium on Geospatial Theory,Processing and Applications,Ottawa,2002.
[8]蘇顯渝,李繼陶.信息光學(xué)[M].北京:科學(xué)出版社,2002:306 -316.
[9]杜立彬,高曉輝,夏進(jìn)軍,等.一種新型激光三維掃描測量儀的研制[J].激光技術(shù),2005,8(4):366 -369.
[10]唐朝偉,梁錫昌,鄒昌平.三維曲面激光精密測量技術(shù)[J].計(jì)量學(xué)報(bào),1994,4(2):99 -103.
[11]Li W X,Mitchell L D.Laser scanning system testing:errors and improvements[J].Measurement,1995,16(2):91 -101.
[12]王學(xué)軍,劉平方.三維曲面激光快速自動測量系統(tǒng)[J].激光雜志,1998(2):46 -50。