[NEW FEATURE] Read tiles catalogue instead of SExtractor detection run

[martinkilbinger]

Desired additional functionality: Reading in a reference (SExtractor) catalogue for each tile to create galaxy candidate catalogue. This catalogue would be the one created by Stephen Gwyn.

This feature should be an additional option to running SExtractor for detection from ShapePipe.

Feature could be implemented as follows:

• in sextractor_runner, to replace call of sextractor, and only use post-processing features of that module.
• to replace sextractor_runner, if post-processing is not required

Potential subsequent changes to ShapePipe are:

• Using segmentation map instead of mask from postage stamp. This would avoid to request postage stamps in the reference catalogue, used in vignet_maker.

[martinkilbinger]

Our current output parameter file is https://github.com/CosmoStat/shapepipe/blob/develop/example/cfis/default.param

[martinkilbinger]

Stephen's current output catalogue parameters are:

#   1 NUMBER                 Running object number                                     
#   2 X_IMAGE                Object position along x                                    [pixel]
#   3 Y_IMAGE                Object position along y                                    [pixel]
#   4 ALPHA_J2000            Right ascension of barycenter (J2000)                      [deg]
#   5 DELTA_J2000            Declination of barycenter (J2000)                          [deg]
#   6 MAG_AUTO               Kron-like elliptical aperture magnitude                    [mag]
#   7 MAGERR_AUTO            RMS error for AUTO magnitude                               [mag]
#   8 MAG_BEST               Best of MAG_AUTO and MAG_ISOCOR                            [mag]
#   9 MAGERR_BEST            RMS error for MAG_BEST                                     [mag]
#  10 MAG_APER               Fixed aperture magnitude vector                            [mag]
#  11 MAGERR_APER            RMS error vector for fixed aperture mag.                   [mag]
#  12 A_WORLD                Profile RMS along major axis (world units)                 [deg]
#  13 ERRA_WORLD             World RMS position error along major axis                  [deg]
#  14 B_WORLD                Profile RMS along minor axis (world units)                 [deg]
#  15 ERRB_WORLD             World RMS position error along minor axis                  [deg]
#  16 THETA_J2000            Position angle (east of north) (J2000)                     [deg]
#  17 ERRTHETA_J2000         J2000 error ellipse pos. angle (east of north)             [deg]
#  18 ISOAREA_IMAGE          Isophotal area above Analysis threshold                    [pixel**2]
#  19 MU_MAX                 Peak surface brightness above background                   [mag * arcsec**(-2)]
#  20 FLUX_RADIUS            Fraction-of-light radii                                    [pixel]

[martinkilbinger]

Not sure we need the following parameters:

• EXT_NUMBER
• BACKGROUND #Background at centroid position [count]
• THRESHOLD #Detection threshold above background [count]
• X_WORLD, Y_WORLD
• X2_IMAGE, Y2_IMAGE, XY_IMAGE + errors
• MU_THRESHOLD #Analysis threshold above background [mag * arcsec**(-2)]
• FWHM_WORLD #FWHM assuming a gaussian core [deg]
• ELONGATION #A_IMAGE/B_IMAGE
• ELLIPTICITY #1 - B_IMAGE/A_IMAGE

To debate:

• VIGNET(51,51) #Pixel data around detection [count]

That leaves the following additional parameters:

FLUXERR_AUTO             #RMS error for AUTO flux                                   [count]
FLUX_WIN                 #Gaussian-weighted flux                                    [count]
FLUXERR_WIN              #RMS error for WIN flux                                    [count]
MAG_WIN                  #Gaussian-weighted magnitude                               [mag]
MAGERR_WIN               #RMS error for MAG_WIN                                     [mag]
FLUX_APER(1)                                                                    
FLUXERR_APER(1)                                                                 

SNR_WIN                  #Gaussian-weighted SNR                                 

XWIN_IMAGE               #Windowed position estimate along x                        [pixel]
YWIN_IMAGE               #Windowed position estimate along y                        [pixel]

XWIN_WORLD               #Windowed position along world x axis                      [deg]
YWIN_WORLD               #Windowed position along world y axis                      [deg]

X2WIN_IMAGE              #Windowed variance along x                                 [pixel**2]
Y2WIN_IMAGE              #Windowed variance along y                                 [pixel**2]
XYWIN_IMAGE              #Windowed covariance between x and y                       [pixel**2]
ERRX2WIN_IMAGE           #Variance of windowed pos along x                          [pixel**2]
ERRY2WIN_IMAGE           #Variance of windowed pos along y                          [pixel**2]
ERRXYWIN_IMAGE           #Covariance of windowed pos between x and y                [pixel**2]

MU_MAX                   #Peak surface brightness above background                  [mag * arcsec**(-2)]

FLAGS_WIN                #Flags for WINdowed parameters                         

IMAFLAGS_ISO             #FLAG-image flags OR'ed over the iso. profile ```

[lbaumo]

this could be useful to convert ascii tables to fits-ldac:astromatic_wrapper from what I can see, this just deals with the OBJECTS table. We also want to be able to access the FIELDS table