Rotating a Map

How to rotate a map.

import matplotlib.pyplot as plt

import astropy.units as u

import sunpy.data.sample
import sunpy.map

We start with the sample data

aia_map = sunpy.map.Map(sunpy.data.sample.AIA_171_IMAGE)

GenericMap provides the rotate method which accepts an angle. This returns a rotated map and does not rotate in place. The data array size is expanded so that none of the original data is lost due to clipping. Note that subsequent rotations are not compounded. The map is only rotated by the specified amount from the original maps orientation.

aia_rotated = aia_map.rotate(angle=30 * u.deg)

Out:

/home/docs/checkouts/readthedocs.org/user_builds/sunpy/conda/stable/lib/python3.8/site-packages/sunpy/image/transform.py:118: SunpyUserWarning: Input data has been cast to float64.
  warnings.warn("Input data has been cast to float64.", SunpyUserWarning)

Let’s now plot the results.

fig = plt.figure()
ax = plt.subplot(projection=aia_rotated)
aia_rotated.plot()
aia_rotated.draw_limb()
aia_rotated.draw_grid()
plt.show()
AIA $171 \; \mathrm{\mathring{A}}$ 2011-06-07 06:33:02

Total running time of the script: ( 0 minutes 3.073 seconds)

Gallery generated by Sphinx-Gallery