function med2img,img

n=size(img)
if n(0) ne 3 then return,img
out=img(*,*,0)
for i=0,n(1)-1 do for j=0,n(2)-1 do out(i,j)=median(reform(img(i,j,*)))
return,out

end