aboutsummaryrefslogtreecommitdiff
path: root/vendor/voclient/voapps/voiminfo.c
diff options
context:
space:
mode:
Diffstat (limited to 'vendor/voclient/voapps/voiminfo.c')
-rw-r--r--vendor/voclient/voapps/voiminfo.c493
1 files changed, 493 insertions, 0 deletions
diff --git a/vendor/voclient/voapps/voiminfo.c b/vendor/voclient/voapps/voiminfo.c
new file mode 100644
index 00000000..fe2baa8d
--- /dev/null
+++ b/vendor/voclient/voapps/voiminfo.c
@@ -0,0 +1,493 @@
+/*
+ * VOIMINFO -- Get the WCS information for a FITS image.
+ *
+ * Usage:
+ * voiminfo [<otps>] image.fits
+ *
+ * @file voiminfo.c
+ * @author Mike Fitzpatrick
+ * @date 11/03/12
+ *
+ * @brief Get the WCS information for a FITS image.
+ */
+
+#include <stdio.h>
+#include <unistd.h>
+#include <stdlib.h>
+#include <string.h>
+#include <ctype.h>
+#include <math.h>
+
+#include "votParse.h" /* keep these in order! */
+#include "voApps.h"
+#include "fitsio.h"
+
+
+
+#define MAX_IMAGES 20480 /* max images to process */
+#define SZ_RESBUF 819200
+
+#define OPT_ALL 0001
+#define OPT_BOX 0002
+#define OPT_CORNERS 0004
+#define OPT_EXTNS 0010
+#define OPT_FRAME 0020
+#define OPT_NAXES 0040
+
+
+static int do_return = 0; /* return result? */
+
+/* Global task declarations. These should all be defined as 'static' to
+ * avoid namespace collisions.
+ */
+static int do_all = 0; /* print global values */
+static int do_sex = 0; /* print sexagesimal values */
+
+static int do_box = 0; /* print box coordinates */
+static int do_corners = 0; /* print image corners */
+static int do_extns = 0; /* print extension values */
+static int do_info = 0; /* print image info */
+static int do_naxis = 0; /* print NAXIS values */
+
+static int debug = 0; /* debug flag */
+static int verbose = 1; /* verbose flag */
+
+
+/* Task specific option declarations. Task options are declared using the
+ * getopt_long(3) syntax.
+ */
+int voiminfo (int argc, char **argv, size_t *len, void **result);
+
+static Task self = { "voiminfo", voiminfo, 0, 0, 0 };
+static char *opts = "%:habcdvnseio:r";
+static struct option long_opts[] = {
+ { "test", 1, 0, '%'}, /* --test is std */
+ { "help", 2, 0, 'h'}, /* --help is std */
+ { "return", 2, 0, 'r'}, /* --return is std */
+ { "debug", 2, 0, 'd'}, /* debug flag */
+ { "verbose", 2, 0, 'v'}, /* verbose flag */
+
+ { "all", 2, 0, 'a'}, /* print global value */
+ { "box", 2, 0, 'b'}, /* print box values */
+ { "corners", 2, 0, 'c'}, /* print image corner */
+ { "extns", 2, 0, 'e'}, /* print each extn */
+ { "info", 2, 0, 'i'}, /* print image info */
+ { "naxes", 2, 0, 'n'}, /* print NEXIS values */
+ { "output", 1, 0, 'o'}, /* output filename */
+ { "sex", 2, 0, 's'}, /* sexagesimal values */
+ { NULL, 0, 0, 0 }
+};
+
+
+/**
+ * Private procedures.
+ */
+static void Usage (void);
+static void Tests (char *input);
+static char *fmt (double pos, int is_ra);
+static char *fmt_naxis (char *imname, ImInfo *im, int do_all);
+static char *fmt_box (char *imname, ImInfo *im, int do_all);
+static char *fmt_corners (char *imname, ImInfo *im, int do_all);
+
+extern int vos_getURL (char *url, char *name);
+
+
+
+/**
+ * Application entry point.
+ */
+int
+voiminfo (int argc, char **argv, size_t *reslen, void **result)
+{
+ /* These declarations are required for the VOApps param interface.
+ */
+ char **pargv, optval[SZ_FNAME], resbuf[SZ_RESBUF];
+ char imname[SZ_LINE];
+
+
+ /* These declarations are specific to the task.
+ */
+ char *oname = NULL, *imlist[MAX_IMAGES], *nimlist[MAX_IMAGES];
+ int i, ch = 0, status = OK, pos = 0, nfiles = 0, narg = 0;;
+ size_t maxlen = SZ_RESBUF;
+ FILE *fd = (FILE *) NULL;
+ ImInfo *im;
+
+
+ /* Initialize result object whether we return an object or not.
+ */
+ *reslen = 0;
+ *result = NULL;
+ memset (imlist, 0, MAX_IMAGES);
+ memset (nimlist, 0, MAX_IMAGES);
+
+
+ /* Parse the argument list. The use of vo_paramInit() is required to
+ * rewrite the argv[] strings in a way vo_paramNext() can be used to
+ * parse them. The programmatic interface allows "param=value" to
+ * be passed in, but the getopt_long() interface requires these to
+ * be written as "--param=value" so they are not confused with
+ * positional parameters (i.e. any param w/out a leading '-').
+ */
+ pargv = vo_paramInit (argc, argv, opts, long_opts);
+ while ((ch = vo_paramNext(opts,long_opts,argc,pargv,optval,&pos)) != 0) {
+ if (ch > 0) {
+ /* If the 'ch' value is > 0 we are parsing a single letter
+ * flag as defined in the 'opts string.
+ */
+ switch (ch) {
+ case '%': Tests (optval); return (self.nfail);
+ case 'h': Usage (); return (OK);
+ case 'r': do_return=1; break;
+ case 'd': debug++; break;
+ case 'v': verbose++; break;
+
+ case 'a': do_all++; break;
+ case 'b': do_box++; break;
+ case 'c': do_corners++; break;
+ case 'e': do_extns++; break;
+ case 'i': do_info++; break;
+ case 's': do_sex++; break;
+ case 'n': do_naxis++; break;
+ case 'o': oname = strdup (optval); break;
+ default:
+ fprintf (stderr, "Invalid option '%s'\n", optval);
+ return (1);
+ }
+
+ } else if (ch == PARG_ERR) {
+ return (ERR);
+
+ } else {
+ /* This code processes the positional arguments. The 'optval'
+ * string contains the value but since this string is
+ * overwritten w/ each arch we need to make a copy (and must
+ * remember to free it later.
+ */
+ imlist[nfiles++] = strdup (optval);
+ narg++;
+ }
+
+ if (narg > MAX_IMAGES) {
+ fprintf (stderr, "ERROR: Too many images to process\n");
+ return (1);
+ }
+ }
+
+
+ /* Sanity checks.
+ */
+ if (imlist[0] == NULL || strcmp (imlist[0], "-") == 0) {
+ free (imlist[0]);
+ imlist[0] = strdup ("stdin");
+ nfiles = 1;
+ }
+ if (do_all && strcasecmp (imlist[0], "stdin") == 0) {
+ fprintf (stderr, "Error: Option not supported with standard input\n");
+ return (ERR);
+ }
+ if (oname == NULL) oname = strdup ("stdout");
+ if (strcmp (oname, "-") == 0) { free (oname), oname = strdup ("stdout"); }
+
+
+ /* Open the output file.
+ */
+ if (strcmp (oname, "stdout") != 0) {
+ if ((fd = fopen (oname, "w+")) == (FILE *) NULL) {
+ fprintf (stderr, "Error: Cannot open output file '%s'\n", oname);
+ return (ERR);
+ }
+ } else
+ fd = stdout;
+
+ /* Default to printing all info if we're not given a specific option.
+ if (!do_box && !do_corners && !do_info && !do_naxis)
+ do_info = 1;
+ */
+
+
+ /**
+ * Main body of task
+ */
+ for (i=0; i < nfiles; i++) {
+
+ memset (imname, 0, SZ_LINE);
+ if (strncmp ("http://", imlist[i], 7) == 0) {
+ if (access ("/tmp/voiminfo.fits", F_OK) == 0)
+ unlink ("/tmp/voiminfo.fits");
+ strcpy (imname, "/tmp/voiminfo.fits");
+ if (vos_getURL (imlist[i], imname) <= 0)
+ continue;
+ } else
+ strcpy (imname, imlist[i]);
+
+ if ((im = vot_imageInfo (imname, do_all)) ) {
+ memset (resbuf, 0, SZ_RESBUF);
+ if (do_info) {
+ vot_printImageInfo (fd, im);
+ } else if (do_naxis) {
+ strcpy (resbuf, fmt_naxis (imname, im, do_all));
+ } else if (do_box) {
+ strcpy (resbuf, fmt_box (imname, im, do_all));
+ } else if (do_corners) {
+ strcpy (resbuf, fmt_corners (imname, im, do_all));
+ } else {
+ sprintf (resbuf, "%s\t%s\t%s\t%s\n", imlist[i],
+ fmt (im->frame.cx,1), fmt (im->frame.cy,0),
+ fmt (im->frame.radius,0));
+ }
+
+ if (do_return)
+ vo_appendResultFromString (resbuf, reslen, result, &maxlen);
+ else
+ fprintf (fd, "%s", resbuf);
+
+ vot_freeImageInfo (im);
+
+ if (strncmp ("http://", imlist[i], 7) == 0)
+ unlink (imname);
+ }
+ }
+
+
+ /* Clean up. Rememebr to free whatever pointers were created when
+ * parsing arguments.
+ */
+ for (i=0; i < MAX_IMAGES; i++) {
+ if (imlist[i])
+ free (imlist[i]);
+ if (nimlist[i])
+ free (nimlist[i]);
+ else
+ break;
+ }
+ if (oname) free (oname);
+ if (fd != stdout)
+ fclose (fd);
+
+ vo_paramFree (argc, pargv);
+
+ return (status); /* status must be OK or ERR (i.e. 0 or 1) */
+}
+
+
+/**
+ * USAGE -- Print task help summary.
+ */
+static void
+Usage (void)
+{
+ fprintf (stderr, "\n Usage:\n\t"
+ "voiminfo [<opts>] votable.xml\n\n"
+ " where\n"
+ " -%%,--test run unit tests\n"
+ " -h,--help this message\n"
+ "\n"
+ " -a,--all print all extns\n"
+ " -i,--info print image info\n"
+ " -b,--box print box values\n"
+ " -c,--corners print image corners\n"
+ " -n,--naxes print NAXIS values\n"
+ " -o,--output output filename\n"
+ " -s,--sex sexagesimal values\n"
+ "\n"
+ " -o,--output=<file> output file\n"
+ "\n"
+ "\n"
+ " Examples:\n\n"
+ " 1) Print the size of each image extension in an MEF file:\n\n"
+ " %% voiminfo -n mef.fits\n"
+ "\n"
+ " 2) Print the LL and UR positions of an image as sexigesimal:\n\n"
+ " %% voiminfo -b -s image.fits\n"
+ "\n"
+ " 3) Print the corner positions of an image as decimal degrees:\n\n"
+ " %% voiminfo -c image.fits\n"
+ "\n"
+ " 4) Print the box values for an entire mosaic MEF file:\n\n"
+ " %% voiminfo -b mef.fits\n"
+ "\n"
+ );
+}
+
+
+/**
+ * Tests -- Task unit tests.
+ */
+static void
+Tests (char *input)
+{
+ Task *task = &self;
+ char *mef = "data/mef.fits";
+ char *sif = "data/sif.fits";
+ char *opts[] = { "-a", "-i", "-b", "-c", "-n", "-s", NULL };
+ register int i;
+
+
+ vo_taskTest (task, "--help", NULL);
+
+ if (access (input, F_OK) == 0) {
+ vo_taskTest (task, "-n", input, NULL); // Ex 1
+ vo_taskTest (task, "-b", "-s", input, NULL); // Ex 2
+ vo_taskTest (task, "-c", input, NULL); // Ex 3
+ vo_taskTest (task, "-b", input, NULL); // Ex 4
+ } else
+ fprintf (stderr, "Warning: input file '%s' not found...\n", input);
+
+ if (access (sif, F_OK) == 0)
+ for (i=0; opts[i]; i++)
+ vo_taskTest (task, opts[i], sif, NULL);
+
+ if (access (mef, F_OK) == 0)
+ for (i=0; opts[i]; i++)
+ vo_taskTest (task, opts[i], mef, NULL);
+
+ vo_taskTestReport (self);
+}
+
+
+
+/***************************************************************************/
+/**** Private Procedures ****/
+/***************************************************************************/
+
+/**
+ * FMT -- Output formatter.
+ */
+static char *
+fmt (double pos, int is_ra)
+{
+ static char strbuf[10][SZ_LINE];
+ static int bufnum = 0;
+
+
+ bufnum = (bufnum + 1) % 10;
+ memset (strbuf[bufnum], 0, SZ_LINE);
+
+ pos = (pos < 0.0 ? -pos : pos);
+
+ if (do_sex) {
+ int d, m;
+ float s, frac;
+ char sign = (pos < 0.0 ? '-' : ' ');
+
+ pos = (is_ra ? (pos / 15.0) : pos);
+ d = (int) pos;
+ frac = (pos - d);
+ m = frac * 60.0;
+ s = ((frac * 60.0) - m) * 60.0;
+
+ sprintf (strbuf[bufnum], "%c%02d:%02d:%04.1f", sign, d, m, s);
+
+ } else
+ sprintf (strbuf[bufnum], "%f", pos);
+
+ return (strbuf[bufnum]);
+}
+
+
+/**
+ * FMT_NAXIS -- Print the size of pixels of an image.
+ */
+static char *
+fmt_naxis (char *imname, ImInfo *im, int do_all)
+{
+ static char line[SZ_RESBUF];
+
+ memset (line, 0, SZ_RESBUF);
+ if (im->nextend && !do_all) {
+ register int i = 0;
+ char buf[SZ_RESBUF];
+
+ for (i=0; i < im->nextend; i++) {
+ memset (buf, 0, SZ_RESBUF);
+ if (im->frame.naxis == 3)
+ sprintf (buf, "%s[%d]\t%d %d %d\n", imname, i,
+ im->extns[i].naxes[0],
+ im->extns[i].naxes[1],
+ im->extns[i].naxes[2]);
+ else
+ sprintf (buf, "%s[%d]\t%d %d\n", imname, i,
+ im->extns[i].naxes[0], im->extns[i].naxes[1]);
+ strcat (line, buf);
+ }
+
+ } else {
+ if (im->frame.naxis == 3)
+ sprintf (line, "%s\t%d %d %d\n", imname,
+ im->frame.naxes[0], im->frame.naxes[1], im->frame.naxes[2]);
+ else
+ sprintf (line, "%s\t%d %d\n", imname,
+ im->frame.naxes[0], im->frame.naxes[1]);
+ }
+
+ return (line);
+}
+
+
+/**
+ * FMT_BOX -- Print the corners of a box and rotation angle.
+ */
+static char *
+fmt_box (char *imname, ImInfo *im, int do_all)
+{
+ static char line[SZ_RESBUF];
+
+ memset (line, 0, SZ_RESBUF);
+ if (im->nextend && !do_all) {
+ register int i = 0;
+ char buf[SZ_RESBUF];
+
+ for (i=0; i < im->nextend; i++) {
+ memset (buf, 0, SZ_RESBUF);
+ sprintf (buf, "%s[%d]\t%s %s\t%s %s\t%s\n", imname, i,
+ fmt (im->extns[i].lx,1), fmt (im->extns[i].ly,0),
+ fmt (im->extns[i].ux,1), fmt (im->extns[i].uy,0),
+ fmt (im->extns[i].rotang,0));
+ strcat (line, buf);
+ }
+
+ } else {
+ sprintf (line, "%s\t%s %s\t%s %s\t%s\n", imname,
+ fmt (im->frame.lx,1), fmt (im->frame.ly,0),
+ fmt (im->frame.ux,1), fmt (im->frame.uy,0),
+ fmt (im->frame.rotang,0));
+ }
+
+ return (line);
+}
+
+
+/**
+ * FMT_CORNERS -- Print the explicit corners of an image footprint.
+ */
+static char *
+fmt_corners (char *imname, ImInfo *im, int do_all)
+{
+ static char line[SZ_RESBUF];
+
+ memset (line, 0, SZ_RESBUF);
+ if (im->nextend && !do_all) {
+ register int i = 0;
+ char buf[SZ_RESBUF];
+
+ for (i=0; i < im->nextend; i++) {
+ memset (buf, 0, SZ_RESBUF);
+
+ sprintf (buf, "%s\t%s %s\t%s %s\t%s %s\t%s %s\n", imname,
+ fmt (im->extns[0].xc[0],1), fmt (im->extns[0].yc[0],0),
+ fmt (im->extns[0].xc[1],1), fmt (im->extns[0].yc[1],0),
+ fmt (im->extns[0].xc[2],1), fmt (im->extns[0].yc[2],0),
+ fmt (im->extns[0].xc[3],1), fmt (im->extns[0].yc[3],0));
+ strcat (line, buf);
+ }
+
+ } else {
+ sprintf (line, "%s\t%s %s\t%s %s\t%s %s\t%s %s\n", imname,
+ fmt (im->frame.xc[0],1), fmt (im->frame.yc[0],0),
+ fmt (im->frame.xc[1],1), fmt (im->frame.yc[1],0),
+ fmt (im->frame.xc[2],1), fmt (im->frame.yc[2],0),
+ fmt (im->frame.xc[3],1), fmt (im->frame.yc[3],0));
+ }
+ return (line);
+}
generated by cgit (git 2.34.1) at 2025-09-21 01:49:46 -0400