IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changeset 15555


Ignore:
Timestamp:
Nov 9, 2007, 2:44:16 PM (19 years ago)
Author:
eugene
Message:

adding new residuals option

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/ippconfig/recipes/psphot.config

    r15079 r15555  
    11
    22# these options turn on different inputs and/or outputs
    3 SAVE.OUTPUT             BOOL    TRUE
    4 SAVE.BACKMDL            BOOL    FALSE
    5 SAVE.BACKMDL.STDEV      BOOL    FALSE
    6 SAVE.BACKGND            BOOL    FALSE
    7 SAVE.BACKSUB            BOOL    FALSE
    8 SAVE.RESID              BOOL    FALSE
    9 SAVE.PSF                BOOL    FALSE
    10 LOAD.PSF                BOOL    FALSE
    11 SAVE.PLOTS              BOOL    FALSE
     3SAVE.OUTPUT                         BOOL  TRUE
     4SAVE.BACKMDL                        BOOL  FALSE
     5SAVE.BACKMDL.STDEV                  BOOL  FALSE
     6SAVE.BACKGND                        BOOL  FALSE
     7SAVE.BACKSUB                        BOOL  FALSE
     8SAVE.RESID                          BOOL  FALSE
     9SAVE.PSF                            BOOL  FALSE
     10LOAD.PSF                            BOOL  FALSE
     11SAVE.PLOTS                          BOOL  FALSE
    1212
    1313# the zero point is used to set a basic scale for DVO
    1414# XXX it may not currently be read : double check this (EAM)
    15 ZERO_POINT          F32  25.000          # zero point used by DVO
    16 ZERO_PT             F32  25.000          # zero point used by DVO
    17 
    18 MASKVAL             STR  SAT,BAD,BLANK   # Mask these types of pixels
    19 
    20 OUTPUT.FORMAT       STR  SMPDATA
     15ZERO_POINT                          F32   25.000          # zero point used by DVO
     16ZERO_PT                             F32   25.000          # zero point used by DVO
     17
     18MASKVAL                             STR   SAT,BAD,BLANK   # Mask these types of pixels
     19
     20OUTPUT.FORMAT                       STR   SMPDATA
    2121
    2222# these parameter govern how the background is measured
    23 BACKGROUND.XBIN     S32   128            # size of background superpixels
    24 BACKGROUND.YBIN     S32   128            # size of background superpixels
    25 IMSTATS_NPIX        S32  10000          # number of pixels to use for sky estimate boxes:
    26 
    27 SKY_BIAS            F32  0.0             # offset applied to measured sky (FOR TESTING)
    28 SKY_FIT_ORDER       S32  0
    29 SKY_FIT_LINEAR      BOOL FALSE
    30 SKY_STAT            STR  FITTED_MEAN_V4  # statistic used to measure background
    31 SKY_CLIP_SIGMA      F32  2.0             # statistic used to measure background
    32 SKY_SIG             F32  1.0             # optional sky error for
     23BACKGROUND.XBIN                     S32   128             # size of background superpixels
     24BACKGROUND.YBIN                     S32   128             # size of background superpixels
     25IMSTATS_NPIX                        S32   10000          # number of pixels to use for sky estimate boxes:
     26
     27SKY_BIAS                            F32   0.0             # offset applied to measured sky (FOR TESTING)
     28SKY_FIT_ORDER                       S32   0
     29SKY_FIT_LINEAR                      BOOL FALSE
     30SKY_STAT                            STR   FITTED_MEAN_V4  # statistic used to measure background
     31SKY_CLIP_SIGMA                      F32   2.0             # statistic used to measure background
     32SKY_SIG                             F32   1.0             # optional sky error for
    3333
    3434# allowed values for SKY_STAT:
     
    3636
    3737# masking parameters (XXX EAM : rework this to use psRegion like ANALYSIS_REGION)
    38 XMIN                F32   0              # minimum valid x-coord
    39 XMAX                F32   0              # maximum valid x-coord
    40 YMIN                F32   0              # minimum valid y-coord
    41 YMAX                F32   0              # maximum valid y-coord
     38XMIN                                F32   0              # minimum valid x-coord
     39XMAX                                F32   0              # maximum valid x-coord
     40YMIN                                F32   0              # minimum valid y-coord
     41YMAX                                F32   0              # maximum valid y-coord
    4242
    4343# peak finding
    44 PEAKS_SMOOTH_SIGMA  F32  1.0            # smoothing kernel sigma in pixels
    45 PEAKS_SMOOTH_NSIGMA F32  2.0            # smoothing kernel width in sigmas
    46 PEAKS_NSIGMA_LIMIT  F32  25.0            # peak significance threshold
    47 PEAKS_NSIGMA_LIMIT_2 F32  5.0            # peak significance threshold
    48 PEAKS_NMAX          S32   0              # on first pass, only keep NMAX peaks (0 == all)
     44PEAKS_SMOOTH_SIGMA                  F32   1.0            # smoothing kernel sigma in pixels
     45PEAKS_SMOOTH_NSIGMA                 F32   2.0            # smoothing kernel width in sigmas
     46PEAKS_NSIGMA_LIMIT                  F32   25.0            # peak significance threshold
     47PEAKS_NSIGMA_LIMIT_2                F32   5.0            # peak significance threshold
     48PEAKS_NMAX                          S32   0               # on first pass, only keep NMAX peaks (0 == all)
    4949
    5050# parameters to control the selection of the peak in the Sx,Sy plane
    51 MOMENTS_SCALE       F32   0.1
    52 MOMENTS_SN_MIN      F32   100.0           # min S/N to measure moments
    53 MOMENTS_SX_MAX      F32   3.0
    54 MOMENTS_SY_MAX      F32   3.0
    55 MOMENTS_AR_MAX      F32   1.5            # maximum axial ratio: 1 / AR < (sx / sy) < AR
     51MOMENTS_SCALE                       F32   0.1       
     52MOMENTS_SN_MIN                      F32   100.0           # min S/N to measure moments
     53MOMENTS_SX_MAX                      F32   3.0
     54MOMENTS_SY_MAX                      F32   3.0
     55MOMENTS_AR_MAX                      F32   1.5            # maximum axial ratio: 1 / AR < (sx / sy) < AR
    5656
    5757# basic object statistics
    58 SKY_INNER_RADIUS    F32  15              # square annulus for local sky measurement
    59 SKY_OUTER_RADIUS    F32  25              # square annulus for local sky measurement
    60 PSF_MOMENTS_RADIUS  F32  3               # calculate initial source moments with this radius
    61 PSF_SN_LIM          F32  50              # minimum S/N for stars used for PSF model
    62 PSF_MAX_NSTARS      S32  200             # limit number of stars used for PSF model
    63 PSF_CLUMP_NSIGMA    F32  1.5             # region of Sx,Sy plane to use for selecting PSF stars
    64 PSF_MIN_DS          F32  0.01
    65 PSF_PARAM_WEIGHTS   BOOL FALSE
     58SKY_INNER_RADIUS                    F32   15              # square annulus for local sky measurement
     59SKY_OUTER_RADIUS                    F32   25              # square annulus for local sky measurement
     60PSF_MOMENTS_RADIUS                  F32   3               # calculate initial source moments with this radius
     61PSF_SN_LIM                          F32   50              # minimum S/N for stars used for PSF model
     62PSF_MAX_NSTARS                      S32   200             # limit number of stars used for PSF model
     63PSF_CLUMP_NSIGMA                    F32   1.5             # region of Sx,Sy plane to use for selecting PSF stars
     64PSF_MIN_DS                          F32   0.01
     65PSF_PARAM_WEIGHTS                   BOOL FALSE
    6666
    6767# PSF model parameters : choose the PSF model
    6868# list as many PSF_MODEL options as desired
    69 # PSF_MODEL           MULTI
    70 PSF_MODEL           STR  PS_MODEL_GAUSS
    71 # PSF_MODEL           STR  PS_MODEL_PGAUSS
    72 # PSF_MODEL           STR  PS_MODEL_QGAUSS
    73 # PSF_MODEL           STR  PS_MODEL_TGAUSS  ## not well tested, not very successful
     69# PSF_MODEL                         MULTI
     70PSF_MODEL                           STR   PS_MODEL_GAUSS
     71# PSF_MODEL                         STR   PS_MODEL_PGAUSS
     72# PSF_MODEL                         STR   PS_MODEL_QGAUSS
     73# PSF_MODEL                         STR   PS_MODEL_TGAUSS # not well tested, not very successful
    7474
    7575# PSF.TREND.MASK must be a 2D polynomial
    7676# the specified values are ignored but define the active components of the polynomial
    77 PSF.TREND.MASK  METADATA 
    78    NORDER_X         S32       0                # number of x orders
    79    NORDER_Y         S32       0                # number of y orders
    80    VAL_X00_Y00      F64       1                # polynomial coefficient
    81    NELEMENTS        S32       1                # number of unmasked components
     77PSF.TREND.MASK                      METADATA 
     78   NORDER_X                         S32   0               # number of x orders
     79   NORDER_Y                         S32   0               # number of y orders
     80   VAL_X00_Y00                      F64   1               # polynomial coefficient
     81   NELEMENTS                        S32   1               # number of unmasked components
    8282END  # folder for 2D polynomial
    8383
    84 PSF.TREND.MODE STR POLY_ORD
    85 PSF.TREND.NX   S32 0
    86 PSF.TREND.NY   S32 0
    87 
    88 PSF_FIT_RADIUS      F32  15.0            # fitting radius for test PSF model
    89 PSF_REF_RADIUS      F32  25.0            # aperture magnitudes are scaled via
    90                                         # curve-of-growth to this radius
     84PSF.TREND.MODE                      STR POLY_ORD         
     85PSF.TREND.NX                        S32   0
     86PSF.TREND.NY                        S32  0
     87
     88PSF_FIT_RADIUS                      F32   15.0            # fitting radius for test PSF model
     89PSF_REF_RADIUS                      F32   25.0            # aperture magnitudes are scaled via
     90                                        # curve-of-growth to this radius
    9191# PSF-like source model parameters
    92 PSF_FIT_NSIGMA       F32  1.0            # significance for pixel included in fit
    93 PSF_FIT_PADDING      F32  2.0            # extra annulus to use for fit
    94 PSF_SHAPE_NSIGMA     F32  3.0            # max significance for shape variation
    95 PSF_MIN_SN           F32  2.0            # reject objects below this significance
    96 PSF_MAX_CHI          F32  50.0           # reject objects worse that this
    97 FULL_FIT_SN_LIM      F32  50.0
    98 
    99 PSF.RESIDUALS        BOOL TRUE
    100 PSF.RESIDUALS.XBIN   S32   1
    101 PSF.RESIDUALS.YBIN   S32   1
    102 PSF.RESIDUALS.NSIGMA F32   3.0 
    103 PSF.RESIDUALS.INTERPOLATION STR BILINEAR
    104 PSF.RESIDUALS.STATISTIC     STR ROBUST_MEDIAN
    105 PSF.RESIDUALS.SPATIAL_ORDER S32 0
     92PSF_FIT_NSIGMA                      F32   1.0             # significance for pixel included in fit
     93PSF_FIT_PADDING                     F32   2.0             # extra annulus to use for fit
     94PSF_SHAPE_NSIGMA                    F32   3.0             # max significance for shape variation
     95PSF_MIN_SN                          F32   2.0             # reject objects below this significance
     96PSF_MAX_CHI                         F32   50.0            # reject objects worse that this
     97FULL_FIT_SN_LIM                     F32   50.0
     98
     99# the RESIDUALS are pixelized psf residual tables
     100PSF.RESIDUALS                       BOOL  TRUE            # generate the residuals?
     101PSF.RESIDUALS.XBIN                  S32   1               # Nx(residual) = Nx(input)*XBIN
     102PSF.RESIDUALS.YBIN                  S32   1               # Ny(residual) = Ny(input)*YBIN
     103PSF.RESIDUALS.NSIGMA                F32   3.0             # clip input stack of NSIGMA outliers
     104PSF.RESIDUALS.INTERPOLATION         STR   BILINEAR        # interpolation to use when reconstructing residual
     105PSF.RESIDUALS.STATISTIC             STR   ROBUST_MEDIAN   # statistic to use for generating the residual
     106PSF.RESIDUALS.SPATIAL_ORDER         S32   0               # fit spatial variations of the residuals at this order (0,1)
     107PSF.RESIDUALS.PIX.SN                F32   0.0             # keep this pixel if residual is more significant than this
    106108 
    107109# EXTended source model parameters
    108 EXT_MODEL            STR  PS_MODEL_PGAUSS
    109 EXT_MIN_SN           F32  50.0           # fit galaxies above this S/N limit
    110 EXT_FIT_NSIGMA       F32  1              # significance for pixel included in fit
    111 EXT_FIT_PADDING      F32  5              # extra annulus to use for fit
    112 EXT_MOMENTS_RADIUS   F32  9
     110EXT_MODEL                           STR   PS_MODEL_PGAUSS
     111EXT_MIN_SN                          F32   50.0            # fit galaxies above this S/N limit
     112EXT_FIT_NSIGMA                      F32   1              # significance for pixel included in fit
     113EXT_FIT_PADDING                     F32   5               # extra annulus to use for fit
     114EXT_MOMENTS_RADIUS                  F32   9
    113115
    114116# Extended source fit parameters
    115 EXTENDED_SOURCE_FITS BOOL FALSE
    116 EXTENDED_SOURCE_SN_LIM F32 20.0
    117 EXTENDED_SOURCE_PSF_CONVOLVED BOOL FALSE
    118 
    119 FITMODE              STR  BLEND
    120 DEBLEND_PEAK_FRACTION   F32 0.1
    121 DEBLEND_SKY_NSIGMA      F32 10.0
    122 
    123 # APTREND            STR  NONE, CONSTANT, SKYBIAS, SKYSAT, XY_LIN, SKY_XY_LIN, SKYSAT_XY_LIN, ALL
    124 APTREND              STR  CONSTANT
    125 AP_MIN_SN            F32  25.0
    126 APTREND.NSTAR.MIN    S32  5
    127 APTREND.ORDER.MAX    S32  5
     117EXTENDED_SOURCE_FITS                BOOL FALSE
     118EXTENDED_SOURCE_SN_LIM              F32  20.0
     119EXTENDED_SOURCE_PSF_CONVOLVED       BOOL FALSE
     120
     121FITMODE                             STR   BLEND 
     122DEBLEND_PEAK_FRACTION               F32  0.1
     123DEBLEND_SKY_NSIGMA                  F32  10.0
     124
     125# APTREND                           STR   NONE, CONSTANT, SKYBIAS, SKYSAT, XY_LIN, SKY_XY_LIN, SKYSAT_XY_LIN, ALL
     126APTREND                             STR   CONSTANT
     127AP_MIN_SN                           F32   25.0
     128APTREND.NSTAR.MIN                   S32   5
     129APTREND.ORDER.MAX                   S32   5
    128130
    129131# test options
    130132# BREAK_POINT may be one of (PEAKS, MOMENTS, PSFMODEL, ENSEMBLE)
    131 BREAK_POINT           STR  NONE
    132 # PEAKS_OUTPUT_FILE   STR  peaks.dat
    133 # MOMENTS_OUTPUT_FILE STR  moments.dat
    134 # ANALYSIS_REGION     STR  [1000:1600,2800:3400]
     133BREAK_POINT                         STR   NONE     
     134# PEAKS_OUTPUT_FILE                 STR   peaks.dat
     135# MOMENTS_OUTPUT_FILE               STR   moments.dat
     136# ANALYSIS_REGION                   STR   [1000:1600,2800:3400]
    135137
    136138# optional parameter to limit the actual analysis to a fraction of the image
    137139# do not uncomment this in the master psphot.config file
    138 # ANALYSIS_REGION     STR  [1000:1600,2800:3400]
    139 
    140 IGNORE_GROWTH BOOL FALSE
    141 CONSTANT_PHOTOMETRIC_WEIGHTS BOOL TRUE # Should the photometric code [currently only ensemblePSF] refuse to weight each pixel by it's significance?
    142 INTERPOLATE_AP BOOL TRUE
    143 
    144 POISSON.ERRORS.PHOT.LMM BOOL TRUE
    145 POISSON.ERRORS.PHOT.LIN BOOL FALSE
    146 POISSON.ERRORS.PARAMS   BOOL TRUE
    147 
    148 PCM_BOX_SIZE        S32 2
    149 
    150 OUTPUT.FORMAT       STR SMPDATA
    151 NOISE.FACTOR        F32 5.0
    152 NOISE.SIZE          F32 2.0
    153 
    154 USE_FOOTPRINTS                  BOOL    F       # use new pmFootprint peak packaging
    155 FOOTPRINT_NPIXMIN               S32     5       # Minimum size of a pmFootprint
    156 FOOTPRINT_NSIGMA_LIMIT          F32     20      # threshold for bright pmFootprint detection
    157 FOOTPRINT_NSIGMA_LIMIT_2        F32     4       # threshold for faint pmFootprint detection
    158 FOOTPRINT_GROW_RADIUS           S32     3       # How much to grow bright footprints
    159 FOOTPRINT_GROW_RADIUS_2         S32     5       # How much to grow faint footprints
    160 FOOTPRINT_CULL_NSIGMA_DELTA     F32     4       # Cull peaks that aren't nsigma above coll to neighbour
    161 FOOTPRINT_CULL_NSIGMA_MIN       F32     1       # Minimum height of colls in units of skyStdev
     140# ANALYSIS_REGION                   STR   [1000:1600,2800:3400]
     141
     142IGNORE_GROWTH                       BOOL FALSE
     143CONSTANT_PHOTOMETRIC_WEIGHTS        BOOL  TRUE            # Should the photometric code [currently only ensemblePSF] refuse to weight each pixel by it's significance?
     144INTERPOLATE_AP                      BOOL TRUE
     145
     146POISSON.ERRORS.PHOT.LMM             BOOL  TRUE   
     147POISSON.ERRORS.PHOT.LIN             BOOL FALSE
     148POISSON.ERRORS.PARAMS               BOOL TRUE
     149
     150PCM_BOX_SIZE                        S32  2
     151
     152OUTPUT.FORMAT                       STR   SMPDATA     
     153NOISE.FACTOR                        F32  5.0
     154NOISE.SIZE                          F32  2.0
     155
     156USE_FOOTPRINTS                      BOOL  F               # use new pmFootprint peak packaging
     157FOOTPRINT_NPIXMIN                   S32   5               # Minimum size of a pmFootprint
     158FOOTPRINT_NSIGMA_LIMIT              F32   20              # threshold for bright pmFootprint detection
     159FOOTPRINT_NSIGMA_LIMIT_2            F32   4               # threshold for faint pmFootprint detection
     160FOOTPRINT_GROW_RADIUS               S32   3               # How much to grow bright footprints
     161FOOTPRINT_GROW_RADIUS_2             S32   5               # How much to grow faint footprints
     162FOOTPRINT_CULL_NSIGMA_DELTA         F32   4               # Cull peaks that aren't nsigma above coll to neighbour
     163FOOTPRINT_CULL_NSIGMA_MIN           F32   1               # Minimum height of colls in units of skyStdev
    162164
    163165# alternate PSPHOT-type recipes
     
    173175END
    174176
    175 TEST_FIT                        BOOL FALSE
    176 TEST_FIT_MODE                   STR  DEFAULT
    177 TEST_FIT_MODEL                  STR  DEFAULT
    178 TEST_FIT_INNER_RADIUS           F32  NAN
    179 TEST_FIT_OUTER_RADIUS           F32  NAN
    180 TEST_FIT_RADIUS                 F32  NAN
    181 TEST_MOMENTS_RADIUS             F32  NAN
    182 TEST_FIT_X                      F32  NAN
    183 TEST_FIT_Y                      F32  NAN
    184 
    185 TEST_FIT_PAR0                   F32  NAN
    186 TEST_FIT_PAR1                   F32  NAN
    187 TEST_FIT_PAR2                   F32  NAN
    188 TEST_FIT_PAR3                   F32  NAN
    189 TEST_FIT_PAR4                   F32  NAN
    190 TEST_FIT_PAR5                   F32  NAN
    191 TEST_FIT_PAR6                   F32  NAN
    192 TEST_FIT_PAR7                   F32  NAN
    193 TEST_FIT_PAR8                   F32  NAN
    194 TEST_FIT_PAR                    F32  NAN
    195 
    196 DIAGNOSTIC.PLOTS                METADATA
    197   IMAGE.BACKGROUND.CELL.HISTOGRAM   BOOL FALSE
    198   IMAGE.BACKGROUND.CELL.HISTOGRAM.X S32 35
    199   IMAGE.BACKGROUND.CELL.HISTOGRAM.Y S32 -1
    200 END
    201 
    202 PSF.FLUXSCALE.NX                S32 5 # number cells to measure flux scale variations
    203 PSF.FLUXSCALE.NY                S32 5 # number cells to measure flux scale variations
     177TEST_FIT                            BOOL FALSE
     178TEST_FIT_MODE                       STR   DEFAULT
     179TEST_FIT_MODEL                      STR   DEFAULT
     180TEST_FIT_INNER_RADIUS               F32   NAN
     181TEST_FIT_OUTER_RADIUS               F32   NAN
     182TEST_FIT_RADIUS                     F32   NAN
     183TEST_MOMENTS_RADIUS                 F32   NAN
     184TEST_FIT_X                          F32   NAN
     185TEST_FIT_Y                          F32   NAN
     186
     187TEST_FIT_PAR0                       F32   NAN
     188TEST_FIT_PAR1                       F32   NAN
     189TEST_FIT_PAR2                       F32   NAN
     190TEST_FIT_PAR3                       F32   NAN
     191TEST_FIT_PAR4                       F32   NAN
     192TEST_FIT_PAR5                       F32   NAN
     193TEST_FIT_PAR6                       F32   NAN
     194TEST_FIT_PAR7                       F32   NAN
     195TEST_FIT_PAR8                       F32   NAN
     196TEST_FIT_PAR                        F32   NAN
     197
     198DIAGNOSTIC.PLOTS                    METADATA
     199  IMAGE.BACKGROUND.CELL.HISTOGRAM   BOOL  FALSE
     200  IMAGE.BACKGROUND.CELL.HISTOGRAM.X S32   35
     201  IMAGE.BACKGROUND.CELL.HISTOGRAM.Y S32   -1
     202END
     203
     204PSF.FLUXSCALE.NX                    S32   5              # number cells to measure flux scale variations
     205PSF.FLUXSCALE.NY                    S32   5              # number cells to measure flux scale variations
Note: See TracChangeset for help on using the changeset viewer.