Changeset 25938 for branches/pap/psModules/src/imcombine/pmStack.c

Timestamp:

Oct 23, 2009, 2:45:37 PM (17 years ago)

Author:

Paul Price

Message:

When TESTING, make output mask = number of valid inputs, variance =
fraction of valid inputs unrejected.

File:

• branches/pap/psModules/src/imcombine/pmStack.c (modified) (7 diffs)

branches/pap/psModules/src/imcombine/pmStack.c

@@ -33 +33 @@
 #define NUM_DIRECT_STDEV 5              // For less than this number of values, measure stdev directly
 
-//#define TESTING                         // Enable test output
-//#define TEST_X 1085                     // x coordinate to examine
-//#define TEST_Y 3371                     // y coordinate to examine
+#define TESTING                         // Enable test output
+#define TEST_X -1                     // x coordinate to examine
+#define TEST_Y -1                     // y coordinate to examine
 
 
@@ -390 +390 @@
         }
     }
+    int numRejected = 0;                // Number of rejected inputs
 #endif
 
@@ -423 +424 @@
               fprintf(stderr, "Single input to combine, safety on, pixel is bad.\n");
           }
+          numRejected = 1;
 #endif
           break;
@@ -468 +470 @@
                               diff, rej, sqrtf(sigma2));
                   }
+                  numRejected = 2;
 #endif
               }
@@ -579 +582 @@
                                       j, diff, sqrtf(pixelLimits->data.F32[j]));
                           }
+                          numRejected++;
 #endif
                       }
@@ -588 +592 @@
                                   j, diff, limit);
                       }
+                      numRejected++;
 #endif
                   }
@@ -632 +637 @@
         variance->data.F32[y][x] = varianceValue;
     }
+
+#ifdef TESTING
+    mask->data.PS_TYPE_IMAGE_MASK_DATA[y][x] = num;
+    if (variance) {
+        variance->data.F32[y][x] = 1.0 - (float)numRejected / (float)num;
+    }
+#endif
 
     return;