圖像信息是人類獲得外界信息的主要來源，數字圖像處理是利用計算機技術對現實生活中獲取的數字圖像進行處理，提取各目标的特征加以分析，最後對不同目标進行相應的自動操作。我們可以通過攝像機、數碼相機、掃描儀來獲取顆粒圖像，采集圖像的過程一般會受多外部因素的影響，如光照強度不均勻，元器件特性不穩定，其他顆粒物體幹擾或顆粒物體本身的光學特性等，使得采集到的圖像不利于後續處理。因此，需要對采集的原始灰度圖像進行一定的預處理來消除外界的幹擾，使得後續的處理能更加方便進行，圖像的預處理包括圖像灰度化、灰度圖像平滑處理、灰度圖像去背景、灰度圖像二值化、形态學相關處理、距離變換這幾部分。這一章對圖像灰度化、平滑處理、去背景操作、二值化等進行介紹。

1.1圖像處理基本概念

1.1.1數字圖像表示

客觀世界在空間上是三維（3-D）的，但一般從客觀景物得到的圖像是二維（2-D）的，一幅M×N的2-D圖像一般用一個二維數組f（x，y）來表示，如式1-1。x和y表示2-D空間XY中一個坐标點的位置，而f則代表圖像在坐标點（x，y）處某種性質F的數值。

                             1-1

灰度圖像中f表示灰度值，常對應客觀景物被觀察到的亮度；二值圖像中f的取值隻有兩個，分别對應黑和白；彩色圖像中f則一般用矢量f來表示，由于在每個圖像點同時具有紅綠藍三個值，可記爲，總之是要根據圖像内不同位置所具有的不同性質來利用圖像的。

1.1.2圖像類型與圖像灰度化

圖像類型中最常用有二值圖像、灰度圖像、彩色圖像這三種不同格式類型。下面就文中圖像處理算法所涉及的三種類型作簡單介紹。

1. 灰度圖像－Gray Scale Image

一幅灰度圖像就是一個數據矩陣，是數字圖像的基本形式，也就是我們所說的黑白照片，它隻有灰度顔色，沒有彩色。灰度即表示圖像像素明暗程度的數值，灰度級表明圖像中不同灰度的最大數量，灰度圖像像素的灰度級通常爲8Bits，範圍大小在0~255。“0”表示純黑色，“255”表示純白色，圖1-1表示了一幅灰度圖像及其對應的數組形式。現實中常用的是8Bits，但不斷的向10Bits、12Bits、16Bits發展。

圖1-1 灰度圖像表示

2. 二值圖像－Binary Image

二值圖像是指圖像的每個像素隻能是黑或者白，沒有中間過渡，二值圖像的像素值爲0、1，0表示黑，1表示白，因此也稱爲1-bit圖像。圖1-2表示了一幅二值圖像及其對應的數組形式。二值圖像隻能反映出圖像中物體的基本形狀，但其數據内容簡單，處理速度快，論文後面章節的顆粒參數提取都是基于二值圖像的。

圖1-2 二值圖像表示

3. 彩色圖像－Color Image

彩色圖像的表示方法有許多種，一般根據描述顔色的三維空間坐标系來進行分類。常用的顔色模型有：RGB（紅、綠、藍）、CMYK（青、品紅、黃、黑）、HSI（色調、色飽和度、亮度），計算機設備中最常用的是RGB彩色圖像模型。RGB彩色空間常用一個RGB彩色立方體加以圖解展示，在立方體的主對角線上，各原色的量相等，産生有暗到亮的白色，即灰度。如果像素的灰度級爲8Bits，則（0，0，0）爲黑，（255，255，255）爲白，正方體的其他6個角點分别爲紅、黃、綠、青、藍和品紅，結果如圖1-3所示。

    

                                       (a)RGB模型                     (b)RGB彩色立方體

圖1-3 RGB彩色空間

本文中最初采集到的顆粒圖像爲RGB彩色圖像，後期的相關處理都是基于灰度圖像和二值圖像的，所以首先要将彩色圖像轉換爲灰度圖像，具體的轉化方法有如下兩種。

1. 一幅RGB圖像就是彩色像素的一個M×N×3的數組，每一個像素點都是在特定空間位置的彩色圖像對應的紅、綠、藍三個分量，同時RGB也可以看成一個由三幅灰度圖像形成的“堆”，所以可直接将分量數據轉化爲對應的灰度圖像。其結果如圖1-4所示。

圖1-4 彩色圖像及其各分量灰度圖像

從上圖觀察出來，每個顔色分量包含白色區域，白色區域表示與每個顔色分量相關的最大值，如R分量中，白色代表純紅色區域部分，G、B分量同理。本文中以玉米顆粒圖像爲主要研究對象，對采集到的玉米顆粒圖像，分别表示出R、G、B分量，每個分量圖像均爲灰度圖像，結果如圖1-5所示。

 

圖1-5玉米圖像與各分量灰度圖像

由圖1-5可以看出，對于R分量的灰度圖像，顆粒目标物與背景對比度最高，利于後續的圖像分割，G分量次之，B分量目标物與背景最接近，效果也是最差的。很容易理解，原始圖像中，玉米顆粒的顔色很接近紅色，藍色分量最少。對于該幅圖像，我們完全可以使用其R分量作爲後續要處理的灰度圖像，但是由于顆粒圖像的多樣性，我們不可能對每幅圖像進行各個分量的判斷，然後選擇合适的分量灰度圖像，需要有一種綜合的方法來進行彩色圖像的灰度化。

2. 将上面提取出的各分量灰度圖像進行組合，通過一定的組合産生所需的灰度圖像，而且所産生的灰度圖像不會因原圖像中顔色的差異而有着明顯的變化。一般來說，将各個分量進行組合時，應滿足式1-2。

                                 

                         1-2

式中，r、g、b分别表示R、G、B各分量的系數，C爲最終的灰度圖像。根據已有的先驗知識，将各顔色分量做如下式1-3組合時，得到的灰度圖像最爲合适。

                                  1-3

将采集到的玉米顆粒圖像各分量做式1-3各比例組合，最終得到我們所需的灰度圖像，實驗結果如圖1-6所示。

圖1-6 灰度圖結果

對比圖1-6和1-5(b),第二種方法的效果要比第一種略差一些，但由于處理目标的多樣性，顔色的多樣性，第一種方法會大大降低算法的應用價值，利用第二種方法轉換的灰度圖像，完全滿足後續的相關算法處理，因此文本采用第二種方法來将彩色圖像轉化爲灰度圖像。

1.1.3圖像分割的定義

在對圖像進行分析和研究中，經常隻是對圖像中的某些部分感興趣，這些部分我們稱爲目标，其餘部分爲背景。爲了對圖像中的目标進行詳細分析，需要将這些目标從圖像中分離提取出來，然後才有可能進一步對目标進行測量和分析。

圖像分割定義爲：令集合R表示整個圖像區域，對R的分割可看作是把R分爲N個滿足以下條件的集合。

圖像分割在圖像處理中有着重要的位置，它是從圖像預處理到圖像分析和識别重要一步。首先圖像分割是表達目标的基礎，對後續的特征參數測量有重要的影響；其次，圖像分割以及後續的目标表達、特征參數提取測量等操作将原始圖像轉化成抽象緊湊的方式，使得後續有關的圖像識别、分析和理解成爲可能。

本文研究的圖像分割首先将顆粒目标物從背景中進行正确分割，對分割後的目标物出現粘連部分進行粘連分割，最終得到單獨的顆粒目标區域，爲統計每個顆粒特征參數打下很好的基礎。

1.2 灰度圖像預處理

灰度圖像預處理屬于圖像處理中底層變換的操作，是根據要求在空域對圖像像素進行的變換，在實際的運算中，常用模闆運算來進行相關處理。圖像預處理相關操作一般是用來抑制或消除對相關圖像處理和分析不利或者無關的信息，爲後續相關處理提供良好的高質量的圖像。這裏對顆粒灰度圖像常用的相關的預處理：圖像平滑處理、圖像背景去除、圖像灰度修正這三個方面進行簡要的介紹。

1.2.1 圖像平滑處理

圖像平滑處理也稱爲圖像平滑濾波，平滑濾波能夠減弱或消除圖像中的高頻部分，但不影響低頻分量（高頻部分對應圖像中的區域邊緣等灰度值有較大較快變化的部分，低頻分量對應圖像中灰度值緩慢變化的區域），在将這些分量濾去後可減少局部灰度起伏，使圖像變得比較平滑。由于噪聲空間的相關性較弱，對應的是較高的空間頻率，所以在實際中，平滑濾波很重要的一個功能就是消除噪聲。噪聲消除在空域可分爲線性平滑濾波和非線性平滑濾波。線性平滑濾波中常用的是均值濾波和高斯濾波；中值濾波是非線性濾波器的代表。實際的很多平滑濾波算法均是基于模闆的卷積進行的，下面就對這三種常用濾波方法進行簡單介紹。

1. 均值濾波

均值濾波是通過一點和鄰域内像素點求平均來去除突變的像素點，從而濾掉一定的噪聲，其數學表達式如式1-4所示。

                                          1-4

其中，N（x,y）是對應f（x，y）中（x，y）處n×n的鄰域，與模闆所覆蓋的範圍對應。常用的3×3和5×5的模闆如下圖1-6所示。

圖1-6 均值濾波模闆

該算法主要的特點是算法簡單，計算速度快，但其代價是會造成圖像一定程度上的模糊。

2. 高斯濾波

高斯濾波是根據高斯分布來設計的濾波器，由于高斯函數有着一些良好的特性，對二維連續高斯分布經采樣、量化，并使模闆歸一化，便可得到二維高斯離散模闆，常用的3×3和5×5的模闆如下圖1-7所示，可以看出來，高斯模闆是一種加權模闆，并且是按二維正态分布進行加權的。

圖1-7 高斯濾波模闆

3. 中值濾波

中值濾波[34]是在非線性濾波中基于排序的一種濾波方法，是依靠模闆排序來實現的，用于圖像的2-D中值濾波的輸出如式1-5所示。

                       1-5

式中，median代表取中值，即對模闆覆蓋的信号序列按照數值大小進行排序，并取排序後處在中間位置的值，N（x,y）對應f（x，y）中（x，y）的n×n的鄰域。對于一個尺寸爲n×n的中值濾波模闆，其輸出值應大于等于模闆中個像素的值，并且小于等于模闆中個像素的值。一般情況下，圖像中尺寸小于模闆一般的過亮或者過暗區域将會在濾波後被消除。可以得出中值濾波的主要作用就是就是讓周圍的像素灰度值的差在比較大的情況下，變換成與周圍像素值接近，這樣就能對單獨孤的噪聲像素有着很強的消除能力，由于中值濾波與均值濾波的原理不同，所以很少産生的模糊情況，也就是說，中值濾波既能消除噪聲又能較好地保持圖像的細節。

中值濾波的效果不僅與所有的模闆尺寸有關，還與模闆中參與運算的像素個數和所用模闆中參與運算的像素構成的圖案有關。一般而言，十字叉模闆保留細的水平線和垂直線，但會濾掉對角線；方形的模闆對圖像的細節最不敏感，會濾除細線并消除邊緣上的點；X形狀的模闆僅保留對角線，常用的一些中值濾波模闆如圖1-8所示。

圖1-8 常用中值濾波模闆

這裏對大米顆粒圖像進行去噪聲實驗，結果如圖1-9所示。圖(d)、(e)、(f)分别表示用3×3模闆進行均值濾波、高斯濾波、中值濾波去除高斯噪聲後的圖像，圖(g)、(h)、(i)分别表示用3×3模闆進行均值濾波、高斯濾波、中值濾波去除椒鹽噪聲後的圖像。

 

圖1-9 噪聲平滑處理

觀察上面實驗結果可知，中值濾波對于噪聲的濾除有着很好地效果，對于采集到的顆粒圖像，一般采取中值濾波來進行相關處理，都會取得很好的平滑結果。

1.2.2 圖像去背景處理

采集圖像時，由于光照的不均勻和背景的反光特性，使得獲得的顆粒圖像往往會出現背景光照不均勻的現象，有的時候背景和目标顆粒有着相同的灰度值，這個時候，如果直接對這樣背景亮度不均勻的圖像做阈值處理會很困難。如圖1-10所示，其中1-10(b)是阈值處理後的圖像。

  

圖1-10 不均勻背景顆粒阈值分割

觀察1-10(b)可以看出，圖像頂部的米粒已被很好地從背景中分離出來，但是圖像底部的米粒并未從背景中正确提取出來。可以通過開運算來補償不均勻的背景亮度，隻要結構元素大到不能完全匹配米粒，對圖像進行開運算可以産生對整個圖像背景的合理估計，産生的背景可以認爲是顆粒圖像的背景，結果如圖1-11(a)所示，然後将圖1-10(a)減去圖1-11(a)，就可以得到具有合适且均勻背景的顆粒圖像，結果如圖1-11(b)所示。

  
    

 

圖1-11 去除不均勻背景

1.2.3 圖像灰度修正

圖像灰度修正是圖像增強技術的一種，适當的修正方法，可以使得原本模糊不清甚至無法分辨的原始圖像處理成清晰并包含大量有用信息的可使用圖像。文中主要使用基于灰度變換的圖像灰度修正。

圖像的灰度變換是基于空間域的圖像處理方法，通過按照一定的變換關系逐點改變原圖中每個像素灰度值的方法，數學表達式如式1-6所示。

                               g(x,y)=T[f(x,y)]                       1-6 

式中，f(x,y)爲變換前原圖像，其灰度範圍爲[a,b]，處理後的圖像爲g(x,y)，其灰度範圍是[c,d]。常見的T有線性和非線性兩種，線性變換可将灰度值區域進行線性放大和縮小，或者是進行分段變換，非線性變換可以使圖像灰度的分布與人的視覺匹配。

這裏以圖1-11(b)爲例進行圖像灰度修正，可以看出來圖1-11(b)灰度對比度比較低，對其進行灰度對比度增強，結果如圖1-12所示，可以看出來，圖像的質量得到了很大的提高，爲方便觀察，将兩幅圖片放置一起做比較。

  

 

圖1-11 圖像灰度修正效果圖

1.3圖像二值化

1.3.1 二值化基本概念

圖像二值化是将灰度圖像轉化爲隻有包含0和1的二值圖像的過程，二值圖像數據量少，處理速度快，能夠很好地保留圖像的外部特征（周長、面積、外觀比等），是顆粒圖像進行分析處理的最常用的圖像格式之一。

灰度圖像二值化一般是通過判斷灰度圖像中每個像素點的特征屬性是否滿足阈值要求，來确定圖像中的該像素點是屬于目标區域還是背景區域，從而将一幅灰度圖像轉換成二值圖像。最簡單的劃分方式是選擇特定的阈值T，按式1-5将灰度圖像f(x,y)轉化爲二值圖像g(x,y)。

                                                  1-5

阈值分割計算簡單，能用封閉且連通的邊界定義不交疊的區域，對目标與背景有較強對比度的圖像顯示出較好的分割效果，關鍵就是阈值T的确定，如果值選取過高，則過多的目标區域将會被劃分成背景；反之，如果阈值選擇過低，則過多的背景就會被劃分爲目标區。前文中，通過對顆粒圖像進行濾波和去背景等處理，已經對圖像質量進行了明顯地改善，采用最大方差阈值法能夠滿足絕大多數圖像阈值化分割，通過最大方差阈值法确定灰度圖像二值化阈值，依此阈值進行圖像二值化，下面對該方法算法進行簡要的介紹。

最大方差阈值分割方法也叫Otsu法，是一種使用類間方差最大的自動确定阈值的方法。設圖像中像素總數是n，是灰度級爲的像素數目，L是圖像中所有可能的灰度級，則的概率密度函數如式1-6所示。

                                      1-6

假設我們現在已經選定了一個阈值k，是一組灰度級爲的像素，是一組灰度級爲的像素。Otsu方法選擇最大類間方差的阈值k，類間方差定義如式1-7所示：

                                          1-7

其中各變量如式1-8所示，

                                                      1-8

令T在[0,L-1]範圍内，以步長1依次遞增取值，當最大時對應的T即爲最佳阈值。确定最佳阈值T後，根據式1-5就可得到所需的二值化圖像。

1.3.2 顆粒圖像二值化

使用最大方差阈值分割方法确定的阈值T具有廣泛的适用性，絕大多數顆粒圖像分析中圖像二值化都是采用該方法，下面對兩組圖片進行二值化處理，其中一幅是上文介紹的去除不均勻背景後再經過灰度值修正的大米顆粒，另外一幅圖片是作爲文中後續主要處理對象的玉米顆粒圖像，二值化結果如圖1-12所示。可以看出，該方法對顆粒圖像有着很好的二值化效果。

  

圖1-12 顆粒圖像二值化

1.4本章小結

本章首先介紹圖像處理中的數字圖像與像素表示，圖像類型等基本概念；另外詳細闡述圖像預處理過程中，兩種彩色圖像轉化爲灰度圖像的方法，對比分析後采用第二種方法來進行彩色圖像灰度化處理。接下來采用中值濾波去除噪聲幹擾，對背景光照不均勻圖像去除其不均勻背景，增加圖像灰度對比度，改善圖像質量。最後利用最大類間方差法選取合适阈值，對圖像進行二值化處理，得到滿意的二值圖像，爲後續處理提供良好的圖像。