Changeset 25489 for trunk/psModules/src/imcombine/pmPSFEnvelope.c

Timestamp:

Sep 22, 2009, 5:28:06 PM (17 years ago)

Author:

Paul Price

Message:

Normalise by the peak, not the centre.

File:

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

trunk/psModules/src/imcombine/pmPSFEnvelope.c

@@ -170 +170 @@
             float y = source->peak->yf + yOffset->data.S32[j]; // y coordinate of source
 
-            double flux = fakeRO->image->data.F32[(int)y][(int)x];
+            // Perhaps I'm being paranoid, but specify a range to check
+            int uMax = PS_MIN(x + radius, numCols - 1), uMin = PS_MAX(x - radius, 0);
+            int vMax = PS_MIN(y + radius, numRows - 1), vMin = PS_MAX(y - radius, 0);
+
+            double flux = -INFINITY;    // Peak flux
+            for (int v = vMin; v <= vMax; v++) {
+                for (int u = uMin; u <= uMax; u++) {
+                    if (fakeRO->image->data.F32[v][u] > flux) {
+                        flux = fakeRO->image->data.F32[v][u];
+                    }
+                }
+            }
+
             if (!isfinite(flux) || flux < 0) {
                 continue;
@@ -240 +252 @@
         int xFake = source->peak->x + xOffset->data.S32[i];
         int yFake = source->peak->y + yOffset->data.S32[i];
-        psRegion region = psRegionSet(xFake - radius, xFake + radius,
-                                      yFake - radius, yFake + radius); // PSF region
+        psRegion region = psRegionSet(xFake - maxRadius, xFake + maxRadius,
+                                      yFake - maxRadius, yFake + maxRadius); // PSF region
         psImage *subImage = psImageSubset(envelope, region); // Subimage of fake PSF