#!/bin/sh
# file: plot3d G. Moody 16 May 1995
# Last revised: 5 May 1999
#
# -----------------------------------------------------------------------------
# plot3d: Make a 3D plot of stdin or a named file using 'gnuplot' in batch mode
# Copyright (C) 1999 George B. Moody
#
# This program is free software; you can redistribute it and/or modify it under
# the terms of the GNU General Public License as published by the Free Software
# Foundation; either version 2 of the License, or (at your option) any later
# version.
#
# This program is distributed in the hope that it will be useful, but WITHOUT
# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
# FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
# details.
#
# You should have received a copy of the GNU General Public License along with
# this program; if not, see <http://www.gnu.org/licenses/>.
#
# You may contact the author by e-mail (wfdb@physionet.org) or postal mail
# (MIT Room E25-505A, Cambridge, MA 02139 USA). For updates to this software,
# please visit PhysioNet (http://www.physionet.org/).
# _____________________________________________________________________________
#
# This shell script is a simple extension of `plot2d', which is itself a
# quick-and-dirty partial replacement for `plt'. `plt', however, does not
# directly support 3D plotting, so 'plot3d' may be useful even for 'plt' users.
#
# See `plot2d.1' (type `man plot2d') for usage and examples.
# Initialize variables.
CFILE=/tmp/gpcommands.$$ # gnuplot command file
TFILE=/tmp/gpdata.$$ # gnuplot data file (used only for piped input)
IFILE='-' # input file: '-' means use standard input
PRINTER='-' # printer name for lpr: '-' means don't print (plot on-screen)
PIPED=no # if 'yes', input is arriving via a pipe
ROTX=60 # viewpoint azimuth in degrees (range 0 to 180)
ROTZ=30 # viewpoint altitude in degrees (range 0 to 360)
TITLE='' # title for the top of the plot (may include spaces if quoted)
XCOL='-2' # column number for abscissas: -2 indicates XCOL not specified
XLABEL='' # label for x-axis (may include spaces if quoted)
XMIN='' # minimum abscissa for data to be included in the plot
XMAX='' # maximum abscissa for data to be included in the plot
XTICK='' # x-axis tick interval (currently ignored)
YCOL='-2' # column number for ordinates: -2 indicates YCOL not specified
YLABEL='' # label for y-axis (may include spaces if quoted)
YMIN='' # minimum y-ordinate for data to be included in the plot
YMAX='' # maximum y-ordinate for data to be included in the plot
YTICK='' # y-axis tick interval (currently ignored)
ZCOL='-2' # column number for ordinates: -2 indicates ZCOL not specified
ZLABEL='' # label for z-axis (may include spaces if quoted)
ZMIN='' # minimum z-ordinate for data to be included in the plot
ZMAX='' # maximum z-ordinate for data to be included in the plot
ZTICK='' # z-axis tick interval (currently ignored)
# Read the argument list and interpret the arguments.
while [ "x$1" != x ]
do
case $1 in
-rx) ROTX=$2; shift ;;
-rz) ROTZ=$2; shift ;;
-t) TITLE=$2; shift ;;
-x) XLABEL=$2; shift ;;
-y) YLABEL=$2; shift ;;
-z) ZLABEL=$2; shift ;;
-xa) XMIN=$2; shift; XMAX=$2; shift; XTICK=$2; shift ;;
-ya) YMIN=$2; shift; YMAX=$2; shift; YTICK=$2; shift ;;
-za) ZMIN=$2; shift; ZMAX=$2; shift; ZTICK=$2; shift ;;
-T) PRINTER=$2;shift ;;
-X) XMIN=$2; shift; XMAX=$2; shift ;;
-Y) YMIN=$2; shift; YMAX=$2; shift ;;
-Z) ZMIN=$2; shift; ZMAX=$2; shift ;;
0|1|2|3|4|5|6|7|8|9)
case $ZCOL in
-2) ZCOL=$1 ;;
*) case $XCOL in
-2) XCOL=$ZCOL; ZCOL=$1 ;;
*) YCOL=$ZCOL; ZCOL=$1 ;;
esac ;;
esac ;;
-*) echo "usage: plot3d [IFILE] [ [ XCOL YCOL ] ZCOL ] [ OPTIONS ...]"
echo " IFILE is the input file. It should contain one or more space- or"
echo " tab-separated columns per line, with each point on a line. Omit"
echo " IFILE to read data from the standard input."
echo " XCOL, YCOL, and ZCOL specify the column numbers of the columns"
echo " that contain the x, y, and z coordinates. The leftmost column"
echo " is column 0. Omit XCOL and YCOL to generate x and y coordinates"
echo " sequentially based on row numbers; a blank line in the input"
echo " resets x and increments y in this case."
echo " OPTIONS include:"
echo " -h Print this help and exit (no plot is made)"
echo " -t TITLE Use TITLE as the plot title"
echo " -x LABEL Use LABEL as the X-axis label"
echo " -y LABEL Use LABEL as the Y-axis label"
echo " -z LABEL Use LABEL as the Z-axis label"
echo " -X XMIN XMAX Plot x-coordinates between XMIN and XMAX only"
echo " -Y YMIN YMAX Plot y-coordinates between YMIN and YMAX only"
echo " -Z ZMIN ZMAX Plot z-coordinates between ZMIN and ZMAX only"
echo " -T PRINTER Plot on the specified PostScript PRINTER (default:"
echo " plot on-screen; use '-T eps' to generate"
echo " encapsulated PostScript on the standard output)"
exit ;;
*) IFILE=$1 ;;
esac
shift
done
# Renumber the columns for gnuplot. Columns are numbered beginning with 1 by
# gnuplot, and beginning with 0 by plot3d. There are four possibilities,
# depending on how many of XCOL, YCOL, and ZCOL were specified:
# Specified Columns plotted
# XCOL YCOL ZCOL X Y Z
# no no no 1 2 3
# no no yes 0 (*) 0 (*) ZCOL+1
# yes no yes XCOL+1 XCOL+2 ZCOL+1
# yes yes yes XCOL+1 YCOL+1 ZCOL+1
#
# (*) In this case, gnuplot generates x- and y-coordinates from the row numbers
# (see the help text above).
case $ZCOL in
-2) XCOL=1; YCOL=2; ZCOL=3 ;;
*) ZCOL=`expr $ZCOL + 1`
case $XCOL in
-2) XCOL=0; YCOL=0 ;;
*) XCOL=`expr $XCOL + 1`
case $YCOL in
-2) YCOL=`expr $XCOL + 1` ;;
*) YCOL=`expr $YCOL + 1` ;;
esac ;;
esac ;;
esac
case $YCOL in
0) PCOL=$ZCOL; PTYPE=noparametric ;;
*) PCOL="$XCOL:$YCOL:$ZCOL"; PTYPE=parametric ;;
esac
# If receiving input from a pipe, save it in a file for plotting. (This
# program treats redirected standard input in the same way.)
case $IFILE in
-) IFILE=$TFILE; PIPED=yes; cat >$TFILE ;;
esac
# Generate the gnuplot command file.
#
# The first section is omitted unless a -T option has been specified. The
# "-T eps" option produces encapsulated PostScript on the standard output;
# use it to generate a figure to be included within another document, or to
# be rendered at a later time. gnuplot's PostScript options follow the
# "set terminal postscript ..." lines below. The commonly-used alternatives
# to those used here are "portrait" (instead of "landscape"), "color" (instead
# of "monochrome"), "Courier" or "Helvetica" (instead of "Times-Roman"), and
# "10", "12", or "16" (instead of "14", the type size in points). See the
# gnuplot manual for details if you wish to customize these options.
case $PRINTER in
-) cp /dev/null $CFILE ;;
eps) cat >$CFILE <<!
set terminal postscript eps
!
;;
*) cat >$CFILE <<!
set terminal postscript landscape monochrome "Times-Roman" 14
set output "| lpr -P$PRINTER"
!
;;
esac
# The next section of the gnuplot command file sets the plot parameters; the
# last line in this section ("splot ...") actually produces the plot itself.
cat >>$CFILE <<!
set tics out
set title '$TITLE'
set xlabel '$XLABEL'
set ylabel '$YLABEL'
set zlabel '$ZLABEL'
set xrange [$XMIN:$XMAX]
set yrange [$YMIN:$YMAX]
set zrange [$ZMIN:$ZMAX]
set $PTYPE
set view $ROTX,$ROTZ
splot '$IFILE' using $PCOL notitle with lines
!
# The last section of the gnuplot command file is omitted if PostScript output
# is to be produced (either on a printer or to a file). Otherwise, the "pause"
# command causes gnuplot to prompt the user after it renders the plot.
case $PRINTER in
-) cat >>$CFILE <<!
pause -1 "Press <RETURN> to dismiss the plot window: "
!
;;
esac
# Make the plot by executing the gnuplot command file. Standard input is
# reattached to the user's keyboard so that the response to the prompt can
# be read.
gnuplot $CFILE </dev/tty
# Clean up and exit.
rm -f $CFILE $TFILE
exit 0