Changeset 19826

Timestamp:

Oct 2, 2008, 8:11:46 AM (18 years ago)

Author:

eugene

Message:

subtract guess position when calculating sums

File:

• trunk/Ohana/src/opihi/lib.data/starfuncs.c (modified) (2 diffs)

trunk/Ohana/src/opihi/lib.data/starfuncs.c

@@ -44 +44 @@
   dsky2 = dsky2 / Npts - sky*sky;
   free (ring);
+
+  float dx, dy;
 
   peak = 0;
@@ -52 +54 @@
       value = gfits_get_matrix_value (matrix, i, j);
       offset = value - sky;
-      x  += i*offset;
-      y  += j*offset;
-      x2 += i*i*offset;
-      y2 += j*j*offset;
-      xy += i*j*offset;
+      dx = i - X;
+      dy = j - Y;
+      x  += dx*offset;
+      y  += dy*offset;
+      x2 += dx*dx*offset;
+      y2 += dy*dy*offset;
+      xy += dx*dy*offset;
       I  += offset;
       Npts ++;