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

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

+/*
+  naptosai.c            NAP to SAI format conversion.
+ ------------
+
+    Read a NAP header and frame.
+    Write a SAI header, (created from the NAP header).
+    Write a succession of SAI frames.
+
+    The options "start" and "length" specify the NAP which is input.
+    The options "width" and "frstep"  specify the SAI frames which are output.
+
+    Special option values are:
+
+      length=max      Input all the NAP from the given start to its end.
+      width=max       Output framewidth is set equal to the given NAP length,
+		      and if this is also "max", then the framewidth is the
+		      remainder of the NAP.
+      frstep=0        Output one frame of given width from the given start of
+		      the NAP.
+
+    For example, use naptosai with width=max and frstep=0 to produce one SAI
+    frame the same size as the input NAP.
+
+Examples:
+
+
+1. To convert a nap or a bmm output to a single frame of sai output:
+
+gennap length=64 output=stdout file | naptosai width=32 frstep=0 > file.sai
+gensai top=1000 useprevious=on file
+
+genbmm output=stdout file | naptosai  > file.sai
+gensai bottom=-100 useprevious=on file
+
+  This uses special case of zero frstep to generate one frame of output with
+  the given width. The default width is the given length, and the default
+  length is the remainder of the input file. If the width=length, then the
+  frstep is zero by default, and therefore the default process, (ie with no
+  arguments), converts the whole of the input to a single SAI frame.
+
+  Note, certain display parameters have different default values for different
+  applications. The SAI display parameters should be set to the appropriate
+  values, in order to plot the SAI on the same scale. For example:
+    When the source application is NAP, set gensai top=1000
+    When the source application is BMM, set gensai bottom=-100
+
+
+2. To convert a nap output to multiple frames of sai output:
+
+gennap length=40 output=stdout file | naptosai width=32 frstep=0.2 > file.sai
+gensai top=1000 useprevious=on file
+
+   The sai output can be plotted as a spiral, but it is best to set the
+   display parameters as follows:
+
+gennap length=40 output=stdout dencf=1 width=500 height=500 file | \
+				naptosai width=32 frstep=0.2 > file.sai
+genspl useprevious=on pensize=2 file
+
+
+3. To convert a nap output to a sai format bitmap, to display an animated
+   stream:
+
+gennap length=40 output=stdout file | naptosai width=32 frstep=0.2 > file.sai
+gensai top=1000 useprevious=on bitmap=on file
+review file
+
+  The same can be done for the spiral display:
+
+gennap length=40 output=stdout dencf=1 width=500 height=500 file | \
+				naptosai width=32 frstep=0.2 > file.sai
+genspl useprevious=on pensize=2 bitmap=on file
+review file
+
+  The plots show the nap moving though 8ms (40-32), in steps of 0.2ms,
+  first in rectangular display, then in spiral-mapped display.
+
+*/
+
+
+#include <stdio.h>
+#include <math.h>
+#include "header.h"
+#include "options.h"
+#include "units.h"
+#include "strmatch.h"
+
+char applic[]     = "NAP to SAI format conversion. " ;
+
+static char *helpstr, *debugstr, *startstr, *lengthstr, *widthstr, *shiftstr, *headerstr ;
+
+static Options option[] = {
+    {   "help"      ,   "off"       ,  &helpstr     ,   "help"                          , DEBUG   },
+    {   "debug"     ,   "off"       ,  &debugstr    ,   "debugging switch"              , DEBUG   },
+    {   "start"     ,   "0"         ,  &startstr    ,   "Start point in i/p NAP"        , VAL     },
+    {   "length"    ,   "max"       ,  &lengthstr   ,   "Length of i/p NAP to process." , VAL     },
+    {   "header"    ,   "on"        ,  &headerstr   ,   "sai header"                    , VAL     },
+    {   "width"     ,   "max"       ,  &widthstr    ,   "Width of o/p SAI frame."       , VAL     },
+    {   "frstep"    ,   "1ms"       ,  &shiftstr    ,   "Step between o/p frames of SAI.\n(frstep=0 gets 1 frame of given width)." , VAL },
+   ( char * ) 0 } ;
+
+
+int     frameheight, framewidth ;       /* Nap parameters read from header */
+int     frames, samplerate ;
+
+int     start, length ;
+
+int     Newframeheight, Newframewidth ; /* Sai parameters to write  */
+int     Newframes ;
+int     Newframeshift ;
+int     Newframebytes ;
+
+
+main(argc, argv)
+int   argc ;
+char *argv[] ;
+{
+    FILE   *fp ;
+    int     i, framebytes, startbytes, frameshift;
+    char   *header, *SaiHeader();
+    short  *nap, *frame, *endframe;
+
+    fp = openopts( option,argc,argv ) ;
+    if ( !isoff( helpstr ) )
+	helpopts( helpstr, argv[0], applic, option ) ;
+
+    if ( ison( debugstr ) ) {
+	printf("start=%s length=%s width=%s shift=%s\n",
+		startstr, lengthstr, widthstr, shiftstr );
+	exit(0);
+    }
+
+    if ( (header = ReadHeader(fp)) == (char *) 0 ) {
+	fprintf(stderr,"naptosai: header not found\n");
+	exit(1);
+    }
+
+    frameheight = HeaderInt( header, "frameheight" );
+    framewidth =  HeaderInt( header, "framewidth"  );
+    frames =      HeaderInt( header, "frames"      );
+    samplerate =  HeaderSamplerate( header, "samplerate"  );
+
+    start =  to_p( startstr,  samplerate );
+    if ( ismax( lengthstr ) )   /* Special case for remainder of input */
+	length = frames - start ;
+    else
+	length = to_p( lengthstr, samplerate );
+
+    if ( start + length > frames ) {
+	fprintf(stderr,"naptosai: nap too small (%d ms) for requested start and length \n", to_ms( frames, samplerate ) );
+	exit(1);
+    }
+
+    framebytes = frameheight * length * sizeof(short) ;
+    startbytes = frameheight * start  * sizeof(short) ;
+
+    Newframeheight = frameheight ;
+    if ( ismax( widthstr ) )    /* Special case for single frame of max width */
+	Newframewidth = length ;
+    else
+	Newframewidth = to_p( widthstr, samplerate ) ;
+    Newframeshift = to_p( shiftstr, samplerate ) ;
+    Newframebytes =  Newframeheight * Newframewidth * sizeof(short) ;
+    if ( Newframeshift > 0 )
+	Newframes = 1 + ( length - Newframewidth ) / Newframeshift ;
+    else {              /* Special case of zero shift */
+	Newframes = 1 ;
+	Newframeshift = 1 ;
+    }
+
+    if ( length < Newframewidth ) {
+	fprintf(stderr,"naptosai: nap too small (%d ms) for requested width\n", to_ms( length, samplerate ) );
+	exit(1);
+    }
+
+    fprintf(stderr,"Output %d sai frames, each %d bytes\n", Newframes, Newframebytes );
+
+    if ( ison( headerstr ) )
+	WriteHeader( SaiHeader(header), stdout );
+
+    /* Allocate space for framebytes of nap data */
+
+    if ( (nap = (short *)malloc( framebytes )) == NULL ) {
+	fprintf(stderr,"naptosai: malloc out of space\n");
+	exit(1);
+    }
+
+    /* Seek to start in blocks of framebytes */
+
+    for (i=framebytes ; i < startbytes ; i += framebytes)
+	if ( fread( nap, framebytes, 1, fp ) == NULL ) {
+	    fprintf(stderr,"naptosai: missing data after header\n");
+	    exit(1);
+	}
+    if ( (startbytes -= (i - framebytes)) > 0 )
+	if ( fread( nap, startbytes, 1, fp ) == NULL ) {
+	    fprintf(stderr,"naptosai: missing data after header\n");
+	    exit(1);
+	}
+
+    /* Read framebytes of i/p nap data */
+
+    if ( fread( nap, framebytes, 1, fp ) == NULL ) {
+	fprintf(stderr,"naptosai: missing data after header\n");
+	exit(1);
+    }
+
+    frameshift = Newframeshift * Newframeheight ;
+    endframe = nap + Newframes * frameshift ;
+
+    for ( frame = nap ; frame < endframe ; frame += frameshift )
+	writeframe( frame, stdout ) ;
+
+    fprintf(stderr,"naptosai done\n" ) ;
+
+}
+
+
+
+/* Write a frame in sai format */
+
+writeframe( frame, fp )
+short  *frame ;
+FILE   *fp ;
+{
+    int    i, row, col;
+
+    for (row=0 ; row < Newframeheight ; row++)
+	for (col=0 , i=row ; col < Newframewidth ; col++, i+=Newframeheight)
+	    fwrite( &frame[i], sizeof(short), 1, fp );
+}
+
+
+/*********************** Read header and build new header *****************/
+
+/*
+   Copy the original nap header to a new sai header, changing in order:
+     frames
+     frameshift
+     framewidth
+     frameheight
+     framebytes
+   Then change applic name [gennap] to [gensai] in the Version string.
+   Finally, update the new header_bytes, and return the new header.
+*/
+
+char *SaiHeader( napheader )
+char *napheader ;
+{
+    char *saiheader;
+    char *p0, *p1, *p2, *s, str[64];
+
+    saiheader = (char *)malloc( strlen(napheader) + 64 ) ;
+
+    p0 = saiheader ;
+    p1 = napheader ;
+
+
+    /** copy up to frames **/
+
+    p2 = HeaderString( napheader , "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 frameshift **/
+
+    p2 = HeaderString( napheader , "frameshift" ) ;
+    while ( p1 < p2 )
+	*p0++ = *p1++ ;
+
+    sprintf(str,"%d\n", Newframeshift);
+    for (s = str ; *s != '\n' ; )
+	*p0++ = *s++;
+    *p0++ = *s;
+    while (*p1 != '\n')
+	*p1++;
+    *p1++;
+
+
+    /** copy up to framewidth **/
+
+    p2 = HeaderString( napheader , "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( napheader , "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( napheader , "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 rest of header **/
+
+    p2 = HeaderString( napheader , "Version" ) ;
+    while ( p1 < p2 )
+	*p0++ = *p1++ ;
+
+    while (*p1 != ']')          /* change applic name [gennap] to [gensai] */
+	*p0++ = *p1++ ;
+    *(p0-3) = 's' ;
+    *(p0-2) = 'a' ;
+    *(p0-1) = 'i' ;
+    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;
+}
+
+