Sample data set overview

An overview of the coordinated sample data set.

import matplotlib.pyplot as plt
import astropy.units as u

import sunpy.timeseries
import as sample_data

On 2011 June 7, various solar instruments observed a spectacular solar eruption from NOAA AR 11226. The event included an M2.5 flare, a filament eruption, a coronal mass ejection, and a global coronal EUV wave (IAU standard: SOL2011-06-07T06:24:00L045C112). This event was spectacular because it features the ejection of a large amount of prominence material, much of which failed to escape and fell back to the solar surface. This event received some press coverage (e.g. National Geographics, Discover Magazine) and the literature contains a number of a papers about it (e.g. Li et al., Inglis et al.)

The following image of the flare is now fairly iconic.

aia_cutout03_map =
fig = plt.figure()
ax = fig.add_subplot(111, projection=aia_cutout03_map)
AIA $193 \; \mathrm{\mathring{A}}$ 2011-06-07 06:45:55

Let’s take a look at the GOES XRS data.

goes = sunpy.timeseries.TimeSeries(sample_data.GOES_XRS_TIMESERIES)
fig = plt.figure()
2011 06 07 sampledata overview


/home/docs/checkouts/ UserWarning: Discarding nonzero nanoseconds in conversion
  return dispatch(args[0].__class__)(*args, **kw)

Next let’s investigate the AIA full disk images that are available. Please note that these images are not at the full AIA resolution.

aia_131_map =
aia_171_map =
aia_211_map =
aia_335_map =
aia_094_map =
aia_1600_map =

fig = plt.figure(figsize=(6, 28))
ax = fig.add_subplot(611, projection=aia_131_map)
aia_131_map.plot(clip_interval=(0.5, 99.9)*u.percent)

ax = fig.add_subplot(612, projection=aia_171_map)
aia_171_map.plot(clip_interval=(0.5, 99.9)*u.percent)

ax = fig.add_subplot(613, projection=aia_211_map)
aia_211_map.plot(clip_interval=(0.5, 99.9)*u.percent)

ax = fig.add_subplot(614, projection=aia_335_map)
aia_335_map.plot(clip_interval=(0.5, 99.9)*u.percent)

ax = fig.add_subplot(615, projection=aia_094_map)
aia_094_map.plot(clip_interval=(0.5, 99.9)*u.percent)

ax = fig.add_subplot(616, projection=aia_1600_map)
aia_1600_map.plot(clip_interval=(0.5, 99.9)*u.percent)

AIA $131 \; \mathrm{\mathring{A}}$ 2011-06-07 06:33:01, AIA $171 \; \mathrm{\mathring{A}}$ 2011-06-07 06:33:02, AIA $211 \; \mathrm{\mathring{A}}$ 2011-06-07 06:33:02, AIA $335 \; \mathrm{\mathring{A}}$ 2011-06-07 06:33:03, AIA $94 \; \mathrm{\mathring{A}}$ 2011-06-07 06:33:05, AIA $1600 \; \mathrm{\mathring{A}}$ 2011-06-07 06:33:05

We also provide a series of AIA cutouts so that you can get a sense of the dynamics of the in-falling material.

aia_cutout01_map =
aia_cutout02_map =
aia_cutout03_map =
aia_cutout04_map =
aia_cutout05_map =

fig = plt.figure(figsize=(6, 28))
ax = fig.add_subplot(511, projection=aia_cutout01_map)

ax = fig.add_subplot(512, projection=aia_cutout02_map)

ax = fig.add_subplot(513, projection=aia_cutout03_map)

ax = fig.add_subplot(514, projection=aia_cutout04_map)

ax = fig.add_subplot(515, projection=aia_cutout05_map)

AIA $193 \; \mathrm{\mathring{A}}$ 2011-06-07 06:33:07, AIA $193 \; \mathrm{\mathring{A}}$ 2011-06-07 06:39:31, AIA $193 \; \mathrm{\mathring{A}}$ 2011-06-07 06:45:55, AIA $193 \; \mathrm{\mathring{A}}$ 2011-06-07 06:52:19, AIA $193 \; \mathrm{\mathring{A}}$ 2011-06-07 06:58:43

There are a number of other data sources available as well, such as SWAP.

swap_map =
fig = plt.figure()
SWAP $174 \; \mathrm{\mathring{A}}$ 2011-06-07 06:33:29


/home/docs/checkouts/ RuntimeWarning: cdelt will be ignored since cd is present
  w2.wcs.cdelt = u.Quantity(self.scale)

And also RHESSI.

rhessi_map =
fig = plt.figure()
RHESSI $[25,~50] \; \mathrm{keV}$ 2011-06-07 06:33:00


/home/docs/checkouts/ FITSFixedWarning: CROTACN1= 0.0 / X Position of Center of Rotation (arcsec)
keyword looks very much like CROTAn but isn't.
  wcsprm = _wcs.Wcsprm(header=header_bytes, key=key,
/home/docs/checkouts/ FITSFixedWarning: CROTACN2= 0.0 / Y Position of Center of Rotation (arcsec)
keyword looks very much like CROTAn but isn't.
  wcsprm = _wcs.Wcsprm(header=header_bytes, key=key,
/home/docs/checkouts/ FITSFixedWarning: CROTA = 0.0 / Rotation Angle (clockwise from N)
keyword looks very much like CROTAn but isn't.
  wcsprm = _wcs.Wcsprm(header=header_bytes, key=key,
/home/docs/checkouts/ SunpyUserWarning: Missing metadata for observer: assuming Earth-based observer.The following sets of keys were checked:
('hgln_obs', 'hglt_obs', 'dsun_obs')
('crln_obs', 'crlt_obs', 'dsun_obs')
  warnings.warn(warning_message, SunpyUserWarning)

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

Gallery generated by Sphinx-Gallery