1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
|
/*
# This file is part of the Astrometry.net suite.
# Licensed under a 3-clause BSD style license - see LICENSE
*/
#include <unistd.h>
#include "log.h"
#include "errors.h"
#include "sip.h"
#include "sip_qfits.h"
#include "fitsioutils.h"
#include "boilerplate.h"
#include "anqfits.h"
const char* OPTIONS = "hv";
static void print_help(char* progname) {
BOILERPLATE_HELP_HEADER(stdout);
printf("\nUsage: %s <input-filename> <extension> <output-filename>\n"
"\n\n", progname);
}
int main(int argc, char** args) {
char* outfn = NULL;
char* infn = NULL;
int ext;
int c;
int loglvl = LOG_MSG;
qfits_header* hdr;
sip_t sip;
while ((c = getopt(argc, args, OPTIONS)) != -1) {
switch (c) {
case '?':
case 'h':
print_help(args[0]);
exit(0);
case 'v':
loglvl++;
break;
}
}
log_init(loglvl);
if (optind != argc-3) {
print_help(args[0]);
exit(-1);
}
infn = args[optind+0];
ext = atoi(args[optind+1]);
outfn = args[optind+2];
logmsg("Reading extension %i from file \"%s\"\n", ext, infn);
hdr = anqfits_get_header2(infn, ext);
if (!hdr) {
ERROR("Failed to read header from extension %i of file \"%s\"\n", ext, infn);
exit(-1);
}
if (!sip_read_header(hdr, &sip)) {
ERROR("Failed to read SIP header.\n");
exit(-1);
}
if (sip.a_order > 0) {
logmsg("Got SIP header.\n");
} else {
logmsg("Got TAN header.\n");
}
logmsg("Writing to file \"%s\"\n", outfn);
if (sip_write_to_file(&sip, outfn)) {
ERROR("Failed to write SIP header to file \"%s\".\n", outfn);
exit(-1);
}
qfits_header_destroy(hdr);
return 0;
}
|
