# import astropy.units as u
# from astropy.coordinates import SkyCoord
# 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
# query = Fido.search(a.Time('2020/06/12 13:20:00', '2020/06/12 13:40:00'),
#                     a.Instrument.hmi, a.Physobs.los_magnetic_field)
# result = Fido.fetch(query[0][0])
#
# # Create map and resample to speed up calculations
# map = Map(result)
# map = map.resample([512, 512] * u.pix)
#
# # Set SOLO as observer
# # observer = get_solo_position(map) # broke on RTD
# observer = SkyCoord(56.18727061*u.deg, 3.5358653*u.deg, 77369481.8542484*u.km,
#                     frame=HeliographicStonyhurst, obstime='2020-06-12 13:20:08.300000')
#
# # Reproject Map
# reprojected_map = reproject_map(map, observer, out_shape=(1024, 1024))
#
# # Run plotting function
# plot_map_reproj(map, reprojected_map)