直方图规定化的编程实现

题目

matlab代码

代码是自己码的&＃xff0c;用最简单粗暴的方法&＃xff08;毫无技巧可言&＃xff09;&＃xff0c;而且还不一定对。。。。仅供参考
里面我设置了一个ratio值&＃xff0c;是用来控制映射的。不同的ratio值结果会不一样噢&＃xff0c;大家注意~

%原始图像&＃xff0c;像素的最高亮度级为160
R&＃61;[1 7 8 9 10 11 14;5 2 6 7 14 12 15;3 4 7 8 6 9 11;2 1 4 7 8 8 9;8 4 5 9 11 12 10;8 10 11 15 16 10 13;13 6 9 16 13 12 10]*10;
%统计像素矩阵中各像素出现的次数
R1&＃61;tabulate(R(:)); %返回&＃xff1a;像素值 出现次数 占比%
R2&＃61;R1(:,1); %I1的第一行&＃xff0c;0到160
R3&＃61;R1(:,2)/49; %占比&＃xff0c;相当于I3&＃61;I1(:,3)/100
[m_r,n_r]&＃61;size(R2); %[160,1]
r&＃61;[];pr&＃61;[]; %for循环后&＃xff0c;筛选出图像中存在的亮度级&＃xff0c;r装像素亮度级&＃xff0c;pr装亮度级对应的占比
for i&＃61;1:m_rif R3(i,1)~&＃61;0r(end&＃43;1)&＃61;R2(i,1);pr(end&＃43;1)&＃61;R3(i,1);end
end
figure
subplot(3,1,1);
stem(r,pr);
axis([0 255 0 0.3]);
title("原图像的灰度直方图");
xlabel("r");ylabel("pr");%参考图像
Z&＃61;[10 9 11 12 10 11 10;9 10 11 12 11 12 11;9 11 10 12 13 14 13;12 12 10 13 14 15 13;11 13 14 16 16 14 15;14 12 12 15 16 15 14;13 12 13 16 13 14 15]*11;
Z1&＃61;tabulate(Z(:));
Z2&＃61;Z1(:,1);
Z3&＃61;Z1(:,2)/49;
[m_z,n_z]&＃61;size(Z2);
z&＃61;[];pz&＃61;[];
for i&＃61;1:m_zif Z3(i,1)~&＃61;0z(end&＃43;1)&＃61; Z2(i,1);pz(end&＃43;1)&＃61; Z3(i,1);end
end
subplot(3,1,2);
stem(z,pz);
axis([0 255 0 0.3]);
title("参考图像的灰度直方图");
xlabel("z");ylabel("pz");[m_i1,n_i1]&＃61;size(pr); %返回[1,16]
[m_j1,n_j1]&＃61;size(pz); %返回[1,8]ratio&＃61;0.9; %控制阈值&＃xff0c;如果当前的temp_v还够放90%的temp_r&＃xff0c;那就整个塞进去&＃xff0c;of course你也可以自己选比例值
v&＃61;[];num_r&＃61;1;num_z &＃61; 1; %“指针”初始化&＃xff0c;指向第一个位置
temp_r&＃61;pr(num_r);temp_v&＃61;pz(num_z); %初始化&＃xff0c;r和z都指向第一个填充位
num_r&＃61;num_r&＃43;1;num_z&＃61;num_z&＃43;1; %指针指向下一个位置while (num_z<&＃61;n_j1)if(temp_v<(temp_r*ratio)) %z在这个填充位的容量太小&＃xff0c;不适合放东西&＃xff0c;跳到下一个v&＃61;[v 0]; %v可以不打分号&＃xff0c;看看输出和自己计算的一不一样temp_v&＃61;temp_v&＃43;pz(num_z);num_z &＃61;num_z&＃43; 1;endif(temp_v>&＃61;(temp_r*ratio))if(temp_v&＃61;&＃61;(temp_r*ratio)) %刚好匹配v&＃61;[v temp_r];temp_r&＃61;pr(num_r);num_r&＃61;num_r&＃43;1;temp_v&＃61;pz(num_z);num_z &＃61;num_z&＃43;1;else %这个填充位置容量大&＃xff0c;看看能不能多放几个load&＃61;temp_r;while(temp_v>&＃61;(load*ratio))temp_r&＃61;load;if(num_r<&＃61;n_i1) %如果r还有下一位load&＃61;temp_r&＃43;pr(num_r);num_r&＃61;num_r&＃43;1;endendnum_r&＃61;num_r-1; %加多了&＃xff0c;回退一位v&＃61;[v temp_r];temp_v&＃61;temp_v-temp_r;if(num_r<&＃61;n_i1)temp_r&＃61;pr(num_r);num_r&＃61;num_r&＃43;1;endtemp_v&＃61;temp_v&＃43;pz(num_z);num_z&＃61;num_z&＃43;1;endendend
v&＃61;[v temp_r];subplot(3,1,3);stem(z,v);
axis([0 255 0 0.3]);
title("原图像经过直方图规定化的灰度直方图");
xlabel("v");ylabel("pv");

运行结果对比

ratio&＃61;0.7

ratio&＃61;0.8

ratio&＃61;0.85

ratio&＃61;0.9