Changeset 41164 for trunk/Ohana/src/opihi/cmd.data/interpolate_presort.c

Timestamp:

Nov 27, 2019, 11:22:02 AM (7 years ago)

Author:

eugene

Message:

extend cut function to apply alternate stats to the output vector; extend extract to allow source region to fall beyond bounds of input image; add commands antialias, kapamemory, match1d, mgaussdev, vsigmoid (fit a sigmoid), virls (irls on vector), vwtmean (weighted mean), vtransitions; extend medimage to report variance; fix error in rotate; update some tests; add vgauss -apply option

File:

• trunk/Ohana/src/opihi/cmd.data/interpolate_presort.c (modified) (3 diffs)

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

@@ -7 +7 @@
   double x0, dx, dy, y0;
   Vector *xout, *yout, *xin, *yin;
+
+  int N;
+  int FillEnds = FALSE;
+  if ((N = get_argument (argc, argv, "-fill-ends"))) {
+    FillEnds = TRUE;
+    remove_argument (N, &argc, argv);
+  }
 
   /** check basic syntax **/
@@ -36 +43 @@
   y0 = yin[0].elements.Flt[0];
   
-  /* in and out vectors are sorted */
-  j = 0;
+  // fill in the start with NANs
+  for (i = 0; (i < xout[0].Nelements) && (xout[0].elements.Flt[i] < xin[0].elements.Flt[0]); i++) {
+    yout[0].elements.Flt[i] = FillEnds ? yin[0].elements.Flt[0] : NAN;
+  }
 
   // every output pixel should get a value unless 
   // we are below the lowest in or above the highest in
-  for (i = 0, j = 0; (i < xout[0].Nelements) && (j < xin[0].Nelements - 1); ) {
+  for (j = 0; (i < xout[0].Nelements) && (j < xin[0].Nelements - 1); ) {
 
     yout[0].elements.Flt[i] = NAN;
@@ -66 +75 @@
 
   // fill in the rest with NANs
+  int NinLast = xin[0].Nelements - 1;
   while (i < yout[0].Nelements) {
-    yout[0].elements.Flt[i] = NAN;
+    yout[0].elements.Flt[i] = FillEnds ? yin[0].elements.Flt[NinLast] : NAN;
     i++;
   }
 
   return (TRUE);
-    
 }