Closed ebuchlin closed 7 months ago
A proposed interface, as part of the FileMetadata
class (corresponding to 1 catalog row):
def plot_fov(self, background_map=None, background_wcs=None, differential_rotation=True):
"""
Plot FOV and a background or blank map
Parameters
-----------------
background_map: sunpy.map.Map
A background map
background_wcs: astropy.wcs.WCS
WCS object to build a blank background map
differential_rotation: bool
Shift FOV edges taking into account differential rotation between background map reference time and SPICE observation
Either background_map or background_wcs must be provided.
"""
Open questions:
Implemented, see PR #43. Remaining linked issues: #39, #40, #41.
The idea is that:
Catalog
object,Catalog.find_files()
...)sunpy.Map
provided by the user, or on a blank map built from WCS headers (e.g. the result ofsunpy.map.make_fitswcs_header()
) of from simple frame parameters.The background/blank map could be in any relevant coordinate frame (including Helioprojective from the same or different viewpoint/time, Carrington, Stonyhurst). The projection should take into account the possibly different point of view between Solar Orbiter and the background/blank map, as well as the time difference between each observation and the reference time for the background/blank map (change of Solar Orbiter's position on its orbit, solar (differential) rotation).