Changeset 9726 for trunk/Ohana/src/opihi/cmd.data

Timestamp:

Oct 24, 2006, 11:58:18 AM (20 years ago)

Author:

eugene

Message:

added mkrgb, fixed shift bug

Location:

trunk/Ohana/src/opihi/cmd.data

Files:

• Makefile (modified) (1 diff)
• init.c (modified) (2 diffs)
• mkrgb.c (added)
• shift.c (modified) (4 diffs)

trunk/Ohana/src/opihi/cmd.data/Makefile

@@ -64 +64 @@
 $(SDIR)/load.$(ARCH).o          \
 $(SDIR)/lookup.$(ARCH).o        \
+$(SDIR)/mkrgb.$(ARCH).o \
 $(SDIR)/mcreate.$(ARCH).o       \
 $(SDIR)/medacc.$(ARCH).o        \

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

@@ -49 +49 @@
 int load             PROTO((int, char **));
 int lookup           PROTO((int, char **));
+int mkrgb            PROTO((int, char **));
 int mcreate          PROTO((int, char **));
 int medacc           PROTO((int, char **));
@@ -160 +161 @@
   {"load",         load,             "load an SAOimage style overlay"},
   {"lookup",       lookup,           "convert vector via lookup table (vector pair)"},
+  {"mkrgb",        mkrgb,            "convert 3 images to rgb jpeg"},
   {"mcreate",      mcreate,          "create a matrix"},
   {"medacc",       medacc,           "accumulate vector values in another vector"},

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

@@ -7 +7 @@
   float *Vin, *Vot;
   double Dx, Dy, fdx, fdy;
-  Buffer *buf;
+  Buffer *in, *out;
 
   ROLL = FALSE;
@@ -15 +15 @@
   }
 
-  if (argc != 4) {
-    gprint (GP_ERR, "USAGE: shift <buffer> dx dy\n");
+  if (argc != 5) {
+    gprint (GP_ERR, "USAGE: shift (input) (output) dx dy\n");
     return (FALSE);
   }
 
-  if ((buf = SelectBuffer (argv[1], OLDBUFFER, TRUE)) == NULL) return (FALSE);
+  // define the input buffer and examine the shift
+  if ((in  = SelectBuffer (argv[1], OLDBUFFER, TRUE)) == NULL) return (FALSE);
 
-  Dx = -atof (argv[2]);
-  Dy = -atof (argv[3]);
+  Dx = atof (argv[3]);
+  Dy = atof (argv[4]);
 
   dx = Dx;
@@ -29 +30 @@
   fdx = Dx - dx;
   fdy = Dy - dy;
-  if (fdx < -0.001) {dx -= 1; fdx += 1;}
-  if (fdy < -0.001) {dy -= 1; fdy += 1;}
+  if (fdx < -0.000001) {dx -= 1; fdx += 1;}
+  if (fdy < -0.000001) {dy -= 1; fdy += 1;}
+  // we always specify a positive fractional shift
+  // the above choice defines a minimum fractional shift of 1e-5
 
-  nx = buf[0].matrix.Naxis[0];
-  ny = buf[0].matrix.Naxis[1];
+  nx = in[0].matrix.Naxis[0];
+  ny = in[0].matrix.Naxis[1];
 
-  if ((dx > buf[0].matrix.Naxis[0]) || (dy > buf[0].matrix.Naxis[1])) {
+  if ((dx > nx) || (dy > ny)) {
     gprint (GP_ERR, "shifting data out of image\n");
     return (FALSE);
   }
   
-  if (dx < 0) {
-    DXin = 0;
-    DXot = -dx;
-  } else {
-    DXin = dx;
-    DXot = 0;
+  // define the output buffer
+  if ((out = SelectBuffer (argv[2], ANYBUFFER, TRUE)) == NULL) return (FALSE);
+
+  gfits_free_matrix (&out[0].matrix);
+  gfits_free_header (&out[0].header);
+  CreateBuffer (out, nx, ny, -32, 0.0, 1.0);
+
+  DXin = (dx < 0) ? -dx : 0;
+  DXot = (dx < 0) ?   0 : dx;
+  DYin = (dy < 0) ? -dy : 0;
+  DYot = (dy < 0) ?   0 : dy;
+  
+  for (j = 0; j < ny - abs(dy); j++) {
+    Vin = (float *)(in[0].matrix.buffer)  + (j + DYin)*nx + DXin;  
+    Vot = (float *)(out[0].matrix.buffer) + (j + DYot)*nx + DXot; 
+    for (i = 0; i < nx - abs(dx); i++, Vin++, Vot++) {
+      *Vot = *Vin;
+    }
+    // fill in the exposed x-border with 0.0
+    Vot = (dx > 0) ? 
+      (float *)(out[0].matrix.buffer) + (j + DYot)*nx : 
+      (float *)(out[0].matrix.buffer) + (j + DYot)*nx + nx - abs(dx);
+    for (i = 0; i < abs(dx); i++, Vot++) {
+      *Vot = 0.0;
+    }   
   }
-  if (dy < 0) {
-    DYin = 0;
-    DYot = -dy;
-  } else {
-    DYin = dy;
-    DYot = 0;
-  }
-  
-  if ((dy > 0) || ((dy == 0) && (dx > 0))) {
-    for (j = 0; j < ny - fabs(dy); j++) {
-      Vin = (float *)(buf[0].matrix.buffer) + (j + DYin)*nx + DXin;  
-      Vot = (float *)(buf[0].matrix.buffer) + (j + DYot)*nx + DXot; 
-      for (i = 0; i < nx - fabs(dx); i++, Vin++, Vot++) {
-        *Vot = *Vin;
-      }
-      if (dx < 0) 
-        Vot = (float *)(buf[0].matrix.buffer) + (j + DYot)*nx; 
-      for (i = 0; i < fabs(dx); i++, Vot++) {
-        *Vot = 0.0;
-      } 
-    }
-    Vot = (float *)(buf[0].matrix.buffer) + nx*(ny - abs(dy));
-    for (j = 0; j < nx * fabs(dy); j++, Vot++) {
-      *Vot = 0.0;
-    }   
-  } else {
-    for (j = ny - fabs(dy) - 1; j >= 0; j--) {
-      Vin = (float *)(buf[0].matrix.buffer) + (j + DYin)*nx + DXin + nx - abs(dx) - 1; 
-      Vot = (float *)(buf[0].matrix.buffer) + (j + DYot)*nx + DXot + nx - abs(dx) - 1; 
-      for (i = 0; i < nx - fabs(dx); i++, Vin--, Vot--) {
-        *Vot = *Vin;
-      }
-      if (dx < 0) 
-        Vot = (float *)(buf[0].matrix.buffer) + (j + DYot)*nx; 
-      for (i = 0; i < fabs(dx); i++, Vot++) {
-        *Vot = 0.0;
-      } 
-    }
-    Vot = (float *)(buf[0].matrix.buffer);
-    for (j = 0; j < nx * fabs(dy); j++, Vot++) {
-      *Vot = 0.0;
-    }
-  } 
 
+  // fill in the exposed y-border with 0.0
+  Vot = (dy > 0) ? 
+    (float *)(out[0].matrix.buffer) :
+    (float *)(out[0].matrix.buffer) + (ny - abs(dy))*nx;
+
+  for (j = 0; j < nx * abs(dy); j++, Vot++) {
+    *Vot = 0.0;
+  }   
+
+  // apply the fractional shift 
   gprint (GP_ERR, "%f %f\n", fdx, fdy);
   if ((fdx > 0) || (fdy > 0)) {
@@ -101 +90 @@
     f11 =    fdx *   fdy;
 
-    Vin = (float *)buf[0].matrix.buffer;
-    Vot = (float *)buf[0].matrix.buffer;
+    Vin = (float *)out[0].matrix.buffer;
+    Vot = (float *)out[0].matrix.buffer;
     for (j = 0; j < ny - 1; j++) {
       for (i = 0; i < nx - 1; i++, Vin++, Vot++) {
         value  = Vin[   0] * f00;
         value += Vin[   1] * f01;
-        value += Vin[ny  ] * f10;
-        value += Vin[ny+1] * f11;
+        value += Vin[nx  ] * f10;
+        value += Vin[nx+1] * f11;
         *Vot = value;
       }