Changeset 29595

Timestamp:

Oct 28, 2010, 5:42:58 PM (16 years ago)

Author:

eugene

Message:

separate the calculation of the kernel normalization from the moments; calculate the moments using fabs(flux), since the kernel may have negative flux

File:

• trunk/psModules/src/imcombine/pmSubtractionAnalysis.c (modified) (1 diff)

trunk/psModules/src/imcombine/pmSubtractionAnalysis.c

@@ -105 +105 @@
         int fullSize = 2 * size + 1;    // Full size of kernel
 
-        /* in the CENTRAL_DELTA case, the final convolution kernel replaces the central delta
-           function which was subtracted to force zero flux for each kernel.  for other cases,
-           we have to include the zero order component of the normalization here.
-         */
-
-        //# if (CENTRAL_DELTA)
-# if (1)
         float norm = 0.0;               // Normalisation (kernel sum)
-# else
-        float norm = 1.0;               // Normalisation (kernel sum)
-# endif
         for (int y = 0; y < fullSize; y++) {
             for (int x = 0; x < fullSize; x++) {