Changeset 28797 for branches/eam_branches/ipp-20100621/psphot/src/psphotSignificanceImage.c

Timestamp:

Jul 30, 2010, 9:39:15 AM (16 years ago)

Author:

eugene

Message:

updates from trunk

Location:

branches/eam_branches/ipp-20100621/psphot

Files:

• . (modified) (1 prop)
• src/psphotSignificanceImage.c (modified) (1 diff)

branches/eam_branches/ipp-20100621/psphot/src/psphotSignificanceImage.c

@@ -76 +76 @@
     // Calculate correction factor for the covariance produced by the (potentially multiple) smoothing
     psKernel *kernel = psImageSmoothKernel(SIGMA_SMTH, NSIGMA_SMTH); // Kernel used for smoothing
-    float factor = 1.0 / psImageCovarianceCalculateFactor(kernel, readout->covariance);
+    double sum2 = 0.0;                                               // Sum of kernel squared
+    for (int y = kernel->yMin; y <= kernel->yMax; y++) {
+        for (int x = kernel->xMin; x <= kernel->xMax; x++) {
+            sum2 += PS_SQR(kernel->kernel[y][x]);
+        }
+    }
+    float factor = 1.0 / (sum2 * psImageCovarianceCalculateFactor(kernel, readout->covariance));
     psFree(kernel);