annotate tools/gate.c @ 0:5242703e91d3 tip

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

/*
  gate.c      gate numbers from input stream
  ------
    Read and write numbers of a given type.
    All numbers within the given gate ranges are replaced by the
    given gate <value>, (a real ie. float number).
    Gate ranges are: xrange (a time range, with optional time units)
    and:             yrange (an amplitude range, in real ie. float numbers).
    If the gate is not a <value>, but is an <operation>, then the numbers
    within the gate ranges are operated on as appropriate.
    Abbreviated forms of the operation names are allowed.

	<operation>     comment
	-----------     -----------------------------
	exclude         exclude numbers from output
	negate          negate numbers
	count           print a count of the numbers found in the gate range
			on the stderr.

Examples:

1. Replace all numbers <=0 by value 0  (ie. half-wave rectification).

gate yrange=min-0 val=0  file

2. Replace all numbers <=0 by their inverse (ie. full-wave rectification).

gate yrange=min-0 val=neg  file

3. Gate the onset of a signal: replace the first 20ms with zeroes.

gate xrange=0-20ms yrange=min-max val=0 file

4. Exclude all numbers <=0

gate yrange=min-0  val=exclude file

5. Exclude all numbers >0

gate yrange=1-max val=exclude file

6. Delete lines 4 to 8 inclusive from ascii input (lines numbered 0,1,2,...)

gate type=ASCII xrange=4-8 yrange=min-max val=exclude  file

7. Replace all numbers in the yrange -1 to +1 inclusive by 0

gate yrange=-1-1 val=0  file

8. Replace all instances of number 10 by -10

gate yrange=10 val=-10  file

9. Print a count of all numbers = 0

gate yrange=0 val=count file

10. Print a count of all numbers < 0

gate yrange=min--1  val=count file

*/


#include <stdio.h>
#include <math.h>
#include "options.h"
#include "units.h"
#include "strmatch.h"

char applic[] = "gate specific numbers in input stream using given value." ;

static char *helpstr, *debugstr, *sampstr, *xrstr, *yrstr, *gstr, *typestr ;

static Options option[] = {
    {   "help"      ,   "off"       ,  &helpstr     ,   "help"                      , DEBUG   },
    {   "debug"     ,   "off"       ,  &debugstr    ,   "debugging switch"          , DEBUG   },
    {   "samplerate",   "20kHz"     ,  &sampstr     ,   "samplerate "               , VAL     },
    {   "xrange"    ,   "min-max"   ,  &xrstr       ,   "time range"                , VAL     },
    {   "yrange"    ,   "min-0"     ,  &yrstr       ,   "amplitude range"           , VAL     },
    {   "value"     ,   "0"         ,  &gstr        ,   "replacement value or <operation>", VAL     },
    {   "type"      ,   "short"     ,  &typestr     ,   "datatype"                  , VAL     },
   ( char * ) 0 } ;


int    samplerate  ;
int    type        ;    /* datatype index */

main(argc, argv)
int    argc;
char **argv;
{
    FILE   *fp ;
    int     i, n=0, ax, bx, real_range(), helpoperations() ;
    float   ay, by, x, val ;

    fp = openopts( option, argc, argv ) ;
    if ( !isoff( helpstr ) )
	helpopts1( helpstr, argv[0], applic, option, helpoperations ) ;

    samplerate = to_Hz( sampstr, 0 ) ;

    if ( ( type = typeindex( typestr ) ) < 0 ) {
	fprintf( stderr, "gate: bad type [%s]\n", typestr ) ;
	exit( 1 ) ;
    }

    if ( range( xrstr, &ax, &bx, samplerate ) == 0 ) {
	fprintf(stderr,"gate: bad xrange [%s]\n", xrstr ) ;
	exit( 1 ) ;
    }

    if ( real_range( yrstr, &ay, &by ) == 0 ) {
	fprintf(stderr,"gate: bad yrange [%s]\n", yrstr ) ;
	exit( 1 ) ;
    }

    if ( iststr( gstr, "count" ) ) {            /* gate=count */

	for ( i = 0 ; readitem( &x, type, 1, fp ) ; i++ )
	    if ( i >= ax && ( i <= bx || bx == (-1) ) )
		if ( x >= ay && x <= by )
		    n++ ;
	fprintf(stderr,"%d numbers in gate\n", n);
    }


    else if ( iststr( gstr, "exclude" ) ) {     /* gate=exclude */

	for ( i = 0 ; readitem( &x, type, 1, fp ) ; i++ ) {
	    if ( i < ax || ( i > bx && bx != (-1) ) )
		writeitem( &x, type, 1, stdout ) ;
	    else if ( x < ay || x > by )
		writeitem( &x, type, 1, stdout ) ;
	}
    }

    else if ( iststr( gstr, "negate" ) ) {      /* gate=negate  */

	for ( i = 0 ; readitem( &x, type, 1, fp ) ; i++ ) {
	    if ( i < ax || ( i > bx && bx != (-1) ) )
		writeitem( &x, type, 1, stdout ) ;
	    else if ( x < ay || x > by )
		writeitem( &x, type, 1, stdout ) ;
	    else {
		x = ( -x ) ;
		writeitem( &x, type, 1, stdout ) ;
	    }
	}
    }

    else {                                      /* gate=<value> */
	val = atof( gstr ) ;
	for ( i = 0 ; readitem( &x, type, 1, fp ) ; i++ ) {
	    if ( i < ax || ( i > bx && bx != (-1) ) )
		writeitem( &x, type, 1, stdout ) ;
	    else if ( x < ay || x > by )
		writeitem( &x, type, 1, stdout ) ;
	    else
		writeitem( &val, type, 1, stdout ) ;
	}
    }

    fclose(fp);
}


/*
Parse a range specifier into real numbers representing amplitude.
*/

int real_range( s, a, b )
char  *s ;
float *a, *b ;
{
    char *val1, *val2 ;

    if ( getvals( s, &val1, &val2, "-" ) == (-3) )
	return ( 0 ) ;
    if      ( ismin( val1 ) ) *a = (-1e10) ;    /* minimum         */
    else if ( ismax( val1 ) ) *a =   1e10  ;    /* maximum         */
    else *a = atof( val1 ) ;

    if ( isempty( val2 ) )    *b = *a ;          /* single object        */
    else if ( ismin( val2 ) ) *b = (-1e10) ;
    else if ( ismax( val2 ) ) *b =   1e10  ;
    else *b = atof( val2 ) ;

    if ( *a > *b ) return 0 ;
    return 1 ;
}


helpoperations()
{
    fprintf(stderr,"\n<operation>:          \n");
    fprintf(stderr,"            exclude     exclude numbers in gate range from output \n");
    fprintf(stderr,"            negate      negate numbers in gate range              \n");
    fprintf(stderr,"            count       print count of numbers in gate range      \n");
    exit( 1 ) ;
}