diff options
Diffstat (limited to 'vendor/voclient/libvoclient/examples/siap4.c')
| -rw-r--r-- | vendor/voclient/libvoclient/examples/siap4.c | 146 |
1 files changed, 146 insertions, 0 deletions
diff --git a/vendor/voclient/libvoclient/examples/siap4.c b/vendor/voclient/libvoclient/examples/siap4.c new file mode 100644 index 00000000..5d5b8def --- /dev/null +++ b/vendor/voclient/libvoclient/examples/siap4.c @@ -0,0 +1,146 @@ +/************************************************************************ + * Call a SIAP search service and download resulting datasets. + * + * Usage: siap4 ra dec size [format [maximages [serviceURL]]] + * + * Or call with no args for the built-in unit test. + * + * M. Fitzpatrick, NOAO, June 2006 + */ + + +#include <stdio.h> +#include <stdlib.h> +#include "VOClient.h" + +double ra = 12.0; /* default values */ +double dec = 0.0; +double size = 0.5; +char *format = "ALL"; +int maximages= 5; + +char *service1 = "http://www-gsss.stsci.edu/gscvo/DSS2.jsp"; +char *service2 = "http://skyview.gsfc.nasa.gov/cgi-bin/vo/sia.pl?"; + +char *server = "6200:localhost"; + +static void callSiapService (char *url, double ra, double dec, double size, + char *format, int maximages); + + +int main (int argc, char *argv[]) +{ + char *service = service2; + + + /* Process command line arguments. + */ + if (argc <= 1) { + /* Use builtin defaults. */ + } else if (argc >= 3) { + int arg = 1; + + /* Look for a server specification. */ + if (strncmp (argv[arg], "-ds",3) == 0) + server = argv[++arg]; + ra = atof (argv[arg++]); + dec = atof (argv[arg++]); + size = atof (argv[arg++]); + if (arg < argc) format = argv[arg++]; + if (arg < argc) maximages = atoi (argv[arg++]); + if (arg < argc) service = argv[arg++]; + + } else { + fprintf (stderr, + "Usage: siap4 [-ds server] ra dec size [fmt [maximg [siapURL]]]\n"); + exit(1); + } + + /* Now call the Siap Service and summarize the results. + */ + callSiapService (service, ra, dec, size, format, maximages); + + return (0); +} + + +/* Simple test routine to call a Siap search service and summarize results. + */ +static void +callSiapService (char *service_url, double ra, double dec, double size, + char *format, int maximages) +{ + char *voc_opts = NULL; + char *acref = NULL; + char fname[SZ_FNAME]; + int i, nattr = 0, nrec = 0; + + DAL siap; /* DAL Connection */ + Query query; /* query handle */ + QResponse qr; /* query response handle */ + QRecord rec; /* result record handle */ + QRAttribute v; /* dataset attribute */ + + + /* Initialize the VOClient code. Error messages are printed by the + * interface so we just quit if there is a problem. + */ + if (voc_initVOClient (voc_opts) == ERR) + return; + + /* Get a new connection to the named service. + */ + siap = voc_openSiapConnection (service_url); /* open a connection */ + + /* Form a query. Here we'll use the one search size we're given for + * both the RA,DEC sizes, and specify a null format. + */ + query = voc_getSiapQuery (siap, ra, dec, size, size, format); + + printf ("Executing Query:\n %s\n\n", + voc_getQueryString (query, SIAP_CONN, 0)); + qr = voc_executeQuery (query); /* execute the query */ + + + /* Summarize response. + */ + if ((nrec = voc_getRecordCount (qr)) <= 0) { + fprintf (stderr, "no records matched"); + exit (1); + + } else { + rec = voc_getRecord (qr, 0); + nattr = (rec != ERR) ? voc_getAttrCount (rec) : 0; + + printf ("# returns %d records containing %d attributes each\n#\n", + nrec, nattr); + } + + + /* Download the first 'maximages' images. + */ + printf ("Downloading images:\n"); + + for (i=0; i < nrec && i < maximages; i++) { + rec = voc_getRecord (qr, i); /* get a row in the table */ + + v = voc_getAttribute (rec, "AccessReference"); + + if (v <= 0) + continue; + + bzero (fname, SZ_FNAME); + sprintf (fname, "dataset.%04d", i); + printf ("Downloading: %s\n", (acref = voc_stringValue (v)) ); + + if ( voc_getDataset (rec, acref, fname) == OK ) + printf ("Downloaded %s\n", fname); + else + printf ("Download failed\n"); + } + + voc_closeConnection (siap); /* close the siap connection */ + voc_closeVOClient (1); /* clean up and shutdown */ + + return; +} |