直方图均衡化函数（彩色图）

% 直方图均衡化函数function [ RGB ] = histEquSelf(image)RGB = image;subplot(121);
imshow(RGB);
title('彩色图');[R, C, K] = size(RGB); % 新增的K表示颜色通道数% 统计每个像素值出现次数
cnt = zeros(K, 256);
for i = 1 : Rfor j = 1 : Cfor k = 1 : Kcnt(k, RGB(i, j, k) + 1) = cnt(k, RGB(i, j, k) + 1) + 1;endend
endf = zeros(3, 256);
f = double(f); cnt = double(cnt);% 统计每个像素值出现的概率， 得到概率直方图
for k = 1 : Kfor i = 1 : 256f(k, i) = cnt(k, i) / (R * C);end
end% 求累计概率，得到累计直方图
for k = 1 : Kfor i = 2 : 256f(k, i) = f(k, i - 1) + f(k, i);end
end% 用f数组实现像素值[0, 255]的映射。 
for k = 1 : Kfor i = 1 : 256f(k, i) = f(k, i) * 255;end
end% 完成每个像素点的映射
RGB = double(RGB);
for i = 1 : Rfor j = 1 : Cfor k = 1 : KRGB(i, j, k) = f(k, RGB(i, j, k) + 1);endend
end% 输出
RGB = uint8(RGB);
subplot(122);
imshow(RGB);
title('彩色直方图均衡化');end

测试代码

A = imread('lena.jpg');
histEquSelf(a);

运行结果

**测试原图：**下载后将原图命名为lena.jpg,与函数和测试代码放在同一个文件夹，即可运行

灰度图像的直方图+均衡化

% 直方图自己绘制
a = imread('lena.jpg');
H=rgb2gray(a);
% H= imread('./image/6.jpg'); 
if length(size(H))>2H=rgb2gray(H);  
end[m,n]=size(H);  
p=zeros(1,256);  
for i=0:255  p(i+1)=length(find(H==i))/(m*n);  
end  
subplot(2,2,1);  
imshow(H);  
title('原图');  
subplot(2,2,2);  
bar(0:255,p,'b');  
title('原图直方图');  s=zeros(1,256);  
for i=1:256  for j=1:i  s(i)=p(j)+s(i);                  end  
end  a=round(s*255);  
b=H;  
for i=0:255  b(H==i)=a(i+1);                
end  
subplot(2,2,3);  
imshow(b)                            
title('均衡化后图像');  for i=0:255  GPeq(i+1)=sum(p(a==i));            
end  
subplot(2,2,4);  
bar(0:255,GPeq,'b'); title('均衡化后的直方图');  

效果演示：