# import astropy.units as u
#
# from astropy.coordinates import SkyCoord
# from parfive import Downloader, SessionConfig
# from aiohttp import ClientTimeout
# from sunpy.net import Fido, attrs as a
# from sunpy.map import Map
# from sunpy.coordinates.frames import HeliographicStonyhurst
#
# from stixpy.visualisation.map_reprojection import reproject_map, plot_map_reproj, get_solo_position
#
# # Search and download map using FIDO
# aia = (a.Instrument.aia &
#     a.Sample(24 * u.hour) &
#     a.Time('2021-04-17', '2021-04-18'))
# wave = a.Wavelength(19.3 * u.nm, 19.3 * u.nm)
# query = Fido.search(wave, aia)
#
# dl = Downloader(config=SessionConfig(ClientTimeout(total=120)))
# results = Fido.fetch(query[0][0], downloader=dl)
#
# # Create map and resample to speed up calculations
# map = Map(results)
# map = map.resample([512, 512]*u.pix)
#
# # Get SOLO position observer
# # observer = get_solo_position(map) # broken on RTD
# observer = SkyCoord(-97.01771373*u.deg, 21.20931737*u.deg, 1.25689226e+08*u.km,
#                     frame=HeliographicStonyhurst, obstime='2021-04-17 00:00:04.840000')
#
# # Reproject Map
# reprojected_map = reproject_map(map, observer, out_shape=(768, 768))
#
# # Run plotting function
# plot_map_reproj(map, reprojected_map)