Skip to content

Instantly share code, notes, and snippets.

@wafels
Created February 23, 2022 20:51
Show Gist options
  • Save wafels/6ca1c463c80a9b27d9c92b6e070411ff to your computer and use it in GitHub Desktop.
Save wafels/6ca1c463c80a9b27d9c92b6e070411ff to your computer and use it in GitHub Desktop.
A small example to show use of spherical screen and composite maps
#
# Use the helioprojective spherical screen to maintain the off-disk emission.
# Use composite map to overlay the two images.
#
import matplotlib.pyplot as plt
from sunpy.map import Map
from sunpy.coordinates import Helioprojective
eit = Map( "/Users/ireland/Data/jp2/TRACE/2008_11_12__08_36_22_272__SOHO_EIT_EIT_284.jp2")
trace = Map("/Users/ireland/Data/jp2/TRACE/2008_11_15__00_11_06_000__TRACE_TRACE_TRACE_195.jp2")
with Helioprojective.assume_spherical_screen(eit.observer_coordinate):
#eit_screen = eit.reproject_to(eit.wcs)
trace_screen = trace.reproject_to(eit.wcs)
comp = Map(eit, trace_screen, composite=True)
fig = plt.figure()
ax = plt.subplot(projection=eit)
comp.plot()
plt.show()
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment