# Issues with overlaying images

**URL:** <https://community.openastronomy.org/t/issues-with-overlaying-images/1301>\
**Category:** SunPy\
**Tags:** sunpy, question\
**Created:** [October 2, 2025, 12:49pm UTC](https://community.openastronomy.org/t/issues-with-overlaying-images/1301 "2025-10-02T12:49:15Z")\
**Posts on this page:** 5\
**Page:** 1

<div class="post-metadata">

**Author:** ![Richardjmorton](https://dub1.discourse-cdn.com/flex005/user_avatar/community.openastronomy.org/richardjmorton/32/671_2.png) [@Richardjmorton](https://community.openastronomy.org/u/Richardjmorton)\
**Post date:** [October 2, 2025, 12:49pm UTC](https://community.openastronomy.org/t/issues-with-overlaying-images/1301/1 "2025-10-02T12:49:15Z")

</div>

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

---

<div class="post-metadata">

**Author:** ![nabobalis](https://dub1.discourse-cdn.com/flex005/user_avatar/community.openastronomy.org/nabobalis/32/14_2.png) [@nabobalis](https://community.openastronomy.org/u/nabobalis)\
**Post date:** [October 2, 2025, 5:00pm UTC](https://community.openastronomy.org/t/issues-with-overlaying-images/1301/2 "2025-10-02T17:00:49Z")

</div>

Can you share the data files you are using?

---

<div class="post-metadata">

**Author:** ![Richardjmorton](https://dub1.discourse-cdn.com/flex005/user_avatar/community.openastronomy.org/richardjmorton/32/671_2.png) [@Richardjmorton](https://community.openastronomy.org/u/Richardjmorton)\
**Post date:** [October 3, 2025, 8:49am UTC](https://community.openastronomy.org/t/issues-with-overlaying-images/1301/3 "2025-10-03T08:49:36Z")

</div>

Hi Nabil,

The suvi data is:

```auto
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

---

<div class="post-metadata">

**Author:** ![nabobalis](https://dub1.discourse-cdn.com/flex005/user_avatar/community.openastronomy.org/nabobalis/32/14_2.png) [@nabobalis](https://community.openastronomy.org/u/nabobalis)\
**Post date:** [October 3, 2025, 6:57pm UTC](https://community.openastronomy.org/t/issues-with-overlaying-images/1301/4 "2025-10-03T18:57:30Z")

</div>

So using the Lasco example, I dropped the custom WCS header and did the following:

```python
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:

 ![reproject](https://europe1.discourse-cdn.com/flex005/uploads/openastronomy/original/1X/b4e8c5940869dd32b43e2cbc2140970bba299361.png) ![mesh_mpl_aligment](https://europe1.discourse-cdn.com/flex005/uploads/openastronomy/original/1X/62dcf54b072ed41118e232671e36aea73990375d.png)

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.

---

<div class="post-metadata">

**Author:** ![Richardjmorton](https://dub1.discourse-cdn.com/flex005/user_avatar/community.openastronomy.org/richardjmorton/32/671_2.png) [@Richardjmorton](https://community.openastronomy.org/u/Richardjmorton)\
**Post date:** [October 3, 2025, 7:46pm UTC](https://community.openastronomy.org/t/issues-with-overlaying-images/1301/5 "2025-10-03T19:46:49Z")

</div>

That’s great. I can work with the first plot! Thanks for coming up with a solution.