obraz_ciemny = imread ("czarne.jpg");
obraz_jasny = imread ("biale.jpg");
obraz_wejsciowy = imread("zdjecie1.jpg");

korekcja = obraz_jasny - obraz_ciemny ;
obraz_wyjsciowy ( : , : , 1 ) = ( obraz_wejsciowy ( : , : , 1 ) - obraz_ciemny ( : , : , 1 ) ) .* (max(max( korekcja ( : , : , 1 ) ) ) ./ korekcja ( : , : , 1 ));
obraz_wyjsciowy ( : , : , 2 ) = (obraz_wejsciowy ( : , : , 2 ) - obraz_ciemny ( : , : , 2 ) ) .* (max(max( korekcja ( : , : , 2 ) ) ) ./ korekcja ( : , : , 2 )) ;
obraz_wyjsciowy ( : , : , 3 ) = ( obraz_wejsciowy (:,:,3) - obraz_ciemny ( : , : , 3 ) ) .* (max(max( korekcja ( : , : , 3 ) ) ) ./ korekcja ( : , : , 3)) ;
imwrite ( obraz_wyjsciowy , "out.jpg" , 'jpg' );