数字图像处理四个实验报告带有源程序文档格式.docx

1、 索引图像（Indexed images） RGB图像（RGB images）（1） 亮度图像一幅亮度图像是一个数据矩阵，其归一化的取值表示亮度。若亮度图像的像素都是uint8类或uint16类，则它们的整数值范围分别是0，255和0，65536。若图像是double类，则像素取值就是浮点数。规定双精度型归一化亮度图像的取值范围是0，1 （2） 二值图像一幅二值图像是一个取值只有0和1的逻辑数组。而一幅取值只包含0和1的uint8类数组，在MATLAB中并不认为是二值图像。使用logical函数可以把数值数组转化为二值数组或逻辑数组。创建一个逻辑图像，其语法为： B=logical（A）其中，

2、B是由0和1构成的数值数组。要测试一个数组是否为逻辑数组，可以使用函数: islogical（c）若C是逻辑数组，则该函数返回1；否则，返回0。（3） 索引图像索引颜色通常也称为映射颜色，在这种模式下，颜色都是预先定义的，并且可供选用的一组颜色也很有限，索引颜色的图像最多只能显示256种颜色。一幅索引颜色图像在图像文件里定义，当打开该文件时，构成该图像具体颜色的索引值就被读入程序里，然后根据索引值找到最终的颜色。（4） RGB图像一幅RGB图像就是彩色像素的一个MN3数组，其中每一个彩色相似点都是在特定空间位置的彩色图像相对应的红、绿、蓝三个分量。按照惯例，形成一幅RGB彩色图像的三个图像常称

3、为红、绿或蓝分量图像。令fR，fG和fB分别代表三种RGB分量图像。一幅RGB图像就利用cat（级联）操作将这些分量图像组合成彩色图像： rgb_image=cat（3,fR,fG,fB）在操作中，图像按顺序放置。2、数据类和图像类型间的转化表1中列出了MATLAB和IPT为表示像素所支持的各种数据类。表中的前8项称为数值数据类，第9项称为字符类，最后一项称为逻辑数据类。工具箱中提供了执行必要缩放的函数（见表2）。以在图像类和类型间进行转化。表1-1 MATLAB和IPT支持数据类型名称描述double双精度浮点数，范围为uint8无符号8比特整数，范围为0 255uint16无符号16比特整

4、数，范围为0 65536uint32无符号32比特整数，范围为0 4294967295int8有符号8比特整数，范围为-128 127int16有符号16比特整数，范围为-32768 32767int32有符号32比特整数，范围为-2147483648 2147483647single单精度浮点数，范围为char字符logical值为0或1表1-2 格式转换函数将输入转化为有效的输入图像数据类im2uint8logical,uint8,uint16和doulbeim2uint16mat2graydouble,范围为0 1im2doubleim2bwuint8,uint16和double 下面给

5、出读取、压缩、显示一幅图像的程序（%后面的语句属于标记语句，编程时可不用输入）I=imread（原图像名.tif）; % 读入原图像,tif格式whos I % 显示图像I的基本信息imshow（I） % 显示图像% 这种格式知识用于jpg格式，压缩存储图像，q是0-100之间的整数imfinfo filename imwrite（I,filename.jpg,quality,q）;imwrite（I,filename.bmp）; % 以位图（BMP）的格式存储图像% 显示多幅图像，其中n为图形窗口的号数figure（n）, imshow（filenamegg=im2bw（ % 将图像转为二值

6、图像figure, imshow（gg） % 显示二值图像三、实验内容及步骤1利用imread（ ）函数读取一幅图像，假设其名为flower.tif，存入一个数组中；I=imread（flower.tif2利用whos 命令提取该读入图像flower.tif的基本信息； whos I % 显示图像I的基本信息 Name Size Bytes Class Attributes I 129x173x3 66951 uint8 3利用imshow（）函数来显示这幅图像； figure（1）,imshow（I） 4利用imfinfo函数来获取图像文件的压缩，颜色等等其他的详细信息； imfinfo f

7、lower.tifans = Filename: FileModDate:31-五月-2012 17:49:15 FileSize: 6777 Format:jpg FormatVersion: Width: 173 Height: 129 BitDepth: 24 ColorType:truecolor FormatSignature: NumberOfSamples: 3 CodingMethod:Huffman CodingProcess:Sequential Comment: 5利用imwrite（）函数来压缩这幅图象，将其保存为一幅压缩了像素的jpg文件,设为flower.jpg；语

8、法：imwrite（原图像，新图像，quality,q）, q取0-100。flower.jpg,10）;figure（2）, imshow（flower2.jpg,50）;figure（3）, imshow（ % 显示多幅图像，其中n为图形窗口的号数6同样利用imwrite（）函数将最初读入的tif图象另存为一幅bmp图像，设为flower.bmp。flower.bmp7用imread（）读入图像：Lenna.jpg 和camema.jpg； P=imread（Lenna.jpgQ=imread（camema.jpg %用imread（）读入图像：8用imfinfo（）获取图像Lenna.j

9、pg和camema.jpg 的大小； imfinfo Lenna.jpg 3957 159 152 8grayscale 1 imfinfo camema.jpg 4857 1549用figure,imshow（）分别将Lenna.jpg和camema.jpg显示出来，观察两幅图像的质量。figure（1）,imshow（P）;figure（2）,imshow（Q）;10用im2bw将一幅灰度图像转化为二值图像，并且用imshow显示出来观察图像的特征。A=imread（lily.tifB=im2bw（A）;% 将图像转为二值图像figure（1）, imshow（A）;figure（2）,i

10、mshow（B）;11将每一步的函数执行语句拷贝下来，写入实验报告，并且将得到第3、9、10步得到的图像效果拷贝下来。四、考核要点1、熟悉在MATLAB中如何读入图像、如何获取图像文件的相关信息、如何显示图像及保存图像等，熟悉相关的处理函数。2、明确不同的图像文件格式，由于其具体的图像存储方式不同，所以文件的大小不同，因此当对同一幅图像来说，有相同的文件大小时，质量不同。五、实验仪器与软件（1） PC计算机（2） MatLab软件/语言包括图像处理工具箱（Image Processing Toolbox） （3） 实验所需要的图片 总结：本次实验的目的是学习MATLAB数字图像处理的一些基本知

11、识，熟悉及掌握在MATLAB中能够处理哪些格式图像，如何读取图像，掌握如何利用MATLAB来获取图像的大小、颜色、高度、宽度等等相关信息，掌握如何在MATLAB中按照指定要求存储一幅图像的方法及图像间如何转化。实验比较容易实验。实验二 图像基本运算一、 实验目的1了解图像的算术运算在数字图像处理中的初步应用。2体会图像算术运算处理的过程和处理前后图像的变化。二、 实验原理图像的代数运算是图像的标准算术操作的实现方法，是两幅输入图像之间进行的点对点的加、减、乘、除运算后得到输出图像的过程。如果输入图像为A（x,y）和B（x,y），输出图像为C（x,y），则图像的代数运算有如下四种形式： C（x,

12、y） = A（x,y） + B（x,y） C（x,y） = A（x,y） - B（x,y） C（x,y） = A（x,y） * B（x,y） C（x,y） = A（x,y） / B（x,y）图像的代数运算在图像处理中有着广泛的应用，它除了可以实现自身所需的算术操作，还能为许多复杂的图像处理提供准备。例如，图像减法就可以用来检测同一场景或物体生产的两幅或多幅图像的误差。使用MATLAB的基本算术符（+、-、*、/ 等）可以执行图像的算术操作，但是在此之前必须将图像转换为适合进行基本操作的双精度类型。为了更方便地对图像进行操作，MATLAB图像处理工具箱包含了一个能够实现所有非稀疏数值数据的算术操

13、作的函数集合。下表列举了所有图像处理工具箱中的图像代数运算函数。 表2-1 图像处理工具箱中的代数运算函数函数名功能描述Imabsdiff两幅图像的绝对差值Imadd两幅图像的加法Imcomplement补足一幅图像Imdivide两幅图像的除法Imlincomb计算两幅图像的线性组合Immultiply两幅图像的乘法imsubtract两幅图像的减法使用图像处理工具箱中的图像代数运算函数无需再进行数据类型间的转换，这些函数能够接受uint8和uint16数据，并返回相同格式的图像结果。虽然在函数执行过程中元素是以双精度进行计算的，但是MATLAB工作平台并不会将图像转换为双精度类型。代数运算

14、的结果很容易超出数据类型允许的范围。例如，uint8数据能够存储的最大数值是255，各种代数运算尤其是乘法运算的结果很容易超过这个数值，有时代数操作（主要是除法运算）也会产生不能用整数描述的分数结果。图像的代数运算函数使用以下截取规则使运算结果符合数据范围的要求：超出数据范围的整型数据将被截取为数据范围的极值，分数结果将被四舍五入。例如，如果数据类型是uint8，那么大于255的结果（包括无穷大inf）将被设置为255。注意：无论进行哪一种代数运算都要保证两幅输入图像的大小相等，且类型相同。三、 实验步骤1图像的加法运算图像相加一般用于对同一场景的多幅图像求平均效果，以便有效地降低具有叠加性质

15、的随机噪声。直接采集的图像品质一般都较好，不需要进行加法运算处理，但是对于那些经过长距离模拟通讯方式传送的图像（如卫星图像），这种处理是必不可少的。在MATLAB中，如果要进行两幅图像的加法，或者给一幅图像加上一个常数，可以调用imadd函数来实现。imadd函数将某一幅输入图像的每一个像素值与另一幅图像相应的像素值相加，返回相应的像素值之和作为输出图像。imadd函数的调用格式如下：Z = imadd（X，Y）其中，X和Y表示需要相加的两幅图像，返回值Z表示得到的加法操作结果。rice.tifJ=imread（rice02.tifK=imadd（I,J）;subplot（2,2,1）;ims

16、how（I）;subplot（2,2,2）;imshow（J）;subplot（2,1,2）;imshow（K）;%图像的加法运算给图像的每一个像素加上一个常数可以使图像的亮度增加。RGB=imread（flower-2.tifRGB2=imadd（RGB,50）;subplot（1,2,1）;imshow（RGB）;subplot（1,2,2）;imshow（RGB2）;%亮度增加 RGB2=imadd（RGB,-50）;%亮度变暗两幅图像的像素值相加时产生的结果很可能超过图像数据类型所支持的最大值，尤其对于uint8类型的图像，溢出情况最为常见。当数据值发生溢出时，imadd函数将数据截取

17、为数据类型所支持的最大值，这种截取效果称之为饱和。为了避免出现饱和现象，在进行加法计算前最好将图像转换为一种数据范围较宽的数据类型。例如，在加法操作前将uint8图像转换为uint16类型。2图像的减法运算图像减法也称为差分方法，是一种常用于检测图像变化及运动物体的图像处理方法。图像减法可以作为许多图像处理工作的准备步骤。例如，可以使用图像减法来检测一系列相同场景图像的差异。图像减法与阈值化处理的综合使用往往是建立机器视觉系统最有效的方法之一。在利用图像减法处理图像时往往需要考虑背景的更新机制，尽量补偿由于天气、光照等因素对图像显示效果造成的影响。在MATLAB中，使用imsubtract函数

18、可以将一幅图像从另一幅图像中减去，或者从一幅图像中减去一个常数。imsubtract函数将一幅输入图像的像素值从另一幅输入图像相应的像素值中减去，再将这个结果作为输出图像相应的像素值。imsubtract函数的调用格式如下：Z = imsubtract（X,Y）；其中，Z是X-Y操作的结果。以下代码首先根据原始图像（如图2.4（a）所示）生成其背景亮度图像，然后再从原始图像中将背景亮度图像减去，从而生成图2.4（b）所示的图像：rice=imread（background=imopen（rice,strel（disk,15）;rice2=imsubtract（rice,background）;

19、imshow（rice）;imshow（rice2）;%原始图像、减去背景图像如果希望从图像数据I的每一个像素减去一个常数，可以将上述调用格式中的Y替换为一个指定的常数值，例如：Z = imsubtract（I,50）；减法操作有时会导致某些像素值变为一个负数，对于uint8或uint16类型的数据，如果发生这种情况，那么imsubtract函数自动将这些负数截取为0。为了避免差值产生负值，同时避免像素值运算结果之间产生差异，可以调用函数imabsdiff。imabsdiff将计算两幅图像相应像素差值的绝对值，因而返回结果不会产生负数。该函数的调用格式与imsubtract函数类似。Z = i

20、msubtract（rice,50）;imshow（Z）;%每一个像素减去一个常数3. 图像的乘法运算两幅图像进行乘法运算可以实现掩模操作，即屏蔽掉图像的某些部分。一幅图像乘以一个常数通常被称为缩放，这是一种常见的图像处理操作。如果使用的缩放因子大于1，那么将增强图像的亮度，如果因子小于1则会使图像变暗。缩放通常将产生比简单添加像素偏移量自然得多的明暗效果，这是因为这种操作能够更好地维持图像的相关对比度。此外，由于时域的卷积或相关运算与频域的乘积运算对应，因此乘法运算有时也被作为一种技巧来实现卷积或相关处理。在MATLAB中，使用immultiply函数实现两幅图像的乘法。immultiply

21、函数将两幅图像相应的像素值进行元素对元素的乘法操作（MATLAB点乘），并将乘法的运算结果作为输出图形相应的像素值。immulitply函数的调用格式如下：Z = immulitply（X,Y）其中，Z=X*Y。例如，以下代码将使用给定的缩放因子对图2.5（a）所示的图像进行缩放，从而得到如图2.5（b）所示的较为明亮的图像：I = imread（room.tifJ = immultiply（I,1.5）;%图像的乘法运算uint8图像的乘法操作一般都会发生溢出现象。Immultiply函数将溢出的数据截取为数据类型的最大值。为了避免产生溢出现象，可以在执行乘法操作之前将uint8图像转换为一

22、种数据范围较大的图像类型，例如uint16。4图像的除法运算除法运算可用于校正成像设备的非线性影响，这在特殊形态的图像（如断层扫描等医学图像）处理中常常用到。图像除法也可以用来检测两幅图像间的区别，但是除法操作给出的是相应像素值的变化比率，而不是每个像素的绝对差异，因而图像除法也称为比率变换。在MATLAB中使用imdivide函数进行两幅图像的除法。imdivide函数对两幅输入图像的所有相应像素执行元素对元素的除法操作（点除），并将得到的结果作为输出图像的相应像素值。imdivide函数的调用格式如下：Z = imdivide（X,Y）其中，Z=X/Y。例如，以下代码将图4所示的两幅图像进

23、行除法运算，请将这个结果和减法操作的结果相比较，对比它们之间的不同之处：A= imread（electric.tifI = double（A）;J= I * 0.43 + 90;B= uint8（J）;Ip = imdivide（B,A）;imshow（A）;imshow（Ip）;%除法操作5图像的四则代数运算可以综合使用多种图像代数运算函数来完成一系列的操作。例如，使用以下语句计算两幅图像的平均值：I2 = imread（imshow（I2）;K = imdivide（imadd（I,I2）,2）;%计算两幅图像的平均值：建议最好不要用这种方式进行图像操作，这是因为，对于uint8或uint

24、16数据，每一个算术函数在将其输出结果传递给下一项操作之前都要进行数据截取，这个截取过程将会大大减少输出图像的信息量。执行图像四则运算操作较好的一个办法就是使用函数imlincomb。函数imlincomb按照双精度执行所有代数运算操作，而且仅对最好的输出结果进行截取，该函数的调用格式如下：Z = imlincomb（A,X,B,Y,C）；其中，Z=A*X+B*Y+C。MATLAB会自动根据输入参数的个数判断需要进行的运算。例如，以下语句将计算Z=A*X+C：Z = imlincomb（A,X,C）而以下语句将计算Z=A*X+B*Y：Z = imlincomb（A,X,B,Y,）本次实验了解并掌握了图像处理的基本运算，通过实验掌握了用MATLAB对图像进行基本的运算。实验三 图像增强空域滤波一、 实验目的 进一步了解MatLab软件/语言，学会使用MatLab对图像作滤波处理，使学生有机会掌握滤波算法，体会滤波效果。了解几种不同滤波方式的使用和使用的场合，培养处理实际图像的能力，并为课堂教学提供配套的实