tomwalters@0: /****************************************************************************
tomwalters@0: *
tomwalters@0: *   compilation:    cc -o xgram xgram.c -lX11
tomwalters@0: *
tomwalters@0: *   "profile" is an array with an element for each pixel across the window.
tomwalters@0: *
tomwalters@0: ****************************************************************************/
tomwalters@0: 
tomwalters@0: #include <stdio.h>
tomwalters@0: #include <math.h>
tomwalters@0: #include "x11coord.h"
tomwalters@0: 
tomwalters@0: char   DEFNAME[]=     "stdin";          /* window name */
tomwalters@0: char   *name =        DEFNAME;
tomwalters@0: 
tomwalters@0: /****************************************************************************
tomwalters@0: *   Defaults for plot
tomwalters@0: ****************************************************************************/
tomwalters@0: #define DEFn            128             /* number of data points on a plot line */
tomwalters@0: #define DEFln           FULLWIDTH>>1    /* length of plot line in pixels */
tomwalters@0: #define DEFl            1000            /* max number of lines to plot */
tomwalters@0: #define DEFs            0               /* number of plot lines offset from start of file */
tomwalters@0: #define DEFf            0.2             /* scale factor for plot y-values */
tomwalters@0: 
tomwalters@0: #define DEFdx           2               /* origin increments for successive lines, in pixels */
tomwalters@0: #define DEFdy           2
tomwalters@0: 
tomwalters@0: /****************************************************************************
tomwalters@0: *   Arguments
tomwalters@0: ****************************************************************************/
tomwalters@0: int    s=DEFs;              /* start position in file */
tomwalters@0: short  n=DEFn;              /* number of data points on a plot line */
tomwalters@0: short  l=DEFl;              /* max number of lines to plot */
tomwalters@0: short  ln=DEFln;            /* length of plot line in pixels */
tomwalters@0: short  dx=DEFdx,dy=DEFdy;   /* increments for successive lines, in pixels */
tomwalters@0: double f=DEFf;              /* scale factor for each y-value */
tomwalters@0: 
tomwalters@0: /****************************************************************************
tomwalters@0: *       main
tomwalters@0: ****************************************************************************/
tomwalters@0: main(argc, argv)
tomwalters@0: int   argc ;
tomwalters@0: char *argv[] ;
tomwalters@0: {
tomwalters@0:     FILE   *fp, *fopen();
tomwalters@0:     Window  w;
tomwalters@0:     short   p;
tomwalters@0:     int     i;
tomwalters@0: 
tomwalters@0: 
tomwalters@0:     while (--argc > 0 && **++argv == '-')
tomwalters@0: 	switch (*++*argv) {
tomwalters@0: 	    case 'n':   n       = atoi(++*argv); break;
tomwalters@0: 	    case 'l':   l       = atoi(++*argv); break;
tomwalters@0: 	    case 's':   s       = atoi(++*argv); break;
tomwalters@0: 	    case 'f':   f      *= atof(++*argv); break;
tomwalters@0: 	    case 'x':   dx      = atoi(++*argv); break;
tomwalters@0: 	    case 'y':   dy      = atoi(++*argv); break;
tomwalters@0: 	    case 'h':
tomwalters@0: 	    case 'H':
tomwalters@0: 	    default:    help();
tomwalters@0: 	}
tomwalters@0:     /* Open data file */
tomwalters@0:     if (argc) {
tomwalters@0: 	name = *argv;
tomwalters@0: 	if ((fp = fopen(name, "r")) == NULL) {
tomwalters@0: 	    fprintf(stderr,"can't open %s\n", *argv);
tomwalters@0: 	    exit(1);
tomwalters@0: 	}
tomwalters@0:     }
tomwalters@0:     else fp = stdin;
tomwalters@0:     /* Seek start of data in file */
tomwalters@0:     if (s > 0) {
tomwalters@0: 	s *= n;         /* convert number of lines to number of points */
tomwalters@0: 	for (i=0 ; i<s && fread(&p, sizeof(short), 1, fp) ; i++)
tomwalters@0: 	    ;
tomwalters@0: 	if (i<s) error("seek overshot end-of-file\n");
tomwalters@0:     }
tomwalters@0: 
tomwalters@0:     set_window_parameters(FULLXORG,FULLYORG, FULLWIDTH,FULLHEIGHT);
tomwalters@0:     set_bordered_box(1);
tomwalters@0: 
tomwalters@0:     /** set data so that points fill window ?? **/
tomwalters@0: 
tomwalters@0:     xopen();
tomwalters@0:     w = xcreate(name, ButtonPressMask | ExposureMask);
tomwalters@0:     xevent_monitor(w,fp);
tomwalters@0:     fclose(fp);
tomwalters@0:     XDestroyWindow(theDisplay,w);
tomwalters@0: }
tomwalters@0: 
tomwalters@0: 
tomwalters@0: /****************************************************************************
tomwalters@0: *       X11 Event Monitor.
tomwalters@0: ****************************************************************************/
tomwalters@0: xevent_monitor(w,fp)
tomwalters@0: Window w;
tomwalters@0: FILE *fp;
tomwalters@0: {
tomwalters@0:     XEvent event;
tomwalters@0: 
tomwalters@0:     for ( ; ; ) {
tomwalters@0: 	XNextEvent(theDisplay,&event);
tomwalters@0: 	switch (event.type) {
tomwalters@0: 	    case ButtonPress:
tomwalters@0: 		switch(event.xbutton.button) {
tomwalters@0: 		    case Button1:   return;     /* Left   */
tomwalters@0: 		    case Button2:   return;     /* Middle */
tomwalters@0: 		    case Button3:   return;     /* Right  */
tomwalters@0: 		    default:
tomwalters@0: 			fprintf(stderr,"button %d not used\n", event.xbutton.button);
tomwalters@0: 		}
tomwalters@0: 		break;
tomwalters@0: 	    case Expose:
tomwalters@0: 		if (event.xexpose.count == 0) {
tomwalters@0: 		    draw_box(w);
tomwalters@0: 		    draw_gram(w,FULLWIDTH,FULLHEIGHT,fp);
tomwalters@0: 		}
tomwalters@0: 		break;
tomwalters@0: 	    default:
tomwalters@0: 		fprintf(stderr,"event type %d not found\n", event.type);
tomwalters@0: 	}
tomwalters@0:     }
tomwalters@0: }
tomwalters@0: 
tomwalters@0: 
tomwalters@0: 
tomwalters@0: /****************************************************************************
tomwalters@0: *       Draw line (spectro)-gram with hidden-line removal.
tomwalters@0: ****************************************************************************/
tomwalters@0: #define     VISIBLE     0
tomwalters@0: #define     HIDDEN      1
tomwalters@0: #define     X0          10              /* origin of plot within window */
tomwalters@0: #define     Y0          10
tomwalters@0: #define     move(x,y)   X=x;Y=y
tomwalters@0: #define     cont(x,y)   XDrawLine(theDisplay, w, theGC, X+X0,height-Y-Y0, x+X0,height-y-Y0);move(x,y)
tomwalters@0: 
tomwalters@0: draw_gram(w,width,height,fp)
tomwalters@0: Window w;
tomwalters@0: int    width,height;    /* window size, in pixels */
tomwalters@0: FILE   *fp;             /* i/p data file */
tomwalters@0: {
tomwalters@0:     int    i;
tomwalters@0:     short  p, *profile;
tomwalters@0:     short  X, Y;        /* current active position */
tomwalters@0:     short  x, y, x0, y0, y1;
tomwalters@0:     short  xorg=0, yorg=0;
tomwalters@0:     float  index, res;
tomwalters@0:     short  rres;
tomwalters@0:     char   stat;
tomwalters@0: 
tomwalters@0:     profile = (short*) malloc(width * sizeof(short));
tomwalters@0:     for (i=0 ; i<width ; i++)
tomwalters@0: 	profile[i] = 0;
tomwalters@0: 
tomwalters@0:     res = (float)ln/n;
tomwalters@0:     rres = (int)(res+0.5);      /* rounded resolution */
tomwalters@0:     /* plot l lines, or until eof */
tomwalters@0:     while (l-->0 && fread(&p, sizeof(short), 1, fp)) {
tomwalters@0: 	/* draw vertical to start of line */
tomwalters@0: 	x = xorg;
tomwalters@0: 	y = yorg + p*f;
tomwalters@0: 	if (y >= (y0=profile[x])) {
tomwalters@0: 	    profile[x] = y;
tomwalters@0: 	    move(x,y0);
tomwalters@0: 	    cont(x,y);
tomwalters@0: 	    stat = VISIBLE;
tomwalters@0: 	}
tomwalters@0: 	else stat = HIDDEN;
tomwalters@0: 	y0 = y;
tomwalters@0: 
tomwalters@0: 	/* draw rest of line */
tomwalters@0: 	index = xorg;
tomwalters@0: 	for (i=1 ; i<n && fread(&p, sizeof(short), 1, fp) ; i++) {
tomwalters@0: 	    /* interpolate y-values for each pixel along line */
tomwalters@0: 	    for (x0=1, y1=yorg+p*f ; x0<=rres && x0+x<width ; x0++) {
tomwalters@0: 		y = (short)( ((float)x0/res) * (y1-y0) + y0 );
tomwalters@0: 		/* a Hidden line becomes Visible */
tomwalters@0: 		if (y >= profile[x+x0]) {
tomwalters@0: 		    profile[x+x0] = y;
tomwalters@0: 		    if (stat == HIDDEN) {
tomwalters@0: 			move(x+x0,y);
tomwalters@0: 			stat = VISIBLE;
tomwalters@0: 		    }
tomwalters@0: 		}
tomwalters@0: 		/* a Visible line becomes Hidden */
tomwalters@0: 		else {
tomwalters@0: 		    if (stat == VISIBLE) {
tomwalters@0: 			cont(x+x0,y);
tomwalters@0: 			stat = HIDDEN;
tomwalters@0: 		    }
tomwalters@0: 		}
tomwalters@0: 	    }
tomwalters@0: 	    y0 = y1;
tomwalters@0: 	    index += res;
tomwalters@0: 	    x = (int)(index+0.5);       /* rounded index */
tomwalters@0: 	    if (stat == VISIBLE)
tomwalters@0: 		cont(x,y);
tomwalters@0: 	}
tomwalters@0: 	cont(x,yorg);
tomwalters@0: 	xorg += dx;
tomwalters@0: 	yorg += dy;
tomwalters@0:     }
tomwalters@0: }
tomwalters@0: 
tomwalters@0: 
tomwalters@0: /****************************************************************************
tomwalters@0: *       Misc.
tomwalters@0: ****************************************************************************/
tomwalters@0: error(s)
tomwalters@0: {
tomwalters@0:     fprintf(stderr,"%s",s);
tomwalters@0:     exit(1);
tomwalters@0: }
tomwalters@0: 
tomwalters@0: help ()
tomwalters@0: {
tomwalters@0:     printf ("\nUsage: a) xgram  [options]  filename\n");
tomwalters@0:     printf ("       b) cat filename  |  xgram  [options]\n");
tomwalters@0:     printf ("    where filename is an input stream of 16-bit binary numbers\n");
tomwalters@0:     printf ("Options:\n");
tomwalters@0:     printf ("-n [int]      = Number of points per plotted line (default=%d).\n", DEFn);
tomwalters@0:     printf ("-l [int]      = Max number of lines plotted (default=%d).\n", DEFl);
tomwalters@0:     printf ("-s [int]      = Offset lines (of n points) from start of file (default=%d).\n", DEFs);
tomwalters@0:     printf ("-f [float]    = Scale factor (default=1).\n");
tomwalters@0:     printf ("-x-y[ints]    = dx/dy changes for successive lines (defaults %d,%d)\n", DEFdx,DEFdy);
tomwalters@0:     exit (0);
tomwalters@0: }