IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changes between Version 1 and Version 2 of Basic_GPC1_Processing_for_End_Users


Ignore:
Timestamp:
Feb 26, 2009, 11:00:23 AM (17 years ago)
Author:
trac
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • Basic_GPC1_Processing_for_End_Users

    v1 v2  
    4242
    4343<code>
    44 cd MD01_20081002_1/36475
    45 psphot -D CAMERA gpc1 -file chip/o4741g0234o.36475.ch.11948.XY55.ch.fits -mask chip/o4741g0234o.36475.ch.11948.XY55.ch.mk.fits -weight chip/o4741g0234o.36475.ch.11948.XY55.ch.wt.fits testout
     44 cd MD01_20081002_1/36475
     45 psphot -D CAMERA gpc1 -file chip/o4741g0234o.36475.ch.11948.XY55.ch.fits \
     46   -mask chip/o4741g0234o.36475.ch.11948.XY55.ch.mk.fits \
     47   -weight chip/o4741g0234o.36475.ch.11948.XY55.ch.wt.fits \
     48   testout
    4649</code>
    4750will create testout.cmf.  If you want to see what psphot is doing as it runs, you may be interested in the '-visual' option
    4851
    49 psphot -file chip/o4741g0234o.36475.ch.11948.XY55.ch.fits -mask chip//o4741g0234o.36475.ch.11948.XY55.ch.mk.fits -weight chip/o4741g0234o.36475.ch.11948.XY55.ch.wt.fits testout -visual
     52<code>
     53 psphot -file chip/o4741g0234o.36475.ch.11948.XY55.ch.fits \
     54        -mask chip//o4741g0234o.36475.ch.11948.XY55.ch.mk.fits \
     55        -weight chip/o4741g0234o.36475.ch.11948.XY55.ch.wt.fits \
     56        -visual \
     57        testout
     58</code>
    5059
    5160The main output from the psphot command will be the .cmf file, which is a binary FITS table that contains information about the objects as gleaned from the image.  One of the fields that generates the most questions is the FLAGS value for objects.  This is a bitmask.  See [wiki:IPP_Detection_Bitmasks IPP Detection Bitmasks] for more details.
     61
     62 * psphot options.
     63
     64There are a variety of options and settings for psphot.  These are controlled by command line optional arguments as well as by the "recipes" file located in the "share/ippconfig/" directory in your IPP installation.  On my machine this is "/Volumes/data/PS1/code/ipp/default.darwin/share/ippconfig".  The default recipes are in "share/ippconfig/recipes", with camera-specific overriding configuration files in the camera subdirectories in "share/ippconfig", e.g. "share/ippconfig/gpc1" for the PS1 GPC1 camera.
     65
     66Here's an example of a recipes file for psphot:
     67
     68<pre>
     69SAVE.OUTPUT                         BOOL  TRUE
     70SAVE.BACKMDL                        BOOL  FALSE
     71SAVE.BACKMDL.STDEV                  BOOL  FALSE
     72SAVE.BACKGND                        BOOL  FALSE
     73SAVE.BACKSUB                        BOOL  FALSE
     74SAVE.RESID                          BOOL  FALSE
     75SAVE.PSF                            BOOL  FALSE
     76LOAD.PSF                            BOOL  FALSE
     77SAVE.PLOTS                          BOOL  FALSE
     78
     79# the zero point is used to set a basic scale for DVO
     80# XXX it may not currently be read : double check this (EAM)
     81ZERO_POINT                          F32   25.000          # zero point used by DVO
     82ZERO_PT                             F32   25.000          # zero point used by DVO
     83
     84OUTPUT.FORMAT                       STR   PS1_DEV_1
     85
     86# these parameter govern how the background is measured
     87BACKGROUND.XBIN                     S32   128             # size of background superpixels
     88BACKGROUND.YBIN                     S32   128             # size of background superpixels
     89IMSTATS_NPIX                        S32   10000           # number of pixels to use for sky estimate boxes:
     90
     91SKY_BIAS                            F32   0.0             # offset applied to measured sky (FOR TESTING)
     92SKY_FIT_ORDER                       S32   0
     93SKY_FIT_LINEAR                      BOOL  FALSE
     94SKY_STAT                            STR   FITTED_MEAN_V4  # statistic used to measure background
     95SKY_CLIP_SIGMA                      F32   2.0             # statistic used to measure background
     96SKY_SIG                             F32   1.0             # optional sky error for
     97
     98# allowed values for SKY_STAT:
     99# SAMPLE_MEAN, SAMPLE_MEDIAN, CLIPPED_MEAN, ROBUST_MEAN, ROBUST_QUARTILE, FITTED_MEAN
     100
     101# masking parameters (XXX EAM : rework this to use psRegion like ANALYSIS_REGION)
     102XMIN                                F32   0               # minimum valid x-coord
     103XMAX                                F32   0               # maximum valid x-coord
     104YMIN                                F32   0               # minimum valid y-coord
     105YMAX                                F32   0               # maximum valid y-coord
     106
     107# peak finding
     108PEAKS_SMOOTH_SIGMA                  F32   1.0             # smoothing kernel sigma in pixels
     109PEAKS_SMOOTH_NSIGMA                 F32   2.0             # smoothing kernel width in sigmas
     110PEAKS_NSIGMA_LIMIT                  F32   25.0            # peak significance threshold
     111PEAKS_NSIGMA_LIMIT_2                F32   5.0             # peak significance threshold
     112PEAKS_NMAX                          S32   0               # on first pass, only keep NMAX peaks (0 == all)
     113PEAKS_MIN_GAUSS                     F32   0.5             # Minimum valid fraction of kernel
     114
     115# parameters to control the selection of the peak in the Sx,Sy plane
     116MOMENTS_SCALE                       F32   0.25       
     117MOMENTS_SN_MIN                      F32   100.0           # min S/N to measure moments
     118MOMENTS_SX_MAX                      F32   50.0
     119MOMENTS_SY_MAX                      F32   50.0
     120MOMENTS_AR_MAX                      F32   1.5             # maximum axial ratio: 1 / AR < (sx / sy) < AR
     121
     122# basic object statistics
     123SKY_INNER_RADIUS                    F32   15              # square annulus for local sky measurement
     124SKY_OUTER_RADIUS                    F32   25              # square annulus for local sky measurement
     125PSF_MOMENTS_RADIUS                  F32   3               # calculate initial source moments with this radius
     126PSF_SN_LIM                          F32   50              # minimum S/N for stars used for PSF model
     127PSF_MAX_NSTARS                      S32   200             # limit number of stars used for PSF model
     128PSF_CLUMP_NSIGMA                    F32   1.5             # region of Sx,Sy plane to use for selecting PSF stars
     129PSF_CLUMP_GRID_SCALE                F32   0.5             # size of Sx,Sy image pixel
     130
     131PSF_CLUMP_NX                        S32   1               # subdivide image in to NX x NY regions for
     132PSF_CLUMP_NY                        S32   1               # selecting PSF stars
     133</pre>
     134
     135
     136
     137Currently 'psphot' does not output it's options or have an help message.  This is something that is being worked on and in fact might even be a result of our efforts to understand psphot better.
    52138
    53139 * DVO.  Let's say we have a set of cmf files and we want to use DVO to combine the detections from these files.  A series of commands such as