Changeset 10175

Timestamp:

Nov 23, 2006, 5:14:54 AM (20 years ago)

Author:

rhl

Message:

Fix smoothing (and normalisation) of weight map --- we're adding them in quadrature

File:

• trunk/psphot/src/psphotFindPeaks.c (modified) (1 diff)

trunk/psphot/src/psphotFindPeaks.c

@@ -23 +23 @@
 
     psImage *smooth_wt = psImageCopy (NULL, readout->weight, PS_TYPE_F32);
-    psImageSmooth (smooth_wt, SIGMA, NSIGMA);
+    psImageSmooth (smooth_wt, SIGMA/sqrt(2), NSIGMA);
+    smooth_wt = (psImage*)psBinaryOp(smooth_wt, smooth_wt, "*",
+                                     psScalarAlloc(1/(4*M_PI*PS_SQR(SIGMA)), PS_TYPE_F32));
     psLogMsg ("psphot", 4, "smooth weight: %f sec\n", psTimerMark ("psphot"));