Changeset 2938 for trunk/Ohana/src/opihi/cmd.data/gridify.c

Timestamp:

Jan 8, 2005, 2:45:18 PM (22 years ago)

Author:

eugene

Message:

minor changes

File:

• trunk/Ohana/src/opihi/cmd.data/gridify.c (modified) (4 diffs)

trunk/Ohana/src/opihi/cmd.data/gridify.c

@@ -3 +3 @@
 int gridify (int argc, char **argv) {
 
-  int i, Nx, Ny, Xb, Yb;
+  int i, Nx, Ny, Xb, Yb, Normalize, N;
   float Xmin, Xmax, dX, Ymin, Ymax, dY;
   float *buf, *val, *x, *y, *z;
@@ -9 +9 @@
   Buffer *bf;
   Vector *vx, *vy, *vz;
+
+  Normalize = TRUE;
+  if ((N = get_argument (argc, argv, "-raw"))) {
+    remove_argument (N, &argc, argv);
+    Normalize = FALSE;
+  }
 
   if (argc != 11) {
@@ -66 +72 @@
     if (Xb >= Nx) continue;
     if (Yb >= Ny) continue;
-    val[Xb + Yb*Nx] = *z;
+    val[Xb + Yb*Nx] += *z;
     Nval[Xb + Yb*Nx]++;
   }
@@ -72 +78 @@
   buf = (float *) bf[0].matrix.buffer;
   for (i = 0; i < Nx*Ny; i++) {
-    if (Nval[i] == 0) {
-      buf[i] = 0;
-      continue;
+    if (Normalize) {
+      if (Nval[i] == 0) {
+        buf[i] = 0;
+        continue;
+      }
+      buf[i] = val[i] / Nval[i];
+    } else {
+      buf[i] = val[i];
     }
-    buf[i] = val[i] / Nval[i];
   }