实验二图像直方图及灰度变换.docxVIP

实验二 图像直方图及灰度变换网络142班 143010 文宗正一、实验目的与要求1．掌握图像灰度直方图的概念及其计算方法，编写灰度直方图统计程序。2．通过对图像直方图的分析，学习应用直方图法解决诸如图像二值化等具体问题。3．熟悉直方图均衡化的计算过程及其应用。4．掌握图像灰度变换技术，通过调整图像的对比度和亮度等参数，改善视觉效果。二、实验相关知识直方图是图像最基本的统计特征，是图像亮度分布的概率密度函数，反映了图像灰度值的分布情况。直方图是多种空间域处理技术的基础。直方图操作能有效地用于图像增强，如通过直方图均衡化处理，可使图像在整个灰度级范围内的分布均匀化，即在每个灰度级上都具有相同的像素点数，从而获得较好的视觉效果。另外，直方图固有的信息也可用在图像分割等其它图像处理的应用中。灰度级变换技术可用g(x,y)=T[f(x,y)]的形式表示，其中f(x,y)为输入图像，g(x,y)为输出图像，T是对图像f进行某种处理的操作。由于(x,y)处的g值仅由f在该点处的亮度决定，T也称为一个亮度或灰度级变换函数，它与位置坐标(x,y)无关，所以通常写成如下的简化形式s=T(r)，其中r和s分别表示图像f和g在相应点(x,y)的亮度。利用灰度变换可调整一幅图像的明暗、对比度等。三、实验内容1、编写一个图像灰度直方图统计函数my_imhist，选择一幅图像利用my_imhist显示其直方图，将结果与MATLAB图像处理工具箱中提供的灰度直方图函数imhist的处理结果进行比较，并在同一窗口中显示出来。源代码：function my_imhist(A)[m,n]=size(A);x=1:256;y=zeros(1,256);for i=1:m*n y(A(i)+1)=y(A(i)+1)+1;endbar(x,y);运行结果：利用以上编写的函数my_imhist或imhist，估算图像iris.tif中瞳孔的半径（以像素为单位）。代码：I=imread(C:\Users\宗正\Desktop\学习\2016图像实验指导书\iris.tif); A=imhist(I); sum=0;for i=1:20 sum=sum+A(i);end display(sqrt(sum));3、利用MATLAB提供的图像直方图修正函数histeq，选择一幅直方图不均匀的图像如pout.tif，对其进行直方图均衡化处理，显示处理前后的图像以及它们的灰度直方图，体会直方图均衡化算法的特点。源代码：function myfun4(A)[m,n]=size(A);I=zeros(m,n);b=zeros(1,256);for i=1:m*n b(A(i)+1)=b(A(i)+1)+1;endc=zeros(1,256);for i=1:256 c(i)=b(i)/(m*n);endfor i=2:256 c(i)=c(i-1)+c(i);endfor i=1:m*n I(i)=256*c(A(i));endimshow(uint8(I));figure(1);subplot(2,2,1);imshow(A);title(原图像);subplot(2,2,2);imshow(uint8(I));subplot(2,2,3);imhist(A);subplot(2,2,4);imhist(floor(uint8(I)));运行结果：4、按照教材68页上的公式（4.1.6） (4.1.6)以及公式，分别编程实现图像的分段线性灰度变换和非线性灰度变换。利用自己编写的代码对图像bottle.tif进行灰度调整，使拍摄主体瓶子的花纹更加清晰地表现出来，以改善视觉效果。将同一个图形窗口分成四个子窗口来分别显示灰度变换前、后的图像以及它们的灰度直方图，注上文字标题。并考虑利用MATLAB提供的imadjust函数是否能完成同样的工作。代码：I=imread(C:\Users\宗正\Desktop\学习\2016图像实验指导书\bottle.tif);[m,n]=size(I);J=zeros(m,n);a=255/3;b=2*a;M=255;c=a+20;d=M-c;for i=1:m*nif(I(i)=0I(i)a) J(i)=(c/a)*I(i);elseif(I(i)=aI(i)b) J(i)=(d-c)/(b-a)*(I(i)-a)+c;else J(i)= (d-b)/(d-a)*(I(i)-b)+d;endendendfigure(1);subplot(2,2,1);imshow(I);title(原图像);subplot(2,2,2);imshow(ui