---

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

---