MEMS/NEMS表面3-D輪廓測量中基于模板的相位解包裹算法

2.2.2　實驗結(jié)果
圖6為圖2所示的圓形薄膜部分的一幅干涉圖像(左)和由五幅干涉圖像合成的亮場圖像(右).使用MATLAB中的edge函數(shù),可以很方便地得到使用不同邊緣檢測算法得到的邊緣提取結(jié)果,如圖7(閾值均為缺省設(shè)置).

圖6　圓形膜部分的干涉圖及亮場圖像

圖7　使用不同算法得到的邊緣檢測結(jié)果
從上圖結(jié)果可以看出,只有使用Canny算子的結(jié)果得到了非理想數(shù)據(jù)區(qū)完整閉合的邊緣.圖8是對Canny算法的閾值進行改變得到的結(jié)果.可見,當(dāng)閾值設(shè)定為(0.3,0.6)時可得到所需剔除區(qū)域的完整邊緣.將相應(yīng)的標(biāo)記有非相容點邊緣的模板代入到生長法程序中,得到解包裹結(jié)果如圖9.

圖8　更改Canny算法的閾值得到的結(jié)果

圖9　使用邊緣檢測法獲取模板得到的相位展開結(jié)果
使用邊緣檢測來獲取相位解包裹模板的方法使用方便,并且可以根據(jù)具體應(yīng)用選擇合適的閾值來進行靈活控制,比較適合應(yīng)用于具有孔洞等非相容區(qū)域的被測表面的相位解包裹運算中.然而,由于干涉測量中噪聲的影響,合成的亮場圖像本身已含有誤差,應(yīng)用于大范圍的比較復(fù)雜的表面時這種方法則較難在控制噪聲和精確定位邊緣之間找到平衡.

2.3　干涉圖灰度差值提取法
2.3.1　方　法
這種方法在原始干涉圖的基礎(chǔ)上得到相位解包裹的模板.在五步相移干涉測量中,驅(qū)動電路驅(qū)動壓電陶瓷使得參考鏡步進移動從而產(chǎn)生步進相移,由此得到時間序列上的五幅干涉圖像.參考光路的微小變化使得光程差發(fā)生相應(yīng)變化,致使各幅干涉圖上被測表面的條紋產(chǎn)生相對位移.而由于非數(shù)據(jù)區(qū)沒有條紋,或者非相容區(qū)域的條紋信息不夠理想,所以可以根據(jù)各幅干涉圖中對應(yīng)像素點灰度的差值來判斷哪些是相容區(qū)域,從而將非相容區(qū)域提取出來,得到相位解包裹的模板.具體實施的方式為:先將各幅干涉圖對應(yīng)像素的灰度兩兩相減,得到它們之間的灰度差值,然后設(shè)定一個合適的閾值,當(dāng)所有的灰度差值都小于這個閾值時,此像素點即被判定為沒有發(fā)生干涉條紋相對變化的非相容點.

2.3.2　實驗結(jié)果
如果想要對如圖10干涉圖所示的圓形薄膜進行全視場范圍的解包裹運算,由于其輪廓較為復(fù)雜,并包含有相當(dāng)部分的非數(shù)據(jù)區(qū)和非理想數(shù)據(jù)區(qū),所以使用干涉圖灰度差值提取的方法來獲得模板.設(shè)定合適的閾值得到的模板結(jié)果如圖11,圖中白色區(qū)域為較為理想的條紋區(qū)域,黑色區(qū)域為解包裹運算中需要繞過的區(qū)域.圖12是將圖11所示的模板以及包裹圖像代入

圖10　圓形膜五幅干涉圖中的一幅(1024x1528像素)

圖11　干涉圖灰度差值提取法獲得的解包裹模板圖

圖12　使用干涉圖灰度差提取法獲取模板得到的相位展開結(jié)果

生長算法得到的相位解包裹后并以圓形膜四周的基底區(qū)域作為基準面調(diào)平后的圓形膜表面3-D輪廓圖像.可以看出,使用干涉圖灰度差值提取法獲得模板并使用生長算法在全視場范圍內(nèi)進行相位展開獲得干涉圖灰度差值提取法的計算相對較為復(fù)雜,但是由于它在最原始的干涉圖像信息上進行判斷,所以成功率高,也可以通過改變閾值靈活地控制模板的靈敏度,并且使用范圍很廣,對于復(fù)雜輪廓、溝槽、噪聲等提取都適用.
3　結(jié)　語
具體應(yīng)用實例的實驗數(shù)據(jù)證明,基于模板的廣度優(yōu)先搜索相位解包裹方法可以根據(jù)不同應(yīng)用的需要標(biāo)記模板,從而繞過非相容區(qū)域準確地實現(xiàn)相位展開.如果有必要,可以根據(jù)被繞過區(qū)域周圍像素點的灰度信息,使用濾波、插值等方法回添這些點的展開相位數(shù)據(jù).此方法能夠克服普通相位展開方法的局限性,并因其簡便、靈活、準確的特點而能被廣泛應(yīng)用于EMS/NEMS結(jié)構(gòu)較為復(fù)雜的輪廓表面測量的相位展開.其不足在于不能應(yīng)用于非連續(xù)表面輪廓的測量.

致謝
感謝北京大學(xué)微米納米加工技術(shù)國家重點實驗室的陳兢副教授和王莎莎同學(xué)等,他們?yōu)楸疚牡难芯抗ぷ魈峁┝瞬糠譁y試結(jié)構(gòu).

參考文獻:
[1]　Denhoff M W. A Measurement of Young’s Modulus and Re-sidual Stress in MEMS Bridges Using a Surface Profiler[J].Journal of Micromechanics and Microengineering, 2003, 13:686-692.
[2]　Lafontan X, Pressecq F, Beaudoin F, et al. The Advent ofMEMS in Space[J]. Microelectronics Reliability, 2003, 43:1061-1083.
[3]　Jensen B D, de Boer M P, Masters N D, et al. Interferometryof Actuated Microcantilevers to Determine Material Propertiesand Test Structure Nonidealities in MEMS[J]. Journal of Mi-croelectromechanical Systems, 2001, 10(3): 336-346.

[4]　Conor O’Mahony, Martin Hill, Magali Brunet et al,“Charac-terization of Micromechanical Structures Using White-Light Interferometry[J]”, Measurement Science and Technology,2003,14: 1807?1814.

[5]　Hu Xiaodong, Hu Chunguang, Guo Tong, et al. Characteri-zing the Behaviour of Micro ElectroMechanical Structures by Optical Interferometry[C] // International Workshop on Mi-crofactories. Shanghai: 2004: 492-497.
[6]　胡春光.利用相移顯微干涉術(shù)和頻閃成像技術(shù)研究MEMS離面運動[J].天津:天津大學(xué),精密儀器與光電子工程學(xué)院,2004.
[7]　Dennis C. Ghiglia and Mark D. Pritt, Two-Dimensional PhaseUnwrapping: Theory, Algorithms and Software[M], John Wiley Sons, Inc., New York 1998: 1-58.
[8]　惠梅,表面微觀形貌測量中相移干涉術(shù)的算法與實驗研究
[D],博士學(xué)位論文,中國科學(xué)院西安光學(xué)精密機械研究所,2001年.
[9]　Hemmert W, Mermelstein M S and Freeman D M. Nanometer Resolution of Three-Dimensional Motions Using Video Inter-ference Microscopy[C]//Twelfth IEEE International Confer-ence on Micro Electro Mechanical Systems. 1999: 302-308.
[10]　李弼程,彭天強,彭　波等.智能圖像處理技術(shù)[M],電子工業(yè)出版社,2004:149-160.
[11]　薛軍平,何繼善.基于人體信息三維可視化的醫(yī)學(xué)圖像邊緣檢測算法研究[J].生物醫(yī)學(xué)工程研究. 2004, 23(1): 7-10.