% BBHE --> Bi-Histogram Equalization:
%dum1 = 'D:\AMIT\MONTHS\8_JAN\ISCE_10\Our paper\Input Images\Couple\';
%dum2 = 'D:\AMIT\MONTHS\4_SEP\Active Contour implementation\MC_image\B & W\';
tic;
o_img = imread('D:\AMIT\MONTHS\8_JAN\ISCE_10\Our paper\Input Images\slide4b.bmp');
imwrite(o_img,'D:\AMIT\MONTHS\8_JAN\ISCE_10\Our paper\Input Images\Man\5.3.01_BPHEME.tiff');
imshow(o_img);
sz = size(o_img);
o_mean = round(mean(o_img(:)));
% HISTOGRAM
h_l = zeros(256,1);
h_u = zeros(256,1);
for i = 1:sz(1)
for j = 1:sz(2)
g_val = o_img(i,j);
if(g_val<=o_mean)
h_l(g_val+1) = h_l(g_val+1) + 1;
else
h_u(g_val+1) = h_u(g_val+1)+ 1;
end
end
end
% NORMALIZED HISTOGRAM OR PDF
nh_l = h_l/sum(h_l);
nh_u = h_u/sum(h_u);
% CDF
hist_l_cdf = double(zeros(256,1));
hist_u_cdf = double(zeros(256,1));
hist_l_cdf(1) = nh_l(1);
hist_u_cdf(1) = nh_u(1);
for k = 2:256
hist_l_cdf(k) = hist_l_cdf(k-1) + nh_l(k);
hist_u_cdf(k) = hist_u_cdf(k-1) + nh_u(k);
end
% IMAGE MODIFICATION
equalized_img = zeros(sz);
range_l = [0 o_mean];
range_u = [(o_mean+1) 255];
for i =1:sz(1)
for j =1:sz(2)
g_val = o_img(i,j);
if(g_val<=o_mean)
equalized_img(i,j) = range_l(1) + round(((range_l(2)-range_l(1))*hist_l_cdf(g_val+1)));
else
equalized_img(i,j) = range_u(1) + round(((range_u(2)-range_u(1))*hist_u_cdf(g_val+1)));
end
end
end
figure,imshow(uint8(equalized_img));
imwrite(uint8(equalized_img),'D:\AMIT\MONTHS\8_JAN\ISCE_10\Our paper\Input Images\City\5.3.02_BBHE.tiff');
toc;
评论0