nray=720*5 nraym1=nray-1 acirc = !PI * 2.0 / float(nray) dp = findgen(nray) * acirc + acirc * (0.25 * (float(nray))) dpx = intarr(nray) dpy = intarr(nray) xcen = 255.5 & ycen = 255.5 & mk3rad = 104.4435 & radius=1.00 rnder = 0.5005 dpx = fix ( cos(dp) * mk3rad * radius + xcen + rnder ) ; For disk size dpy = fix ( sin(dp) * mk3rad * radius + ycen + rnder ) ; For disk size mask = bytarr(512,512) mask(*,*) = 1 for radius = 0.0 , 1.1223, 0.001 do begin dpx = fix ( cos(dp) * mk3rad * radius + xcen + rnder ) ; For disk size dpy = fix ( sin(dp) * mk3rad * radius + ycen + rnder ) ; For disk size mask (dpx,dpy) = 0 endfor tv,bytscl(mask,min=0,max=1) openw,1,'mk3_innermask' writeu,1,mask close,1 end