數(shù)字圖像處理及實現(xiàn)

數(shù)字圖像處理及實現(xiàn)第一頁，共七十頁，編輯于2023年，星期三Contents

8.1圖像分割定義

8.2使用閾值進行圖像分割

8.3基于梯度的圖像分割方法8.4邊緣檢測和連接第二頁，共七十頁，編輯于2023年，星期三Contents

8.5區(qū)域增長(RegionGrowing)

8.6二值圖像處理(BinaryImageProcessing)

8.7分割圖像的結構小結第三頁，共七十頁，編輯于2023年，星期三8.1圖像分割定義

圖像分割處理定義：將數(shù)字圖像劃分成互不相交（不重疊）區(qū)域的過程.區(qū)域(region)定義:像素的連通集。連通(connectedness)定義:在一個連通集中的任意兩個像素之間，存住一條完全由這個集合的元素構成的連通路徑。第四頁，共七十頁，編輯于2023年，星期三

8.2.1全局閾值化思想：整個圖像中將灰度閾值的值設置為常數(shù)。前提條件：如果背景的灰度值在整個圖像中可合理地看作為恒定，而且所有物體與背景都具有幾乎相同的對比度，那么，只要選擇了正確的閾值，使用一個固定的全局閾值一般會有較好的效果。第五頁，共七十頁，編輯于2023年，星期三8.2.2自適應閾值改進方法：在許多的情況下，背景的灰度值并不是常數(shù)，物件和背景的對比度在圖像中也有變化，這時，一個在圖像中某一區(qū)域效果良好的閾值在其它區(qū)域卻可能效果很差。在這種情況下，把灰度閾值取成一個隨圖像中位置緩慢變化的函數(shù)值是適宜的。第六頁，共七十頁，編輯于2023年，星期三8.2.3最佳閾值的選擇除非圖像中的物體有陡峭的邊沿，否則灰度閾值的取值對所抽取物體的邊界的定位和整體的尺寸有很大的影響。這意味著后續(xù)的尺寸（特別是面積）的測量對于灰度閾值的選擇很敏感。由于這個原因，我們需要一個最佳的，或至少是具有一致性的方法確定閾值。第七頁，共七十頁，編輯于2023年，星期三1.直方圖技術含有一個與背景明顯對比的物體的圖像其有包含雙峰的灰度直方圖第八頁，共七十頁，編輯于2023年，星期三直方圖生成a=imread('d:\pic\i_boat_gray.bmp');imshow(a)figureimhist(a)利用灰度閾值T對物體面積進行計算的定義是:第九頁，共七十頁，編輯于2023年，星期三2.最大類間方差法(OTSU)OTSU算法定義：該算法是在灰度直方圖的基礎上用最小二乘法原理推導出來的，具有統(tǒng)計意義上的最佳分割閾值。第十頁，共七十頁，編輯于2023年，星期三OTSU基本原理：以最佳閾值將圖像的灰度直方圖分割成兩部分，是兩部分之間的方差取最大值，即分離性最大。第十一頁，共七十頁，編輯于2023年，星期三3.迭代法求閾值原理：圖像中前景與背景之間的灰度分布為相互不重疊，在該前提下，實現(xiàn)對兩類對象的閾值分割方法。第十二頁，共七十頁，編輯于2023年，星期三8.3基于梯度的圖像分割方法思路對比：

區(qū)域分割方法：通過將圖像劃分為內(nèi)部點集和外部點集來實現(xiàn)分割。邊界方法：利用邊界具有高梯度值的性質(zhì)直接把邊界找出來。第十三頁，共七十頁，編輯于2023年，星期三8.3.1邊界跟蹤

算法步驟1：我們從一個梯度幅值圖像著手進行處理，這個圖像是從一幅處于和物體具有反差的背景中的單一物體的圖像進行計算得來的。2:搜索以邊界起始點為中心的3×3鄰域，找出具有最大灰度級的鄰域點作為第2個邊界點。第十四頁，共七十頁，編輯于2023年，星期三8.3.2梯度圖像二值化如果用適中的閾值對一幅梯度圖像進行二值化，Kirsch的分割法利用了這種現(xiàn)象。算法步驟用一個中偏低的灰度閾值對梯度圖像進行二值化從而檢測出物體和背景，物體與背景被處于閾值之上的邊界點帶分開。隨著閾值逐漸提高，就引起物體和背景的同時增長。當它們接觸上而又不至于合并時，可用接觸點來定義邊界。這是分水嶺算法在梯度圖像中的應用。

第十五頁，共七十頁，編輯于2023年，星期三Kirsch的分割算法第十六頁，共七十頁，編輯于2023年，星期三8.3.3拉普拉斯邊緣檢測

拉普拉斯算于是對二維函數(shù)進行運算的二階導數(shù)標量算子。它定義為:第十七頁，共七十頁，編輯于2023年，星期三一個無噪聲圖像具有陡峭的邊緣，可用拉普拉斯算子將它們找出來。對經(jīng)拉普拉斯算子濾波后的圖像用零灰度值進行二值化會產(chǎn)生閉合的、連通的輪廓并消除了所有的內(nèi)部點。但是由于噪聲的存在，在運用拉普拉斯算子之前需要先進行低通濾波。第十八頁，共七十頁，編輯于2023年，星期三選用高斯低通濾波器進行預先平滑是很合適的。由卷積的結合律可以將拉普拉斯算子和高斯脈沖響應組合成一個單一的高斯拉普拉斯核：這個脈沖響應對x和y是可分離的，因此可以有效地加以實現(xiàn)。第十九頁，共七十頁，編輯于2023年，星期三8.4邊緣檢測和連接邊緣點：確定圖像中的物體邊界的另一種方法是先檢測每個像素和其直接鄰域的狀態(tài)，以決定該像素是否確實處于一個物體的邊界上。具有所需特性的像素被標為邊緣點。邊緣圖像:當圖像中各個像素的灰度級用來反映各像素符合邊緣像素要求的程度時，這種圖像被稱為邊緣圖像。第二十頁，共七十頁，編輯于2023年，星期三8.4.1邊緣檢測邊緣檢測算子檢查每個像素的鄰域并對灰度變化率進行量化，通常也包括方向的確定。有若干種方法可以使用，其中大多數(shù)是基于方向?qū)?shù)掩模求卷積的方法。第二十一頁，共七十頁，編輯于2023年，星期三Roherts邊緣算子

其中，f（x，y）是具有整數(shù)像素坐標的輸入圖像。其中的平方根運算使該處理類似于人類視覺系統(tǒng)中發(fā)生的過程。第二十二頁，共七十頁，編輯于2023年，星期三Sobel邊緣算子

兩個卷積核形成了Sobel邊緣算子。圖像中的每個點都用這兩個核做卷積。一個核對通常的垂直邊緣響應最大而另一個對水平邊緣響應最大。兩個卷積的最大值作為該點的輸出值。運算結果是一幅邊緣幅度圖像。第二十三頁，共七十頁，編輯于2023年，星期三Sobel邊緣算子圖第二十四頁，共七十頁，編輯于2023年，星期三Prewitt邊緣算子

Prewitr邊緣算子第二十五頁，共七十頁，編輯于2023年，星期三Kirsch邊緣算子

圖像中的每個點均與這8個模板進行卷積，每個掩模對某個特定邊緣方向作出最大響應。所有8個方向中的最大值作為邊緣幅度圖像的輸出。最大響應掩模的序號構成了對邊緣方向的編碼。Kirsch算子的梯度幅度值第二十六頁，共七十頁，編輯于2023年，星期三

Kirsch邊緣算子第二十七頁，共七十頁，編輯于2023年，星期三

Kirsch邊緣算子第二十八頁，共七十頁，編輯于2023年，星期三邊緣檢測器性能:使用兩個掩模板組成邊緣檢測器時，通常取較大得幅度作為輸出值。這使得它們對邊緣的走向有些敏感。取它們的平方和的開方可以獲得性能更‘致的全方位響應。這與真實的梯度值更接近。第二十九頁，共七十頁，編輯于2023年，星期三邊緣提取方法原圖第三十頁，共七十頁，編輯于2023年，星期三

邊緣提取方法邊緣提取后第三十一頁，共七十頁，編輯于2023年，星期三8.4.2邊緣連接問題：如果邊緣很明顯，而且噪聲級低，那么可以將邊緣圖像二值化并將其細化為單像素寬的閉合連通邊界圖。然而在非理想條件下，這種邊緣圖像會有間隙出現(xiàn)，需要加以填充。第三十二頁，共七十頁，編輯于2023年，星期三填充小的間隙可以簡單地實現(xiàn)，通過搜索一個以某端點為中心的5×5或更大的鄰域，在鄰域中找出其它端點并填充上必要的邊界像素，從而將它們連接起來。對具有許多邊緣點的復雜場景，這種方法可能會對圖像過度分割。為了避免過度的分割，可以規(guī)定：兩個端點只有在邊緣強度和走向相近的情況下才能連接。第三十三頁，共七十頁，編輯于2023年，星期三解決方法啟發(fā)式搜索曲線擬合Hough變換第三十四頁，共七十頁，編輯于2023年，星期三(1)啟發(fā)式搜索

假定在一幅邊緣圖像的某條邊界上有一個像間隙的缺口，但是這個缺口太長而不能僅用一條直線填充，它還可能不是同一條邊界上的缺口，可能在兩條邊界上。作為質(zhì)量的度量，我們可以建立一個可以在任意連接兩端點（稱為A，B）的路徑上進行計算的函數(shù)。。

第三十五頁，共七十頁，編輯于2023年，星期三

如果邊緣質(zhì)量函數(shù)很復雜而且要評價的缺口既多又長，啟發(fā)式搜索技術的計算會很復雜。這樣的技術在相對簡單的圖像中性能很好，但不一定能找出兩端點間的全局最佳路徑。第三十六頁，共七十頁，編輯于2023年，星期三(2)曲線擬合假定有一組散布在兩個特定邊緣點A和B之間的邊緣點，我們希望從中選取一個子集作為從A到B一條分段線性路徑上的結點集。首先：從A到B引一條直線。其次：接著計算其它的每個邊緣點到該直線的垂直距離。第三十七頁，共七十頁，編輯于2023年，星期三(3)Hough變換

直線y=mx+b可用極坐標表示為其中定義了一個從原點到線上最近點的向量。這個向量與該直線垂直。第三十八頁，共七十頁，編輯于2023年，星期三如果有一組位于由參數(shù)確定的直線上的邊緣點，則每個邊緣點對應了空間的一條正弦型曲線。所有這些曲線必交于點，因為這是它們共享的一條直線的參數(shù)。建立一個在空間的二維直方圖。對每個邊緣點，我們將給所有與該點的Hough變換（正弦曲線）對應的空間的直方圖方格一個增量。當對所有邊緣點施行完這種操作后，包含的方格將具有局部最大值。然后對空間的直方圖進行局部最大值搜索可以獲得邊界線段的參數(shù)。第三十九頁，共七十頁，編輯于2023年，星期三hough直線檢測結果

原圖直線檢測結果第四十頁，共七十頁，編輯于2023年，星期三8.5區(qū)域增長

方法：從把一幅圖像分成許多小區(qū)域開始的。這些初始的區(qū)域可能是小的鄰域甚至是單個像素。在每個區(qū)域中，對經(jīng)過適當定義的能反映一個物體內(nèi)成員隸屬程度的性質(zhì)（度量）進行計算。用于區(qū)分不同物體內(nèi)像素的性質(zhì)（度量）包括平均灰度值，紋理，或顏色信息。第四十一頁，共七十頁，編輯于2023年，星期三區(qū)域增長算法主要分成兩類(1)簡單連接這是基于單個像素的區(qū)域增長法，它從滿足的檢測的點（連接核）開始，考察其周圍(4鄰域或8鄰域)的不屬于任何一個區(qū)域的點，如果其特性符合接收準則，就把它作為同一個區(qū)域加以合并，形成連接核，繼而檢測周圍的點，并把符合接入準則的點并入，產(chǎn)生新的核。重復上述過程，直到?jīng)]有可并入的點時，生產(chǎn)過程結束。

第四十二頁，共七十頁，編輯于2023年，星期三（2）子區(qū)合并法合并過程：首先：將圖像分割成個，大小為的小區(qū)域（簡稱子區(qū)）。其次：從左上角第一個子區(qū)開始，分別計算子區(qū)和相鄰子區(qū)的灰度統(tǒng)計量，然后做相似性判定。若兩者的灰度分布相似且符合接收準則。相鄰子區(qū)并入當前子區(qū)，形成下一輪判定合并時的當前子區(qū)。第四十三頁，共七十頁，編輯于2023年，星期三如果某個相鄰子區(qū)不符合接收準則，就打上“未分割標記”。繼續(xù)新一輪判定，使當前子區(qū)不斷“生長”，知道沒有可歸并的子區(qū)為止，一個區(qū)域分割完畢。最后：搜索圖像全域，對凡具有“未分割標記”的子區(qū)重復上述步驟。第四十四頁，共七十頁，編輯于2023年，星期三相似性判定準則一般是：第四十五頁，共七十頁，編輯于2023年，星期三8.6二值圖像處理

二值圖像也就是只具有兩個灰度級的圖像，它是數(shù)字圖像的一個重要子集。一個二值圖像（例如，一個剪影像或一個輪廓圖）通常是由一個圖像分割操作產(chǎn)生的。如果初始的分割不夠令人滿意，對二值圖像的某些形式的處理通常能提高其質(zhì)量。第四十六頁，共七十頁，編輯于2023年，星期三形態(tài)學圖像處理第四十七頁，共七十頁，編輯于2023年，星期三術語1.集合論術語(Definition)形態(tài)學處理語言中，二值圖像B和結構元素S都是定義在二維笛卡兒網(wǎng)格上的集合，“1”是這些集合中的元素。當一個結構元素的原點位移到點(x，y)處時，我們將其記作。形態(tài)學運算的輸出是另一個集合，這個運算可用一個集合論方程來確定。第四十八頁，共七十頁，編輯于2023年，星期三2.腐蝕和膨脹(DilationandErosion)第四十九頁，共七十頁，編輯于2023年，星期三(1)腐蝕一般意義的腐蝕概念定義為：

E=BS={x，y|B}如果S的原點位移到點(x，y)，那么S將完全包含于B中。使用基本的3×3結構元素時，一般意義的腐蝕簡化為簡單腐蝕。第五十頁，共七十頁，編輯于2023年，星期三(2)膨脹一般膨脹定義為：

E=BS={x，y|B}S對B膨脹產(chǎn)生的二值圖像D是由這樣的點（x，y）組成的集合，如果S的原點位移到(x，y)，那么它與B的交集非空。采用基本的3×3結構造元素時，一般膨脹簡化為簡單膨脹。第五十一頁，共七十頁，編輯于2023年，星期三8.6.2開運算和閉運算開運算：先腐蝕后膨脹的過程稱為開運算。它具有消除細小物體、在纖細點處分離物體、和平滑較大物體的邊界時又不明顯改變其面積的作用。開運算定義為：第五十二頁，共七十頁，編輯于2023年，星期三閉運算：先膨脹后腐蝕的過程稱為閉運算。它具有填充物體內(nèi)細小空洞、連接鄰近物體、在不明顯改變物體面積的情況下平滑其邊界的作用。閉運算定義為：第五十三頁，共七十頁，編輯于2023年，星期三8.6.3腐蝕和膨脹的變體

通常反復施以腐蝕運算，將使一個物體變得不存在。類似地，反復膨脹將把一幅圖像中的所有物體合并為一個。然而，這些過程可以改變一下，以便在一些應用中產(chǎn)生更合適的結果。第五十四頁，共七十頁，編輯于2023年，星期三(1)收縮定義：當腐蝕以一種不觸及單像素物件的方式執(zhí)行時，這個過程稱為收縮。收縮可以迭代方式為一個包含近似圓形物體的二值圖像生成物體尺寸的分布。收縮時會使非常不圓的物體（如啞鈴狀的物體）分解，因此這種技術存它的局限性。第五十五頁，共七十頁，編輯于2023年，星期三(2)細化定義:將一個曲線形物體細化為一條單像素寬的線，從而圖形化地顯示出其拓撲性質(zhì)。腐蝕可編程過程：第一步：足一個正常的腐蝕。第二步：只將那些消除后并不破壞連通性的點消除，否則保留。第五十六頁，共七十頁，編輯于2023年，星期三(3)抽骨架定義：一個與細化有關的運算，也稱為中軸變換(Medialaxistransform)或焚燒草地技術(grass-firetechnigue)。抽骨架的實現(xiàn)與細化相似．可采用一個兩步有條件腐蝕實現(xiàn)，但是刪除像素的規(guī)則略有不同。第五十七頁，共七十頁，編輯于2023年，星期三(4)剪裁細化和抽骨架過程會在所生成的阿中留下毛刺。毛刺是由邊界上單像素尺寸的起伏造成的，這些起伏產(chǎn)生了小的分支。它們可通進一系列的消除端點的3×3運算（導致所有的分支縮短）除去，然后再重建那些留下的分支。第五十八頁，共七十頁，編輯于2023年，星期三(5)粗化一些分割技術傾向于用緊貼的邊界擬臺物體以避免錯誤地合并它們。通常孤立物體的最佳邊界太緊貼并不利于后續(xù)測量。粗化可在不合并彼此分離的物體的前提下擴大邊界，從而修正了這種不足。第五十九頁，共七十頁，編輯于2023年，星期三

原圖腐蝕第六十頁，共七十頁，編輯于2023年，星期三膨脹開操作第六十一頁，共七十頁，編輯于2023年，星期三

閉操作第六十二頁，共七十頁，編輯于2023年，星期三8.7分割圖像的結構希望制作一幅新圖像，以顯示物體是如何調(diào)整的，或者用單獨的圖像顯示每個物體。甚至還可能希望對單個物體逐個地進行進一步的測量或其它處理。在這些情況下，就需要抽取并以更方便的形式存儲各個物體。通常，每個物體在被檢測時都應該標以一個序號。這個物體編號可用來識別和跟蹤景物中的物體。在這一節(jié)，我們討論三種對分割圖像進行結構化的方法。第六十三頁，共七十頁，編輯于2023年，星期三8.7.1分水嶺算法

分水嶺算法第六十四頁，共七十頁，編輯于2023年，星期三分水嶺算法不是簡單地將圖像在最佳灰度級進行閾值處理，而是從一個偏低但仍然能正確分割各個物體的閾值開始。然后隨著閾值逐漸上升到最佳值，使各個物體不會被合并。這個方法可以解決那些由于物體靠得太近而不能用全局閾值解決的問題。如果初始的閾值太低，那么低對比度的物體開始時會被丟失，然后隨著閾值的增加就會和相鄰的物體合并。如果初始閾值太高，物體一開始便會被合并。最終的閾值決定了最后的邊界與實際