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

Timestamp:

Dec 2, 2021, 8:23:14 AM (5 years ago)

Author:

eugene

Message:

oops, need to use nGoodClip not nGood

File:

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

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

@@ -1657 +1657 @@
             int Nhi = nGood; 
 
+            if ((x == 3161) && (y == 450)) {
+              fprintf (stderr, "test point\n");
+            }
+
             if (nGood >= 5) {
               int midPoint = nGood / 2;
@@ -1665 +1669 @@
               int P75 = 0.75*nGood;
               float rawSigma = 0.74*(pixelData->data.F32[P75] - pixelData->data.F32[P25]);
-              float minThresh = rawMedian - 7.0*rawSigma;
-              float maxThresh = rawMedian + 7.0*rawSigma;
+              float minThresh = rawMedian - 5.0*rawSigma;
+              float maxThresh = rawMedian + 5.0*rawSigma;
 
               // find the entries which are in the range
@@ -1691 +1695 @@
 # if (AT_LEAST)
             int Ns = MIN(MAX(1, 0.5*rejectFraction * nGoodClip), nGoodClip) + Nlo;
-            int Ne = nGood - Ns + Nlo;
+            int Ne = nGoodClip - Ns + Nlo;
             int Npt = Ne - Ns;
 # else
-            int Ns = MIN(MAX(0, 0.5*rejectFraction * nGood), nGood) + Nlo;
-            int Ne = nGood - Ns + Nlo;
+            int Ns = MIN(MAX(0, 0.5*rejectFraction * nGoodClip), nGoodClip) + Nlo;
+            int Ne = nGoodClip - Ns + Nlo;
             int Npt = Ne - Ns;
 # endif