Hi,
I am trying to overlay an AIA image on a SUVI image. I have used the SDO/LASCO example in the docs as a guide. Wherever it mentions lasco_map, i swapped this for suvi_map.
However, when I make the plot I can only see the AIA data and the both axis shows ranges from 0” to 8000”.
Any help would be appreciated.
Thanks, Richard
Can you share the data files you are using?
Hi Nabil,
The suvi data is:
time = a.Time(“2025-03-23 16:00”, “2025-03-23 17:00”)
query = Fido.search(time,a.Instrument.suvi,a.goes.SatelliteNumber(16),a.Wavelength(195*u.angstrom), a.Level.two)
Any file there will do. And the other is standard SDO data from the same time frame
Cheers,
Richard
So using the Lasco example, I dropped the custom WCS header and did the following:
import matplotlib.pyplot as plt
import astropy.units as u
from astropy.coordinates import SkyCoord
from sunpy.coordinates import SphericalScreen
import sunpy.map
from sunpy.net import Fido
from sunpy.net import attrs as a
time = a.Time("2025-03-23 16:00", "2025-03-23 16:00:15")
suvi_query = Fido.search(time,a.Instrument.suvi,a.goes.SatelliteNumber(16),a.Wavelength(195*u.angstrom), a.Level.two)
aia_query = Fido.search(time,a.Instrument.aia,a.Wavelength(171*u.angstrom))
aia_file = Fido.fetch(aia_query, site="NSO")
suvi_file = Fido.fetch(suvi_query)
aia_map = sunpy.map.Map(aia_file)
suvi_map = sunpy.map.Map(suvi_file)
with SphericalScreen(suvi_map.observer_coordinate):
aia_reprojected = aia_map.reproject_to(suvi_map.wcs)
fig = plt.figure()
ax = fig.add_subplot(projection=suvi_map)
suvi_map.plot(axes=ax)
aia_reprojected.plot(axes=ax, clip_interval=(1, 99.9)*u.percent, alpha=0.5)
plt.show()
fig = plt.figure()
ax = fig.add_subplot(projection=suvi_map)
suvi_map.plot(axes=ax)
aia_map.plot(axes=ax, clip_interval=(1, 99.9)*u.percent, alpha=0.5)
plt.show()
So the last two images are:
Not sure if this is helpful. Or what you wanted. I am unsure I can explain the difference in behaviour between your output and mine.
That’s great. I can work with the first plot! Thanks for coming up with a solution.
