Add img_size guide parameter to control guide star image readout size

proseco/catalog.py

@@ -62,6 +62,8 @@ def get_aca_catalog(obsid=0, **kwargs):
     :param exclude_ids_fid: list of fiducial lights to exclude by index
     :param include_ids_guide: list of AGASC IDs of stars to include in guide catalog
     :param exclude_ids_guide: list of AGASC IDs of stars to exclude from guide catalog
+    :param img_size_guide: readout window size for guide stars (6, 8, or ``None``).
+                           For default value ``None`` use 8 for no fids, 6 for fids.
     :param optimize: optimize star catalog after initial selection (default=True)
     :param verbose: provide extra logging info (mostly calc_p_safe) (default=False)
     :param print_log: print the run log to stdout (default=False)
@@ -102,6 +104,7 @@ def _get_aca_catalog(**kwargs):
     assembling the merged aca catalog.
     """
     raise_exc = kwargs.pop('raise_exc')
+    img_size_guide = kwargs.pop('img_size_guide', None)
 
     aca = ACATable()
     aca.call_args = kwargs.copy()
@@ -141,7 +144,8 @@ def _get_aca_catalog(**kwargs):
     aca.acqs.fid_set = aca.fids['id']
 
     aca.log('Starting get_guide_catalog')
-    aca.guides = get_guide_catalog(stars=aca.acqs.stars, fids=aca.fids, **kwargs)
+    aca.guides = get_guide_catalog(stars=aca.acqs.stars, fids=aca.fids,
+                                   img_size=img_size_guide, **kwargs)
 
     # Make a merged starcheck-like catalog.  Catch any errors at this point to avoid
     # impacting operational work (call from Matlab).
@@ -439,7 +443,14 @@ def merge_cats(fids=None, guides=None, acqs=None):
         fids['p_acq'] = 0
         fids['sz'] = '8x8'
 
-    guide_size = '8x8' if len(fids) == 0 else '6x6'
+    # Set guide star image readout size, either from explicit user input or
+    # using the default rule that OR's get 6x6 and ER's (no fids) get 8x8.
+    img_size = guides.img_size
+    if img_size is None:
+        guide_size = '8x8' if len(fids) == 0 else '6x6'
+    else:
+        guide_size = f'{img_size}x{img_size}'
+
     if len(guides) > 0:
         guides['slot'] = 0  # Filled in later
         guides['type'] = 'GUI'

proseco/guide.py

@@ -67,6 +67,13 @@ def get_guide_catalog(obsid=0, **kwargs):
     return guides
 
 
+class ImgSizeMetaAttribute(MetaAttribute):
+    def __set__(self, instance, value):
+        if value not in (4, 6, 8, None):
+            raise ValueError('img_size must be 4, 6, 8, or None')
+        instance.meta[self.name] = value
+
+
 class GuideTable(ACACatalogTable):
     # Define base set of allowed keyword args to __init__. Subsequent MetaAttribute
     # or AliasAttribute properties will add to this.
@@ -88,6 +95,7 @@ class GuideTable(ACACatalogTable):
 
     cand_guides = MetaAttribute(is_kwarg=False)
     reject_info = MetaAttribute(default=[], is_kwarg=False)
+    img_size = ImgSizeMetaAttribute()
 
     def reject(self, reject):
         """

proseco/tests/test_catalog.py

@@ -766,3 +766,40 @@ def test_dark_date_warning():
 
     assert len(warns) == 1
     assert 'Unexpected dark_date: dark_id nearest dark_date' in str(warns[0].message)
+
+
+def test_img_size_guide():
+    """
+    Test img_size_guide for setting guide star readout image size
+    """
+    dark = DARK40.copy()
+    stars = StarsTable.empty()
+    stars.add_fake_constellation(mag=8.0, n_stars=8)
+
+    # Confirm that for an inferred ER, boxes are 8x8
+    aca = get_aca_catalog(**mod_std_info(stars=stars, dark=dark, n_fid=0))
+    assert np.all(aca.guides['sz'] == '8x8')
+    assert aca.guides.img_size is None
+
+    # Confirm that for an inferred OR, boxes are 6x6
+    aca = get_aca_catalog(**mod_std_info(stars=stars, dark=dark, n_fid=3))
+    assert np.all(aca.guides['sz'] == '6x6')
+    assert aca.guides.img_size is None
+
+    # Confirm that for explicit img_size_guide of 8 boxes are '8x8'
+    aca = get_aca_catalog(**mod_std_info(stars=stars, dark=dark, n_fid=3, img_size_guide=8))
+    assert np.all(aca.guides['sz'] == '8x8')
+    assert aca.guides.img_size == 8
+
+    # Confirm that for explicit img_size_guide of 6 boxes are '6x6'
+    aca = get_aca_catalog(**mod_std_info(stars=stars, dark=dark, n_fid=0, img_size_guide=6))
+    assert np.all(aca.guides['sz'] == '6x6')
+    assert aca.guides.img_size == 6
+
+    # Confirm that for explicit img_size_guide of 6 boxes are '6x6'
+    aca = get_aca_catalog(**mod_std_info(stars=stars, dark=dark, n_fid=0, img_size_guide=4))
+    assert np.all(aca.guides['sz'] == '4x4')
+    assert aca.guides.img_size == 4
+
+    with pytest.raises(ValueError, match='img_size must be 4, 6, 8, or None'):
+        get_aca_catalog(**mod_std_info(stars=stars, dark=dark, img_size_guide=3))