Saving and loading sunpy Maps with FITS#
In this example we are going to look at how to save and load a
GenericMap as FITS files.
FITS is a modern file format designed to meet the needs of the astronomy community. It has deep integration with Python, SunPy and Astropy as well as IDL and many other languages.
import astropy.units as u import sunpy.data.sample import sunpy.map
We begin by creating an
AIAMap object using the
We can now save this object to a FITS file to use later. Saving it like this
allows us to preserve all of the metadata of the object along with the actual
array data. When we load the FITS file again, we get an identical
All changes to
GenericMap are saved within the FITS file. For example,
if we rotate the image by 45 degrees and then save it to a FITS file, these changes
are reflected in both the image data and metadata when we load the data back in.
aia_map = aia_map.rotate(45*u.deg) # Please be aware that if you try to save this twice, # it will throw an exception rather than overwriting the file. aia_map.save('aia_map.fits')
This FITS file is portable and can be safely loaded by any FITS reader.
Map, we can easily read our rotated image back in:
aia_map_from_fits = sunpy.map.Map('aia_map.fits') aia_map_from_fits.peek(clip_interval=(1, 99.99)*u.percent)
Total running time of the script: (0 minutes 1.398 seconds)