Changeset 41933

Timestamp:

Nov 24, 2021, 7:11:32 AM (5 years ago)

Author:

eugene

Message:

more restrictive clipping for percentile method: MIN_GOOD = 3, exclude suspect pixels, use AT_LEAST method: these should all be user options

File:

• branches/eam_branches/ipp-20211108/psModules/src/imcombine/pmStack.c (modified) (4 diffs)

branches/eam_branches/ipp-20211108/psModules/src/imcombine/pmStack.c

@@ -1515 +1515 @@
 }
 
-# define MIN_GOOD_PERCENTILE 2
 # define SUSPECT_FRACTION 0.65
 
@@ -1610 +1609 @@
                 if (fabs(image->data.F32[yIn][xIn]) > 1e5) continue; // XXX this cut is a bit arbitrary..
                 if (mask->data.PS_TYPE_IMAGE_MASK_DATA[yIn][xIn] & badMaskBits) continue;
+                if (mask->data.PS_TYPE_IMAGE_MASK_DATA[yIn][xIn] & suspectMaskBits) continue;
 
                 // count the number of times a given mask bit is set in the input pixels.
@@ -1629 +1629 @@
                 // accumulate exposure times if required
                 if (expTime)  { expTime->data.F32[nGood] = data->exp; }
-
                 nGood ++;
             }
@@ -1641 +1640 @@
             // rejection needs to be symmetric
 
-# define AT_LEAST 0
+            // 'AT_LEAST' means we reject at least 'rejectFraction' of values, but it could
+            // be a higher percentage for smaller numbers of inputs.
+# define AT_LEAST 1
+# define MIN_GOOD_PERCENTILE 3
 # if (AT_LEAST)
             int Ns = MIN(MAX(1, 0.5*rejectFraction * nGood), nGood);