磁共振图像后处理算法设计

地理与生物信息学院

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 有了更深层次的认识，比如说很多编程的技巧问题，程序的结构设计问题，对于程序的运行效率非常有帮助。有的时候，编出来的程序，能够运行，但是耗时太长，也就是说程序没有错，但是不适合实际运行。或者说，对于规模小的问题能够解决，但是规模大一点的问题就需要很长很长的时间，这就需要对程序的结构和算法问题进行改进。