; ; Make a radial mask for SwRI eclipse images ; x0 = 851.4 y0 = 829.3 pixX = 1024 pixY = 256 rsun = 728.3 img = fltarr(pixX,pixY) loadct,3 for ii=0,pixX-1 do begin x2 = (1.*ii)-x0 x2 = x2*x2 for jj=0,pixY-1 do begin y2 = (1.*jj)-y0 y2 = y2*y2 r = sqrt(y2+x2)/rsun img(ii,jj)=r^7 end print,ii end tv,img*256. openw,1,'r7mask' writeu,1,img close,1 end