[rm]=rmedian(im2,L,U)

%Relaxed memdian filter
%im2 - input image
%L - Lower limit
%U-Upper limit
function [rm]=rmedian(im2,L,U)
im2=double(im2);
n3=im2;
[x y]=size(im2);
%set defaults
if nargin<3
    L=3;
    U=5;
end
nt=medfilt2(im2);
n1=medfilt2(im2,[L L]);
n2=medfilt2(im2,[U U]);
for i=1:x
    for j=1:y
        if (n3(i,j)~=n1(i,j) & (n3(i,j)~=n2(i,j)))
           n3(i,j)=nt(i,j);
        end
    end
end
rm=uint8(n3);