def aperture_ring(radius_arcsec,wcs):
a_in=radius_arcsec-1.5
a_out=radius_arcsec
b_out=a_out*incl
ring_sky=SkyEllipticalAnnulus(position,a_in=a_in*u.arcsec,a_out=a_out*u.arcsec,b_out=b_out*u.arcsec,theta=PA*u.degree)
ring_pix=ring_sky.to_pixel(wcs=wcs)
return ring_pix
def aperture_ring(radius_in,radius_out,wcs):
a_in=radius_in
a_out=radius_out
b_out=a_out*incl
ring_sky=SkyEllipticalAnnulus(center,a_in=a_in*u.arcsec,a_out=a_out*u.arcsec,b_out=b_out*u.arcsec,theta=PA*u.degree)
ring_pix=ring_sky.to_pixel(wcs=wcs)
return ring_pix
def aperture_rings(a_in, a_out, wcs, cosi, pa):
ring_sky = SkyEllipticalAnnulus(positions=center,
a_in=a_in * u.arcsec,
a_out=a_out * u.arcsec,
b_out=a_out * cosi * u.arcsec,
theta=pa * u.degree)
ring_pix = ring_sky.to_pixel(wcs=wcs)
return ring_pix
fitsimage = '../test3/run3/maps/NGC_5258_local_1mom.fits'
model = fits_import(fitsimage)[1]
model_cut = cut_2d(model, center, size, wcs)[1]
# draw rings
a_in = 4 * u.arcsec
a_out = 8 * u.arcsec
b_out = a_out * 0.5
theta = 204 * u.deg
ring_sky = SkyEllipticalAnnulus(positions=center,
a_in=a_in,
a_out=a_out,
b_out=b_out,
theta=theta)
ring_pix = ring_sky.to_pixel(wcs_cut)
phi = 204
xmax = 360
xmin = 0
ymax = 360
ymin = 0
xcen = 180
ycen = 180
x = np.arange(0, xmax - xmin, 0.1)
y = np.tan(np.radians(phi - 90)) * (x - xcen) + ycen
fig = plt.figure()
ax = plt.subplot(121)
ax.imshow(data_cut, origin='lower', cmap='gist_rainbow', vmin=-150, vmax=200)
ring_pix.plot()