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

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图像投影滤波重建
中 南 大 学 医学图像处理实验报告 CT 图像投影域加噪与滤波去噪 中南大学生物医学工程系 生物医学工程 1101 班 指导老师:喻罡 文勋喆、苏猛、龚书滔、王...
CT反投影滤波重建算法设计
地理与生物信息学院 2012 / 2013 学年 第二 学期 实验报告 课程名称: 医学图像处理和成像技术 实验名称:CT投影滤波重建算法设计 班学指 级生导 学姓教 号...
文档CT图像滤波反投影重建算法的研究
文档CT图像滤波投影重建算法的研究_数学_自然科学_专业资料。CT图像滤波投影重建算法的研究西北工业大学 学位论文知识产权声明书 本人完全了解学校有关保护知识产权...
用图像的投影数据进行重建程序
图像的投影数据进行重建程序_IT/计算机_专业资料。今日推荐 160份文档 四季...CT投影滤波重建算法设... 8页 2下载券©2014 Baidu 使用百度前必读 | 文库...
滤波反投影
滤波投影重建 CT 图像实验指导书 一、 实验目的 1. 了解傅立叶变换法、直接反投影重建 CT 图像的原理; 2. 掌握滤波投影重建 CT 图像的原理和基本...
基于投影的图像重建
(CT)为代表的各种层析成像方法对被测对象产生的并不是其 直接三维图像,而是从...因此可将这样一个重建的过程命名为滤波投影。 3 有限频带的作用 3.1 对一...
CT影像重建技术
CT影像重建技术_医药卫生_专业资料。CT是一种先进的成像技术,它通过扫描和图像重建...3 图像重建算法 图像重建算法主要分为两大类,一类是滤波投影算法,以傅里叶...
CT图像重建
年月日 学号 CT 图像重建 开课实验室: 3208 姓名 指导教师 2008 年 12 月...第二步:卷积(滤波)将第一步修改了的投影与响应函数 g (? ) 进行卷积 (2....
CT断层图像重建算法研究
CT 断层图像重建算法研究专业:通信工程 姓名:刘明帅 指导教师:骆岩红 摘要 CT ...重点对解析法中的滤波投影算法从 平行束与扇束投影方式进行研究,最后通过 Visua...
CT技术的图像重建
从反投影算法到滤波投影算法和迭代算法,从平行光束到扇形光束、锥形光束的重建 算法,图像重建问题在数学方法上的进展,为 CT 技术在各个领域成功的和不断拓广的...
更多相关标签:
滤波反投影重建算法 | 基于投影的图像重建 | ct图像重建 | ct图像三维重建 | matlab ct 图像重建 | ct图像重建算法 | ct图像重建算法 pdf | ct 图像重建 原理 |