Changeset 41379 for trunk/Ohana/src/opihi/cmd.astro

Timestamp:

Jun 22, 2020, 11:00:22 AM (6 years ago)

Author:

eugene

Message:

add -copy option to coords; add mjdtolst; add option to control major tick spacing; add -copy option to header; add usage info to avperiodomatch

Location:

trunk/Ohana/src/opihi/cmd.astro

Files:

• coords.c (modified) (2 diffs)
• init.c (modified) (2 diffs)
• jdtolst.c (modified) (2 diffs)

trunk/Ohana/src/opihi/cmd.astro/coords.c

@@ -12 +12 @@
   Vector *xvec = NULL;
   Vector *yvec = NULL;
+
+  if ((N = get_argument (argc, argv, "-copy"))) {
+    Buffer *src = NULL;
+
+    remove_argument (N, &argc, argv);
+    if ((src = SelectBuffer (argv[N], OLDBUFFER, TRUE)) == NULL) goto escape;
+
+    if (!GetCoords (&coords, &src[0].header)) {
+      gprint (GP_ERR, "error getting WCS elements from src buffer %s\n", argv[N]);
+      return (FALSE);
+    }
+
+    remove_argument (N, &argc, argv);
+
+    if (argc != 2) goto syntax;
+
+    if ((buf = SelectBuffer (argv[1], OLDBUFFER, TRUE)) == NULL) goto escape;
+    if (!PutCoords (&coords, &buf[0].header)) {
+      gprint (GP_ERR, "error getting WCS elements from src buffer %s\n", argv[1]);
+      return (FALSE);
+    }
+    return TRUE;
+  }
 
   MOSAIC = NULL;
@@ -116 +139 @@
  syntax:
   gprint (GP_ERR, "USAGE: coords [buffer] (-c R D) | (-p X Y)\n");
-  gprint (GP_ERR, "only one of -p or -c can be used\n");
+  gprint (GP_ERR, "USAGE: coords [tgtbuffer] (-copy srcbuffer)\n");
+  gprint (GP_ERR, "only one of -p or -c or -copy can be used\n");
   gprint (GP_ERR, " -p : from pixels to ra/dec\n");
   gprint (GP_ERR, " -c : from ra/dec to pixels\n");
   gprint (GP_ERR, " coordinates are in degrees\n");
+  gprint (GP_ERR, " -copy : copy coordinate WCS keywords from srcbuffer to tgtbuffer\n");
  escape:
   if (MOSAIC != NULL) free (MOSAIC);

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

@@ -34 +34 @@
 int imsub                   PROTO((int, char **));
 int jdtolst                 PROTO((int, char **));
+int mjdtolst                PROTO((int, char **));
 int medianmap               PROTO((int, char **));
 int galsectors              PROTO((int, char **));
@@ -105 +106 @@
   {1, "imsub",       imsub,        "subtract function"},
   {1, "jdtolst",     jdtolst,      "JD to LST conversion"},
+  {1, "mjdtolst",    mjdtolst,     "MJD to LST conversion"},
   {1, "medianmap",   medianmap,    "small median image"},
   {1, "mkgauss",     mkgauss,      "generate a 2-D gaussian centered in image"},

trunk/Ohana/src/opihi/cmd.astro/jdtolst.c

@@ -23 +23 @@
   char *long_str = get_variable ("LONGITUDE");
   if (!long_str) {
-    gprint (GP_ERR, "please define $LONGITUDE\n");
+    gprint (GP_ERR, "please define $LONGITUDE (decimal hours, west positive)\n");
     FREE (OutputName);
     return FALSE;
   }
   double longitude = atof (long_str);
-  if (VERBOSE) gprint (GP_ERR, "using longitude of %f\n", longitude);
+  if (VERBOSE) gprint (GP_ERR, "using longitude of %f (decimal hours, west positive)\n", longitude);
   free (long_str);
 
@@ -87 +87 @@
   gprint (GP_ERR, "  (JD) may be a number or a vector.\n");
   gprint (GP_ERR, "  if (JD) is a vector, (output) must be supplied and is used as the name of the output vector \n");
+  gprint (GP_ERR, "  $LONGITUDE must be set (decimal hours, west positive)\n");
   return (FALSE);
 }
+
+# define MJD_OFFSET 2400000.5
+
+int mjdtolst (int argc, char **argv) {
+
+  int N;
+  Vector *mjdvec, *outvec;
+
+  int VERBOSE = FALSE;
+  VERBOSE = FALSE;
+  if ((N = get_argument (argc, argv, "-v"))) {
+    VERBOSE = TRUE;
+    remove_argument (N, &argc, argv);
+  }
+
+  char *OutputName = NULL;
+  if ((N = get_argument (argc, argv, "-output"))) {
+    remove_argument (N, &argc, argv);
+    OutputName = strcreate (argv[N]);
+    remove_argument (N, &argc, argv);
+  }
+  if (argc != 2) goto syntax;
+
+  char *long_str = get_variable ("LONGITUDE");
+  if (!long_str) {
+    gprint (GP_ERR, "please define $LONGITUDE (decimal hours, west positive)\n");
+    FREE (OutputName);
+    return FALSE;
+  }
+  double longitude = atof (long_str);
+  if (VERBOSE) gprint (GP_ERR, "using longitude of %f (decimal hours, west positive)\n", longitude);
+  free (long_str);
+
+  if (ISNUM(argv[1][0])) {
+
+    char *endptr;
+    double mjd = strtod (argv[1], &endptr);
+    if (endptr == argv[1]) {
+      gprint (GP_ERR, "error interpretting %s as a number\n", argv[1]);
+      FREE (OutputName);
+      return FALSE;
+    }
+
+    double lst = ohana_lst (mjd + MJD_OFFSET, longitude);
+
+    if (OutputName) {
+      set_variable (OutputName, lst);
+      FREE (OutputName);
+      return (TRUE);
+    }
+    
+    gprint (GP_ERR, "lst: %f\n", lst);
+    return (TRUE);
+  }
+
+  if (!OutputName) {
+    gprint (GP_ERR, " if (MJD) is a vector, -output must be supplied\n");
+    return FALSE;
+  }
+
+  if ((mjdvec = SelectVector (argv[1], OLDVECTOR, TRUE)) == NULL) {
+    FREE (OutputName);
+    return FALSE;
+  }
+
+  if (mjdvec->type != OPIHI_FLT) {
+    gprint (GP_ERR, " (MJD) must be of type FLOAT\n");
+    FREE (OutputName);
+    return FALSE;
+  }
+
+  if ((outvec = SelectVector (OutputName, ANYVECTOR, TRUE)) == NULL) {
+    FREE (OutputName);
+    return FALSE;
+  }
+
+  ResetVector (outvec, OPIHI_FLT, mjdvec[0].Nelements);
+
+  for (int i = 0; i < mjdvec[0].Nelements; i++) {
+    double lst = ohana_lst (mjdvec->elements.Flt[i] + MJD_OFFSET, longitude);
+    outvec->elements.Flt[i] = lst;
+  }
+  FREE (OutputName);
+  return TRUE;
+
+ syntax:
+  gprint (GP_ERR, "USAGE: mjdtolst (MJD) [-output name]\n");
+  gprint (GP_ERR, "  (MJD) may be a number or a vector.\n");
+  gprint (GP_ERR, "  if (MJD) is a vector, (output) must be supplied and is used as the name of the output vector \n");
+  gprint (GP_ERR, "  $LONGITUDE must be set (decimal hours, west positive)\n");
+  return (FALSE);
+}