数字图像本操作及灰度调整.doc

 数字图像处理实验（一） 专 业 班 级 学 号 姓 名 实 验 名 称 实验一 数字图像基本操作及灰度调整 一．实验目的 1.掌握读、写图像的基本方法； 2.掌握MATLAB语言中图像数据与信息的读取方法； 3.理解图像灰度变换处理在图像增强的作用； 4．掌握绘制灰度直方图的方法，理解灰度直方图的灰度变换及均衡化的方法。 二．实验基本原理 灰度变换 灰度变换是图像增强的一种重要手段，它常用于改变图象的灰度范围及分布，是图象数字化及图象显示的重要工具。 图像反转 灰度级范围为[0, L-1]的图像反转可由下式获得 对数运算：有时原图的动态范围太大，超出某些显示设备的允许动态范围，如直接使用原图，则一部分细节可能丢失。解决的方法是对原图进行灰度压缩，如对数变换： s = clog(1 + r)，c为常数，r≥ 0 幂次变换： 对比拉伸：在实际应用中,为了突出图像中感兴趣的研究对象，常常要求局部扩展拉伸某一范围的灰度值，或对不同范围的灰度值进行不同的拉伸处理，即分段线性拉伸： 其对应的数学表达式为： 直方图均衡化 灰度直方图的横坐标是灰度级，纵坐标是该灰度级出现的频度，它是图像最基本的统计特征。依据定义，在离散形式下， 用rk代表离散灰度级，用pr(rk)代表pr(r)，并且有下式成立： 式中：nk为图像中出现rk级灰度的像素数，n是图像像素总数，而nk/n即为频数。 直方图均衡化处理是以累积分布函数变换法为基础的直方图修正法。假定变换函数为 (a) Lena图像 (b) Lena图像的直方图 图1-1 Lena图像及直方图 当灰度级是离散值时，可用频数近似代替概率值，即 式中：l是灰度级的总数目，pr(rk)是取第k级灰度值的概率，nk是图像中出现第k级灰度的次数，n是图像中像素总数。 所以积分可以表示为下列累计分布函数(cumulative distribution function, CDF) 三．实验结果与分析 1. 熟悉MATLAB语言中对图像数据读取，显示等基本函数。 a.将MATLAB目录下work文件夹中的forest.tif图像文件读出. 对应源代码如下： [I ,map] = imread(forest.tif), Imshow(I,map) 结果如下图所示： 图3-1 尝试修改map颜色矩阵的值，再将图像显示出来，相关代码如下： [I ,map] = imread(forest.tif),Imshow(I,map/4) 对应图像颜色的变化如下： 图3-2 Imfinfo函数的使用，显示图片相关信息,其对应源代码如下： imfinfo(forest.tif,tif) 对应结果如下 ans = Filename: forest.tif FileModDate: 04-Dec-2000 13:57:58 FileSize: 124888 Format: tif FormatVersion: [] Width: 447 Height: 301 BitDepth: 8 ColorType: indexed FormatSignature: [73 73 42 0] ByteOrder: little-endian NewSubfileType: 0 BitsPerSample: 8 Compression: PackBits PhotometricInterpretation: RGB Palette StripOffsets: [17x1 double] SamplesPerPixel: 1 RowsPerStrip: 18 StripByteCounts: [17x1 double]