二維數(shù)字圖像相關(guān)方法位移應(yīng)變測量中成像性因素的影響

二維數(shù)字圖像相關(guān)方法位移應(yīng)變測量中成像性因素的影響

1維數(shù)字圖像相關(guān)測量原理由于測量設(shè)備和測量過程簡單，測量模型極其廣泛（從模型到大變形、從微規(guī)模到宏觀規(guī)模、從傳統(tǒng)材料到建筑材料、從實驗室到室外和高溫氣候），整個測量過程的精度非常高。因此，光學(xué)鏡頭和單相機二維數(shù)字圖像的相關(guān)方法（2d-cic）已廣泛應(yīng)用于各種材料和結(jié)構(gòu)表面的位移和變形測量。在測量過程中,安裝了成像鏡頭的數(shù)字相機放置在被測物體的正前方,并盡可能使鏡頭光軸與被測平面試樣表面垂直,記錄不同狀態(tài)下試樣表面的散斑圖像。用基于灰度相似原理的模板匹配算法,可跟蹤“參考圖像”中各感興趣點在“變形后圖像”中亞像素精度的像面位移和全場應(yīng)變,所獲得的像面位移和應(yīng)變可進一步轉(zhuǎn)換為物體表面的真實變形。在二維數(shù)字圖像相關(guān)測量中絕大多數(shù)文獻都采用簡單的針孔成像模型描述像面位移和物面位移之間的對應(yīng)關(guān)系,即認(rèn)為利用數(shù)字圖像相關(guān)方法獲得的圖像位移(u,v)與真實位移(U,V)之間存在著簡單的線性對應(yīng)關(guān)系:u=MU,v=MV,這里M(等于針孔成像模型的像距L與物矩Z之比)為成像系統(tǒng)的放大倍數(shù)。雖然物面、像面位移間的這種簡單的線性比例關(guān)系已被廣泛地使用,但在實際測量中,由于各種不利因素的存在,這種理想的線性對應(yīng)關(guān)系并不嚴(yán)格成立或在變形后發(fā)生不同程度的改變。首先,由于實際被測平面與相機傳感器平面(像面)之間不會嚴(yán)格平行、加載方向也會與像平面存在一定偏差、被測材料會發(fā)生泊松效應(yīng)等各種因素的存在,試樣加載過程中其表面各點將不可避免地出現(xiàn)或多或少的離面位移。離面位移會改變針孔成像模型的物距Z,導(dǎo)致像面圖像放大倍數(shù)的改變,從而引起額外的壓縮或拉伸虛應(yīng)變。其次,相機通電后產(chǎn)生的自熱或測量過程中環(huán)境溫度的變化會使其內(nèi)部機械支撐部件產(chǎn)生熱變形,從而使相機傳感器的位置發(fā)生微小改變,進而改變針孔成像模型中的像距L。近期研究顯示一般商業(yè)相機自熱會在記錄圖像中引起約200με的虛應(yīng)變。此外,由于光學(xué)透鏡固有的像差、加工誤差和鏡片裝配誤差等原因,任何鏡頭都或多或少地存在畸變。鏡頭畸變會使實際像點偏離其理想位置,并引起額外的位移和應(yīng)變測量誤差。由于鏡頭畸變具有非均勻分布特性,因此由其引起的位移和應(yīng)變測量誤差也是非均勻分布的。使用普通鏡頭的二維數(shù)字圖像相關(guān)測量系統(tǒng)對物距、像距的微小改變非常敏感,某些鏡頭還存在較大的、不可忽略的畸變,因此會引起較大的位移和應(yīng)變測量誤差。客觀地說,要準(zhǔn)確測量材料在彈性范圍內(nèi)的小變形,以上這些不利因素所引起的測量誤差往往是不可忽略的。如何避免或消除以上3種不利因素引起的測量誤差以獲得高精度的變形測量結(jié)果是數(shù)字圖像相關(guān)方法研究人員和使用者最為關(guān)心的問題?；谶@一目的,本文分析了二維數(shù)字圖像相關(guān)測量中由被測物體的離面位移、相機自熱和鏡頭畸變引起的測量誤差。使用高質(zhì)量的雙遠(yuǎn)心鏡頭建立了高精度二維數(shù)字圖像相關(guān)測量系統(tǒng),并通過實際實驗仔細(xì)研究了使用3種不同鏡頭(標(biāo)像鏡頭、物方遠(yuǎn)心鏡頭和雙遠(yuǎn)心鏡頭)的二維數(shù)字圖像相關(guān)測量系統(tǒng)在3種不利因素出現(xiàn)情況下的位移和應(yīng)變測量結(jié)果。實驗結(jié)果顯示高質(zhì)量的雙遠(yuǎn)心鏡頭不僅對被測物體表面的離面位移和相機自熱不敏感,并且鏡頭畸變也非常小,因此是高精度二維數(shù)字圖像相關(guān)測量中必不可少的光學(xué)元件。2相關(guān)方法和二維數(shù)字圖像的誤差分析2.1維數(shù)字圖像相關(guān)測量系統(tǒng)的設(shè)計圖1為二維數(shù)字圖像相關(guān)方法測量系統(tǒng)示意圖以及一幅典型的散斑圖。被測平面物體表面需具有隨機的灰度分布(通常稱為散斑場),該散斑場作為變形信息載體隨試樣表面一起變形。實驗過程中,不同加載條件下被測物體表面的散斑圖由放置在被測物體正前方的數(shù)字相機采集并存入計算機,隨后用數(shù)字圖像相關(guān)方法分析可提取像面位移和應(yīng)變信息,所獲得的像面變形可進一步轉(zhuǎn)換為物體表面的真實變形。成像鏡頭是二維數(shù)字圖像相關(guān)測量系統(tǒng)中最重要的組成部分,其作用是將變形前后被測平面物體表面的散斑圖案成像到相機傳感器平面上以獲得數(shù)字圖像。有關(guān)二維數(shù)字圖像相關(guān)方法的文獻多用簡單的針孔成像模型定量地描述像點坐標(biāo)(x,y)和物點坐標(biāo)(X,Y)的對應(yīng)關(guān)系,認(rèn)為變形前后的物距和像距為常數(shù),而且不考慮成像鏡頭的畸變。在這種理想條件下,由數(shù)字圖像相關(guān)方法獲得的像面位移[u(x,y),v(x,y)]顯然與物體表面的位移[U(X,Y),V(X,Y)]呈線性比例關(guān)系,即u(x,y)=MU(X,Y),v(x,y)=MV(X,Y)。2.2維和二維數(shù)字圖像相關(guān)測量結(jié)果的基本過程分析盡管在眾多文獻中,物面、像面位移的線性比例關(guān)系被廣泛地使用,但在實際的位移、應(yīng)變測量中,由于以下各種不利因素的存在,這種理想的線性對應(yīng)關(guān)系并不嚴(yán)格成立,在某些情況下會引起不可忽略的較大測量誤差。首先,由于實際被測物面和相機傳感器平面(像面)之間不會嚴(yán)格平行、加載方向也與像面存在不同程度的偏差、被測材料會產(chǎn)生泊松效應(yīng)等各種因素的存在,試樣加載過程中其表面各點將不可避免地出現(xiàn)或多或少的離面位移。離面位移會改變物距,如果物距減小,會使圖像放大;反之,圖像縮小。因此,離面位移即意味著圖像中會出現(xiàn)拉伸或壓縮正應(yīng)變。圖2顯示了物距為Z的平面試樣出現(xiàn)了剛體離面位移ΔZ(假設(shè)物體移向成像鏡頭的離面位移ΔZ為正),可見高度為Y的物點在像平面中的高度由y變成y′。由簡單的推導(dǎo)可知,對坐標(biāo)(X,Y)的物點,剛體離面位移引起的面內(nèi)位移和正應(yīng)變分別為顯然,當(dāng)物體表面靠近相機時,剛體離面位移將帶來雙向拉伸虛應(yīng)變;當(dāng)物體遠(yuǎn)離相機時,剛體離面位移將帶來雙向壓縮虛應(yīng)變,并且虛應(yīng)變在數(shù)值上隨著離面位移量ΔZ的增加而增大。如果被測物體產(chǎn)生了非均勻的離面位移(如加載過程中出現(xiàn)了離面轉(zhuǎn)動或試樣局部出現(xiàn)了頸縮),則在像面產(chǎn)生非均勻分布的虛應(yīng)變。為減小離面位移對二維數(shù)字圖像相關(guān)方法的影響,可使成像系統(tǒng)盡可能遠(yuǎn)離被測物體,即通過增加物距Z來減小離面位移ΔZ的影響。其次,相機作為一種電子器件,在通電后其內(nèi)部的電子元器件會因電阻發(fā)熱效應(yīng)產(chǎn)生熱量,使相機內(nèi)部溫度升高,即所謂的“相機自熱”。最近Ma等對6種不同型號相機的測試顯示,相機自熱會在通電后的1.5h后產(chǎn)生約10℃的溫升,隨后進入較穩(wěn)定的熱平衡階段。相機自熱(或測量過程中環(huán)境溫度的變化)會使相機內(nèi)部的支撐機械元件產(chǎn)生熱變形(如圖3所示),導(dǎo)致相機傳感器位置(像距L)出現(xiàn)微小改變,從而改變原來的線性對應(yīng)關(guān)系,在像平面引起100~250με的虛應(yīng)變。為減小相機自熱的影響,文獻中建議在相機達(dá)到熱平衡后再使用。最后,由于光學(xué)透鏡固有的像差、加工誤差和鏡片裝配誤差等原因,任何成像鏡頭都或多或少地存在畸變。鏡頭畸變會使實際像點偏離其理想位置,并引起額外的位移測量誤差。由于應(yīng)變場是通過對位移場差分獲得,因此由畸變引起的應(yīng)變誤差將更為明顯。圖4顯示鏡頭畸變的存在對數(shù)字圖像相關(guān)方法位移測量結(jié)果的影響,圖中(Cx,Cy)為假想的畸變中心;(xu,yu)、(xd,yd)分別為無畸變和畸變的像點坐標(biāo);δr為徑向畸變,δt為切向畸變。圖4中的插圖顯示,由于徑向畸變的存在,實際測量獲得的位移矢量dm與真實位移矢量da存在一定的偏差。文獻顯示,徑向畸變在參考圖像中某點引起的位移測量誤差(Δu,Δv)與徑向畸變系數(shù)k1、該點的圖像坐標(biāo)(xd,yd)和該點由數(shù)字圖像相關(guān)方法直接計算獲得的畸變位移(u′,v′)有關(guān):推導(dǎo)顯示,由徑向畸變引起的應(yīng)變測量誤差(Δεx,Δεy,Δγxy)也與以上3個參數(shù)有關(guān),即為了保證數(shù)字圖像相關(guān)測量結(jié)果的準(zhǔn)確性,必須要研究成像鏡頭的畸變程度。不同研究人員也提出不同的畸變校正方法,有興趣的讀者可參考文獻[6~10]。3圖像大小的變化在機器視覺領(lǐng)域,特殊設(shè)計的遠(yuǎn)心鏡頭常用來避免傳統(tǒng)鏡頭的透視畸變。在物方遠(yuǎn)心鏡頭中,孔徑光闌被放置在鏡頭像方的焦平面上,僅使與光軸平行的光線通過焦平面上的孔徑光闌在相機靶面成像[見圖5(a)],因此所有的光線可看作來自無窮遠(yuǎn)處。這種特殊的平行光路設(shè)計使其可以在一定的物距(遠(yuǎn)心深度)范圍內(nèi),使得到的圖像放大倍率不會隨物距的變化而變化。在圖5(b)像方遠(yuǎn)心鏡頭中,孔徑光闌被放置在鏡頭物方焦平面上,使像方主光線平行于光軸。因此像方遠(yuǎn)心鏡頭對相機圖像傳感器(像平面)位置的微小變化不敏感,即像距的改變不會影響圖像的大小。圖5(c)所示的雙遠(yuǎn)心鏡頭綜合了物方遠(yuǎn)心鏡頭和像方遠(yuǎn)心鏡頭的優(yōu)點,對在其遠(yuǎn)心深度內(nèi)變化的物距和像距變化均不敏感。將雙遠(yuǎn)心鏡頭用在二維數(shù)字圖像相關(guān)測量系統(tǒng)中,其對物距和像距變化均不敏感的優(yōu)點將有助于在被測物體變形前后獲得放大倍數(shù)恒定的圖像,對因加載出現(xiàn)的離面位移和因相機自熱或環(huán)境溫度變化引起的靶面位置的微小移動不敏感,從而有助于實現(xiàn)高精度的面內(nèi)位移和應(yīng)變測量。4實驗與研究4.1維數(shù)字圖像為了定量研究使用雙遠(yuǎn)心鏡頭的二維數(shù)字圖像相關(guān)測量系統(tǒng)在以上3種不利因素存在情況下的測量精度并與其他鏡頭比較,本實驗中使用了如下3種具有代表性的鏡頭:1)日本CBC集團公司生產(chǎn)的焦距為55mm的定焦鏡頭(型號:Tec-55);2)北京大恒光電生產(chǎn)的物方遠(yuǎn)心鏡頭(型號:GCO-2301+GCO-235);3)德國施耐德光學(xué)公司生產(chǎn)的雙遠(yuǎn)心鏡頭(型號:Xenoplan1…5)。這3種鏡頭分別與空間分辨率為1280pixel×1024pixel的CMOS數(shù)字?jǐn)z像機(DH-HV1351UM,大恒圖像)組成3種不同的二維數(shù)字圖像相關(guān)測量系統(tǒng)。這3種不同鏡頭如圖6所示,表1給出了3種鏡頭在實驗中的實際工作距離、遠(yuǎn)心深度以及放大倍數(shù)。圖7是使用施耐德雙遠(yuǎn)心鏡頭的二維數(shù)字圖像相關(guān)測量系統(tǒng)照片。實驗中的被測平板玻璃試件表面制作了隨機散斑,該平板玻璃固定在一個定位精度為5μm的精密三維平移臺上,可沿著水平(X軸)、豎直(Y軸)和光軸(Z軸)3個方向準(zhǔn)確平移。4.2圖像采集和測點布置為全面準(zhǔn)確地分析使用不同鏡頭的二維數(shù)字圖像相關(guān)測量系統(tǒng)的測量精度,設(shè)計了3個實驗以重現(xiàn)上述3種不利因素結(jié)果的影響。1)離面剛體位移實驗。在該實驗中,平板玻璃試樣放置在各種鏡頭的工作距離處清晰成像,并保存一幅參考圖像。隨后用精密平移臺使試樣沿Z軸正(靠近鏡頭)、負(fù)(遠(yuǎn)離鏡頭)方向每次平移1mm后保存1幅圖像,最大離面剛體位移為5mm,共保存10幅沿Z軸剛體平移的圖像。為了盡可能地減小相機自熱對數(shù)字圖像相關(guān)方法測量的位移場和應(yīng)變場的影響,根據(jù)文獻的建議對相機事先預(yù)熱了2h后才進行圖像采集。2)相機自熱實驗。在該實驗中,平板玻璃試樣始終牢靠地固定在平移臺上并保存靜止,相機通電后首先采集一幅圖像作為參考圖像,隨后每隔2min保存1幅圖像,圖像采集持續(xù)5h,共保存150幅不同時間的靜態(tài)圖像。為盡量避免環(huán)境溫度變化的影響,實驗在環(huán)境溫度變化極小的密閉室內(nèi)進行。因為試樣保持靜止?fàn)顟B(tài),所以由數(shù)字圖像相關(guān)方法獲得的圖像位移只可能是由相機自熱引起的傳感器靶面移動導(dǎo)致的。3)面內(nèi)平移實驗。該實驗用以定量分析各鏡頭的畸變造成的測量誤差。平板玻璃試樣放置在各種鏡頭的工作距離處清晰成像,并保存1幅參考圖像。隨后用精密平移臺使試樣沿X軸正、負(fù)方向每次平移0.5mm后保存1幅剛體位移后的圖像,最大面內(nèi)剛體位移為2mm,共保存8幅平移后的圖像。同樣,相機也是事先預(yù)熱了2h后再進行實驗。4.3計算方法的選取以上實驗所采集的變形前后的數(shù)字圖像用自編的數(shù)字圖像相關(guān)分析軟件進行處理可獲得像面位移和應(yīng)變。具體計算時,選擇位于參考圖像中間大小為820pixel×680pixel的矩形區(qū)域作為計算區(qū)域,計算位移場時所用圖像子區(qū)大小為41pixel×41pixel,步長為10pixel,共計算了5727個計算點。位移計算采用了作者最近提出的快速、高精度Newton-Raphson算法。計算應(yīng)變場則采用作者之前提出的逐點最小二乘應(yīng)變估計方法,應(yīng)變計算窗口大小為21pixel×21pixel。此外,在計算整個計算區(qū)域的平均應(yīng)變時,采用線性平面擬合的方法,即對兩個位移場分別用雙線性函數(shù)進行最小二乘擬合,擬合函數(shù)的系數(shù)可直接作為計算區(qū)域的平均應(yīng)變。5結(jié)果5.1大恒光電的物方遠(yuǎn)心鏡頭的拉伸和壓縮虛應(yīng)變圖8是平板玻璃沿Z軸正方向平移5mm后由3種不同類型鏡頭拍攝獲得的水平和豎直方向的位移場。從圖8(a)中等間距均勻分布的位移等值線可以看出普通鏡頭對離面剛體位移非常敏感,離面位移引起的像面位移場呈明顯的雙向“拉伸”狀態(tài),這與第3節(jié)分析吻合。從圖8(b)對于大恒光電的物方遠(yuǎn)心鏡頭位移場也能看出雙向拉伸狀態(tài),但因拉伸位移較小,圖像噪聲的影響較為明顯。這表明離面位移對大恒光電的物方遠(yuǎn)心鏡頭影響相對較小。而施耐德雙遠(yuǎn)心鏡頭采集的圖像幾乎完全不受離面位移的影響,圖8(c)中的位移場不僅數(shù)值較小,而且看不出任何規(guī)律性的分布。圖9是使用不同鏡頭的二維數(shù)字圖像相關(guān)測量系統(tǒng)獲得的不同離面剛體位移時的兩個方向(X,Y方向)面內(nèi)平均虛應(yīng)變。對該圖中的離散數(shù)據(jù)進行線性擬合,擬合直線的斜率可定量地表征各鏡頭對離面剛體位移的敏感程度?？梢钥闯?普通鏡頭對離面位移敏感,單位毫米離面位移約引起2579με的拉伸或壓縮虛應(yīng)變。大恒光電的物方遠(yuǎn)心鏡頭,單位毫米離面位移引起360~380με的拉伸或壓縮虛應(yīng)變。而離面位移對施耐德雙遠(yuǎn)心鏡頭幾乎沒有任何影響,單位毫米的離面位移引起最大僅為4με的拉伸或壓縮虛應(yīng)變。5.2虛應(yīng)變與時間的關(guān)系由不同鏡頭所記錄的序列圖像與初始時刻的參考圖像用4.3節(jié)介紹的計算方法進行分析可獲得由相機自熱引起的X、Y方向位移場和平均應(yīng)變。圖10繪出了相機自熱引起的虛應(yīng)變與時間的關(guān)系曲線,該圖顯示普通鏡頭、大恒光電物方遠(yuǎn)心鏡頭對相機自熱引起的靶面位移都很敏感,在相機通電后的1.5~2h內(nèi),由相機自熱會在水平和豎直方向引起拉伸正應(yīng)變,但對切應(yīng)變無影響(切應(yīng)變在零值附近波動)。拉伸正應(yīng)變隨時間線性增加,隨后進入熱平衡階段,最大熱應(yīng)變穩(wěn)定200~250με之間。相比之下,施耐德雙遠(yuǎn)心鏡頭則對相機自熱也不敏感,相機自熱會在水平和豎直方向引起最大僅為25με的壓縮正應(yīng)變。5.3其他性別角色圖像的應(yīng)變分布圖v場,如圖11和圖12所示。在理想情況下,面內(nèi)剛體平移試樣表面各點的位移分量應(yīng)該為常數(shù),位移場應(yīng)該為一平面。然而,由于非線性鏡頭畸變和圖使用三種鏡頭的二維數(shù)字圖像相關(guān)測量系統(tǒng)獲得試樣沿水平方向平移1mm后的畸變位移u場和像噪聲的影響,所有的位移場均為非均勻分布。其中大恒光電物方遠(yuǎn)心鏡頭所獲取圖像的水平和豎直方向位移場分布呈現(xiàn)對稱的二次拋物面和反對稱的馬鞍形分布,這種規(guī)律性的位移分布可用文獻中介紹的徑向畸變模型很好地解釋。對于Computar鏡頭,也呈現(xiàn)類似的規(guī)律性分布,只是徑向畸變相對較小,圖像噪聲的影響更為明顯。使用施耐德雙遠(yuǎn)心鏡頭的測量系統(tǒng)所測量的位移場規(guī)律性不明顯,其兩個方向位移的標(biāo)準(zhǔn)差均小于0.02pixel,顯示該雙遠(yuǎn)心鏡頭畸變極小。用逐點最小二乘法對圖11和圖12中的位移場進行差分處理,進一步獲得由不同鏡頭所得圖像的應(yīng)變分布圖,如圖13