实验二像直方图及灰度变换q络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、编写一个图像灰度直方图统计函数myjmhist ,选择一幅图像利用 my」mhist显示其直方图，将结果与MATLAB图像处理工具箱中提供的灰度直函数imhist的处理结果进行比较，并在同一窗口中显示出来源代码：function myJmhist(A) [m#n]=size(A);x=l:256;y=zeros(l,256);y(A(i)+l)=y(A(i) + l) + l;end bar(xfy);运行结果：lure 1文件 賴 查看(Y)插入⑴工具①桌面Q)窗帮助Od)q a a 久久 a 逼 x ▼ I s 口18001600140012001000800600400200050100150200250300t，就~ - fr\ m •2、利用以上编写的函数myjmhist或imhist，估算半径(以像素为单位X代码：I=imread (丨 C: \Users\^』l:\Desktop\学习\2016 图像实验指导 |5\iris . tif 1);A=imhist(工)； sum=0；for i=l:20sum=sum+A(i)；enddisplay(sqrt(sum))；ans =66.86553、利用MATLAB提供的图像直方图修正函数histeq ,选择一幅直方图不 均匀的图像如pouttif ,对其进行直方图均衡化处理，显示处理前后的图像以及 它们的灰度直方图，体会直方图均衡化算法的特点。
源代码：function myfun4(A) [m/n]=size(A);I=zeros(m,n);b=zeros(l,256);for i = l:m*nb(A(i) + l) = b(A(i) + l)+l;end c=zeros(l,256);for i = l:256c(i) = b(i)/(m*n);endfor i=2:256c(i)=c(i-l)+c(i);endfor i = l:m*nI(i)=256*c(A(i));end imshow(uint8(I)); figured);subplot(2,2,l);imshow(A);title(•原图像•);subplot(2,2,2);imshow(uint8(I));subplot(2/2,3)；imhist(A);subplot(2,2,4);imhist(floor(uint8(I)));运行结果：Q 1MkL •騙X) 7*0 助必 WOQfiO 麵咖□ u J 焱 •、•、O X • a □ 63 ■ off田a4、按照教材68页上的公式（4.1.6 ）(4.1.6)(c/a)f(x,y) 0< /(x, y)
利用自己编写的代码对像bottle.tif进行灰度调整，使拍摄主体瓶子的花纹更加清晰地表现出来，以改善视觉效果将同一个形窗口分成四个子并考虑利用MATLAB提供的imadjust函数是否能完成同样的工作 代码：I=imreadCC:\Users\宗正\Desktop\学习\2016 图像实验指导[m,n]=size ⑴；J=zeros(m,n);a = 255/3;b=2*a;M=255;c=a+20;d = M-c;for i=l:m*nelseJ(i)=(d-c)/(b-a)*(I(i)-a)+c;elseJ(i)= (d-b)/(d-a)*0(i)-b)+d;endend end figure(l); subplot(2,2，l);imshow(I);titleC原图像•);subplot(2,2,2);imshow(uint8(J));titleC变换后D;subplot(2,2,3)；imhist ⑴；subplot(2,2,4);imhist(uint8(J));运行结果：BIft斯D 曹釁⑦癦1*0助必費□⑽□ □ y ct •、••、o w z • a □ u ■ aa12000100008000GOOO40002000。