AGCWD Code

clear all; close all; dbstop if error img_path = './imgs/2.bmp'; I = im2double(imread(img_path)); %imshow(imread(img_path)); % rgb to HSV hsv_I = rgb2hsv(I); V = hsv_I(:,:,3); %figure;imshow(V); % calculate pdf [p, x] = imhist(V); pdf = p./sum(p); %lmax = find(p > 0, 1, 'last'); %lmin = find(p > 0, 1, 'first'); lmax = max(x(p > 0)); lmin = min(x(p > 0)); %% eq(6) weighting distribution adjusted_parameter = 0.75; pdf_max = max(pdf); pdf_min = min(pdf); pdfw = pdf_max * ((pdf - pdf_min)./(pdf_max - pdf_min)).^adjusted_parameter; pdf_sum = sum(pdfw); pdfw = pdfw ./ pdf_sum; % plot the adjusted pdf if 1 horz = linspace(0, 255, 256); figure; b1 = bar(horz, pdf, 'FaceColor',[1 0 0], 'EdgeColor',[1 0 0]); set(b1,'facealpha',0.8) %alpha(1); hold on b2 = bar(horz, pdfw, 'FaceColor',[0 1 0], 'EdgeColor',[0 1 0]); set(b2,'facealpha',0.2) axis([0 255 0 0.02]) %set(gca, 'xtick', 0:50:255); %set(gca, 'ytick',); end %% eq(7) cdfw = cumsum(pdfw); gamma = 1 - cdfw; gamma_mapping = lmax*(x./lmax).^gamma; % plot the curve if 1 figure; plot(0:255, 255*gamma_mapping); end %% restore image out_V = zeros(size(V)); for i=1:256 l = x(i); out_V(V == l) = gamma_mapping(i); end hsv_I(:,:,3) = out_V; rI = hsv2rgb(hsv_I); figure;imshow([I rI]);