Changeset 42382 for trunk/ppImage/src/ppImageDetrendReadout.c

Timestamp:

Feb 8, 2023, 11:57:40 AM (3 years ago)

Author:

eugene

Message:

merge from eam_branches/ipp-20220316. add code for PatternDeadCells, fix organization of DetrendPattern, handle psphotModelBackground failure

Location:

trunk/ppImage

Files:

• . (modified) (1 prop)
• src/ppImageDetrendReadout.c (modified) (2 diffs)

trunk/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?)
+    // 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);
+
+        // 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);
+        psMetadataAddF32(hdu->header, PS_LIST_TAIL, "BACK_ERR", PS_META_REPLACE, "Stdev of cell background", stats->robustStdev);
+        psFree (stats);
+        psFree (rng);
     }