annotate saitools/napgraph.c @ 0:5242703e91d3 tip

Initial checkin for AIM92 aimR8.2 (last updated May 1997).

/*
    Copyright (c) Applied Psychology Unit, Medical Research Council. 1993
    ===========================================================================

    Permission to use, copy, modify, and distribute this software without fee 
    is hereby granted for research purposes, provided that this copyright 
    notice appears in all copies and in all supporting documentation, and that 
    the software is not redistributed for any fee (except for a nominal 
    shipping charge). Anyone wanting to incorporate all or part of this 
    software in a commercial product must obtain a license from the Medical 
    Research Council.

    The MRC makes no representations about the suitability of this 
    software for any purpose.  It is provided "as is" without express or 
    implied warranty.
 
    THE MRC DISCLAIMS ALL WARRANTIES WITH REGARD TO THIS SOFTWARE, INCLUDING 
    ALL IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS, IN NO EVENT SHALL 
    THE A.P.U. BE LIABLE FOR ANY SPECIAL, INDIRECT OR CONSEQUENTIAL DAMAGES 
    OR ANY DAMAGES WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, 
    WHETHER IN AN ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, 
    ARISING OUT OF OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS 
    SOFTWARE.
*/
/* napgraph.c 
*
* M. Akeroyd  11th June 1993. Revised Winter 1994.
*/





#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include "tip.h"


/* Function declarations */
/*---------------------------------------------------------------------------*/

void parsecommandline (int argc, char *argv[], char inputfn[], char outputfn[]);

int readheader_nap(FILE *inputfp);
void checkheader();
void writeheader(FILE *outputfp);
int fakeheader_sai();
void changeheader(int no_frames, int new_pwidth, int new_nwidth, int framewidth_samples_output);

void find_intervals_nap(int channel, int maximum_interval, int sample);

void writedata_output (FILE *outputfp, int samples_to_write);
FILE *open_file (char filefn[], FILE *dir_default, int streamtype);

void write_matrix_hori_time(int framewidth, FILE *outputfigurefp);
void write_matrix_hori_freq(int framewidth, FILE *outputfigurefp);




/* Data arrays: global */
/*---------------------------------------------------------------------------*/


short freqdata[MAX_CHANNELS];                /* input .nap: a single column */
short previousfreqdata[MAX_CHANNELS];
short previouspreviousfreqdata[MAX_CHANNELS];
short outputdata[MAX_DATA]; 
short outputfiguredata[MAX_DATA]; 
short outputgrounddata[MAX_DATA]; 

short interval[MAX_CHANNELS][MAX_DATA];      /* histogram */ 
short summation[MAX_CHANNELS];

short clip[MAX_CHANNELS];                    /* set to ON of clipping occured*/


/* other variables */
/*---------------------------------------------------------------------------*/
struct Peak peak[MAX_PEAKS];
struct Peak peak_nap_right[MAX_CHANNELS];
struct Peak peak_nap_left[MAX_CHANNELS];


/* variables read from header */
/*---------------------------------------------------------------------------*/

char header[MAX_LINES_HEADER][MAX_LINE_LENGTH];
int header_lines;

int no_frames;
int frameheight;             /* number of channels */
int framewidth_samples;      /* pwidth + nwidth * samplerate */
int frameshift_samples;      /* frstep_aid * samplerate */
int pwidth;                  /* in msecs */
int nwidth;                  /* in msecs: NEGATIVE */
int width_win;               /* pixels */
int height_win;              /* pixels */
long samplerate;             /* samples per sec */
int mincf;                   /* Hz */
int maxcf;                   /* Hz */

int no_columns;
int nap_height;
int greyscaleflag = OFF;


/* misc */
/*---------------------------------------------------------------------------*/

char progname[MAX_STRING_LENGTH];
char outputfigurefn[MAX_STRING_LENGTH];


int verboseflag = OFF;               /* -v */
int removeflag = OFF;                /* -r */
double removevalue = 0.1;
int widthflag = OFF;                 /* -w */     /* if OFF, use DEFAULT */
int total_width;                     /* -w */     /* total width of new .sai */
int asciiflag = OFF;
int asciireverseflag = OFF;
int summationflag = OFF;
int oppositearchflag = OFF;          /* -oparch : see saigraph.c for info */



/* ...................           Main           .............................*/
/* ..........................................................................*/
/* ..........................................................................*/



void main (int argc, char *argv[])
{
  int sample;
  int n;

  int column, channel;
  int framewidth_samples_input,
      framewidth_samples_output;
  int new_pwidth, 
         new_nwidth;
  int header_bytes = 0;
 
  char inputfn[MAX_STRING_LENGTH], 
       outputbasefn[MAX_STRING_LENGTH];

  FILE *inputfp, *outputfigurefp;


  strcpy(progname, argv[0]);
  strcpy(inputfn, "");
  strcpy(outputbasefn, "");
  strcpy(outputfigurefn, "");

  /* parse command line */
  parsecommandline(argc, argv, inputfn, outputbasefn);

  if (strcmp(outputbasefn, "") == 0) 
    strcpy(outputfigurefn, "");
  else {
    strcpy(outputfigurefn, outputbasefn);
/*    strcat(outputfigurefn, OUTPUT_EXT);*/
  }


  /* open files */
  /* default directions are:
   * input = stdin 
   * figure = stdout
   */

  inputfp = open_file(inputfn, stdin, READ);
  outputfigurefp = open_file(outputfigurefn, stdout, WRITE);

  /* read Header */
  header_bytes = readheader_nap(inputfp);

  /* do some error-checking on the header, and set the fread pointer 
   * to the right place */
  checkheader();

  /* fake the .sai header, but first copy the important variables:
   *   no of columns (no_frames)
   *   height (frameheight)
   */
  no_columns = no_frames;
  nap_height = frameheight;
  header_lines = fakeheader_sai();

  /* reset the variables */
  framewidth_samples_input = framewidth_samples;
  if (widthflag == OFF) {
    new_pwidth = + PWIDTH;
    new_nwidth = - NWIDTH;}
  else {
    new_pwidth = + NEW_PWIDTH;
    new_nwidth = - (new_pwidth + total_width);}

  framewidth_samples_output = (int) ((int) (abs(new_pwidth) + abs(new_nwidth)) * samplerate ) / 1000;
  
  if (framewidth_samples_output >= MAX_DATA) {
    fprintf(stderr, "%s: new frame is  too wide at %i: only alllowed %i samples\n", progname, framewidth_samples_output, MAX_DATA);
    exit(-1); }

  
  /* change header */
  changeheader(1, new_pwidth, new_nwidth, framewidth_samples_output);

  /* write header */
  if (asciiflag == OFF) 
    writeheader(outputfigurefp);
  else {
    /* begining of ascii header. The rest depends on the matrix format ... */
    fprintf(outputfigurefp, "# napgraph output.  %s\n", inputfn);
    fprintf(outputfigurefp, "# \n");
    fprintf(outputfigurefp, "# samplerate=%i  pwidth=%d  nwidth=%d  columns=%i ", samplerate, new_pwidth, new_nwidth, no_columns);
    fprintf(outputfigurefp, " mincf=%i  maxcf=%i  channels=%i \n", mincf, maxcf, nap_height);
    fprintf(outputfigurefp, "# \n");
  }

/* --------------------------------------------------------------------------*/


  if  (verboseflag==ON) { 
    fprintf(stderr, " channels "); 
    fflush(stderr);}

  
  /* clear arrays */
  
  for (channel=0; channel<MAX_CHANNELS; channel++) {
    previousfreqdata[channel] = previouspreviousfreqdata[channel] = 0;
    peak_nap_left[channel].sample = 0;
    peak_nap_right[channel].sample = 0;
    for(sample=0; sample<MAX_DATA; sample++)
      interval[channel][sample] = 0;
    clip[channel] = OFF;
  }


  /* Main loop */


/* .....  */

  for (column=1; column<=no_columns; column++) {

    if ((column % 100) == 1 ) { 
      if (verboseflag==ON) {
	fprintf(stderr, " %i ", column); 
	fflush(stderr);}
    }

    /* load a column's worth of data  */
    fread (freqdata, 2, (size_t) nap_height, inputfp);
    
    /* move everything UP one byte.
     * Don't know why ... */
    for (channel=nap_height; channel >=1; channel--)
      freqdata[channel] = freqdata[channel-1];

    /* This next is a simple input check: if any numbers < 0 , then say so */
    for(channel=0; channel<nap_height; channel++)
      if (freqdata[channel] < 0 )  {
	fprintf(stderr, "%s: something's gone wrong: the data is negative.\n", progname);
	exit(-1); }


    /* the important bits ... */
    for (channel=1; channel <=nap_height; channel++)
      find_intervals_nap(channel, framewidth_samples_output, column);
    

  } /* column */


/* --------------------------------------------------------------------------*/
/* --------------------------------------------------------------------------*/


  fflush(stdout);

  /* If required, sum the values */
  if (summationflag==ON) {
    /* clear ... */
    for (channel = 1; channel <= frameheight; channel ++)
	summation[channel] = 0;
    /* set ... */
    for (channel = 1; channel <= frameheight; channel ++)
      for (sample = 0; sample < framewidth_samples_output; sample++)
	summation[channel] += interval[channel][sample];
  }


  /* write output */

  if (asciiflag == OFF) {
    for (channel = 1; channel <= frameheight; channel ++){
      for (sample = 0; sample < framewidth_samples_output; sample++)
	outputfiguredata[sample] = interval[channel][sample];
      writedata_output(outputfigurefp, framewidth_samples_output);
    }
  }  
  else {
    frameheight = nap_height;
    if (asciireverseflag == OFF) 
      write_matrix_hori_freq(framewidth_samples_output, outputfigurefp);
    else
      write_matrix_hori_time(framewidth_samples_output, outputfigurefp);
}    


/* --------------------------------------------------------------------------*/
/* --------------------------------------------------------------------------*/


  /* Tidy up and exit */
  fprintf(stderr, "\n");
  fclose(inputfp); 
  fclose(outputfigurefp); 
 
  if ( ferror(inputfp) != 0) {
    fprintf(stderr, " %s : error closing input file.\n", progname);
    exit(-1);}

  if ( ferror(outputfigurefp) != 0) {
    fprintf(stderr, " %s : error closing figure file.\n", progname);
    exit(-1);}

  exit(0);

} /* Main */








/* ..........................................................................*/
/* ..........................................................................*/
/* ..........................................................................*/
/* ..........................................................................*/










void parsecommandline(int argc, char *argv[], char inputfn[], char outputbasefn[])
{
  int x=1, helpflag = OFF;
  int group;

  while (x < argc){
    if      (!strcmp(argv[x], "-o")) {strcpy(outputbasefn, argv[x+1]); x+=2;}
    else if (!strcmp(argv[x], "-output")) {strcpy(outputbasefn, argv[x+1]); x+=2;}
    else if (!strcmp(argv[x], "-i")) {strcpy(inputfn, argv[x+1]); x+=2;}
    else if (!strcmp(argv[x], "-input")) {strcpy(inputfn, argv[x+1]); x+=2;}
    else if (!strcmp(argv[x], "-help")) {helpflag = ON; x+=1;}
    else if (!strcmp(argv[x], "-h")) {helpflag = ON; x+=1;}
    else if (!strcmp(argv[x], "-verbose")) {verboseflag = ON; x+=1;}
    else if (!strcmp(argv[x], "-v")) {verboseflag = ON; x+=1;}
    else if (!strcmp(argv[x], "-w")) {widthflag = ON; total_width=atof(argv[x+1]); x+=2;}
    else if (!strcmp(argv[x], "-width")) {widthflag = ON; total_width=atof(argv[x+1]); x+=2;}
    else if (!strcmp(argv[x], "-v")) {verboseflag = ON; x+=1;}
    else if (!strcmp(argv[x], "-r")) { removeflag = ON; removevalue = atof(argv[x+1]); x+=2;}
    else if (!strcmp(argv[x], "-a")) { asciiflag = ON;x+=1;}
    else if (!strcmp(argv[x], "-ar")) { asciiflag = ON; asciireverseflag=ON; x+=1;}
    else if (!strcmp(argv[x], "-s")) { summationflag=ON; x+=1;}
    else if (!strcmp(argv[x], "-g")) { greyscaleflag=ON; x+=1;}
    else if (!strcmp(argv[x], "-grey")) { greyscaleflag=ON; x+=1;}
    else if (!strcmp(argv[x], "-gray")) { greyscaleflag=ON; x+=1;}
    else if (!strcmp(argv[x], "-oparch")) {oppositearchflag = ON; x+=1;}
    else {fprintf(stderr, "%s: unknown option %s\n", progname, argv[x]);
	  exit(-1);}
  }
  if (helpflag == ON)
    {
      fprintf(stderr, "\n");
      fprintf(stderr, "                         %s\n", progname);
      fprintf(stderr, " -----------------------------------------------------------\n");
      fprintf(stderr, " Works out a histogram of peak intervals, from a .nap\n");
      fprintf(stderr, " Default output is a .sai file: if -a or -ar specified,\n");
      fprintf(stderr, " then output is a ASCII matrix.\n");
      fprintf(stderr, "\n");
      fprintf(stderr, " The -oparch option is REQUIRED if reading a Sun .sai files \n");
      fprintf(stderr, " on a DEC cpu or vice-versa.\n");
      fprintf(stderr, " Add .nap to input filename, .sai to output filename.\n");
      fprintf(stderr, " -----------------------------------------------------------\n");
      fprintf(stderr, " -i <.nap file>      input file             default = stdin\n");
      fprintf(stderr, " -o <.sai file>      output file            default = stdout\n");
      fprintf(stderr, " -oparch          assume input was generated on an opposing architecture cpu\n");
      fprintf(stderr, "\n");
      fprintf(stderr, " -a                  ASCII output: rows = channels\n");
      fprintf(stderr, " -ar                 ASCII output: rows = time intervals\n");
      fprintf(stderr, " -s                  include summations (ASCII only)\n");
      fprintf(stderr, "\n");
      fprintf(stderr, " -w <int>            width of graph.sai (msecs).\n");
      fprintf(stderr, " -r <fraction>       remove any peaks F(n) < F(n+1) & F(n-1) * fract \n");
      fprintf(stderr, "                                            default = %.3f\n", removevalue);
      fprintf(stderr, "\n");
      fprintf(stderr, " -v                  verbose output (to stderr)\n");
      fprintf(stderr, " -g                  Change .sai format to 'grayscale'.\n");
      fprintf(stderr, "\n\n");
      exit(-1);}

}