# HMI images have been downloaded for you and are already in the folder (using the online search form https://vso.nascom.nasa.gov/cgi-bin/search).
# Let's add the fits in our folder to the database.
db.add_from_dir("./data/", ignore_already_added=True,time_string_parse_format="%d/%m/%Y")
# Let's see what's now in our database:
for database_entry in db:
if database_entry.observation_time_start is None and database_entry.observation_time_end is None:
db.remove(database_entry)
print(display_entries(db,['id', 'observation_time_start','instrument', 'wavemin']))
# We can now search the database. Here we are searching for images in wavelengths between 1-2 nm or 10-20 Angstrom.
# We have also chosen to sort our results by wavlength!
print(display_entries(db.search(a.Wavelength(1.0*u.nm, 2.0*u.nm)),['id', 'observation_time_start', 'instrument', 'wavemin'], sort=True))
# Another way to fetch data is to use:
# entries = db.fetch(a.Time('2019/05/06 12:00', '2019/05/06 12:01'),a.Instrument('aia'),a.vso.Sample(2*u.minute),a.Wavelength(94*u.angstrom))
# This automatically adds the observations to our database and downloads files checking for duplicates.
# Now we have our images let's plot them!
AIA = db.search(a.Wavelength(.1*u.nm, 60.0*u.nm))
AIAplotDic = {}
for obs in AIA:
AIAplotDic[int(obs.wavemin*10)] = [obs.path,"sdoaia"+str(int(obs.wavemin*10))]
# fig = plt.figure(num=3, figsize=(len(AIA)*2,2))
# count = 1
# for im in sorted(AIAplotDic.keys()):
