Changeset 17351

Timestamp:

Apr 6, 2008, 2:42:31 PM (18 years ago)

Author:

eugene

Message:

move psf normalization into psphotKernelFromPSF

File:

• branches/eam_branch_20080324/psphot/src/psphotKernelFromPSF.c (modified) (1 diff)

branches/eam_branch_20080324/psphot/src/psphotKernelFromPSF.c

@@ -12 +12 @@
   psKernel *psf = psKernelAlloc (-nPix, +nPix, -nPix, +nPix);
 
+  double sum = 0.0;
+
   for (int j = psf->yMin; j <= psf->yMax; j++) {
     for (int i = psf->xMin; i <= psf->xMax; i++) {
-      psf->kernel[j][i] = source->psfFlux->data.F32[y0 + j][x0 + i];
+        double value = source->psfFlux->data.F32[y0 + j][x0 + i];
+        psf->kernel[j][i] = value;
+        sum += value;
     }
   }
+  assert (sum > 0.0);
+
+  // psf must be normalized (integral = 1.0)
+  for (int i = 0; i < psf->image->numRows; i++) {
+      for (int j = 0; j < psf->image->numCols; j++) {
+          psf->image->data.F32[i][j] /= sum;
+      }
+  }
+
   return psf;
 }