pro clean_cosmic,a on_error,2 dim=size(a) HS=dim(1) VS=dim(2) frameout=fltarr(HS+2,VS+2) frameout(1:HS,1:VS)=a frameout(1:HS,0)=a(*,1) frameout(1:HS,VS+1)=a(*,VS-1) frameout(0,1:VS)=a(1,*) frameout(VS+1,1:VS)=a(VS-1,*) frameout=sigma_filter(frameout) a=frameout(1:HS,1:VS) end