Changeset 26341

Timestamp:

Dec 5, 2009, 5:03:10 AM (16 years ago)

Author:

eugene

Message:

show fractional residual image (chisq contribution)

File:

• branches/eam_branches/20091201/psModules/src/imcombine/pmSubtractionVisual.c (modified) (2 diffs)

branches/eam_branches/20091201/psModules/src/imcombine/pmSubtractionVisual.c

@@ -352 +352 @@
     for (int y = -footprint; y <= footprint; y++) {
         for (int x = -footprint; x <= footprint; x++) {
-            fresidualImage->data.F32[y + NYpix][x + NXpix] = residualImage->data.F32[y + NYpix][x + NXpix] / target->kernel[y][x];
+            fresidualImage->data.F32[y + NYpix][x + NXpix] = PS_SQR(residualImage->data.F32[y + NYpix][x + NXpix]) / (target->kernel[y][x] + 25.0);
         }
     }
@@ -364 +364 @@
     if (!pmVisualInitWindow(&kapa2, "ppSub:Misc")) return false;
 
+    KiiEraseOverlay (kapa1, "red");
+    KiiEraseOverlay (kapa2, "red");
+
     pmVisualScaleImage(kapa1, targetImage, "Target Stamps", 0, true);
     pmVisualScaleImage(kapa1, sourceImage, "Source Stamps", 1, true);
     pmVisualScaleImage(kapa1, convolutionImage, "Convolution Stamps", 2, true);
 
-    // pmVisualScaleImage(kapa2, fresidualImage, "Frac Residual Stamps", 2, true);
+    pmVisualScaleImage(kapa2, fresidualImage, "Frac Residual Stamps", 2, true);
     pmVisualScaleImage(kapa2, residualImage, "Residual Stamps", 1, true);
     pmVisualScaleImage(kapa2, differenceImage, "Difference Stamps", 0, true);