基于自適應(yīng)灰度分段線性變換的醫(yī)學(xué)C T序列圖像增強

西安思源學(xué)院工學(xué)院 王浩軍

第四軍醫(yī)大學(xué)口腔醫(yī)院 張 圃

西安思源學(xué)院工學(xué)院 楊 燕 屈瑞娜

基于自適應(yīng)灰度分段線性變換的醫(yī)學(xué)C T序列圖像增強

西安思源學(xué)院工學(xué)院 王浩軍

第四軍醫(yī)大學(xué)口腔醫(yī)院 張 圃

西安思源學(xué)院工學(xué)院 楊 燕 屈瑞娜

為了有助于CT圖像的臨床診斷，本文提出了一個自適應(yīng)的灰度分段線性變換方法，對醫(yī)學(xué)CT序列圖像進行了增強。該方法首先基于直方圖多尺度變換和EM(Expectation-Maximization)算法自適應(yīng)地選取分段線性變換的分段點，然后對背景和不感興趣的區(qū)域進行灰度范圍壓縮，對感興趣目標區(qū)域進行灰度范圍拉伸處理。實驗結(jié)果顯示感興趣目標區(qū)域?qū)Ρ榷缺缓侠碓鰪姡C明了該方法的有效性。

灰度分段線性變換；圖像增強；CT序列圖像；直方圖分析

醫(yī)學(xué)CT成像在臨床醫(yī)學(xué)中是一種無損傷的成像技術(shù)，CT圖像目前已成為醫(yī)學(xué)診斷中重要的手段之一。由于醫(yī)學(xué)CT圖像可針對不同的人體組織和病變呈現(xiàn)不同的像素灰度分布，從而有效地分辨人體組織內(nèi)微小的差別，使影像診斷的范圍大大擴大。特別是醫(yī)生通過CT序列圖像可了解和觀察組織器官的立體解剖結(jié)構(gòu)，分析和診斷病變組織及其定位。和其他醫(yī)學(xué)成像技術(shù)一樣，CT成像技術(shù)由于受成像設(shè)備、獲取條件等多種因素的影響，獲得的圖像質(zhì)量會出現(xiàn)退化，造成結(jié)構(gòu)邊界模糊不清，給醫(yī)生診斷帶來一些困難，因此圖像增強和其他圖像處理技術(shù)可作為輔助的圖像優(yōu)化技術(shù)，在醫(yī)學(xué)輔助診斷中起到了很大的作用。

圖像增強技術(shù)是針對給定圖像的應(yīng)用場合，改善圖像的視覺效果，有目的地強調(diào)圖像的整體或局部特性，將原來不太清晰的圖像變得清晰或更加強調(diào)某些感興趣目標特征，擴大圖像不同目標之間的差別，抑制不感興趣的區(qū)域，實現(xiàn)圖像質(zhì)量的改善，加強圖像的判讀和識別效果。通過對CT圖像的增強，突出了感興趣組織和病變的視覺特征，更適合醫(yī)生的明確診斷。圖像增強方法很多，根據(jù)處理過程所在的空間不同，可分為基于空域的算法和基于頻域的算法，前者直接對圖像所在空間進行處理，后者在圖像的變換域-頻域空間上對圖像進行處理。空域的圖像增強又分為點運算和鄰域運算，點運算的增強方法主要包括灰度變換，直方圖處理等，目的是使圖像成像均勻，擴大圖像的動態(tài)范圍，擴展圖像對比度。鄰域運算主要是以圖像平滑和銳化為目的，即消除圖像噪聲、突出目標邊緣輪廓，主要有均值濾波、中值濾波、高通濾波等算法。

本文針對CT序列圖像的特點，在基本的灰度分段線性變換方法的基礎(chǔ)上，提出一個基于直方圖多尺度變換和EM算法相結(jié)合的方法，自適應(yīng)選取灰度分段線性變換中的分段點，從而實現(xiàn)不同灰度段的自動變換，抑制了背景，提高了感興趣目標區(qū)域的對比度，最后還驗證了該算法的有效性。

一、基本的灰度分段線性變換方法

所謂圖像的灰度分段線性變換，是將圖像的灰度區(qū)間分成兩段以上，分別對這些分段的灰度區(qū)域進行相應(yīng)的灰度變換，例如將其中的某一段拉伸或壓縮到一定的灰度范圍。圖1為一典型的三分段的灰度線性變換示意圖，橫坐標S代表了原始的灰度，通過一個線性函數(shù)，將原始的灰度映射到縱軸T，即代表經(jīng)過分段線性變換后的灰度。可以看出通過這樣的變換，原圖中灰度值在0到s1和s2到L-1間的灰度范圍減少到了0到t1和t2到L-1，從而抑制了不感興趣區(qū)域的灰度級；而原圖中在s1和s2間的灰度范圍增加到了t1到t2，作用是拉伸了這個范圍的灰度級，增強了特征目標的對比度，這里L-1一般指最大的灰度值255。從圖中可以得出相應(yīng)的變換函數(shù)表達式：

圖1 灰度分段變換示意圖

在灰度分段線性變換的應(yīng)用中，確定圖1中折線的分段點是關(guān)鍵，圖中（s1，t1）和（s2，t2）的位置控制著折線的形狀。在該傳統(tǒng)方法的應(yīng)用中，分段點多采用人工確定，為了達到最佳的變換效果，往往需要反復(fù)地試驗，才能找到適當?shù)姆侄吸c。通過上述分析，對于CT圖像，如果能夠通過某種方式自動地找到背景、不感興趣區(qū)域和感興趣區(qū)域的灰度分段點，就可以通過分段線性變換進行合適的灰度壓縮和拉伸等，從而達到抑制背景，突出目標的目的。

二、自適應(yīng)灰度分段線性變換法

在本文增強的CT序列圖像中，背景分布在灰度低段，要突出的組織和病變主要分布在圖1的第2段，灰度值的高段主要由不感興趣的區(qū)域組成，如果采用一個自適應(yīng)的確定分段點的方法，則可以彌補傳統(tǒng)的灰度分段線性變換方法的不足。 自適應(yīng)分段點的計算方法很多，但針對CT序列圖像的還比較少見，考慮到CT圖像中感興趣的病變目標對比度低，為了有助于醫(yī)生的進一步診斷，文中提出了一個基于直方圖多尺度變換和EM算法相結(jié)合的方法，自適應(yīng)選取灰度分段變換中的分段點，從而實現(xiàn)不同灰度段的自動變換。該方法對CT序列圖像進行灰度線性變換方法流程如圖2所示。

（一）直方圖分析

圖像直方圖雖然不包含圖像目標的空間信息，但它可以提供圖像目標的灰度分布。通常最直接和最常用的方法，是把圖像的直方圖考慮成為一個高斯混合模型（GMM）(如公式（2）,（3）)，其中每一高斯分量代表圖像中一個目標類的灰度分布，后面可用EM算法去估計和優(yōu)化高斯混合模型參數(shù)，其中被估計的直方圖高斯混合模型表示為：

這里是被估計的歸一化直方圖，g變量代表圖像的灰度值，一般取值范圍是[0,255]，C是模型中高斯分量的數(shù)目，分別是第i個高斯分量的權(quán)值、均值和方差。第i個高斯分量的公式如下：

圖2 CT序列圖像自適應(yīng)灰度分段變換流程圖

（二）方圖多尺度變換和初始化參數(shù)確定

使用直方圖多尺度變換的目的是獲得高斯混合模型的初始參數(shù)，然后利用EM算法優(yōu)化這些參數(shù)，從而自適應(yīng)地找到分段線性變換的分段點。多尺度變換可以產(chǎn)生灰度直方圖的多尺度空間表達，以便提取相應(yīng)的特征點。一般來說，可利用一系列尺度不同的高斯函數(shù)與直方圖卷積去實現(xiàn)。在尺度空間中，高斯函數(shù)的尺度大小是由高斯函數(shù)的標準差σ體現(xiàn)的，σ隨尺度的增大而逐漸增大，卷積后的直方圖越平滑。在每個尺度下，通過直方圖二階導(dǎo)數(shù)零交叉點的性質(zhì)，確定直方圖的峰值特征。隨著尺度的增加，直方圖的峰點逐漸變少，有意義的峰點被保留，最后在尺度空間中，可自上而下精確的確定這些峰點位置。

在多尺度變換中，尺度的選擇是非常重要的。文中可通過分析CT序列圖像的相應(yīng)目標區(qū)域類別先驗的得到，按照灰度值范圍有小到大的次序，可以分為黑背景區(qū)，不感興趣的低灰度區(qū)，感興趣的目標區(qū)域和不感興趣的高灰度區(qū)四大類。這四類決定了直方圖中有意義的峰點數(shù)目，以及直方圖高斯混合模型的高斯分量數(shù)目，我們用該數(shù)目去控制最大的尺度大小σ。在多尺度變換中，一旦峰點數(shù)目＜=4尺度變換就停止，然后較容易的獲得上述4類區(qū)域的大致范圍及它們的高斯分量初始化參數(shù)。直方圖高斯混合模型的初始化參數(shù)計算如下：

式中是高斯混合模型第i個分量的均值，g是灰度值(0～255),LT和RT分別為峰點左右兩邊的谷底拐點對應(yīng)的灰度值，該拐點同樣可用類似于峰點查找的方法確定，不同的是其二階導(dǎo)數(shù)零交叉點的變化方向與峰點不同。

這里LI和HI分別是第i分量的最低灰度值和最高灰度值，它們可利用與均值點位置最近鄰的原則確定。

（三）EM算法優(yōu)化參數(shù)和分段點的確定

Duda和Hart強調(diào)，在EM算法中確定好的初始化參數(shù)是十分必要的，主要是為了更加逼近全局最優(yōu)解，而不是陷入局部最優(yōu)，因此上述的直方圖多尺度變換方法，就是為了幫助我們自動地獲得好的高斯分量的初始參數(shù)。在測試實驗中，上述4類區(qū)域利用2.2節(jié)的方法都能較好的被確定和定位，從而方便確定分段線性變換的分段點，完成相應(yīng)區(qū)域的圖像對比度加強，和其他不重要區(qū)域的對比度削弱。

EM算法優(yōu)化參數(shù)是一個迭代的過程，主要根據(jù)已有的數(shù)據(jù)獲得高斯混合模型的最大似然估計，當滿足一定的收斂準則，表明當前模型參數(shù)是最優(yōu)解。具體的原理和迭代實現(xiàn)可參考文獻。最后根據(jù)EM算法估計的最優(yōu)參數(shù)確定各高斯分量的分布及分界點，從而確定2.4節(jié)的分段點。圖3結(jié)果顯示一CT序列圖像的歸一化直方圖，以及相應(yīng)的高斯混合模型和各類分界點，其中自動的選取s1和s2灰度值作為我們分段線性變換的兩分段點的橫坐標值。

圖3 CT序列圖像歸一化直方圖及其對應(yīng)的高斯混合模型

（四）具體的分段線性灰度變換

在CT序列圖像中，感興趣的目標主要分布在圖3中第三個高斯分量所在的灰度區(qū)域，因此該區(qū)域的灰度將被拉伸，其他都將被壓縮。在實際變換時，只要給出對應(yīng)段的拉伸系數(shù)或壓縮系數(shù)，然后變換區(qū)域的動態(tài)范圍就可以被確定。圖1中兩個分段點將整個灰度范圍分成三段，其中第一段和第三段范圍將被壓縮，當這兩段被壓縮后，中間的灰度范圍自然被拉伸了。為了實現(xiàn)這樣的變換，就是要找到兩個分段點的橫坐標和縱坐標值，其中橫坐標s1和s2通過上節(jié)所述方法獲得，接下來就是計算它們的縱坐標t1和t2：

其中L-1是圖像灰度的最大值，一般取為255，k為圖一第一段和第三段的壓縮系數(shù)，該取值要小于1。這兩段壓縮系數(shù)取為相同的的目的，是為了簡化運算，同時不影響變換效果。當?shù)谝欢魏偷谌蔚幕叶确秶粔嚎s后，中間那段的灰度范圍自然被拉伸了，對應(yīng)的區(qū)域就實現(xiàn)了對比度加強。總之，當分段點確定后根據(jù)公式(1)，就可完成基于灰度的分段線性變換了。

三、實驗結(jié)果

利用提出的方法，針對連續(xù)的的4幅CT序列圖像進行了灰度分段線性變換實驗。圖4分別給出原始圖像，見圖中第一行，以及在不同壓縮系數(shù)取值k=0.2，k=0.4和k=0.7下的變換結(jié)果(分別顯示在第二行，第三行和第四行)，k值越小，壓縮的相對較大。這里變換空間的兩個分段點分別取值為（82，t1）和（198，t2）。

圖4 CT序列圖像及不同壓縮系數(shù)下的變換結(jié)果

該圖像中的病人為左側(cè)腮腺低分化癌，箭頭所指處為病變部位，原始圖像病變顯示部位和右側(cè)相同部位對比，圖像差異不是太明顯。通過該方法變換后，明顯可以看出，病變部位的灰度與相應(yīng)的正常組織相比，灰度值偏高，灰度分布不均勻。另外比較圖中第二行、第三行和第四行可以看出，k值過小，高低灰度區(qū)域的抑制過大，目標細節(jié)丟失過多,見圖中第二行；k值過大，高低灰度值區(qū)域的抑制不明顯,目標對比度增強不明顯，見圖中第四行；當k值選取合適時（見圖中第三行），在抑制背景的同時目標的細節(jié)得以保留，該區(qū)域?qū)Ρ榷缺挥行г鰪姡Y(jié)果將有助于臨床的影像診斷。

四、總結(jié)

本文提出了一個基于直方圖多尺度變換和EM相結(jié)合的自適應(yīng)灰度分段變換方法。根據(jù)應(yīng)用的圖像不同，可自適應(yīng)地選取變換分段點，對感興趣目標段、不感興趣區(qū)域段、背景可進行相應(yīng)的拉伸、壓縮等處理。另外還分析了方法中壓縮系數(shù)k的選取對結(jié)果的影響；實驗證明了該方法的有效性和實用性。今后在更廣泛的應(yīng)用中，如果被處理圖像灰度變化復(fù)雜，還可考慮對各分段區(qū)域?qū)嵤┓蔷€性灰度變換等處理。

[1]余建明．實用醫(yī)學(xué)影像技術(shù)[M]．人民軍醫(yī)出版社,北京,2015(10)．

[2]章毓晉．圖象處理和分析[M]．清華大學(xué)出版社,北京,1999(10)．

[3]張秀君,孫曉麗．分段線性變換增強的自適應(yīng)方法[J]．電子科技,2005,3:13-16．

[4]M．J．Carlotto．Histogram analysis using a scale-space approach．IEEE trans on PAMI,1987,9(1):121-129．

[5]R．O．Duda,P．E．Hart,D．G．Stork． Pattern Classification-2nd ed．[M]．Wiley&SONS,2001

[6]S．Theodoridis,K．Koutroumbas．Pattern Recognition-2nd ed．[M]．Elservier Science,2003．

王浩軍（1968-），女，河北寧晉人，博士，副教授，研究方向為圖像處理和模式識別。

張圃（1964-），男，陜西西安人，博士，副教授，研究方向為口腔頜面外科學(xué)。

楊燕（1971-），女，陜西西安人，碩士，高工，研究方向為控制工程。

屈瑞娜（1978-），女，河南平頂山人，本科，講師，研究方向為控制工程。

陜西省教育廳自然科學(xué)專項項目“基于增強CT的口腔癌計算機識別和輔助分析”（項目編號：2013JK1173）。