实验三 灰度变换和直方图处理.doc

实验三 灰度变换与直方图处理 一、实验目的 1、理解图像直方图的含义； 2、掌握MATLAB中灰度转换和直方图处理相关函数； 3、了解直方图均衡化的实现方法。 二、实验步骤及内容 1、灰度转换 灰度变换是图像增强技术中最简单的一类。常用的三类基本函数：线性函数、对数函数、幂次函数。 （1）图像反转：对教材中61页的图像Fig3.04(a)的图像进行反转变换 Image=imread(Fig3.04(a).jpg); imshow(Image); Image=255-double(Image); %反转 figure(2),imshow(uint8(Image)); （2）对数变换：对教材中62页的图像Fig3.05(a)进行对数变换 close all Image=imread(Fig3.05(a).jpg); imshow(Image); Image=log(1+double(Image)); %对数变换 figure(2),imshow(Image,[]); （3）幂次变换：对教材中66页的图像Fig3.09(a)进行幂次变换 close all Image=imread(Fig3.09(a).jpg); imshow(Image); Image2=double(Image).^3.0; %γ＝3 figure(2),imshow(Image2,[]); Image2=double(Image).^4.0; figure(3),imshow(Image2,[]); Image2=double(Image).^5.0; figure(4),imshow(Image2,[]); 2、灰度直方图的显示 imhist imhist是MATLAB提供的显示图像灰度直方图的函数。 (1) 显示图像cameraman.tif 的直方图。 close all x=imread(cameraman.tif); figure subplot(2,1,1) imshow(x) subplot(2,1,2) imhist(x) （2）还可以将直方图显示成其他形式，比如柱形图（续一）或者折线图（续二）。 （续一）： h=imhist(x); h1=h(1:10:256); horz=1:10:256; bar(horz,h1) axis([0 255 0 1500]) set(gca, xtick,0:50:255) set(gca, ytick,0:200:1500) （续二）： h=imhist(x); plot(h) axis([0 255 0 1500]) set(gca, xtick,0:50:255) set(gca, ytick,0:200:1500) 思考题1： 显示教材中71页，暗、亮、低对比度和高对比度图以及它们对应的直方图，即显示实验文件夹中Fig3.15(a)1.jpg ～Fig3.15(a)4.jpg各图的直方图，观察并说明其直方图的分布情况，进一步了解直方图的含义。 图2-1 暗、亮、低对比度和高对比度图 close all x=imread(Fig3.15(a)1.jpg); figure subplot(2,1,1) imshow(x) subplot(2,1,2) imhist(x) 分析：此图像较暗，观察图形，我们发现直方图分布在灰度值较小处。 close all x=imread(Fig3.15(a)2.jpg); figure subplot(2,1,1) imshow(x) subplot(2,1,2) imhist(x) 此图像较亮，观察图形，我们发现直方图分布在灰度值较大处。 close all x=imread(Fig3.15(a)3.jpg); figure subplot(2,1,1) imshow(x) subplot(2,1,2) imhist(x) 此图像对比度较小，观察图形，我们发现直方图比较集中地分布在灰度值的中心范围。 close all x=imread(Fig3.15(a)4.jpg); figure subplot(2,1,1) imshow(x) subplot(2,1,2) imhist(x) 此图像对比度较大，观察图形，我们发现直方图分布整个灰度值范围。 灰度直方图是灰度级的函数，它表示图象中具有每种灰度级的象素的个数，反映图象中每种灰度出现的频率imadjust （1）imadjust是MATLAB提供的图像强度调整函数。利用它也可以实现图像反转。 close all I=imread(Fig3.04(a).jpg); imshow(I,[]); JJ=imadjust(I,[0 1],