---

tclSadie_ROI.c

Go to the documentation of this file.

#include "sadie.h"



void CREATE_ROI (
IMAGE *in,      /*  I   Pointer to the image                                  */
char  *name,    /*  I   Name of the region to create                          */
int   xbgn,     /*  I   Starting x coord of region of interest                */
int   ybgn,     /*  I   Starting y coord of region of interest                */
int   xend,     /*  I   Ending x coord of region of interest                  */
int   yend,     /*  I   Ending y coord of region of interest                  */
short matrix,   /*  I   Indicates whether to compute cov/cor matrix           */
short inc,      /*  I   Pixel subsample increment                             */
short nlev      /*  I   > 0:    Number of bins in the histogram               */
                /*      = 0:    Histogram is suppressed                       */
/*----------------------------------------------------------------------------*/
) {   char msg[SLEN];
    short i, j, k;
    long  nval,n;
    PIXEL *gval=NULL;
    ROIPtr r=NULL, s=NULL;
    


    if (NAMES) {
        MESSAGE('I',"");
        MESSAGE('I',"CREATE_ROI");
        MESSAGE('I',"");
        sprintf(msg," Region Name:                %s",name);
        MESSAGE('I',msg);
        sprintf(msg," Region Start Point:         ( %d , %d )",xbgn+1,ybgn+1);
        MESSAGE('I',msg);
        sprintf(msg," Region End Point:           ( %d , %d )",xend+1,yend+1);
        MESSAGE('I',msg);
        sprintf(msg," Region Height:              %d pixels",(yend-ybgn)+1);
        MESSAGE('I',msg);
        sprintf(msg," Region Width:               %d pixels",(xend-xbgn)+1);
        MESSAGE('I',msg);
        MESSAGE('I'," ...............");
    }
    
    /* check input */
    if (!CHECKIMG(in)) {
        MESSAGE('E'," Can't identify input image.");
        goto the_end;
    }
    
    /* Determine the number of pixels in the region */
    nval = (long)((xend-xbgn)+1) * (long)((yend-ybgn)+1);
    
    /* Find the end of the ROI list */
    for (i = 1, s = in->regions; s != NULL; s = s->next, i++);

    /* Allocate space for the array of pixel values */
    if ((gval    = (PIXEL  *)malloc(in->nbnd*nval*sizeof(PIXEL))) != NULL  &&
        (r       = (ROIPtr  )malloc(sizeof(ROI))) != NULL  &&
        (r->name = (char   *)malloc((strlen((char *)name)+1)*sizeof(char))) != NULL  &&
        (r->zero = (long   *)malloc(in->nbnd*sizeof(long))) != NULL  &&
        (r->gmin = (PIXEL  *)malloc(in->nbnd*sizeof(PIXEL))) != NULL  &&
        (r->gmax = (PIXEL  *)malloc(in->nbnd*sizeof(PIXEL))) != NULL  &&
        (r->mean = (double *)malloc(in->nbnd*sizeof(double))) != NULL  &&
        (r->dev  = (double *)malloc(in->nbnd*sizeof(double))) != NULL  &&
        (!matrix  ||  (r->cov  = (double *)malloc(in->nbnd*in->nbnd*sizeof(double))) != NULL)  &&
        (!matrix  ||  (r->cor  = (double *)malloc(in->nbnd*in->nbnd*sizeof(double))) != NULL)) {

      strcpy(r->name,name);
      r->nval = nval;

      for (n = i = 0; i< in->nbnd; i++) {
        for (j=ybgn; j<= yend; j++) {
          for (k=xbgn; k<= xend; k++,n++) {
            gval[n] = in->data[i][j][k];
          }
        }
      }
      
      if (!matrix) {
        r->cov = NULL;
        r->cor = NULL;
      }
      ROISTATS(r->name,in->nbnd,nval,gval,inc,nlev,r->gmin,r->gmax,r->mean,r->dev,r->zero,r->cov,r->cor);
      r->next   = NULL;
      if (in->regions == NULL) {
        in->regions = r;
      } else {
        for (s = in->regions; s->next != NULL; s = s->next);
        s->next = r;
      }
      
      

    } else {
      MESSAGE('E',"Memory request failed.");                            
      goto the_end;                                                             
    }                                                                          
    


    the_end:
    if (gval) free(gval);
}

---