#include <ossimHistogram.h>

Inheritance diagram for ossimHistogram:

Classes
class	ossimProprietaryHeaderInformation

Public Types
enum	FillAlgorithmType { HISTOGRAM_FILL_DEFAULT = 0, HISTOGRAM_FILL_THIN_PLATE = 1 }

Public Member Functions
	ossimHistogram ()

	ossimHistogram (int xres, float min, float max)

	ossimHistogram (float , float , int)

	ossimHistogram (const ossimHistogram &his)

	ossimHistogram (const ossimHistogram *, float width)

	ossimHistogram (const double *samples, ossim_uint32 size, ossim_uint32 numBins)
	Uses samples array to establish a histogram with numBins: More...

virtual int	GetIndex (float) const

ossimHistogram *	fillInteriorEmptyBins (int type=HISTOGRAM_FILL_THIN_PLATE) const

ossimHistogram *	Scale (float scale_factor)

ossimHistogram *	CumulativeGreaterThanEqual () const

ossimHistogram *	CumulativeLessThanEqual () const

ossimHistogram *	NonMaximumSupress (int radius=1, bool cyclic=false)

void	create (int xres, float val1, float val2)

void	UpCount (float newval, float occurences=1)

float	GetCount (float uval) const

float	SetCount (float pixelval, float count)

float	GetMinVal () const

float	GetMaxVal () const

float	GetMaxCount () const

float	GetRangeMin () const

float	GetRangeMax () const

float *	GetVals ()

const float *	GetVals () const

float *	GetCounts ()

const float *	GetCounts () const

int	GetRes () const

float	GetBucketSize () const

float *	GetMinValAddr ()

float *	GetMinCountAddr ()

const float *	GetMinValAddr () const

const float *	GetMinCountAddr () const

float	ComputeArea (float low, float high) const

float	ComputeArea () const

float	getLowFractionFromValue (float val) const

float	getHighFractionFromValue (float val) const

float	LowClipVal (float clip_fraction) const

float	HighClipVal (float clip_fraction) const

float	GetValFromIndex (ossim_uint32 idx) const

float	GetMinValFromIndex (ossim_uint32 idx) const

float	GetMaxValFromIndex (ossim_uint32 idx) const

int	GetValIndex (float val) const

float	GetMean () const

float	GetStandardDev () const

void	Print () const

void	Dump (char *) const

int	WritePlot (const char *fname) const

virtual	~ossimHistogram ()

virtual bool	importHistogram (const ossimFilename &inputFile)

virtual bool	importHistogram (istream &in)

virtual bool	saveState (ossimKeywordlist &kwl, const char *prefix=0) const

virtual bool	loadState (const ossimKeywordlist &kwl, const char *prefix=0)

virtual bool	saveState (ossimRefPtr< ossimXmlNode > xmlNode) const

virtual bool	loadState (const ossimRefPtr< ossimXmlNode > xmlNode)

Public Member Functions inherited from ossimObject
	ossimObject ()

virtual	~ossimObject ()

virtual ossimObject *	dup () const

virtual ossimString	getShortName () const

virtual ossimString	getLongName () const

virtual ossimString	getDescription () const

virtual ossimString	getClassName () const

virtual RTTItypeid	getType () const

virtual bool	canCastTo (ossimObject *obj) const

virtual bool	canCastTo (const RTTItypeid &id) const

virtual bool	canCastTo (const ossimString &parentClassName) const

virtual std::ostream &	print (std::ostream &out) const
	Generic print method. More...

virtual bool	isEqualTo (const ossimObject &obj, ossimCompareType compareType=OSSIM_COMPARE_FULL) const

virtual void	accept (ossimVisitor &visitor)

Public Member Functions inherited from ossimReferenced
	ossimReferenced ()

	ossimReferenced (const ossimReferenced &)

ossimReferenced &	operator= (const ossimReferenced &)

void	ref () const
	increment the reference count by one, indicating that this object has another pointer which is referencing it. More...

void	unref () const
	decrement the reference count by one, indicating that a pointer to this object is referencing it. More...

void	unref_nodelete () const
	decrement the reference count by one, indicating that a pointer to this object is referencing it. More...

int	referenceCount () const

Protected Member Functions
virtual void	deleteAll ()

Protected Member Functions inherited from ossimReferenced
virtual	~ossimReferenced ()

Protected Attributes
float *	m_vals

float *	m_counts

int	m_num

float	m_delta

float	m_vmin

float	m_vmax

float	m_mean

float	m_standardDev

Private Attributes
int	m_statsConsistent

Detailed Description

Definition at line 28 of file ossimHistogram.h.

Member Enumeration Documentation

◆ FillAlgorithmType

enum ossimHistogram::FillAlgorithmType

Enumerator
HISTOGRAM_FILL_DEFAULT
HISTOGRAM_FILL_THIN_PLATE

Definition at line 79 of file ossimHistogram.h.

    {
       HISTOGRAM_FILL_DEFAULT    = 0,
       HISTOGRAM_FILL_THIN_PLATE = 1
    };

Constructor & Destructor Documentation

◆ ossimHistogram() [1/6]

ossimHistogram::ossimHistogram ( )

Definition at line 46 of file ossimHistogram.cpp.

Referenced by CumulativeGreaterThanEqual(), CumulativeLessThanEqual(), fillInteriorEmptyBins(), NonMaximumSupress(), and Scale().

    :
    m_statsConsistent(MEAN_FLAG | SD_FLAG),
    m_vals(new float [1]),
    m_counts(new float [1]),
    m_num(0),
    m_delta(0.0),
    m_vmin(0),
    m_vmax(0),
    m_mean(0.0),
    m_standardDev(0.0)
 {
    m_vals[0] = 0.0;
    m_counts[0] = 0.0;
 }

◆ ossimHistogram() [2/6]

ossimHistogram::ossimHistogram	(	int	xres,
		float	min,
		float	max
	)

Definition at line 62 of file ossimHistogram.cpp.

    :
    m_statsConsistent(MEAN_FLAG | SD_FLAG),
    m_vals(new float [xres]),
    m_counts(new float [xres]),
    m_num(xres),
    m_delta(0.0),
    m_vmin(0),
    m_vmax(0),
    m_mean(0.0),
    m_standardDev(0.0)
 {
    m_vmax = MAX(val1, val2);
    m_vmin = MIN(val1, val2);
 
    //---
    // Set the delta which is used to index the bins.
    // Note: that using "(m_vmax - m_vmin) / xres" was dropping the
    // last bin on integer data.
    //---
    if ( (m_vmax - m_vmin + 1) == xres )
    {
       m_delta = 1.0;
    }
    else
    {
       m_delta = (m_vmax - m_vmin) / xres;
    }
 
    m_mean = (float)((m_vmax + m_vmin)/2.0);
    m_standardDev = (float)((m_vmax - m_vmin)/(2.0*sqrt(3.0)));
    int i = 0;
 
    if (m_vals == NULL || m_counts == NULL)
    {
       fprintf(stderr, "Histogram : Ran out of memory for arrays.\n");
       m_vals = NULL;
       m_counts = NULL;
       m_num = 0;
       m_vmin = 0;
       m_vmax = 0;
       m_delta = 0.0;
    }
    else
    {
       //std::cout << std::setprecision(15) << m_vmin << ", " << m_vmax <<  ", " <<m_delta <<", "<< xres << std::endl; 
       for(i = 0; i < xres; i++)
       {
          m_vals[i] = m_vmin + m_delta * (float)(i + 0.5);
          //std::cout << m_vals[i] << std::endl;
          m_counts[i] = 0.0;
       }
    }
 }

◆ ossimHistogram() [3/6]

ossimHistogram::ossimHistogram	(	float *	uvals,
		float *	ucounts,
		int	xres
	)

Definition at line 117 of file ossimHistogram.cpp.

    :
    m_statsConsistent(MEAN_FLAG | SD_FLAG),
    m_vals(uvals),
    m_counts(ucounts),
    m_num(xres),
    m_delta(0.0),
    m_vmin(0),
    m_vmax(0),
    m_mean(0.0),
    m_standardDev(0.0)
 {
    if ( ( xres >= 2 ) && uvals && ucounts )
    {
       m_delta = m_vals[1] - m_vals[0]; // Changed this from delta = 1.0
       //  m_vmax = GetMaxVal();
       //  m_vmin = GetMinVal(); JAF version
       m_vmin = uvals[0] - .5f*m_delta;
       m_vmax = uvals[m_num-1] + .5f*m_delta;
       m_mean = GetMean();
       m_standardDev = GetStandardDev();
    }   
 }

◆ ossimHistogram() [4/6]

ossimHistogram::ossimHistogram ( const ossimHistogram & his )

Definition at line 186 of file ossimHistogram.cpp.

References GetBucketSize(), GetCounts(), GetMaxVal(), GetMean(), GetMinVal(), GetRes(), GetStandardDev(), GetVals(), m_counts, m_delta, m_mean, m_num, m_standardDev, m_statsConsistent, m_vals, m_vmax, and m_vmin.

 :
 m_statsConsistent(0),
 m_vals(0),
 m_counts(0),
 m_num(0),
 m_delta(0.0),
 m_vmin(0),
 m_vmax(0),
 m_mean(0.0),
 m_standardDev(0.0)
 {
 
    int i = 0;
    m_num = his.GetRes();
 
    m_vals = new float[m_num];
    const float* his_vals = his.GetVals();
 
    m_counts = new float[m_num];
    const float* his_counts = his.GetCounts();
 
    if (m_vals == NULL || m_counts == NULL)
    {
       fprintf(stderr, "Histogram : Ran out of memory for arrays.\n");
       m_vals = NULL;
       m_counts = NULL;
       m_num = 0;
       m_vmin = 0;
       m_vmax = 0;
       m_delta = 0.0;
       m_statsConsistent = 0;
       return;
    }
 
    m_mean = his.GetMean();
    m_standardDev = his.GetStandardDev();
 
    for(i=0; i<m_num; i++)
    {
       m_vals[i] = his_vals[i];
       m_counts[i] = his_counts[i];
    }
    m_vmax = his.GetMaxVal();
    m_vmin = his.GetMinVal();
    m_delta = his.GetBucketSize();
 
    m_statsConsistent = 0;
    m_statsConsistent |= (MEAN_FLAG | SD_FLAG);
 }

◆ ossimHistogram() [5/6]

ossimHistogram::ossimHistogram	(	const ossimHistogram *	his,
		float	width
	)

Definition at line 241 of file ossimHistogram.cpp.

References GetBucketSize(), GetRes(), GetVals(), m_counts, m_delta, m_mean, m_num, m_standardDev, m_statsConsistent, m_vals, m_vmax, and m_vmin.

 :
 m_statsConsistent(0),
 m_vals(0),
 m_counts(0),
 m_num(0),
 m_delta(0.0),
 m_vmin(0),
 m_vmax(0),
 m_mean(0.0),
 m_standardDev(0.0)
 {
 
    m_statsConsistent =0;
 
 // Attributes of original histogram
 
    float del = his->GetBucketSize();
    int max_index = his->GetRes() - 1;
    float minvalue = his->GetVals()[0] - del*.5f;
    float maxvalue = his->GetVals()[max_index] + del*.5f;
 
 
 // Intialize a new histogram
    if(width == del) m_num = his->GetRes();
    else if(!(width == 0.0))
       m_num = (int)ceil((maxvalue - minvalue)/width);
    else
       m_num = 1; // This shouldn't happen anyway.
 
    m_vals = new float [m_num];
    m_counts = new float [m_num];
    m_delta = width;
    float mean_val = (maxvalue + minvalue)/2.0f;
    float half_range = (m_num * m_delta)/2.0f;
    m_vmax =  mean_val + half_range;
    m_vmin =  mean_val - half_range;
    int i = 0;
 
    if (m_vals == NULL || m_counts == NULL)
    {
       fprintf(stderr,
               "Histogram : Ran out of memory for arrays.\n");
       m_vals = NULL;
       m_counts = NULL;
       m_num = 0;
       m_vmin = 0;
       m_vmax = 0;
       m_delta = 0.0;
       m_mean = 0.0;
       m_standardDev = 0.0;
       m_statsConsistent |= (MEAN_FLAG | SD_FLAG);
       return;
        
    }
    
    else
    {
       for(i = 0; i < m_num; i++)
       {
          m_vals[i] = m_vmin + m_delta * (i + 0.5f);
          m_counts[i] = 0.0;
       }
    }
 
 
 // Cases:
 
 
    if(width == del)    // Then just copy his
    {
       const float* his_counts = his->GetCounts();
       for(i=0; i<m_num; i++)
          m_counts[i] = his_counts[i];
       m_mean = GetMean();
       m_standardDev = GetStandardDev();
       m_statsConsistent |= (MEAN_FLAG | SD_FLAG);
       return;
    }
 
 
    if(del > width)     // Then interpolate his m_counts.
    {
 
 // Boundary conditions:
 //    Start
       float his_start = minvalue + .5f*del;
       float start = m_vmin + .5f*m_delta;
       float c0 = his->GetCount(his_start);
       float c1 = his->GetCount(his_start + del);
       float s0 = (c1 - c0)/del;
 
       for(float x = start; x <= (his_start + del + m_delta);)
       {
          float interp = s0 * (x - his_start) + c0;
          if(interp < 0) interp = 0; //Can be negative
          SetCount(x,interp);
          x += width;
       }
 //    End
       float his_end = maxvalue - .5f*del;
       float end = m_vmax - .5f*m_delta;
       float cn = his->GetCount(his_end);
       float cn_1 = his->GetCount(his_end - del);
       float sn = (cn_1 - cn)/del;
 
       for(float y = end; y >= (his_end - del + m_delta);)
       {
          float interp = sn * (his_end - y) + cn;
          if(interp < 0) interp = 0; //Can be negative
          SetCount(y, interp);
          y -= m_delta;
       }
 // Interior Loop
 
       for(float z = his_start + del; z <= (his_end - del);)
       {
          float ci = his->GetCount(z);
          float ci_1 = his->GetCount(z-del);
          float cip1 = his->GetCount(z+del);
          float deriv = (cip1 - ci_1)/(2.0f*del);
          float second_drv =
             ((cip1 + ci_1)/2.0f - ci)/(del*del);
          int fine_x_index = GetIndex(z);
          if (fine_x_index < 0)
          {
             if (z<m_vmin) fine_x_index = 0;
             else fine_x_index = m_num-1;
          }
          float fine_x = m_vals[fine_x_index];
          for(float xfine = fine_x; xfine < z + del;)
          {
             float interp = ci + deriv*(xfine -z) +
                second_drv*(xfine - z)*(xfine - z);
 
             if(interp < 0) interp = 0; //Can be negative
             SetCount(xfine, interp);
             xfine += width;
          }
          z += del;
       }
    }
 
 
    if(del < width)    //Just accumulate samples from his into larger bins
    {
       if( del != 0.0){
          float his_start = minvalue + .5f*del;
          float his_end = maxvalue - .5f*del;
          for(float x = his_start; x <= his_end;)
          {
             SetCount(x, (GetCount(x) + his->GetCount(x)));
             x += del;
          }
       }
    }
    m_mean = GetMean();
    m_standardDev = GetStandardDev();
    m_statsConsistent =0;
    m_statsConsistent |= (MEAN_FLAG | SD_FLAG);
 }

◆ ossimHistogram() [6/6]

ossimHistogram::ossimHistogram	(	const double *	samples,
		ossim_uint32	size,
		ossim_uint32	numBins
	)

Uses samples array to establish a histogram with numBins:

Definition at line 140 of file ossimHistogram.cpp.

References GetMean(), GetStandardDev(), m_counts, m_delta, m_num, m_vals, m_vmax, m_vmin, size, and UpCount().

    :
    m_statsConsistent(0),
    m_vals(0),
    m_counts(0),
    m_num((int)xres),
    m_delta(0.0),
    m_vmin(0),
    m_vmax(0),
    m_mean(0.0),
    m_standardDev(0.0)
 {
    if ((size == 0) || (xres == 0))
       return;
 
    // scan the dataset for min/max:
    m_vmin=(float)(data[0]);
    m_vmax=(float)(data[0]);
    for (ossim_uint32 i=1; i<size; ++i)
    {
       if ((float)(data[i]) < m_vmin)
          m_vmin = (float)(data[i]);
       else if ((float)(data[i]) > m_vmax)
          m_vmax = (float)(data[i]);
    }
 
    // Allocate histogram:
    m_delta = (m_vmax - m_vmin) / m_num;
    m_vals = new float [m_num];
    m_counts = new float [m_num];
    for (ossim_int32 i=0; i<m_num; ++i)
    {
       m_vals[i] = m_vmin + m_delta * (i + 0.5);
       m_counts[i] = 0.0;
    }
 
    // compute histogram:
    for (ossim_uint32 i=0; i<size; i++)
       UpCount((float)(data[i]));
 
    GetMean();
    GetStandardDev();
 }

◆ ~ossimHistogram()

ossimHistogram::~ossimHistogram ( )

virtual

Definition at line 1300 of file ossimHistogram.cpp.

References deleteAll().

 {
    deleteAll();
 }

Member Function Documentation

◆ ComputeArea() [1/2]

float ossimHistogram::ComputeArea	(	float	low,
		float	high
	)		const

Definition at line 1051 of file ossimHistogram.cpp.

References GetIndex(), GetMaxVal(), GetMinVal(), m_counts, m_num, and m_vmin.

 {
    float sum = 0.0;
    float maxval = GetMaxVal();
    float minval = GetMinVal();
 
    if (low < minval) low = minval;
    if (high > maxval) high = maxval;
 
    if (low <= high)
    {
       int indexlow, indexhigh;
       indexlow = (int) GetIndex(low);
       if (indexlow < 0)
       {
          if (low<m_vmin) indexlow = 0;
          else indexlow = m_num-1;
       }
       indexhigh = (int) GetIndex(high);
       if (indexhigh < 0)
       {
          if (high<m_vmin) indexhigh = 0;
          else indexhigh = m_num-1;
       }
       int i=indexlow;
 
       while (i<=indexhigh)
       {
          sum+= m_counts[i];
          i++;
       }
    }
 
    return sum;
 }

◆ ComputeArea() [2/2]

float ossimHistogram::ComputeArea ( ) const

Definition at line 1089 of file ossimHistogram.cpp.

References GetMaxVal(), GetMinVal(), m_vmax, and m_vmin.

Referenced by HighClipVal(), LowClipVal(), and Print().

 {
    float m_vmin = this->GetMinVal();
    float m_vmax = this->GetMaxVal();
    if(m_vmin>m_vmax)
    {
       float temp = m_vmin;
       m_vmin = m_vmax;
       m_vmax = temp;
    }
    return this->ComputeArea(m_vmin, m_vmax);
 }

◆ create()

void ossimHistogram::create	(	int	xres,
		float	val1,
		float	val2
	)

Definition at line 403 of file ossimHistogram.cpp.

References deleteAll(), m_counts, m_delta, m_mean, m_num, m_standardDev, m_statsConsistent, m_vals, m_vmax, m_vmin, MAX, and MIN.

Referenced by importHistogram(), and loadState().

 {
    // clear all the data
    deleteAll();
 
    // now set it up and initialize;
    xres = xres >0? xres:1;
    
    m_vals   = new float [xres];
    m_counts = new float [xres];
    m_num = xres;
    m_vmax = MAX(val1, val2);
    m_vmin = MIN(val1, val2);
 
    m_delta = (m_vmax - m_vmin) / xres;
    m_mean = (float)((m_vmax + m_vmin)/2.0);
    m_standardDev = (float)((m_vmax - m_vmin)/(2.0*sqrt(3.0)));
    m_statsConsistent = 0;
    m_statsConsistent |= (MEAN_FLAG | SD_FLAG);
    int i = 0;
    if (m_vals == NULL || m_counts == NULL)
    {
       ossimNotify(ossimNotifyLevel_FATAL) << "Histogram : Ran out of memory for arrays.\n";
       m_vals = NULL;
       m_counts = NULL;
       m_num = 0;
       m_vmin = 0;
       m_vmax = 0;
       m_delta = 0.0;
    }
    else
    {
       for(i = 0; i < xres; i++)
       {
          m_vals[i] = m_vmin + m_delta * (float)(i + 0.5);
          m_counts[i] = 0.0;
       }
    }   
 }

◆ CumulativeGreaterThanEqual()

ossimHistogram * ossimHistogram::CumulativeGreaterThanEqual ( ) const

Definition at line 573 of file ossimHistogram.cpp.

References GetCounts(), GetRes(), and ossimHistogram().

 {
    ossimHistogram* cum_his = new ossimHistogram(*this);
    const float* density_counts = this->GetCounts();
    int res = this->GetRes();
 
    // Intitialize cumulative m_counts
    float* cum_counts = cum_his->GetCounts();
    int i = 0;
    for(i=0; i < res; i++)
       cum_counts[i] = 0;
   
    cum_counts[res-1] = density_counts[res-1];
    for(i = res-2; i>=0; --i)
    {
       cum_counts[i] += (density_counts[i] + cum_counts[i+1]);
    }
 
    return cum_his;
 }

◆ CumulativeLessThanEqual()

ossimHistogram * ossimHistogram::CumulativeLessThanEqual ( ) const

Definition at line 594 of file ossimHistogram.cpp.

References GetCounts(), GetRes(), and ossimHistogram().

 {
    ossimHistogram* cum_his = new ossimHistogram(*this);
    const float* density_counts = this->GetCounts();
    int res = this->GetRes();
 
    // Intitialize cumulative m_counts
    float* cum_counts = cum_his->GetCounts();
    int i = 0;
    for(i=0; i < res; i++)
       cum_counts[i] = 0;
 
    cum_counts[0] = density_counts[0];
    for(i = 1; i < res; i++)
    {
       cum_counts[i] += (density_counts[i] + cum_counts[i-1]);
    }
   
    return cum_his;
 }

◆ deleteAll()

void ossimHistogram::deleteAll ( )

protectedvirtual

Definition at line 1286 of file ossimHistogram.cpp.

References m_counts, and m_vals.

Referenced by create(), and ~ossimHistogram().

 {
    if (m_vals)
    {
       delete []m_vals;
       m_vals = NULL;
    }
    if (m_counts)
    {
       delete []m_counts;
       m_counts = NULL;
    }  
 }

◆ Dump()

void ossimHistogram::Dump ( char * dumpfile ) const

Definition at line 1248 of file ossimHistogram.cpp.

References m_counts, m_num, and m_vals.

 {
    FILE *dumpfp = fopen(dumpfile, "w");
 
    if (!dumpfp)
    {
       fprintf(stderr, "Error opening histogram data file.\n");
       return;
    }
    int i = 0;
 
    for(i = 0; i < m_num; i++)
       fprintf(dumpfp, "%f %f\n", m_vals[i], m_counts[i]);
   
    fclose(dumpfp);
    return;
 }

◆ fillInteriorEmptyBins()

ossimHistogram * ossimHistogram::fillInteriorEmptyBins ( int type = HISTOGRAM_FILL_THIN_PLATE ) const

Definition at line 442 of file ossimHistogram.cpp.

References ossimThinPlateSpline::addPoint(), GetCounts(), ossimThinPlateSpline::getPoint(), HISTOGRAM_FILL_DEFAULT, HISTOGRAM_FILL_THIN_PLATE, m_num, m_statsConsistent, ossimHistogram(), and ossimThinPlateSpline::solve().

 {
    if(m_num < 1) return 0;
    ossimHistogram* result = new ossimHistogram(*this);
    switch(type)
    {
       case HISTOGRAM_FILL_THIN_PLATE:
       case HISTOGRAM_FILL_DEFAULT:
       {
          ossimThinPlateSpline spline(1);
          double pvars[1];
          float* new_counts = result->GetCounts();
          ossim_int32 idxLeft = 0;
          ossim_int32 idxRight = m_num-1;
          while((idxLeft < m_num) && (new_counts[idxLeft]  < 1))++idxLeft;
          while((idxRight > -1) && (new_counts[idxRight] < 1))--idxRight;
          if(idxLeft < idxRight)
          {
             ossim_int32 idx = idxLeft;
             while(idx <= idxRight)
             {
                if(new_counts[idx]>0)
                {
                   pvars[0] = new_counts[idx];
                   spline.addPoint(idx, 0, pvars);
                }
                ++idx;
             }
             if(spline.solve())
             {
                idx = idxLeft;
                while(idx <= idxRight)
                {
                   if(spline.getPoint(idx, 0, pvars))
                   {
                      new_counts[idx] = pvars[0];
                   }
                   ++idx;
                }
                m_statsConsistent = 0;
             }
             else
             {
             }
          }
          
          break;
       }
    }
    
    return result;
 }

◆ GetBucketSize()

float ossimHistogram::GetBucketSize ( ) const

inline

Definition at line 149 of file ossimHistogram.h.

Referenced by ossimHistogram(), and Print().

149 { return m_delta; }

ossimHistogram::m_delta

float m_delta

Definition: ossimHistogram.h:47

◆ GetCount()

float ossimHistogram::GetCount ( float uval ) const

Definition at line 971 of file ossimHistogram.cpp.

References GetIndex(), and m_counts.

 {
    int index = GetIndex(pixelval);
 
    if (index < 0)
       return -1;
    else
       return m_counts[index];
 }

◆ GetCounts() [1/2]

float* ossimHistogram::GetCounts ( )

inline

Definition at line 134 of file ossimHistogram.h.

Referenced by ossimShorelineTool::autoComputeThreshold(), ossimKMeansFilter::computeKMeans(), CumulativeGreaterThanEqual(), CumulativeLessThanEqual(), fillInteriorEmptyBins(), HighClipVal(), ossimHistogramEqualization::initializeLuts(), loadState(), LowClipVal(), NonMaximumSupress(), ossimHistogram(), ossimImageData::populateHistogram(), Print(), and Scale().

     { 
       m_statsConsistent = 0; // m_counts might change.
       return m_counts; 
     }

◆ GetCounts() [2/2]

const float* ossimHistogram::GetCounts ( ) const

inline

Definition at line 140 of file ossimHistogram.h.

     { 
       //m_statsConsistent = 0; // m_counts might change.
       return m_counts; 
     }

◆ getHighFractionFromValue()

float ossimHistogram::getHighFractionFromValue ( float val ) const

Returns the fraction of accumulation down to and including "val" bucket from max divided by the total accumulation. returns OSSIM_FLT_NAN if "val" is not between GetMin() and GetMax().

Definition at line 1133 of file ossimHistogram.cpp.

References GetMaxVal(), GetMinVal(), GetRes(), GetValIndex(), m_counts, and ossim::nan().

Referenced by ossimHistogramRemapper::setHighClipPoint().

 {
 //   float min = floor(GetMinVal());
 //   float max = ceil(GetMaxVal());
    float minValue = GetMinVal();
    float maxValue = GetMaxVal();
    if (val < minValue || val > maxValue)
    {
       return ossim::nan();
    }
 
    int total_buckets = GetRes();
    int cutoff_bucket = GetValIndex(val);
    float partial_sum = 0.0;
    float total_sum   = 0.0;
    
    for(int i = (total_buckets-1); i >= 0; --i)
    {
       total_sum += m_counts[i];
       if (i >= cutoff_bucket)
       {
          partial_sum += m_counts[i];
       }
    }
 
    return (partial_sum/total_sum);
 }

◆ GetIndex()

int ossimHistogram::GetIndex ( float pixelval ) const

virtual

Definition at line 865 of file ossimHistogram.cpp.

References GetRes(), m_delta, m_num, m_vals, m_vmax, and m_vmin.

Referenced by ComputeArea(), GetCount(), ossimHistogramEqualization::initializeLuts(), ossimHistogramEqualization::runEqualizationAlgorithm(), Scale(), SetCount(), and UpCount().

 {
 #if 1
    if ((pixelval > m_vmax) || (pixelval < m_vmin) || (m_num==0) )
    {
       return -1;
    }
 //   ossim_float32 d = m_vmax-m_vmin;
    int bandIdx = (ossim_int32)((pixelval-m_vmin)/m_delta);
    return bandIdx<GetRes()?bandIdx:-1;
 //    if(bandIdx == m_num)
 //    {
 //       return m_num-1;
 //    }
 //    else if(bandIdx < m_num)
 //    {
 //       return bandIdx;
 //    }
 //    return -1;
 #else
    if ((pixelval > m_vmax) || (pixelval < m_vmin))
       return -1;
 
    int idx = 0;
    int i = 0;
 
    for(i = 0; i < m_num; i++)
    {
       //std::cout << std::setprecision(15) << m_vals[i] << std::endl;
       // RWMC: This is very dangerous - might get an intermediate
       // value which is between m_vals[i]+0.5*m_delta and
       // m_vals[i+1]-0.5*m_delta, which would then return index of 0.
       // Changed to check range one-sided, which is safe because of
       // previous check on range.
       //       if ((pixelval > (m_vals[i] - 0.5 * m_delta)) &&
       //           (pixelval <= (m_vals[i] + 0.5 * m_delta)))
       if (pixelval <= (m_vals[i] + 0.5 * m_delta))
       {
          idx = i;
          break;
       }
    }
 //std::cout << idx << std::endl;
    return idx;
 #endif
 }

◆ getLowFractionFromValue()

float ossimHistogram::getLowFractionFromValue ( float val ) const

Returns the fraction of accumulation up to and including "val" bucket from min divided by the total accumulation. returns OSSIM_FLT_NAN if "val" is not between GetMinVal and GetMaxVal.

Definition at line 1102 of file ossimHistogram.cpp.

References GetMaxVal(), GetMinVal(), GetRes(), GetValIndex(), m_counts, and ossim::nan().

Referenced by ossimHistogramRemapper::setLowClipPoint().

 {
    // std::cout << "ossimHistogram::getLowFractionFromValue(float val)\n";
 //   float minValue = floor(GetMinVal());
 //   float maxValue = ceil(GetMaxVal());
    float minValue = GetMinVal();
    float maxValue = GetMaxVal();
    if (val < minValue || val > maxValue)
    {
       return ossim::nan();
    }
 // std::cout << "VAL: " << val << "\n"
 //           << "MIN: " << minValue << "\n"
 //           << "MAX: " << maxValue << "\n"; 
    int total_buckets = GetRes();
    int cutoff_bucket = GetValIndex(val);
    float partial_sum = 0.0;
    float total_sum   = 0.0;
    // std::cout << "CUTOFF BUCKET ===" << cutoff_bucket << "\n";
    for(int i = 0; i < total_buckets; ++i)
    {
       total_sum += m_counts[i];
       if (i <= cutoff_bucket)
       {
          partial_sum += m_counts[i];
       }
    }
    // std::cout << "FRACTION ==== " << (partial_sum/total_sum) << "\n";
    return (partial_sum/total_sum);
 }

◆ GetMaxCount()

float ossimHistogram::GetMaxCount ( ) const

Definition at line 1010 of file ossimHistogram.cpp.

References m_counts, m_num, and max.

Referenced by ossimHistogramEqualization::initializeLuts().

 {
    int i=0;
    float max;
    max = 0.0;
    for (i=0; i < m_num; i++)
       if (m_counts[i] > max)
          max = m_counts[i];
    return max;
 }

◆ GetMaxVal()

float ossimHistogram::GetMaxVal ( ) const

Definition at line 996 of file ossimHistogram.cpp.

References m_counts, m_num, and m_vals.

Referenced by ComputeArea(), getHighFractionFromValue(), getLowFractionFromValue(), HighClipVal(), LowClipVal(), ossimHistogram(), and Print().

 {
    int i=m_num-1;
 
    while (i>0 && !m_counts[i])
       i--;
 
    if (i < 0)
       return 0.0;
 
    return m_vals[i];
 }

◆ GetMaxValFromIndex()

float ossimHistogram::GetMaxValFromIndex ( ossim_uint32 idx ) const

Definition at line 924 of file ossimHistogram.cpp.

References m_delta, m_num, m_vals, and m_vmax.

 {
    float result = 0.0;
    
    if((int)idx < m_num)
    {
       result = (m_vals[idx]+(0.5 * m_delta));
       if(result > m_vmax) result = m_vmax;
    }
 
    return result;  
 }

◆ GetMean()

float ossimHistogram::GetMean ( ) const

Definition at line 802 of file ossimHistogram.cpp.

Referenced by ossimHistogram(), and Scale().

 {
    float xsum = 0.0;
 
    if(MEAN_FLAG&m_statsConsistent)
       return m_mean;
    else
    {
       if( this->GetBucketSize() > 0.0){
          for(float x=this->GetMinVal(); x<= this->GetMaxVal(); x +=this->GetBucketSize())
             xsum += x*GetCount(x);
       }
 
       float area = ComputeArea(m_vmin, m_vmax);
       if(area <= 0.0)
       {
          //       fprintf(stderr, "Histogram : Area <= 0.0\n");
          return 0.0;
       }
       else
       {
          m_statsConsistent |=1;
          m_mean = xsum/area;
          return m_mean;
       }
    }
 }

◆ GetMinCountAddr() [1/2]

float* ossimHistogram::GetMinCountAddr ( )

inline

Definition at line 154 of file ossimHistogram.h.

155 { return m_counts+GetIndex(GetMinVal()); }

ossimHistogram::GetIndex

virtual int GetIndex(float) const

Definition: ossimHistogram.cpp:865

ossimHistogram::GetMinVal

float GetMinVal() const

Definition: ossimHistogram.cpp:983

ossimHistogram::m_counts

float * m_counts

Definition: ossimHistogram.h:41

◆ GetMinCountAddr() [2/2]

const float* ossimHistogram::GetMinCountAddr ( ) const

inline

Definition at line 160 of file ossimHistogram.h.

161 { return m_counts+GetIndex(GetMinVal()); }

ossimHistogram::GetIndex

virtual int GetIndex(float) const

Definition: ossimHistogram.cpp:865

ossimHistogram::GetMinVal

float GetMinVal() const

Definition: ossimHistogram.cpp:983

ossimHistogram::m_counts

float * m_counts

Definition: ossimHistogram.h:41

◆ GetMinVal()

float ossimHistogram::GetMinVal ( ) const

Definition at line 983 of file ossimHistogram.cpp.

References m_counts, m_num, and m_vals.

Referenced by ComputeArea(), getHighFractionFromValue(), getLowFractionFromValue(), HighClipVal(), LowClipVal(), ossimHistogram(), and Print().

 {
    int i=0;
 
    while (i<m_num-1 && !m_counts[i])
       i++;
 
    return m_vals[i];
 }

◆ GetMinValAddr() [1/2]

float* ossimHistogram::GetMinValAddr ( )

inline

Definition at line 151 of file ossimHistogram.h.

152 { return m_vals+GetIndex(GetMinVal()); }

ossimHistogram::GetIndex

virtual int GetIndex(float) const

Definition: ossimHistogram.cpp:865

ossimHistogram::GetMinVal

float GetMinVal() const

Definition: ossimHistogram.cpp:983

ossimHistogram::m_vals

float * m_vals

Definition: ossimHistogram.h:39

◆ GetMinValAddr() [2/2]

const float* ossimHistogram::GetMinValAddr ( ) const

inline

Definition at line 157 of file ossimHistogram.h.

158 { return m_vals+GetIndex(GetMinVal()); }

ossimHistogram::GetIndex

virtual int GetIndex(float) const

Definition: ossimHistogram.cpp:865

ossimHistogram::GetMinVal

float GetMinVal() const

Definition: ossimHistogram.cpp:983

ossimHistogram::m_vals

float * m_vals

Definition: ossimHistogram.h:39

◆ GetMinValFromIndex()

float ossimHistogram::GetMinValFromIndex ( ossim_uint32 idx ) const

Definition at line 911 of file ossimHistogram.cpp.

References m_delta, m_num, m_vals, and m_vmin.

 {
    float result = 0.0;
    
    if((int)idx < m_num)
    {
       result = (m_vals[idx]-(0.5 * m_delta));
       if(result < m_vmin) result = m_vmin;
    }
 
    return result;
 }

◆ GetRangeMax()

float ossimHistogram::GetRangeMax ( ) const

inline

Definition at line 119 of file ossimHistogram.h.

Referenced by ossimKMeansFilter::initialize(), and ossimHistogramEqualization::initializeLuts().

     {
        return m_vmax;
     }

◆ GetRangeMin()

float ossimHistogram::GetRangeMin ( ) const

inline

Definition at line 115 of file ossimHistogram.h.

Referenced by ossimKMeansFilter::initialize(), and ossimHistogramEqualization::initializeLuts().

     {
        return m_vmin;
     }

◆ GetRes()

int ossimHistogram::GetRes ( ) const

inline

◆ GetStandardDev()

float ossimHistogram::GetStandardDev ( ) const

Definition at line 832 of file ossimHistogram.cpp.

Referenced by ossimHistogram(), and Scale().

 {
    float sum = 0.0;
 
    if(SD_FLAG&m_statsConsistent)
       return m_standardDev;
    else
    {
       float xm = this -> GetMean(); // Force an Update of m_mean
 
       if( this->GetBucketSize() > 0.0){
          for(float x=this->GetMinVal();
              x<= this->GetMaxVal();
              x +=this->GetBucketSize())
 
             sum += (x-xm)*(x-xm)*GetCount(x);
       }
 
       float area = ComputeArea(m_vmin, m_vmax);
       if(area <= 0.0)
       {
          //       fprintf(stderr, "Histogram : Area <= 0.0\n");
          return 0.0;
       }
       else
       {
          m_statsConsistent |= 2;
          m_standardDev = (float)sqrt(sum/area);
          return m_standardDev;
       }
    }
 }

◆ GetValFromIndex()

float ossimHistogram::GetValFromIndex ( ossim_uint32 idx ) const

Definition at line 937 of file ossimHistogram.cpp.

References m_num, and m_vals.

 {
    float result = 0.0;
    if((int)idx < m_num)
    {
       result = m_vals[idx];
    }
 
    return result;
 }

◆ GetValIndex()

int ossimHistogram::GetValIndex ( float val ) const

Definition at line 948 of file ossimHistogram.cpp.

References m_delta, m_num, m_vals, m_vmax, and m_vmin.

Referenced by getHighFractionFromValue(), and getLowFractionFromValue().

 {
    if ((pixelval > m_vmax) || (pixelval < m_vmin))
       return -1;
 
    int idx = 0;
    int i = 0;
 
    for(i = 0; i < m_num; i++)
    {
       if ((pixelval > (m_vals[i] - 0.5 * m_delta)) &&
           (pixelval <= (m_vals[i] + 0.5 * m_delta)))
       {
          idx = i;
          break;
       }
    }
 
    return idx;
 }

◆ GetVals() [1/2]

float* ossimHistogram::GetVals ( )

inline

Definition at line 123 of file ossimHistogram.h.

Referenced by ossimShorelineTool::autoComputeThreshold(), ossimKMeansFilter::computeKMeans(), HighClipVal(), LowClipVal(), ossimHistogram(), and Print().

     {
       m_statsConsistent = 0; // Values might change.
       return m_vals; 
     }

◆ GetVals() [2/2]

const float* ossimHistogram::GetVals ( ) const

inline

Definition at line 128 of file ossimHistogram.h.

     {
       m_statsConsistent = 0; // Values might change.
       return m_vals; 
     }

◆ HighClipVal()

float ossimHistogram::HighClipVal ( float clip_fraction ) const

Definition at line 1194 of file ossimHistogram.cpp.

References ComputeArea(), GetCounts(), GetMaxVal(), GetMinVal(), GetRes(), GetVals(), m_counts, and m_vals.

Referenced by ossimHistogramRemapper::getHighClipPoint(), ossimMultiBandHistogramTileSource::runLinearStretchAlgorithm(), and ossimHistogramThreshholdFilter::runThreshholdStretchAlgorithm().

 {
    if(clip_fraction<0) clip_fraction=0.0;
    if(clip_fraction>1.0) clip_fraction=1.0;
    float area = this->ComputeArea();
    if(area==0.0) return this->GetMaxVal();
    if(clip_fraction==0.0) return this->GetMaxVal();
    if(clip_fraction==1.0) return this->GetMinVal();
    float clip_area = area*clip_fraction;
    const float* m_counts = this->GetCounts();
    const float* m_vals = this->GetVals();
    int res = this->GetRes();
    float sum = 0;
    int i = (res-1);
    for(; i>=0; i--)
    {
       sum+=m_counts[i];
       if(sum>=clip_area)
          break;
    }
    return m_vals[i];
 }

◆ importHistogram() [1/2]

bool ossimHistogram::importHistogram ( const ossimFilename & inputFile )

virtual

Definition at line 1352 of file ossimHistogram.cpp.

References ossimString::c_str(), and ossimFilename::exists().

Referenced by ossimMultiBandHistogram::importHistogram().

 {
    if(inputFile.exists())
    {
       ifstream input(inputFile.c_str());
 
       return importHistogram(input);
    }
 
    return false;  
 }

◆ importHistogram() [2/2]

bool ossimHistogram::importHistogram ( istream & in )

virtual

Definition at line 1306 of file ossimHistogram.cpp.

References ossimString::c_str(), create(), getline(), ossimHistogram::ossimProprietaryHeaderInformation::getNumberOfBins(), ossimHistogram::ossimProprietaryHeaderInformation::parseStream(), SetCount(), and ossimString::toDouble().

 {
    ossimProprietaryHeaderInformation header;
    bool binsCreated = false;
    
    if(header.parseStream(in))
    {
       long numberOfBins = header.getNumberOfBins();
       
       if(numberOfBins)
       {
          create(numberOfBins, 0, numberOfBins - 1);
          binsCreated = true;
 
          if(binsCreated)
          {
             ossimString buffer;
             ossimString binNumber;
             ossimString count;
 
             while(in.good() &&
                   !in.eof() &&
                   *binNumber.c_str() != '.')
             {
                
                getline(in, buffer);
 
                istringstream s(buffer);
 
                s >> binNumber >> count;
                if(*binNumber.c_str() != (char)'.')
                {
                   SetCount((float)binNumber.toDouble(),
                            (float)count.toDouble());
                }
             }
          }
       }
       else
       {
          return false;
       }
    }  
    return true;
 }

◆ loadState() [1/2]

bool ossimHistogram::loadState	(	const ossimKeywordlist &	kwl,
		const char *	prefix = `0`
	)

virtual

Method to the load (recreate) the state of the object from a keyword list. Return true if ok or false on error.

Reimplemented from ossimObject.

Definition at line 1514 of file ossimHistogram.cpp.

References create(), ossimString::empty(), ossimKeywordlist::extractKeysThatMatch(), ossimKeywordlist::find(), GetCounts(), ossimKeywordlist::getMap(), ossimKeywordlist::getSubstringKeyList(), size, ossimString::toDouble(), ossimString::toString(), ossimString::toUInt32(), and ossim::toVector().

Referenced by ossimMultiBandHistogram::loadState().

 {
 //   std::cout << "ossimHistogram::loadState!!!!\n";
    const char* number_of_bins = kwl.find(prefix, "number_of_bins");
 
 //   std::cout << "NBINS = " << number_of_bins << std::endl;
    if(number_of_bins)
    {
       ossim_uint32 bins = ossimString(number_of_bins).toUInt32();
 
 //      std::cout << "BINS ======== " << bins << std::endl;
       if(bins > 0)
       {
          // setup some defaults
          float minValue = 0;
          float maxValue = bins - 1;
 
          // see if there is a range set for the data
          const char* min_value = kwl.find(prefix, "min_value");
          const char* max_value = kwl.find(prefix, "max_value");
 
          if(min_value)
          {
             minValue = (ossim_float32)ossimString(min_value).toDouble();
          }
          if(max_value)
          {
             maxValue = (ossim_float32)ossimString(max_value).toDouble();
          }
 
          create((int)bins, minValue, maxValue);
          float* countsPtr = GetCounts();
          memset(countsPtr, '\0', bins*sizeof(float));
          // this is new style histogram creation
          //
          ossimString binsString = kwl.find(prefix, "bins");
          if(!binsString.empty())
          {
             std::vector<ossimDpt> result;
             ossim::toVector(result, binsString);
             if(!result.empty())
             {
                ossim_uint32 idx = 0;
                for(idx = 0; idx < result.size();++idx)
                {
                   ossim_uint32 binIdx = static_cast<ossim_uint32>(result[idx].x);
                   if(binIdx < bins)
                   {
                      countsPtr[binIdx] = result[idx].y;
                   }
                }
             }
          }
          else
          {
             ossimKeywordlist binsKwl;
             ossim_uint32 offset = (ossim_uint32)(ossimString(prefix)+"bin").size();
             ossimString regExpression =  ossimString("^(") + ossimString(prefix) + "bin[0-9]+)";
             kwl.extractKeysThatMatch(binsKwl,regExpression);
             const ossimKeywordlist::KeywordMap& kwlMap = binsKwl.getMap();
             ossimKeywordlist::KeywordMap::const_iterator iter = kwlMap.begin();
             while(iter != kwlMap.end())
             {
                ossimString numberStr(iter->first.begin() + offset,
                                      iter->first.end());
                countsPtr[numberStr.toUInt32()] = ossimString(iter->second).toDouble();
                ++iter;
             }
          }
 
 //         ossimKeywordlist kwl;
 //         this->saveState(kwl);
 //         std::cout << kwl << std::endl;
          
          return true;
 #if 0
          // create the bins
          ossimString binNumber = "";
          ossimString regExpression =  ossimString("^(") + ossimString(prefix) + "bin[0-9]+)";
          vector<ossimString> keys = kwl.getSubstringKeyList( regExpression );
          ossim_uint32 numberOfBins = (ossim_uint32)keys.size();
          ossim_uint32 offset = (ossim_uint32)(ossimString(prefix)+"bin").size();
 
          std::vector<ossim_uint32> theNumberList(numberOfBins);
          ossim_uint32 idx = 0;
          for(idx = 0; idx < theNumberList.size();++idx)
          {
             ossimString numberStr(keys[idx].begin() + offset,
                                   keys[idx].end());
             theNumberList[idx] = numberStr.toUInt32();
 
          }
          
          float* countsPtr = GetCounts();
          memset(countsPtr, '\0', bins*sizeof(float));
          for(idx = 0; idx < numberOfBins;++idx)
          {
             const char* binCount = kwl.find(prefix, ossimString("bin") + ossimString::toString(theNumberList[idx]));
             countsPtr[theNumberList[idx]] = (float)ossimString(binCount).toDouble();
          }
 #endif
       }
    }
    return true;
 }

◆ loadState() [2/2]

bool ossimHistogram::loadState ( const ossimRefPtr< ossimXmlNode > xmlNode )

virtual

Definition at line 1621 of file ossimHistogram.cpp.

References create(), ossimXmlNode::findFirstNode(), GetCounts(), ossimXmlNode::getText(), ossimString::push_back(), ossimString::toFloat32(), and ossimRefPtr< T >::valid().

 {
    ossimRefPtr<ossimXmlNode> binValues =  xmlNode->findFirstNode("binValues");
    ossimRefPtr<ossimXmlNode> minValueNode  =  xmlNode->findFirstNode("minValue");
    ossimRefPtr<ossimXmlNode> maxValueNode  =  xmlNode->findFirstNode("maxValue");
 
    if(binValues.valid())
    {
       ossim_uint32 count = 0;
       float minValue = 0.0;
       float maxValue = 0.0;
       std::vector<float> floatValues;
       std::istringstream in(binValues->getText());
       ossimString vString;
       while(!in.fail())
       {
          in>>vString;
          if(!in.fail())
          {
             floatValues.push_back(vString.toFloat32());
          }
       }
       count = (ossim_uint32)floatValues.size();
      
       if(count)
       {
          minValue = 0;
          maxValue = count - 1;
 
          if(minValueNode.valid())
          {
             minValue = minValueNode->getText().toFloat32();
          }
          if(maxValueNode.valid())
          {
             maxValue = maxValueNode->getText().toFloat32();
          }
 
          create(count, minValue, maxValue);
          float* countsPtr = GetCounts();
          ossim_uint32 idx = 0;
          for(idx = 0; idx < count; ++idx)
          {
             countsPtr[idx] = floatValues[idx];
          }
          return true;
       }
    }
 
    return false;
 }

◆ LowClipVal()

float ossimHistogram::LowClipVal ( float clip_fraction ) const

Definition at line 1165 of file ossimHistogram.cpp.

References ComputeArea(), GetCounts(), GetMaxVal(), GetMinVal(), GetRes(), GetVals(), m_counts, and m_vals.

Referenced by ossimHistogramRemapper::getLowClipPoint(), ossimMultiBandHistogramTileSource::runLinearStretchAlgorithm(), and ossimHistogramThreshholdFilter::runThreshholdStretchAlgorithm().

 {
    if(clip_fraction<0) clip_fraction=0.0;
    if(clip_fraction>1.0) clip_fraction=1.0;
    float area = this->ComputeArea();
    if(area==0.0) return this->GetMinVal();
    if(clip_fraction==0.0) return this->GetMinVal();
    if(clip_fraction==1.0) return this->GetMaxVal();
    float clip_area = area*clip_fraction;
    const float* m_counts = this->GetCounts();
    const float* m_vals = this->GetVals();
    int res = this->GetRes();
    float sum = 0;
    int i=0;
 
    for(; i<res; i++)
    {
       sum+=m_counts[i];
       if(sum>=clip_area)
          break;
    }
 
    return m_vals[i];
 }

◆ NonMaximumSupress()

ossimHistogram * ossimHistogram::NonMaximumSupress	(	int	radius = `1`,
		bool	cyclic = `false`
	)

Definition at line 762 of file ossimHistogram.cpp.

References GetCounts(), GetExtendedCount(), GetRes(), m_num, ossimHistogram(), ossimNotify(), ossimNotifyLevel_WARN, and RemoveFlatPeaks().

 {
    if((2*radius +1)> m_num/2)
    {
       ossimNotify(ossimNotifyLevel_WARN)<<"ossimHistogram::NonMaximumSupress the radius is too large \n";
       return NULL;
    }
    //Get the m_counts array of "this"
    ossimHistogram* h_new = new ossimHistogram(*this);
    int n_buckets = h_new->GetRes();
    float* counts_old = this->GetCounts();
 
    //Make a new Histogram for the suppressed version
    float* counts_new = h_new->GetCounts();
    int i;
    for( i =0; i < n_buckets; i++)
       counts_new[i] = 0;
   
    //Find local maxima
    for( i = 0; i<  n_buckets; i++)
    {
       //find the maxium value in the current kernel
       float max_count = counts_old[i];
       int k = 0;
       for(k = -radius; k <= radius ;k++)
       {
          int index = i+k;
          float c = GetExtendedCount(index, n_buckets, counts_old, cyclic);
          if( c > max_count)
         max_count = c;
       }
       //Is position i a local maxium?
       if(max_count == counts_old[i])
          counts_new[i] = max_count;//Yes. So set the m_counts to the max value
    }
    RemoveFlatPeaks(n_buckets, counts_new, cyclic);
    return h_new;
 }

◆ Print()

void ossimHistogram::Print ( ) const

Definition at line 1219 of file ossimHistogram.cpp.

References ComputeArea(), GetBucketSize(), GetCounts(), GetMaxVal(), GetMinVal(), GetRes(), GetVals(), m_counts, m_vals, ossimNotify(), and ossimNotifyLevel_INFO.

 {
    ostream& out = ossimNotify(ossimNotifyLevel_INFO);
    const float* m_vals = this->GetVals();
    const float* m_counts = this->GetCounts();
    int res = this->GetRes();
    int width = 0;
    int i = 0;
    for(i =0; i < res; i++)
    {
       if(width++ > 5)
       {
          width = 0;
          out << "\n";
       }
       out << m_vals[i] << " " << m_counts[i] << " | " ;
    }
    out << "\n MaxVal " << this->GetMaxVal() << "\n";
    out << " MinVal " << this->GetMinVal() << "\n";
    out << " BucketSize " << this->GetBucketSize() << "\n";
    out << " Resolution " << this->GetRes() << "\n";
    out << " Area "
        << this->ComputeArea(this->GetMinVal(),this->GetMaxVal()) << "\n";
    out << "------------------------------------------------\n\n";
 }

◆ saveState() [1/2]

bool ossimHistogram::saveState	(	ossimKeywordlist &	kwl,
		const char *	prefix = `0`
	)		const

virtual

Method to save the state of the object to a keyword list. Return true if ok or false on error.

Reimplemented from ossimObject.

Definition at line 1449 of file ossimHistogram.cpp.

References ossimKeywordlist::add(), ossimString::c_str(), FLT_EPSILON, m_counts, m_num, m_vmax, m_vmin, and ossimString::toString().

 {
    kwl.add(prefix,
            "type",
            "ossimHistogram",
            true);
    kwl.add(prefix,
            "number_of_bins",
            m_num,
            true);
    kwl.add(prefix,
            "min_value",
            m_vmin,
            true);
    kwl.add(prefix,
            "max_value",
            m_vmax,
            true);
    
 
 
    ossimString binArrayList = "(";
    bool firstValue = true;
 
    for(ossim_int32 index = 0; index < m_num; ++index)
    {
       if(fabs(m_counts[index]) > FLT_EPSILON)
       {
 
          if(!firstValue)
          {
             binArrayList += ",";
          }
          else
          {
             firstValue = false;
          }
          binArrayList += "("+ossimString::toString(index)+","+ossimString::toString(m_counts[index])+")";
      }
    }
 
    binArrayList += ")";
 
    kwl.add(prefix, "bins", binArrayList, true);
 #if 0
    ossimString binValue = "";
    for(ossim_int32 index = 0; index < m_num; ++index)
    {
       if(fabs(m_counts[index]) > FLT_EPSILON)
       {
          //     binValue = prefix;
          binValue = "bin";
          binValue += ossimString::toString(index);
          
          kwl.add(prefix,
                  binValue.c_str(),
                  m_counts[index],
                  true);
       }
    }
 #endif   
    return true;
 }

◆ saveState() [2/2]

bool ossimHistogram::saveState ( ossimRefPtr< ossimXmlNode > xmlNode ) const

virtual

Definition at line 1673 of file ossimHistogram.cpp.

References ossimXmlNode::addChildNode(), ossimRefPtr< T >::get(), m_counts, m_mean, m_num, m_standardDev, m_vmax, m_vmin, ossimXmlNode::setTag(), ossimXmlNode::setText(), and ossimString::toString().

 {
    ossimRefPtr<ossimXmlNode> binValues = new ossimXmlNode;
    xmlNode->setTag("ossimHistogram");
    xmlNode->addChildNode("minValue", ossimString::toString(m_vmin));
    xmlNode->addChildNode("maxValue", ossimString::toString(m_vmax));
    xmlNode->addChildNode("standardDeviation", ossimString::toString(m_standardDev));
    xmlNode->addChildNode("m_mean", ossimString::toString(m_mean));
    binValues->setTag("binValues");
    std::ostringstream out;
 
    ossim_int32 idx = 0;
    if(m_num > 0)
    {
       for(idx = 0; idx < m_num;++idx)
       {
          out << ossimString::toString(m_counts[idx], 8) << " ";
       }
       binValues->setText(out.str());
    }
    xmlNode->addChildNode(binValues.get());
    
    return true;
 }

◆ Scale()

ossimHistogram * ossimHistogram::Scale ( float scale_factor )

Definition at line 499 of file ossimHistogram.cpp.

References GetCounts(), GetIndex(), GetMean(), GetStandardDev(), m_counts, m_delta, m_mean, m_num, m_standardDev, m_vals, m_vmin, ossimHistogram(), and x.

 {
 
 // Extract attributes of self
 
 //    float lowvalue = m_vals[0];
    float highvalue = m_vals[m_num-1];
 
 // Construct a new histogram
 
    ossimHistogram* scaled_his = new ossimHistogram(this, m_delta);
    float* new_counts = scaled_his->GetCounts();
    int i = 0;
    for(i=0; i < m_num; i++)  // Initialize
       new_counts[i] = 0.0;
 
 // Compute scaled values
 // We assume that the new histogram is to be scaled down from his
 
    float scale = scale_factor;
    if(scale_factor > 1.0) scale = 1.0;
 
    for(float x = highvalue; x > m_vmin;)
    {
       float trans_x = (x-m_vmin)*scale + m_vmin; // Scaled x.
       int index = GetIndex(trans_x);
       if (index < 0)
       {
          if (trans_x<m_vmin) index = 0;
          else index = m_num-1;
       }
       float fraction = (trans_x - m_vals[index])/m_delta;
       float abs_fraction = (float)fabs(fraction);
       int x_index = GetIndex(x);
       if (x_index < 0)
       {
          if (x<m_vmin) x_index = 0;
          else x_index = m_num-1;
       }
 
 // Distribute the m_counts in proportion
 
       new_counts[index] += (1.0f - abs_fraction)*m_counts[x_index];
       if(fraction > 0)
          if(index < (m_num-1))
             new_counts[index + 1] +=
                abs_fraction*m_counts[x_index];
          else
             new_counts[index] +=
                abs_fraction*m_counts[x_index];
       else
          if(index > 0)
             new_counts[index - 1] +=
                abs_fraction*m_counts[x_index];
          else
             new_counts[index] +=
                abs_fraction*m_counts[x_index];
       x -= m_delta;
    }
 
 // Compute new Histogram attributes
 
    m_mean = scaled_his->GetMean();
    m_standardDev = scaled_his->GetStandardDev();
    return scaled_his;
 }

◆ SetCount()

float ossimHistogram::SetCount	(	float	pixelval,
		float	count
	)

Definition at line 1024 of file ossimHistogram.cpp.

References GetIndex(), m_counts, and m_statsConsistent.

Referenced by importHistogram(), and ossimHistogramRemapper::setNullCount().

 {
    m_statsConsistent = 0;
 
    int index = GetIndex(pixelval);
 
    if (index < 0)
       return -1;
    else
    {
       m_counts[index] = count;
       return count;
    }
 }

◆ UpCount()

void ossimHistogram::UpCount	(	float	newval,
		float	occurences = `1`
	)

Definition at line 1040 of file ossimHistogram.cpp.

References GetIndex(), m_counts, and m_statsConsistent.

Referenced by ossimHistogram(), and ossimImageData::populateHistogram().

 {
 
    m_statsConsistent = 0;
    int idx = GetIndex(pixelval);
    if (idx >= 0)  // Originally (index > 0)
    {
       m_counts[idx] += occurences;
    }
 }

◆ WritePlot()

int ossimHistogram::WritePlot ( const char * fname ) const

Definition at line 1269 of file ossimHistogram.cpp.

References m_counts, m_num, and m_vals.

 {
    FILE *fp = fopen(fname, "w");
 
    if (!fp)
    {
       fprintf(stderr, "Error opening histogram plot file.\n");
       return 0;
    }
 
    for(int j = 0; j < m_num; j++)
       fprintf(fp, "%f %f\n", m_vals[j], m_counts[j]);
 
    fclose(fp);
    return 1;
 }

Member Data Documentation

◆ m_counts

float* ossimHistogram::m_counts

protected

Definition at line 41 of file ossimHistogram.h.

Referenced by ComputeArea(), create(), deleteAll(), Dump(), GetCount(), getHighFractionFromValue(), getLowFractionFromValue(), GetMaxCount(), GetMaxVal(), GetMinVal(), HighClipVal(), LowClipVal(), ossimHistogram(), Print(), saveState(), Scale(), SetCount(), UpCount(), and WritePlot().

◆ m_delta

float ossimHistogram::m_delta

protected

Definition at line 47 of file ossimHistogram.h.

Referenced by create(), GetIndex(), GetMaxValFromIndex(), GetMinValFromIndex(), GetValIndex(), ossimHistogram(), and Scale().

◆ m_mean

float ossimHistogram::m_mean

mutableprotected

Definition at line 49 of file ossimHistogram.h.

Referenced by create(), ossimHistogram(), saveState(), and Scale().

◆ m_num

int ossimHistogram::m_num

protected

Definition at line 45 of file ossimHistogram.h.

Referenced by ComputeArea(), create(), Dump(), fillInteriorEmptyBins(), GetIndex(), GetMaxCount(), GetMaxVal(), GetMaxValFromIndex(), GetMinVal(), GetMinValFromIndex(), GetValFromIndex(), GetValIndex(), NonMaximumSupress(), ossimHistogram(), saveState(), Scale(), and WritePlot().

◆ m_standardDev

float ossimHistogram::m_standardDev

mutableprotected

Definition at line 50 of file ossimHistogram.h.

Referenced by create(), ossimHistogram(), saveState(), and Scale().

◆ m_statsConsistent

int ossimHistogram::m_statsConsistent

mutableprivate

Definition at line 32 of file ossimHistogram.h.

Referenced by create(), fillInteriorEmptyBins(), ossimHistogram(), SetCount(), and UpCount().

◆ m_vals

float* ossimHistogram::m_vals

protected

Definition at line 39 of file ossimHistogram.h.

Referenced by create(), deleteAll(), Dump(), GetIndex(), GetMaxVal(), GetMaxValFromIndex(), GetMinVal(), GetMinValFromIndex(), GetValFromIndex(), GetValIndex(), HighClipVal(), LowClipVal(), ossimHistogram(), Print(), Scale(), and WritePlot().

◆ m_vmax

float ossimHistogram::m_vmax

protected

Definition at line 48 of file ossimHistogram.h.

Referenced by ComputeArea(), create(), GetIndex(), GetMaxValFromIndex(), GetValIndex(), ossimHistogram(), and saveState().

◆ m_vmin

float ossimHistogram::m_vmin

protected

Definition at line 48 of file ossimHistogram.h.

Referenced by ComputeArea(), create(), GetIndex(), GetMinValFromIndex(), GetValIndex(), ossimHistogram(), saveState(), and Scale().

The documentation for this class was generated from the following files:

Classes

Public Types

Public Member Functions

Protected Member Functions

Protected Attributes

Private Attributes

Detailed Description

Member Enumeration Documentation

◆ FillAlgorithmType

Constructor & Destructor Documentation

◆ ossimHistogram() [1/6]

◆ ossimHistogram() [2/6]

◆ ossimHistogram() [3/6]

◆ ossimHistogram() [4/6]

◆ ossimHistogram() [5/6]

◆ ossimHistogram() [6/6]

◆ ~ossimHistogram()

Member Function Documentation

◆ ComputeArea() [1/2]

◆ ComputeArea() [2/2]

◆ create()

◆ CumulativeGreaterThanEqual()

◆ CumulativeLessThanEqual()

◆ deleteAll()

◆ Dump()

◆ fillInteriorEmptyBins()

◆ GetBucketSize()

◆ GetCount()

◆ GetCounts() [1/2]

◆ GetCounts() [2/2]

◆ getHighFractionFromValue()

◆ GetIndex()

◆ getLowFractionFromValue()

◆ GetMaxCount()

◆ GetMaxVal()

◆ GetMaxValFromIndex()

◆ GetMean()

◆ GetMinCountAddr() [1/2]

◆ GetMinCountAddr() [2/2]

◆ GetMinVal()

◆ GetMinValAddr() [1/2]

◆ GetMinValAddr() [2/2]

◆ GetMinValFromIndex()

◆ GetRangeMax()

◆ GetRangeMin()

◆ GetRes()

◆ GetStandardDev()

◆ GetValFromIndex()

◆ GetValIndex()

◆ GetVals() [1/2]

◆ GetVals() [2/2]

◆ HighClipVal()

◆ importHistogram() [1/2]

◆ importHistogram() [2/2]

◆ loadState() [1/2]

◆ loadState() [2/2]

◆ LowClipVal()

◆ NonMaximumSupress()

◆ Print()

◆ saveState() [1/2]

◆ saveState() [2/2]

◆ Scale()

◆ SetCount()

◆ UpCount()

◆ WritePlot()

Member Data Documentation

◆ m_counts

◆ m_delta

◆ m_mean

◆ m_num

◆ m_standardDev

◆ m_statsConsistent

◆ m_vals

◆ m_vmax

◆ m_vmin