Aligning AIA and HMI Data with Reproject¶
This example shows how to reproject one image to be aligned with another image.
Here we use the reproject package to transform one image to the reference frame of another image. The reference frames are defined by the respective World Coordinate System (WCS) information. This is a very generic way of aligning data, and can be very accurate.
You will need
reproject v0.6 or higher installed.
See Auto-Aligning AIA and HMI Data During Plotting for an alternate approach to image alignment, where the images are aligned as they are plotted, with no modified map created.
import matplotlib.pyplot as plt import astropy.units as u import sunpy.data.sample import sunpy.map
We use the AIA image and HMI image from the sample data. For the HMI map, we use the special HMI color map, which expects the plotted range to be -1500 to 1500.
Plot both images side by side.
<matplotlib.image.AxesImage object at 0x7f22750fbe20>
We can now reproject the HMI image to the WCS of the AIA image. We are using
reproject_interp, however the slower but most accurate
reproject_exact would also work well here. The
reproject_exact function only works when reprojecting between
two WCSes with the same observer, which makes it well suited to aligning
out_hmi = map_hmi.reproject_to(map_aia.wcs)
Plot the images side by side.
<matplotlib.image.AxesImage object at 0x7f22723d81c0>
As both of these images are now on the same pixel grid we can directly plot them over one another, by setting the transparency of the HMI plot.
Total running time of the script: ( 0 minutes 3.575 seconds)