Changeset 42158

Timestamp:

Apr 4, 2022, 4:01:04 PM (4 years ago)

Author:

eugene

Message:

test code to measure and record the background level (robust median) for each cell

File:

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

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

@@ -63 +63 @@
     // psLogMsg ("ppImage", 6, "check video: %f sec\n", psTimerMark ("detrend.readout"));
 
+    // XXX for GPC2, there are some bizzare cell levels that need to be filtered out
+    // Here I measure the median background before overscan subtraction, but an alternative
+    // is to use the overscan-subtracted image.
+    // XX if (0) {
+    // XX       psStats *stats = psStatsAlloc (PS_STAT_ROBUST_MEDIAN);
+    // XX       psRandom *rng = psRandomAlloc(PS_RANDOM_TAUS);
+    // XX       psImageBackground (stats, NULL, input->image, NULL, 0xffff, rng);
+    // XX       readoutBackground = stats->robustMedian;
+    // XX }
+
     // Masking on the basis of pixel value needs to be done before anything else, so the values are pristine.
     if (options->doMaskBuild) {
@@ -113 +123 @@
       }
       // psLogMsg ("ppImage", 6, "subtract overscan: %f sec\n", psTimerMark ("detrend.readout"));
+    }
+
+    // measure the overscan-subtracted readoutBackground here (or subtract the overscan value?)
+    if (1) {
+        psStats *stats = psStatsAlloc (PS_STAT_ROBUST_MEDIAN);
+        psRandom *rng = psRandomAlloc(PS_RANDOM_TAUS);
+        psImageBackground (stats, NULL, input->image, NULL, 0xffff, rng);
+
+        // save this value somewhere
+        pmHDU *hdu = pmHDUFromReadout(input);  // HDU of interest
+        psMetadataAddF32(hdu->header, PS_LIST_TAIL, "BACK_VAL", PS_META_REPLACE, "Median cell background", stats->robustMedian);
+        psFree (stats);
+        psFree (rng);
     }