Changeset 18190

Timestamp:

Jun 18, 2008, 5:03:53 PM (18 years ago)

Author:

Paul Price

Message:

Need to correct for the non-unity scaling of the kernel.

File:

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

trunk/psModules/src/imcombine/pmSubtraction.c

@@ -4 +4 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.94 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2008-06-17 22:16:38 $
+ *  @version $Revision: 1.95 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2008-06-19 03:03:53 $
  *
  *  Copyright 2004-2007 Institute for Astronomy, University of Hawaii
@@ -642 +642 @@
 
     double sumKernel2 = 0.0;            // Sum of the kernel squared
+    double sumKernel = 0.0;             // Sum of the kernel
     for (int y = kernel->yMin; y <= kernel->yMax; y++) {
         for (int x = kernel->xMin; x <= kernel->xMax; x++) {
+            sumKernel += kernel->kernel[y][x];
             sumKernel2 += PS_SQR(kernel->kernel[y][x]);
         }
@@ -650 +652 @@
     psFree(kernel);
 
-    return sumKernel2;
+    return sumKernel2 / PS_SQR(sumKernel);
 }