LAB03

clear all;
close all;

a = imread("lena1.bmp");
b = imread("lena2.bmp");
c = imread("lena3.bmp");
d = imread("lena4.bmp");
n = 256;

figure('Name', 'Lena - histogramy');
subplot(4,2,1);
imshow(a);
subplot(4,2,2);
imhist(a, n);
subplot(4,2,3);
imshow(b);
subplot(4,2,4);
imhist(b, n);
subplot(4,2,5);
imshow(c);
subplot(4,2,6);
imhist(c, n);
subplot(4,2,7);
imshow(d);
subplot(4,2,8);
imhist(d, n);

hist1 = imread("hist1.bmp");
rozc = imadjust(hist1);

figure('Name', 'Hist1');
subplot(2,2,1)
imshow(hist1)
subplot(2,2,2)
imhist(hist1,256)
subplot(2,2,3)
imshow(rozc)
subplot(2,2,4)
imhist(rozc,256)
%obraz stał sie znacznie bardziej przejrzysty, teraz można dostrzec co tak
%naprawdę kryło się w jakby się mogło wydawać jednym odcieniu szarości

[H,x] = imhist(hist1);
c = cumsum(H);
k = max(c)/max(H);
C2=c/k;
figure(3)
subplot(2,2,1)
bar(x,H)
subplot(2,2,2)
plot(c)
subplot(2,2,3)
plot(x,H)
hold on
plot(x,C2)

k2 = max(c)/255;
C3=c/k2;
C4=uint8(C3);
lut = intlut(hist1, C4);

[H2,x2] = imhist(lut);
F = cumsum(H2);
m = max(F)/max(H2);
F2=F/m;

figure(4)
subplot(2,2,1)
imshow(lut)
subplot(2,2,2)
imhist(lut)
subplot(2,2,3)
plot(x2,H2)
hold on
plot(x2,F2)

wyr = histeq(hist1,256);
clahe = adapthisteq(hist1);

figure(5)
subplot(3,2,1)
imshow(lut)
subplot(3,2,2)
imhist(lut)
subplot(3,2,3)
imshow(wyr)
subplot(3,2,4)
imhist(wyr)
subplot(3,2,5)
imshow(clahe)
subplot(3,2,6)
imhist(clahe)


hist2 = imread("hist2.bmp");
rozc2 = imadjust(hist2);
hist3 = imread("hist3.bmp");
rozc3 = imadjust(hist3);
hist4 = imread("hist4.jpg");
rozc4 = imadjust(hist4);

figure(6)
subplot(2,4,1)
imshow(hist2)
title("Oryginalny")
subplot(2,4,2)
imhist(hist2)
subplot(2,4,3)
imshow(rozc2)
title("Rozciaganie")
subplot(2,4,4)
imhist(rozc2)
subplot(2,4,5)
imshow(histeq(hist2))
title("Wyrównywanie HE")
subplot(2,4,6)
imhist(histeq(hist2))
subplot(2,4,7)
imshow(adapthisteq(hist2))
title("Wyrównywanie CLAHE")
subplot(2,4,8)
imhist(adapthisteq(hist2))

figure(7)
subplot(2,4,1)
imshow(hist3)
title("Oryginalny")
subplot(2,4,2)
imhist(hist3)
subplot(2,4,3)
imshow(rozc3)
title("Rozciaganie")
subplot(2,4,4)
imhist(rozc3)
subplot(2,4,5)
imshow(histeq(hist3))
title("Wyrównywanie HE")
subplot(2,4,6)
imhist(histeq(hist3))
subplot(2,4,7)
imshow(adapthisteq(hist3))
title("Wyrównywanie CLAHE")
subplot(2,4,8)
imhist(adapthisteq(hist3))

figure(8)
subplot(2,4,1)
imshow(hist4)
title("Oryginalny")
subplot(2,4,2)
imhist(hist4)
subplot(2,4,3)
imshow(rozc4)
title("Rozciaganie")
subplot(2,4,4)
imhist(rozc4)
subplot(2,4,5)
imshow(histeq(hist4))
title("Wyrównywanie HE")
subplot(2,4,6)
imhist(histeq(hist4))
subplot(2,4,7)
imshow(adapthisteq(hist4))
title("Wyrównywanie CLAHE")
subplot(2,4,8)
imhist(adapthisteq(hist4))

=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=--=-=-=-=-=-=-=-=

=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=--=-=-=-=-=-=-=-=

clear all;
close all;

load histogramZadany;
fobos=imread("phobos.bmp");

wyr1=histeq(fobos);
wyr2=histeq(fobos,histogramZadany);
rozc=imadjust(fobos);
CLAHE=adapthisteq(fobos);


figure;
subplot(2,2,1);
imshow(fobos);
title("Oryginal");
subplot(2,2,2);
imhist(fobos);
subplot(2,2,3);
imshow(wyr1);
title("Histeq");
subplot(2,2,4);
imhist(wyr1);

figure;
subplot(2,2,1);
imshow(fobos);
title("Oryginal");
subplot(2,2,2);
imhist(fobos);
subplot(2,2,3);
imshow(wyr2);
title("Histogram Zadany");
subplot(2,2,4);
imhist(wyr2);

figure;
subplot(3,2,1);
imshow(fobos);
title("Oryginal");
subplot(3,2,2);
imhist(fobos);
subplot(3,2,3);
imshow(rozc);
title("Imadjust");
subplot(3,2,4);
imhist(rozc);
subplot(3,2,5);
imshow(CLAHE);
title("CLAHE");
subplot(3,2,6);
imhist(CLAHE);