ossimHistogram.cpp File Reference

#include <ossim/base/ossimCommon.h>
#include <ossim/base/ossimHistogram.h>
#include <ossim/base/ossimNotifyContext.h>
#include <ossim/base/ossimThinPlateSpline.h>
#include <ossim/base/ossimDpt.h>
#include <cmath>
#include <cstdio>
#include <fstream>
#include <iostream>
#include <iomanip>
#include <sstream>

Go to the source code of this file.

Macros
#define	MAX(x, y)   ((x)>(y)?(x):(y))
#define	MIN(x, y)   ((x)>(y)?(y):(x))

Functions
	RTTI_DEF1 (ossimHistogram, "ossimHistogram", ossimObject)
float	GetExtendedCount (int bin, int n_bins, float *cnts, bool cyclic)
void	RemoveFlatPeaks (int nbins, float *cnts, bool cyclic)

Macro Definition Documentation

MAX

#define MAX	(		x,
			y
	)		((x)>(y)?(x):(y))

Definition at line 39 of file ossimHistogram.cpp.

Referenced by ossimHistogram::create(), ossimChipperUtil::getIrect(), ossimPixelFlipper::getMaxPixelValue(), ossimPixelFlipper::getMinPixelValue(), and ossimFfL7::readRadiomRecord().

MIN

#define MIN	(		x,
			y
	)		((x)>(y)?(y):(x))

Definition at line 40 of file ossimHistogram.cpp.

Referenced by ossimHistogram::create(), ossimPixelFlipper::getMaxPixelValue(), and ossimPixelFlipper::getMinPixelValue().

Function Documentation

GetExtendedCount()

float GetExtendedCount ( int  bin, int  n_bins, float *  cnts, bool  cyclic  )

inline

Definition at line 620 of file ossimHistogram.cpp.

Referenced by ossimHistogram::NonMaximumSupress(), and RemoveFlatPeaks().

{
   int nbm = n_bins-1;
   if(!cyclic)
   {
      if(bin < 0)
         return cnts[0];
      if(bin >= n_bins)
         return cnts[nbm];
   }
   else
   {
      if(bin<0)
         return cnts[nbm+bin];
      if(bin >= n_bins)
         return cnts[bin-n_bins];
   }
   return cnts[bin];
}

RemoveFlatPeaks()

void RemoveFlatPeaks ( int  nbins, float *  cnts, bool  cyclic  )

inline

Definition at line 643 of file ossimHistogram.cpp.

References GetExtendedCount().

Referenced by ossimHistogram::NonMaximumSupress().

{
   int nbm = nbins-1;

   //Here we define a small state machine - parsing for runs of peaks
   //init is the state corresponding to an initial run (starting at i ==0)
   bool init=(GetExtendedCount(0, nbins, cnts, cyclic) !=0 ) ? true : false;
   int init_end =0;

   //start is the state corresponding to any other run of peaks
   bool start=false;  
   int start_index=0; 
   int i = 0;

   //The scan of the state machine
   for(i = 0; i < nbins; i++)
   {
      float v = GetExtendedCount(i, nbins, cnts, cyclic);

      //State init: a string of non-zeroes at the begining.
      if(init&&v!=0)
         continue;

      if(init&&v==0)
      {
         init_end = i;
         init = false;
         continue;
      }

      //State !init&&!start: a string of "0s"
      if(!start&&v==0)
         continue;

      //State !init&&start: the first non-zero value
      if(!start&&v!=0)
      {
         start_index = i;
         start = true;
         continue;
      }
      //State ending flat peak: encountered a subsequent zero after starting
      if(start&&v==0)
      {
         int peak_location = (start_index+i-1)/2;//The middle of the run
         int k = 0;
         for(k = start_index; k<=(i-1); k++)
        if(k!=peak_location)
               cnts[k] = 0;
         start = false;
      }
   }
   //Now handle the boundary conditions
   //The non-cyclic case
   if(!cyclic)
   {
      if(init_end!=0)  //Was there an initial run of peaks?
      {
         int init_location = (init_end-1)/2;
         int k = 0;
         for(k = 0; k<init_end; k++)
        if(k!=init_location)
               cnts[k] = 0;
      }
      if(start)       // Did we reach the end of the array in a run of pks?
      {
         int end_location = (start_index + nbm)/2;
         int k = 0;
         for(k = start_index; k<nbins; k++)
        if(k!=end_location)
               cnts[k] = 0;
      }
   }
   else  //The cyclic case
   {
      if(init_end!=0)  //Is there a run which crosses the cyclic cut?
      {
         if(start)    
         { //Yes, so define the peak location accordingly
        int peak_location = (start_index + init_end - nbm -1)/2;
        int k;
        if(peak_location < 0) //Is the peak to the left of the cut?
            {// Yes, to the left
               peak_location += nbm; 
               for( k = 0; k< init_end; k++)
          cnts[k]=0;
               for( k= start_index; k <nbins; k++)
          if(k!=peak_location)
                     cnts[k] = 0;
            }
        else   
            {//No, on the right.
               for( k = start_index; k< nbins; k++)
          cnts[k]=0;
               for( k= 0; k < init_end; k++)
          if(k!=peak_location)
                     cnts[k] = 0;
            }
         }
         else  
         {//There wasn't a final run so just clean up the initial run
        int init_location = (init_end-1)/2;
        int k = 0;
        for(k = start_index; k<init_end; k++)
               if(k!=init_location)
                  cnts[k] = 0;
         }
      }
   }
}

RTTI_DEF1()

RTTI_DEF1	(	ossimHistogram	,
		"ossimHistogram"	,
		ossimObject
	)