---

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

---