Tutorial:How to contribute to the flare wiki pages

From stix
Revision as of 11:37, 14 May 2021 by Hualin (talk | contribs) (Created page with "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.ba...")
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

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' )