diff options
Diffstat (limited to 'src/analysis/modhead.c')
| -rw-r--r-- | src/analysis/modhead.c | 114 |
1 files changed, 114 insertions, 0 deletions
diff --git a/src/analysis/modhead.c b/src/analysis/modhead.c new file mode 100644 index 0000000..0ffc40b --- /dev/null +++ b/src/analysis/modhead.c @@ -0,0 +1,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); +} + |