matlab数字图像处理编码实现直方图均衡化

matlab中直方图均衡化有对应的内置函数histeq函数，其具体形式为：

g=histeq(f,nlev)

参数f为输入图像，nlev是为输出图像设置的灰度级（默认为64），g为直方图均衡化后的输出图像。

下面是自写matlab代码实现histeq相同功能的函数及测试：

函数定义代码：

function img_out=hist_equal(img,nlev)%直方图均衡化函数，其中参数img是输入图像，nlev是设定输出图像灰度值数量
    img_out=img;
    img=uint16(img);%unit8类型转为uint16类型，避免uint8的最大数255加1时出现溢出错误
    h=imhist(img)/numel(img);%直方图数据归一化
    s=0:length(h)-1;%创建输出灰度值向量
    rk=0.0;%输入图像灰度概率求和的变量
    for i=1:length(h)
        rk=rk+h(i);
        s(i)=round(rk*(nlev-1));
    end
    [m,n,p]=size(img);
    for i=1:m
        for j=1:n
            for k=1:p
                img_out(i,j,k)=s(img(i,j,k)+1);
            end
        end
    end
end

相关测试代码：

img=imread('e.jpg');
img=rgb2gray(img);%将rgb图像转化为灰度图像

figure
subplot(321)
imshow(img)%原始图像
title('原始图像')

subplot(322)
imhist(img)%原始图像直方图
title('原始图像直方图')

img_out1=hist_equal(img,256);%直方图均衡化,调用自写函数
subplot(323)
imshow(img_out1)%直方图均衡化图像
title('直方图均衡化图像(自写函数）')

subplot(324)
imhist(img_out1)%直方图均衡化后的图像的直方图
title('直方图均衡化后的图像的直方图（自写函数）')

img_out2=histeq(img,256);%直方图均衡化,调用matlab内置函数histeq函数
subplot(325)
imshow(img_out2)%直方图均衡化图像
title('直方图均衡化图像（histeq函数）')

subplot(326)
imhist(img_out2)%直方图均衡化后的图像的直方图
title('直方图均衡化后的图像的直方图（histeq函数）')

运行结果：

 可以看到，自写的直方图均衡化函数与matlab内置的histeq函数的测试结果一模一样。

感谢阅读！