casacore::LatticeStatistics< T > Class Template Reference

Compute and display various statistics from a lattice. More...

#include <LatticeStatistics.h>

Inheritance diagram for casacore::LatticeStatistics< T >:

Classes
struct	AlgConf
Public Types
typedef NumericTraits< T > ::PrecisionType	AccumType
typedef std::pair< String, String >	stat_element
typedef std::list< stat_element >	stat_list
Public Member Functions
	LatticeStatistics (const MaskedLattice< T > &lattice, LogIO &os, Bool showProgress=True, Bool forceDisk=False)
	Constructor takes the lattice and a `LogIO` object for logging.
	LatticeStatistics (const MaskedLattice< T > &lattice, Bool showProgress=True, Bool forceDisk=False)
	Constructor takes the lattice only.
	LatticeStatistics (const LatticeStatistics< T > &other)
	Copy constructor.
virtual	~LatticeStatistics ()
	Destructor.
LatticeStatistics< T > &	operator= (const LatticeStatistics< T > &other)
	Assignment operator.
Bool	setAxes (const Vector< Int > &cursorAxes)
	Set the cursor axes (0 relative).
Bool	setInExCludeRange (const Vector< T > &include, const Vector< T > &exclude, Bool setMinMaxToInclude=False)
	You may specify a pixel intensity range as either one for which all pixels in that range are included or one for which all pixels in that range are excluded.
Bool	setList (const Bool &doList)
	This function allows you to control whether the statistics are written to the output stream if you are also making a plot.
Bool	display ()
	Display the statistics by listing and/or plotting them.
Bool	getLayerStats (String &stats, Double area, Int zAxis=-1, Int zLayer=-1, Int hAxis=-1, Int hLayer=-1)
Bool	getLayerStats (stat_list &stats, Double area, Int zAxis=-1, Int zLayer=-1, Int hAxis=-1, Int hLayer=-1)
Vector< Int >	displayAxes () const
	Return the display axes.
Bool	getStatistic (Array< AccumType > &stat, LatticeStatsBase::StatisticsTypes type, Bool dropDeg=True)
	Recover the desired Statistic into an array.
Bool	getConvertedStatistic (Array< T > &stat, LatticeStatsBase::StatisticsTypes type, Bool dropDeg=True)
Bool	getMinMaxPos (IPosition &minPos, IPosition &maxPos)
	Recover position of min and max.
Bool	getStats (Vector< AccumType > &, const IPosition &pos, const Bool posInLattice=False)
	This function gets a vector containing all the statistics for a given location.
void	resetError ()
	Reset argument error condition.
Bool	getFullMinMax (T &dataMin, T &dataMax)
	Get full lattice min and max only.
String	errorMessage () const
	Recover last error message.
Bool	setNewLattice (const MaskedLattice< T > &lattice)
	Set a new MaskedLattice object.
Bool	hasLogger () const
	Did we construct with a logger ?
void	configureClassical ()
	configure object to use Classical Statistics The time, t_x, it takes to compute classical statistics using algorithm x, can be modeled by t_x = n_sets(a_x + b_xn_el) where n_sets is the number of independent sets of data to compute stats on, each containing n_el number of elements.
void	configureClassical (Double aOld, Double bOld, Double aNew, Double bNew)
void	configureFitToHalf (FitToHalfStatisticsData::CENTER centerType=FitToHalfStatisticsData::CMEAN, FitToHalfStatisticsData::USE_DATA useData=FitToHalfStatisticsData::LE_CENTER, AccumType centerValue=0)
	configure to use fit to half algorithm.
void	configureHingesFences (Double f)
	configure to use hinges-fences algorithm
void	configureChauvenet (Double zscore=-1, Int maxIterations=-1)
	configure to use Chauvenet's criterion
std::map< String, uInt >	getChauvenetNiter () const
	get number of iterations associated with Chauvenet criterion algorithm
Protected Member Functions
virtual Bool	_canDoFlux () const
	Virtual Functions.
virtual Quantum< AccumType >	_flux (Bool &, AccumType, Double) const
virtual void	listMinMax (ostringstream &osMin, ostringstream &osMax, Int oWidth, DataType type)
virtual Bool	listStats (Bool hasBeam, const IPosition &dPos, const Matrix< AccumType > &ord)
	List the statistics to the logger.
virtual Bool	listLayerStats (const Matrix< AccumType > &ord, ostringstream &rslt, Int zLayer)
virtual void	getLabels (String &higherOrderLabel, String &xAxisLabel, const IPosition &dPos) const
*IPosition	locInLattice (const IPosition &storagePosition, Bool relativeToParent=True) const
	Given a location in the storage lattice, convert those locations on the non-statistics axis (the last one) and optionally account for the lattice subsectioning.
void	setStream (ostream &os, Int oPrec)
	Non-virtual functions.
IPosition	_storageLatticeShape () const
	get the storage lattice shape
virtual Bool	_computeFlux (Array< AccumType > &flux, const Array< AccumType > &npts, const Array< AccumType > &sum)
virtual Bool	_computeFlux (Quantum< AccumType > &flux, AccumType sum, const IPosition &pos, Bool posInLattice)
void	_latticePosToStoragePos (IPosition &storagePos, const IPosition &latticePos)
	convert a position in the input lattice to the corresponding position in the stats storage lattice.
Protected Attributes
LogIO	os_p
Vector< Int >	cursorAxes_p
Vector< Int >	displayAxes_p
Bool	goodParameterStatus_p
Bool	haveLogger_p
Bool	fixedMinMax_p
Bool	doRobust_p
	doRobust means that when the storage lattice is generated, the robust statistics are generated as well
Bool	doList_p
IPosition	minPos_p
IPosition	maxPos_p
IPosition	blcParent_p
String	error_p
Private Member Functions
void	_setDefaultCoeffs ()
virtual void	summStats ()
	Summarize the statistics found over the entire lattice.
virtual void	displayStats (AccumType nPts, AccumType sum, AccumType median, AccumType medAbsDevMed, AccumType quartile, AccumType sumSq, AccumType mean, AccumType var, AccumType rms, AccumType sigma, AccumType dMin, AccumType dMax, AccumType q1, AccumType q3)
Bool	calculateStatistic (Array< AccumType > &slice, LatticeStatsBase::StatisticsTypes type, Bool dropDeg)
	Calculate statistic from storage lattice and return in an array.
void	generateRobust ()
	Find the median per cursorAxes chunk.
Bool	generateStorageLattice ()
	Create a new storage lattice.
IPosition	locInStorageLattice (const IPosition &latticePosition, LatticeStatsBase::StatisticsTypes type) const
	Given a location in the lattice and a statistic type, work out where to put it in the storage lattice.
void	minMax (Bool &none, AccumType &dMin, AccumType &dMax, const Vector< AccumType > &d, const Vector< AccumType > &n) const
	Find min and max of good data in arrays specified by pointers.
Bool	retrieveStorageStatistic (Array< AccumType > &slice, const LatticeStatsBase::StatisticsTypes type, const Bool dropDeg)
	Retrieve a statistic from the storage lattice and return in an array.
Bool	retrieveStorageStatistic (Vector< AccumType > &slice, const IPosition &pos, const Bool posInLattice)
	Retrieve a statistic from the storage lattice at the specified location and return in an array.
IPosition	statsSliceShape () const
	Find the shape of slice from the statistics lattice at one spatial pixel.
Bool	someGoodPoints ()
	See if there were some valid points found in the storage lattice.
void	stretchMinMax (AccumType &dMin, AccumType &dMax) const
	Stretch min and max by 5%.
CountedPtr < StatisticsAlgorithm < AccumType, const T , const Bool > >	_createStatsAlgorithm () const
void	_configureDataProviders (LatticeStatsDataProvider< T > &lattDP, MaskedLatticeStatsDataProvider< T > &maskedLattDP) const
void	_doStatsLoop (uInt nsets, CountedPtr< LattStatsProgress > progressMeter)
Private Attributes
const MaskedLattice< T > *	pInLattice_p
CountedPtr< TempLattice < AccumType > >	pStoreLattice_p
Vector< Int >	nxy_p
Vector< Int >	statsToPlot_p
Vector< T >	range_p
Bool	noInclude_p
Bool	noExclude_p
Bool	needStorageLattice_p
Bool	doneSomeGoodPoints_p
Bool	someGoodPointsValue_p
Bool	showProgress_p
Bool	forceDisk_p
T	minFull_p
T	maxFull_p
Bool	doneFullMinMax_p
AlgConf	_algConf
std::map< String, uInt >	_chauvIters
Double	_aOld
Double	_bOld
Double	_aNew
Double	_bNew

Detailed Description

template<class T>
class casacore::LatticeStatistics< T >

Compute and display various statistics from a lattice.

Intended use:

Public interface

Review Status

Date Reviewed:: yyyy/mm/dd

Prerequisite

Etymology

This is a class designed to display and retrieve statistics from lattices

Synopsis

This class enable you to display and/or retrieve statistics evaluated over specified regions of a lattice. The dimension of the region is arbitrary, but the size of each dimension is always the shape of the corresponding lattice axis. The statistics are displayed as a function of location of the axes not used to evaluate the statistics over. The axes which you evaluate the statistics over are called the cursor axes, the others are called the display axes.

For example, consider a lattice cube (call the axes xyz or [0,1,2]). You could display statistics from xy planes (cursor axes [0,1]) as a function of z (display axes [2]). Or you could retrieve statistics from the z axis (cursor axes [2]) for each [x,y] location (display axes [0,1]).

This class inherits from LatticeStatsBase This base class provides an enum defining allowed statistics types and a helper function to convert between a String and a Vector<Int> describing the desired statistics to plot. An example is shown below.

This class can list, plot and retrieve statistics. When it lists statistics, it always lists all the available statistics. When you plot statistics, you must specify which ones you would like to see.

This class generates a "storage lattice" into which it writes the accumulated statistical sums. It is from this storage lattice that the plotting and retrieval arrays are drawn. The storage lattice is either in core or on disk depending upon its size (if > 10% of memory given by .aipsrc system.resources.memory then it goes into a disk-based PagedArray). If on disk, the storage lattice is deleted when the LatticeStatistics class object destructs. However, currently, if the process is terminated ungracefully, the storage lattice will be left over.

Tip: This class has a few virtual functions; they are not part of a nice general polymorphic interface; rather they have specialized functionality; The idea of these is that you can derive a class from LatticeStatistics, such as ImageStatistics which provides you with a little more information when displaying/logging the statistics (such as world coordinates) The virtual functions are

getBeamArea can be used to return the synthesized beam area so that the FLUX statistic can be computed

listStats is used to list the statistics to the logger

getLabelsM find the X-axis label and the title label for the plotting;

Tip: If you ignore return error statuses from the functions that set the state of the class, the internal status of the class is set to bad; This means it will just keep on returning error conditions until you explicitly recover the situation; A message describing the last error condition can be recovered with function errorMessage;

Example

// Construct PagedImage (which isA MaskedLattice) from file name
 
 PagedImage<Float> inImage(inName);
 
// Construct statistics object
 
 LogOrigin or("myClass", "myFunction(...)", WHERE);
 LogIO os(or);
 LatticeStatistics<Float> stats(SubImage<FLoat>(inImage), os);
 
// Set cursor axes to see statistics of yz planes (0 relative)
 
 Vector<Int> cursorAxes(2)
 cursorAxes(0) = 1;
 cursorAxes(1) = 2;
 if (!stats.setAxes(cursorAxes)) return 1;
 
// Set to list and plot mean, sigma and rms
 
 if (!stats.setList(True)) return 1;
 String device = "/xs";
 Vector<Int> nxy(2);
 nxy(0) = 1;
 nxy(1) = 1;
 Vector<Int> statsToPlot = LatticeStatsBase::toStatisticTypes("mean,rms,sigma");
 if (!stats.setPlotting(statsToPlot, device, nxy)) return 1;
 
// Now activate actual listing and plotting
 
 if (!stats.display ()) return 1;
 
// Retrieve statistics into array
 
 Array<Double> sum;
 if (!stats.getStatistic(sum, LatticeStatsBase::SUM)) return 1;

In this example, a PagedImage is constructed (which isA MaskedLattice) with . We set the cursor axes to be the y and z axes, we specify to list the statistics if we plot them, and we ask to plot the mean, standard deviation, and root mean square of each yz plane as a function of x location on the PGPLOT device "/xs" with 1 subplot per page (there will be only one in this case). After the plotting and listing, we also retrieve the sum of the selected pixels as a function of x location into an array.

Motivation

The generation of statistical information from a lattice is a basic and necessary capability.

To Do

Implement plotting for complex lattices
Retrieve statistics at specified location of display axes

Definition at line 196 of file LatticeStatistics.h.

Member Typedef Documentation

template<class T>

typedef NumericTraits<T>::PrecisionType casacore::LatticeStatistics< T >::AccumType

Reimplemented in casacore::ImageStatistics< T >.

Definition at line 202 of file LatticeStatistics.h.

template<class T>

typedef std::pair<String,String> casacore::LatticeStatistics< T >::stat_element

Definition at line 294 of file LatticeStatistics.h.

template<class T>

typedef std::list<stat_element> casacore::LatticeStatistics< T >::stat_list

Definition at line 295 of file LatticeStatistics.h.

Constructor & Destructor Documentation

template<class T>

casacore::LatticeStatistics< T >::LatticeStatistics	(	const MaskedLattice< T > &	lattice,
		LogIO &	os,
		Bool	showProgress = `True`,
		Bool	forceDisk = `False`
	)

Constructor takes the lattice and a LogIO object for logging.

You can specify whether you want to see progress meters or not. You can force the storage lattice to be disk based, otherwise the decision for core or disk is taken for you.

template<class T>

casacore::LatticeStatistics< T >::LatticeStatistics	(	const MaskedLattice< T > &	lattice,
		Bool	showProgress = `True`,
		Bool	forceDisk = `False`
	)

Constructor takes the lattice only.

In the absence of a logger you get no messages. This includes error messages and potential listing of the statistics. You can specify whether you want to see progress meters or not. You can force the storage lattice to be disk based, otherwise the decision for core or disk is taken for you.

template<class T>

casacore::LatticeStatistics< T >::LatticeStatistics ( const LatticeStatistics< T > & other )

Copy constructor.

Copy semantics are followed. Therefore any storage lattice that has already been created for other is copied to *this

template<class T>

virtual casacore::LatticeStatistics< T >::~LatticeStatistics ( ) [virtual]

Destructor.

Member Function Documentation

template<class T>

virtual Bool casacore::LatticeStatistics< T >::_canDoFlux ( ) const [inline, protected, virtual]

Virtual Functions.

See implementation to figure it all out !

FIXME The indirect dependence of this class on ImageInterface related issues (eg flux density) breaks encapsulation. All the ImageInterface related code should be encapsulated in ImageStatistics. Unfortunately, that requires significantly more time than I have atm. A return value of False means that the object in question cannot compute flux density values. The default implementation returns False.

Reimplemented in casacore::ImageStatistics< T >.

Definition at line 424 of file LatticeStatistics.h.

References casacore::False.

template<class T>

virtual Bool casacore::LatticeStatistics< T >::_computeFlux	(	Quantum< AccumType > &	flux,
		AccumType	sum,
		const IPosition &	pos,
		Bool	posInLattice
	)			`[protected, virtual]`

Reimplemented in casacore::ImageStatistics< T >.

template<class T>

virtual Bool casacore::LatticeStatistics< T >::_computeFlux	(	Array< AccumType > &	flux,
		const Array< AccumType > &	npts,
		const Array< AccumType > &	sum
	)			`[protected, virtual]`

Reimplemented in casacore::ImageStatistics< T >.

template<class T>

void casacore::LatticeStatistics< T >::_configureDataProviders	(	LatticeStatsDataProvider< T > &	lattDP,
		MaskedLatticeStatsDataProvider< T > &	maskedLattDP
	)			const `[private]`

template<class T>

CountedPtr<StatisticsAlgorithm<AccumType, const T*, const Bool*> > casacore::LatticeStatistics< T >::_createStatsAlgorithm ( ) const [private]

template<class T>

void casacore::LatticeStatistics< T >::_doStatsLoop	(	uInt	nsets,
		CountedPtr< LattStatsProgress >	progressMeter
	)			`[private]`

template<class T>

virtual Quantum<AccumType> casacore::LatticeStatistics< T >::_flux	(	Bool &	,
		AccumType	,
		Double
	)			const `[inline, protected, virtual]`

Reimplemented in casacore::ImageStatistics< T >.

Definition at line 426 of file LatticeStatistics.h.

References ThrowCc.

template<class T>

void casacore::LatticeStatistics< T >::_latticePosToStoragePos	(	IPosition &	storagePos,
		const IPosition &	latticePos
	)			`[protected]`

convert a position in the input lattice to the corresponding position in the stats storage lattice.

The number of elements in storagePos will not be changed and only the first N elements will be modified where N = the number of elements in latticePos. storagePos must therefore have at least as many elements as latticePos. Returns False if latticePos is inconsistent with the input lattice.

template<class T>

void casacore::LatticeStatistics< T >::_setDefaultCoeffs ( ) [inline, private]

coefficients from timings run on PagedImages on etacarinae.cv.nrao.edu (dmehring's development machine)

Definition at line 512 of file LatticeStatistics.h.

References casacore::LatticeStatistics< T >::_aNew, casacore::LatticeStatistics< T >::_aOld, casacore::LatticeStatistics< T >::_bNew, and casacore::LatticeStatistics< T >::_bOld.

template<class T>

IPosition casacore::LatticeStatistics< T >::_storageLatticeShape ( ) const [inline, protected]

get the storage lattice shape

Definition at line 471 of file LatticeStatistics.h.

References casacore::LatticeStatistics< T >::pStoreLattice_p.

template<class T>

Bool casacore::LatticeStatistics< T >::calculateStatistic	(	Array< AccumType > &	slice,
		LatticeStatsBase::StatisticsTypes	type,
		Bool	dropDeg
	)			`[private]`

Calculate statistic from storage lattice and return in an array.

template<class T>

void casacore::LatticeStatistics< T >::configureChauvenet	(	Double	zscore = `-1`,
		Int	maxIterations = `-1`
	)

configure to use Chauvenet's criterion

template<class T>

void casacore::LatticeStatistics< T >::configureClassical	(	Double	aOld,
		Double	bOld,
		Double	aNew,
		Double	bNew
	)

template<class T>

void casacore::LatticeStatistics< T >::configureClassical ( )

configure object to use Classical Statistics The time, t_x, it takes to compute classical statistics using algorithm x, can be modeled by t_x = n_sets*(a_x + b_x*n_el) where n_sets is the number of independent sets of data to compute stats on, each containing n_el number of elements.

a_x is the time it takes to compute stats a a single set of data, and b_x is the time it takes to accumulate a single point. The old algorithm was developed in the early history of the project, I'm guessing by Neil Kileen, while the new algorithm was developed in 2015 by Dave Mehringer as part of the stats framework project. The old algorithm is faster in the regime of large n_sets and small n_el, while the new algorithm is faster in the regime of small n_sets and large n_el. If one always wants to use one of these algorithms, that algorithm's coefficients should be set to 0, while setting the other algorithm's coefficients to positive values. Note that it's the relative, not the absolute, values of these coeffecients that is important The version that takes no parameters uses the default values of the coefficients;

template<class T>

void casacore::LatticeStatistics< T >::configureFitToHalf	(	FitToHalfStatisticsData::CENTER	centerType = `FitToHalfStatisticsData::CMEAN`,
		FitToHalfStatisticsData::USE_DATA	useData = `FitToHalfStatisticsData::LE_CENTER`,
		AccumType	centerValue = `0`
	)

configure to use fit to half algorithm.

template<class T>

void casacore::LatticeStatistics< T >::configureHingesFences ( Double f )

configure to use hinges-fences algorithm

template<class T>

Bool casacore::LatticeStatistics< T >::display ( )

Display the statistics by listing and/or plotting them.

If you don't call this function then you won't see anything ! A return value of False indicates an invalid plotting device, or that the internal state of the class is bad.

template<class T>

Vector<Int> casacore::LatticeStatistics< T >::displayAxes ( ) const [inline]

Return the display axes.

The returned vector will be valid only if setAxes has been called, or if one of the active "display" or "get*" methods has been called.

Definition at line 302 of file LatticeStatistics.h.

References casacore::LatticeStatistics< T >::displayAxes_p.

template<class T>

virtual void casacore::LatticeStatistics< T >::displayStats	(	AccumType	nPts,
		AccumType	sum,
		AccumType	median,
		AccumType	medAbsDevMed,
		AccumType	quartile,
		AccumType	sumSq,
		AccumType	mean,
		AccumType	var,
		AccumType	rms,
		AccumType	sigma,
		AccumType	dMin,
		AccumType	dMax,
		AccumType	q1,
		AccumType	q3
	)			`[private, virtual]`

Reimplemented in casacore::ImageStatistics< T >.

template<class T>

String casacore::LatticeStatistics< T >::errorMessage ( ) const [inline]

Recover last error message.

Definition at line 351 of file LatticeStatistics.h.

References casacore::LatticeStatistics< T >::error_p.

template<class T>

void casacore::LatticeStatistics< T >::generateRobust ( ) [private]

Find the median per cursorAxes chunk.

template<class T>

Bool casacore::LatticeStatistics< T >::generateStorageLattice ( ) [private]

Create a new storage lattice.

template<class T>

std::map<String, uInt> casacore::LatticeStatistics< T >::getChauvenetNiter ( ) const [inline]

get number of iterations associated with Chauvenet criterion algorithm

Definition at line 400 of file LatticeStatistics.h.

References casacore::LatticeStatistics< T >::_chauvIters.

template<class T>

Bool casacore::LatticeStatistics< T >::getConvertedStatistic	(	Array< T > &	stat,
		LatticeStatsBase::StatisticsTypes	type,
		Bool	dropDeg = `True`
	)

template<class T>

Bool casacore::LatticeStatistics< T >::getFullMinMax	(	T &	dataMin,
		T &	dataMax
	)

Get full lattice min and max only.

Returns False if no unmasked data, else returns True. Honours any include or exclude range if set.

template<class T>

virtual void casacore::LatticeStatistics< T >::getLabels	(	String &	higherOrderLabel,
		String &	xAxisLabel,
		const IPosition &	dPos
	)			const `[protected, virtual]`

Reimplemented in casacore::ImageStatistics< T >.

template<class T>

Bool casacore::LatticeStatistics< T >::getLayerStats	(	stat_list &	stats,
		Double	area,
		Int	zAxis = `-1`,
		Int	zLayer = `-1`,
		Int	hAxis = `-1`,
		Int	hLayer = `-1`
	)

template<class T>

Bool casacore::LatticeStatistics< T >::getLayerStats	(	String &	stats,
		Double	area,
		Int	zAxis = `-1`,
		Int	zLayer = `-1`,
		Int	hAxis = `-1`,
		Int	hLayer = `-1`
	)

template<class T>

Bool casacore::LatticeStatistics< T >::getMinMaxPos	(	IPosition &	minPos,
		IPosition &	maxPos
	)

Recover position of min and max.

Only works if there are no display axes (i.e. statistics found over entire image), otherwise, the returned values are resized to 0 shape. A return value of False indicates that the internal state of the class is bad.

template<class T>

Bool casacore::LatticeStatistics< T >::getStatistic	(	Array< AccumType > &	stat,
		LatticeStatsBase::StatisticsTypes	type,
		Bool	dropDeg = `True`
	)

Recover the desired Statistic into an array.

If you choose to use the T version, be aware that the values in the AccumType version of the Array may not be representable in the T version (e.g. large values for SumSq). The shape of the array is the shape of the display axes (e.g. if the shape of the lattice is [nx,ny,nz] and you ask for the mean of the y axis the shape of the returned array would be [nx,nz]. A returned array of zero shape indicates that there were no good values. A return value of False indicates that the internal state of the class is bad.

template<class T>

Bool casacore::LatticeStatistics< T >::getStats	(	Vector< AccumType > &	,
		const IPosition &	pos,
		const Bool	posInLattice = `False`
	)

This function gets a vector containing all the statistics for a given location.

If posInLattice=True then the location is a location in the input lattice. Any positions on the display axes are ignored. Otherwise, you should just give locations for the display axes only. Use can use the enum in class LatticeStatsBase to find out which locations in the vector contain which statistics. A returned vector of zero shape indicates that there were no good values. A return value of False indicates that the internal state of the class is bad.

template<class T>

Bool casacore::LatticeStatistics< T >::hasLogger ( ) const [inline]

Did we construct with a logger ?

Definition at line 358 of file LatticeStatistics.h.

References casacore::LatticeStatistics< T >::haveLogger_p.

template<class T>

virtual Bool casacore::LatticeStatistics< T >::listLayerStats	(	const Matrix< AccumType > &	ord,
		ostringstream &	rslt,
		Int	zLayer
	)			`[protected, virtual]`

template<class T>

virtual void casacore::LatticeStatistics< T >::listMinMax	(	ostringstream &	osMin,
		ostringstream &	osMax,
		Int	oWidth,
		DataType	type
	)			`[protected, virtual]`

Reimplemented in casacore::ImageStatistics< T >.

template<class T>

virtual Bool casacore::LatticeStatistics< T >::listStats	(	Bool	hasBeam,
		const IPosition &	dPos,
		const Matrix< AccumType > &	ord
	)			`[protected, virtual]`

List the statistics to the logger.

The implementation here is adequate for all lattices. See ImageStatistics for an example of where extra information is provided. hasBeam is the return value of getBeamArea. If it is true, that means that the FLUX statistics will be available in the storage lattice. dPos is the location of the start of the cursor in the storage image for this row. stats(j,i) is the statistics matrix. for the jth point and the ith statistic. The return value is False if something goes wrong ! Have a look at the implementation to see what you really have to do.

Reimplemented in casacore::ImageStatistics< T >.

template<class T>

* IPosition casacore::LatticeStatistics< T >::locInLattice	(	const IPosition &	storagePosition,
		Bool	relativeToParent = `True`
	)			const `[protected]`

Given a location in the storage lattice, convert those locations on the non-statistics axis (the last one) and optionally account for the lattice subsectioning.

template<class T>

IPosition casacore::LatticeStatistics< T >::locInStorageLattice	(	const IPosition &	latticePosition,
		LatticeStatsBase::StatisticsTypes	type
	)			const `[private]`

Given a location in the lattice and a statistic type, work out where to put it in the storage lattice.

template<class T>

void casacore::LatticeStatistics< T >::minMax	(	Bool &	none,
		AccumType &	dMin,
		AccumType &	dMax,
		const Vector< AccumType > &	d,
		const Vector< AccumType > &	n
	)			const `[private]`

Find min and max of good data in arrays specified by pointers.

template<class T>

LatticeStatistics<T>& casacore::LatticeStatistics< T >::operator= ( const LatticeStatistics< T > & other )

Assignment operator.

Deletes any storage lattice associated with the object being assigned to and copies any storage lattice that has already been created for "other".

Reimplemented in casacore::ImageStatistics< T >.

template<class T>

void casacore::LatticeStatistics< T >::resetError ( ) [inline]

Reset argument error condition.

If you specify invalid arguments to one of the above set functions, an internal flag will be set which will prevent the work functions from doing anything (should you have chosen to ignore the Boolean return values of the set functions). This function allows you to reset that internal state to good.

Definition at line 344 of file LatticeStatistics.h.

References casacore::LatticeStatistics< T >::goodParameterStatus_p, and casacore::True.

template<class T>

Bool casacore::LatticeStatistics< T >::retrieveStorageStatistic	(	Vector< AccumType > &	slice,
		const IPosition &	pos,
		const Bool	posInLattice
	)			`[private]`

Retrieve a statistic from the storage lattice at the specified location and return in an array.

template<class T>

Bool casacore::LatticeStatistics< T >::retrieveStorageStatistic	(	Array< AccumType > &	slice,
		const LatticeStatsBase::StatisticsTypes	type,
		const Bool	dropDeg
	)			`[private]`

Retrieve a statistic from the storage lattice and return in an array.

template<class T>

Bool casacore::LatticeStatistics< T >::setAxes ( const Vector< Int > & cursorAxes )

Set the cursor axes (0 relative).

A return value of False indicates you have asked for an invalid axis. The default state of the class is to set the cursor axes to all axes in the lattice.

template<class T>

Bool casacore::LatticeStatistics< T >::setInExCludeRange	(	const Vector< T > &	include,
		const Vector< T > &	exclude,
		Bool	setMinMaxToInclude = `False`
	)

You may specify a pixel intensity range as either one for which all pixels in that range are included or one for which all pixels in that range are excluded.

One or the other of include and exclude must therefore be a zero length vector if you call this function. If you are setting an include range, then if you set setMinMaxToInclude=True, the minimum and maximum values that this class returns will always be the minimum and maximum of the include range, respectively. A return value of False indicates that you have given both an include and an exclude range. A vector of length 1 for include and/or exclude means that the range will be set to (say for include) -abs(include(0)) to abs(include(0)). A return value of False indicates that both an inclusion and exclusion range were given or that the internal state of the class is bad. If you don't call this function, the default state of the class is to include all pixels.

template<class T>

Bool casacore::LatticeStatistics< T >::setList ( const Bool & doList )

This function allows you to control whether the statistics are written to the output stream if you are also making a plot.

A return value of False indicates that the internal state of the class is bad. If you have created the LatticeStatistics object without a LogIO object, you won't see any listings, but no error conditions will be generated. The default state of the class is to not list the output when making a plot.

template<class T>

Bool casacore::LatticeStatistics< T >::setNewLattice ( const MaskedLattice< T > & lattice )