ana_ehdr.c

Go to the documentation of this file.

/*****************************************************************************/
#include "tpcclibConfig.h"
/*****************************************************************************/
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <string.h>
#include <time.h>
/*****************************************************************************/
#include "libtpcmisc.h"
#include "libtpcimgio.h"
/*****************************************************************************/
 
/*****************************************************************************/
static char *info[] = {
  "Edits the information in Analyze 7.5 header file.",
  "Note that frame times and isotope information are not stored in Analyze,",
  "header but in SIF.",
  " ",
  "Usage: @P [Options] headerfile fieldname value ",
  " ",
  "Options:",
  " -stdoptions", // List standard options like --help, -v, etc
  " ",
  "Example:",
  "   @P s2345dy1.hdr patient_id 'PD013' ",
  " ",
  "See also: ana_lhdr, nii_lhdr, iftlist, eframe, ecat2ana",
  " ",
  "Keywords: image, Analyze, header, IFT",
  0};
/*****************************************************************************/
 
/*****************************************************************************/
/* Turn on the globbing of the command line, since it is disabled by default in
   mingw-w64 (_dowildcard=0); in MinGW32 define _CRT_glob instead, if necessary;
   In Unix&Linux wildcard command line processing is enabled by default. */
/*
#undef _CRT_glob
#define _CRT_glob -1
*/
int _dowildcard = -1;
/*****************************************************************************/
 
/*****************************************************************************/
int main(int argc, char **argv)
{
  int          ai, help=0, version=0, verbose=1;
  int          ret;
  char         hdrfile[FILENAME_MAX], *cptr;
  char         fieldname[31], fieldvalue[256];
  ANALYZE_DSR  dsr;
 
 
  /*
   *  Get arguments
   */
  if(argc==1) {tpcPrintUsage(argv[0], info, stderr); return(1);}
  hdrfile[0]=fieldname[0]=fieldvalue[0]=(char)0;
  /* Options */
  for(ai=1; ai<argc; ai++) if(*argv[ai]=='-') { /* options */
    if(tpcProcessStdOptions(argv[ai], &help, &version, &verbose)==0) continue;
    fprintf(stderr, "Error: invalid option '%s'.\n", argv[ai]);
    return(1);
  } else break;
  
  /* Print help or version? */
  if(help==2) {tpcHtmlUsage(argv[0], info, ""); return(0);}
  if(help) {tpcPrintUsage(argv[0], info, stdout); return(0);}
  if(version) {tpcPrintBuild(argv[0], stdout); return(0);}
 
  /* Process other arguments, starting from the first non-option */
  for(; ai<argc; ai++) {
    if(!hdrfile[0]) {strcpy(hdrfile, argv[ai]); continue;}
    if(!fieldname[0]) {strcpy(fieldname, argv[ai]); continue;}
    if(strcmp(argv[ai], ":=")==0) continue;
    if(!fieldvalue[0]) {strcpy(fieldvalue, argv[ai]); continue;}
    fprintf(stderr, "Error: invalid argument '%s'.\n", argv[ai]);
    return(1);
  }
 
  /* Is something missing? */
  if(!fieldvalue[0]) {
    fprintf(stderr, "Error: missing command-line argument.\n"); 
    return(1);
  }
 
  /* In verbose mode print arguments and options */
  if(verbose>1) {
    printf("hdrfile := %s\n", hdrfile);
    printf("fieldname := %s\n", fieldname);
    printf("fieldvalue := %s\n", fieldvalue);
  }
  if(verbose>2) ANALYZE_TEST=verbose-2; else ANALYZE_TEST=0;
 
 
  /*
   *  Read Analyze header file
   */
  /* Check that filename contains .hdr extension */
  cptr=filenameGetExtension(hdrfile);
  if(cptr==NULL || strcmp(cptr, ".hdr")!=0) strcat(hdrfile, ".hdr");
  /* Read header file */
  if(verbose>1) printf("reading %s\n", hdrfile);
  ret=anaReadHeader(hdrfile, &dsr);
  if(ret) {
    fprintf(stderr, "Error: cannot read %s.\n", hdrfile);
    if(verbose>1) fprintf(stderr, "  ret := %d\n", ret);
    return(2);
  }
 
 
  /*
   *  Change header field
   */
  ret=anaEditHeader(&dsr, fieldname, fieldvalue);
  if(ret==1) {
    fprintf(stderr, "Error: invalid field name: %s\n", fieldname);
    return(1);
  } else if(ret==2) {
    fprintf(stderr, "Error: invalid field value: %s\n", fieldvalue);
    return(1);
  }
  if(verbose>8) anaPrintHeader(&dsr, stdout);
 
 
  /*
   *  Write header
   */
  if(verbose>1) printf("writing %s\n", hdrfile);
  ret=anaWriteHeader(hdrfile, &dsr);
  if(ret) {
    fprintf(stderr, "Error: cannot write %s (%d).\n", hdrfile, ret);
    return(11);
  }
 
  return(0);
}
/*****************************************************************************/
 
/*****************************************************************************/