Tutorial:How to contribute to the flare wiki pages
import os import io import astropy.units as u from matplotlib import pyplot as plt from sunpy.net import Fido, attrs as a from sunpy.map import Map from matplotlib.backends.backend_agg import FigureCanvasAgg as FigureCanvas from matplotlib.figure import Figure
def plot(folder,flare_id , start_utc, end_utc, peak_utc, wavelen=131, overwrite=False):
sdo_query = Fido.search(a.Time(utc_start, utc_end), a.Instrument('AIA'),
a.Wavelength(wavelen* u.angstrom))
sdo_res = Fido.fetch(sdo_query[0], progress=False, path=folder)
if not sdo_res:
print('AIA data not available')
return None
if sdo_res:
sdo = Map(sdo_res[0])
fig = plt.figure(figsize=(6, 6), dpi=100)
ax = fig.add_subplot(projection=sdo)
sdo.plot(clip_interval=[1, 100] * u.percent, axes=ax)
fname=os.path.join(folder, f'AIA_{wavelen}_{flare_id}.png')
plt.savefig(fname, dpi=100)
return fname
return None
if __name__=='__main__':
#for test only
plot('.',0,'2021-04-01T00:00:00', '2021-04-01T00:10:00', '2021-04-01T00:05:00' )