majorityfilter.m

function newm = majorityfilter(img, amt)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%takes an image of classifications and smooths based on the majority within
%a neighborhood amt.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%img = imread(img);
tx = size(img,1);
ty = size(img,2);

newm = zeros(tx,ty);

clusters = max(max(img));

for i = 1:tx;
    for j = 1:ty;
        m = zeros(clusters,1);
        for dx = -amt:amt;
            for dy= -amt:amt;
                if abs(tx/2 - (i + dx) + .5) < tx/2 && abs(ty/2 + .5 - (j + dy)) < ty/2;
                    m(img(i + dx, j + dy),1) = m(img(i + dx, j + dy), 1) + 1;
                end
            end
        end
        max2 = 1;
        for it = 1 : clusters;
            if m(it, 1) > m(max2,1);
                max2 = it;
            end
        end
        newm(i,j) = max2;
    end
end
end