摘  要： 精確的圖像配準(zhǔn)是超分辨率重建取得成功的前提。為此，針對大多數(shù)圖像配準(zhǔn)算法存在精度低、抗噪聲能力差等缺點，提出一種基于紋理特征的圖像自動配準(zhǔn)方法。首先用Canny邊緣檢測算子提取裂縫紋理圖像，然后用Photoshop圖層工具，通過對待配準(zhǔn)幀關(guān)于參考幀作變換，得到空間幾何變換參數(shù)，最后用Matlab編程實現(xiàn)待配準(zhǔn)幀的縮放、旋轉(zhuǎn)和平移。對配準(zhǔn)后的酸蝕巖板裂縫圖像，進行超分辨率重建，結(jié)果證明，該方法具有抗噪聲能力強、配準(zhǔn)精度高、易于編程實現(xiàn)等優(yōu)點，取得了較好的效果。
關(guān)鍵詞： 圖像配準(zhǔn)；酸蝕巖板裂縫圖像；圖像紋理；插值重建

　裂縫型油氣藏是21世紀(jì)石油增儲上產(chǎn)的重要領(lǐng)域之一，但裂縫儲層本身的復(fù)雜性在一定程度上增加了研究難度[1]。以往對儲層巖石的研究，大多是憑借個人經(jīng)驗直觀地觀察鉆取巖心樣品來獲得裂縫信息。但由于巖心樣品取出后，地應(yīng)力釋放，得到的裂縫信息有誤差。為獲取更真實的信息，借助自主開發(fā)的儲層巖石微觀裂縫可視化裝置，模擬地下高壓環(huán)境，用工業(yè)電子內(nèi)窺鏡采集巖石裂縫圖像。但由于受觀察孔內(nèi)光線不足、耐壓玻璃及釜體內(nèi)液體反光和鏡頭抖動等因素的影響，往往得到的是多幀低分辨率模糊形變圖像。為了獲得高分辨率的巖石裂縫圖像，需要對多幀不同角度下順序拍得的圖像進行超分辨率重建，其中極為關(guān)鍵的一步就是圖像配準(zhǔn)。圖像配準(zhǔn)不僅是重建的前提，而且決定了重建的準(zhǔn)確度[2]。
　目前研究較多的配準(zhǔn)算法是基于仿射變換或圖像信息熵的方法[3]。但是由于參考幀與待配準(zhǔn)幀的拍攝時間與條件的變化，以及成像設(shè)備的限制，使得參考幀與待配準(zhǔn)幀之間存在著嚴(yán)重的灰度差異，因而這些配準(zhǔn)算法有著計算量大、配準(zhǔn)精度低等不足，不適合實際工程應(yīng)用中抗噪聲和高精確度的要求。為此，本文用邊緣檢測算子Canny提取紋理特征圖像，然后利用Photoshop圖層工具對待配準(zhǔn)幀關(guān)于參考幀作變換，得到空間幾何變換參數(shù)，最后通過Matlab編程實現(xiàn)待配準(zhǔn)圖像的縮放、旋轉(zhuǎn)和平移。對配準(zhǔn)后的酸蝕巖板裂縫圖像進行超分辨率重建，獲得了較好的效果，表明本配準(zhǔn)方法是行之有效的。

2 邊緣檢測獲取紋理圖像
　邊緣是圖像的最基本特征，它具有勾畫出物體輪廓、能被局部定義以及傳遞大部分圖像信息等優(yōu)點[6]，通過邊緣檢測獲取含有紋理特征的圖像，進而實現(xiàn)圖像匹配，是一種簡單而又高效的方法。常用的邊緣檢測算子有微分算子、拉普拉斯算子(Laplacian)和Canny算子。微分算子(如Robert算子、Sobel算子、Prewitt算子)對邊緣定位不是很準(zhǔn)確，圖像的邊緣不止一個像素，對灰度突變十分敏感；Laplacian算子對噪聲比較敏感，所以實際應(yīng)用中很少用該算子檢測邊緣；Canny邊緣檢測方法抗噪聲能力強，它的突出優(yōu)點是，使用兩種不同的閾值分別檢測強邊緣和弱邊緣，并且當(dāng)弱邊緣與強邊緣相連時，才會將弱邊緣包含在輸出圖像中，因此能夠檢測到真正的邊緣。鑒于此，本文選用了具有很好邊緣檢測性能的Canny算子。
2.1 Canny邊緣檢測方法
　Canny算子的實現(xiàn)過程如下：
　(1)首先用高斯濾波器平滑圖像以消除噪聲。
　(2)利用一階偏導(dǎo)的有限差分計算圖像灰度沿著兩個方向的偏導(dǎo)數(shù)(Gx，Gy)，并求出梯度的大?。?br />
　(4)把邊緣的梯度方向大致分為水平、豎直、45°和135°四種。通過梯度的方向，就可以找到這個像素梯度方向的鄰接像素。
　(5)遍歷圖像，若某個像素的灰度值與其梯度方向上前后兩個像素的灰度值相比不是最大時，則將該像素置于0，即不是邊緣。
　(6)使用累計直方圖計算兩個閾值，大于高閾值的一定是邊緣，小于低閾值的一定不是邊緣，介于兩者之間的，則判斷這個像素的臨界像素中有沒有超過高閾值的邊緣像素，如果有，則是邊緣；否則不是邊緣[7]。
2.2 邊緣檢測的實現(xiàn)
　在Matlab圖像處理工具箱中，提供了edge函數(shù)利用Canny邊緣檢測算法來檢測灰度圖像的邊緣。
調(diào)用Canny算子檢測巖板裂縫邊緣圖像的程序如下：
　r=imread(′E：\image\uvs100413-011.bmp′)；
　g=imread(′E：\image\uvs100413-034.bmp′)； %讀入圖像
　r=rgb2gray(r)；
　g=rgb2gray(g)； %圖像灰度化
　R=edge(r，′canny′)；
　G=edge(g，′canny′)； %Canny算子提取圖像邊緣
　R_edge=imwrite(′E：\image\R_edge.bmp′)；
　G_edge=imwrite(′E：\image\G_edge.bmp′)； %保存圖像
　應(yīng)用Canny邊緣檢測算子檢測灰度圖像的邊緣效果如圖1~4所示。對比圖1和圖3以及圖2和圖4，可以看到裂縫的邊緣和方向都比較準(zhǔn)確，基本保留了原始圖像的紋理特征。

3 求取空間幾何變換參數(shù)
　為簡化計算，通過Photoshop對獲取的參考幀和待配準(zhǔn)幀相應(yīng)的邊緣輪廓圖進行反相、去除冗余信息、增色操作，得到參考幀和待配準(zhǔn)幀僅包含裂縫信息的圖像。
　圖層是Photoshop最基本也是最重要的工具，在Photoshop下打開如圖5和圖6所示的圖像，選中并復(fù)制圖6，以圖5為背景新建圖層，默認(rèn)的新建圖層名稱為圖層1，將圖6粘貼于選中的圖層1上。把圖層1的不透明度選項設(shè)置為50%，則可以看到如圖7所示的一幅兩條裂縫交叉的圖像。這時以參考幀裂縫圖像為背景(其屬性為不可修改)，在編輯欄中先選擇自由變換，記錄下初始待配準(zhǔn)圖像的位置信息，順時針旋轉(zhuǎn)待配準(zhǔn)圖像直至兩裂縫平行，然后根據(jù)裂縫紋理特征，平移待配準(zhǔn)圖像，直至兩裂縫基本重合。兩圖像之間若有縮放，可選擇編輯欄中變換-縮放、縮小或放大待配準(zhǔn)圖像，直至最后得到如圖8所示的兩圖像的紋理特征達到最佳重合，記錄下變換后待配準(zhǔn)圖像的位置信息。最終得到的幾何變換參數(shù)如下：縮放因子k=1，水平方向平移量Δx=0，垂直方向平移量Δy=0.4，順時針旋轉(zhuǎn)角度θ=2.9°。

4 圖像配準(zhǔn)和重建結(jié)果分析
4.1 待配準(zhǔn)幀空間幾何變換的實現(xiàn)
　空間幾何變換s的實現(xiàn)流程圖如圖9所示。

4.2 重建實驗結(jié)果分析
　圖像配準(zhǔn)的目的是獲得超分辨重建后的高分辨率圖像，為此，本文實現(xiàn)了基于插值的超分辨率重建算法，得到用紋理特征配準(zhǔn)的重建圖像如圖12所示，用全參數(shù)配準(zhǔn)的重建圖像如圖13所示。對比圖12和圖13，可以明顯看出，圖13比圖12效果好。其原因是：以往大部分的圖像配準(zhǔn)算法是基于圖像灰度信息進行的配準(zhǔn)，也就是說，僅適用于當(dāng)兩幅圖像像素的空間分布改變，而各處亮度值不變的情況[8]，而且對噪聲極其敏感，魯棒性差，對待配準(zhǔn)幀圖像要求很苛刻。而本文提出的基于紋理特征的圖像自動配準(zhǔn)方法，有很強的抗噪聲能力和較高的配準(zhǔn)精度。

　隨著圖像處理在工程應(yīng)用中的廣泛使用，圖像的超分辨率重建技術(shù)，因其能在不增加硬件成本的前提下，通過軟件提高圖像的分辨率而得到人們越來越多的關(guān)注，而精確的圖像配準(zhǔn)是重建成功的先決條件。在實際應(yīng)用中，特別是工程現(xiàn)場，順序獲取的同一場景圖像之間的亮度值不可能完全相同，因此，本文提出的基于紋理特征圖像自動配準(zhǔn)方法，也更加符合于工程實際的需要。所得到的高精度配準(zhǔn)圖像，也為后續(xù)的超分辨率重建研究奠定了堅實的基礎(chǔ)。
參考文獻
[1] 周新桂，操成杰，袁嘉音.儲層構(gòu)造裂縫定量預(yù)測與油氣滲流規(guī)律研究現(xiàn)狀和進展[J].地球科學(xué)進展，2003，18(3)：398-402.
[2] BANNORE V. Iterative-interpolation super-resolution image reconstruction： A Computationally Efficient Techique[D]. Spring， 2009： 77-91.
[3] ZITOVA B， FLUSSER J. Image registration mathods： a survey[J]. Image and Vision Computing， 2003，21：977-1000.
[4] 惠文華.多傳感器遙感圖像配準(zhǔn)方法研究[J].西北大學(xué)學(xué)報，2006(8)：648-650.
[5] 羅俊，黃本雄.基于全參數(shù)的分層搜索圖像配準(zhǔn)[J].現(xiàn)代計算機，2005(11)：40-43.
[6] 任文杰.圖像邊緣檢測方法的研究[D].濟南：山東大學(xué)，2008:9-12.
[7] 江笑嬋，萬振凱，陳利.基于Matlab邊緣提取的幾種方法的比較[J].電腦知識與技術(shù)，2005(10)：138-141.
[8] BERGEN J R， ANANDAN P. Hierachical model-based motion estimation[C]. Proceeding of the Second European Conference on Computer Vision， ECCV′92，1992.