clust_util.cpp

高斯混合模型的C语言实现

/*
* All questions regarding the software should be addressed to
* 
*       Prof. Charles A. Bouman
*       Purdue University
*       School of Electrical and Computer Engineering
*       1285 Electrical Engineering Building
*       West Lafayette, IN 47907-1285
*       USA
*       [phone]
*       [phone] (fax)
*       email:  [e-mail]
*       http://www.ece.purdue.edu/~bouman
* 
* Copyright (c) 1995 The Board of Trustees of Purdue University.
*
* Permission to use, copy, modify, and distribute this software and its
* documentation for any purpose, without fee, and without written agreement is
* hereby granted, provided that the above copyright notice and the following
* two paragraphs appear in all copies of this software.
*
* IN NO EVENT SHALL PURDUE UNIVERSITY BE LIABLE TO ANY PARTY FOR DIRECT,
* INDIRECT, SPECIAL, INCIDENTAL, OR CONSEQUENTIAL DAMAGES ARISING OUT OF THE
* USE OF THIS SOFTWARE AND ITS DOCUMENTATION, EVEN IF PURDUE UNIVERSITY HAS
* BEEN ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
*
* PURDUE UNIVERSITY SPECIFICALLY DISCLAIMS ANY WARRANTIES, INCLUDING, BUT NOT
* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A
* PARTICULAR PURPOSE.  THE SOFTWARE PROVIDED HEREUNDER IS ON AN "AS IS" BASIS,
* AND PURDUE UNIVERSITY HAS NO OBLIGATION TO PROVIDE MAINTENANCE, SUPPORT,
* UPDATES, ENHANCEMENTS, OR MODIFICATIONS.
*/

#include "stdafx.h"

#include "clust_defs.h"
#include "alloc_util.h"
#include "clust_util.h"


int
I_SigSetNClasses(struct SigSet *S)
{
  int i, count;

  for (i = 0, count = 0; i < S->nclasses; i++)
    if (S->ClassSig[i].used) count++;

  return count;
}


struct ClassData *
I_AllocClassData(struct SigSet *S, struct ClassSig *C, int npixels)
{
  struct ClassData *Data;

  Data = (ClassData *)&(C->ClassData);
  Data->npixels = npixels;
  Data->x = G_alloc_matrix (npixels, S->nbands);
  Data->p = G_alloc_matrix (npixels, C->nsubclasses);
  Data->w = G_alloc_vector (npixels*sizeof(double));
  Data->ml_label = G_alloc_ivector (npixels*sizeof(double));
  return Data;
}


void
I_InitSigSet(struct SigSet *S)
{
  S->nbands = 0;
  S->nclasses = 0;
  S->ClassSig = NULL;
  S->title = NULL;
}


void
I_SigSetNBands(struct SigSet *S, int nbands)
{
  S->nbands = nbands;
}


struct ClassSig *
I_NewClassSig(struct SigSet *S)
{
  struct ClassSig *Sp;
  if (S->nclasses == 0)
	  S->ClassSig = (struct SigSet::ClassSig *) G_malloc (sizeof(struct ClassSig));
  else
    S->ClassSig = (struct SigSet::ClassSig *) G_realloc((char *)S->ClassSig,
                  sizeof(struct ClassSig)*(S->nclasses+1));

  Sp = (struct ClassSig *)&S->ClassSig[S->nclasses++];
  Sp->classnum = 0;
  Sp->nsubclasses = 0;
  Sp->used = 1;
  Sp->type = SIGNATURE_TYPE_MIXED;
  Sp->title = NULL;
  Sp->ClassData.npixels = 0;
  Sp->ClassData.SummedWeights = 0.0;
  Sp->ClassData.x = NULL;
  Sp->ClassData.p = NULL;
  Sp->ClassData.w = NULL;
  Sp->ClassData.ml_label = NULL;  
  return Sp;
}


struct SubSig *
I_NewSubSig(struct SigSet *S, struct ClassSig *C)
{
	struct ClassSig::SubSig *Sp;
  int i;

  if (C->nsubclasses == 0)
	  C->SubSig = (struct ClassSig::SubSig *) G_malloc (sizeof(struct SubSig));
  else
    C->SubSig = (struct ClassSig::SubSig *) G_realloc ((char *)C->SubSig,
		sizeof(struct SubSig) * (C->nsubclasses+1));

  Sp = &C->SubSig[C->nsubclasses++];
  Sp->used = 1;
  Sp->R = (double **) G_calloc (S->nbands, sizeof(double *));
  Sp->R[0] = (double *) G_calloc (S->nbands * S->nbands, sizeof(double));
  for (i = 1; i < S->nbands; i++)
    Sp->R[i] = Sp->R[i-1] + S->nbands;
  Sp->Rinv = (double **) G_calloc (S->nbands, sizeof(double *));
  Sp->Rinv[0] = (double *) G_calloc (S->nbands * S->nbands, sizeof(double));
  for (i = 1; i < S->nbands; i++)
    Sp->Rinv[i] = Sp->Rinv[i-1] + S->nbands;
  Sp->means = (double *) G_calloc (S->nbands, sizeof(double));
  Sp->N = 0;
  Sp->pi = 0;
  Sp->cnst = 0;
  return (struct SubSig *)Sp;
}


void
I_SetSigTitle(struct SigSet *S, char *title)
{
  if (title == NULL) title = "";
  if (S->title) free (S->title);
  S->title = G_malloc (strlen (title)+1);
  strcpy(S->title, title);
}


char *
I_GetSigTitle(struct SigSet *S)
{
  if (S->title) return S->title;
  else return "";
}


void
I_SetClassTitle(struct ClassSig *C, char *title)
{
  if (title == NULL) title = "";
  if (C->title) free (C->title);
  C->title = G_malloc (strlen (title)+1);
  strcpy(C->title, title);
}

char *
I_GetClassTitle(struct ClassSig *C)
{
  if (C->title) return C->title;
  else return "";
}


/* Deallocators */

void
I_DeallocClassData(struct SigSet *S, struct ClassSig *C)
{
  struct ClassData * Data;

  Data = (struct ClassData *)&(C->ClassData);
  G_free_matrix(Data->x);
  Data->x = NULL;
  G_free_matrix(Data->p);
  Data->p = NULL;
  G_free_vector(Data->w);
  Data->w = NULL;
  G_free_ivector(Data->ml_label);
  Data->ml_label = NULL;
  Data->npixels = 0;
  Data->SummedWeights = 0.0;
}

    
void
I_DeallocSubSig(struct ClassSig *C)
{
  struct SubSig *Sp;

  Sp = (struct SubSig *)&C->SubSig[--C->nsubclasses];

  G_dealloc( (char *) Sp->R[0] );
  G_dealloc( (char *) Sp->R );
  G_dealloc( (char *) Sp->Rinv[0] );
  G_dealloc( (char *) Sp->Rinv );
  G_dealloc( (char *) Sp->means );

  C->SubSig = (struct ClassSig::SubSig *) G_realloc((char *)C->SubSig,
		    sizeof(struct SubSig) * (C->nsubclasses));
}


void
I_DeallocClassSig(struct SigSet *S)
{
  struct ClassSig *Sp;
    
  Sp = (struct ClassSig *)&(S->ClassSig[--S->nclasses]);
    
  I_DeallocClassData( S, Sp ); 
    
  while(Sp->nsubclasses>0) {
    I_DeallocSubSig(Sp);
  }
  S->ClassSig = (struct SigSet::ClassSig *) G_realloc ((char *)S->ClassSig,
                sizeof(struct ClassSig) * (S->nclasses));
}     
    

void
I_DeallocSigSet(struct SigSet *S)
{
  while(S->nclasses>0) {
    I_DeallocClassSig( S );
  }
}