对于一些对比度较低的图像，使用直方图算法处理该类图像，可以提高图像的对比度。

主要的操作步骤：

计算每个灰度的个数n_k

计算每个灰度出现的概率，并绘制直方图 P(K) = n_k/(h*w)

_{计算每个灰度的累计概率 : P(K)的和S(K)}

_{直方图的均衡化 S1(K) = S(K) * 255}

MATLAB代码如下：

1. [filename,pathname]=uigetfile('*.*','select an image'); input_I=imread([pathname filename]);%读入图像 figure(1); imshow(input_I); title('原始图像');%显示原始图像 I_gray = rgb2gray(input_I); figure(2); subplot(251);imshow(I_gray); title('灰度图像'); subplot(252);imhist(I_gray); title('灰度图像的直方图');%显示灰度图像及直方图 [height,width]=size(I_gray); %---------------------------------------------------------------------------------- Num = zeros(1,256); for i = 1 : height for j = 1 : width Num(I_gray(i,j)+1) = Num(I_gray(i,j)+1) + 1; end end subplot(253);stem(Num); xlim([0 280]);%把横坐标的范围限定在0--280之间 title('计算所得灰度图像的直方图'); %------------------------------------------------------------------------------------ %计算每个灰阶出现的概率，并绘制直方图(归一化后) P = zeros(1,256); for i = 1:256 P(i) = Num(i)/(height*width); end subplot(254); bar(0:255,P,'b'); %绘制直方图 xlim([0 255]); %将横坐标限定在0~255 title('归一化直方图'); xlabel('灰阶'); ylabel('概率 '); %--------------------------------------------------------------------------------------- %在P后计算累计的归一化直方图 S = zeros(1,256); for i = 1 : 256 for j = 1 : i S(i) = S(i) + P(j); end end subplot(255); bar(0:255,S,'k'); %绘制直方图 xlim([0 255]); title('累计后的归一化直方图 '); xlabel('灰阶'); ylabel('概率 '); %直方图均衡化 IMG = zeros(height,width); for i = 1 : height for j = 1 : width IMG(i,j) = floor(S(I_gray(i,j)+1)*255); end end subplot(256); imshow(uint8(IMG)); %显示最终的图像 title('均衡后图像 '); subplot(257);imhist(IMG); title('均衡后的图像直方图');%显示灰度图像及直方图 %-------------------------------------------------------------------------------------------------------- M = zeros(1,256); for i = 1 : height for j = 1 : width M(IMG(i , j ) + 1 ) = M(IMG(i , j ) + 1 ) + 1; end end subplot(258); bar(0:255,M,'k'); %绘制直方图 xlim([0 255]); title('均衡后的每个灰阶的个数 '); Q = zeros(1,256) ; for i = 1:256 Q(i) = M(i) / (height * width); end subplot(259); bar(0:255,P,'b'); %绘制直方图 xlim([0 255]); %将横坐标限定在0~255 title('均衡后的归一化直方图'); xlabel('灰阶'); ylabel('概率 '); QP = zeros(1,256); for i = 1 : 256 for j = 1 : i QP(i) = QP(i) + Q(j); end end subplot(2,5,10); bar(0:255,QP,'k'); %绘制直方图 xlim([0 255]); title('累计后的归一化直方图 '); xlabel('灰阶'); ylabel('概率 ');

   效果图因为网速不好，发不出来！