Changeset 42167 for branches/eam_branches/ipp-20220316/ppImage/src/ppImageDetrendReadout.c

Timestamp:

Apr 13, 2022, 10:53:25 AM (4 years ago)

Author:

eugene

Message:

add PATTERN.DEAD.CELLS option; split ppImageDetrendPattern into separate functions for the 4 different options; functions like ppImageDetrendFringeApply and ppImageDetrendPatternApply now test the relevant option and return true if not selected (no need to check the option in the calling function)

File:

• branches/eam_branches/ipp-20220316/ppImage/src/ppImageDetrendReadout.c (modified) (2 diffs)

branches/eam_branches/ipp-20220316/ppImage/src/ppImageDetrendReadout.c

@@ -126 +126 @@
 
     // measure the overscan-subtracted readoutBackground here (or subtract the overscan value?)
-    if (1) {
-        psStats *stats = psStatsAlloc (PS_STAT_ROBUST_MEDIAN);
+    // XXX this is the measurements and should be independent of the pattern masking
+    if (options->doPatternDeadCells) {
+        psStats *stats = psStatsAlloc (PS_STAT_ROBUST_MEDIAN | PS_STAT_ROBUST_STDEV);
         psRandom *rng = psRandomAlloc(PS_RANDOM_TAUS);
         psImageBackground (stats, NULL, input->image, NULL, 0xffff, rng);
@@ -134 +135 @@
         pmHDU *hdu = pmHDUFromReadout(input);  // HDU of interest
         psMetadataAddF32(hdu->header, PS_LIST_TAIL, "BACK_VAL", PS_META_REPLACE, "Median cell background", stats->robustMedian);
+        psMetadataAddF32(hdu->header, PS_LIST_TAIL, "BACK_ERR", PS_META_REPLACE, "Stdev of cell background", stats->robustStdev);
         psFree (stats);
         psFree (rng);