直方图规定化——GML和SML映射规则
GML和SML原理
GML规定化（组映射） 先列出图像灰度值i,j,计算原始直方图Pr(i),列出规定直方图Pz(j),计算原始累计直方图Pi和规定累计直方图Pj，然后按照Pj到Pi寻找最接近的值存储为j，当满足给定的值多于一个时，按惯例取最小值。确定变换关系为i->j,最后列出匹配直方图P(j)。 SML规定化（单映射） 先列出图像灰度值i,j,计算原始直方图Pr(i),列出规定直方图Pz(j),计算原始累计直方图Pi和规定累计直方图Pj，然后按照Pi到Pj寻找最接近的值存储为j，当满足给定的值多于一个时，按惯例取最小值。确定变换关系为i->j,最后列出匹配直方图P(j)。
映射结果与分析
GML（组映射）映射结果 GML映射分析 先列出图像灰度值i,j,计算原始直方图Pr(i),列出规定直方图Pz(j),计算原始累计直方图Pi和规定累计直方图Pj，然后按照Pj到Pi寻找最接近的值存储为j，当满足给定的值多于一个时，按惯例取最小值。确定变换关系为i->j,最后列出匹配直方图P(j)。 SML（单映射）映射结果 SML映射分析 先列出图像灰度值i,j,计算原始直方图Pr(i),列出规定直方图Pz(j),计算原始累计直方图Pi和规定累计直方图Pj，然后按照Pi到Pj寻找最接近的值存储为j，当满足给定的值多于一个时，按惯例取最小值。确定变换关系为i->j,最后列出匹配直方图P(j)。 两种映射方法对比分析
代码
GML映射
%Write by 长安 Rjex
clc ,clear all
A = [0.1 0.05 0.15 0.2 0.2 0.15 0.05 0.1]; %原直方图矩阵
C = [0 0.3 0 0.45 0 0 0.25 0]; %规定直方图矩阵
F=zeros(1,size(A,2)); %规定化后直方图矩阵
for i = 1:size(A,2) %B为A的累计直方图矩阵
if i == 1
B(i) = A (i);
else
B(i) = B(i-1) + A(i);
end
end
for i = 1:size(C,2)
if i == 1
D(i) = C (i);
else
D(i) = D(i-1) + C(i);
end
end
min=1;
flag2 = 1;
for i = 1:size(D,2)
if C(i)~=0
for j = flag2:size(B,2)
t=abs(D(i)-B(j));
if(t min=t; %min为求差最小绝对值 flag=j; %flag为所求的最小绝对值下标 end end for k = flag2:flag E(k) = i; %E为GML映射规则矩阵 end flag2 = flag+1; min=1; end end for i = 1:size(E,2) F(E(i))=F(E(i))+A(i); end j=1:size(A,2); subplot(1,2,1),stem(j-1,A(j),'fill','black'),title('原始直方图'),xlabel('r_{k}'),ylabel('p_{r}(r_{k})'); subplot(1,2,2),stem(j-1,F(j),'fill','black'),title('GML映射规则'),xlabel('s_{k}'),ylabel('p_{s}(s_{k})'); SML映射 %Write by 长安 Rjex clc ,clear all A = [0.1 0.05 0.15 0.2 0.2 0.15 0.05 0.1]; %原直方图矩阵 C = [0 0.3 0 0.45 0 0 0.25 0]; %规定直方图矩阵 F=zeros(1,size(A,2)); %规定化后直方图矩阵 for i = 1:size(A,2) %B为A的累计直方图矩阵 if i == 1 B(i) = A (i); else B(i) = B(i-1) + A(i); end end for i = 1:size(C,2) %D为C的累计直方图矩阵 if i == 1 D(i) = C (i); else D(i) = D(i-1) + C(i); end end min=1; flag = 1; for i = 1:size(B,2) for j = 1:size(D,2) if C(j)~=0 t=abs(B(i)-D(j)); if(t min=t; %min为求差最小绝对值 flag=j; %flag为所求的最小绝对值下标 end end end E(i) = flag; %E为CML映射规则矩阵 min=1; end for i = 1:size(E,2) %F为规定化后直方图矩阵 F(E(i))=F(E(i))+A(i); end j=1:size(A,2); subplot(1,2,1),stem(j-1,A(j),'fill','black'),title('原始直方图'),xlabel('r_{k}'),ylabel('p_{r}(r_{k})'); subplot(1,2,2),stem(j-1,F(j),'fill','black'),title('CML映射规则'),xlabel('s_{k}'),ylabel('p_{s}(s_{k})');