MaskedArray.h

Classes

MaskedArray -- Class for masking an Array for operations on that Array. (full description)
Global Functions -- General global functions for MaskedArrays, and MaskedArrays and Arrays. (full description)

template<class T> class MaskedArray

Interface

Public Members
MaskedArray()
void setData(const Array<T> & data, const LogicalArray & mask, Bool isReadOnly=False)
void setData(const MaskedArray<T> & array, Bool isReadOnly=False)
MaskedArray(const Array<T> &inarray, const LogicalArray &inmask, Bool isreadonly)
MaskedArray(const Array<T> &inarray, const LogicalArray &inmask)
MaskedArray(const MaskedArray<T> &inarray, const LogicalArray &inmask, Bool isreadonly)
MaskedArray(const MaskedArray<T> &inarray, const LogicalArray &inmask)
MaskedArray(const Array<T> &inarray, const MaskedLogicalArray &inmask, Bool isreadonly)
MaskedArray(const Array<T> &inarray, const MaskedLogicalArray &inmask)
MaskedArray(const MaskedArray<T> &inarray, const MaskedLogicalArray &inmask, Bool isreadonly)
MaskedArray(const MaskedArray<T> &inarray, const MaskedLogicalArray &inmask)
MaskedArray(const MaskedArray<T> &other, Bool isreadonly)
MaskedArray(const MaskedArray<T> &other)
~MaskedArray()
MaskedArray<T> operator() (const LogicalArray &mask) const
MaskedArray<T> operator() (const MaskedLogicalArray &mask) const
MaskedArray<T> operator()(const IPosition &start, const IPosition &end)
MaskedArray<T> operator()(const IPosition &start, const IPosition &end, const IPosition &inc)
MaskedArray<T> operator()(const Slicer&)
MaskedArray<T> copy(Bool isreadonly) const
MaskedArray<T> copy() const
const Array<T> & getArray() const
Array<T> & getRWArray() const
const LogicalArray & getMask() const
uInt ndim() const
uInt nelements() const
uInt nelementsValid() const
Bool ok() const
Bool conform(const Array<T> &other) const
Bool conform(const MaskedArray<T> &other) const
IPosition shape() const
Bool isReadOnly() const
void setReadOnly() const
MaskedArray<T> &operator=(const Array<T> &inarray)
MaskedArray<T> &operator=(const MaskedArray<T> &other)
MaskedArray<T> &operator=(const T &value)
Array<T> getCompressedArray () const
Array<T> getCompressedArray (const IPosition & shape) const
void getCompressedArray (Array<T> & inarr) const
void setCompressedArray (const Array<T> & inarr)
const T * getArrayStorage (Bool &deleteIt) const
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
void freeMaskStorage(const LogicalArrayElem *&storage, Bool deleteIt) const
See Also
LogicalArray forwards -- Forward declarations for LogicalArrays.
LogicalArray -- Logical valued Arrays.
LogicalCube -- Logical valued Cubes.
LogicalMatrix -- Logical valued Matrices.
LogicalVector -- Logical valued Vectors.
MaskedArray IO -- Ascii input/output operations for MaskedArrays.
MaskedArray logical operations -- Logical operations for MaskedArrays, and between MaskedArrays and Arrays.
MaskedArray mathematical operations -- Mathematical operations for MaskedArrays, and between MaskedArrays and Arrays.
MaskedArray general global functions -- General global functions for MaskedArrays, and between MaskedArrays and Arrays.
MaskedLogicalArray forwards -- Forward declarations for MaskedLogicalArrays.
MaskedLogicalArray -- Masked LogicalArrays.

Description

Review Status

Reviewed By:
UNKNOWN
Date Reviewed:
before2004/08/25
Programs:
Tests:

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

Member Description

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.

void 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.

void setData(const MaskedArray<T> & array, Bool isReadOnly=False)

MaskedArray(const Array<T> &inarray, const LogicalArray &inmask, Bool isreadonly)
MaskedArray(const Array<T> &inarray, const LogicalArray &inmask)

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

MaskedArray(const MaskedArray<T> &inarray, const LogicalArray &inmask, Bool isreadonly)
MaskedArray(const MaskedArray<T> &inarray, const LogicalArray &inmask)

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

MaskedArray(const Array<T> &inarray, const MaskedLogicalArray &inmask, Bool isreadonly)
MaskedArray(const Array<T> &inarray, const MaskedLogicalArray &inmask)

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

MaskedArray(const MaskedArray<T> &inarray, const MaskedLogicalArray &inmask, Bool isreadonly)
MaskedArray(const MaskedArray<T> &inarray, const MaskedLogicalArray &inmask)

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

MaskedArray(const MaskedArray<T> &other, Bool isreadonly)
MaskedArray(const MaskedArray<T> &other)

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.

Thrown Exceptions

~MaskedArray()

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.

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.

MaskedArray<T> operator()(const IPosition &start, const IPosition &end, const IPosition &inc)

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

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

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 Slicer&)

Get a reference to an array using a Slicer.

MaskedArray<T> copy(Bool isreadonly) const
MaskedArray<T> copy() 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.

const Array<T> & getArray() const

Return the internal Array.

Array<T> & getRWArray() const

Return the internal Array, writeable.

Thrown Exceptions

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. This is the number of elements in the underlying Array.

uInt nelementsValid() const

The number of valid elements of this masked array. This is the number of elements of the mask which are TRUE.

Bool 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.

Bool conform(const Array<T> &other) const
Bool conform(const MaskedArray<T> &other) const

Are the shapes identical?

IPosition shape() 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.

Thrown Exceptions

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

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

Array<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.

Array<T> getCompressedArray (const IPosition & shape) 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 the input shape. This shape must give the returned Array the required number of elements.

Thrown Exceptions

void 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

void 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

T * getRWArrayStorage (Bool &deleteIt) const

Manipulate the storage for the underlying Array. See the description of the corresponding Array functions for more information.

Thrown Exceptions

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

Manipulate the storage for the underlying Array. See the description of the corresponding Array functions for more information.

Thrown Exceptions

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

Manipulate the storage for the underlying Array. See the description of the corresponding Array functions for more information.

const LogicalArrayElem *getMaskStorage (Bool &deleteIt) const
void freeMaskStorage(const LogicalArrayElem *&storage, Bool deleteIt) const

Manipulate the storage for the underlying Mask. See the description of the corresponding Array functions for more information.


General global functions for MaskedArrays, and MaskedArrays and Arrays. (source)

Interface

Protected Members
Bool conform2 (const MaskedArray<T> &left, const Array<U> &right)
Bool conform2 (const Array<T> &left, const MaskedArray<U> &right)
Bool conform2 (const MaskedArray<T> &left, const MaskedArray<U> &right)

Description

Review Status

Reviewed By:
UNKNOWN
Date Reviewed:
before2004/08/25
Programs:
Tests:
  • tMaskedArray

Prerequisite

Synopsis

These are generally useful global functions which operate on all MaskedArrays, or on MaskedArrays and Arrays.

Member Description

Bool conform2 (const MaskedArray<T> &left, const Array<U> &right)
Bool conform2 (const Array<T> &left, const MaskedArray<U> &right)
Bool conform2 (const MaskedArray<T> &left, const MaskedArray<U> &right)

Test conformance for masked arrays and arrays of different types. Are the shapes identical?