Changeset 23773 for trunk/psModules/src/imcombine/pmStack.c

Timestamp:

Apr 9, 2009, 12:13:16 PM (17 years ago)

Author:

Paul Price

Message:

Weighted Olympic mean wasn't working as intended: high pixels weren't being discarded correctly.

File:

• trunk/psModules/src/imcombine/pmStack.c (modified) (3 diffs)

trunk/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 2318                     // x coordinate to examine
-//#define TEST_Y 2306                     // y coordinate to examine
+#define TESTING                         // Enable test output
+#define TEST_X 1085                     // x coordinate to examine
+#define TEST_Y 3371                     // y coordinate to examine
 
 
@@ -265 +265 @@
 
     int numBad = frac * numGood + 0.5;  // Number of bad values
-    int low = numBad / 2, high = low + numGood; // Indices (modulo masked pixels) delimiting range of interest
+    int low = numBad / 2, high = low + numGood - numBad; // Indices (modulo masked pixels)
 
     sortBuffer = psVectorSortIndex(sortBuffer, values);
@@ -536 +536 @@
 #ifdef TESTING
                   if (x == TEST_X && y == TEST_Y) {
-                      fprintf(stderr, "Rejection limit: %f\n", limit);
+                      fprintf(stderr, "Rejecting without variance; rejection limit: %f\n", limit);
                   }
 #endif
               } else {
+#ifdef TESTING
+                  if (x == TEST_X && y == TEST_Y) {
+                      fprintf(stderr, "Rejecting with variance...\n");
+                  }
+#endif
                   median = combinationWeightedOlympic(pixelData, pixelWeights, pixelMasks, discard, sort);
               }