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
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
|
/*****************************************************************************
* Johns Hopkins University
* Center For Astrophysical Sciences
* FUSE
*****************************************************************************
*
* Synopsis: modhead file.fit[ext] KEYWORD NEW_VALUE
*
* Description: Returns or modifies the value of a header keyword.
* If NEW_VALUE is not specified, current value of KEYWORD is
* returned. If KEYWORD is not found and NEW_VALUE is given,
* KEYWORD is added to header. If an extension beyond the
* primary HDU is requested, file name must be given in
* parentheses: "file.fit[3]" or an alternate construction
* used: file.fit+3 Note that the primary HDU is extension 0.
*
* History: 05/15/05 tc started work
* 08/30/05 wvd Allow multi-word strings.
*
****************************************************************************/
#include <string.h>
#include <stdio.h>
#include "fitsio.h"
int main(int argc, char *argv[])
{
fitsfile *fptr; /* FITS file pointer, defined in fitsio.h */
char card[FLEN_CARD], newcard[FLEN_CARD];
char oldvalue[FLEN_VALUE], comment[FLEN_COMMENT];
int status = 0; /* CFITSIO status value MUST be initialized to zero! */
int i, iomode, keytype;
if (argc == 3)
iomode = READONLY;
else if (argc >= 4)
iomode = READWRITE;
else if (argc < 3) {
printf("Usage: modhead filename[ext] keyword new_value\n");
printf("\n");
printf("Write or modify the value of a header keyword.\n");
printf("If 'newvalue' is not specified, then just print \n");
printf("the current value. \n");
printf("\n");
printf("Examples: \n");
printf(" modhead file.fits dec - list the DEC keyword \n");
printf(" modhead file.fits dec 30.0 - set DEC = 30.0 \n");
printf("\n");
printf(" modhead \"file.fits[2]\" OBSTIME - list OBSTIME from extention 2 \n");
printf(" modhead file.fits+2 OBSTIME - list OBSTIME from extention 2\n");
return(0);
}
if (!fits_open_file(&fptr, argv[1], iomode, &status))
{
if (fits_read_card(fptr,argv[2], card, &status))
{
printf("Keyword does not exist\n");
card[0] = '\0';
comment[0] = '\0';
status = 0; /* reset status after error */
}
else
printf("%s\n",card);
if (argc >= 4) /* write or overwrite the keyword */
{
/* check if this is a protected keyword that must not be changed */
if (*card && fits_get_keyclass(card) == TYP_STRUC_KEY)
{
printf("Protected keyword cannot be modified.\n");
}
else
{
/* get the comment string */
if (*card)fits_parse_value(card, oldvalue, comment, &status);
/* construct template for new keyword */
strcpy(newcard, argv[2]); /* copy keyword name */
strcat(newcard, " = "); /* '=' value delimiter */
if (argc == 4) strcat(newcard, argv[3]); /* new value */
else {
strcat(newcard, "'");
for (i = 3; i < argc-1; i++) {
strcat(newcard, argv[i]); /* new value */
strcat(newcard, " ");
}
strcat(newcard, argv[i]); /* new value */
strcat(newcard, "'");
}
if (*comment) {
strcat(newcard, " / "); /* comment delimiter */
strcat(newcard, comment); /* append the comment */
}
/* reformat the keyword string to conform to FITS rules */
fits_parse_template(newcard, card, &keytype, &status);
/* overwrite the keyword with the new value */
fits_update_card(fptr, argv[2], card, &status);
printf("Keyword has been changed to:\n");
printf("%s\n",card);
}
} /* if argc == 4 */
fits_close_file(fptr, &status);
} /* open_file */
/* if error occured, print out error message */
if (status) fits_report_error(stderr, status);
return(status);
}
|
