Overplotting HEK feature/event polygons on a map#
How to overplot HEK outlines on a map.
import matplotlib.pyplot as plt
import numpy as np
import astropy.units as u
from astropy.coordinates import SkyCoord
from astropy.time import TimeDelta
from sunpy.coordinates import frames
from sunpy.net import attrs as a
from sunpy.net import hek
from sunpy.physics.differential_rotation import solar_rotate_coordinate
from sunpy.time import parse_time
We start with the sample data.
aia_map = sunpy.map.Map(sunpy.data.sample.AIA_171_IMAGE)
Look for coronal holes detected using the SPoCA feature recognition method:
Let’s find the biggest coronal hole within 80 degrees north/south of the equator:
Next let’s get the boundary of the coronal hole.
The coronal hole was detected at different time than the AIA image was taken so we need to rotate it to the map observation time.
/home/docs/checkouts/readthedocs.org/user_builds/sunpy/conda/stable/lib/python3.10/site-packages/sunpy/physics/differential_rotation.py:206: SunpyUserWarning: Using 'time' assumes an Earth-based observer.
warn_user("Using 'time' assumes an Earth-based observer.")
Now let’s plot the rotated coronal hole boundary on the AIA map, and fill it with hatching.
Total running time of the script: (0 minutes 1.196 seconds)