Procedura di cleaning per cosmic ray

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