Changeset 26007 for branches/pap/psModules/src/imcombine/pmSubtraction.c

Timestamp:

Nov 2, 2009, 5:08:31 PM (17 years ago)

Author:

Paul Price

Message:

Serious rework of pmStackCombine for clarity and some additional functionality. Split the combinePixels function into two operating modes (test and combine) to make clear what we're doing; pmStackCombine could be split into two separate functions (e.g., pmStackTest and pmStackCombine) in the future, since the testing part doesn't need to produce any output image, but this involves moving a lot of stuff around in ppStack, so not doing it yet. Verified that variance scaling is correct, and softening of variance is correct. Big improvement in minimising loss of good pixels comes from throwing out only the most variant pixel on each rejection pass (consequently changed iterations parameter to iterations per input); this involved another rework of the core combinePixels function (moving the switch on the number of valid inputs inside the iteration). Much happier with the outlier rejection now.

File:

• branches/pap/psModules/src/imcombine/pmSubtraction.c (modified) (1 diff)

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

@@ -832 +832 @@
     psTrace("psModules.imcombine", 1, "RMS deviation: %f\n", rms);
 
+    fprintf(stderr, "Mean = %f ; chi^2 = %f for %ld d.o.f.\n",
+            mean,
+            mean * PS_SQR(2 * stamps->footprint + 1) * numStamps,
+            PS_SQR(2 * stamps->footprint + 1) * numStamps - kernels->solution1->n);
+
+
     kernels->mean = mean;
     kernels->rms = rms;