CT重建实验报告.docx

平行光照射的CT重建
艾文书220171807
一、实验目的
认识并熟悉CT重建的原理和成像基础。通过编写CT图像重建程序，进一步了解CT重建过程，同时加强图像处理程序的编程训练。
二、实验软件：MATLAB
三、实验要求
.递
15
16
17
13
19
20
21
22
23
342
26
27
28
2930丸
五、MATLAB程序算法及结果
第辑器发布视圄
lradon+mrIprojectm:+
■
—clc»clear：
，1adonS城前的图像%
—P=phantoin('UodLfiedsh*pp-LcganT256)；
-figure,imshow(?):
-titled原图,)；
%Radcj^换■
-theta=1inspace(0s179r360)：
R=radQn(PT;
为利用imd。口变抵用36。个角度来重构国保
Il=irado3i(RTthfeta)：
figurefinshow(ll)360个角度重构图憬/);
用1radon%
-p=phantoii(^Modifiedshepp-logari,R256):
-A-linspace(0,L79,130)：
—xp.-radontp,A);
%方法L频域避波法一一一一一
—*leruAN]=sise(Kj:
—JenFilt-2"nextpov2(2*size(Rj1.)),4疆;皮器的长度
n=0;(LenFiLt/2);
v=2*pi*n/lenFiit；
filt-Xend:-1:2)^;t(1：1：end-l)?]:
,堂接到颜域、源液、傅里叶反笠摭
sl-aeros(siae(R>l)9size(K?
fori=l:sisetK2)
t€»p=^(:,i);
3?
33
34瑞
33
33
40
41
42君
44石的
47
43
49口。=1
3Z
S3
54
—一8
5S
57
53
□3
temp'.lenFilt|=0:
t^mpF=fftshift(fft(t&tp))：
terapF=teiiipF,*fiIt:小滤波
tenp*Teal(ifft(fftshift(teapF))):sl(;,i)=temp(l;len);
Lemi

h=ssras(I>lenFilt/3>l):
n=Q：LenJilt/2;
h(l)=1/4:
h(2:2:end)=-I./((n(2:ZendO.*2)*(pi'Z));
li-［匕'；h(end—1：-1：2).11；
"Ifori=L:AN
S=;
S(lanFdlt)=0：
S=ccwvCS,h,lenFilt);
醴2,0-SU；1«1)：
Lend
\反报殷
N-2*flooi(si_2e(XD〃2返中七圾)))；%空间分辨窣拳数
■iproject(Kl»I);
=ipjQjcctCs2,N);
士显不
figure*imshow(I2)：titleC频博魂装反投题图像,〕：
ficuEB,imshowt13);tithC时域卷"反投影图最'):
1L
1