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

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

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

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

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

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

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

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

灰度分段線性變換；圖像增強(qiáng)；CT序列圖像；直方圖分析

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

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

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

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

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

圖1 灰度分段變換示意圖

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

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

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

（一）直方圖分析

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

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

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

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

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

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

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

這里L(fēng)I和HI分別是第i分量的最低灰度值和最高灰度值，它們可利用與均值點(diǎn)位置最近鄰的原則確定。

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

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

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

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

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

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

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

三、實(shí)驗(yàn)結(jié)果

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

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

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

四、總結(jié)

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

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

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

[3]張秀君,孫曉麗．分段線性變換增強(qiáng)的自適應(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-），女，河北寧晉人，博士，副教授，研究方向?yàn)閳D像處理和模式識(shí)別。

張圃（1964-），男，陜西西安人，博士，副教授，研究方向?yàn)榭谇活M面外科學(xué)。

楊燕（1971-），女，陜西西安人，碩士，高工，研究方向?yàn)榭刂乒こ獭?/p>

屈瑞娜（1978-），女，河南平頂山人，本科，講師，研究方向?yàn)榭刂乒こ獭?/p>

陜西省教育廳自然科學(xué)專項(xiàng)項(xiàng)目“基于增強(qiáng)CT的口腔癌計(jì)算機(jī)識(shí)別和輔助分析”（項(xiàng)目編號(hào)：2013JK1173）。