---

maxlike.c

Go to the documentation of this file.

#include        "sadie.h"

/*-Copyright Information------------------------------------------------------*/
/* Copyright (c) 1994 by the University of Arizona Digital Image Analysis Lab */
/*----------------------------------------------------------------------------*/
/*-General Information--------------------------------------------------------*/
/*                                                                            */
/*   This function performs a maximum likelihood classification on an image.  */
/*                                                                            */
/*----------------------------------------------------------------------------*/
/*-Interface Information------------------------------------------------------*/
void
MAXLIKE (IMAGE * in,            /*  I  Pointer to the input image.                             */
         short nlev,            /*  I  Number of classes                                       */
         /*     (and number of graylevels in the output image).         */
         double *cov,           /*  I  Pointer to the array                                    */
         /*     of covariance matrices[nlev][nbnd][nbnd].               */
         double *prob,          /*  I  Pointer to the array[nlev] of a priori probabilities.   */
         IMAGE ** out           /*  O  Address of a pointer to the output image.               */
/*--------------------------------------------------------------------------- */
  )
{
  register short i, j, k, m, n;
  char msg[SLEN];

  if (TIMES)
    TIMING (T_MAXLIKE);
  if (NAMES)
    {
      MESSAGE ('I', "");
      MESSAGE ('I', "MAXLIKE");
      MESSAGE ('I', "");
      sprintf (msg, " Input image:          %s", in->text);
      MESSAGE ('I', msg);
      sprintf (msg, " Number of classes:    %d", nlev);
      MESSAGE ('I', msg);
      for (i = 0; i < nlev; i++)
        {
          sprintf (msg, " Class%3d:          ", i + 1);
          MESSAGE ('I', msg);
          MESSAGE ('I', "  Covariance matrix:");
          for (j = 0; j < in->nbnd; j++)
            {
              sprintf (msg, "  Band%4d:          ", j + 1);
              for (k = 0, m = strlen (msg); k < in->nbnd;
                   k++, m = strlen (msg))
                {
                  sprintf (msg + m, "%12.4e",
                           cov[i * in->nbnd * in->nbnd + j * in->nbnd + k]);
                }
              MESSAGE ('I', msg);
            }
          sprintf (msg, "  Probability:       %12.4e", prob[i]);
          MESSAGE ('I', msg);
        }
      MESSAGE ('I', " ...............");
    }

the_end:
  if (TIMES)
    TIMING (T_EXIT);
}

---