Patch fixes for v1.0.1

docs/source/installation.md

@@ -33,7 +33,7 @@ cd shapepipe
 
 where `<VERSION>` is a
 [tagged release](https://github.com/CosmoStat/shapepipe/releases) of ShapePipe
-(e.g. `v1.0.0`). It is recommend to use the
+(e.g. `v1.0.1`). It is recommend to use the
 [latest release](https://github.com/CosmoStat/shapepipe/releases/latest)
 unless you want to reproduce an older set of results.
 

example/cfis/config_Rc.ini

@@ -65,3 +65,12 @@ N_RANDOM = 50000
 
 # N_RANDOM is per square degrees if True
 DENSITY = True
+
+# Output healpix mask if True
+SAVE_MASK_AS_HEALPIX = True
+
+# Healpix mask file base name (used if SAVE_MASK_AS_HEALPIX is True)
+HEALPIX_OUT_FILE_BASE = mask_hp
+
+# Healpix mask nside (used if SAVE_MASK_AS_HEALPIX is True)
+HEALPIX_OUT_NSIDE = 1024

example/cfis/config_make_cat_mccd.ini

@@ -70,7 +70,7 @@ FILE_EXT = .fits, .fits, .sqlite, .fits
 NUMBERING_SCHEME = -000-000
 
 SM_DO_CLASSIFICATION = True
-SM_STAR_STRESH = 0.003
+SM_STAR_THRESH = 0.003
 SM_GAL_THRESH = 0.01
 
 SHAPE_MEASUREMENT_TYPE = ngmix

example/cfis/config_tile_MiViSmVi.ini

@@ -90,7 +90,7 @@ ME_LOG_WCS = $SP_RUN/output/log_exp_headers.sqlite
 
 [VIGNETMAKER_RUNNER_RUN_1]
 
-# Create vignets for tiles weights
+# Create vignets for tile weights
 
 INPUT_DIR = last:sextractor_runner_run_1, last:uncompress_fits_runner
 
@@ -177,5 +177,5 @@ PREFIX =
 # Additional parameters for path and file pattern corresponding to single-exposure
 # run outputs
 ME_IMAGE_DIR = split_exp_runner, split_exp_runner, split_exp_runner, sextractor_runner_run_2
-ME_IMAGE_PATTERN = flag, image, weight, sexcat_background
+ME_IMAGE_PATTERN = flag, image, weight, background
 ME_LOG_WCS = $SP_RUN/output/log_exp_headers.sqlite

shapepipe/modules/mask_package/mask.py

@@ -364,7 +364,7 @@ def make_mask(self):
                         self._config[_type]['cat_path'],
                         size_plus=self._config[_type]['size_plus'],
                         flag_value=self._config[_type]['flag'],
-                        typ=_type,
+                        obj_type=_type,
                     )
 
         if not self._err:
@@ -535,7 +535,13 @@ def mask_border(self, width=100, flag_value=4):
 
         return flag
 
-    def mask_dso(self, cat_path, size_plus=0.1, flag_value=8, typ='Messier'):
+    def mask_dso(
+        self,
+        cat_path,
+        size_plus=0.1,
+        flag_value=8,
+        obj_type='Messier',
+    ):
         """Mask DSO.
 
         Create a circular patch for deep-sky objects (DSOs), e.g.

shapepipe/modules/random_cat_package/__init__.py

@@ -24,8 +24,13 @@
 DENSITY : bool, optional
     Option to interpret the number of random objects per square degree; the
     default is ``False``
-TILE_LIST : str, optional
-    Path to tile IDs for overlap flagging
+SAVE_MASK_AS_HEALPIX : bool
+    Output healpix mask if ``True``
+HEALPIX_OUT_FILE_BASE : str, optional
+    Output halpix mask file base name; used only if SAVE_MASK_AS_HEALPIX is
+    ``True``
+HEALPIX_OUT_NSIDE : int, optional
+    Output healpix mask nside; used only if SAVE_MASK_AS_HEALPIX is ``True``
 
 """
 

shapepipe/modules/random_cat_package/random_cat.py

@@ -72,38 +72,42 @@ def __init__(
         self._w_log = w_log
         self._healpix_options = healpix_options
 
-    def save_as_healpix(self, hdu_mask, header):
+    def save_as_healpix(self, mask, header):
         """Save As Healpix.
 
         Save mask as healpix FITS file.
 
         Parameters
         ----------
-        hdu_mask : class HDUList
-            HDU with 2D pixel mask image
+        mask : numpy.ndarray
+            2D pixel mask image
+
         header : class Header
             Image header with WCS information
 
         """
         if not self._healpix_options:
             return
 
+        # Tranform config entry from str to int
+        nside = int(self._healpix_options['NSIDE'])
+
         mask_1d, footprint = reproject_to_healpix(
-            (hdu_mask, header),
+            (mask, header),
             'galactic',
-            nside=self._healpix_options['OUT_NSIDE']
+            nside=nside,
         )
 
         t = Table()
         t['flux'] = mask_1d
         t.meta['ORDERING'] = 'RING'
         t.meta['COORDSYS'] = 'G'
-        t.meta['NSIDE'] = self._healpix_options['OUT_NSIDE']
+        t.meta['NSIDE'] = nside
         t.meta['INDXSCHM'] = 'IMPLICIT'
 
         output_path = (
-            f'{output_dir}/{self._healpix_options["FILE_BASE"]}-'
-            + f'{file_number_string}.fits'
+            f'{self._output_dir}/{self._healpix_options["FILE_BASE"]}-'
+            + f'{self._file_number_string}.fits'
         )
         t.write(output_path)
 
@@ -135,15 +139,15 @@ def process(self):
         mask = hdu_mask[0].data
 
         # Save mask in healpix format (if option is set)
-        self._save_as_healpix(hdu_mask, header)
+        self.save_as_healpix(mask, header)
 
         # Number of pixels
         n_pix_x = mask.data.shape[0]
         n_pix_y = mask.data.shape[1]
         n_pix = n_pix_x * n_pix_y
 
         # Number of non-masked pixels
-        n_unmasked = len(np.where(hdu_mask[0].data == 0)[0])
+        n_unmasked = len(np.where(mask == 0)[0])
 
         # Compute various areas
 
@@ -209,10 +213,10 @@ def process(self):
 
         # Tile ID
         output_path = (
-            f'{output_dir}/{output_file_pattern}-'
-            + f'{file_number_string}.fits'
+            f'{self._output_dir}/{self._output_file_pattern}-'
+            + f'{self._file_number_string}.fits'
         )
-        file_name = os.path.split(self._output_path)[1]
+        file_name = os.path.split(output_path)[1]
         file_base = os.path.splitext(file_name)[0]
         tile_ID_str = re.split('-', file_base)[1:]
         tile_id = float('.'.join(tile_ID_str))
@@ -230,7 +234,7 @@ def process(self):
 
         # TODO: Add units to header
         output = file_io.FITSCatalogue(
-            self._output_path,
+            output_path,
             open_mode=file_io.BaseCatalogue.OpenMode.ReadWrite
         )
         output.save_as_fits(cat_out, names=column_names)

shapepipe/modules/random_cat_runner.py

@@ -56,7 +56,7 @@ def random_cat_runner(
     )
     if save_mask_as_healpix:
         healpix_options = {}
-        for option_trunc in ['FILE_BASE', 'OUT_NSIDE']:
+        for option_trunc in ['FILE_BASE', 'NSIDE']:
             option = f'HEALPIX_OUT_{option_trunc}'
             healpix_options[option_trunc] = config.get(
                 module_config_sec,