当前位置:首页 >> 信息与通信 >>

CT图像投影滤波重建










医学图像处理实验报告

CT 图像投影域加噪与滤波去噪

中南大学生物医学工程系 生物医学工程 1101 班 指导老师:喻罡 文勋喆、苏猛、龚书滔、王正果

2014 年 4 月 27 日

1.

实验目的
A.熟悉 CT 图像重建原理; B.对比几种重建图像方法,掌握反投影重建法; C.在投影域加噪声,熟悉去噪算法。

2. 实验原理
根据著名的 Radon 变换,能到体膜在不同角度的

3. 实验平台
Windows 7 操作系统 Matlab 2013b 软件

4. 实验步骤
a. b. c. d. e. 构造体膜 Radon 变换,加上高斯噪声 利用算法进行去噪声 iRadon 变换得到去噪后的图像 得到结果,提出结论

5. 实验结果

A. 体膜图像
I=phantom(256);

B. Radon 变换
theta=0:179; [R]=radon(I,theta)

C. 图像加噪处理
MatLab 加噪函数: 高斯噪声: J1=10*imnoise(1.5e-2*R,'gaussian',0,0.003);

椒盐噪声:

J2=10*imnoise(1.5e-2*R,'salt & pepper',0.03);

随机噪声:

Jr=5*randn(size(R));

D. 投影域去噪滤波
MatLab滤波函数:
高斯滤波器: Lg=fspecial('gaussian'); L1=imfilter(J,Lg)

自适应维纳滤波器:

L2=wiener2(J,[3 3]);

中值滤波器:

L3=medfilt2(J);

均值平滑滤波:

LL=fspecial('average'); L4imfilter(J,LL);

E. 图像反投影重建

Reimage=iradon(L,[0:179],'linear','Shepp-Logan');

6. 结论分析
A. 在投影域里的加入各种噪声,如高斯噪声,椒盐噪声等,它们都比较大程度的影响 了重建图像的质量,因此需要除噪声。 B. 除噪声既能在重建图上进行,也能在投影图上进行。此次试验在投影图上进行,及 时去除噪声,以免噪声信号进入重建环节而使图像质量变差。这样使图像的除噪性 能更好。 C. 中值滤波和均值滤波能有效的去除噪声,但当随着噪声方差加大,噪声变大,效果 会变得不理想。 D. 加噪声函数 imnoise 会影响图像重建效果,如实验结果中可以看出,去除用 imnoise 函数写的噪声图像经过去噪再重建的效果不好,而用 randn 函数加上的随机高斯噪 声经过滤波在反投影重建图像时,效果明显好于前者。

附录
程序源代码:
I=phantom(256); theta=0:179; [R]=radon(I,theta); J1=10*imnoise(1.5e-2*R,'gaussian',0,0.003); J2=10*imnoise(1.5e-2*R,'salt & pepper',0.03); Jr=5*randn(size(R)); J3=R+Jr figure,imshow(R,[]) figure,imshow(J1,[]) figure,imshow(J2,[]) figure,imshow(J3,[]) Lg=fspecial('gaussian'); L1_1=imfilter(J1,Lg); L1_2=imfilter(J2,Lg); L1_3=imfilter(J3,Lg); figure,imshow(L1_1,[]),xlabel('???????¨?÷?????????ù'); figure,imshow(L1_2,[]),xlabel('???????¨?÷???·?????ù'); figure,imshow(L1_3,[]),xlabel('???????¨?÷?????ú???ù'); L2_1=wiener2(J1,[3 3]); L2_2=wiener2(J2,[3 3]); L2_3=wiener2(J3,[3 3]); figure,imshow(L2_1,[]),xlabel('???????¨?÷?????????ù'); figure,imshow(L2_2,[]),xlabel('???????¨?÷???·?????ù'); figure,imshow(L2_3,[]),xlabel('???????¨?÷?????ú???ù'); L3_1=medfilt2(J1); L3_2=medfilt2(J2); L3_3=medfilt2(J3); figure,imshow(L3_1,[]),xlabel('???????¨?÷?????????ù'); figure,imshow(L3_2,[]),xlabel('???????¨?÷???·?????ù'); figure,imshow(L3_3,[]),xlabel('???????¨?÷?????ú???ù'); LL=fspecial('average'); L4_1=imfilter(J1,LL); L4_2=imfilter(J2,LL); L4_3=imfilter(J3,LL); figure,imshow(L4_1,[]),xlabel('?ù?????????¨?????????ù'); figure,imshow(L4_2,[]),xlabel('?ù?????????¨???·?????ù');

figure,imshow(L4_3,[]),xlabel('?ù?????????¨?÷?????ú???ù'); Reimage_1=iradon(L1_3,[0:179],'linear','Shepp-Logan'); Reimage_2=iradon(L2_3,[0:179],'linear','Shepp-Logan'); Reimage_3=iradon(L3_3,[0:179],'linear','Shepp-Logan'); Reimage_4=iradon(L4_3,[0:179],'linear','Shepp-Logan'); figure,imshow(Reimage_1),xlabel('?????????ú???ù·????°???¨'); figure,imshow(Reimage_2),xlabel('?????????ú???ù·????°???¨'); figure,imshow(Reimage_3),xlabel('?????????ú???ù·????°???¨'); figure,imshow(Reimage_4),xlabel('?ù???????????ú???ù·????°???¨');


相关文章:
CT反投影滤波重建算法设计
地理与生物信息学院 2012 / 2013 学年 第二 学期 实验报告 课程名称: 医学图像处理和成像技术 实验名称:CT投影滤波重建算法设计 班学指 级生导 学姓教 号...
CT成像模拟
X-CT 三维图像重建的基本算法 X-CT 三维图像重建的基本算法包括:矩阵法、迭代法、傅里叶变换法、反投影法,其 中反投影法包括滤波投影法和卷积反投影法。 ...
CT仿真实验报告
二、实验要求用 CTSim 程序完成滤波投影重建实验,分析评价结果。 三、实验内容 1、对人脑体模仿真数据与实际的 CT 图像数据进行滤波投影重 建实验, 比较直接...
第六讲-CT成像技术(3)
第六讲-CT 成像技术(3) X-CT 反投影重建成像原理 X-CT 图像的质量评价一...? 现在应用比较多的是投影法、二维傅里叶重建法和滤波投影法(filtering back...
CT图像重建采用的基本方法是 A.滤波反投影 B.傅立叶变...
CT图像重建采用的基本方法是 A.滤波投影 B.傅立叶变换 C.迭代重建D.重复反投影 E.二维傅立叶变换 正确答案及相关解析 正确答案 A 解析 暂无解析 ...
CT成像的原理和方法
一般情况下, 我们需要得到 360 个全方位角度的投影数据。 3、 CT 图像重建原理: CT 重建算法有很多:直接反投影重建算法、滤波投影重建算法、卷积反 投影重建...
CT实验平行束投影与反投影重建报告
投影,滤波重建程序 (2 周) 3.将二维平行束 CT 程序融于数字图像平台软件(1 周) 4.学习二维扇束计算机层析成像原理(2 周) 5.编制二维扇束投影,滤波与...
CT成像技术模拟试题(16)
CT 成像技术模拟试题 成像技术模拟试题(16) 1、CT 图像重建采用的基本方法是 A. 滤波投影 B. 傅立叶变换 C. 迭代重建 D. 重复反投影 E. 二维傅立叶...
压缩感知在有限投影角度CT重建中的应用 (1)
1.3 有限角度 CT 重建的发展现状 CT 图像重建算法是 CT 技术中较为关键的部分,常用的 CT 重建算法有这两 种: 滤波投影重建算法和迭代重建算法。在投影数据...
CT实验一
CT实验一_临床医学_医药卫生_专业资料。实验课程:现代医学影像设备 生物医学工程...图2-4滤波投影的结果 实验结果: 由上图可知,反投影重建滤波投影重建都...
更多相关标签:
滤波反投影重建算法 | 二维图像灰度投影滤波 | 由投影重建图像 | 图像投影重建 | ct图像重建 | ct图像重建算法 | ct图像三维重建 | matlab ct 图像重建 |