pro replace im = readfits("F:dataRMS.fits",Header) for i=0,n_elements(im[*,1])-1 do begin for j=0,n_elements(im[1,*])-1 do begin if finite(im[i,j]) eq 0 then begin im[i,j]=-10 endif endfor endfor WRITEFITS,"F:dataRMS1.fits", im, Header end
pro replace n=1l im = readfits("F:dataRMS.fits",Header) for i=0,n_elements(im[*,1])-1 do begin for j=0,n_elements(im[1,*])-1 do begin if finite(im[i,j]) ne 0 then begin n=n+1
endif endfor endfor ;n是用来统计所有值不为NULL的个数,下面会建立个数为n的数组a a=dblarr(n) m=0l im = readfits("F:dataRMS.fits",Header) for i=0,n_elements(im[*,1])-1 do begin for j=0,n_elements(im[1,*])-1 do begin if finite(im[i,j]) ne 0 then begin a[m]=im[i,j] m++