Changeset 2826 for trunk/Ohana/src/opihi/lib.data/starfuncs.c

Timestamp:

Dec 26, 2004, 10:22:20 AM (22 years ago)

Author:

eugene

Message:

added imsmooth, vload, findpeaks

File:

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

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

@@ -6 +6 @@
   double x, y, x2, y2, I, sky, FWHMx, FWHMy, value, mag;
   int i, j, n, Npix2, Nring, Nmax;
-  double Npts, gain, dsky2, dmag;
+  double Npts, gain, dsky2, dmag, peak, offset;
   char *string;
   
@@ -45 +45 @@
   free (ring);
 
+  peak = 0;
   Npts = Nmax = 0;
   x = y = x2 = y2 = I = 0;
@@ -50 +51 @@
     for (j = Y - Npix2; j < Y + Npix2 + 1; j++) {
       value = fits_get_matrix_value (matrix, i, j);
-      x  += i*(value - sky);
-      y  += j*(value - sky);
-      x2 += i*i*(value - sky);
-      y2 += j*j*(value - sky);
-      I  += (value - sky);
+      offset = value - sky;
+      x  += i*offset;
+      y  += j*offset;
+      x2 += i*i*offset;
+      y2 += j*j*offset;
+      I  += offset;
       Npts ++;
       if (value > max) {
         Nmax ++;
       }
+      if (value > peak) peak = value;
     }
   }
@@ -80 +83 @@
   set_variable ("dZg", dmag);
   set_variable ("Zcg", I);
+  set_variable ("Zpk", peak);
   set_int_variable ("Nsat", Nmax);