IrregularHistogram1D Class Reference

#include <histogram.h>

Collaboration diagram for IrregularHistogram1D:

Public Member Functions
void	init (const Histogram1D &oldHistogram, const MultidimArray< int > &bins)
	Initialize class. More...
int	val2Index (double value) const
	Return the index corresponding to a certain value. More...
void	selfNormalize ()
	Normalize to be a probability density function. More...
double	operator() (int i) const
	Get value. More...
const Histogram1D &	getHistogram () const
	Get histogram. More...

Public Attributes
Histogram1D	__hist
MultidimArray< double >	__binsRightLimits

Friends
std::ostream &	operator<< (std::ostream &_out, const IrregularHistogram1D &_h)
	Show. More...

Detailed Description

Irregular histogram.

This class is initialized with a regular histogram and a set of bins. The class builds a new histogram in which the values within bins of the original histogram are gathered. The class provides the standard operations with histograms.

Definition at line 386 of file histogram.h.

Member Function Documentation

getHistogram()

const Histogram1D & IrregularHistogram1D::getHistogram

(

)

const

Get histogram.

Definition at line 467 of file histogram.cpp.

{
    return __hist;
}

init()

void IrregularHistogram1D::init	(	const Histogram1D &	oldHistogram,
		const MultidimArray< int > &	bins
	)

Initialize class.

Definition at line 393 of file histogram.cpp.

{
    int steps_no = XSIZE(bins);
    __binsRightLimits.initZeros(steps_no);
    __hist.initZeros(steps_no);

    int k = 0;
    for (int i = 0; i < steps_no; ++i)
    {
        hist.index2val(A1D_ELEM(bins,i), A1D_ELEM(__binsRightLimits,i));
        for (int j = k; j <= A1D_ELEM(bins,i); ++j)
            A1D_ELEM(__hist,i) += A1D_ELEM(hist,j);
        k = A1D_ELEM(bins,i) + 1;
    }
}

operator()()

double IrregularHistogram1D::operator() ( int  i ) const

inline

Get value.

Definition at line 406 of file histogram.h.

    {
        return DIRECT_A1D_ELEM(__hist,i);
    }

selfNormalize()

void IrregularHistogram1D::selfNormalize

(

)

Normalize to be a probability density function.

Definition at line 451 of file histogram.cpp.

{
    __hist *= 1.0/__hist.sum();
}

val2Index()

int IrregularHistogram1D::val2Index

(

double

value

)

const

Return the index corresponding to a certain value.

Definition at line 411 of file histogram.cpp.

{
    int binsNo = XSIZE(__binsRightLimits);
    /* Binary search is not interesting for small vectors
     * However, we leave the code, just in case it is needed in the future
     int ileft=1, iright=binsNo-1;
     int iret=0;
     if (value<=DIRECT_A1D_ELEM(__binsRightLimits,0))
     return =0;
     else if (value>DIRECT_A1D_ELEM(__binsRightLimits,iright))
     return iright;
     else
     do
     {
     iret = (iright+ileft)/2;
     if (value>DIRECT_A1D_ELEM(__binsRightLimits,iret))
     ileft=iret+1;
     else
     iright=iret;
     if (value>DIRECT_A1D_ELEM(__binsRightLimits,iret-1) &&
     value<=DIRECT_A1D_ELEM(__binsRightLimits,iret))
     break;
     else if (ileft == iright)
     {
     iret=ileft;
     break;
     }
     }
     while (true);
     */

    for (int i = 0; i < binsNo; ++i)
        if (value <= DIRECT_A1D_ELEM(__binsRightLimits,i))
            return i;

    //In case the value is greater, we return the last bin
    return binsNo - 1;
}

Friends And Related Function Documentation

operator<<

std::ostream& operator<< ( std::ostream &  _out, const IrregularHistogram1D &  _h  )

friend

Show.

Definition at line 457 of file histogram.cpp.

{
    for (size_t i = 0; i < XSIZE(_hist.__binsRightLimits); i++)
        _out << "\t" << _hist.__binsRightLimits(i) << "\t\t" << _hist.__hist(i)
        << std::endl;
    return _out;
}

Member Data Documentation

__binsRightLimits

MultidimArray<double> IrregularHistogram1D::__binsRightLimits

Definition at line 390 of file histogram.h.

__hist

Histogram1D IrregularHistogram1D::__hist

Definition at line 389 of file histogram.h.

---

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

• xmippCore/core/histogram.h
• xmippCore/core/histogram.cpp