matlab 直方图匹配（规定）代码实现

% I为输入灰度图像，href为目标图像直方图，
function I = hist_Match(I,href)
[m1,n1,c]=size(I);
if c==3
    error('请输入灰度图像');
end
% 累计概率
I_cum_pro_ = cumsum(imhist(I))/(m1*n1);
href_cum_pro = cumsum(href/sum(href));
% 遍历寻找映射灰度值,寻找原图与目标图累计概率最接近的灰度值，存储在match中
match=uint8(zeros(256,1));
for i=1:256
    a = I_cum_pro_(i);
    min = 1;
    index=1;
    for j = 1:256
        b = href_cum_pro(j);
       if abs(a-b) < min  
           min=abs(a-b);
           index = j;
       end
    end
    match(i)=index - 1;  %灰度级=直方图统计下标 - 1 
end
% 遍历原图像，更新规定化的灰度值
for i=1:m1
    for j=1:n1
        t=I(i,j);
        I(i,j) = match(t+1);
    end
end