Note

Click here to download the full example code

Plotting a difference image

How to compute and plot a difference image. The example uses sunpy.map.MapSequence to compute a difference image and then plot it. This basic method works for base difference or running difference. Just change whether you’re subtracting the previous image or the first image in a sequence.

import matplotlib.pyplot as plt
import matplotlib.colors as colors
import astropy.units as u
from sunpy.net import Fido, attrs as a
import sunpy.map

First we’ll download a couple images and store them in a sunpy.map.MapSequence. This could be from any channel of any imager. Here, we use SDO/AIA 304 Å.

instrument = a.Instrument('AIA')
wave = a.Wavelength(30 * u.nm, 31 * u.nm)
result = Fido.search(a.Time('2015-06-18T00:00:00', '2015-06-18T00:00:10') |
                     a.Time('2015-06-18T01:03:30', '2015-06-18T01:03:35'),
                     instrument,
                     wave)
downloaded_files = Fido.fetch(result)
maps = sunpy.map.Map(downloaded_files, sequence=True)

Out:

/home/docs/checkouts/readthedocs.org/user_builds/sunpy/conda/stable/lib/python3.8/site-packages/parfive/downloader.py:86: DeprecationWarning: The loop argument is deprecated since Python 3.8, and scheduled for removal in Python 3.10.
  self.http_queue = asyncio.Queue(loop=self.loop)
/home/docs/checkouts/readthedocs.org/user_builds/sunpy/conda/stable/lib/python3.8/asyncio/queues.py:48: DeprecationWarning: The loop argument is deprecated since Python 3.8, and scheduled for removal in Python 3.10.
  self._finished = locks.Event(loop=loop)
/home/docs/checkouts/readthedocs.org/user_builds/sunpy/conda/stable/lib/python3.8/site-packages/parfive/downloader.py:87: DeprecationWarning: The loop argument is deprecated since Python 3.8, and scheduled for removal in Python 3.10.
  self.http_tokens = asyncio.Queue(maxsize=self.max_conn, loop=self.loop)
/home/docs/checkouts/readthedocs.org/user_builds/sunpy/conda/stable/lib/python3.8/asyncio/queues.py:48: DeprecationWarning: The loop argument is deprecated since Python 3.8, and scheduled for removal in Python 3.10.
  self._finished = locks.Event(loop=loop)
/home/docs/checkouts/readthedocs.org/user_builds/sunpy/conda/stable/lib/python3.8/site-packages/parfive/downloader.py:88: DeprecationWarning: The loop argument is deprecated since Python 3.8, and scheduled for removal in Python 3.10.
  self.ftp_queue = asyncio.Queue(loop=self.loop)
/home/docs/checkouts/readthedocs.org/user_builds/sunpy/conda/stable/lib/python3.8/asyncio/queues.py:48: DeprecationWarning: The loop argument is deprecated since Python 3.8, and scheduled for removal in Python 3.10.
  self._finished = locks.Event(loop=loop)
/home/docs/checkouts/readthedocs.org/user_builds/sunpy/conda/stable/lib/python3.8/site-packages/parfive/downloader.py:89: DeprecationWarning: The loop argument is deprecated since Python 3.8, and scheduled for removal in Python 3.10.
  self.ftp_tokens = asyncio.Queue(maxsize=self.max_conn, loop=self.loop)
/home/docs/checkouts/readthedocs.org/user_builds/sunpy/conda/stable/lib/python3.8/asyncio/queues.py:48: DeprecationWarning: The loop argument is deprecated since Python 3.8, and scheduled for removal in Python 3.10.
  self._finished = locks.Event(loop=loop)

Files Downloaded:   0%|          | 0/2 [00:00<?, ?file/s]/home/docs/checkouts/readthedocs.org/user_builds/sunpy/conda/stable/lib/python3.8/site-packages/aiohttp/connector.py:964: DeprecationWarning: The loop argument is deprecated since Python 3.8, and scheduled for removal in Python 3.10.
  hosts = await asyncio.shield(self._resolve_host(
/home/docs/checkouts/readthedocs.org/user_builds/sunpy/conda/stable/lib/python3.8/site-packages/aiohttp/connector.py:964: DeprecationWarning: The loop argument is deprecated since Python 3.8, and scheduled for removal in Python 3.10.
  hosts = await asyncio.shield(self._resolve_host(
/home/docs/checkouts/readthedocs.org/user_builds/sunpy/conda/stable/lib/python3.8/site-packages/aiohttp/locks.py:21: DeprecationWarning: The loop argument is deprecated since Python 3.8, and scheduled for removal in Python 3.10.
  self._event = asyncio.Event(loop=loop)


aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:   0%|          | 0.00/7.92M [00:00<?, ?B/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:   0%|          | 29.2k/7.92M [00:00<00:39, 201kB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:   0%|          | 0.00/7.95M [00:00<?, ?B/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:   1%|          | 69.2k/7.92M [00:00<00:35, 218kB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:   0%|          | 29.3k/7.95M [00:00<00:39, 201kB/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:   2%|1         | 119k/7.92M [00:00<00:31, 245kB/s] 


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:   1%|          | 72.2k/7.95M [00:00<00:35, 222kB/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:   2%|2         | 179k/7.92M [00:00<00:27, 278kB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:   2%|1         | 127k/7.95M [00:00<00:30, 253kB/s] 

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:   3%|3         | 251k/7.92M [00:00<00:23, 320kB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:   2%|2         | 192k/7.95M [00:00<00:26, 291kB/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:   4%|4         | 334k/7.92M [00:00<00:20, 368kB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:   3%|3         | 272k/7.95M [00:00<00:22, 338kB/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:   5%|5         | 431k/7.92M [00:01<00:17, 425kB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:   5%|4         | 370k/7.95M [00:00<00:19, 397kB/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:   7%|6         | 554k/7.92M [00:01<00:14, 499kB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:   6%|6         | 490k/7.95M [00:01<00:15, 469kB/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:   9%|8         | 697k/7.92M [00:01<00:12, 585kB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:   8%|7         | 630k/7.95M [00:01<00:13, 554kB/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:  11%|#         | 863k/7.92M [00:01<00:10, 684kB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:  10%|#         | 796k/7.95M [00:01<00:10, 654kB/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:  13%|#3        | 1.06M/7.92M [00:01<00:08, 803kB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:  12%|#2        | 986k/7.95M [00:01<00:09, 769kB/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:  16%|#6        | 1.28M/7.92M [00:01<00:07, 934kB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:  15%|#5        | 1.20M/7.95M [00:01<00:07, 897kB/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:  19%|#9        | 1.53M/7.92M [00:01<00:05, 1.08MB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:  18%|#8        | 1.45M/7.95M [00:01<00:06, 1.04MB/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:  23%|##2       | 1.80M/7.92M [00:02<00:04, 1.23MB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:  22%|##1       | 1.73M/7.95M [00:01<00:05, 1.20MB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:  26%|##5       | 2.03M/7.95M [00:02<00:04, 1.42MB/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:  26%|##6       | 2.10M/7.92M [00:02<00:04, 1.39MB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:  28%|##7       | 2.21M/7.95M [00:02<00:03, 1.53MB/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:  31%|###       | 2.42M/7.92M [00:02<00:03, 1.56MB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:  32%|###1      | 2.54M/7.95M [00:02<00:03, 1.76MB/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:  35%|###4      | 2.77M/7.92M [00:02<00:02, 1.75MB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:  35%|###4      | 2.77M/7.95M [00:02<00:02, 1.90MB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:  40%|####      | 3.18M/7.95M [00:02<00:02, 2.19MB/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:  40%|###9      | 3.16M/7.92M [00:02<00:02, 2.01MB/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:  44%|####4     | 3.50M/7.92M [00:02<00:01, 2.29MB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:  44%|####3     | 3.46M/7.95M [00:02<00:01, 2.34MB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:  50%|####9     | 3.94M/7.95M [00:02<00:01, 2.67MB/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:  48%|####8     | 3.84M/7.92M [00:02<00:01, 2.43MB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:  53%|#####3    | 4.25M/7.95M [00:02<00:01, 2.79MB/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:  53%|#####2    | 4.19M/7.92M [00:02<00:01, 2.68MB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:  60%|######    | 4.81M/7.95M [00:02<00:00, 3.18MB/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:  59%|#####8    | 4.65M/7.92M [00:03<00:01, 2.95MB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:  65%|######4   | 5.16M/7.95M [00:03<00:00, 3.27MB/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:  64%|######4   | 5.10M/7.92M [00:03<00:00, 3.28MB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:  73%|#######3  | 5.80M/7.95M [00:03<00:00, 3.70MB/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:  71%|#######1  | 5.64M/7.92M [00:03<00:00, 3.58MB/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:  78%|#######7  | 6.15M/7.92M [00:03<00:00, 3.93MB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:  78%|#######8  | 6.21M/7.95M [00:03<00:00, 3.77MB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:  87%|########6 | 6.91M/7.95M [00:03<00:00, 4.24MB/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:  86%|########5 | 6.79M/7.92M [00:03<00:00, 4.27MB/s]


aia_lev1_304a_2015_06_18t01_03_31_13z_image_lev1.fits:  93%|#########2| 7.37M/7.95M [00:03<00:00, 4.31MB/s]

aia_lev1_304a_2015_06_18t00_00_07_12z_image_lev1.fits:  93%|#########3| 7.40M/7.92M [00:03<00:00, 4.66MB/s]

                                                                                                           
Files Downloaded:  50%|#####     | 1/2 [00:05<00:05,  5.47s/file]


                                                                                                           
Files Downloaded: 100%|##########| 2/2 [00:05<00:00,  2.74s/file]
/home/docs/checkouts/readthedocs.org/user_builds/sunpy/conda/stable/lib/python3.8/asyncio/base_events.py:654: ResourceWarning: unclosed event loop <_UnixSelectorEventLoop running=False closed=False debug=False>
  _warn(f"unclosed event loop {self!r}", ResourceWarning, source=self)

Now we’ll do a standard plot of the second image just to see it.

plt.figure()
ax = plt.subplot(projection=maps[1])
maps[1].plot()
../../../_images/sphx_glr_difference_images_001.png

Out:

<matplotlib.image.AxesImage object at 0x7f3898a03c10>

And now we can do take the actual difference.

diff = maps[1].data - maps[0].data

But we have to decide what to do with the metadata. For example, what time does this difference image correspond to? The time of the first or second image? The mean time? You’ll have to decide what makes most sense for your application. Here we’ll just use the metadata from the second image. Then we can store the difference and header back in a Map.

meta = maps[1].meta
diff_map = sunpy.map.Map(diff, meta)

Finally, we’ll plot it. We’ll apply a colormap and renormalize the intensity so that it shows up well.

plt.figure()
ax_diff = plt.subplot(projection=diff_map)
diff_map.plot(cmap='Greys_r',
              norm=colors.Normalize(vmin=-50, vmax=50))
plt.show()

#sphinx_gallery_thumbnail_number = 2
../../../_images/sphx_glr_difference_images_002.png

Total running time of the script: ( 0 minutes 17.355 seconds)

Gallery generated by Sphinx-Gallery