地理与生物信息学院2012/2013学年第二学期实验报告课程名称：医学成像技术实验名称：磁共振图像后处理算法设计班级学号:B10090405学生姓名:陈洁指导教师:戴修斌日期：2013年5月实验题目：磁共振图像后处理算法设计二、实验内容：1．对图像进行去除噪声操作；2．对图像进行灰度变换操作；三、实验目的：1.加强下同学们实际的动手编程能力；2.重在体验和过程；实验过程：实验1：对图像进行去除噪声操作：1．操作步骤：对图像加入高斯噪声使用中值滤波对图像进行去噪处理模板尺寸设为5×5，也可自己设定图像边缘缺失部分使用对称方法补足算法实现流程：读入图像函数：imread()，中值滤波函数：medfilt2()；实验2：对图像进行灰度变换操作1．操作步骤：原图像灰度范围[50150]内的像素灰度值转成[10250]范围；原图像灰度范围[50150]内的像素灰度值转成[20200]范围；2.算法实现流程：源代码：clear;clc;iptsetpref('ImshowBorder','tight');I=imread('C:\DocumentsandSettings\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('中值滤波');a=50;b=150;m=10;n=250;c=20;d=200;g=zeros(512,512);f=zeros(512,512);forx=1:512fory=1:512ifI(x,y)<ag(x,y)=m;elseifI(x,y)>bg(x,y)=n;elseg(x,y)=(n-m)/(b-a)*(I(x,y)-a)+m;endendendendforx=1:512fory=1:512ifI(x,y)<af(x,y)=c;elseifI(x,y)>bf(x,y)=d;elsef(x,y)=(d-c)/(b-a)*(I(x,y)-a)+c;endendendendfigure,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有了更深层次的认识，比如说很多编程的技巧问题，程序的结构设计问题，对于程序的运行效率非常有帮助。有的时候，编出来的程序，能够运行，但是耗时太长，也就是说程序没有错，但是不适合实际运行。或者说，对于规模小的问题能够解决，但是规模大一点的问题就需要很长很长的时间，这就需要对程序的结构和算法问题进行改进。