CT迭代重建算法的加速方法研究.pdf

1、CT成像作為放射性斷層成像的主要手段,在醫(yī)學(xué)診斷、輔助治療和工業(yè)無損檢測等方面都有廣泛應(yīng)用。經(jīng)過幾十年的研究,傳統(tǒng)的濾波反投影方法已漸趨成熟,目前的CT技術(shù)對完整均勻的無噪聲投影數(shù)據(jù)能以較快的速度完成重建。抵抗噪聲的CT成像重建方法主要是迭代重建法,它的缺陷是重建速度很慢,遠(yuǎn)遠(yuǎn)低于應(yīng)用需求,對該類算法的加速非常值得探討。本論文研究兩種對CT迭代重建算法進(jìn)行加速的方法,并通過模擬實驗來分析它們的主觀效果和客觀誤差。 第一部分工作是

2、基于統(tǒng)計迭代投影估計的快速CT重建算法?？紤]到傳統(tǒng)的圖像空間迭代法中對龐大的系統(tǒng)矩陣的計算和處理都極為耗時,嘗試在投影空間對噪聲投影進(jìn)行統(tǒng)計估計和迭代求解,獲得所有方向上的理想投影數(shù)據(jù),并結(jié)合具有速度優(yōu)勢的解析法來完最后的重建。首要的問題是對噪聲建立統(tǒng)計模型,采用了針對校準(zhǔn)后的數(shù)據(jù)建立的不同于傳統(tǒng)泊松描述的高斯模型,在該模型下,用貝葉斯方法推導(dǎo)了最大后驗概率(MAP)估計的能量函數(shù),同時還推導(dǎo)了懲罰加權(quán)最小均方(PWLS)估計的能量函數(shù)

3、,在采用相同的二次懲罰項的情況下,得到了類似的結(jié)果。 考慮到本文所處理的估計問題中,初始值可信度高,又是唯一可能的初始點,采用了迭代條件模式(ICM)優(yōu)化算法對MAP目標(biāo)函數(shù)進(jìn)行迭代求解。完成投影估計步驟之后,再進(jìn)行濾波反投影步驟的實驗。實驗中產(chǎn)生兩組不同信噪比的隨機噪聲數(shù)據(jù),分別進(jìn)行50次迭代估計和50次圖像重建,并對每次的投影估計和圖像重建結(jié)果都進(jìn)行主觀比較和客觀誤差分析。同時分別選取一種常用的線性濾波器和非線性濾波器進(jìn)行對

4、照。 實驗的圖像結(jié)果以及總的誤差趨勢分析和比較表明,在選擇合理的系統(tǒng)模型與優(yōu)化算法的情況下,該方法能獲得理想的估計投影,信噪比得到有效提升,誤差很快地減小,進(jìn)而有效地改善了圖像重建質(zhì)量。由于不需要處理系統(tǒng)矩陣,并且選取了收斂相當(dāng)快的優(yōu)化算法,該統(tǒng)計方法的時間代價極小。從誤差圖中發(fā)現(xiàn),選擇較少的迭代次數(shù)可避免少數(shù)點的誤差累積。正則因子beta的選取對投影估計的迭代后期影響較大,正則因子較小時,估計結(jié)果較好,而對重建結(jié)果來說,由于選

5、擇了較少的迭代次數(shù),其影響并不顯著。 第二部分工作是基于傅立葉方法的快速重投影。重投影是一般迭代法中必不可少、反復(fù)進(jìn)行的步驟,也是迭代法的速度瓶頸所在,對這一過程實現(xiàn)加速很有意義。考慮到CT成像的基本原理與直接傅里葉重建法有關(guān),討論了近年來非均勻傅里葉變換(NUFFT)技術(shù)在直接傅里葉重建中的研究進(jìn)展,引入該類技術(shù)對迭代算法的重投影步驟進(jìn)行加速,并與標(biāo)準(zhǔn)的系統(tǒng)矩陣重投影法進(jìn)行比較。同時用傳統(tǒng)的雙線性插值法實現(xiàn)重投影作為傅立葉方法

6、的對比。在NUFFT重投影算子的實現(xiàn)中,利用Kaiser-Bessel卷積核優(yōu)化方面的新進(jìn)展進(jìn)行前向投影算子和反向投影算子的加速,與最近的一種從最小化最壞近似誤差(Min-Max)判據(jù)直接導(dǎo)出的卷積核的重投影實現(xiàn)進(jìn)行比較,給出了實驗結(jié)果和誤差對比。 由實驗結(jié)果可以看到,第一,與傳統(tǒng)的系統(tǒng)矩陣計算方式相比,基于NUFFT的傅里葉方法在前向投影和反投影算子的實現(xiàn)上有明顯的速度優(yōu)勢。第二,在采用優(yōu)化Kaiser-Bessel插值器以及