Implements SOLOHi imager

Author: AthKouloumvakos

PyThea/PyThea_app.py

@@ -85,7 +85,7 @@ def footer_text():
             More imaging data have been added:
             - SDO/AIA images from 211A channel.
             - PSP/WISPR inner and outer telescope images.
-            - SolO/EUI and METIS images.
+            - SolO/EUI, METIS, and HI(tiles 1&2) images.
             ''', icon='ℹ️')
     st.warning('''
            **NOTE: From PyThea >0.8.1 the JSON fitting files will be slightly different from the old ones.**

PyThea/config/selected_imagers.py

@@ -80,3 +80,10 @@
 imager_dict['METIS'] = {'fido': (a.Instrument('METIS'), a.soar.Product('METIS-VL-TB'), a.Level(2)),
                         'process': {'dimensions': (1024*u.pixel, 1024*u.pixel), 'superpixel': 2},
                         'source': 'SOLO', 'instrument': 'METIS', 'detector': 'VLD', 'wavelength': 'TB'}
+
+for tile in range(1, 5):
+    z = 'T' if tile in [1, 2] else 'G'
+    imager_dict[f'SOLOHI-T{tile}'] = {'fido': (a.Instrument('SOLOHI'),
+                                               a.soar.Product(f'SOLOHI-{tile}F{z}'), a.Level(2)),
+                                      'process': {'superpixel': 2},
+                                      'source': 'SOLO', 'instrument': 'SOLOHI', 'detector': f'T{tile}'}

PyThea/sunpy_dev/map/maputils.py

@@ -317,6 +317,9 @@ def maps_clims(images):
             return [14., 14.]
         elif images[i].detector == 'Inner':
             return [13., 13.]
-    if images[i].instrument == 'Metis':
+    elif images[i].instrument == 'Metis':
         return [12.80, 13.00]
+    elif images[i].instrument.startswith('SoloHI'):
+        return [13.95, 14.05]
+
     return [np.nanquantile(images[i].data, 0.20)-10, np.nanquantile(images[i].data, 0.80)+10]

PyThea/utils.py

@@ -128,10 +128,7 @@ def make_figure(map, cmap='Greys_r', clim=[-20, 20], clip_model=True, **kwargs):
     if median_filter_value != 1:
         map = sunpy.map.Map(median_filter(map.data, size=int(median_filter_value)), map.meta)
 
-    if map.instrument == 'WISPR':
-        clim = [-10**-clim[0], 10**-clim[1]]
-
-    if map.instrument == 'Metis':
+    if map.instrument in ['WISPR', 'Metis'] or map.instrument.startswith('SoloHI'):
         clim = [-10**-clim[0], 10**-clim[1]]
 
     if cmap == 'default':
@@ -155,9 +152,19 @@ def make_figure(map, cmap='Greys_r', clim=[-20, 20], clip_model=True, **kwargs):
     if cref.Ty > 0:
         axis.invert_yaxis()
 
-    axis.set_title(re.sub(r'\d{4}-\d{2}-\d{2} \d{2}:\d{2}:\d{2}',
-                          ' $T_{AGV}:$' + parse_time(map.date_average).strftime('%Y-%m-%d %H:%M:%S'),
-                          map.latex_name.replace('VLD', 'METIS-VDL')),
+    if map.instrument == 'SoloHI':
+        title = 'SoloHI' + f' Tile-{map.detector}' ' $T_{AGV}:$' + parse_time(map.date_average).strftime('%Y-%m-%d %H:%M:%S')
+    elif map.instrument == 'Metis':
+        title = re.sub(r'\d{4}-\d{2}-\d{2} \d{2}:\d{2}:\d{2}',
+                       ' $T_{AGV}:$' + parse_time(map.date_average).strftime('%Y-%m-%d %H:%M:%S'),
+                       map.latex_name.replace('VLD', 'METIS-VDL'))
+    else:
+        title = re.sub(r'\d{4}-\d{2}-\d{2} \d{2}:\d{2}:\d{2}',
+                       ' $T_{AGV}:$' + parse_time(map.date_average).strftime('%Y-%m-%d %H:%M:%S'),
+                       map.latex_name)
+    print(map.instrument)
+    print(map.latex_name)
+    axis.set_title(title,
                    fontsize=10, pad=8)
 
     return fig, axis