% 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
版权声明:本文为wayne8088原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。