图像增强—灰度变换增强

成 绩 评 定 表 学生姓名 赵诗雅 班级学号 1303030403 专 业 电子信息工程 课程设计题目 图像空域增强算法设计— 灰度变换增强 评 语 组长签字： 成绩   日期 2016年7月18日           课程设计任务书 学 院 信息科学与工程 专 业 电子信息工程 学生姓名 赵诗雅 班级学号 1303030403 课程设计题目 图像空域增强算法设计—灰度变换增强 实践教学要求与任务： 本设计要求利用Matlab进行编程及仿真，仿真内容为图像空域增强算法设计——灰度变换增强。拟利用所学数字图象处理技术知识，在Matlab软件系统上来实现灰度变换增强，并且对程序进行测试。要求如下： （1）掌握课程设计的相关知识、概念、思路及目的。 （2）程序设计合理、能够正确运行且操作简单，可实施性强。 （3）掌握灰度变换增强的。 （4）能够利用线性变换，非线性变换对灰度图像进行变换。 工作与进度安排： 第一阶段（1-2天）：熟悉matlab编程环境，查阅相关资料； 第二阶段（2-3天）：算法设计； 第三阶段（2-3天）：编码与调试； 第四阶段（1-2天）：实验与； 第五阶段（1-2天）：编写文档。 指导教师： 2016年7月3日 专业负责人： 2016年7月4日 学院教学副院长： 2016年7月4日             摘 要 空域增强在数字图像处理中起到对图像灰度的拉伸、压缩变换的作用，目前这种方法在处理图像灰度值方面得到广泛的运用。图像变换增强是利用一系列的变换方法使图像的对比度得到提升，以达到增强图像的目的，便于观察，易于区分不同灰度的图像。 灰度变换的方法有线性灰度变换、分段线性灰度变换、非线性灰度变换。本文采用灰度变换技术，分别对灰度图像进行了线性灰度变换增强、非线性灰度增强和分段线性灰度增强。仿真结果表明灰度变换可使图像动态范围增大，图像对比度扩展，图像变清晰，特征明显；对于灰度局限在某一个很小范围内的数字图像，采用线性函数对图像的每一个像素进行线性扩展，扩大像素的对比度，将有效地改善视觉效果，达到了提高图像对比度，增强图像效果的目的，证明了图像变换增强在数字图像处理中的重要作用。 关键词：空域增强；灰度变换；线性灰度变换；非线性灰度变换 目 录 1 设计目的与要求    1 1.1 设计目的    1 1.2 设计要求    1 2 设计原理    1 2.1 线性灰度变换    1 2.2分段线性灰度变换    2 2.3非线性灰度变换    3 3 设计    3 3.1设计思想    3 3.2 设计流程    4 4 软件实现    5 4.1 线性灰度变换增强    5 4.2 分段线性灰度变换增强设    6 4.3 非线性灰度变换增强    7 5 仿真与结果分析    7 5.1 仿真    8 5.2 结果分析    10 结  论    11 参考文献    12 1 设计目的与要求 1.1 设计目的 空域增强在数字图像处理中起到对图像灰度的拉伸、压缩变换的作用，目前这种方法在处理图像灰度值方面得到广泛的运用。灰度变换是空域变换的一种重要方式，可使图像对比度扩展，图像清晰，特征明显。本实验中，我们就用灰度变换来对灰度图像进行处理。 1.2 设计要求 （1）熟练掌握 MATLAB软件的基本操作 （2）学会掌握 MATLAB软件的程序编程 （3）查阅相关资料并分析，掌握灰度变换增强的主要思路 （4）培养独立分析和解决问题的能力，学会撰写课程设计的总结 （5）善于总结和改进方案，提高可实施性和高效性。 2 设计原理 2.1 线性灰度变换 线性变换简单的线性灰度变换法可表示为g(x,y)=(d-c/b-a)[f(x,y)-a]+c，其中：b和a分别是输入图像亮度分量的最大值和最小值，d和c分别是输出图像亮度分量的最大值和最小值。经过线性灰度变化法，图像亮度分量的线性范围从[a,b]变化到[c,d]，如图2.1所示： 图2.1 线性拉伸示意图 若图像中大部分像素的灰度级分布在区间[a,b]内,M为原图的最大灰度级,只有很小一部分的灰度级超过了此区间，则为了改善增强效果，可以令： （2.1） 由于人眼对灰度级别的分辨能力有限，只有当相邻像素的灰度值相差到一定程度时才能被辨别出来。通过上述变换，图像中相邻像素灰度的差值增加，例如在曝光不足或过度的情况下，图像的灰度可能会局限在一个很小的范围内，这时得到的图像可能是一个模糊不清，似乎没有灰度层次的图像。采用线性变换对图像中每一个像素灰度作线性拉伸，将有效改善图像视觉效果。[1] 2.2分段线性灰度变换 为了突出图像中感兴趣的目标或灰度区间，相对抑制那些不感兴趣的灰度区间，可采用分段线性变换，它将图像灰度区间分成两段乃至多段分别作线性变换。进行变换时，把0～255整个灰度值区间分为若干线段，每一个直线段都对应一个局部的线性变换关系。 常用的三段线性变换如图2.2所示： 图2.2 分段线性拉伸示意图 图中，中间段的灰度得到拉伸，两端灰度被压缩公式如2.2所示： （2.2） 2.3非线性灰度变换 非线性拉伸不是对图像的整个灰度范围进行扩展，而是有选择的对某一灰度范围进行扩展，其他范围的灰度值则有可能被压缩。非线性拉伸利用变换函数的数学性质实现对不同灰度值区间的扩展与压缩。下面介绍其中一种非线性扩展法：对数变换。 对数变换，是指输出图像的像素点的灰度值与对应的输出图像的像素灰度值之间为对数关系,其一般公式为：g(x,y)=a+ln[f(x,y)+1]/blnc，式中a,b,c 都是可以选择的参数，式中f(x,y)+1是为了避免对0求对数，确保ln[f(x,y)+1]≥0。当f(x,y)=0时，ln[f(x,y0+1)，则y=a,则a为y轴上的截距,确定了变换曲线的初始位置的变换关系b,c，两个参数确定变换曲线的变换速率。对数变换扩展了低灰度区，压缩了高灰度区，能使低灰度区的图像较清晰地显示出来。[2] 3 设计方案 3.1设计思想 灰度变换可使图像对比度扩展，图像清晰，特征明显。它是图像增强的重要手段。在图像空间所进行的灰度变换是一种点处理方法，它将输入图像中的每个像素（x，y）的灰度值f（x，y），通过映射函数变换成输出图像中的灰度g（x，y），即 g（x，y）=T[ f（x，y）]                      （3.1） 根据函数的性质，灰度变换的方法有以下几种： （1）线性灰度变换。 （2）分段线性灰度变换。 （3）非线性灰度变换。 对于线性灰度变换和非线性灰度变换，是直接用确定的公式依次对每个像素进行处理，也称为直接线性变换。[3] 3.2 设计流程 如下图3.1所示为线性变换流程图。 图3.1 线性变换流程 如图3.1所示，先将输入的灰度图像f(x,y)的位置及标签设置好之后，选择线性函数对原灰度图像进行灰度拉伸的处理，增强其对比度，得到变换之后的图像并输出。图3.2为分段线性变换增强流程图。 图3.2 分段线性变换流程 由图3.2可知，先将输入的灰度图像f(x,y)的位置及标签设置好之后，再设置灰度值变化的范围，分别为0-70，70-180和180-255。这三个范围的变化函数不同，从而达到使图像的灰度值拉伸程度不同的目的。经过分段线性变换后，得到变换后的图像g(x,y)。图3.3为非线性变换增强流程。[4] 图 3.3 非线性变换流程 由图3.3可知，先将输入图像f(x,y)的位置和标签设置好之后，确定图像灰度拉伸和压缩的范围，按照要求设置参数，利用对数函数对原始灰度图像进行非线性变换得到输出图像g(x,y)。 4 软件实现 4.1 线性灰度变换增强 通过调用imread函数和imshow函数来实现线性灰度变换增强，其核心代码如下： I = imread('fish.jpg');            % 读入原图像 I = im2double(I);                  % 转换数据类型为double [M,N] = size(I);                  % 计算图像面积 figure(1);subplot(121);            % 打开新窗口，设置图像位置 imshow(I);                        % 显示原图像 title('原图像');                  % 定义标签 Fa = 6; Fb = -160; O = Fa .* I + Fb/255;              % 线性变换公式 figure(1);subplot(122); imshow(O); title('灰度变换后的图像'); 4.2 分段线性灰度变换增强 通过调用imread函数和imshow函数来实现分段线性灰度变换增强，其核心代码如下： X1=imread('girl.jpg');% 读入原图像 subplot(221),imshow(X1)；          %输出图像 title('原图像');                %在原始图像中加标题 f0=0;g0=0;                      %进行灰度级线性分段增强； f1=70;g1=30; f2=180;g2=230; f3=255;g3=255; subplot(223),plot([f0,f1,f2,f3],[g0,g1,g2,g3]); axis tight,xlabel('f'),ylabel('g'),title('灰度变换曲线');      %绘制灰度变换线 %增强对比度 r1=(g1-g0)/(f1-f0); b1=g0-r1*f0; r2=(g2-g1)/(f2-f1); b2=g1-r2*f1; r3=(g3-g2)/(f3-f2); b3=g2-r3*f2; [m,n]=size(X1);