function im_out = Contrast_Enhance_CUT( im,prog )

im = uint8(im);
H = imhist(im);
ind = find(H'>prog);
Lmin = ind(1)-1;
n=size(ind,2);
Lmax = ind(n)-1;
im = double(im);
wsp_a = 255/(Lmax-Lmin);
wsp_b = (-wsp_a*Lmin);
im_out = wsp_a*im+wsp_b;
im_out = uint8(im_out);

end

-----------------


function [ k1,k2,k3,k4 ] = Quality( im )

Sx = size(im,1);
Sy = size(im,2);
Lmin = double(min(min(im)));
Lmax = double(max(max(im)));
Lsr = sum(sum(im,2),1)/(Sx*Sy);
k1 = (Lmax-Lmin)/255;
k2 = (Lmax-Lmin)/Lsr;
k3 = (Lmax-Lmin)/(Lmax+Lmin);
k4 = 4*sum(sum((im-Lsr).^2,2),1)/(255^2*Sx*Sy);

end



--------------------


[ k1,k2,k3,k4 ] = Quality( im1 );
[ kCE_cut1,kCE_cut2,kCE_cut3,kCE_cut4 ] = Quality( imCE_cut1 );
disp(['k1 = ' num2str(k1)]);
disp(['k2 = ' num2str(k2)]);
disp(['k3 = ' num2str(k3)]);
disp(['k4 = ' num2str(k4)]);

disp(['kCE_cut1 = ' num2str(kCE_cut1)]);
disp(['kCE_cut2 = ' num2str(kCE_cut2)]);
disp(['kCE_cut3 = ' num2str(kCE_cut3)]);
disp(['kCE_cut4 = ' num2str(kCE_cut4)]);