Changeset 24295 for trunk/psModules/src/imcombine/pmSubtractionKernels.c

Timestamp:

Jun 1, 2009, 5:45:42 PM (17 years ago)

Author:

Paul Price

Message:

Fixing normalisation: was normalising both x and y by the sum of the products.

File:

• trunk/psModules/src/imcombine/pmSubtractionKernels.c (modified) (2 diffs)

trunk/psModules/src/imcombine/pmSubtractionKernels.c

@@ -145 +145 @@
 
                 // Calculate moments
-                double sum = 0.0;       // Sum of kernel component, for normalisation
                 double moment = 0.0;    // Moment, for penalty
                 for (int v = -size, y = 0; v <= size; v++, y++) {
                     for (int u = -size, x = 0; u <= size; u++, x++) {
                         double value = xKernel->data.F32[x] * yKernel->data.F32[y]; // Value of kernel
-                        sum += value;
                         moment += value * (PS_SQR(u) + PS_SQR(v));
                     }
@@ -163 +161 @@
                         }
                     }
-                    sum = 1.0 / sum;
+                    sum = 1.0 / sqrt(sum);
                     psBinaryOp(xKernel, xKernel, "*", psScalarAlloc(sum, PS_TYPE_F32));
                     psBinaryOp(yKernel, yKernel, "*", psScalarAlloc(sum, PS_TYPE_F32));