casa  $Rev:20696$
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines
Public Member Functions | Protected Attributes
casa::MaskedArray< T > Class Template Reference

Class for masking an Array for operations on that Array. More...

#include <MaskedArray.h>

List of all members.

Public Member Functions

 MaskedArray ()
 Default constructor for a MaskedArray does not allocate any memory for the Data array or Mask.
void setData (const Array< T > &data, const LogicalArray &mask, Bool isReadOnly=False)
 Reset the data and mask of the the MaskedArray.
void setData (const MaskedArray< T > &array, Bool isReadOnly=False)
 MaskedArray (const Array< T > &inarray, const LogicalArray &inmask, Bool isreadonly)
 Create a MaskedArray from an Array and a LogicalArray.
 MaskedArray (const Array< T > &inarray, const LogicalArray &inmask)
 MaskedArray (const MaskedArray< T > &inarray, const LogicalArray &inmask, Bool isreadonly)
 Create a MaskedArray from a MaskedArray and a LogicalArray.
 MaskedArray (const MaskedArray< T > &inarray, const LogicalArray &inmask)
 MaskedArray (const Array< T > &inarray, const MaskedLogicalArray &inmask, Bool isreadonly)
 Create a MaskedArray from an Array and a MaskedLogicalArray.
 MaskedArray (const Array< T > &inarray, const MaskedLogicalArray &inmask)
 MaskedArray (const MaskedArray< T > &inarray, const MaskedLogicalArray &inmask, Bool isreadonly)
 Create a MaskedArray from a MaskedArray and a MaskedLogicalArray.
 MaskedArray (const MaskedArray< T > &inarray, const MaskedLogicalArray &inmask)
 MaskedArray (const MaskedArray< T > &other, Bool isreadonly)
 Copy constructor.
 MaskedArray (const MaskedArray< T > &other)
 ~MaskedArray ()
 
     

MaskedArray< T > operator() (const LogicalArray &mask) const
 Return a MaskedArray.
MaskedArray< T > operator() (const MaskedLogicalArray &mask) const
 Return a MaskedArray.
MaskedArray< T > operator() (const IPosition &start, const IPosition &end)
 Get a reference to an array part which extends from "start" to end.
MaskedArray< T > operator() (const IPosition &start, const IPosition &end, const IPosition &inc)
 Along the ith axis, every inc[i]'th element is chosen.
MaskedArray< T > operator() (const Slicer &)
 Get a reference to an array using a Slicer.
MaskedArray< T > copy (Bool isreadonly) const
 Make a copy of the masked array.
MaskedArray< T > copy () const
const Array< T > & getArray () const
 Return the internal Array.
Array< T > & getRWArray () const
 Return the internal Array, writeable.
const LogicalArray & getMask () const
 Return the (const) internal Mask.
uInt ndim () const
 The dimensionality of this masked array.
uInt nelements () const
 The number of elements of this masked array.
uInt size () const
uInt nelementsValid () const
 The number of valid elements of this masked array.
Bool ok () const
 Check to see if the masked array is consistent.
Bool conform (const Array< T > &other) const
 Are the shapes identical?
Bool conform (const MaskedArray< T > &other) const
const IPositionshape () const
 The length of each axis.
Bool isReadOnly () const
 Is the array read only?
void setReadOnly () const
 Set the array to be read only.
MaskedArray< T > & operator= (const Array< T > &inarray)
 Copy the values in inarray to this, only copying those elements for which the corresponding mask element is True.
MaskedArray< T > & operator= (const MaskedArray< T > &other)
 Copy the values in other to this, only copying those elements for which the logical AND of the corresponding mask elements of both MaskedArrays is True.
MaskedArray< T > & operator= (const T &value)
 Set every element of this array to "value", only setting those elements for which the corresponding mask element is True.
Array< T > getCompressedArray () const
 Return a "compressed" Array containing only the valid elements of the MaskedArray.
Array< T > getCompressedArray (const IPosition &shape) const
 The returned Array will have the input shape.
void getCompressedArray (Array< T > &inarr) const
 Fill the argument "compressed" Array with only the valid elements of the MaskedArray.
void setCompressedArray (const Array< T > &inarr)
 Set only the valid elements of the MaskedArray from the argument "compressed" Array.
const T * getArrayStorage (Bool &deleteIt) const
 Manipulate the storage for the underlying Array.
T * getRWArrayStorage (Bool &deleteIt) const
 
   

void freeArrayStorage (const T *&storage, Bool deleteIt) const
void putArrayStorage (T *&storage, Bool deleteAndCopy) const
 
   

const LogicalArrayElem * getMaskStorage (Bool &deleteIt) const
 Manipulate the storage for the underlying Mask.
void freeMaskStorage (const LogicalArrayElem *&storage, Bool deleteIt) const

Protected Attributes

Array< T > * pArray
 
      

LogicalArray * pMask
 The mask.
uInt nelemValid
 Cache the number of valid elements.
Bool nelemValidIsOK
 Is the number of valid elements cache OK? i.e.
Bool isRO
 Is the array read only?

Detailed Description

template<class T>
class casa::MaskedArray< T >

Class for masking an Array for operations on that Array.

Review Status

Reviewed By:
UNKNOWN
Date Reviewed:
before2004/08/25
Test programs:
tMaskedArray tMaskArrExcp

Prerequisite

Etymology

MaskedArray is a class for masking elements of an Array while performing operations on that Array.

Synopsis

A MaskedArray is an association between an Array and a mask. The mask selects elements of the Array. Only elements of the Array where the corresponding element of the mask is True are defined. Thus, operations on a MaskedArray only operate on those elements of the Array where the corresponding element of the mask is True.

A MaskedArray should be thought of as a manipulator for an Array, analogous to an iterator. It allows one to perform whole Array operations on selected elements of the Array.

The mask used in the constructor for the MaskedArray must conform to the Array, thus have the same shape. The internal mask is (will be) copy constructed with reference semantics from the input mask. Therefore, it is (will be) possible to change the internal mask by changing values in the input mask after the MaskedArray has been constructed. To ensure that the internal mask is independent of the input mask after construction, use mask.copy() as the input argument.

One can explicitly construct a MaskedArray from an Array and a mask or a MaskedArray and a mask. One can also use operator() on an Array or a MaskedArray to implicitly create a MaskedArray.

One can create a MaskedArray from a MaskedArray and a mask. The resulting MaskedArray has as its Array the Array from the original MaskedArray. The mask for the resulting MaskedArray is the AND of the mask from the original MaskedArray and the input mask.

Any operation involving a MaskedArray or a set of MaskedArrays is only performed for those elements where the AND of the masks is True.

Any operation involving a MaskedArray or a set of MaskedArrays results in a MaskedArray whose mask is the AND of the masks of the original MaskedArrays. The only exception to this is assignment, where the mask determines which elements of the underlying Array are assigned.

Masks, which are LogicalArrays, can be constructed by logical operations involving Arrays. They can also, of course, be constructed by individually setting individual elements of an LogicalArray.

MaskedArrays constructed directly from Arrays are by default writeable. MaskedArrays constructed indirectly from Arrays by operator() are writeable if the Array is non-const and are readonly if the Array is const. MaskedArrays constructed from other MaskedArrays, either directly by constructors or indirectly by operator(), are by default writeable if the input MaskedArray is writeable, and readonly if the input MaskedArray is readonly.

A given MaskedArray can be set to be readonly. One specifies this in the constructor with the Bool argument isreadonly, or calls the setReadOnly() member function. A MaskedArray which would default to be readonly cannot be forced to be writeable. It will remain readonly even if the Bool argument isreadonly is set to be False.

The isReadOnly(), member function is used to test whether the MaskedArray is readonly.

Member functions which change the MaskedArray test to see whether the MaskedArray is readonly, and throw an ArrayError exception if it is. These member functions are:

The copy() member function makes a deep copy of a MaskedArray. By default it returns a writeable MaskedArray, but the MaskedArray returned can be made readonly by using the Bool argument "isreadonly" to copy() (or by calling setReadOnly() on the new MaskedArray).

The valid elements of the MaskedArray can be manipulated as a "compressed" Array which contains only the valid elements. The number of elements in this "compressed" Array is the number of valid elements in the MaskedArray, nelementsValid(). The "compressed" Array can have any shape which meets this requirement. The MaskedArray can have any shape.

The getCompressedArray() member functions get a compressed Array from the valid members of the MaskedArray, while the setCompressedArray() member function sets the valid members of the MaskedArray from the input compressed Array.

Many mathematical and logical global operators and functions which operate on MaskedArrays are defined. Typically, they are defined for all sensible combinations of MaskedArrays, Arrays, and scalars.

Mathematical global operators and functions are defined in Arrays/MaskArrMath.h . The list is:

Logical global operators and functions are defined in Arrays/MaskArrLogi.h . The list is:

Example

Use an explicit MaskedArray to limit the maximum value of an Array.

      Vector<Int> arr (20);
         . . .
      MaskedArray<Int> marr (arr, (arr > 5));
      marr = 5;

This sets all elements of arr which are greater than 5 to 5.

Example

Use an implicit MaskedArray to limit the minimum value of an Array.

      Vector<Int> arr (20);
         . . .
      arr (arr < 0) = 0;

This sets all elements of arr which are less than 0 to 0.

Example

It does not matter where in an expression the MaskedArrays are located. The operation is only performed on those elements where the AND of the masks is True.

The following expressions are all equivalent. The first (and second) expressions are the most efficient, since the sum is only performed for those elements where ((a > 0) && (b > 0)). The third example is less efficient, since the sum is performed for all elements of a and b, and then the assignment is only performed for those elements where ((a > 0) && (b > 0)).

      Vector<Int> arr (20);
      Vector<Int> a (20);
      Vector<Int> b (20);
         . . .
      arr = a(a > 0) + b(b > 0);
   
      arr = a(b > 0) + b(a > 0);
   
      arr ((a > 0) && (b > 0)) = a + b;
   
      arr = (a + b) ((a > 0) && (b > 0));
   
      arr (a > 0) = a + b(b > 0);

All of these expressions set those elements of arr where ((a > 0) && (b > 0)) to a + b. Those elements of arr where the condition is False are unchanged.

Example

This example extracts the valid elements of the MaskedArray as a "compressed" Vector, manipulates this Vector, and then puts the result back into the MaskedArray.

      Matrix<Int> arr (20,5);
         . . .
      MaskedArray<Int> marr (arr, (arr>0) && (arr<10));
      Vector<Int> vec (marr.getCompressedArray());
         . . .
      marr.setCompressedArray (vec);

Motivation

A MaskedArray is an association between an Array and a LogicalArray which masks the Array. It allows one to perform whole Array manipulations with a single expression, selecting those elements of the Array to modify based either on a logical expression, typically involving some of the Arrays involved in the expression, or based on a specifically set mask.

To Do

Definition at line 325 of file MaskedArray.h.


Constructor & Destructor Documentation

template<class T>
casa::MaskedArray< T >::MaskedArray ( )

Default constructor for a MaskedArray does not allocate any memory for the Data array or Mask.

Hence the masked array should not be used until some data is allocated to the object using one of the set functions.

template<class T>
casa::MaskedArray< T >::MaskedArray ( const Array< T > &  inarray,
const LogicalArray &  inmask,
Bool  isreadonly 
)

Create a MaskedArray from an Array and a LogicalArray.

The internal mask is a total copy of the input mask, and is completely independent of the input mask.

The Array is copy constructed, which means that it is a really smart pointer to the underlying Block, and shares this Block with the input Array.

By default, the MaskedArray constructed is writeable. If isreadonly is True, then the MaskedArray returned is readonly.

Thrown Exceptions

template<class T>
casa::MaskedArray< T >::MaskedArray ( const Array< T > &  inarray,
const LogicalArray &  inmask 
)
template<class T>
casa::MaskedArray< T >::MaskedArray ( const MaskedArray< T > &  inarray,
const LogicalArray &  inmask,
Bool  isreadonly 
)

Create a MaskedArray from a MaskedArray and a LogicalArray.

The internal mask is the AND of the input mask and the mask of the input MaskedArray.

The Array from the input MaskedArray is copy constructed, which means that it is a really smart pointer to the underlying Block, and shares this Block with the Array from the input MaskedArray.

By default, the MaskedArray constructed is writeable if the input MaskedArray is writeable, and readonly if the input MaskedArray is readonly. If isreadonly is True, then the MaskedArray returned is readonly. If isreadonly is False and the input MaskedArray is readonly, then the constructed MaskedArray is readonly.

Thrown Exceptions

template<class T>
casa::MaskedArray< T >::MaskedArray ( const MaskedArray< T > &  inarray,
const LogicalArray &  inmask 
)
template<class T>
casa::MaskedArray< T >::MaskedArray ( const Array< T > &  inarray,
const MaskedLogicalArray &  inmask,
Bool  isreadonly 
)

Create a MaskedArray from an Array and a MaskedLogicalArray.

The internal mask is the AND of the internal LogicalArray and the internal mask of the MaskedLogicalArray.

The Array is copy constructed, which means that it is a really smart pointer to the underlying Block, and shares this Block with the input Array.

By default, the MaskedArray constructed is writeable. If isreadonly is True, then the MaskedArray returned is readonly.

Thrown Exceptions

template<class T>
casa::MaskedArray< T >::MaskedArray ( const Array< T > &  inarray,
const MaskedLogicalArray &  inmask 
)
template<class T>
casa::MaskedArray< T >::MaskedArray ( const MaskedArray< T > &  inarray,
const MaskedLogicalArray &  inmask,
Bool  isreadonly 
)

Create a MaskedArray from a MaskedArray and a MaskedLogicalArray.

The internal mask is the AND of the internal LogicalArray and the internal mask of the MaskedLogicalArray, ANDed with the mask of the input MaskedArray.

The Array from the input MaskedArray is copy constructed, which means that it is a really smart pointer to the underlying Block, and shares this Block with the Array from the input MaskedArray.

By default, the MaskedArray constructed is writeable if the input MaskedArray is writeable, and readonly if the input MaskedArray is readonly. If isreadonly is True, then the MaskedArray returned is readonly. If isreadonly is False and the input MaskedArray is readonly, then the constructed MaskedArray is readonly.

Thrown Exceptions

template<class T>
casa::MaskedArray< T >::MaskedArray ( const MaskedArray< T > &  inarray,
const MaskedLogicalArray &  inmask 
)
template<class T>
casa::MaskedArray< T >::MaskedArray ( const MaskedArray< T > &  other,
Bool  isreadonly 
)

Copy constructor.

The internal mask is a total copy of the mask from the input MaskedArray, and is completely independent of this input mask.

The Array from the input MaskedArray is copy constructed, which means that it is a really smart pointer to the underlying Block, and shares this Block with the Array from the input MaskedArray.

By default, the MaskedArray constructed is writeable if the input MaskedArray is writeable, and readonly if the input MaskedArray is readonly. If isreadonly is True, then the MaskedArray returned is readonly. If isreadonly is False and the input MaskedArray is readonly, then the constructed MaskedArray is readonly.

template<class T>
casa::MaskedArray< T >::MaskedArray ( const MaskedArray< T > &  other)
template<class T>
casa::MaskedArray< T >::~MaskedArray ( )

     


Member Function Documentation

template<class T>
Bool casa::MaskedArray< T >::conform ( const Array< T > &  other) const

Are the shapes identical?

template<class T>
Bool casa::MaskedArray< T >::conform ( const MaskedArray< T > &  other) const
template<class T>
MaskedArray<T> casa::MaskedArray< T >::copy ( Bool  isreadonly) const

Make a copy of the masked array.

This is a deep copy. The Array and mask components of the returned MaskedArray are deep copies of the Array and mask in the input MaskedArray pointed to by this. In other words, the Array and mask in the output MaskedArray are completely independent of those in the input MaskedArray.

By default, the MaskedArray returned is writeable. If isreadonly is True, then the MaskedArray returned is readonly.

template<class T>
MaskedArray<T> casa::MaskedArray< T >::copy ( ) const
template<class T>
void casa::MaskedArray< T >::freeArrayStorage ( const T *&  storage,
Bool  deleteIt 
) const
template<class T>
void casa::MaskedArray< T >::freeMaskStorage ( const LogicalArrayElem *&  storage,
Bool  deleteIt 
) const
template<class T>
const Array<T>& casa::MaskedArray< T >::getArray ( ) const

Return the internal Array.

template<class T>
const T* casa::MaskedArray< T >::getArrayStorage ( Bool deleteIt) const

Manipulate the storage for the underlying Array.

See the description of the corresponding Array functions for more information.

template<class T>
Array<T> casa::MaskedArray< T >::getCompressedArray ( ) const

Return a "compressed" Array containing only the valid elements of the MaskedArray.

The number of elements in the Array will be nelementsValid() for the MaskedArray. The MaskedArray can have any shape.

The returned Array will have dimension one.

template<class T>
Array<T> casa::MaskedArray< T >::getCompressedArray ( const IPosition shape) const

The returned Array will have the input shape.

This shape must give the returned Array the required number of elements.

Thrown Exceptions

template<class T>
void casa::MaskedArray< T >::getCompressedArray ( Array< T > &  inarr) const

Fill the argument "compressed" Array with only the valid elements of the MaskedArray.

The size of the Array must be nelementsValid() for the MaskedArray. The Array can have any shape which meets this requirement. The MaskedArray can have any shape.

Thrown Exceptions

template<class T>
const LogicalArray& casa::MaskedArray< T >::getMask ( ) const

Return the (const) internal Mask.

template<class T>
const LogicalArrayElem* casa::MaskedArray< T >::getMaskStorage ( Bool deleteIt) const

Manipulate the storage for the underlying Mask.

See the description of the corresponding Array functions for more information.

template<class T>
Array<T>& casa::MaskedArray< T >::getRWArray ( ) const

Return the internal Array, writeable.

Thrown Exceptions

template<class T>
T* casa::MaskedArray< T >::getRWArrayStorage ( Bool deleteIt) const

   

Thrown Exceptions

template<class T>
Bool casa::MaskedArray< T >::isReadOnly ( ) const [inline]

Is the array read only?

Definition at line 559 of file MaskedArray.h.

References casa::MaskedArray< T >::isRO.

template<class T>
uInt casa::MaskedArray< T >::ndim ( ) const

The dimensionality of this masked array.

template<class T>
uInt casa::MaskedArray< T >::nelements ( ) const

The number of elements of this masked array.

This is the number of elements in the underlying Array.

Referenced by casa::MaskArrMath_global_functions_MaskedArray_mathematical_operations::median(), and casa::MaskedArray< T >::size().

template<class T>
uInt casa::MaskedArray< T >::nelementsValid ( ) const

The number of valid elements of this masked array.

This is the number of elements of the mask which are TRUE.

template<class T>
Bool casa::MaskedArray< T >::ok ( ) const

Check to see if the masked array is consistent.

This is about the same thing as checking for invariants. If AIPS_DEBUG is defined, this is invoked after construction and on entry to most member functions.

template<class T>
MaskedArray<T> casa::MaskedArray< T >::operator() ( const LogicalArray &  mask) const

Return a MaskedArray.

The new MaskedArray is masked by the input LogicalArray "anded" with the mask of the original MaskedArray. This mask must conform to the array.

The MaskedArray constructed is writeable if the input MaskedArray is writeable, and readonly if the input MaskedArray is readonly.

template<class T>
MaskedArray<T> casa::MaskedArray< T >::operator() ( const MaskedLogicalArray &  mask) const

Return a MaskedArray.

The new MaskedArray is masked by the input MaskedLogicalArray "anded" with the mask of the original MaskedArray. This mask must conform to the array.

The MaskedArray constructed is writeable if the input MaskedArray is writeable, and readonly if the input MaskedArray is readonly.

template<class T>
MaskedArray<T> casa::MaskedArray< T >::operator() ( const IPosition start,
const IPosition end 
)

Get a reference to an array part which extends from "start" to end.

"

template<class T>
MaskedArray<T> casa::MaskedArray< T >::operator() ( const IPosition start,
const IPosition end,
const IPosition inc 
)

Along the ith axis, every inc[i]'th element is chosen.

template<class T>
MaskedArray<T> casa::MaskedArray< T >::operator() ( const Slicer )

Get a reference to an array using a Slicer.

template<class T>
MaskedArray<T>& casa::MaskedArray< T >::operator= ( const Array< T > &  inarray)

Copy the values in inarray to this, only copying those elements for which the corresponding mask element is True.

Thrown Exceptions

template<class T>
MaskedArray<T>& casa::MaskedArray< T >::operator= ( const MaskedArray< T > &  other)

Copy the values in other to this, only copying those elements for which the logical AND of the corresponding mask elements of both MaskedArrays is True.

Thrown Exceptions

template<class T>
MaskedArray<T>& casa::MaskedArray< T >::operator= ( const T &  value)

Set every element of this array to "value", only setting those elements for which the corresponding mask element is True.

In other words, a scalar behaves as if it were a constant conformant array.

Thrown Exceptions

template<class T>
void casa::MaskedArray< T >::putArrayStorage ( T *&  storage,
Bool  deleteAndCopy 
) const

   

Thrown Exceptions

template<class T>
void casa::MaskedArray< T >::setCompressedArray ( const Array< T > &  inarr)

Set only the valid elements of the MaskedArray from the argument "compressed" Array.

The size of the Array must be nelementsValid() for the MaskedArray. The Array can have any shape which meets this requirement. The MaskedArray can have any shape.

Thrown Exceptions

template<class T>
void casa::MaskedArray< T >::setData ( const Array< T > &  data,
const LogicalArray &  mask,
Bool  isReadOnly = False 
)

Reset the data and mask of the the MaskedArray.

There should perhaps be a whole family of setData functions with different arguements, analogous to the constructors. However these are sufficient for the moment.

template<class T>
void casa::MaskedArray< T >::setData ( const MaskedArray< T > &  array,
Bool  isReadOnly = False 
)
template<class T>
void casa::MaskedArray< T >::setReadOnly ( ) const

Set the array to be read only.

template<class T>
const IPosition& casa::MaskedArray< T >::shape ( ) const [inline]

The length of each axis.

Definition at line 555 of file MaskedArray.h.

References casa::MaskedArray< T >::pArray.

template<class T>
uInt casa::MaskedArray< T >::size ( ) const [inline]

Definition at line 534 of file MaskedArray.h.

References casa::MaskedArray< T >::nelements().


Member Data Documentation

template<class T>
Bool casa::MaskedArray< T >::isRO [protected]

Is the array read only?

Definition at line 691 of file MaskedArray.h.

Referenced by casa::MaskedArray< T >::isReadOnly().

template<class T>
uInt casa::MaskedArray< T >::nelemValid [protected]

Cache the number of valid elements.

Definition at line 684 of file MaskedArray.h.

template<class T>
Bool casa::MaskedArray< T >::nelemValidIsOK [protected]

Is the number of valid elements cache OK? i.e.

has it been calculated?

Definition at line 688 of file MaskedArray.h.

template<class T>
Array<T>* casa::MaskedArray< T >::pArray [protected]

      

The array.

Definition at line 678 of file MaskedArray.h.

Referenced by casa::MaskedArray< T >::shape().

template<class T>
LogicalArray* casa::MaskedArray< T >::pMask [protected]

The mask.

Definition at line 681 of file MaskedArray.h.


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