def test_draw_focal_planes(self):
""" Draw a DECam focal planes """
ra, dec = 45, -45
radius = 1.5
delta = 1.0
plt.figure()
# Use the Cassini projection (because we can)
m = Skymap(projection='cass',
lon_0=ra,
lat_0=dec,
celestial=False,
llcrnrlon=ra + 2 * radius,
urcrnrlon=ra - 2 * radius,
llcrnrlat=dec - 2 * radius,
urcrnrlat=dec + 2 * radius)
# Can plot individual fields
m.draw_focal_planes([ra + delta / 2], [dec - delta / 2], color='g')
# Or as arrays
m.draw_focal_planes([ra, ra - delta, ra - delta],
[dec, dec + delta, dec - delta],
color='r')
# Draw the grid lines
m.draw_parallels(np.linspace(dec - 2 * radius, dec + 2 * radius, 5),
labelstyle='+/-',
labels=[1, 0, 0, 0])
m.draw_meridians(np.linspace(ra - 2 * radius, ra + 2 * radius, 5),
labelstyle='+/-',
labels=[0, 0, 0, 1])
plt.title('DECam Focal Planes')
