mcd.c

speech signal process tools

/*
 * This material contains unpublished, proprietary software of 
 * Entropic Research Laboratory, Inc. Any reproduction, distribution, 
 * or publication of this work must be authorized in writing by Entropic 
 * Research Laboratory, Inc., and must bear the notice: 
 *
 *    "Copyright (c) 1986-1990  Entropic Speech, Inc. 
 *    "Copyright (c) 1990-1991  Entropic Research Laboratory, Inc. 
 *                   All rights reserved"
 *
 * The copyright notice above does not evidence any actual or intended 
 * publication of this source code.     
 *
 * Written by:  Rod Johnson
 * Checked by:
 * Revised by: Ken Nelson
 *
 * Brief description: Sun filter for plot data in Unix GPS format.
 *
 */

static char *sccs_id = "@(#)mcd.c	3.5	6/14/91	ESI/ERL";
int    debug_level = 0;

/*+-----------------------------------------------------------------------+
|									|
|   mcd - Sun filter for plot data in Unix GPS format			|
|									|
|   Rodney W. Johnson, Entropic Speech, Inc.				|
|									|
+----------------------------------------------------------------------*/
#define VERSION "3.1"
#define DATE "10/7/87"

#include <stdio.h>
#include <signal.h>
#include <esps/unix.h>
#include <esps/esps.h>
#include <esps/spsassert.h>
#include <esps/constants.h>

#define SYNTAX	    USAGE("mcd [-r region][-u][-x level] [file.gps]")
#define EOF_ERR	    {Fprintf(stderr, "%s: unexpected end of file.\n", \
		    ProgName); exit(1);}
#define READ	    fread((char *) &input_word, sizeof input_word, 1, infile)
#define SAVE	    (void) fwrite((char *) &input_word, sizeof input_word, 1, tfile)
#define GET1(code, length) (READ && (code = (input_word>>12) & 0xf, \
		    length = input_word & 0xfff, 1))
#define GET2(coord) {if (READ) coord = input_word; else EOF_ERR}
#define GET3(m,n)   {if (READ) {m = (input_word>>8) & 0xff; \
		    n = input_word & 0xff;} else EOF_ERR}
#define GETCHRS(ch1,ch2) {if (READ) {ch1 = input_chr[0]; \
		    ch2 = input_chr[1];} else EOF_ERR}

#define L_ROUND(n) ((n)>=0 ? (long)((n)+0.5) : -(long)(-(n)+0.5))

#define LINES	    0
#define ARC	    3
#define TEXT	    2
#define HARDWARE    4
#define COMMENT	    15

#define BUFSIZE	    2046
#define TXTSIZE	    8185
#define TEXT_OFFSET   20	/* space to save for text when using range */
extern int text_offset;

double          sin(), cos();

int 
init_dev(), erase(), label(), line(), circle(),
arc(), move(), cont(), point(), linemod(),
space(), closepl();
int             text();
int             plotline();
void            do_plot();
void            start_plot();

char           *ProgName = "mcd";
char           *Version = VERSION;
char           *Date = DATE;

short           input_word;
char           *input_chr = (char *) &input_word;
int             c;
extern          optind;
extern char    *optarg;
int             e_flag = NO;
int             u_flag = NO;
int             region = 0;
FILE           *infile;
char           *inname;
FILE           *tfile;
char           *tmpname;
int             code;
int             length;
static int      xdata[BUFSIZE], ydata[BUFSIZE];
static char     chdata[TXTSIZE];
static int      region_x[] = {-32767, -20479, -8191, 4096, 16384};
static int      region_y[] = {-32767, -19660, -6553, 6553, 19660};
int             x, y;
int             bundle;
int             unused;
int             size;
int             rotation;
double          cosrot, sinrot;
int             i;
int             x0 = 32767, y0 = 32767, x1 = -32767, y1 = -32767;
int		R_flag = NO;
int		P_flag = NO;
char		*play_options = NULL;
char 		*title = "ESPS Plot";
char		*icon_file = NULL;
int		Argc;
char		**Argv;
void		convert_options();


main(argc, argv)
	int             argc;
	char          **argv;
{
	int		onintr();
	(void)signal(SIGTERM, onintr);      /* set up to catch SIGTERM */
    	(void)signal(SIGINT, onintr);       /* set up to catch SIGINT */

	Argc = argc;
	Argv = argv;	
        convert_options(Argv, Argc);
	init_device();

	while ((c = getopt(Argc, Argv, "T:r:ux:t:RPi:p:")) != EOF)
		switch (c) {
		case 'T':
			break;
		case 'e':
			e_flag = YES;
			break;
		case 'r':
			region = atoi(optarg);
			u_flag = NO;
			break;
		case 'u':
			region = 0;
			u_flag = YES;
			break;
		case 'x':
			debug_level = atoi(optarg);
			break;
		case 't':
			title = optarg;
			break;
		case 'R':
			R_flag = YES;
			text_offset = TEXT_OFFSET;
			break;
		case 'i':
			icon_file = optarg;
			break;
		case 'P':
			P_flag = YES; /* enable play from range mode */
			break;
		case 'p':
			play_options = optarg;
			break;
		default:
			SYNTAX;
			break;
		}

	if (optind >= Argc) {
		inname = "<stdin>";
		infile = stdin;
	} else {
		inname = argv[optind++];
		if (strcmp(inname, "-") == 0) {
			inname = "<stdin>";
			infile = stdin;
		} else
			TRYOPEN(ProgName, inname, "r", infile);
	}

	if (!e_flag)
		erase();


	/*
	 * read the gps stream into a temporary file
	 */


	tmpname = mktemp("/usr/tmp/mcdXXXXXX");
	TRYOPEN(ProgName, tmpname, "w", tfile);
	while (GET1(code, length)) {
		SAVE;
		if (debug_level >= 1)
			Fprintf(stderr, "Writing temporary file.\n");
		switch (code) {
		case LINES:
			if (debug_level >= 1)
				Fprintf(stderr, "lines\t(%d words)\n", length);
			spsassert(length % 2 == 0,
			 "can't draw line with odd number of coordinates.");
			spsassert(length >= 4,
				  "cant't draw line with no coordinates.");

			for (i = 0; i < (length - 2) / 2; i++) {
				GET2(x);
				SAVE;
				GET2(y);
				SAVE;
				if (x < x0)
					x0 = x;
				if (y < y0)
					y0 = y;
				if (x > x1)
					x1 = x;
				if (y > y1)
					y1 = y;
			}

			GET3(bundle, unused);
			SAVE;

			break;
		case ARC:
			if (debug_level >= 1)
				Fprintf(stderr, "arc\t(%d words)\n", length);
			spsassert(length % 2 == 0,
			  "can't draw arc with odd number of coordinates.");
			spsassert(length != 0,
				  "cant't draw arc with no coordinates.");
			for (i = 0; i < (length - 2) / 2; i++) {
				GET2(xdata[i]);
				SAVE;
				GET2(ydata[i]);
				SAVE;
			}

			GET3(bundle, unused);
			SAVE;

			break;
		case TEXT:
			if (debug_level >= 1)
				Fprintf(stderr, "text\t(%d words)\n", length);
			spsassert(length >= 6, "incomplete text command");
			GET2(x);
			SAVE;
			GET2(y);
			SAVE;
			GET3(bundle, unused);
			SAVE;
			GET3(size, rotation);
			SAVE;

			for (i = 0; i < 2 * (length - 5); i += 2) {
				GETCHRS(chdata[i], chdata[i + 1]);
				SAVE;
			}
			chdata[i] = '\0';

			cosrot = cos(PI * rotation / 128.0);
			sinrot = sin(PI * rotation / 128.0);
			{
				int xhi = x, xlo = x, yhi = y, ylo = y;
				double rgt = 
				 4.0 * size * (strlen(chdata) - 0.5), 
				 lft = 2.0 * size, top = 2.5 * size, 
				 bot = 4.5 * size;

				if (cosrot > 0) {
					xhi += rgt * cosrot;
					xlo -= lft * cosrot;
					yhi += top * cosrot;
					ylo -= bot * cosrot;
				} else {
					xhi -= lft * cosrot;
					xlo += rgt * cosrot;
					yhi -= bot * cosrot;
					ylo += top * cosrot;
				}
				if (sinrot > 0) {
					xhi += bot * sinrot;
					xlo -= top * sinrot;
					yhi += rgt * sinrot;
					ylo -= lft * sinrot;
				} else {
					xhi -= top * sinrot;
					xlo += bot * sinrot;
					yhi -= lft * sinrot;
					ylo += rgt * sinrot;
				}
				if (xlo < x0)
					x0 = xlo;
				if (ylo < y0)
					y0 = ylo;
				if (xhi > x1)
					x1 = xhi;
				if (yhi > y1)
					y1 = yhi;
			}

			break;
		case HARDWARE:
			if (debug_level >= 1)
				Fprintf(stderr, 
				 "hardware text\t(%d words)\n", length);
			spsassert(length >= 4, "incomplete hardware text command");
			GET2(x);
			SAVE;
			GET2(y);
			SAVE;
			for (i = 0; i < 2 * (length - 3); i += 2) {
				GETCHRS(chdata[i], chdata[i + 1]);
				SAVE;
			}
			chdata[i] = '\0';

			break;
		case COMMENT:
			if (debug_level >= 1)
				Fprintf(stderr,"comment\t(%d words)\n", length);
			for (i = 0; i < length - 1; i++) {
				GET2(unused);
				SAVE;
			}
			break;
		default:
			Fprintf(stderr, 
			 "%s: unimplemented plotting command (code %d).\n",
			 ProgName, code);
			exit(1);
			break;
		}
	}

	(void)fclose(tfile);

	if (u_flag) {
		space(-32767, -32767, 32767, 32767);
		if (debug_level >= 1)
			Fprintf(stderr, "space(%d, %d, %d, %d)\n",
				-32767, -32767, 32767, 32767);
	} else if (region) {
		x = region_x[(region - 1) % 5];
		y = region_y[(region - 1) / 5];
		space(x, y, x + 16383, y + 13107);
		if (debug_level >= 1)
			Fprintf(stderr, 
			 "space(%d, %d, %d, %d)\n", x, y, x + 16383, y + 13107);
	}

/* call Sunviews stuff to start the action */
	start_plot();

/* if you want to use this file as a GPS front-end to some other
   device or window system, all you have to do, is to have the start_plot
   routine call back the do_plot function.
*/
}

void
do_plot()
{
	TRYOPEN(ProgName, tmpname, "r", tfile);
	infile = tfile;
	if (debug_level >= 1)
		Fprintf(stderr, "Reading temporary file.\n");
	space(x0, y0, x1, y1);
	if (debug_level >= 1)
		Fprintf(stderr, "space(%d, %d, %d, %d)\n", x0, y0, x1, y1);
	while (GET1(code, length)) {
		switch (code) {
		case LINES:
			if (debug_level >= 1)
				Fprintf(stderr, "lines\t(%d words)\n", length);
			spsassert(length % 2 == 0,
			 "can't draw line with odd number of coordinates.");
			spsassert(length >= 4,
				  "cant't draw line with no coordinates.");

			for (i = 0; i < (length - 2) / 2; i++) {
				GET2(xdata[i])
					GET2(ydata[i])
					if (debug_level >= 2)
					Fprintf(stderr, 
					 "\t%d %d\n", xdata[i], ydata[i]);
			}

			GET3(bundle, unused)
				if (debug_level >= 2)
				Fprintf(stderr, 
				 "\tbundle %d, unused %d.\n", bundle, unused);

			move(xdata[0], ydata[0]);
			for (i = 1; i < (length - 2) / 2; i++)
				cont(xdata[i], ydata[i]);

			break;
		case ARC:
			if (debug_level >= 1)
				Fprintf(stderr, "arc\t(%d words)\n", length);
			spsassert(length % 2 == 0,
			  "can't draw arc with odd number of coordinates.");
			spsassert(length != 0,
				  "cant't draw arc with no coordinates.");
			for (i = 0; i < (length - 2) / 2; i++) {
				GET2(xdata[i])
					GET2(ydata[i])
					if (debug_level >= 2)
					Fprintf(stderr, 
					 "\t%d %d\n", xdata[i], ydata[i]);
			}

			GET3(bundle, unused)
				if (debug_level >= 2)
				Fprintf(stderr, 
				 "\tbundle %d, unused %d.\n", bundle, unused);
			break;
		case TEXT:
			if (debug_level >= 1)
				Fprintf(stderr, "text\t(%d words)\n", length);
			spsassert(length >= 6, "incomplete text command");
			GET2(x)
				GET2(y)
				if (debug_level >= 2)
				Fprintf(stderr, "\tx %d, y %d\n", x, y);
			GET3(bundle, unused)
				if (debug_level >= 2)
				Fprintf(stderr, 
				 "\tbundle %d, unused %d.\n", bundle, unused);
			GET3(size, rotation)
				if (debug_level >= 2)
				Fprintf(stderr, 
				 "\tsize %d, rotation %d.\n", size, rotation);
			for (i = 0; i < 2 * (length - 5); i += 2)
				GETCHRS(chdata[i], chdata[i + 1])
					chdata[i] = '\0';
			if (debug_level >= 2)
				Fprintf(stderr, "\t%s\n", chdata);

			cosrot = cos(PI * rotation / 128.0);
			sinrot = sin(PI * rotation / 128.0);
			text(L_ROUND(x - (2.0 * cosrot - 2.5 * sinrot) * size),
			  L_ROUND(y - (2.0 * sinrot + 2.5 * cosrot) * size),
			   (long) size * 4, ROUND(rotation * 360.0 / 256.0),
			     chdata, plotline);

			break;
		case HARDWARE:
			if (debug_level >= 1)
				Fprintf(stderr, 
				 "hardware text\t(%d words)\n", length);
			spsassert(length >= 4, 
			 "incomplete hardware text command");
			GET2(x)
				GET2(y)
				if (debug_level >= 2)
				Fprintf(stderr, "\tx %d, y %d\n", x, y);
			for (i = 0; i < 2 * (length - 3); i += 2)
				GETCHRS(chdata[i], chdata[i + 1])
					chdata[i] = '\0';
			if (debug_level >= 2)
				Fprintf(stderr, "\t%s\n", chdata);

			move(x, y);
			label(chdata);

			break;
		case COMMENT:
			if (debug_level >= 1)
				Fprintf(stderr, 
				 "comment\t(%d words)\n", length);
			for (i = 0; i < length - 1; i++) {
				GET2(unused)
					if (debug_level >= 2)
					Fprintf(stderr, "%x\n", unused);
			}
			break;
		default:
			Fprintf(stderr, 
			 "%s: unimplemented plotting command (code %d).\n",
				ProgName, code);
			exit(1);
			break;
		}
	}
	(void)fclose(tfile);
}


int
plotline(n, u, v)
	long            n;
	long           *u, *v;
{
	long            i;

	spsassert(n > 0, "cant't draw arc with no coordinates.");
	move((int) u[0], (int) v[0]);
	for (i = 1; i < n; i++)
		cont((int) u[i], (int) v[i]);
}

void
convert_options(Argv,Argc)
char **Argv;
int  Argc;
/* find -i xx and change it to -WI xx,
   and  -t xx and change it to -Wl xx   */

{
	int i;
	
	for (i=0; i<Argc; i++) {
		if (strcmp(Argv[i], "-i") == 0) Argv[i] = "-WI";
		if (strcmp(Argv[i], "-t") == 0) Argv[i] = "-Wl";
	}
	return;
}

/* on interrupt or kill, remove the temp file
*/
onintr()
{
        (void)unlink(tmpname);
        exit(1);
}