Changeset 41954

Timestamp:

Dec 2, 2021, 3:58:54 PM (5 years ago)

Author:

eugene

Message:

add deactivated test pixels, allow values > 1e5

File:

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

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

@@ -1617 +1617 @@
                 // skip obviously bad input data
                 if (!isfinite(image->data.F32[yIn][xIn])) continue;
-                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;
@@ -1657 +1656 @@
             int Nhi = nGood; 
 
-            if ((x == 3161) && (y == 450)) {
-              fprintf (stderr, "test point\n");
-            }
-
             if (nGood >= 5) {
               int midPoint = nGood / 2;
@@ -1682 +1677 @@
             int nGoodClip = Nhi - Nlo;
 
+            int isTest = false;
+            isTest = isTest || ((x == 4743) && (y == 2903));
+            isTest = isTest || ((x == 4953) && (y == 2919));
+            isTest = isTest || ((x == 4751) && (y == 2725));
+            isTest = false;
+
+            if (isTest) {
+              char testname[256];
+              snprintf (testname, 256, "testpix.%04d.%04d.txt", (int) x, (int) y);
+              FILE *f = fopen (testname, "w");
+              int fd = fileno (f);
+
+              fprintf (f, "# nGood: %d, Nlo: %d, Nhi: %d, nGoodClip: %d\n", nGood, Nlo, Nhi, nGoodClip);
+              p_psVectorPrint (fd, pixelData, "pixelData");
+              fclose (f);
+            }
             // rather than define a min and max value,
             // what we really want is a symmetric selection about the middle,