《精》数字图像处理实验指导书2011.docVIP

《数字图像处理》 实 验 指 导 书 计算机科学与技术系 魏佳 2011年9月 实验一、图像直方图统计（验证性） 一、实验目的 1掌握MATLAB处理数字图像的基本方法； 2掌握MATLAB绘制图像直方图的方法。 二、实验原理 1 MATLAB窗口构成:在缺省的情况下，由三个窗口组成。命令窗口(command window)、命令历史（command history)、工作空间（workbench)。 注意：缺省窗口的设置为：MATLAB菜单/view选项/Desktop layout/default 调入一个文件：i=imread(.bmp)（观察此时的i出现在什么窗口？是什么类型的变量？大小是多少？显示这幅图：imshow(i);imhist(i)； 4 将变量i转置成j,即j=i;显示j即imshow(j); 写入到一个新的图像文件abc.tif中即imwrite(j,abc.tif)清除变量命令：clear清除用户开设的窗口命令：close all‘lena.bmp’); %读入原图像文件 imshow(I); %显示原图像 a=0; b=255; c=50; d=205; m=uint8(zeros(255,255)); for j=1:size(i,1) for k=1:size(i,2) if(i(j,k)=ai(j,k)b) m(j,k)=(d-c)/(b-a)*(i(j,k)-a)+c; else m(j,k)=i(j,k); end end end imshow(m); 三、实验步骤 1打开计算机，安装和启动MATLAB程序；程序组中“work”文件夹中应有待处理的图像文件； 2利用MatLab工具箱中的函数编制灰度变换文件FG; 3 运行FG文件； 4记录和整理实验报告。 四、实验仪器 1计算机； 2 MATLAB程序； 3移动式存储器（U盘等）。 4记录用的笔、纸。 五、实验报告内容 1叙述实验过程； 2提交实验的原始图像和结果图像。 六、思考题 1如何将两幅图像显示在一个窗口内？ 实验三、图像去噪（验证性） 一、实验目的 1掌握灰度图像简单的去噪方法； 2熟练掌握MATLAB实现模板操作的方法； 3熟练掌握空域滤波中常用的平滑滤波器； 4利用MATLAB程序进行图像去噪。 二、实验原理 见教材P65,4.4去除噪声。 下面给出图像去噪的MATLAB程序： I=imread(lenan.bmp); % 读入原图像 subplot(1,3,1); imshow(I); title(原图); J=fspecial(average); J1=uint8(filter2(J,I)); subplot(1,3,2); imshow(J1); title(3*3滤波); K=fspecial(average,9); K1=uint8(filter2(K,I)); subplot(1,3,3); imshow(K1); title(9*9滤波); 处理后的图像去除了噪声，但是会变得模糊。 三、实验步骤 1打开计算机，启动MATLAB程序；程序组中“work”文件夹中应有待处理的图像文件； 2给加噪图片“lenan.bmp”去噪； 4记录和整理实验报告。 四、实验仪器 1计算机； 2 MATLAB程序； 3移动式存储器（U盘等）。 4记录用的笔、纸。 五、实验报告内容 1叙述实验过程； 2提交实验的原始图像和结果图像。 六、思考题 1 试一试用fspecial函数生成其他滤波算子，观察去噪效果RGB = imread(原图像名.tif); I = rgb2gray(RGB); J = dct2(I); imshow(log(abs(J)),[]), colormap(jet(64)), colorbar J(abs(J) 10) = 0; K = idct2(J); figure,imshow(I) figure,imshow(K,[0 255]) 2）利用离散余弦变换进行JPEG图像压缩 I=imread(‘原图像名.tif’); %读入原图像； I=im2double(I); %将原图像转为双精度数据类型； T=dctmtx(8); %产生二维DCT变换矩阵 B=blkproc(I,[8 8],’P1*x*P2’,T,T’); %计算二维DCT，矩阵T及其转置T’是DCT函数P1*x*P2的参数 Mask=[ 1 1 1 1 0 0 0 0 1 1 1 0 0 0 0 0 1 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0