Changeset 26894 for trunk/psphot/doc/stack.txt
- Timestamp:
- Feb 10, 2010, 7:36:29 PM (16 years ago)
- File:
-
- 1 edited
-
trunk/psphot/doc/stack.txt (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/doc/stack.txt
r24583 r26894 1 2 20100126: 3 4 * watch out for psphotSetMomentsWindow & MOMENTS_SX_MAX,etc 5 * watch out for psphotSignificanceImage.c:,psphotEfficiency using the FWHM_MAJ from psphotChoosePSF 6 * ppSimDetections.c : SIGMA_SMOOTH 7 ppSmooth/src/ppSmoothReadout.c: psMetadataAddF32(recipe, PS_LIST_TAIL, "EFFECTIVE_AREA", PS_META_REPLACE, "Effective Area", effArea); 8 ppSmooth/src/ppSmoothReadout.c: psMetadataAddF32(recipe, PS_LIST_TAIL, "SIGNIFICANCE_SCALE_FACTOR", PS_META_REPLACE, "Signicance scale factor", factor); 9 10 11 20100120 : more stack processing mods: 12 13 there are a number of data collections generated during psphotReadout: 14 15 * detections (currently a structure containing multiple arrays) 16 * sources (a psArray) 17 * psf (the psf model structure) 18 19 * new sources vs old sources: 20 there is a sequence on the second pass in which we need to distinguish the sources 21 from the first pass from those on the second pass: 22 23 - add noise (old sources only) 24 - find detections 25 - subtract noise (old sources only) 26 - generate sources (new detections only) 27 - source classifications (new sources only) 28 - guess models (new sources only) 29 - replace sources (old sources only) 30 - merge sources (new + old -> sources) 31 32 currently we are distiguishing the old vs new based on different arrays. 33 can we use the processing flags to distinguish the these cases and carry around 34 only a single source list? 35 36 * detections->peaks holds the most recently detected set of peaks 37 detections->oldPeaks holds the previous collection 38 39 detections->footprints holds the full set of merged footprints, 40 including assigned peaks from the old and new set. 41 42 20100107 : updates for stack processing 43 44 I am building a new version of psphot that operates on a set of 45 images (eg, multi-filter or multi-epoch). The assumptions are: 46 47 * each input image represents the same sky pixels : they are warped to a common frame. 48 * each image has been previously processed, with the background 49 subtracted (and the psf model determined?) 50 * all planes are / can be loaded into memory at once (otherwise I 51 need to add a lot of I/O layers) 52 53 psphot currently expects the pmFPAfile of interest to be available in 54 config->files with the name PSPHOT.INPUT. To update the code for 55 stacks, I am extending this concept with the pmFPAfileSelectSingle 56 API: config->files may contain multiple PSPHOT.INPUT entries, and 57 functions which access this file / these files need to specify 58 *which* image they want. The number of PSPHOT.INPUT entries is 59 saved in config->arguments as PSPHOT.INPUT.NUM (this can be 60 generated from the filerule PSPHOT.INPUT for extension). 61 62 As a result, we need a number of wrapper functions which loop over 63 all PSPHOT.INPUT.NUM entries and perform a particular operation on 64 one of the entries. Here are the functions which I have modified 65 in this way (function -> child) 66 67 * psphotAddPhotcode -> psphotAddPhotcodeReadout 68 * note that the photcode is now saved on readout->analysis 69 * psphotSetMaskAndVariance -> psphotSetMaskAndVarianceReadout 70 * psphotModelBackground -> psphotModelBackgroundReadoutFileIndex 71 * psphotSubtractBackground -> psphotSubtractBackgroundReadout 72 73 side notes on 74 75 * psphotModelBackground vs psphotBackgroundModel : I've renamed 76 psphotBackgroundModel (used only by ppStack for now) to 77 psphotModelBackgroundReadoutNoFile. I've also named the 78 single-readout version used by psphotModelBackground to 79 psphotModelBackgroundReadoutFileIndex. These two functions do the 80 same thing, but psphotModelBackgroundReadoutNoFile does not save 81 a pmFPAfile on config->files 82 1 83 2 84 20090606 : design notes on the multi-image photometry analysis
Note:
See TracChangeset
for help on using the changeset viewer.
