Closed keflavich closed 2 years ago
We don't have any documentation on how to show PV diagram figures.
Super short version is:
pvfits = fits.open('extracted_pv.fits') ww = wcs.WCS(pvfits[0].header) ax = pl.subplot(111, projection=ww)
but there's a lot of subtlety in how to make it look good and change units etc.
https://docs.astropy.org/en/stable/visualization/wcsaxes/controlling_axes.html is part of the answer, https://github.com/radio-astro-tools/tutorials/blob/6810376c0353f0bdf3be2b9b7231c388e886adba/PVDiagramPlotting.ipynb is another part
also https://docs.astropy.org/en/stable/visualization/wcsaxes/slicing_datacubes.html
closed by #102
We don't have any documentation on how to show PV diagram figures.
Super short version is:
but there's a lot of subtlety in how to make it look good and change units etc.