实验一 灰度变换及直方图变换

《实验一 灰度变换及直方图变换》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、实验一图像空域增强——灰度变换一．实验内容：图像灰度变换二．实验目的：1、理解、掌握对数变换、幂次变换、对比度拉伸的原理方法2、了解图像灰度级的概率分布情况，掌握直方图均衡化原理和方法3、编程实现图像的对数变换、幂次变换、对比度拉伸、直方图均衡化三．实验步骤：1．获取实验用图像：Fig3.05(a).jpg.使用imread函数将图像读入Matlab；并对其进行对数变换（令c=1）；应用sublpot在同一个图像窗口中显示原图像及对数变换后的图像，分别用title标注。原理：对数变换的一般表达式为：其中是一个常数，并假设。该种变换

2、是一窄带低灰度输入图像映射为一宽带输出值。相应的是输入灰度的高调整值。可以利用这种变换来扩展被压缩的高值图像中的暗像素。我们用matlab编程获取了图像Fig3.05(a).jpg，并对其进行对数变换，所得结果如下图1所示：图1分析图像：比较上面两幅图像，我们会发现经过对数变换后的图像的细节可见程度比原始图像更明显些，在原始图像中最亮的像素显示为重点，灰度低的像素（恰恰是重要的）却观察不到。经过对数变换扩展灰度低的像素从而使暗像素观察的很清楚。实验结果分析：经过matlab编程很好地对原始图像进行了对数变换，并取得了很好的效果，达

3、到了实验地目的。程序如下：clearall;closeall;clc;f=imread('D:实验一灰度变换及直方图变换Fig3.05(a).jpg');s=log(1+double(f));subplot(1,2,1);imshow(f);title('原始图像');subplot(1,2,2);imshow(s);title('对数变换');1．产生灰度变换函数T1，使得：0.3rr<0.35s=0.105+2.6333(r–0.35)0.35≤r≤0.651+0.3(r–1)r>0.65用T1对原图像Fig3.10(b)

4、.jpg进行处理，在同一窗口显示处理前后的图像。原理：灰度拉伸又叫对比度拉伸，它是最基本的一种灰度变换，使用的是最简单的分段线性变换函数，它的主要思想是提高图像处理时灰度的动态范围，适用于低对比度图像的处理，增强对比度。我们用matlab编程获取了图像Fig3.10(b).jpg，并对其进行对比度拉伸，所得结果如下图2所示：图2分析图像：比较图2中的两幅图像，我们发现经过对比度拉伸处理后的图像明显比原始图像清晰。在原始图像中，花粉的对比度较低，经过灰度变换后的图像的对比度比较高。这是因为变换函数把灰度级有原范围线性地拉伸至饱和范围

5、，使得变换后的图像对比度显著增强。实验结果分析：运用matlab，对原始图像进行灰度变换，显著增强了图片的对比度，达到了实验目的。程序如下：clearall;closeall;clc;f=imread('D:实验一灰度变换及直方图变换Fig3.10(b).jpg');r=f;ifr<0.35s=0.3*r;elseifr<=0.65s=0.105+2.6333*(r-0.35);elses=1+0.3*(r-1);endendsubplot(1,2,1);imshow(f);title('原始图像');subplot(1,2,

6、2);imshow(s,[]),title('处理后图像');3.分别用s=r0.6;s=r0.4;s=r0.3对Fig3.08(a).jpg图像进行幂次变换。为简便起见，可使用Matlab中的imadjust函数；在同一图像窗口显示原图像，伽玛为0.6，0.4，0.3的图像原理：幂次变换的基本形式为：其中和为正常数。有时考虑到偏移量（即当输入为0时可测量输出），上式也可写为。当时把输入窄带暗值映射到宽带输出值，增强图像的对比度。当时，正好与时相反。根据此原理，我们用matlab进行编程来实现对图像Fig3.08(a).jpg进行

7、的幂次变换处理，其结果图像如下图3所示：图3分析图像：由于原始图像整体呈现暗色，我们知道当伽马小于1时，幂次变换把窄带低灰度输入映射到宽带输出值，此时，幂次变换相当于对数变换，只不过比对数变换更灵活些，因此，我们通过不断尝试伽马值来达到增强图像的目的。从图3中我们可观察到：在伽马值从0.6到0.4时，可看到更多的细节，当伽马值进一步减小到0.3时，背景中的细节进一步增强，但图像看上去开始有轻微的“冲淡”外观的那一点的对比度开始减少，尤其是背景。比较所有结果，我们发现对比度和可辨别细节的最好增强效果在伽马为0.4时。程序如下：cle

8、arall;closeall;clc;f=imread('D:实验一灰度变换及直方图变换Fig3.08(a).jpg');r=im2double(f);s1=power(r+10.^(-100),0.6);s2=power(r+10.^(-10