Changeset 30266 for branches/eam_branches/ipp-20101205/psModules/src/imcombine/pmSubtractionStamps.c

Timestamp:

Jan 12, 2011, 9:32:08 PM (16 years ago)

Author:

eugene

Message:

add errors to kernel coeffs; calculate and plot chisq and moments

File:

• branches/eam_branches/ipp-20101205/psModules/src/imcombine/pmSubtractionStamps.c (modified) (1 diff)

branches/eam_branches/ipp-20101205/psModules/src/imcombine/pmSubtractionStamps.c

@@ -812 +812 @@
     stamps->normWindow2 = 2.0*R2;
     psLogMsg ("psModules.imcombine", PS_LOG_DETAIL, "Kron Radii: %f for 1, %f for 2\n", stamps->normWindow1, stamps->normWindow2);
+
+    // Generate a weighting window based on the kron radii
+    { 
+        float radius = 1.5 * PS_MAX(R1, R2);
+
+        psImageInit(stamps->window->image, 0.0);
+
+        for (int y = -size; y <= size; y++) {
+            for (int x = -size; x <= size; x++) {
+                if (hypot(x,y) > radius) continue;
+                stamps->window->kernel[y][x] = 1.0;
+            }
+        }
+    }
 
 # else