【文章內容簡介】 或相似準則。生長準則一般可分為三種：基于區(qū)域灰度差準則、基于區(qū)域內灰度分布統(tǒng)計性質準則和基于區(qū)域形狀準則。分裂合并法是先將圖像分割成很多的一致性較強的小區(qū)域，再按一定的規(guī)則將小區(qū)域融合成大區(qū)域，達到分割圖像的目的。區(qū)域提取法的缺點是往往會造成過度分割，即將圖像分割成過多的區(qū)域，因此近年來針對這種方法的研究較少。閾值分割法是一種簡單高效的圖像分割技術。它通過設定不同的特征閾值，把圖像像素點分為若干類。在一幅圖像中用灰度等級表示各像素點的特征，許多閾值分割方法根據一維灰度直方圖或者二維灰度直方圖從背景中提取感興趣的目標。閾值的選取是閾值分割技術的關鍵，如果閾值選取過高，過多的目標點被誤歸為背景；閾值選3取過低，則會出現(xiàn)相反的情況。現(xiàn)有的大部分算法都是集中在閾值確定的研究上。閾值分割方法根據圖像本身的特點，可分為單閾值分割方法和多閾值分割方法：也可分為基于像素值的閾值分割方法、基于區(qū)域性質的閾值分割方法和基于坐標位置的閾值分割方法。2 基于 MATLAB 的圖像分割 MATLAB 的優(yōu)點 MATLAB 使用方便MATLAB 允許用戶以數學形式的語言編寫程序，用戶在命令窗口中輸入命令即可直接得出結果，這比 C++、Fortran 和 Basic 等等該機語言都要方便的多。而且它是用C 語言開發(fā)的，其流程控制語句與 C 語言中的相應語句幾乎一致。這給使用上帶來了方便，使我能較快的適應與使用 MATLAB 這門語言。MATLAB 內部函數豐富MATLAB 的內部函數提供了相當豐富的函數，這些函數解決許多基本問題，如矩陣的輸入。在其它語言中(比如 C 語言中)，要輸入一個矩陣，先要編寫一個矩陣的子函數，而 MATLAB 語言則提供了一個人機交互的數學系統(tǒng)環(huán)境，該系統(tǒng)的基本數據結構是矩陣，在生成矩陣對象時，不要求做明確的維數說明。與利用 C 語言或 Fortran1等等高級語言編寫數值計算的程序相比，利用 MATLAB 可以節(jié)省大量的編程時間。這就給用戶節(jié)省了很多的時間，使用戶可以把自己的精力放到創(chuàng)造方面，而把繁瑣的問題交給內部函數來解決。除了這些數量巨大的基本內部函數外，MATLAB 還有為數不少的工具箱。這些工具箱用于解決某些領域的復雜問題。MATLAB 強大的圖形和符號功能MATLAB 具有強大的圖形處理功能，它本身帶有許多繪圖的庫函數，可以很輕松地畫出各種復雜的二維和多維圖形。這些圖形可以在與運行該程序的計算機連接的任何打印機設備上打印出來，這使得 MATLAB 成為技術數據可視化的杰出代表。1譯為“公式編譯器” ， 它是世界上最早出現(xiàn)的計算機高級程序設計語言，廣泛應用于科學和工程計算領域。43 圖像分割的主要研究方法 圖像分割定義圖像分割是圖像處理中的一項關鍵技術，自 20 世紀 70 年代起一直受到人們的高度重視，至今已提出了上千種各種類型的分割算法，現(xiàn)提出的分割算法大都是針對具體問題的，并沒有一種適合于所有圖像的通用分割算法，而且近年來每年都有上百篇相關研究報道發(fā)表。然而，還沒有制定出選擇合適分割算法的標準，這給圖像分割技術的應用帶來許多實際問題。因此，對圖像分割的研究還在不斷深入之中，是目前圖像處理中研究的熱點之一 [1]。多年來人們對圖像分割提出了不同的解釋和表述，借助集合概念對圖像分割可給出如下定義：令 集 合 R 代 表 整 個 圖 像 區(qū) 域 ， 對 R 的 圖 像 分 割 可 以 看 做 是 將 R 分 成N 個 滿 足 以 下 條 件 的 非 空 子 集 R1， R2， R3， …， RN；（1）在分割結果中，每個區(qū)域的像素有著相同的特性；（ 2） 在 分 割 結 果 中 ， 不 同 子 區(qū) 域 具 有 不 同 的 特 性 ， 并 且 它 們 沒 有 公 共 特 性 ；（ 3） 分 割 的 所 有 子 區(qū) 域 的 并 集 就 是 原 來 的 圖 像 ；（ 4） 各 個 子 集 是 連 通 的 區(qū) 域 ；圖像分割是把圖像分割成若干個特定的、具有獨特性質的區(qū)域并提取出感興趣目標的技術和過程，這些特性可以是像素的灰度 、顏色、紋理等提取的目標可以是對應的單個區(qū)域，也可以是對應的多個區(qū)域。 圖像分割方法綜述圖像分割是指將圖像劃分為與其中含有的真實世界的物體或區(qū)域有強相關性的組成部分的過程。圖像分割是圖像處理和分析中的重要問題,也是計算機視覺研究中的一個經典難題。盡管它一直受到科研人員的重視,但是它的發(fā)展很慢,被認為是計算機視覺的一個瓶頸。迄今為止,還沒有一種圖像分割方法適用于所有的圖像,也沒有一類圖像所有的方法都適用于它。近幾年來,研究人員不斷改進原有方法并將其它學科的新理論和新方法引入圖像分割,提出了不少新的分割方法。本文對傳統(tǒng)的圖像分割方法進行分析。5 典型的圖像分割方法有閾值法，邊緣檢測法，區(qū)域法。分析各種圖像分割方法可以發(fā)現(xiàn)，它們分割圖像的基本依據和條件有以下 4 方面: (l)分割的圖像區(qū)域應具有同質性，如灰度級別相近、紋理相似等。 (2)區(qū)域內部平整，不存在很小的小空洞。 (3)相鄰區(qū)域之間對選定的某種同質判據而言，應存在顯著差異性。 (4)每個分割區(qū)域邊界應具有齊整性和空間位置的準確性。 現(xiàn)有的大多數圖像分割方法只是部分滿足上述判據。如果加強分割區(qū)域的同性質約束，分割區(qū)域很容易產生大量小空洞和不規(guī)整邊緣: 若強調不同區(qū)域間性質差異的顯著性，則極易造成非同質區(qū)域的合并和有意義的邊界丟失。不同的圖像分割方法總有在各種約束條件之間找到適當的平衡點 [2]。 邊緣檢測法 邊緣檢測原理邊緣(或邊沿 )是指其周圍像素灰度有階躍變化或“ 屋頂” 變化的那些像素的集合，也即邊緣是灰度值不連續(xù)的結果，這種不連續(xù)?？梢杂们髮捣奖愕臋z測到，一般常用一階導數和二階導數來檢測邊緣。邊緣廣泛存在于物體與背景之間、物體與物體之間、基元與基元之間。因此，它是圖像分割所依賴的重要特征，而邊緣信息是一種圖像的緊描述，所包含的往往是圖像中最重要的信息，故對圖像提取邊緣能極大地降低我們要處理的數據量 [3]。在討論邊緣算子之前，首先給出下列術語的定義。邊緣點：圖像中亮度顯著變化的點。邊緣段：邊緣點坐標[i,j]及其方向 θ 的綜合，邊緣的方向是可以是梯度角。邊緣檢測器：從圖像中抽取邊緣（邊緣點或邊緣段）集合的算法。輪廓：邊緣列表，或是一條邊緣列表的曲線模型。邊緣連接：從無序邊緣表形成有序邊緣表的過程。習慣上邊緣的表示采用順時針方向來排序。邊緣跟蹤：一個用來確定輪廓圖像（指濾波后的圖像）的搜索過程。邊緣點的坐標可以是邊緣位置像素點的行、列整數標號，也可以在子像素分辨率水平上表示。邊緣坐標可以在原始圖像坐標系上表示，但大多數情況下是在邊緣檢測濾波器的是輸出圖像的坐標系上表示，因為濾波過程可能導致圖像坐標平移或者縮放。邊緣段可以用像素點尺寸大小的小線段定義，或用具有方向屬性的一個點定義。在實6際應用中，邊緣點和邊緣段都稱為邊緣。雖然圖像邊緣點產生的原因不同，但他們都是圖像上灰度不連續(xù)點，或是灰度變化劇烈的地方。 經典的邊緣提取方法是考察圖像的每個像素在某個鄰域內灰度的變化，利用邊緣鄰近一階或二階方向導數變化規(guī)律，用簡單的方法檢測邊緣。 這種方法稱為邊緣檢測局部算子法。邊緣檢測算子 檢查每個像素的鄰域并對灰度變化率進行量化，也包括方向的確定。常用的邊緣檢測方法的基礎是微分運算，邊緣點對應于一階微分幅度大的點，或對應于二階微分的過零點。 傳統(tǒng)的邊緣檢測算法通過梯度算子來實現(xiàn)，在求邊緣的梯度時，需要對每個象素位置計算。經典的梯度算子模板有 Prewitt 模板、Canny 模板、Sobel 模板、Log 模板等[4]。 Canny 算子 Canny 邊緣檢測利用高斯函數的一階微分，在噪聲抑制和邊緣檢測之間尋求較好的平衡，其表達式近似于高斯函數的一階導數。Canny 邊緣檢測算子對受加性噪聲影響的邊緣檢測是最優(yōu)的。在高斯噪聲中，一個典型的邊緣代表一個階躍的強度變化。根據這個模型，好的邊緣檢測算子應該有 3 個指標：（1）低失誤概率，即真正的邊緣點盡可能少的丟失又要盡可能避免將非邊緣點檢測為邊緣：（2）高位置精度，檢測的邊緣應盡可能接近真實的邊緣；（3）對每一個邊緣點有惟一的響應，得到單像素寬度的邊緣?？材崴阕犹岢隽诉吘壦阕拥娜缦?3 個準則：信噪比準則：信噪比越大，提取的邊緣質量越高。信噪比 定義為：SNR （式 )??????wdxhGSNR2?其中 G（x）代表邊緣函數，h(x)代表寬度為 W 的濾波器的脈沖響應。定位精確度準則：7邊緣定位精度 L 如下定義： （式 )??dxhGLw????239。39。39。?其中 和 分別是 和 的導數。L 越大表明定位精度越高。??XG39。H39。 ??X單邊緣響應準則：為了保證單邊緣只有一個響應，檢測算子的脈沖響應導數的零交叉點平均距離應滿足：??39。fD （式 ）????2139。239。39。 ?????????????dxhfD?是 的二階導數??xh39。以上述指標和準則為基礎，利用泛函數求導的方法可導出坎尼邊緣檢測器是信噪比與定位之乘積的最優(yōu)逼近算子，表達式近似于高斯函數的一階導數。將坎尼 3 個準則相結合可以獲得最優(yōu)的檢測算子 [5]。 Prewitt 算子Prewitt 與 Sobel 算子的方程完全一樣，只是常系數 c=1。所以 和 可分別用卷積xsy模板表示為：1 0 1 1 0 1 1 0 1 1 1 10 0 01 1 18下面使用 MATLAB 圖像處理工具箱中的 edge 函數利用以上算子來檢測邊緣。Edge 函數提供許多微分算子模板，對于某些模板可以指定其是對水平邊緣還是對垂直邊緣（或者二者都有）敏感（即主要檢測是水平邊緣還是垂直邊緣）。Edge 函數在檢測邊緣時可以指定一個灰度閾值，只有滿足這個閾值條件的點才視為邊界點。Edge 函數的基本調用格式如下： （式 ）???,39。39。parmettyIedgBW?其中，I 表示輸入圖像，type 表示使用的算子類型，parameter 則是與具體算子有關的參數。 Sobel 算子 采用 33 鄰域可以避免在像素之間內插點上計算梯度。考慮下圖中所示的點周??ji,圍點的排列。Sobel 算子也是一種梯度幅值: （式 ）2yxsM??其中的偏導數用下式計算: （式 ）????670432 acacsx ?? （式 ）y??其中常系數 c=2。和其他的梯度算子一樣 和 可分別用卷積模板表示為:xsy0a12a7??ji,36a54a1 0 12 0 29 這一算子把重點放在接近于模板中心的像素點。Sobel 算子是邊緣檢測中最常用的算子之一。 Log 算子 Log 算子也就是 LaplacianGauss 算子，它把 Gauss 平滑濾波器和 Laplacian 銳化濾波器結合了起來，先平滑掉噪聲，再進行邊緣檢測。 這種方法的特點是圖像首先與高斯濾波器 1進行卷積，這一步既平滑了圖像又降低了噪聲，孤立的噪聲點和較小的結構組織將被濾除。由于平滑會導致邊緣的延展，因此邊緣檢測器只考慮那些具有局部梯度最大值的點為邊緣點。這一點可以用二階導數的零交叉點來實現(xiàn)。拉普拉斯函數用作二維二階的近似，是因為它是一種無方向算子。為了避免檢測出非顯著邊緣，應選擇一階導數大于某一閾值的零交叉點作為邊緣點 [6]。LoG 算子對圖像 進行邊緣檢測，輸出 是通過卷積運算得到的，即：??yxf, ??yxh, （式 )??yxfeyxh ,*2, 242?????????????????濾波（或平滑）、增強和檢測 3 個邊緣檢測的步驟對 LoG 算子邊緣檢測依然成立。其中高斯濾波器對圖像進行平滑，拉普拉斯算子將邊緣點轉換成零交叉點來實現(xiàn)，邊緣檢測通過零交叉點的檢測實現(xiàn)。1 高 斯 濾 波 實 質 上 是 一 種 信 號 的 濾 波 器 ， 其 用 途 是 信 號 的 平 滑 處 理 ， 我 們 知 道 數 字 圖 像 用 于 后 期 應 用 ， 其 噪聲 是 最 大 的 問 題 。1 0 11 2 10 0 01 2 110 區(qū)域生長法 區(qū)域生長原理區(qū)域生長的基本思想是將具有相似性質的像素集合起來構成區(qū)域。具體先對每個需要分割的區(qū)域找一個種子像素作為生長的起點，然后將種子像素周圍鄰域中與種子像素有相同或相似性質的像素（根據某些事先確定的生長或相似準則來判定）合并到種子像素所在的區(qū)域中。將這些新像素當作新的種子像素繼續(xù)進行上面的過程，直到再沒有滿足條件的像素可被包括進來。這樣一個區(qū)域就長成了 [7]。 區(qū)域生長的一個關鍵是選擇合適的生長或相似準則，大部分區(qū)域生長準則使用圖像的局部性質。生長準則可根據不同原則制定，而使用不同的生長準則會影響生長的過程。下面介紹 2 種基本的生長準則和方法。 灰度差準則區(qū)域生長方法將圖像以像素為基本單位來進行操作，基于區(qū)域灰度差的方法主要有如下步驟：1. 設灰度差的閾值 1為零，用上述方法進行區(qū)域擴張，使灰度相同像素合并。2．求出所以鄰接區(qū)域之間的平均灰度差，并合并具有最小灰度差的鄰接區(qū)域。3. 設定終止準則，通過反復進行上述步驟（2）中的操作將區(qū)域依次合并直到中指準則滿足為止