---

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;
}

---