腦外科CT圖像的綜合邊緣提取算法

寧 祎，高紅偉

（河南工業大學 機器人研究所，河南 鄭州 450001）

眾所周知，腦外科手術是最復雜而危險的外科手術，只要一個微小的錯誤動作，就足以造成中樞神經系統的永久損害。因此若需進行腦外科手術，事先必須縝密的規劃，并做好充分的準備工作，包括分毫不差地確認腫瘤位置、決定最佳的切除方式以及了解腦細胞與其他重要大腦部位的相對位置等。 隨 著計算機 技術 、CT （Computed Tomography）、MRI（Magnetic Resonance Imaging）等醫學影像技術的不斷發展，虛擬現實技術也越來越多地應用到現代醫療領域[1]。利用計算機圖像處理和數據可視化技術，根據醫學影像設備提供的二維斷層圖像，進行人體器官的三維重建已是現代醫學重要發展方向之一。圖像最基本的特征是邊緣，邊緣是區域屬性發生突變的地方，也是圖像中不確定性最大的地方和圖像信息最集中的地方，包含著圖像的豐富信息[2－3]。圖像的邊緣對大腦的三維重建具有重要意義，是三維重建的重要基礎，圖像邊緣提取的正確性和可靠性將直接影響到三維重建后的效果。CT二維圖像的邊緣提取作為器官三維重建的第一步，如何快速、準確地提取圖像邊緣信息一直受到國內外學者的關注。1959年，Julez[4]最早提出了邊緣檢測；1965年，Roberts[5]最早開始系統地研究邊緣檢測；1986年，Canny提出基于最優化算法的邊緣檢測算子。經過50多年的研究，已經出現了許多不同的邊緣檢測方法，各自有其特點和局限性，但針對腦外科CT圖像的研究報道并不多。因此，本研究針對腦外科CT圖像紋理復雜、細節豐富等特點，設計了一種采用圖像增強、Canny算子與數學形態學相結合的綜合邊緣提取算法，達到了理想的邊緣提取效果。

1 本文綜合算法的提出與圖像增強

1.1 本文綜合算法的提出

經典的邊緣檢測算法是對原始圖像中像素的某小鄰域來構造邊緣檢測算子，常用的邊緣檢測算子有Prewitt、Log、Roberts和Canny算子等。近些年，隨著數學形態學理論不斷完善與發展，數學形態學在圖像邊緣檢測中也得到了廣泛的應用[6]。Prewitt算子的特點是抑制噪聲效果好，但對邊緣的定位不準確；Log算子的特點是邊緣定位較準確，但對噪聲較敏感，對圖像中的某些邊緣容易產生雙重響應；Roberts算子的特點是計算時采用較少的像素，邊緣定位準確，但是對噪聲高度敏感。Canny算子是在圖像處理中人們公認的功能較為強大的檢測算子之一，但有時對于灰度分布不均勻的圖像，容易將非邊緣點檢測為邊緣。數學形態學是利用一個結構元素去探測一個圖像，看是否能夠將這個結構元素很好地填放在圖像的內部，同時驗證填放結構元素的方法是否有效。數學形態學的腐蝕算法會導致圖像區域的分裂，而膨脹算法容易產生零星的像素點。

從目前的研究情況來看，單純采用某種邊緣檢測算法很難實現腦外科CT圖像準確的邊緣提取。因而本文針對腦外科圖像紋理復雜、細節豐富等特點，考慮到醫學圖像數據多層切片間的相關性，提出了一種綜合的圖像邊緣檢測算法，該算法集成了用平均值濾波器對圖像增強，Canny邊緣檢測算子與數學形態學相結合的算法。

該算法的實現過程是在Matlab 7.5軟件上對圖像進行處理分析，運用了Matlab圖像處理工具箱 （Image Processing Toolbox）中的大量函數。本文邊緣提取算法的主要步驟為：

1）對原始圖像運用平均值濾波器來對圖像增強；

2）對增強后的圖像用Canny算子進行邊緣檢測；

3）利用數學形態學中的膨脹算法進行邊緣處理。

1.2 圖像的增強

圖像增強是數字圖像處理過程中經常采用的一種方法，這種根據圖像的特點或存在的問題采取的簡單改善方法或者加強特征的措施稱為圖像增強（Image Enhancement）。圖像增強的目標有2個：第一是改善圖像的視覺效果，提高圖像的清晰度；第二是使圖像變得更利于計算機處理。所有的圖像處理就是對這些組成圖像的象素灰度值進行增或減的計算處理。從增強處理的作用域出發，圖像增強的方法分為兩大類：空間域方法和頻域方法。由于CT圖像存在大量的噪聲，本文采用了平滑線性濾波器增強的方法，它是空間域的一種。平滑線性濾波器的輸出是包含在濾波掩模鄰域內像素的簡單平均值，因此也稱為均值濾波器[7]。

平均值濾波器定義如下，對于給定的圖像 f（i，j）中的每個像素點（m，n），取其鄰域S。設S含有M個像素，取其平均值作為處理后所得圖像像素點（m，n）處的灰度值，即用一像素鄰域內各像素灰度平均值來代替該像素原來的灰度值。鄰域S的形狀和大小根據圖像特點確定，一般取的形狀是正方形、矩形及十字形等，鄰域S的形狀和大小可以在全圖處理過程中保持不變，也可根據圖像的局部統計特性而變化，點（m，n）一般位于鄰域 S的中心。

文中根據腦外科圖像的特點和大量的實驗比較，得出比較理想的效果是采用鄰域S為正方形5×5鄰域，像素點（m，n）位于鄰域S的中心，共25個像素點，通過公式（1）計算該點的灰度平均值：

式中，fˉ（m，n）為像素灰度平均值；m，n分別為像素點的橫、縱坐標值；i，j分別為鄰域S的橫、縱坐標值。

原始圖像如圖1所示。在Matlab命令窗口中，首先應用imread（filename）函數將原始圖像導入到Matlab中，然后應用filter2（fspecial（‘average’，5），I）/255 語句對原始圖像進行增強，增強后圖像如圖2所示。

圖1 原始圖像Fig.1 Original image

圖2 圖像增強后Fig.2 Image enhancement

從圖2可以看出，經過平滑線性濾波器后的圖像，去除了CT圖像中因掃描存在的噪聲、平滑了圖像信號，為后面能夠更好地進行邊緣檢測打下了基礎。

2 Canny算子的邊緣檢測

邊緣檢測是對圖像的邊緣進行處理獲得需要的邊緣信息，常用的邊緣檢測算子有Prewitt、Log、Roberts和Canny算子等。每種邊緣檢測算子都有各自的優缺點，Prewitt算子抑制噪聲效果好，但對邊緣的定位不準確，Log和Roberts算子都是邊緣定位較準確，但對噪聲較敏感，因此本文選用了具有高信噪比、高定位精度和單邊緣響應等優良性能的Canny算子。

Canny算子采用二維高斯函數的任意方向上的一階方向導數為噪聲濾波器，通過與圖像卷積進行濾波，然后對濾波后的圖像尋找圖像梯度的局部最大值，以此來確定圖像邊緣[8]。由于噪聲的影響，一個閾值的判斷往往是不夠的，邊緣信號的響應只有近一半是大于這個閾值的，由此造成了邊緣斷裂。如果降低這個閾值，又會發現一些錯誤的“邊緣”。為了解決這個問題，Canny提出了雙閾值方法[9]。Canny算子進行邊緣檢測的具體步驟如下：

1）用二維高斯濾波器對圖像進行濾波，圖像通過二維高斯濾波函數。

進行平滑，抑制圖像噪聲，其中σ為高斯濾波器參數，它控制著平滑程度。

2）用導數算子找到圖像灰度沿著兩個方向的偏導數（Gx，Gy），并計算出梯度的大小和方向。

3）對梯度幅值進行非極大值抑制。僅僅得到全局的梯度并不足以確定邊緣，因此為確定邊緣，必須保留局部梯度最大的點，抑制非極大值。若某個像素的灰度值與其梯度方向上前后兩個像素的灰度值相比不是最大的，那么該像素值置為0，即不是邊緣，這個過程稱為“非極大值抑制”。

4）用雙閾值法檢測和連接邊緣。對梯度取兩次閾值Tl和T2，T1=T2，0＜α＜1，可根據噪聲的大小適當調整。 用 T1、T2兩個閾值對非極大值抑制圖像進行雙閾值化后，可得兩個檢測結果，分別記為圖像1和圖像2。圖像l閾值較低，則保留了較多信息；圖像2閾值較高，所以噪聲較少，但會造成邊緣信息的損失。于是以圖像2為基礎，以圖像l為補充，連接圖像的邊緣。

應用 Matlab 圖像處理工具箱中 edge（I，’ ’，thresh）函數，幾種經典的邊緣檢測算子對增強后的圖像進行邊緣檢測，效果如圖3所示。

通過對圖3的對比分析，可以清楚地看到：圖 3（a）Prewitt算子對噪聲有抑制作用，檢測的圖像沒有較多孤立的邊緣像素點，但對邊緣的定位不是很準確，圖像的邊界損失較多；圖3（b）Log算子檢測的圖像邊緣信息損失較少，但對圖像中的某些邊緣產生雙重響應；圖3（c）Roberts算子對邊緣定位比較準確，但存在較多孤立的像素點，抑制噪聲能力差；圖3（d）用Canny算子檢測后的效果最好，由于Canny算子具有非極大值抑制和雙閾值法連接邊緣等優點，比前3種提取的邊緣清晰而且連續，虛假邊緣少，但是提取出的邊緣仍然存在不光滑和斷裂的現象。為了實現更好的邊緣檢測效果，下一步就要用本文提出的數學形態學算法來具體解決該問題。

3 數學形態學邊緣檢測

3.1 數學形態學的基本原理

數學形態學是一門綜合了多學科知識的交叉科學，它建立在嚴格的數學理論基礎上。其基本思想是用具有一定形態的結構元素 SE（structuring element）來度量和提取圖像中的對應形狀，以達到對圖像的分析和識別的目的[10]。

膨脹和腐蝕是形態學的兩個基本操作，膨脹過程是在圖像中對圖象的邊界添加像素點，向外擴張，而腐蝕是它的逆過程。膨脹和腐蝕這兩種變換都對灰度值變化明顯的圖像邊緣較為敏感，膨脹可以填充圖像中的小孔及圖像邊緣上較小的凹陷部分，腐蝕可以消除圖像中細小的成分。下面介紹2個基本運算的定義：

膨脹運算定義為：

設A是原始圖像，B是結構元素，式中A、B為Z2中的集合，φ為空集，B?為B的映射，A被B膨脹，記為A⊕B，⊕為膨脹算子。該式表明膨脹過程為B首先做關于原點的映射，然后平移x。A被B膨脹是被所有x平移后與A至少有一個非零的公共元素。

腐蝕運算定義為：

圖3 邊緣檢測Fig.3 Edge detection

設A是原始圖像，B是結構元素，式中A、B為Z2中的集合，A被B腐蝕，記為AΘB。該式表明A被B腐蝕的結果是所有使B被x平移后包含于A的點x的集合。

3.2 數學形態學方法在圖像處理中的運用

應用 Matlab圖像處理工具箱中 imdilate（I，SE）函數，對經過Canny算子邊緣檢測后的圖像進行數學形態法中的膨脹運算。利用數學形態學進行圖像處理時，選擇簡單、表現力強的結構元素是關鍵，是形態變換中最重要的參數；其次，還要綜合考慮目標體的清晰度和噪聲的大小來選取結構元素的大小。一般目標體輪廓不清晰時，選擇較小的結構元素；噪聲顆粒較大時，選擇較大的結構元素。本文根據腦外科圖像的特點和大量的實驗比較，采用平面結構元素為5×5的正方形，處理后的效果如圖4所示。

圖4 膨脹后效果圖Fig.4 Dilated image

從圖4中可以看到，經過膨脹后得到的邊緣清晰而且光滑連續，基本上消除了圖3（d）中所存在的邊緣模糊和斷裂現象，取得了較好的效果。

4 結 論

從以上的實驗結果和分析可以看出：針對腦外科CT圖像的紋理復雜、細節豐富等特點，本文提出的用平均值濾波器對圖像增強、Canny邊緣檢測算子與數學形態學相結合的綜合邊緣提取算法，不僅具有簡單快捷的處理能力，而且具有很好的邊緣提取、較強的抗噪性能，提取出的邊緣清晰連續，能夠比較準確的提取腦外科CT圖像的邊緣，有效地保留了圖像中的重要信息，為腦外科CT圖像三維重建和虛擬手術的研究等打下了堅實的基礎。

[1]M_engen S B， Marjunath B S， KenneyC.Image segmentation using multi-region stability and edge strength[C]//IEEE International Conference on Image Processing.Barcelona：Spain，2003：429-432.

[2]黃金國，戴志明，周培源.一種改進的基于梯度的圖像邊緣檢測算法[J].武漢科技學院學報，2010，23（3）：33-35.

HUANG Jin-guo， DAI Zhi-ming， ZHOU Pei-yuan.An image edge detection algorithm based on gradient[J].Journal of Wuhan University of Science and Engineering，2010，23（3）：33-35.

[3]蘇恒陽，袁先珍.一種改進的Canny的圖像邊緣檢測算法[J].計算機仿真，2010，27（10）：242-245.

SU Heng-yang， YUAN Xian-zhen.An improved image edge detectionalgorithm ofCannyoperator[J].Computersimulation，2010，27（10）：242-245.

[4]Julez B.A method of coding TV signals based on edge detection[J].BellSystem Tech，1959，38（4）:1001-1020.

[5]Roberts L D.Machine perception of three-dimension solids in optica and electro-optimal information processing[C]//Massachusetts Institute of Technology Press,1966:157-197.

[6]Huang C P，Wang R Z.An intergrated edge detection method using mathematical morphology[J].Pattern Recgnition and Image Analysis，2006，16（3）：406-412.

[7]孫兆林.MATLAB 6.x圖像處理[M].北京：清華大學出版社，2002.

[8]張斌,賀賽先.基于Canny算子的邊緣提取改善方法[J].紅外技術，2006，28（3）：165-169.

ZHANG Bin,HE Sai-xian.An improved method of edge detection based on Canny operator[J].Journal of Infrared technology，2006，28（3）：165-169.

[9]王家文，曹宇.MATLAB 6.5圖形圖像處理[M].北京：國防工業出版社，2004.

[10]張小琴,楊闊,郝葆青.基于Canny算子與數學形態學的細胞邊緣提取[J].生物醫學研究，2009，28（4）：275-279.

ZHANG Xiao-qin,YANG Kuo,HAO Bao-qing.Cell edge detection based on Canny operatorand Mathematical morphology[J].Journal of Biomedical Research，2009，28（4）：275-279.