图像显示实验报告

1、医学图像处理评分实验报告 课程名称 医学图像处理 实验名称 图像显示 专业班级 姓 名 学 号 实验日期 实验地点 20152016学年度第 2 学期一、 实验目的MATLAB的图像处理工具箱提供了多种图像显示和颜色、光照技术。本实验将具体介绍MATLAB中用于图像显示、颜色转换的各个函数的具体用法。二、实验环境 1、硬件配置：处理器：Intel（R） Core(TM) i7-3770 CPU 3.40GHz 3.40GHz 安装内存（RAM）：4.00GB 系统类型：64位操作系统 2、软件环境：MATLAB R2013b软件三、实验内容1. 读取图像使用函数imread使用函数imread

2、可以将图像读入MATLAB环境，imread的语法为imread(filename)其中，filename是一个含有图像文件全名的字符串（包括任何可用的扩展名）。注意：此处的单引号用以界定filename字符串。如命令行：f=imread(d:imagerice.tif);将TIFF图像rice.tif读入图像数组f。命令行结尾处的分号用于取消输出。若命令行未包含分号，则MATLAB会立即显示该行命令的运算结果。MATLAB桌面工具条上的当前目录窗口会显示MATLAB的当前工作路径。当filename中不包含任何路径信息时，imread会从当前目录中寻找并读取图像文件。函数size可以给出一幅

3、图像的行数和列数：size(f)ans= 256 256也可以使用如下形式自动确定一幅图像的大小：M,N=size(f);将返回图像的行数(M)和列数(N)。2. 显示图像使用函数imshow3. 保存图像使用函数imwrite4. 特殊图像显示技术添加颜色条显示多图像序列显示多帧图像动画显示z把真彩图像映射到柱面上5. 图像颜色转换技术（1） 图像的褪色 褪色函数imapprox，rgb2ind（2） 练习使用颜色模型转换函数：rgb2hsv,hsv2rgb,ntsc2rgb,rgb2ntsc,rgb2ycbcr,ycbcr2rgb,rgb2gray,im2bw四、实验结果与分析 （包括实验

4、原理、数据的准备、运行过程分析、源程序（代码）、图形图象界面等）注：本项可以增加页数例1J=imread('eight.tif');%将TIFF图像eight.tif读入图像数组Jsubplot(121);%把绘图窗口分成一行两列两块区域，在第一块区域绘图imshow(J,100,200); %在灰度范围100,200内显示图像subplot(122);%在第二块区域绘图imshow(J,20);%以灰度等级20显示图图1原eight图 图2灰度图像 因为在机房无法以灰度等级20显示图，所以下面用增大灰度范围显示图像J=imread('eight.tif');

5、%读取图像subplot(121); %生成一行二列的子图，这是第一个子图imshow(J,100,200); %在灰度范围100,200内显示图像subplot(122); %生成一行二列的子图，这是第二个子图imshow(J,0，100); %在灰度范围0,100内显示图像图3原eight图 图4灰度图像分析：显示灰度图像例2BW=imread('eight.tif '); %读取图像imshow(BW) %显示二值图像figure(2) %用figure函数，保持第一幅图像同时显示第二幅图像imshow(BW) %显示颜色反转后的二值图像图5 二值图像 图6颜色反转后的二

6、值图像分析：显示二值图像例3I=imread('cell.tif'); %读取图像imshow(I); %显示图像H=1 2 1;0 0 0;-1 -2 -1; %生成x方向上的Sobel算子模板I2=filter2(H,I); %以H模板进行filter傅里叶变换imshow(I2,),colorbar %在I2图像旁边垂直方向，加上颜色条图7 cell.Tif 图8 锐化图和颜色条分析：添加颜色条；Sobel算子主要是滤波作用，加强边缘和细节，削弱图像中的低频分量，具有高通滤波特性。例4mri=uint8(zeros(128,128,1,27);%产生8位无符号的4维整数数

7、组，其维数128*128*1*27，用来存储27帧核磁共振图像for frame=1:27 mri(:,:,:,frame),map=imread('mri.tif',frame);%循环读出多帧图像中的每一幅图像，共有27帧endimshow(mri(:,:,:,3),map); %显示第3帧图像图9 第3帧核磁共振图像分析：显示多帧图像中的第3帧；map是调色板矩阵，用于存储该索引色图像相关的颜色映射表例5mri=uint8(zeros(128,128,1,27); %产生8位无符号的4维整数数组，其维数128*128*1*27，用来存储27帧核磁共振图像for frame

8、=1:27 %循环过程从1至27 mri(:,:,:,frame),map=imread('mri.tif',frame); %读入图像至fram中，每次循环都存储在mri中，并给出其色谱值endmontage(mri,map); %拼接索引图像图10 27帧核磁共振图像分析：显示多帧图像，montage可以同时显示多帧图像的所有帧例6mri=uint8(zeros(128,128,1,27); %产生8位无符号的4维整数数组，其维数128*128*1*27，用来存储27帧核磁共振图像for frame=1:27 mri(:,:,:,frame),map=imread('

9、;mri.tif',frame); %循环读出多帧图像中的每一幅图像，共有27帧endmov=immovie(mri,map); %构建动画文件movie(mov); %播放动画图11 动画显示的最后一帧核磁共振图像分析：按顺序动画显示27帧图像例7x,y,z=cylinder; %纹理映射方式观察，向量X，Y，Z轴定义圆柱面中心线，矩阵X，Y，Z轴确定圆柱面RGB=imread('LungRGB.jpg'); %RGB=imread('9.jpg') %从图像文件中读取图像数据;warp(x,y,z,RGB); %将真彩色图像映射到柱面上图12 真彩图

10、的映射分析：把真彩色图像映射到柱面上，三维图像卷起生成例8X,map=imread('canoe.tif'); %读取索引图像Y,newmap=imapprox(X,map,20); %利用彩色映射newmap返回一个数组Y，该数组最多有20种颜色；用20种颜色代替原来的256种颜色。X可以是uint8类，uint16类或double类；若n小于256，则Y是uint8类，若n大于256，则Y是double类。imshow(X,map) %显示索引图像，X为数据图像矩阵，map是调色板figure,imshow(Y,newmap) %用figure函数，保持第一幅图像同时显示第

11、二幅图像对比原图图13原canoe图 图14 褪色后的图像分析:图像的褪色，用20种颜色代替原来的256种颜色。 颜色模型转换函数rgb2hsv函数:rgbmap=imread('canoe.tif')cmap=rgb2hsv(rgbmap)rgb_image=imread('canoe.tif')hsv_image=rgb2hsv(rgb_image)或者%将RGB模型转换为HSV模型J=imread('LungRGB.jpg');%读入文件名为LungRGB.jpg的图像HSV=rgb2hsv(J);%真彩色图像转换为HSV色彩空间图像ims

12、how(J)%显示图像Jfigure(2)%建立图形窗口2imshow(HSV)%显示经转换后的HSV色彩空间图像 图15 RGB图像 图16 HSV图像rgb2ntsc函数：%将RGB模型转换为NTSC模型J=imread('LungRGB.jpg');%读入文件名为LungRGB.jpg的图像YIQ=rgb2ntsc(J);%真彩色图像转换为YIQ色彩空间图像imshow(J)%显示图像Jfigure(2)%建立图形窗口2imshow(YIQ)%显示经转换后的YIQ色彩空间图像rgb2ycbcr函数：%将RGB模型转换为YCbCr模型J=imread('LungRG

13、B.jpg');%读入文件名为LungRGB.jpg的图像Ycbcr=rgb2ycbcr(J);%真彩色图像转换为YCbCr空间图像imshow(J)%显示图像Jfigure(2)%建立图形窗口2imshow(Ycbcr)%显示经转换后的YCbCr空间图像hsv2rgb函数：%将HSV模型转换为RGB模型J=imread('LungRGB.jpg');%读入文件名为LungRGB.jpg的图像HSV=rgb2hsv(J);%真彩色图像转换为HSV色彩空间图像RGB=hsv2rgb(HSV);%HSV色彩空间图像转换为真彩色图像imshow(HSV)%显示HSV色彩空间图

14、像figure(2)%建立图形窗口2imshow(RGB)%显示HSV色彩空间图像ntsc2rgb函数：%将NTSC模型转换为RGB模型J=imread('LungRGB.jpg');%读入文件名为LungRGB.jpg的图像NTSC=rgb2ntsc(J);%真彩色图像转换为YIQ色彩空间图像RGB=ntsc2rgb(NTSC);%YIQ色彩空间图像转换为真彩色图像imshow(NTSC)%显示YIQ色彩空间图像figure(2)%建立图形窗口2imshow(RGB)%显示真彩色图像ycbcr2rgb函数：%将YCbCr模型转换为RGB模型J=imread('Lung

15、RGB.jpg');%读入文件名为LungRGB.jpg的图像YCbCr=rgb2ycbcr(J);%真彩色图像转换为YCbCr空间图像RGB=ycbcr2rgb(YCbCr);%YCbCr空间图像转换为真彩色图像imshow(YCbCr)%显示YCbCr空间图像figure(2)%建立图形窗口2imshow(RGB)%显示真彩色图像五、实验小结： （包括主要实验问题的最终结果描述、详细的收获体会，待解决的问题等）在此次实验中，由于实验内容更贴近生活应用，因此比起上学期，我们更容易领悟该程序的表达，只是在细节方面还是很容易出错，甚至不容易拐过弯来。但经过此次实验，我们懂得要学着从复杂的程序中剥茧抽丝，把程序尽可能的简单化。 在实验中应注意的点：在编写运行文件的时候要选择正确的路径，图片保存在