实验指导书2016汇总

1、数字信号处理实验指导书实验一、图像的傅立叶变换实验三、图像压缩实验四、图像融合实验一、图像的傅立叶变换实验目的1了解图像变换的意义和手段；2熟悉傅里叶变换的性质；3热练掌握FFT酌方法反应用；4通过实验了解二维频谱的分布特点；5通过本实验掌握利用 MATLAB编程实现数字图像的傅立叶变换。实验原理1应用傅立叶变换进行图像处理傅里叶变换是线性系统分析的一个有力工具，它能够定量地分析诸如数字化系统、 采样点、电子放大器、卷积滤波器、噪音和显示点等的作用。通过实验培养这项技能，将有助于 解决大多数图像处理问题。对任何想在工作中有效应用数字图像处理技术的人来说，把时间用在学习和掌握博里叶变换上是很有必

2、要的。2傅立叶（Fourier ）变换的定义对于二维信号，二维 Fourier变换定义为：u uJfBi ®逆变换of-二维离散傅立叶变换为:六、思考题1曲杓叽r靠逆变换图像的傅立叶变换与一维信号的傅立叶变换变换一样，有快速算法，具体参见参考书目，有关傅立叶变换的快速算法的程序不难找到。实际上，现在有实现傅立叶变换的芯片，可以 实时实现傅立叶变换。3利用MATLAB软件实现数字图像傅立叶变换的程序：gif '%读入原图像文件%显示原图像%二维离散傅立叶变换%直流分量移到频谱中心%取傅立叶变换的实部%取傅立叶变换的虚部l=imread(原图像名imshow(I);fftI=ff

3、t2(I);sfftl=fftshift(fftl);RR=real(sfftI); ll=imag(sfftl);A=sqrt(RRA2+l|A2);% 计算频谱幅值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并利用

4、自编的函数显示其频谱；C)讨论不同的图像内容与FFT频谱之间的对应关系。4记录和整理实验报告。四、实验仪器1计算机；2 MATLAB 程序；3移动式存储器(软盘、U盘等)。4记录用的笔、纸。五、实验报告内容1傅里叶变换有哪些重要的性质 ?2图像的二维频谱在显示和处理时应注意什么实验目的1 掌握灰度直方图的概念及其计算方法；2 熟练掌握直力图均衡化和直方图规定化的计算过程；3 熟练掌握空域滤波中常用的平滑和锐化滤波器；4 掌握色彩直方图的概念和计算方法5利用MATLAB序进行图像增强。实验原理图像增强是指按特定的需要突出一幅图像中的某些信息， 同时，消弱或去除某些不需要 的信息的处理方法。其主要

5、目的是处理后的图像对某些特定的应用比原来的图像更加有效。 图像增强技术主要有直方图修改处理、 图像平滑化处理、 图像尖锐化处理和彩色处理技术等。 本实验以直方图均衡化增强图像对比度的方法为主要内容， 其他方法同学们可以在课后自行 联系。直方图是多种空间城处理技术的基础。 直方图操作能有效地用于图像增强。 除了提供有 用的图像统计资料外， 直方图固有的信息在其他图像处理应用中也是非常有用的， 如图像压 缩与分割。 直方图在软件中易于计算，也适用于商用硬件设备，因此， 它们成为了实时图像 处理的一个流行工具。直方图是图像的最基本的统计特征， 它反映的是图像的灰度值的分布情况。 直方图均衡化的目的是

6、使图像在整个灰度值动态变化范围内的分布均匀化，改善图像的亮度分布状态， 增强图像的视觉效果。灰度直方图是图像预处理中涉及最广泛的基本概念之一。图像的直方图事实上就是图像的亮度分布的概率密度函数， 是一幅图像的所有象素集合 的最基本的统计规律。 直方图反映了图像的明暗分布规律， 可以通过图像变换进行直方图调 整，获得较好的视觉效果。即在每个灰度直方图均衡化是通过灰度变换将一幅图像转换为另一幅具有均衡直方图， 级上都具有相同的象素点数的过程。下面给出直方图均衡化增强图像对比度的 MATLAB 程序：I=imread( 原图像名 .gif); % 读入原图像 J=histeq(I); %对原图像进行

7、直方图均衡化处理 Imshow(I); % 显示原图像 Title( 原图像 '); %给原图像加标题名 Figure； imshow(J); %对原图像进行屏幕控制；显示直方图均衡化后的图像 Title( 直方图均衡化后的图像 '); %给直方图均衡化后的图像加标题名 Figure; subplot(1,2,1) ;% 对直方图均衡化后的图像进行屏幕控制；作一幅子图作为并排 两幅图的第 1 幅图Imhist(I,64); %将原图像直方图显示为 64 级灰度 Title( 原图像直方图 '); %给原图像直方图加标题名 Subplot(1,2,2); % 作第 2 幅

8、子图Imhist(J,64) ; % 将均衡化后图像的直方图显示为 64 级灰度 Title( 均衡变换后的直方图 '); %给均衡化后图像直方图加标题名 处理后的图像直方图分布更均匀了， 图像在每个灰度级上都有像素点。 从处理前后的图 像可以看出， 许多在原始图像中看不清楚的细节在直方图均衡化处理后所得到的图像中都变 得十分清晰。实验步骤1 打开计算机，启动 MA TLAB 程序；程序组中“ work ”文件夹中应有待处理的图像文 件；2 调入“实验一”中获取的数字图像，并进行计算机均衡化处理；3 显示原图像的直方图和经过均衡化处理过的图像直方图。4 记录和整理实验报告四、实验仪器1

9、 计算机；2 MATLAB 程序；3 移动式存储器(软盘、 U 盘等)。4 记录用的笔、纸。五、实验报告内容121234直方图是什么概念？它反映了图像的什么信息？ 直方图均衡化是什么意思？它的主要用途是什么？实验三、图像压缩实验目的理解有损压缩和无损压缩的概念； 理解图像压缩的主要原则和目的； 了解几种常用的图像压缩编码方式。利用MATLAB序进行图像压缩。实验原理1.图像压缩原理图像压缩主要目的是为了节省存储空间， 增加传输速度。 图像压缩的理想标准是信息丢 失最少， 压缩比例最大。 不损失图像质量的压缩称为无损压缩， 无损压缩不可能达到很高的 压缩比； 损失图像质量的压缩称为有损压缩， 高

10、的压缩比是以牺牲图像质量为代价的。 压缩 的实现方法是对图像重新进行编码，希望用更少的数据表示图像。信息的冗余量有许多种，如空间冗余，时间冗余，结构冗余，知识冗余，视觉冗余等， 数据压缩实质上是减少这些冗余量。高效编码的主要方法是尽可能去除图像中的冗余成分，从而以最小的码元包含最大的图像信息。 编码压缩方法有许多种， 从不同的角度出发有不同的分类方法， 从信息论角度出发可分 为两大类。（1） .冗余度压缩方法，也称无损压缩、信息保持编码或嫡编码。具体说就是解码图像 和压缩编码前的图像严格相同，没有失真，从数学上讲是一种可逆运算。（ 2）信息量压缩方法，也称有损压缩、失真度编码或烟压缩编码。也就

11、是说解码图像 和原始图像是有差别的，允许有一定的失真。3 类：编码， Lempel zev 编码。应用在多媒体中的图像压缩编码方法，从压缩编码算法原理上可以分为以下 （1）无损压缩编码种类 哈夫曼（ Huffman ）编码，算术编码，行程（ RLE ）（ 2）有损压缩编码种类 预测编码， DPCM ，运动补偿； 频率域方法：正交变换编码 （如 DCT） ，子带编码； 空间域方法：统计分块编码； 模型方法：分形编码，模型基编码； 基于重要性：滤波，子采样，比特分配，向量量化；（ 3）混合编码。DCT ）压缩和行程编码（Run Length有 JBIG，H261，JPEG，MPEG 等技术标准。

12、本实验主要利用 MA TLAB 程序进行离散余弦变换（Encoding， RLE）。1）离散余弦变换 （DCT） 图像压缩原理JPEG、 MPEG 等数据压缩标准的离散余弦变换 DCT 在图像压缩中具有广泛的应用，它是 重要数学基础。和相同图像质量的其他常用文件格式（如GIF（可交换的图像文件格式）,TIFF（标签图像文件格式 ）， PCX（ 图形文件格式 ）相比， JPEG 是目前静态图像中压缩比最高的。 JPEG 比其 他几种压缩比要高得多，而图像质量都差不多 （JPEG 处理的图像只有真彩图和灰度图）。正是由于其高压缩比，使得 JPEG 被广泛地应用于多媒体和网络程序中。 JPEG 有几

13、种模式， 其中最常用的是基于 DCT 变换的顺序型模式，又称为基本系统 （Baseline） 。用 DCT 压缩图像的过程为：（1）首先将输入图像分解为 8X 8或16X 16的块，然后对每个子块进行二维DCT变换。（2）将变换后得到的量化的 DCT 系数进行编码和传送，形成压缩后的图像格 式。用 DCT 解压的过程为：（1）对每个8X 8或16X 16块进行二维 DCT反变换。（2）将反变换的矩阵的块合成一个单一的图像。 余弦变换具有把高度相关数据能量集中的趋势， DCT 变换后矩阵的能量集中在矩阵的左上角，右下的大多数的 DCT 系数值非常接近于 0。对于通常的图像来说，舍弃这些接近 于

14、0 的 DCT 的系数值，并不会对重构图像的画面质量带来显著的下降。所以，利用DCT变换进行图像压缩可以节约大量的存储空间。 压缩应该在最合理地近似原图像的情况下使用 最少的系数。使用系数的多少也决定了压缩比的大小。在压缩过程的第 2 步中，可以合理地舍弃一些系数，从而得到压缩的目的。在压缩 过程的第 2 步，还可以采用 RLE 和 Huffman 编码来进一步压缩。2) 行程编码(RLE原理： 例如如下这幅的二值图像,0000111 000001110 MX t tooQOOO111000 moLnoooo OILIOOOOO OllLOOOOOO lllOOOOOOO如果采用行程编码可以按

15、如下格式保存10807130d1306130613061306130fl1306130e对于灰对使BMP其中10和8表示图像的宽和高。在这个小例子中行程编码并没有起到压缩图像的作用。 这是由于这个图的尺寸过小，当图像尺寸较大时行程编码还是不错的无损压缩方法。:度图像和二值图像，用行程编码一般都有很高的压缩率。行程编码方法实现起来很容易，于具有长重复值的串的压缩编码很有效，例如：对于有大面积的阴影或颜色相同的图像，用这种方法压缩效果很好。很多位图文件格式都采用行程编码，如TIFF，PCX GEM等。3.图像压缩编码的 MATLAB程序语句1)利用DCT变换进行图像压缩的MATLAB程序RGB =

16、 imreadC 原图像名.tif);I = rgb2gray(RGB);J = dct2(l);imshow(log(abs(J),), colorma p(jet(64), colorbarJ(abs(J) < 10) = 0;K = idct2(J);figure,imshow(I)2)利用离散余弦变换进行JPEG图像压缩I=imread(原图像名.tif '%读入原图像；I=im2double(I);T=dctmtx(8);B=blk proc(l,8figure,imshow(K,0 255)%将原图像转为双精度数据类型；%产生二维DCT变换矩阵8, ' 1*x

17、*P2 ',T,T' %计算二维 DCT，矩阵T及其转置 T'是DCT函数P1*x*P2的参数Mask= 1 1 1 11 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 00 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)

18、; %逆 DCT ，重构图像 Subplot(1,2,1);Imshow(I);title( Subplot(1,2,2);Imshow(I2);title(原图像 '）； %显示原图像压缩图像 '）； %显示压缩后的图像。对比原始图像和压缩后的图像， 虽然舍弃了 85%的 DCT 系数， 但图像仍然清晰 （当然 有一些质量损失）3)利用行程编码(I=checkerboard(10,2); %调入原图像 m n=size(I);J=;for i=1:mvalue=I(i,1);num=1;for j=2:nif I(i,j)=value num=num+1;elseRLE ）进

19、行图像压缩J=J num value; num=1; value=I(i,j);endendI=J num value 0 0;%添加的行判断位 0 0enddisp( 原图像大小： ')whos(I');disp( 压缩图像大小： ')whos(J');disp( 图像的压缩比： ')disp(m*n/length(J)实验步骤1 打开计算机，启动 MA TLAB 程序；2 调入“实验一”中获取的数字图像，并进行数据的行程（ RLE ）编码压缩处理；3 将原图像在 Photoshop 软件中打开，分别以不同的位图文件格式进行“另保存” ，比 较它们的数

20、据量。4 记录和整理实验报告四、实验仪器1 计算机；U 盘等）。2 MATLAB 、Photoshop 等程序；3 移动式存储器（软盘、4 记录用的笔、纸。五、实验报告内容1 叙述实验过程；2 提交实验的原始图像和结果图像。六、思考题1图像中哪些信息是主要的，哪些信息是次要的？2简述离散余弦变换（ DCT ）和行程编码（ RLE ）的原理。实验四、图像融合实验目的1 熟悉图像融合的意义和用途，理解图像融合的原理；2 掌握图像融合的一般方法；3 掌握运用 MA TLAB 软件进行图像融合的操作。实验原理图像融合的目的把来自多传感器的数据互补信息合并成一幅新的图像 , 以改善图像的质 量。图像融合

21、最简单的理解就是两个 (或多个)图像间的相加运算。这一技术广泛应用于多 频谱图像理解和医学图像处理等领域。主要分为空域和频域相加。本实验主要应用 MA TLAB 软件进行两幅图像的融合。方法有：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); % 在空域内直接融合title( 两图像直接