Changeset 40792

Timestamp:

Jun 12, 2019, 3:49:40 PM (7 years ago)

Author:

eugene

Message:

ignore inf/nan values in cut

File:

• branches/eam_branches/ohana.20190329/src/opihi/cmd.data/cut.c (modified) (2 diffs)

branches/eam_branches/ohana.20190329/src/opihi/cmd.data/cut.c

@@ -59 +59 @@
       /* accumulate values */
       Vin = (float *)(buf[0].matrix.buffer) + sy*Nx + sx + i; 
+      int Npix = 0;
       for (j = 0; j < ny; j++, Vin += Nx) {
+        if (!isfinite(*Vin)) continue;
         Vbuf[j] = *Vin;
+        Npix ++;
+      }
+      if (Npix == 0) {
+        value = NAN;
+        goto set_yvalue;
       }
       /* apply stat of choice */
       if (Mode == MEDIAN) {
-        fsort (Vbuf, ny);
-        value = Vbuf[(int)(0.5*ny)];
+        fsort (Vbuf, Npix);
+        value = Vbuf[(int)(0.5*Npix)];
       } else {
         value = 0;
-        for (j = 0; j < ny; j++) {
+        for (j = 0; j < Npix; j++) {
           value += Vbuf[j];
         }
-        if (Mode == MEAN) { value /= ny; }
+        if (Mode == MEAN) { value /= Npix; }
       }
+    set_yvalue:
       yvec[0].elements.Flt[i] = value;
     }
@@ -91 +99 @@
       /* accumulate values */
       Vin = (float *)(buf[0].matrix.buffer) + (sy + i)*Nx + sx; 
+      int Npix = 0;
       for (j = 0; j < nx; j++, Vin ++) {
+        if (!isfinite(*Vin)) continue;
         Vbuf[j] = *Vin;
+        Npix ++;
+      }
+      if (Npix == 0) {
+        value = NAN;
+        goto set_xvalue;
       }
       /* apply stat of choice */
       if (Mode == MEDIAN) {
-        fsort (Vbuf, nx);
-        value = Vbuf[(int)(0.5*nx)];
+        fsort (Vbuf, Npix);
+        value = Vbuf[(int)(0.5*Npix)];
       } else {
         value = 0;
-        for (j = 0; j < nx; j++) {
+        for (j = 0; j < Npix; j++) {
           value += Vbuf[j];
         }
-        if (Mode == MEAN) { value /= nx; }
+        if (Mode == MEAN) { value /= Npix; }
       }
+    set_xvalue:
       yvec[0].elements.Flt[i] = value;
     }