---

Sadie_Plot.c

Go to the documentation of this file.

/*
##########################################
# Sadie_Plot.c -
#   Set of routines for creating plots of SADIE data in the
#   session log through tcl/tk.
#
# RCS: $Id: Sadie_Plot.c,v 2.2 1999/03/20 00:34:26 conner Exp $ 
# 
# Digital Image Analysis Lab
# Dept. of Electrical and Computer Engineering
# University of Arizona
##########################################
*/

#include <tcl.h>
#include <tk.h>
#include <stdio.h>
#include <string.h>
#include "sadie.h"

/* RCS Indentification information */
static char rcsid[] =
  "$Id: Sadie_Plot.c,v 2.2 1999/03/20 00:34:26 conner Exp $";

#define FALSE 0
#define TRUE  1

extern Tcl_Interp *MainTclInterp;

int **x_array_buffer;           /* The data where info is stored      */
int point;                      /* Counter for points used            */
int count_point;                /* Flag decide to count or fill array */
int plotcount = 0;

int
Sadie_Plot_ProfileCmd (ClientData client_data, Tcl_Interp * interp, int argc,
                       char *argv[])
{

  int i, j;
  int startx, endx;
  int starty, endy;
  int option;
  int inimgaddr;
  IMAGE *inimg;
  char msg[SLEN];
  PIXEL **arrays = NULL;
  char minlbl[20], maxlbl[20];

  if (argc <= 6)
    {
      return TCL_ERROR;
    }

  sscanf (argv[1], "%x", &inimgaddr);
  inimg = (IMAGE *) inimgaddr;

  if (Tcl_GetInt (interp, argv[2], &startx) != TCL_OK)
    {
      return TCL_ERROR;
    }

  if (Tcl_GetInt (interp, argv[3], &starty) != TCL_OK)
    {
      return TCL_ERROR;
    }

  if (Tcl_GetInt (interp, argv[4], &endx) != TCL_OK)
    {
      return TCL_ERROR;
    }

  if (Tcl_GetInt (interp, argv[5], &endy) != TCL_OK)
    {
      return TCL_ERROR;
    }

  if (Tcl_GetInt (interp, argv[6], &option) != TCL_OK)
    {
      return TCL_ERROR;
    }

  /* check input */
  if (option != UNIT && option != NORM && option != LOG)
    {
      MESSAGE ('E', " Option must be one of UNIT, NORM, or LOG.");
      return TCL_ERROR;
    }

  /* Count points */
  count_point = TRUE;
  point = 0;
  brshnm (startx, starty, endx, endy);

  /* Allocate memory */
  x_array_buffer = (int **) malloc (point * sizeof (int *));
  if (x_array_buffer == (int **) NULL)
    {
      sprintf (msg, "No memory for x_array_buffer");
      MESSAGE ('E', msg);
      goto the_end;
    }

  for (i = 0; i < point; i++)
    {
      x_array_buffer[i] = (int *) calloc (2, sizeof (int));
      if (x_array_buffer[i] == (int) NULL)
        {
          sprintf (msg, "No memory for element %d", x_array_buffer[i]);
          MESSAGE ('E', msg);
          goto the_end;
        }
    }

  /* Fill array */
  count_point = FALSE;
  point = 0;
  brshnm (startx, starty, endx, endy);

  /* Allocate memory for plotting */
  arrays = (PIXEL **) malloc (inimg->nbnd * sizeof (PIXEL *));
  if (!arrays)
    {
      MESSAGE ('E', " Error allocating memory for plotting!");
      goto the_end;
    }
  for (i = 0; i < inimg->nbnd; i++)
    {
      arrays[i] = NULL;
      arrays[i] = (PIXEL *) malloc (point * sizeof (PIXEL));
      if (!arrays[i])
        {
          MESSAGE ('E', " Error allocating memory for plotting!");
          goto the_end;
        }
    }

  /* Fill the plot array */
  for (i = 0; i < inimg->nbnd; i++)
    {
      for (j = 0; j < point; j++)
        {
          arrays[i][j] =
            inimg->data[i][x_array_buffer[j][1]][x_array_buffer[j][0]];
        }
    }

  sprintf (minlbl, "(%d,%d)", startx, starty);
  sprintf (maxlbl, "(%d,%d)", endx, endy);
  PLOT (arrays, inimg->nbnd, point, minlbl, maxlbl, option);

  if (NAMES)
    {
      MESSAGE ('I', " ");
      MESSAGE ('I', " ...............");
    }

the_end:

  /* Free Array */
  for (i = 0; i < point; i++)
    {
      free (x_array_buffer[i]);
      x_array_buffer[i] = NULL;
    }

  free (x_array_buffer);
  x_array_buffer = NULL;

  if (arrays)
    {
      for (i = 0; i < inimg->nbnd; i++)
        {
          if (arrays[i])
            free (arrays[i]);
          arrays[i] = NULL;
        }
      free (arrays);
      arrays = NULL;
    }

  return TCL_OK;

}

int
Sadie_Plot_StatisticsCmd (ClientData client_data, Tcl_Interp * interp,
                          int argc, char *argv[])
{
  int imgaddr;
  IMAGE *in;
  int incr, nbins, matrix;
  double *mean = 0, *dev = 0, *cov = 0, *cor = 0;
  PIXEL *gmin = 0, *gmax = 0, *grng = 0;
  long *zero = 0;
  int rangeoption;
  PIXEL rangemin, rangemax;
  double tempval;

  if (argc <= 7)
    {
      return TCL_ERROR;
    }

  sscanf (argv[1], "%x", &imgaddr);
  in = (IMAGE *) imgaddr;

  if (Tcl_GetInt (interp, argv[2], &incr) != TCL_OK)
    {
      return TCL_ERROR;
    }

  if (Tcl_GetInt (interp, argv[3], &nbins) != TCL_OK)
    {
      return TCL_ERROR;
    }

  if (Tcl_GetInt (interp, argv[4], &matrix) != TCL_OK)
    {
      return TCL_ERROR;
    }

  if (Tcl_GetInt (interp, argv[5], &rangeoption) != TCL_OK)
    {
      return TCL_ERROR;
    }

  if (Tcl_GetDouble (interp, argv[6], &tempval) != TCL_OK)
    {
      return TCL_ERROR;
    }
  rangemin = (PIXEL) tempval;

  if (Tcl_GetDouble (interp, argv[7], &tempval) != TCL_OK)
    {
      return TCL_ERROR;
    }
  rangemax = (PIXEL) tempval;

  if (!CHECKIMG (in))
    {
      Tcl_AppendResult (interp,
                        "couldn't identify SADIE image from input string.",
                        NULL);
      return TCL_ERROR;
    }
  if (rangeoption != IGNORE && rangeoption != INCLUDE
      && rangeoption != EXCLUDE)
    {
      MESSAGE ('E', " Option must be one of IGNORE, INCLUDE, or EXCLUDE.");
      return TCL_ERROR;
    }

  RANGE (in);

  /* allocate space */
  if (!(mean = (double *) malloc (in->nbnd * sizeof (double))))
    {
      MESSAGE ('E', " Memory request failed.");
      goto the_end;
    }
  if (!(dev = (double *) malloc (in->nbnd * sizeof (double))))
    {
      MESSAGE ('E', " Memory request failed.");
      goto the_end;
    }
  if (!(gmin = (PIXEL *) malloc (in->nbnd * sizeof (PIXEL))))
    {
      MESSAGE ('E', " Memory request failed.");
      goto the_end;
    }
  if (!(gmax = (PIXEL *) malloc (in->nbnd * sizeof (PIXEL))))
    {
      MESSAGE ('E', " Memory request failed.");
      goto the_end;
    }
  if (!(grng = (PIXEL *) malloc (in->nbnd * sizeof (PIXEL))))
    {
      MESSAGE ('E', " Memory request failed.");
      goto the_end;
    }
  if (!(zero = (long *) malloc (in->nbnd * sizeof (long))))
    {
      MESSAGE ('E', " Memory request failed.");
      goto the_end;
    }
  if (!(cov = (double *) malloc (in->nbnd * in->nbnd * sizeof (double))))
    {
      MESSAGE ('E', " Memory request failed.");
      goto the_end;
    }
  if (!(cor = (double *) malloc (in->nbnd * in->nbnd * sizeof (double))))
    {
      MESSAGE ('E', " Memory request failed.");
      goto the_end;
    }

  STATS (in, incr, nbins, rangeoption, rangemin, rangemax, mean, dev, gmin,
         gmax, grng, zero);

  if (matrix)
    COVAR (in, incr, cov, cor);

the_end:
  if (cov)
    free (cov);
  if (cor)
    free (cor);
  if (mean)
    free (mean);
  if (dev)
    free (dev);
  if (gmin)
    free (gmin);
  if (gmax)
    free (gmax);
  if (grng)
    free (grng);
  if (zero)
    free (zero);

  return TCL_OK;
}

int
Sadie_Plot_CreateROICmd (ClientData client_data, Tcl_Interp * interp,
                         int argc, char *argv[])
{
  int imgaddr;
  IMAGE *in;
  char *name;
  int incr, nlev, matrix;
  int point1x, point2x, point1y, point2y;
  int xbgn, ybgn, xend, yend;

  if (argc <= 9)
    {
      return TCL_ERROR;
    }

  sscanf (argv[1], "%x", &imgaddr);
  in = (IMAGE *) imgaddr;

  if (Tcl_GetInt (interp, argv[2], &incr) != TCL_OK)
    {
      return TCL_ERROR;
    }

  if (Tcl_GetInt (interp, argv[3], &nlev) != TCL_OK)
    {
      return TCL_ERROR;
    }

  if (Tcl_GetInt (interp, argv[4], &matrix) != TCL_OK)
    {
      return TCL_ERROR;
    }

  if (Tcl_GetInt (interp, argv[5], &point1x) != TCL_OK)
    {
      return TCL_ERROR;
    }
  if (Tcl_GetInt (interp, argv[6], &point1y) != TCL_OK)
    {
      return TCL_ERROR;
    }
  if (Tcl_GetInt (interp, argv[7], &point2x) != TCL_OK)
    {
      return TCL_ERROR;
    }
  if (Tcl_GetInt (interp, argv[8], &point2y) != TCL_OK)
    {
      return TCL_ERROR;
    }

  name = argv[9];

  if (!CHECKIMG (in))
    {
      Tcl_AppendResult (interp,
                        "couldn't identify SADIE image from input string.",
                        NULL);
      return TCL_ERROR;
    }

  RANGE (in);

  xbgn = min (point1x, point2x);
  xend = max (point1x, point2x);
  ybgn = min (point1y, point2y);
  yend = max (point1y, point2y);

  CREATE_ROI (in, name, xbgn, ybgn, xend, yend, matrix, incr, nlev);

the_end:

  return TCL_OK;
}

int
Sadie_Plot_PlotArrayCmd (PIXEL * array, int arraysize)
{
  int i, j;
  PIXEL gmin, gmax;
  int minval, maxval;
  char minlbl[20], maxlbl[20];

  /* Find min and max plot values */
  gmax = gmin = array[0];
  for (i = 0; i < arraysize; i++)
    {
      gmin = min (gmin, array[i]);
      gmax = max (gmax, array[i]);
    }

  minval = (int) rnd (gmin);
  maxval = (int) rnd (gmax);

  sprintf (minlbl, "%d", minval);
  sprintf (maxlbl, "%d", maxval);

  PLOT (&array, 1, arraysize, (char *) minlbl, (char *) maxlbl, NORM);
}

/*************************************************************************
 *
 * This procedure initializes the the Sadie_Plot procedures when
 * this package is loaded.
 *
 *************************************************************************/
int
Sadie_Plot_Init (Tcl_Interp * interp)
{
  Tcl_CreateCommand (interp, "Sadie_Plot_Profile", Sadie_Plot_ProfileCmd,
                     (ClientData) NULL, NULL);
  Tcl_CreateCommand (interp, "Sadie_Plot_Statistics",
                     Sadie_Plot_StatisticsCmd, (ClientData) NULL, NULL);
  Tcl_CreateCommand (interp, "Sadie_Plot_CreateROI", Sadie_Plot_CreateROICmd,
                     (ClientData) NULL, NULL);

  return TCL_OK;
}

---