(2016版)图像处理实验1.docVIP

实验一、空域图像增强技术 一、实验目的 1结合实例学习如何在视频显示程序中增加图像处理算法； 2理解和掌握图像的线性变换和直方图均衡化的原理和应用； 3了解平滑处理的算法和用途，学习使用均值滤波、中值滤波和拉普拉斯锐化进行图像增强处理的程序设计方法； 4 了解噪声模型及对图像添加噪声的基本方法。 二、实验原理 1 灰度线性变换就是将图像中所有点的灰度按照线性灰度变换函数进行变换。 2 直方图均衡化通过点运算将输入图像转换为在每一级上都有相等像素点数的输出图像。按照图像概率密度函数PDF的定义： 通过转换公式获得： 3 均值（中值）滤波是指在图像上，对待处理的像素给定一个模板，该模板包括了其周围的临近像素。将模板中的全体像素的均值（中值）来代替原来像素值的方法。 4 拉普拉斯算子如下： 拉普拉斯算子首先将自身与周围的8个像素相减，表示自身与周围像素的差异，再将这个差异加上自身作为新像素的灰度。 三、实验步骤 1 启动MATLAB程序，对图像文件分别进行灰度线性变换（参考教材50页，例4.1）、直方图均衡化（参考教材56页，例4.6）、均值滤波（参考教材60页，例4.9）、中值滤波（参考教材64页，例4.11）和梯度锐化操作（参考教材66页，例4.12）。添加噪声，重复上述过程观察处理结果。 2记录和整理实验报告 四、实验仪器 1计算机； 2 MATLAB程序； 3记录用的笔、纸。 五、实验报告内容 实验前将准备好的灰度图像复制到计算机上（为了方便调用，我选择复制图片在MATLAB的work文件下）。学习课本线性灰度变换、直方图均衡化、均值滤波、中值滤波和梯度锐化等知识内容，参考书上的相关MATLAB程序代码，编写相应的程序。程序执行完后，保存输出的结果和图片，写成实验报告。 （1）直接对原图像进行处理 1、线性灰度变换 源程序 %读入并显示原始图像 I=imread(123.jpg); subplot(121); imshow(I); I=double(I); [M,N]=size(I); %进行线性灰度变换 for i=1:M for j=1:N if I(i,j)=30 I(i,j)=I(i,j); elseif I(i,j)=150 I(i,j)=(200-30)/(150-30)*(I(i,j)-30)+30; else I(i,j)=(255-200)/(255-150)*(I(i,j)-150)+200; end end end %显示变换后的结果 subplot(122);imshow(uint8(I)); 2、直方图均衡化 源程序 %直方图均衡化 I=imread(456.jpg); subplot(231);imshow(I);title(原图); I0=rgb2gray(I); %灰度化 subplot(232);imshow(I0) title(灰度图); subplot(233);imhist(I0) title(直方图); I1=histeq(I0); %均衡化 subplot(234);imshow(I1) subplot(235);imhist(I1) 3、均值滤波 源程序 %均值滤波 I=imread(456.jpg); I0=rgb2gray(I); %灰度图 [M,N]=size(I0); II1=zeros(M,N); for i=1:16 II(:,:, i)=imnoise(I0,gaussian,0,0.01); II1=II1+double(II(:,:, i)); if or (or(i==1,i==4),or(i==8,i==16)); figure; imshow(uint8(II1/i)); end end 4、中值滤波 源程序 %中值滤波 I=imread(456.jpg); I0=rgb2gray(I); J=imnoise(I0,salt pepper,0.02); subplot(231),imshow(I0);title(原图像); subplot(232),imshow(J);title(添加椒盐噪声图像)