數(shù)字圖像匹配

1、數(shù)字圖像匹配一.需求分析題目:數(shù)字圖像匹配問題描述:所謂圖像匹配，就是指圖像之間的比較、得到不同圖像之間的相似度?；跀?shù)字圖像，編寫對兩副數(shù)字圖像進行匹配的算法及演示程序?；疽筮M行匹配的兩幅圖像為格式或格式。能夠進行對兩幅數(shù)字圖像的匹配。（3采）用.交互式程序對圖像進行匹配。提高要求:能夠對數(shù)字進行簡單處理，例如放大，縮小，翻轉，灰度處理，圖象二值化。開發(fā)環(huán)境：自帶的界面編輯器基本功能:通過分析題目的基本要求，我將此軟件的基本功能主要分為2大模塊:一個是數(shù)字圖像處理模塊，另一個是數(shù)字圖像匹配模塊。在數(shù)字圖像處理模塊中，用戶可以對數(shù)字圖像進行簡單的處理，可以對圖像進行放大，縮小，翻轉，灰度處

2、理。在數(shù)字圖像匹配模塊中，用戶可以對兩張圖像進行匹配并顯示匹配結果。性能、接口：輸入輸出形式此軟件以編輯器開發(fā)出的界面作為載體對相映的圖像行相應的操作，所以輸入輸出形式主要是通過編輯器開發(fā)出的界面來實現(xiàn)的。輸入形式：輸入任何一幅格式或格式的數(shù)字圖像。輸出形式：將經(jīng)過相應操作處理后的圖片顯示出來。測試數(shù)據(jù)要求任何一幅格式或格式的數(shù)字圖像。二算法設計1.概念解釋：數(shù)字圖像數(shù)字圖像是由被稱做像素的小塊區(qū)域組成的二維像素矩陣。一般把圖像分成3種形式：單色圖像，灰度圖像和彩色圖像。像素：表示圖像顏色的最小單位灰度圖像：灰度圖是指只含亮度信息，不含色彩信息的圖像，就像平時看到的黑白照片：亮度由暗到明，變化

3、是連續(xù)的?；叶葓D的每個像素的亮度用一個數(shù)值來表示，通常數(shù)值范圍在025之5間，即可用一個字節(jié)來表示，0表示黑，25表示白，而其他表示灰度。點陣圖:顯示器的屏幕由可以發(fā)光的像素點組成.并且從幾何位置看,所用這些像素點構成一個矩形的陣列.利用計算機控制各像素點按我們指定的要求發(fā)光,就構成了我們需要的圖形.這種方式構成的圖形我們可稱之為點陣圖形.點陣圖形的坐標系統(tǒng)：各像素點有一個坐標唯一指定了它的位置如果點陣圖形的大小是義那么它的點陣共有行列每個像素點的位置就由它所在的行和列的位置所唯一確定.這個行和列的位置就給出了點陣圖形的坐標系統(tǒng).按照前面的順序第行第列的像素點順序數(shù)就是反之順序數(shù)為的像素點在第

4、行第列這里的是除以后的余數(shù)是的整數(shù)部分需要注意的是第行第列的像素點的坐標可能不是而是這是因為有時為了在計算機中處理的方便像素點的行列的排序不是從1,而是從0開始的.我們常用的顯示器的像素坐標就是如此.2數(shù)字圖像匹配算法設計：在此軟件中我采用了兩種圖像匹配算法：基于灰度的模板匹配算法基于灰度的快速匹配算法。由于各種各樣的原因如（成象條件的差異）圖象預處理，引入的誤差等，參與圖象匹配的模板與潛在的匹配子圖象間通常存在著程度不同的不一致，因此根據(jù)模板在一幅陌生圖象中檢測出潛在的匹配對象并得出它在圖象中的位置是一件復雜的工作?；诨叶鹊拿嫫ヅ渌惴０迤ヅ涫侵赣靡粋€較小的圖像，即模板與源圖像進行比較，

5、以確定在源圖像中是否存在與該模板相同或相似的區(qū)域，若該區(qū)域存在，還可確定其位置并提取該區(qū)域。模板匹配常用的一種測度為模手術臺與源圖像對應區(qū)域的誤差平方和。設為義的源圖像，為義WW的模板圖像則誤差平方和測度定義為Dx,y)J1Kf(xj,y)(j,k)Rj10k10由上式展開可得:D(X,y)K11f(Xj,yk)22.t(j,k)If(xj,yk)t(j,k)2jIDkIDjkjIDkID令DS(x,y)f(xj,yk)RjkDST(x,y)2.t(j,k)胃(xj,yk)jkDT(x,y)t(j,k)2jknDS(x,y)稱為源圖像中與模板對應區(qū)域的能量，它與像素位置(x,y)有關，但隨像素

6、位置(x,y)的變化，DS(x,y)變化緩慢。DST(x,y)模板與源圖像對應區(qū)域的互相關，它隨像素位置(x,y)的變化而變化，當模板t(j,k)和源圖像中對應區(qū)域相匹配時取最大值。DT(x,y)稱為模板的能量，它與圖像像素位置(x,y)無關，只用一次計算便可。顯然，計算誤差平方和測度可以減少計算量。基于上述分析，若設DS(x,y)也為常數(shù)，則用DST(x,y)便可進行圖像匹配，當DST(x,y)取最大值時，便可認為模板與圖像是匹配的。但假設DS(x,y)為常數(shù)會產(chǎn)生誤差，嚴重時將無法下確匹配，因此可用歸一化互相關作為誤差平方和測度，其定義為：t(j,k)If(xj,yk)R(x,y)./.k

7、卬_產(chǎn)f(xj,yk)t(j,k)2j*k*j*k*下圖給出了模板匹配的示意圖，其中假設源圖像f(x,y)和模板圖像t(k,l)的原點都在左上角。對任何一個f(x,y)中的(x,y),根據(jù)上式都可以算得一個R(x,y).當x和y變化時，t(j,k)在源圖像區(qū)域中移動并得出R(x,y)所有值。R(x,y)的最大值指出了與t(j,k)匹配的最佳位置，若從該位置開始在源圖像中取出與模板大小相同的一個區(qū)域，便可得到匹配圖像?；诨叶鹊目焖倨ヅ渌惴?局部灰度特征的編碼與計算首先將整幅圖像劃分為X尺寸且互不重疊的方塊，可根據(jù)問題任意選擇,稱該方塊為塊.如果圖像的邊長不是的整數(shù)倍，則將最底部與最右邊剩余的幾

8、行、幾列裁剪掉下文將說明這并不影響最終的匹配結果對邊長為的圖像，共可得到個塊.對于塊，表示所包含像素的灰度值之和.定義1.塊如圖1中的所示與其周圍個相鄰的塊如圖1中的,所示組成塊的鄰域.將塊的鄰域分為個部分，分別為12346789,如圖1所示稱為塊的鄰域.塊分別屬于個鄰域，即1234511uUUUUUUUUUUU.245245783568942356對于每個鄰域中的個塊,可規(guī)定一個順序如圖中所取的逆時針序對所包含的個塊的像素灰度值之和，RR做排序，顯然jj1j2j3j4共有種可能，每種排序結果可以用位的二進制編碼來表示，記作0000,00001，10111.圖將塊所在的個塊的做位串拼接得到即,

9、TOC o 1-5 h z其中為所在的鄰域的二進編碼為移位操作其后面的數(shù)字表示移位jij位數(shù).定義為塊的位二進制編碼特征表示簡稱塊的編碼iii對一幅圖像提取它所有塊的編碼需要計算各個塊的灰度值和、計ii算各個鄰域的編碼、計算各個塊的編碼等共步圖像最外一圈的jii塊的編碼無定義對于邊長為的圖像上述運算的時間復雜度為i顯然表示塊的灰度與相鄰個塊灰度的分布序關系體現(xiàn)了圖ii像灰度的相對值因此對整體灰度值的變化具有相對的穩(wěn)定性通過對塊尺寸的選擇,可以改變圖像處理粒度的大小,以改變抵抗不同頻率噪聲的能力.2特征的匹配過程定義在待搜索圖上以模板的長、寬為橫向、縱向步長從的左上角開始按模板的大小劃分得到的子

10、圖稱為限制塊記作其中為限制塊左上角頂點在搜索圖上的坐標這樣劃分后如果在搜索圖的右側或底部有剩余部分則相應地從的最右側開始向左或從最底部開始向上劃分出一列或一行限制塊使得全部限制塊可以完全覆蓋搜索圖這樣得到的圖上的限制塊的數(shù)量為H定義限制塊與模板都是尺寸為義的圖像各自的塊特征集合用階方陣與表示稱為特征編碼矩陣這里為塊的邊長在與作特征比較時即比較每一個元素與中每一個元素是否相等如果相等則記下矩陣中的行號、列號三應用程序設計：1數(shù)字圖像匹配相關函數(shù)：函數(shù)功能：從圖像文件等格式中讀入圖像數(shù)據(jù)。格式：文件名，（圖像文件格式說明:文件名為指定圖像文件名稱的字符串。圖像文件格式為圖像文件格式的字符串。文件名

11、必須在當前目錄或路徑中，如果找不到則錄找文件名圖像文件格式為無符號位整數(shù)。如果文件灰度圖像詳細解釋見“算法設計”中的“概念解釋”）則為一個二維數(shù)組；如果文件是一個真彩色圖像，則是一個三維數(shù)組lXX）。舉例：讀取存放在盤下的名為的圖像，并將其整個點陣數(shù)據(jù)傳給注：在中申請一個變量時不需要指定類型，其具體類型由系統(tǒng)自動判別函數(shù)功能：交圖像寫入圖像文件（以等格式）。格式：文件名，圖像文件格式說明：文件名為指定圖像文件名稱的字符串。圖像文件格式是指定圖像文件的保存格式的字符串。如果是一個無符號位整數(shù)表示的灰度圖像或真彩色圖像，直接交數(shù)組中的值寫入文件。如果為雙精度浮點數(shù)，首先使用自動將數(shù)組中的值變換為無

12、符號位整數(shù)，即交0，1范圍內的浮點數(shù)變換為0，255范圍內8位整數(shù)，然后寫入文件。函數(shù)功能：圖像顯示格式：說明：在圖形窗口顯示，其中為用函數(shù)賦值的變量是一個圖形矩陣。函數(shù)功能：轉換圖像或顏色映像表為灰度圖像。格式：I=rgb2gray(RGB)說明：為用函數(shù)賦值的變量是一個圖形矩陣；為一個表示灰度圖像的二維數(shù)組。函數(shù)功能：圖像剪裁格式：III說明：格式I為交互方式，顯示輸入圖像，等待用戶用鼠標定義要剪裁的矩形。格式I為非交互方式，對灰度圖像進行非交互方式的剪裁操作，通過四元素向量指定剪裁矩陣，這些值是定義在坐標系中的，和是圖像矩陣的一處坐標，是以此坐標為準向右的延伸的長度，是以此坐標為準向下延

13、伸的長度，由此構成要剪裁的矩形區(qū)域。函數(shù)功能：用來求出某一個變量的的大小參數(shù)格式說明：用函數(shù)賦值的變量是一個圖形矩陣且為灰度圖像，用來保存的高度值，用來保存的寬度值。2數(shù)字圖像匹配函數(shù)：(1)基于灰度的摸板匹配算法函數(shù)功能：求圖像均值格式：說明：函數(shù)可用來計算圖像矩陣的均值是雙精度標量。舉例：m=mean2(I)運行結果：函數(shù)功能：求圖像間相關系數(shù)格式：說明：函數(shù)計算圖像矩陣與的相關系數(shù)，矩陣與的大小相同，如均為：X。計算公式如下：(ABA)(BB)TOC o 1-5 h zmnmn HYPERLINK l bookmark16 r,mn(AA)2).(BB)2) HYPERLINK l bo

14、okmark18 mnmnmnmn其中，AB。相關系數(shù)是表示矩陣與的線性聯(lián)系密切程度的一個統(tǒng)諦睛系數(shù)值是小于等于的正數(shù)。值為時，表示矩陣與的線性聯(lián)系最為密值為時，表示矩陣與不相關，線性聯(lián)系最最弱。舉例：，運行結果如下：函數(shù)功能：在線性坐標系中繪制二維數(shù)據(jù)格式：（顏色說明：可以為任意一個失量。舉例：在縱坐標為處畫一條平行于x由，長度為的紅色的線段其中表示從到以默認值為的速度遞增表示紅色命令功能：可以交圖形添加到現(xiàn)有的圖形中，從而實現(xiàn)圖形的疊加。如果交設置為，則在繪圖時交不再清除已經(jīng)存在的圖形，而是交新的數(shù)據(jù)直接加到當前圖形中。當新的數(shù)據(jù)落在原來的坐標范圍之外時，交自動進行坐標比例的調整。功能：實

15、現(xiàn)兩幅數(shù)字圖像的匹配實現(xiàn)細節(jié)：I讀取原圖像搜索圖和匹配圖像（模板圖）讀取原圖像讀取匹配圖像II將原圖像和匹配圖像灰度化，以便于計算將原圖灰度化將匹配圖灰度化III獲取原圖像和匹配圖像的尺寸獲取匹配圖象的大小尺寸，其值為高，寬將匹配圖象的長度值賦給將匹配圖象的寬度值賦給獲取原圖象的大小尺寸，其值為高，寬將原圖象的長度值賦給將原圖象的寬度值賦給w顯示出原圖，并調用w顯示出原圖，并調用命令，以實現(xiàn)以后在原圖像上做標記顯示原圖像V在原圖像上進行掃描，橫坐標到原圖像的寬度減去匹配圖的寬度，縱坐標到原圖像高度減去匹配圖的高度；對于原圖像上掃描到的每個像素點，以此像素點的坐標為準，剪裁參見函數(shù)說明出一個大小

16、和匹配圖大小一樣的矩陣，此像素點為該矩陣的左上角第一個像素；對這個臨時矩陣和匹配圖像的矩陣取相關系數(shù)值參見函數(shù)說明規(guī)定相關系數(shù)的值(我設為)如果取得的相關系數(shù)值大于規(guī)定值，可認為在以此坐標為左上角第一個像素的臨時矩陣和匹配圖矩陣是相匹配的；如果在此坐標下是匹配的話，則以此為左上角，將匹配區(qū)域標記出來。取得相關系數(shù)規(guī)定值為下面用函數(shù)在原圖的坐標系上畫出匹配區(qū)域(2)基于灰度的快速模板匹配算法函數(shù)功能：十進制數(shù)轉換成二進制數(shù)格式：位數(shù)說明：為十進制數(shù)，為轉換后的二進制序列，位數(shù)為的二進制位數(shù)。舉例：將轉換為位的二進制運行結果：01101函數(shù)功能：二進制數(shù)轉換成十進制數(shù)格式：說明：為二進制序列，為轉

17、換后的十進制數(shù)。舉例：將二進制為的序列轉化為十進制運行結果：45函數(shù)功能：根據(jù)輸入的四項參數(shù)的大小順序，返回一個5位二進制的編碼格式：shunxu=getOrder(temp_shunxu_1,temp_4)舉例：運行結果：10000函數(shù)功能：根據(jù)特征塊的左上角坐標得出該特征塊的特征值，特征值的定義詳見基于灰度的快速匹配算法中的定義2。格式：說明：為特征塊左上角的橫坐標與衷縱坐標，為圖像矩陣，為最佳分塊標準值函數(shù)功能：求向量中的最小值格式：說明：為一個向量，中保存的是中的最小值舉例：運算結果2函數(shù)功能：求向量中的最大值格式：說明：為一個向量，中保存的是中的最大值舉例：運算結果66函數(shù)功能：求向

18、量中和格式：說明：為一個向量，中保存的是中的和舉例：運算結果函數(shù)實現(xiàn)細節(jié)：I讀取原圖像搜索圖和匹配圖像（模板圖）讀取原圖像讀取匹配圖像II將原圖像和匹配圖像灰度化，以便于計算將原圖灰度化將匹配圖灰度化III獲取原圖像和匹配圖像的尺寸及分塊標準值獲取匹配圖象的大小尺寸，其值為高，寬將匹配圖象的長度值賦給將匹配圖象的寬度值賦給獲取原圖象的大小尺寸,其值為高，寬將原圖象的長度值賦給將原圖象的寬度值賦給尺兩個邊中最短的一條邊以分成義特征塊為基準，得出最佳分塊標準值W計算及記錄匹配圖的個特征塊的信息；顯示出原圖，并調用命令，以實現(xiàn)以后在原圖像上做標記去掉外圍的一圈取順時針為正方向temp_pipeitu

19、(t,s)=bin2dec將二進制編碼值(特征塊值)以十進制的形式保存在義的矩陣中的對應位置V在原圖像上進行掃描，橫坐標到原圖像的寬度減去匹配圖的寬度，縱坐標到原圖像高度減去匹配圖的高度；對于原圖像上掃描到的每個像素點，以此像素點的坐標為準，剪裁參見函數(shù)說明出一個大小和匹配圖大小一樣的矩陣，此像素點為該矩陣的左上角第一個像素；對這個臨時矩陣進行特征值提取，并和匹配圖的特征值進行比較，如果特征塊的匹配數(shù)大于一定數(shù)(我設定為大于等于8),可認為在以此坐標為左上角第一個像素的臨時矩陣和匹配圖矩陣是相匹配的；如果在此坐標下是匹配的話，則以此為左上角，將匹配區(qū)域標記出來。:去掉3外圍的一圈=1%:取3順

20、時針為正方向如果特征值相同，就加下面用函數(shù)在原圖的坐標系上畫出匹配區(qū)域plot(j:j+pipei_width,i,r);plot(j:j+pipei_width,i+pipei_height,r);plot(j,i:i+pipei_heightr,);plot(j+pipei_width,i:i+pipei_height,r);endendendD3界面設計：本程序采用交互式來演示圖像匹配過程，采用自帶的界面編輯器。主要涉及圖像顯示，各個功能按鈕對顯示的圖像所進行的操作(如：放大，縮小，旋轉和圖像匹配)，文件對話框(用來讀取和保存圖像)。uigetfiie函數(shù)口功能：打開磁盤內的文件格式：p

21、name,adrname口uigetfiie(格式,名稱)口說明：pname為選擇的文件名，adrname為路徑；uigetfiie通常和if口exist(文件名連用，如果存在文件的話，則把相應的值傳給pname和口口口口口口口口口舉例：pname,adrname口uigetfiie(*.jpg,*.bmp)口運行結果:Duiputfiie功能：保存文件格式：temp_pname,temp_adrname=uiputfiie(strcat(adrname,pname),保存圖象);axes函數(shù)：口功能：固定地點，此程序用此函數(shù)在固定區(qū)域內顯示圖像舉例：axes（handles.axesl）;%handles.a