基于malab圖像處理的振動(dòng)位移測(cè)試

基于malab圖像處理的振動(dòng)位移測(cè)試

在實(shí)驗(yàn)室進(jìn)行的各種測(cè)試中，許多機(jī)器輸出的結(jié)果只是一個(gè)曲線，而沒(méi)有曲線中點(diǎn)的坐標(biāo)數(shù)據(jù)。例如:使用測(cè)振儀測(cè)量構(gòu)件振動(dòng)位移的測(cè)試結(jié)果就是一條繪制在坐標(biāo)紙上的曲線。在隨后的數(shù)據(jù)處理過(guò)程中,通常的做法是手工標(biāo)出曲線上的若干點(diǎn),然后根據(jù)坐標(biāo)紙讀出這些點(diǎn)的坐標(biāo)值,作為實(shí)驗(yàn)得到的初始數(shù)據(jù),以此為基礎(chǔ)進(jìn)行計(jì)算和分析。這種處理方法不僅繁瑣、工作量大,而且在曲線上取點(diǎn)的多少和坐標(biāo)值的判讀都因人而異,給測(cè)試結(jié)果帶來(lái)了較大的誤差。作為一種應(yīng)用廣泛的編程工具,MATLAB在圖形處理方面有著明顯的優(yōu)勢(shì):具有強(qiáng)大的矩陣運(yùn)算功能,在進(jìn)行圖像處理時(shí)可以避免許多繁瑣的運(yùn)算;圖形顯示方便,可以在調(diào)試過(guò)程中隨時(shí)觀察圖形的變化;帶有豐富的圖像處理函數(shù)庫(kù),其圖像處理工具箱(imageprocessingtoolbox)幾乎涵蓋了所有常用的圖像處理函數(shù)。筆者利用MATLAB的圖像處理工具箱,對(duì)實(shí)驗(yàn)輸出的圖形進(jìn)行處理,獲得實(shí)驗(yàn)圖線的二進(jìn)制圖像矩陣,進(jìn)而通過(guò)對(duì)二進(jìn)制矩陣的分析獲得實(shí)驗(yàn)圖線的坐標(biāo)值。實(shí)踐證明,該過(guò)程簡(jiǎn)單實(shí)用,獲取坐標(biāo)值精度高,是一種解決上述問(wèn)題的有效方法。1原始輸出流預(yù)覽的處理1.1基于原始圖像的生成實(shí)驗(yàn)中輸出的圖形,一般可通過(guò)掃描儀掃入計(jì)算機(jī),使其成為可進(jìn)一步編輯的計(jì)算機(jī)圖像,見(jiàn)圖1。在MATLAB中,通過(guò)imread函數(shù)可以將圖形讀入內(nèi)存,返回一個(gè)大小與原圖形大小(圖像的大小以像素計(jì))相同的矩陣,輸入語(yǔ)句:通常情況下,獲得的原始圖像是一個(gè)RGB圖像,由于涉及的圖像處理函數(shù)多針對(duì)灰度圖像,因此需要將RGB圖像ori轉(zhuǎn)換為灰度圖像bw,輸入語(yǔ)句:灰度圖像矩陣中每個(gè)元素代表一個(gè)像素點(diǎn),如果矩陣為unit8類型,其數(shù)據(jù)范圍為,不同的元素?cái)?shù)值代表不同的亮度或灰度級(jí),數(shù)值0代表黑色,數(shù)值255代表白色,數(shù)值越大,表明該位置像素點(diǎn)亮度越高。通過(guò)imshow函數(shù),可以將矩陣代表的圖形顯示在圖形窗口上,可以看到,顯示出來(lái)的圖形與原圖形是完全一樣的。1.2圖像的二維特征為了消除圖形中無(wú)關(guān)像素的干擾,需要除去圖形的背景,僅僅保留輸出結(jié)果圖線。將圖形轉(zhuǎn)換為二進(jìn)制圖形,再施以相應(yīng)的變換,就可以達(dá)到這個(gè)目標(biāo),滿足數(shù)字化的需要。筆者首先利用MATLAB中雙閾值選擇的方法,有效地清除圖形中的次要區(qū)域。用小波分析的方法處理掃描圖形,根據(jù)圖形掃描強(qiáng)度和明暗的變化,可以確定2個(gè)閾值,用這2個(gè)閾值對(duì)原圖形矩陣進(jìn)行篩選,如閾值threshold1=100,threshold2=210,輸入語(yǔ)句:篩選后返回的bw1和bw2矩陣是2個(gè)二進(jìn)制矩陣,與原圖形矩陣大小相同,其中的元素僅有0和1組成。數(shù)字0代表原圖形矩陣該元素所在位置像素小于等于設(shè)定的閾值(即“否”),而數(shù)字1代表大于設(shè)定的閾值(即“是”)。與灰度圖像相同,在顯示為圖形時(shí),數(shù)值0代表黑色,數(shù)值1代表白色,因此二進(jìn)制圖像可以看作是一種特殊的灰度圖像,見(jiàn)圖2。由于閾值設(shè)定的不同,bw1中黑色像素點(diǎn)的數(shù)量遠(yuǎn)遠(yuǎn)少于bw2中的數(shù)量,人們可以選擇合適的閾值,使bw1中包含的全是結(jié)果曲線上的像素點(diǎn),再以bw1為基礎(chǔ)對(duì)bw2進(jìn)行選擇,可以很好地清除圖形中實(shí)驗(yàn)曲線之外的背景像素,如圖形記錄紙上的網(wǎng)格線等均可被清除。首先確定bw1中像素點(diǎn)的位置,輸入語(yǔ)句:r和c是2個(gè)數(shù)組,分別代表bw1矩陣中0元素所在的行值和列值,也就是黑像素點(diǎn)所在位置。對(duì)bw2進(jìn)行選擇,輸入語(yǔ)句:函數(shù)bwselect的作用是尋找bw2中與指定位置(c,r)相重疊的物體。需要指出的是,MATLAB二進(jìn)制圖像中的“物體”,指的是白色像素點(diǎn)組成的區(qū)域,因此函數(shù)bwselect中的選擇源應(yīng)該是bw2圖形的“負(fù)片”,即矩陣中的元素0和1互換。如前所述,(c,r)指示的都是曲線上點(diǎn)的位置,因此返回的bwclean就是原圖形中的結(jié)果曲線,見(jiàn)圖3。通過(guò)imshow(bwclean)可以看到,獲得的bwclean曲線還并不完善,曲線較粗且不光滑,需要對(duì)它作進(jìn)一步的處理,才能滿足數(shù)字化的需要。MATLAB提供了一個(gè)強(qiáng)大的二進(jìn)制圖像處理函數(shù)bwmorph,其用法為bwmorph(pic,operation,n),通過(guò)指定“operation”為不同的參數(shù),可以完成對(duì)二進(jìn)制圖像不同的處理。這里首先使用“skel”參數(shù),將曲線進(jìn)一步細(xì)化,使像素點(diǎn)的位置更加精確,見(jiàn)圖4。從圖4可以看出返回的bwskel曲線已經(jīng)很細(xì),只是仍不夠光滑,一些微小的“毛刺”需要進(jìn)一步去除。再次使用bwmorph函數(shù),選用“spur”參數(shù),將bwskel曲線上的“毛刺”去除。觀察返回的圖形,曲線已經(jīng)非常精確,圖形中不包含無(wú)關(guān)的像素點(diǎn),見(jiàn)圖5。2圖形的數(shù)字化2.1bwpruned矩陣大小經(jīng)過(guò)圖像處理,最終得到的bwpruned矩陣為二進(jìn)制格式,矩陣中元素為0和1,原圖片中每一個(gè)像素點(diǎn)對(duì)應(yīng)矩陣中的一個(gè)元素。由于處理后的圖形已經(jīng)沒(méi)有無(wú)關(guān)像素點(diǎn),因此矩陣中0(顯示為黑色像素)所在位置就是原圖片中曲線所在位置,而圖片的空白處在矩陣中則全部顯示為1。bwpruned矩陣的大小與原圖片的大小總是相同的,所以矩陣中某一個(gè)0所在相對(duì)位置和其對(duì)應(yīng)的曲線上像素點(diǎn)在原圖片中的相對(duì)位置也是相同的。確定了矩陣中元素0的位置,就相當(dāng)于確定了曲線在原圖片中的位置。將確定的相對(duì)位置在圖片的真實(shí)坐標(biāo)中定位,即可知道曲線上每個(gè)點(diǎn)的坐標(biāo)值。例如:原曲線上一點(diǎn)a,其對(duì)應(yīng)的數(shù)值0出現(xiàn)在二進(jìn)制矩陣的i行j列,原圖形的大小(即矩陣大小)為m×n,而原圖形對(duì)應(yīng)的真實(shí)坐標(biāo)尺度為x×y,則可求出a點(diǎn)在真實(shí)坐標(biāo)中的位置為x·(j-1)/(n-1),y·(i-1)/(m-1),再根據(jù)真實(shí)坐標(biāo)的起點(diǎn)進(jìn)行坐標(biāo)平移,即可得到a點(diǎn)的坐標(biāo)值。2.2tempy位置進(jìn)行轉(zhuǎn)換在MATLAB中,判斷一特定元素在矩陣中的位置是非常簡(jiǎn)單的,完全不必像其他語(yǔ)言一樣使用循環(huán)嵌套判斷,只需使用一個(gè)find命令即可,輸入語(yǔ)句:返回的tempy和tempx是2個(gè)數(shù)組,分別記錄了矩陣中每個(gè)0元素所在的行與列,由于矩陣的行數(shù)計(jì)算是自上而下的,而實(shí)際的縱坐標(biāo)通常是自下而上的,輸入下列語(yǔ)句將tempy轉(zhuǎn)換為相對(duì)于矩陣最后一行的相對(duì)值,可表示曲線的縱向相對(duì)位置:其中matixheight為矩陣的總行數(shù)。將獲得的橫向位置tempx和縱向位置tempy與圖片的真實(shí)坐標(biāo)比較,進(jìn)行相應(yīng)的縮放和平移,即可獲得曲線上每點(diǎn)的坐標(biāo)值,見(jiàn)表1。比較圖6與原始曲線圖1可以看出,通過(guò)獲得的坐標(biāo)值重新繪制的圖線與原掃描圖線符合得非常好,應(yīng)用坐標(biāo)值計(jì)算的