1 簡(jiǎn)介

FDK 重建算法是一種近似的重建方法，即它是把所有不通過(guò)幾何中心平面的錐束投影數(shù)據(jù)近似地看做是幾何中心平面的扇束經(jīng)過(guò)傾斜一個(gè)角度得到的投影數(shù)據(jù)，然后對(duì)投影數(shù)據(jù)進(jìn)行修正，最后使用扇束濾波反投影算法進(jìn)行重建。實(shí)際上FDK 重建方法也是一種濾波反投影算法。

?

根據(jù)Feldkamp、Davis和Kress等學(xué)者所做的研究工作，可知實(shí)現(xiàn)FDK重建算法，必須滿足以下條件：探測(cè)器-射線源兩者固定不動(dòng)，被掃描物體繞垂直軸作自旋運(yùn)動(dòng)；被掃描物體必須完全處于 X 射線源的錐束范圍內(nèi)；X 光射線源必須為點(diǎn)光源；重建結(jié)果采用體素模型表示。 為了方便分析和實(shí)現(xiàn) FDK 重建算法，人們提出了虛擬探測(cè)器這個(gè)概念。虛擬探測(cè)器與實(shí)際的探測(cè)器平行，且位于載物臺(tái)的旋轉(zhuǎn)中心軸上。

2 部分代碼

clear allclose allhead =phantom(256);%獲得扇形束等距離投影數(shù)據(jù)D = 512; ?%射線源到圖象中心的距離[shade] = fanbeam(head,D,'FanSensorGeometry','line','FanSensorSpacing',1); ? ?%%%%獲得投影數(shù)據(jù) 線形檢測(cè)器 通道數(shù)sp ?旋轉(zhuǎn)角度360。figure;imshow(shade,[]);title('360度投影數(shù)據(jù)');clear Xangle;figuresubplot(121);imshow(re,[])title('重建后的圖象');subplot(122);imshow(head,[])title('原始圖象');p1 = re(128,:);p2 = head(128,:);figureplot(p1,'r');hold onplot(p2,'b')title('紅---重建圖象????藍(lán)---原始圖象')

3 仿真結(jié)果

4 參考文獻(xiàn)

[1]張斌, 潘晉孝. 三維圖像重建P—FDK算法的一種改進(jìn)方法[J]. 期刊論文, 2009.

博主簡(jiǎn)介：擅長(zhǎng)智能優(yōu)化算法、神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)、信號(hào)處理、元胞自動(dòng)機(jī)、圖像處理、路徑規(guī)劃、無(wú)人機(jī)等多種領(lǐng)域的Matlab仿真，相關(guān)matlab代碼問(wèn)題可私信交流。

部分理論引用網(wǎng)絡(luò)文獻(xiàn)，若有侵權(quán)聯(lián)系博主刪除。