地理与生物信息学院
2012/ 2013 学年第二学期
实验报告
课程名称：医学成像技术
实验名称：磁共振图像后处理算法设计
班级学号: B10090405
学生姓名: 陈洁
指导教师: 戴修斌
日期：2013 年 5 月
一、实验题目：磁共振图像后处理算法设计 二、实验内容：
1．对图像进行去除噪声操作 ；
2．对图像进行灰度变换操作 ；
三、实验目的：
1.加强下同学们实际的动手编程能力 ；
2.重在体验和过程 ；
四、 实验过程：
实验1：对图像进行去除噪声操作： 1．操作步骤：
1) 对图像加入高斯噪声
2) 使用中值滤波对图像进行去噪处理 3) 模板尺寸设为5×5，也可自己设定 4) 图像边缘缺失部分使用对称方法补足
51141671
81
91
71819151141611
21
31
1121311121511471
81
71
51113121161481
311691
91
1471
81
51718171
51711481
91
1691811691
91
2． 算法实现流程：
1) 读入图像函数：imread()，中值滤波函数：medfilt2()；
实验2：对图像进行灰度变换操作
1．操作步骤：
1) 原图像灰度范围[50 150]内的像素灰度值转成[10 250]范围；
2) 原图像灰度范围[50 150]内的像素灰度值转成[20 200]范围；
2.算法实现流程：
源代码：
clear;clc;
iptsetpref('ImshowBorder','tight');
I = imread('C:\Documents and Settings\nupt\桌面\4.bmp'); J = imnoise(I,'gaussian',0.02,0.02); K = medfilt2(J,[5,5]); figure,imshow(I),title('原图'); figure,imshow(J),title('高斯噪声'); figure,imshow(K),title('中值滤波');
f (x , y )
a
m
b
n g (x , y )
⎪⎪
⎩⎪
⎪⎨⎧>≤≤+---<=b y x f n b y x f a m a y x f a b m
n a y x f m y x g ),(
),( ]),([),(
),(
a=50;
b=150;
m=10;
n=250;
c=20;
d=200;
g=zeros(512,512);
f=zeros(512,512);
for x=1:512
for y=1:512
if I(x,y) g(x,y)=m; else if I(x,y)>b g(x,y)=n; else g(x,y)=(n-m)/(b-a)*(I(x,y)-a)+m; end end end end for x=1:512 for y=1:512 if I(x,y) f(x,y)=c; else if I(x,y)>b f(x,y)=d; else f(x,y)=(d-c)/(b-a)*(I(x,y)-a)+c; end end end end figure,imshow(g,[]),xlabel('原灰度图像50,150--10,250'); figure,imshow(g,[]),xlabel('原灰度图像50,150--20,200'); 3.运行结果： 原图：高斯噪声： 中值滤波'：原灰度图像50,150--10,250： 原灰度图像50,150--20,200： 五．实验小结： 通过这次实验，我了解了磁共振图像后处理算法设计的内容，同时也对matlab 有了更深层次的认识，比如说很多编程的技巧问题，程序的结构设计问题，对于程序的运行效率非常有帮助。有的时候，编出来的程序，能够运行，但是耗时太长，也就是说程序没有错，但是不适合实际运行。或者说，对于规模小的问题能够解决，但是规模大一点的问题就需要很长很长的时间，这就需要对程序的结构和算法问题进行改进。
