如上图所示在fits数据中存在NaN值,不利于接下来的数据处理,下面的小程序,可以把NaN替换掉,例如替换为-10 pro replace im = readfits("F:dataRMS.fits",Header) for i=0,n_elements(im )-1 do begin for j=0,n_elements(im )-1 do begin   ...
pro zhang READCOL,'E:test.dat',F='F,F,F,F',x, y, C, D FITEXY, x, y, A, B, X_SIG= C, Y_SIG= D, sigma_A_B, chi_sq, q print,'Intercept=',A print,'Slope=',B print,'Intercept standard deviation=',sigma_A_B print,'Slope standard deviation=',sigma_A_B pri ...