数字图像处理（MATLAB）（第二版） 冈萨雷斯 学习笔记

数字图像处理（MATLAB）（第二版） 冈萨雷斯
主要用来记录自己的学习过程，相关实验内容，便于日后查找相关知识点及复习所用
第二章 灰度变换和空间滤波
本章重点讨论两类空间域处理方法：亮度（或灰度）变换与空间滤波。空间滤波有时又称为邻域处理或空间卷积。
2.1 背景
空间域技术直接对图像像素进行操作。 本章空间域处理表达式：g(x,y)=T[f(x,y)]
f(x,y)为输入图像，g(x,y)为输出图像 T是在点(x,y)的指定邻域上定义的对图像f进行处理的算子。 此外，T还可以对一组图像进行处理，如为了降低噪声而叠加K幅图像。

2.2灰度变换函数
由于输出值仅取决于某一点的灰度值，而不取决于该点的邻域，因此灰度变换函数通常写成如下形式：s=T(r)
式中，r表示图像f中相应点(x,y)的灰度，s表示图像g中相应点(x,y)的灰度。

2.2.1 函数imadjust 和stretchlim
函数imadjust是一个基本的图像处理工具箱函数，用于对灰度级图像进行灰度变换，其语法格式为：
g=imadjust(f,[low_in high_in],[low_out high_out],gamma)

对于函数imadjust，除了f和gamma,其所有输入值都被限定在0和1之间，而与f的类别无关。
函数mat2gray可用于将一幅图像转换为double类，并将其灰度标定到范围[0,1]，这与输入图像的类别无关。

例2.1 使用函数imadjust
f=imread(‘Fig0203(a).tif’);
imshow(f)

上图（a）是一幅数字乳房图像f,图中显示了一处病灶。
g1=imadjust(f,[0 1],[1 0]);
imshow(g1)

上图是使用g1命令得到的明暗反转图像（负片图像），这种过程对于增强嵌入在一大片黑色区域中的白色或灰色细节非常有用。例如在图b中就很容易分析乳房的组织。
图片的负片同样可以利用工具箱函数imcomplement得到：
g=imcomplement(f);
imshow(g)

命令：g2=imadjust(f,[0.5 0.75],[0 1]);
将0.5到0.75之间的灰度扩展到整个[0.1]范围，这种类型的处理对于强调感兴趣灰度区非常有用。

命令:g3=imadjust(f,[ ],[ ],2);
通过压缩灰度级的低端并扩展高端[见图d]，得到类似于图c的结果（增加了更多的灰色调）。

有时能够自动地使用函数imadjust而不必关心上面讨论的低参数或高参数是非常有用的，可使用函数stretchlim，基本语法为
low_High=stretcchlim(f)

g=imadjust(f,stretchlim(f),[ ]);
imshow(g)

g=imadjust(f,stretchlim(f),[1 0]);
imshow(g)

观察图片对比度的提升。

2.2.2对数及对比度拉伸变换
对数和对比度拉伸变换是动态范围操作的基本工具，对数通过以下表达式实现：
g=c*log(1+f)
其中，c为常数，f是浮点数。伽马曲线的形状是可变的，对数函数的形状是固定的，对数变换的一项主要应用是压缩动态范围。

执行对数变换时，我们希望使得压缩值出现在显示的完整范围内，对8比特来说，最简方法是使用语句：
gs=im2uint8(mat2gray(g));
通过使用函数mat2gray会将值限定在[0,1]范围内，使用函数im2uint8会将值限定在[0,255]范围内，把该图像转化成uint8类。

s=T(r)=11+(mr)E