comparison tools/bufframe.c @ 0:5242703e91d3 tip

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

+/*
+    bufframe.c    shifting AIM output frame buffer.
+    ----------
+
+
+
+Examples:
+
+1. To convert gennap output to multiple animated frames:
+
+gennap len=16ms display=off output=stdout file | \
+bufframe width=8ms frstep=0.2ms anim=on > file.sai
+gensai useprev=on top=1000 cegc             -(for landscape plot)
+genspl useprev=on top=1000 pensize=2 cegc   -(for spiral plot)
+
+    (Note: spirals look better in a square box, so you might use options:
+	   dencf=1 width=500 height=500 )
+
+2. To view the basilar membrane from a cross section, animating the waves on it.
+
+genbmm mincf=220 maxcf=660 len=8ms output=stdout display=off file | \
+bufframe width=1p frstep=1p Tran=on display=on anim=on > foo.sai
+gensai bott=-100 useprev=on mag=.2 foo
+
+or:
+
+genbmm mincf=220 maxcf=660 len=32ms output=stdout display=off file | \
+bufframe width=1p frstep=1p Tran=on Header=off > foo
+
+x11play -n75 -a500 foo
+
+*/
+
+
+
+#include <stdio.h>
+#include <math.h>
+#include "header.h"
+#include "options.h"
+#include "units.h"
+#include "strmatch.h"
+
+
+char applic[]     = "shifting AIM output frame buffer. " ;
+
+
+static char *helpstr,   *debugstr, *dispstr,    *anistr,    *transtr  ;
+static char *widthstr,  *shiftstr, *headerstr,  *framestr             ;
+
+static Options option[] = {
+    {   "help"      ,   "off"       ,  &helpstr     ,   "help"                                  , DEBUG   },
+    {   "debug"     ,   "off"       ,  &debugstr    ,   "debugging switch"                      , DEBUG   },
+    {   "width"     ,   "32ms"      ,  &widthstr    ,   "Width of image frame."                 , VAL     },
+    {   "frstep"    ,   "16ms"      ,  &shiftstr    ,   "Step between image frames."            , VAL     },
+    {   "frames"    ,   "min-max"   ,  &framestr    ,   "Select frames inclusively"             , VAL     },
+    {   "Transpose" ,   "off"       ,  &transtr     ,   "Transpose image frames"                , SETFLAG },
+    {   "Header"    ,   "on"        ,  &headerstr   ,   "Header (for gensai useprevious=on)."   , VAL     },
+    {   "display"   ,   "on"        ,  &dispstr     ,   "Display output image."                 , SETFLAG },
+    {   "animate"   ,   "off"       ,  &anistr      ,   "Animate cartoon."                      , SETFLAG },
+   ( char * ) 0 } ;
+
+
+int        frameheight,    framewidth,    frames,    framebytes ;
+int     Newframeheight, Newframewidth, Newframes, Newframebytes ;
+
+int        rows,    cols,    numframes ;
+int     Newrows, Newcols, Newnumframes ;
+
+int     samplerate  ;
+int     frstep      ;
+int     application ;
+int     startbytes  ;
+int     frshift     ;
+int     Transpose   ;
+
+main(argc, argv)
+int   argc ;
+char *argv[] ;
+{
+    FILE   *fp ;
+    char   *header, *Header() ;
+    short  *frame, *endframe ;
+    char   *val1, *val2 ;
+    int     i,j,k, a,b ;
+
+    fp = openopts( option,argc,argv ) ;
+    if ( !isoff( helpstr ) )
+	helpopts( helpstr, argv[0], applic, option ) ;
+
+    if ( ( header = ReadHeader( fp ) ) == (char *) 0 ) {
+	fprintf( stderr, "edframe: header not found\n" ) ;
+	exit( 1 ) ;
+    }
+    frameheight = HeaderInt( header, "frameheight" );
+    framewidth  = HeaderInt( header, "framewidth"  );
+    frames      = HeaderInt( header, "frames"      );
+    samplerate  = HeaderSamplerate( header );
+
+    if ( ( application = Applic( header) ) < 0 ) {
+	fprintf(stderr,"bufframe: application name not found in header\n" ) ;
+	exit( 1 ) ;
+    }
+    if ( Format( application ) != NAP ) {
+	fprintf(stderr,"bufframe: buffer for %s format not implemented\n", gen_applics[application] ) ;
+	exit( 1 ) ;
+    }
+    frame_to_matrix( NAP, &rows, &cols, &numframes, frameheight, framewidth, frames ) ;
+
+    Transpose    = ison( transtr ) ;
+    frstep       = to_p( shiftstr, samplerate ) ;
+    Newrows      = frameheight ;
+    Newcols      = to_p( widthstr, samplerate ) ;
+    Newnumframes = ( cols - Newcols ) / frstep + 1 ;    /* since: Newcols + ( Newnumframes - 1 ) * frstep <= cols */
+
+
+    /* Get limits on specified number of frames */
+
+    if ( getvals( framestr, &val1, &val2, "-" ) == BADVAL ) {
+	fprintf(stderr,"bufframe: bad frame selector [%s]\n", framestr ) ;
+	exit( 1 ) ;
+    }
+    if      ( ismin( val1 ) ) a = 0 ;
+    else if ( ismax( val1 ) ) a = Newnumframes - 1 ;
+    else if ( Units( val1 ) ) a = (int)to_p( val1, samplerate ) / frstep ;
+    else                      a = atoi( val1 ) ;
+
+    if ( isempty( val2 ) )    b = a ;
+    else if ( ismin( val2 ) ) b = 0 ;
+    else if ( ismax( val2 ) ) b = Newnumframes - 1 ;
+    else if ( Units( val2 ) ) b = (int)to_p( val2, samplerate ) / frstep ;
+    else                      b = atoi( val2 ) ;
+
+    if (a<0 || b>=Newnumframes || a>b) {
+	if ( b>=Newnumframes ) fprintf(stderr,"bufframe: bad frame selector [%d frames needed but %d available]\n", b+1, Newnumframes );
+	else                fprintf(stderr,"bufframe: bad frame selector \n" );
+	exit(1);
+    }
+
+    Newnumframes = b - a + 1 ;
+
+    framebytes   = Newrows * ( Newcols + ( Newnumframes - 1 ) * frstep ) * sizeof(short) ;
+
+    if ( Transpose )
+	matrix_to_frame( SAI, Newcols, Newrows, Newnumframes, &Newframeheight, &Newframewidth, &Newframes ) ;
+    else
+	matrix_to_frame( SAI, Newrows, Newcols, Newnumframes, &Newframeheight, &Newframewidth, &Newframes ) ;
+
+    Newframebytes  = Newframeheight * Newframewidth * sizeof(short) ;
+
+    if ( ison( headerstr ) )
+	WriteHeader( Header( header ), stdout ) ;
+
+
+    /* Allocate space for framebytes of data */
+
+    if ( ( frame = (short *)malloc( framebytes ) ) == NULL ) {
+	fprintf( stderr, "bufframe: malloc out of space\n" ) ;
+	exit( 1 ) ;
+    }
+
+    frshift = frstep * rows ;
+
+    /* Seek past `a' frames in blocks of framebytes */
+
+    startbytes = a * frshift * sizeof(short) ;
+    for (i=framebytes ; i < startbytes ; i += framebytes)
+	if ( fread( frame, framebytes, 1, fp ) == NULL ) {
+	    fprintf(stderr,"bufframe: missing data after header\n");
+	    exit(1);
+	}
+    if ( (startbytes -= (i - framebytes)) > 0 )
+	if ( fread( frame, startbytes, 1, fp ) == NULL ) {
+	    fprintf(stderr,"bufframe: missing data after header\n");
+	    exit(1);
+	}
+
+    /* Read framebytes of data */
+
+    if ( fread( frame, framebytes, 1, fp ) == NULL ) {
+	fprintf(stderr,"bufframe: insufficient data after header\n");
+	exit(1);
+    }
+    fclose( fp ) ;
+
+    endframe = frame + Newframes * frshift ;
+
+    fprintf(stderr,"bufframe: %d frames   [height=%d  width=%d]\n", Newnumframes, Newframeheight, Newframewidth );
+
+    if ( ison( debugstr ) ) {
+	fprintf( stderr,"   frameheight=%3d    framewidth=%3d       frames=%3d    framebytes=%d\n",     frameheight,    framewidth,    frames,    framebytes ) ;
+	fprintf( stderr,"Newframeheight=%3d Newframewidth=%3d    Newframes=%3d Newframebytes=%d\n",  Newframeheight, Newframewidth, Newframes, Newframebytes ) ;
+	fprintf( stderr,"          rows=%3d          cols=%3d    numframes=%d\n",     rows,    cols,    numframes ) ;
+	fprintf( stderr,"       Newrows=%3d       Newcols=%3d Newnumframes=%d\n",  Newrows, Newcols, Newnumframes ) ;
+	fprintf( stderr,"frstep=%d  startbytes=%d frshift=%d\n",  frstep, startbytes, frshift ) ;
+	exit( 1 ) ;
+    }
+
+    if ( Transpose ) {
+
+	for (  ; frame < endframe ; frame += frshift )
+	    for ( j=0 ; j < Newcols ; j++ )
+		fwrite( frame + j*frameheight, sizeof(short), Newrows, stdout ) ;
+    }
+
+    else {
+	for (  ; frame < endframe ; frame += frshift )
+	    for (j=0 ; j < Newrows ; j++)
+		for (k=0 , i=j ; k < Newcols ; k++, i+=frameheight)
+		    fwrite( &frame[i], sizeof(short), 1, stdout );
+    }
+    fprintf(stderr,"bufframe: done\n" );
+}
+
+
+
+/* Return 1 if the str has appended units. Otherwise return 0.             */
+/* (A freqency specifier with no units is interpreted as a channel number) */
+
+Units(str)
+char *str ;
+{
+    char *eptr = str + strlen( str ) ;
+
+    if( isdigit( *--eptr ) ) return 0 ; /* last char is digit, so no units */
+    else return 1 ;
+}
+
+
+/*
+   Copy the original (NAP format) header to a new sai header, changing in order:
+     frstep
+     frames
+     framewidth
+     frameheight
+     framebytes
+     animate
+     display
+   Copy application name "sai" into the version string.
+   Finally, update the new header_bytes, and return the new header.
+*/
+
+char *Header( oldheader )
+char *oldheader ;
+{
+    char *saiheader;
+    char *p0, *p1, *p2, *s, str[64];
+
+    saiheader = (char *)malloc( strlen(oldheader) + 256 ) ;
+
+    p0 = saiheader ;
+    p1 = oldheader ;
+
+    /* copy to first item after the header_bytes */
+
+    p2 = HeaderString( oldheader , "header_bytes" ) ;
+    while( p1 < p2 )
+	*p0++ = *p1++ ;
+    while (*p1 != '\n')
+	*p0++ = *p1++ ;
+    *p0++ = *p1++;
+
+    /* insert frstep_aid at top of new header */
+
+    sprintf( str,"frstep_aid=%s\n", shiftstr ) ;
+    for (s = str ; *s != '\n' ; )
+	*p0++ = *s++;
+    *p0++ = *s;
+
+
+    /** copy up to channels_afb (if transpose changes apparent number of channels) **/
+
+    if ( Transpose ) {
+	p2 = HeaderString( oldheader , "channels_afb" ) ;
+	while( p1 < p2 )
+	    *p0++ = *p1++ ;
+
+	sprintf(str,"%d\n", Newframeheight);
+	for (s = str ; *s != '\n' ; )
+	    *p0++ = *s++;
+	*p0++ = *s;
+	while (*p1 != '\n')
+	    *p1++;
+	*p1++;
+    }
+
+
+    /** copy up to frames **/
+
+    p2 = HeaderString( oldheader , "frames" ) ;
+    while( p1 < p2 )
+	*p0++ = *p1++ ;
+
+    sprintf(str,"%d\n", Newframes);
+    for (s = str ; *s != '\n' ; )
+	*p0++ = *s++;
+    *p0++ = *s;
+    while (*p1 != '\n')
+	*p1++;
+    *p1++;
+
+
+    /** copy up to framewidth **/
+
+    p2 = HeaderString( oldheader , "framewidth" ) ;
+    while ( p1 < p2 )
+	*p0++ = *p1++ ;
+
+    sprintf(str,"%d\n", Newframewidth);
+    for (s = str ; *s != '\n' ; )
+	*p0++ = *s++;
+    *p0++ = *s;
+    while (*p1 != '\n')
+	*p1++;
+    *p1++;
+
+
+    /** copy up to frameheight **/
+
+    p2 = HeaderString( oldheader , "frameheight" ) ;
+    while ( p1 < p2 )
+	*p0++ = *p1++ ;
+
+    sprintf(str,"%d\n", Newframeheight);
+    for (s = str ; *s != '\n' ; )
+	*p0++ = *s++;
+    *p0++ = *s;
+    while (*p1 != '\n')
+	*p1++;
+    *p1++;
+
+
+    /** copy up to framebytes **/
+
+    p2 = HeaderString( oldheader , "framebytes" ) ;
+    while ( p1 < p2 )
+	*p0++ = *p1++ ;
+
+    sprintf(str,"%d\n", Newframebytes);
+    for (s = str ; *s != '\n' ; )
+	*p0++ = *s++;
+    *p0++ = *s;
+    while (*p1 != '\n')
+	*p1++;
+    *p1++;
+
+
+    /** copy up to animate_ctn **/
+
+    p2 = HeaderString( oldheader , "animate_ctn" ) ;
+    while ( p1 < p2 )
+	*p0++ = *p1++ ;
+
+    sprintf(str,"%s\n", anistr );
+    for (s = str ; *s != '\n' ; )
+	*p0++ = *s++;
+    *p0++ = *s;
+    while (*p1 != '\n')
+	*p1++;
+    *p1++;
+
+
+    /** copy up to display **/
+
+    p2 = HeaderString( oldheader , "display" ) ;
+    while ( p1 < p2 )
+	*p0++ = *p1++ ;
+
+    sprintf(str,"%s\n", dispstr );
+    for (s = str ; *s != '\n' ; )
+	*p0++ = *s++;
+    *p0++ = *s;
+    while (*p1 != '\n')
+	*p1++;
+    *p1++;
+
+
+    /** copy rest of header **/
+
+    p2 = HeaderString( oldheader , "Version" ) ;
+    while ( p1 < p2 )
+	*p0++ = *p1++ ;
+
+       /** ( copy application name "sai" into header ) **/
+
+    if ( ( p2 = ApplicString( oldheader ) ) == (char *)0 )  {
+	fprintf(stderr,"bufframe: application name not found in header\n");
+	exit( 1 ) ;
+    }
+    while ( p1 < p2 )
+	*p0++ = *p1++ ;
+
+    *p0++ = 's' ; *p0++ = 'a' ; *p0++ = 'i' ;
+    p1 += 3 ;
+
+    while (*p1 != '\n')
+	*p0++ = *p1++ ;
+    *p0++ = *p1++ ;
+
+
+    /** update header_bytes **/
+
+    sprintf(str, "%0*d", 7, p0-saiheader);
+    p0 = HeaderString( saiheader , "header_bytes" ) ;
+    s = str;
+    while(*p0 != '\n')
+	*p0++ = *s++ ;
+
+
+    return saiheader;
+}
+
+
+
+