WCSAxes is a framework for making plots of Astronomical data in Matplotlib. It was previously distributed as a standalone package (wcsaxes), but is now included in astropy.visualization.
The following is a very simple example of plotting an image with the WCSAxes package:
import matplotlib.pyplot as plt
from astropy.wcs import WCS
from astropy.io import fits
from astropy.utils.data import get_pkg_data_filename
filename = get_pkg_data_filename('galactic_center/gc_msx_e.fits')
hdu = fits.open(filename)[0]
wcs = WCS(hdu.header)
fig, ax = plt.subplots(subplot_kw=dict(projection=wcs))
ax.imshow(hdu.data, vmin=-2.e-5, vmax=2.e-4, origin='lower')
ax.grid(color='white', ls='solid')
ax.set(xlabel='Galactic Longitude', ylabel='Galactic Latitude')
The above example uses the partially-object oriented matplotlib interface (keeping explicit references to figure and axes objects), as opposed to the pyplot interface. Note that this needed to access some of the more advanced functionality of WCSAxes. An example of this usage is:
ax = fig.add_subplot(projection=wcs, label='overlays')
ax.imshow(hdu.data, vmin=-2.e-5, vmax=2.e-4, origin='lower')
ax.coords.grid(True, color='white', ls='solid')
ax.coords[0].set_axislabel('Galactic Longitude')
ax.coords[1].set_axislabel('Galactic Latitude')
overlay = ax.get_coords_overlay('fk5')
overlay.grid(color='white', ls='dotted')
overlay[0].set_axislabel('Right Ascension (J2000)')
overlay[1].set_axislabel('Declination (J2000)')
In the rest of this documentation we will assume that you have kept a reference to the axes object, which we will refer to as ax. However, we also note when something can be done directly with the pyplot interface.
WCSAxes supports a number of advanced plotting options, including the ability to control which axes to show labels on for which coordinates, overlaying contours from data with different coordinate systems, overlaying grids for different coordinate systems, dealing with plotting slices from data with more dimensions than the plot, and defining custom (non-rectangular) frames.
Using WCSAxes# Reference/API# astropy.visualization.wcsaxes Package# Functions# Classes# astropy.visualization.wcsaxes.frame Module# Classes#RetroSearch is an open source project built by @garambo | Open a GitHub Issue
Search and Browse the WWW like it's 1997 | Search results from DuckDuckGo
HTML:
3.2
| Encoding:
UTF-8
| Version:
0.7.4
