《图像处理》实验

实验一、图像的傅立叶变换 、实验目的 1 了解图像变换的意义和手段； 2 熟悉傅里叶变换的孩本性质； 3 热练掌握 FFT 酌方法反应用； 4 通过实验了解二维频谱的分布特点； 5 通过本实验掌握利用 MATLAB 编程实现数字图像的傅立叶变换。 、实验原理 1 应用傅立叶变换进行图像处理 傅里叶变换是线性系统分析的一个有力工具，它能够定量地分析诸如数字化系统、采 样点、电子放大器、卷积滤波器、噪音和显示点等的作用。通过实验培养这项技能，将有 助于解决大多数图像处理问题。对任何想在工作中有效应用数字图像处理技术的人来说， 把时间用在学习和掌握博里叶变换上是很有必要的。 2 傅立叶（Fourier）变换的定义 对于二维信号，二维 Fourier 变换定义为： 二维离散傅立叶变换为： 图像的傅立叶变换与一维信号的傅立叶变换变换一样，有快速算法，具体参见参考书 目，有关傅立叶变换的快速算法的程序不难找到。实际上，现在有实现傅立叶变换的芯片， 可以实时实现傅立叶变换。 3 利用 MATLAB 软件实现数字图像傅立叶变换的程序： I=imread(原图像名.gif); %读入原图像文件 imshow(I); %显示原图像 fftI=fft2(I); %二维离散傅立叶变换 sfftI=fftshift(fftI); %直流分量移到频谱中心 RR=real(sfftI); %取傅立叶变换的实部 II=imag(sfftI); %取傅立叶变换的虚部 A=sqrt(RR.2+II.2);%计算频谱幅值 A=（A-min(min(A)）/(max(max(A)-min(min(A)*225; %归一化 figure; %设定窗口 imshow(A); %显示原图像的频谱 、实验步骤 1 打开计算机，安装和启动 MATLAB 程序；程序组中“work ”文件夹中应有待处理的 图像文件； 2 利用 MatLab 工具箱中的函数编制 FFT 频谱显示的函数; 3 a).调入、显示“实验一”获得的图像；图像存储格式应为 “.gif”; b)对这三幅图像做 FFT 并利用自编的函数显示其频谱; c)讨论不同的图像内容与 FFT 频谱之间的对应关系。 4 记录和整理实验报告。 、实验仪器 1 计算机； 2 MATLAB 程序； 3 移动式存储器（软盘、U 盘等） 。 4 记录用的笔、纸。 、实验报告内容 1 叙述实验过程； 2 提交实验的原始图像和结果图像。 、思考题 1傅里叶变换有哪些重要的性质? 2图像的二维频谱在显示和处理时应注意什么? 实验二、图像增强 、实验目的 1 掌握灰度直方图的概念及其计算方法； 2 熟练掌握直力图均衡化和直方图规定化的计算过程； 3 熟练掌握空域滤波中常用的平滑和锐化滤波器； 4 掌握色彩直方图的概念和计算方法 5 利用 MATLAB 程序进行图像增强。 、实验原理 图像增强是指按特定的需要突出一幅图像中的某些信息，同时，消弱或去除某些不需 要的信息的处理方法。其主要目的是处理后的图像对某些特定的应用比原来的图像更加有 效。图像增强技术主要有直方图修改处理、图像平滑化处理、图像尖锐化处理和彩色处理 技术等。本实验以直方图均衡化增强图像对比度的方法为主要内容，其他方法同学们可以 在课后自行联系。 直方图是多种空间城处理技术的基础。直方图操作能有效地用于图像增强。除了提供 有用的图像统计资料外，直方图固有的信息在其他图像处理应用中也是非常有用的，如图 像压缩与分割。直方图在软件中易于计算，也适用于商用硬件设备，因此，它们成为了实 时图像处理的一个流行工具。 直方图是图像的最基本的统计特征，它反映的是图像的灰度值的分布情况。直方图均 衡化的目的是使图像在整个灰度值动态变化范围内的分布均匀化，改善图像的亮度分布状 态，增强图像的视觉效果。灰度直方图是图像预处理中涉及最广泛的基本概念之一。 图像的直方图事实上就是图像的亮度分布的概率密度函数，是一幅图像的所有象素集 合的最基本的统计规律。直方图反映了图像的明暗分布规律，可以通过图像变换进行直方 图调整，获得较好的视觉效果。 直方图均衡化是通过灰度变换将一幅图像转换为另一幅具有均衡直方图，即在每个灰 度级上都具有相同的象素点数的过程。 下面给出直方图均衡化增强图像对比度的 MATLAB 程序： I=imread(原图像名.gif); % 读入原图像 J=histeq(I); %对原图像进行直方图均衡化处理 Imshow(I); %显示原图像 Title(原图像); %给原图像加标题名 Figure；imshow(J); %对原图像进行屏幕控制；显示直方图均衡化后的图像 Title(直方图均衡化后的图像) ; %给直方图均衡化后的图像加标题名 Figure; subplot(1,2,1) ;%对直方图均衡化后的图像进行屏幕控制；作一幅子图作为并排 两幅图的第 1 幅图 Imhist(I,64); %将原图像直方图显示为 64 级灰度 Title(原图像直方图) ; %给原图像直方图加标题名 Subplot(1,2,2); %作第 2 幅子图 Imhist(J,64) ; %将均衡化后图像的直方图显示为 64 级灰度 Title(均衡变换后的直方图) ; %给均衡化后图像直方图加标题名 处理后的图像直方图分布更均匀了，图像在每个灰度级上都有像素点。从处理前后的 图像可以看出，许多在原始图像中看不清楚的细节在直方图均衡化处理后所得到的图像中 都变得十分清晰。 、实验步骤 1 打开计算机，启动 MATLAB 程序；程序组中“work ”文件夹中应有待处理的图像文 件； 2 调入“实验一”中获取的数字图像，并进行计算机均衡化处理； 3 显示原图像的直方图和经过均衡化处理过的图像直方图。 4 记录和整理实验报告 、实验仪器 1 计算机； 2 MATLAB 程序； 3 移动式存储器（软盘、U 盘等） 。 4 记录用的笔、纸。 、实验报告内容 1 叙述实验过程； 2 提交实验的原始图像和结果图像。 、思考题 1 直方图是什么概念？它反映了图像的什么信息？ 2 直方图均衡化是什么意思？它的主要用途是什么？ 实验三、图像压缩 一、 实验目的 1 理解有损压缩和无损压缩的概念； 2 理解图像压缩的主要原则和目的； 3 了解几种常用的图像压缩编码方式。 4 利用 MATLAB 程序进行图像压缩。 二、 实验原理 1.图像压缩原理 图像压缩主要目的是为了节省存储空间，增加传输速度。图像压缩的理想标准是信息 丢失最少，压缩比例最大。不损失图像质量的压缩称为无损压缩，无损压缩不可能达到很 高的压缩比；损失图像质量的压缩称为有损压缩，高的压缩比是以牺牲图像质量为代价的。 压缩的实现方法是对图像重新进行编码，希望用更少的数据表示图像。 信息的冗余量有许多种，如空间冗余，时间冗余，结构冗余，知识冗余，视觉冗余等， 数据压缩实质上是减少这些冗余量。高效编码的主要方法是尽可能去除图像中的冗余成分， 从而以最小的码元包含最大的图像信息。 编码压缩方法有许多种，从不同的角度出发有不同的分类方法，从信息论角度出发可 分为两大类。 （1）.冗余度压缩方法，也称无损压缩、信息保持编码或嫡编码。具体说就是解码图 像和压缩编码前的图像严格相同，没有失真，从数学上讲是一种可逆运算。 （2）信息量压缩方法，也称有损压缩、失真度编码或烟压缩编码。也就是说解码图像 和原始图像是有差别的，允许有一定的失真。 应用在多媒体中的图像压缩编码方法，从压缩编码算法原理上可以分为以下 3 类： （1）无损压缩编码种类 哈夫曼（Huffman）编码，算术编码，行程（ RLE）编码，Lempel zev 编码。 （2）有损压缩编码种类 预测编码，DPCM，运动补偿； 频率域方法：正交变换编码(如 DCT)，子带编码； 空间域方法：统计分块编码； 模型方法：分形编码，模型基编码； 基于重要性：滤波，子采样，比特分配，向量量化； （3）混合编码。 有 JBIG，H261，JPEG ，MPEG 等技术标准。 本实验主要利用 MATLAB 程序进行离散余弦变换（ DCT）压缩和行程编码（Run Length Encoding， RLE） 。 1) 离散余弦变换(DCT)图像压缩原理 离散余弦变换 DCT 在图像压缩中具有广泛的应用，它是 JPEG、MPEG 等数据压缩标准的 重要数学基础。 和相同图像质量的其他常用文件格式(如 GIF(可交换的图像文件格式)，TIFF(标签图像 文件格式)，PCX( 图形文件格式)相比，JPEG 是目前静态图像中压缩比最高的。 JPEG 比其 他几种压缩比要高得多，而图像质量都差不多(JPEG 处理的图像只有真彩图和灰度图) 。正 是由于其高压缩比，使得 JPEG 被广泛地应用于多媒体和网络程序中。JPEG 有几种模式， 其中最常用的是基于 DCT 变换的顺序型模式，又称为基本系统(Baseline)。 用 DCT 压缩图像的过程为： (1)首先将输入图像分解为 88 或 1616 的块，然后对每个子块进行二维 DCT 变换。 (2)将变换后得到的量化的 DCT 系数进行编码和传送，形成压缩后的图像格 式。 用 DCT 解压的过程为： (1)对每个 88 或 1616 块进行二维 DCT 反变换。 (2)将反变换的矩阵的块合成一个单一的图像。 余弦变换具有把高度相关数据能量集中的趋势，DCT 变换后矩阵的能量集中在矩阵的 左上角，右下的大多数的 DCT 系数值非常接近于 0。对于通常的图像来说，舍弃这些接近 于 0 的 DCT 的系数值，并不会对重构图像的画面质量带来显著的下降。所以，利用 DCT 变换进行图像压缩可以节约大量的存储空间。压缩应该在最合理地近似原图像的情况下使 用最少的系数。使用系数的多少也决定了压缩比的大小。 在压缩过程的第 2 步中，可以合理地舍弃一些系数，从而得到压缩的目的。在压缩 过程的第 2 步，还可以采用 RLE 和 Huffman 编码来进一步压缩。 2)行程编码（RLE）原理： 例如如下这幅 的二值图像， 如果采用行程编码可以按如下格式保存 其中 10 和 8 表示图像的宽和高。在这个小例子中行程编码并没有起到压缩图像的作用。 这是由于这个图的尺寸过小，当图像尺寸较大时行程编码还是不错的无损压缩方法。对于 灰度图像和二值图像，用行程编码般都有很高的压缩率。行程编码方法实现起来很容易， 对于具有长重复值的串的压缩编码很有效，例如：对于有大面积的阴影或颜色相同的图像， 使用这种方法压缩效果很好。很多位图文件格式都采用行程编码，如 TIFF，PCX，GEM，BMP 等。 3. 图像压缩编码的 MATLAB 程序语句 1） 利用 DCT 变换进行图像压缩的 MATLAB 程序 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 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0; %二值掩膜，用来压缩 DCT 系数，只留下 DCT 系数中左上角的 10 个 B2=blkproc(B,8 8, P1.*x.,mask); %只保留 DCT 变换的 10 个系数 I2= blkproc(B2,8,8,P1*x*P2,T,T); %逆 DCT，重构图像 Subplot(1,2,1); Imshow(I);title(原图像)； %显示原图像 Subplot(1,2,2); Imshow(I2);title(压缩图像)；% 显示压缩后的图像。对比原始图像和压缩后的图像， 虽然舍弃了 85%的 DCT 系数，但图像仍然清晰（当 然有一些质量损失） 3）利用行程编码（RLE）进行图像压缩 I=checkerboard(10,2); %调入原图像 m n=size(I); J=; for i=1:m value=I(i,1); num=1; for j=2:n if I(i,j)=value num=num+1; else J=J num value; num=1; value=I(i,j); end end I=J num value 0 0; %添加的行判断位 0 0 end disp(原图像大小：) whos(I); disp(压缩图像大小：) whos(J); disp(图像的压缩比：) disp(m*n/length(J) 、实验步骤 1 打开计算机，启动 MATLAB 程序； 2 调入“实验一”中获取的数字图像，并进行数据的行程（RLE）编码压缩处理； 3 将原图像在 Photoshop 软件中打开，分别以不同的位图文件格式进行“另保存” ，比 较它们的数据量。 4 记录和整理实验报告 、实验仪器 1 计算机； 2 MATLAB、Photoshop 等程序； 3 移动式存储器（软盘、U 盘等） 。 4 记录用的笔、纸。 、实验报告内容 1 叙述实验过程； 2 提交实验的原始图像和结果图像。 、思考题 1图像中哪些信息是主要的，哪些信息是次要的？ 2简述离散余弦变换（DCT）和行程编码（RLE）的原理。 实验四、图像融合 、实验目的 1 熟悉图像融合的意义和用途，理解图像融合的原理； 2 掌握图像融合的一般方法； 3 掌握运用 MATLAB 软件进行图像融合的操作。 、实验原理 图像融合的目的把来自多传感器的数据互补信息合并成一幅新的图像,以改善图像的质 量。图像融合最简单的理解就是两个（或多个）图像间的相加运算。这一技术广泛应用于 多频谱图像理解和医学图像处理等领域。主要分为空域和频域相加。 本实验主要应用 MATLAB 软件进行两幅图像的融合。方法有： 1 图像直接融合； 2 图像傅立叶变换融合； 3 图像小波变换融合。 图像融合的 MATLAB 程序如下： 1）调入、显示两幅图像的程序语句 load tartan; X1=X;map1=map; Load sinsin; X2=X;map2=map; %打开图像 Subplot(1 2 1) Image(X1),colormap(map1); Title(图像 map1) Subplot(1 2 2) Image(X2),colormap(map2); Title(图像 map2) %显示两幅图像 2）两幅图像直接融合的程序语句 figure,subplot(1 3 1) image(X1+X2)/2),colormap(map2); %在空域内直接融合 ti