function f = adpmedian(g, Smax)

%ADPMEDIAN performs adaptive median filtering
%  F = ADPMEDIAN (G, SMAX) performs adaptive median filtering of
%  image G. The median filter starts at size 3 by 3 and iterates up
%  to size SMAX by SMAX. SMAX must be an odd integer greater than 1.
%SMAX must be an odd, positive integer grater than 1.

if (Smax <= 1) | (Smax/2 == round(Smax/2)) | (Smax ~= round(Smax))
    error('SMAX harus bilangan ganjil > 1')
end

[M, N] = size(g);

%Initial setup
f = zeros(M, N);
setelahDiproses = false(M, N);

%Begin filtering
for k = 3:2:Smax
    zmin = ordfilt2(g, 1, ones(k, k), 'symmetric');
    zmax = ordfilt2(g, k*k, ones(k, k), 'symmetric');
    zmed = medfilt2(g, [k, k], 'symmetric');
    prosesLevelB = (zmed > zmin) & (zmax > zmed) & ~setelahDiproses;
    zB = (g > zmin) & (zmax > g);
    outputZxy = prosesLevelB & zB;
    outputZmed = prosesLevelB & ~zB;
    f(outputZxy) = g(outputZxy);
    f(outputZmed) = zmed(outputZmed);
    
    setelahDiproses = setelahDiproses | prosesLevelB;
    if all(setelahDiproses(:))
        break;
    end
end

%Output zmed for any remaining unprocessed pixels. Note that this
%zmed was computed using a window of size Smax by Smax, which is
%the final value of k in the loop.

f(~setelahDiproses) = zmed(~setelahDiproses);