图像灰度直方图处理..docVIP

第15周上机实验报告 学号：912104520120 姓名：鲍春霖 1、 内容 （1）图像线性对比度增强，测试图Fig69.bmp （2）图像非线性调整（对数，指数），侧视图Fig69.bmp （3）图像直方图均衡化，测试图Fig61.bmp 2、 原理 1. 图像线性对比度增强程原理 把分界数值定在100，所有0- 2. 图像非线性调整（对数，指数）原理： 将图片进行指数变换的算子是：B(i,j)=C*log(A(i,j)+1)，通过这个算子可以对图像进行指数上的变换。其中 指数非线性调整算子：1.1^(C*(A1(i,j)-a))-1，方法同对数非线性调整正好是逆运算关系。其中C也是一个常数。 通过计算每一个灰度值的像素点个数，进而算出每一个灰度值像素点的分布概率，然后进行直方图均衡的 3、 方法及程序 1. 图像线性对比度增强程序 clear all; A=imread(D:\Fig69.bmp); [row,col]=size(A); A1=double(A); fmax=max(max(A1)); fmin=min(min(A1)); for i=1:row for j=1:col if (A1(i,j)=100) B(i,j)=A1(i,j)*2; else if(A1(i,j)100A1(i,j)=255) B(i,j)=(255-200)/(255-100)*(A1(i,j)-100)+200; end end end end B=uint8(B); figure(1); subplot(1,2,1),imshow(A),title(?-í?); subplot(1,2,2),imshow(B),title(??D???±è?è????í?); figure(2); subplot(1,2,1),imhist(A),title(?-í?); subplot(1,2,2),imhist(B),title(??D???±è?è????í?); 2. 图像非线性调整（对数，指数） A=imread(D:\Fig69.bmp) A1=double(A); C=30; [row,col]=size(A1); for i=1:row for j=1:col B(i,j)=C*log(A1(i,j)+1); end end subplot(2,2,1),imshow(A),title(?-í?); subplot(2,2,2),imshow(B),title(·???D?±???oóí???); subplot(2,2,3),imhist(A); subplot(2,2,4),imhist(B); A=imread(D:\Fig69.bmp) C=0.8; a=1; A1=double(A); [row,col]=size(A1); for i=1:row for j=1:col B(i,j)=1.1^(C*(A1(i,j)-a))-1; end end subplot(2,2,1),imshow(A),title(?-í?); subplot(2,2,2),imshow(B),title(??êy±???oóí???); subplot(2,2,3),imhist(A); subplot(2,2,4),imhist(B); 3. 图像直方图均衡化 %?±·?í??ùoa?ˉ clear all; I = imread(Fig61.bmp); [row,col] = size(I); figure(1) subplot(2,2,1) imshow(I)%??ê??-ê?í??? subplot(2,2,2) imhist(I)%??ê??-ê?í????±·?í? %??DD?????ò?èí3??; s = zeros(1,256);%í3???÷?ò?èêy??￡?12256???ò?è?? for i = 1:row for j = 1:col s(I(i,j) + 1) = s(I(i,j) + 1) + 1;%??ó|?ò?è?μ????μ?êyá????óò? end end %?????ò?è·?2??ü?è p = zeros(1,256); for i = 1:256 p(i) = s(i) / (row * col ); end