资源描述:
《医学影像处理 实验报告.doc》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库。
1、医学成像技术与图像处理实验报告学院(系):电子信息与电气工程学部专业:生物医学工程学生姓名:学号:指导教师:刘惠完成日期:2013.06.25实验一:数字图像读取及色彩、亮度、对比度变化一、实验目的了解数字图像的存储格式,并学会对图像的某些视觉特征作简单处理。二、实验要求1.从最常用的“.BMP”图像格式中读取图像数据;2.对数字图像的表示方式(如RGB、YUV)及各种表示方式之间的转换有初步了解;3.根据输入参数改变数字图像的色彩、亮度、对比度。三、实验步骤1.利用工具(如ACDSee、Pho
2、toShop)将Sample.1.jpg转换为Sample1.bmp;2.借助imread命令将图像内容读入内存数组;3.通过访问数字图像RGB三个通道的对应矩阵,改变数字图像的色彩;4.将数字图像的RGB表示转换为YUV表示;Y=0.30R+0.59G+0.11BU=0.70R-0.59G-0.11BV=-0.30R-0.59G+0.89B5.通过访问Y(亮度)通道,改变数字图像的亮度;6.通过Y(亮度)通道作灰度的线性变换,改变数字图像的对比度。四、实验程序及图像处理结果pic=imread
3、('E:实验1.jpg');pic=rgb2gray(pic);subplot(1,2,1);imshow(pic);title('GRAY');INFO=imfinfo('E:实验1.jpg');INFO.Width,INFO.Heightpic1=pic(1:2:end,1:2:end);subplot(1,2,2);imshow(pic1);title('RE');functionfA=imread('E:实验Sample1.bmp');%将原图像读入内存[r,c,d]=siz
4、e(A);red=A;red(:,:,1)=A(:,:,1);red(:,:,2)=zeros(r,c);red(:,:,3)=zeros(r,c);red=uint8(red);subplot(3,3,1);imshow(red);%通过访问R通道来观察图像title('Red');green=zeros(r,c);green(:,:,2)=A(:,:,2);green(:,:,1)=zeros(r,c);green(:,:,3)=zeros(r,c);green=uint8(green);s
5、ubplot(3,3,2);imshow(green);%通过访问G通道来观察图像title('Green');blue=zeros(r,c);blue(:,:,3)=A(:,:,3);blue(:,:,1)=zeros(r,c);blue(:,:,2)=zeros(r,c);blue=uint8(blue);subplot(3,3,3);imshow(blue);%通过访问B通道来观察图像title('Blue');Y=0.3*red+0.59*green+0.11*blue;%将图像的RGB
6、转换为YUVU=0.70*red-0.59*green-0.11*blue;V=-0.30*red-0.59*green+0.89*blue;subplot(3,3,4);imshow(Y);title('原图亮度');%原图像亮度显示Y=2*Y;subplot(3,3,5);imshow(Y);%增大图像亮度并显示title('增大亮度');Y=Y/4;subplot(3,3,6);imshow(Y);%减小图像亮度并显示title('减少亮度');GreyR(:,:,1)=100+100/2
7、55*A(:,:,1);GreyR(:,:,2)=100+100/255*A(:,:,2);GreyR(:,:,3)=100+100/255*A(:,:,3);subplot(3,3,7);imshow(GreyR);title('线性变换');实验二:数字图像的噪声去除一、实验目的:完善程序,对图像进行去噪操作。二、实验要求1.用均值滤波器去除图像中的噪声;2.用中值滤波器去除图像中的噪声;3.比较两种方法的处理结果三、实验步骤1.将图像‘Lena-1.bmp’和‘Lena-2.bmp’读入内
8、存数组;2.用均值滤波器去除图像中的噪声;3.用中值滤波器去除图像中的噪声;将两种处理方法的结果与原图比较;4.注意两种处理方法对边缘的影响。四、实验程序及图像结果clci=imread('E:实验34.bmp');subplot(1,3,1);imshow(i);i1=avefiltm(i,5);%进行均值滤波subplot(1,3,2),imshow(i1);title('均值滤波');i2=midfiltm(i,5);subplot(1,3,3),imshow(i2);title('中