前言：本站為你精心整理了小波分析法在醫(yī)學(xué)中的運用范文，希望能為你的創(chuàng)作提供參考價值，我們的客服老師可以幫助你提供個性化的參考范文，歡迎咨詢。

1基于小波變換的醫(yī)學(xué)圖像融合的提出

傳統(tǒng)的融合方法多是在時間域?qū)τ跋襁M行算術(shù)運算,沒對處理影像時相應(yīng)頻率的變化進行考慮.不同的方法在頻率域有不同的表現(xiàn),這種頻率域表現(xiàn)是否合適,在一定程度上說明了融合方法是否有效.如果沒有充分考慮時間域和頻率域的互動關(guān)系,融合的效果會受到限制.對于影像融合,在頻率域進行比在時間域進行更為有效,融合算法的設(shè)計必須把融合技術(shù)目的和影像的頻率域表現(xiàn)(即融合理論基礎(chǔ))結(jié)合起來考慮.小波變換的多分辨率分析是當(dāng)前信息和圖像處理等領(lǐng)域的研究熱點,小波變換可將原始圖像分解成一系列具有不同空間分辨率和頻域特性的子圖像,可充分反映原始圖像的局部化特征,為醫(yī)學(xué)圖像數(shù)據(jù)融合提供了有利條件.CT和MRI圖像是臨床用于顱骨疾病診斷最有效的方法.對于同一身體部位在CT成像過程中,骨組織將產(chǎn)生強信號,也就對應(yīng)著灰度圖像上的較亮區(qū)域,而在MRI成像過程中,由于骨組織缺乏活動的質(zhì)子,因此產(chǎn)生的信號較其他弱,也就對應(yīng)著灰度圖像上的較暗區(qū)域,而軟組織卻清晰.目前臨床多使用T1加權(quán)和T2加權(quán)脈沖序列進行MRI成像.在CT,MRI的T1加權(quán)像和T2加權(quán)像中,圖像的灰度按從高到低的順序排列.小波分析方法為CT提供的骨信息和MR提供的關(guān)于軟組織的信息都能在最后的融合圖像中有機地結(jié)合在一起提供了有利的數(shù)學(xué)工具.

2基于小波變換的融合和結(jié)果定量評價

Mallat提出了小波變換的快速分解與重構(gòu)算法[4,5],利用兩個一維濾波器實現(xiàn)對二維圖像的快速小波分解,再利用兩個一維重構(gòu)濾波器實現(xiàn)圖像的重構(gòu).小波變換的目的是將原始圖像分別分解到一系列頻率信道中,利用其分解后的塔形結(jié)構(gòu),對不同分解層,不同頻帶分別進行融合處理,可有效地將來自不同圖像的細節(jié)融合在一起.融合時,將被融合圖像各自攜帶的不同特征與細節(jié)融合在多個分解層,多個頻帶上分別從不同景象進行融合.Campbell和Robson[6]的實驗表明,人的視網(wǎng)膜是在不同的頻率信道中進行處理的.基于小波分解的圖像融合也是在不同的頻率信道上進行融合處理的,可獲得與人的視覺特性更為接近的融合效果.這里采用了CT與MRI的融合實驗.設(shè)A,B分別為兩幅原始圖像(CT和MRI),F為融合后的圖像.在圖像融合過程中,融合規(guī)則及融合算子的選擇對于融合的質(zhì)量至關(guān)重要.也是圖像融合至今尚未能解決的難點問題.這里給出對不同融合規(guī)則和融合算子及其融合結(jié)果的比較:⑴最大值融合規(guī)則對小波分解中的低頻分量LL,高頻分量LHi,HLi,HHi取兩幅圖像中相應(yīng)系數(shù)矩陣中對應(yīng)項的(i=1,2,3,…,N)的最大值[7].⑵局部平均梯度融合規(guī)則分別計算對CT和MRI圖像的小波分解后的高頻子圖像上每個像素的局部平均梯度,以將像素的局部平均梯度作為準(zhǔn)則,確定融合后的高頻子圖像上的像素值.設(shè)A(x,y)和B(x,y)分別為CT和MRI圖像,不同分辨率高頻子圖像分別為Akj(x,y)和Bkj(x,y)(k=1,2,3,j為尺度參數(shù)),不同分辨率的高頻子圖像的梯度圖像為GAkj(x,y)和GBkj(x,y),若GAkj(x,y)≥GBkj(x,y),則小波分解后融合的高頻信息為Fkj(x,y)=K1×GAkj(x,y)+K2×GBkj(x,y),(3)若GAkj(x,y)<GBkj(x,y),則K1<K2.⑶局部區(qū)域能量的融合規(guī)則①對于分解后圖像的低頻部分(位于最高分解層)采取平均算子;②對于高頻帶,采用基于區(qū)域矩形特性選擇加權(quán)平均算子;③對于3個方向的高頻帶,分別選用不同的特性選擇算子.對于高頻分量,基于區(qū)域特性選擇的融合算子的確定方法如下:先分別計算CT與MRI圖像對應(yīng)局部區(qū)域的能量Ej,A,Ej,B.即Ej(n,m)=∑n′∈L,m′∈KWε(n′,m′)[Dεj(n+n′,m+m′)]2,ε=1,2,3.其中Ej(n,m)表示2j分辨率下,以(n,m)為中心位置的局部區(qū)域能量;Dεj表示2-j分辨率下3個方向的高頻分量,Wε(n′,m′)為與Dεj對應(yīng)的權(quán)系數(shù);L,K定義了局部區(qū)域的大小(3×3),n′,m′的變化范圍在L,K內(nèi),最后確定融合算子.Dεj,F=Wεj,maxDεj,A+Wεj,minDεj,B,Ej,A≥Ej,B,Dεj,F=Wεj,minDεj,A+Wεj,maxDεj,B,Ej,A<Ej,B,其中Wεj,min+Wεj,max=1.不同融合方法的結(jié)果,可用目視判別:優(yōu)點是直接、簡單,可直接根據(jù)圖像處理前后的對比做出定性評價,缺點是主觀性較強.為了進一步客觀定量評價融合效果,從融合圖像包含的信息進行分析.對于不同的融合方法的結(jié)果,可以用熵、交叉熵、平均梯度、標(biāo)準(zhǔn)偏差等來定量描述.熵和標(biāo)準(zhǔn)偏差越大,說明融合后的圖像上的信息量增加得越多,達到富集表示的目的.圖像交叉熵直接反映了兩幅圖像對應(yīng)像素的差異,可用來確定各種融合效果的優(yōu)劣.交叉熵值越小,說明融合結(jié)果圖像與兩幅原始圖像的差異越小.平均梯度是敏感反映圖像對微小細節(jié)反差和紋理變化特征表達的能力,同時也反映了圖像的清晰度.一般平均梯度越大,圖像層次越多,融合后圖像越清晰,融合達到了提高空間分辨率的目的.

3融合實驗

實驗所采用的切片來自于美國國家醫(yī)學(xué)圖書館使用的可視化人體的橫截面解剖圖像(Cross-sectionalAnatomy)部分CT和MRI數(shù)據(jù)集.在圖像融合實驗中用的是已經(jīng)配準(zhǔn)的圖像。顱骨MRI醫(yī)學(xué)圖像,軟組織卻清晰.用提出的基于小波變換的區(qū)域能量融合規(guī)則,局部平均梯度融合規(guī)則,最大值融合規(guī)則算法分別實驗.圖1(~(分別是CT,MRI使用上述3種融合規(guī)則,分解2層和分解3層的融合結(jié)果.從圖1(和(,(和(,(和(可知,分解2層的融合圖像不如分解3層的融合圖像對微小細節(jié)反差和紋理變化特征表達的敏感能力高,這在融合規(guī)則的融合結(jié)果圖像統(tǒng)計參數(shù)對比表中反映為平均梯度增大,同時融合圖像的熵,標(biāo)準(zhǔn)偏差也增大.但是由于同時分解層數(shù)增多,交叉熵也增大,其融合目視效果不如分解2層的結(jié)果.因為圖像交叉熵直接反映了兩幅圖像對應(yīng)像素的差異,是確定各種融合效果優(yōu)劣的一個重要標(biāo)準(zhǔn).綜合圖1(~(,分解2層的基于小波變換的區(qū)域能量融合規(guī)則的融合圖像圖1(的融合效果最好.從均值、方差來看,此融合結(jié)果與原圖像(CT/MRI)更為接近.在分解2層的融合結(jié)果中,從熵值增大,交叉熵減小來看,融合結(jié)果攜帶的信息量較大,損失小達到富集的目的,又與兩幅原圖像的差異較小.從平均梯度增大來看,融合后的圖像清晰度提高.在融合圖像中,可看出CT提供的骨信息和MR提供的關(guān)于軟組織的信息都非常清晰.CT,MRI醫(yī)學(xué)圖像采用小波分析方法融合處理后,可有效地保留原始多源圖像的邊緣和紋理特征,避免了融合圖像因平均化出現(xiàn)模糊現(xiàn)象.充分利用了小波變換的時頻性,對于CT和MRT的圖像的不同頻率,不同空間細節(jié)做不同處理從而使各自的原始圖像的高低頻分量同時做了綜合提取,保存了邊緣信息避免了融合后的2D圖像的模糊,而且CT提供的骨信息和MR提供的關(guān)于軟組織的信息都能在最后的融合圖像中有機的結(jié)合在一起.這在放射治療計劃中,很好地利用了能確定目標(biāo)療區(qū)輪廓的MR圖像和能提供計算劑量分布的放射衰減系數(shù)的CT圖像,二者的基于小波變換的有效融合能夠提高放射治療計劃的準(zhǔn)確性.此外,采用的雙正交小波基,在小波分解時有效的提取圖像邊緣特征,重建圖像時可減少圖像方塊效應(yīng).

4實驗結(jié)論

為了融合信息的良好表達,圖像融合不是圖像的簡單加減運算,它是一個復(fù)雜的處理過程,如何在一幅圖像上提供更多的信息,不僅提供骨組織、病變的信息,還要盡可能地把其他組織的信息清楚地表達出來,這樣信息結(jié)構(gòu)、層次越完整,融合圖像的效果越好.同時融合圖像統(tǒng)計參數(shù)對比表明任何一種算法有其局限性,在具體的融合過程中應(yīng)根據(jù)處理對象的特點,通過反復(fù)實驗選擇最佳方法,實現(xiàn)多源圖像的優(yōu)勢互補.小波分析方法在醫(yī)學(xué)圖像融合中的應(yīng)用廣泛.它不僅可用于CT圖像上觀察到的骨組織結(jié)構(gòu)和MR圖像上對照軟組織信息的融合,而且還用于來源于CT或MR圖像的解剖信息與來源于PET或SPECT圖像的功能信息融合.在外科手術(shù)導(dǎo)航系統(tǒng)中,將手術(shù)前所得的CT和MR的病灶三維圖像與手術(shù)中所得到的實時X熒光圖像或超聲圖像進行融合[8,9],有利于實時地指導(dǎo)和觀察,確保手術(shù)順利準(zhǔn)確地進行.