get_dark_channel.m

%ÇóÍźĎń°ľÍ¨ľŔ Jdark = min(min(Jc))
function Jdark = get_dark_channel(image, win_size)

[m, n, ~] = size(image);

Ir = image( : , : , 1);
Ig = image( : , : , 2);
Ib = image( : , : , 3);

%ŃĄČĄ´°żÚś¨ŇĺłőĘźťŻ
Irr = 1./zeros(m + (win_size-1), n + (win_size-1));
Igg = 1./zeros(m + (win_size-1), n + (win_size-1));
Ibb = 1./zeros(m + (win_size-1), n + (win_size-1));

radius_size = floor(win_size / 2);
Irr(radius_size : (m + radius_size-1), radius_size :(n + radius_size-1)) = Ir;
Igg(radius_size : (m + radius_size-1), radius_size :(n + radius_size-1)) = Ig;
Ibb(radius_size : (m + radius_size-1), radius_size :(n + radius_size-1)) = Ib;

Jdark = zeros(m, n);

for i = 1 : 1 : m
    for j = 1 : 1 : n
        Rmin = min(min( Irr(i : i + (win_size-1), j : j + (win_size-1)) ));
        Gmin = min(min( Igg(i : i + (win_size-1), j : j + (win_size-1)) ));
        Bmin = min(min( Ibb(i : i + (win_size-1), j : j + (win_size-1)) ));
        
        Jdark(i,j) = min(min(Rmin,Gmin),Bmin);
    end
end