casacore
|
Interface for converting between world and pixel coordinates. More...
#include <Coordinate.h>
Public Types | |
enum | Type { LINEAR, DIRECTION, SPECTRAL, STOKES, TABULAR, QUALITY, COORDSYS } |
This enum lists the types of the derived classes. More... | |
enum | formatType { DEFAULT, SCIENTIFIC, FIXED, MIXED, TIME } |
This enum is used for formatting world values into Strings. More... | |
Public Member Functions | |
virtual | ~Coordinate () |
Destructor. More... | |
virtual Type | type () const =0 |
List the type of this Coordinate object. More... | |
virtual String | showType () const =0 |
virtual uInt | nPixelAxes () const =0 |
How many world/pixel axes are there in this Coordinate? While the number of world and pixel axes will generally be the same, it is not a requirement. More... | |
virtual uInt | nWorldAxes () const =0 |
virtual Bool | toWorld (Vector< Double > &world, const Vector< Double > &pixel, Bool useConversionFrame=True) const =0 |
Convert an absolute pixel position to an absolute world position or vice versa. More... | |
virtual Bool | toPixel (Vector< Double > &pixel, const Vector< Double > &world) const =0 |
virtual Bool | toMix (Vector< Double > &worldOut, Vector< Double > &pixelOut, const Vector< Double > &worldIn, const Vector< Double > &pixelIn, const Vector< Bool > &worldAxes, const Vector< Bool > &pixelAxes, const Vector< Double > &worldMin, const Vector< Double > &worldMax) const |
Mixed absolute pixel/world coordinate conversion. More... | |
virtual Bool | setWorldMixRanges (const IPosition &shape) |
Set the world min and max ranges, for use in function toMix , for a lattice of the given shape for this coordinate. More... | |
virtual void | setDefaultWorldMixRanges () |
Vector< Double > | worldMixMin () const |
Vector< Double > | worldMixMax () const |
virtual Bool | toWorldMany (Matrix< Double > &world, const Matrix< Double > &pixel, Vector< Bool > &failures) const |
Batch up a lot of transformations. More... | |
virtual Bool | toPixelMany (Matrix< Double > &pixel, const Matrix< Double > &world, Vector< Bool > &failures) const |
virtual void | makePixelRelative (Vector< Double > &pixel) const |
Make absolute coordinates relative and vice-versa (with respect to the reference value). More... | |
virtual void | makePixelAbsolute (Vector< Double > &pixel) const |
virtual void | makeWorldRelative (Vector< Double > &world) const |
virtual void | makeWorldAbsolute (Vector< Double > &world) const |
virtual void | makeWorldAbsoluteRef (Vector< Double > &world, const Vector< Double > &refVal) const |
Make absolute coordinates relative and vice versa with respect to the given reference value. More... | |
virtual void | makePixelRelativeMany (Matrix< Double > &pixel) const |
Batch up a lot of absolute/relative transformations. More... | |
virtual void | makePixelAbsoluteMany (Matrix< Double > &pixel) const |
virtual void | makeWorldRelativeMany (Matrix< Double > &world) const |
virtual void | makeWorldAbsoluteMany (Matrix< Double > &world) const |
virtual Vector< String > | worldAxisNames () const =0 |
Return the requested attributed. More... | |
virtual Vector< Double > | referencePixel () const =0 |
virtual Matrix< Double > | linearTransform () const =0 |
virtual Vector< Double > | increment () const =0 |
virtual Vector< Double > | referenceValue () const =0 |
virtual Vector< String > | worldAxisUnits () const =0 |
virtual Bool | setWorldAxisNames (const Vector< String > &names)=0 |
Set the requested attribute. More... | |
virtual Bool | setReferencePixel (const Vector< Double > &refPix)=0 |
virtual Bool | setLinearTransform (const Matrix< Double > &xform)=0 |
virtual Bool | setIncrement (const Vector< Double > &inc)=0 |
virtual Bool | setReferenceValue (const Vector< Double > &refval)=0 |
virtual Bool | setWorldAxisUnits (const Vector< String > &units)=0 |
Change the units. More... | |
virtual Coordinate * | makeFourierCoordinate (const Vector< Bool > &axes, const Vector< Int > &shape) const |
Find the Coordinate for when we Fourier Transform ourselves. More... | |
const String & | errorMessage () const |
If the last conversion to world or pixel coordinates resulted in an error, report that error. More... | |
virtual Bool | near (const Coordinate &other, Double tol=1.0e-6) const =0 |
Comparison to fractional tolerance (for floating point values). More... | |
virtual Bool | near (const Coordinate &other, const Vector< Int > &excludeAxes, Double tol=1.0e-6) const =0 |
virtual void | getPrecision (Int &precision, Coordinate::formatType &format, Bool showAsAbsolute, Int defPrecScientific, Int defPrecFixed, Int defPrecTime) const |
Provide a common interface to getting formatted representations of coordinate values. More... | |
virtual String | format (String &units, Coordinate::formatType format, Double worldValue, uInt axis, Bool isAbsolute=True, Bool showAsAbsolute=True, Int precision=-1, Bool usePrecForMixed=False) const |
String | formatQuantity (String &units, Coordinate::formatType format, const Quantum< Double > &worldValue, uInt axis, Bool isAbsolute=True, Bool showAsAbsolute=True, Int precision=-1) |
virtual Bool | save (RecordInterface &container, const String &fieldName) const =0 |
Used for persistence. More... | |
virtual Coordinate * | clone () const =0 |
Make a copy of ourself. More... | |
virtual Bool | doNearPixel (const Coordinate &other, const Vector< Bool > &thisAxes, const Vector< Bool > &otherAxes, Double tol=1.0e-6) const |
Comparison only made for specified axes in this and other Coordinate The default implementation should be ok for all Coordinate types except Stokes and Quality... More... | |
virtual Coordinate * | rotate (const Quantum< Double > &angle) const |
return the result of rotating the coordinate clockwise through the specified angle. More... | |
Static Public Member Functions | |
static String | typeToString (Coordinate::Type type) |
static void | set_wcs (::wcsprm &wcs) |
Call wcsset on the wcs structure. More... | |
static void | init_wcs (::wcsprm &wcs, int naxis) |
Call wcsini on the wcs structure. More... | |
static void | sub_wcs (const ::wcsprm &src, int &nsub, int axes[],::wcsprm &dst) |
Call wcssub on the src/dst pair. More... | |
static void | copy_wcs (const ::wcsprm &src,::wcsprm &dst) |
Call wcssub on the src/dst pair with null nsub/axes. More... | |
Protected Member Functions | |
Coordinate () | |
Default constructor. More... | |
Coordinate (const Coordinate &other) | |
Copy constructor (copy semantics) More... | |
Coordinate & | operator= (const Coordinate &other) |
Assignment (copy semantics) More... | |
void | set_error (const String &errorMsg) const |
Set error message. More... | |
Bool | find_scale_factor (String &error, Vector< Double > &factor, const Vector< String > &units, const Vector< String > &oldUnits) |
void | fourierUnits (String &nameOut, String &unitOut, String &unitInCanon, Coordinate::Type type, Int axis, const String &unitIn, const String &nameIn) const |
Tries to find a canonical unit for input unit (e.g. More... | |
Bool | toWorldWCS (Vector< Double > &world, const Vector< Double > &pixel, wcsprm &wcs) const |
Functions to interconvert pixel<->world via wcs. More... | |
Bool | toPixelWCS (Vector< Double > &pixel, const Vector< Double > &world, wcsprm &wcs) const |
Bool | toWorldManyWCS (Matrix< Double > &world, const Matrix< Double > &pixel, Vector< Bool > &failures, wcsprm &wcs) const |
Bool | toPixelManyWCS (Matrix< Double > &pixel, const Matrix< Double > &world, Vector< Bool > &failures, wcsprm &wcs) const |
void | toCurrentMany (Matrix< Double > &world, const Vector< Double > &toCurrentFactors) const |
Functions for handling conversion between the current units and the wcs units. More... | |
void | fromCurrentMany (Matrix< Double > &world, const Vector< Double > &toCurrentFactors) const |
virtual void | convertTo (Vector< Double > &) const |
Functions for handling conversion between the current reference frame and the native one. More... | |
virtual void | convertFrom (Vector< Double > &) const |
void | convertToMany (Matrix< Double > &world) const |
Functions for handling conversion between the current reference frame and the native one for many conversions. More... | |
void | convertFromMany (Matrix< Double > &world) const |
void | pcToXform (Matrix< Double > &xForm, const wcsprm &wcs) const |
Interconvert between wcs PC cards and Matrix xForm format. More... | |
void | xFormToPC (wcsprm &wcs, const Matrix< Double > &xForm) const |
Protected Attributes | |
Vector< Double > | worldMin_p |
toMix ranges. More... | |
Vector< Double > | worldMax_p |
Private Member Functions | |
void | checkFormat (Coordinate::formatType &format, const Bool absolute) const |
Check format type. More... | |
void | makeWorldAbsRelMany (Matrix< Double > &value, Bool toAbs) const |
void | makePixelAbsRelMany (Matrix< Double > &value, Bool toAbs) const |
Private Attributes | |
String | error_p |
Interface for converting between world and pixel coordinates.
Public interface
The Coordinate class defines the generic interface whereby a pixel position is converted to a world (sky, frequency, stokes,...) position and vice versa. The pixel and world coordinates are in general multi-dimensional values. In general there need not be the same number of pixel and world axes, although this will normally be the case.
The fundamental model is that a pixel is first converted into a relative physical coordinate by:
After this linear stage, the final coordinate value is computed from this relative physical unit and a reference value, and possibly some other parameters. In the case of a sky position, these latter include at least the projection type. In the case of a purely linear coordinate, the reference value is merely added to the relative physical coordinate. The interface also allows the axes to be assigned names (reasonable defaults will be selected), and for physical units.
Both absolute and relative coordinates are supported. The main interface supports conversion between absolute pixel and absolute world coordinate. There are then functions to convert absolute coordinates to relative and vice versa. A relative pixel coordinate is defined according to
relative = absolute - reference
A relative world coordinate is similar, although there may be deviations from this formula (e.g. for DirectionCoordinate a cos(latitude) term is incorporated and for StokesCoordinate relative world coordinates are defined to be the same as absolute world coordinates.
Caution: All absolute pixels coordinates are zero relative;
This is a base class so there is no direct example, but see the example in Coordinates.h for use of the derived classes.
Encapsulate the common interface to coordinate conversion so that it may be used polymorphically.
Definition at line 139 of file Coordinate.h.
This enum is used for formatting world values into Strings.
Enumerator | |
---|---|
DEFAULT |
Default; formatter decides. |
SCIENTIFIC |
Scientific format (e.g. -1.2397E+03) |
FIXED |
Fixed floating format (e.g. 12.134) |
MIXED |
Either scientific or floating point, auto-selected by the C++ STL formatting routines. May not be available for all Coordinate types. |
TIME |
HHH:MM:SS.SSS style formatting. |
Definition at line 162 of file Coordinate.h.
This enum lists the types of the derived classes.
It is primarly used in the CoordinateSystem class.
Enumerator | |
---|---|
LINEAR |
Linear axes. |
DIRECTION |
A direction. Usually RA/DEC. |
SPECTRAL |
A spectral axis. |
STOKES |
A Stokes axis. |
TABULAR |
A one-dimensional Cooordinate system, usually created from a table although it can also be purely linear. |
QUALITY |
to mark DATA and ERROR values |
COORDSYS |
A CoordinateSystem (a collection of Coordinates). |
Definition at line 144 of file Coordinate.h.
|
virtual |
Destructor.
Needs to be public so the user can delete Coordinate* objects
|
protected |
Default constructor.
Make an empty coordinate. Used by derived classes.
|
protected |
Copy constructor (copy semantics)
|
private |
Check format type.
|
pure virtual |
Make a copy of ourself.
This pointer has been allocated with new
and must be deleted by the caller.
Implemented in casacore::CoordinateSystem, casacore::SpectralCoordinate, casacore::DirectionCoordinate, casacore::TabularCoordinate, casacore::StokesCoordinate, casacore::LinearCoordinate, and casacore::QualityCoordinate.
Reimplemented in casacore::DirectionCoordinate, and casacore::SpectralCoordinate.
Definition at line 549 of file Coordinate.h.
Functions for handling conversion between the current reference frame and the native one.
The default implementations do nothing. They should be over-ridden in the derived classes.
Reimplemented in casacore::DirectionCoordinate, and casacore::SpectralCoordinate.
Definition at line 547 of file Coordinate.h.
Functions for handling conversion between the current reference frame and the native one for many conversions.
These functions just call the virtual functions for single conversions.
|
static |
Call wcssub on the src/dst pair with null nsub/axes.
|
virtual |
Comparison only made for specified axes in this and other Coordinate The default implementation should be ok for all Coordinate types except Stokes and Quality...
Reimplemented in casacore::StokesCoordinate, and casacore::QualityCoordinate.
|
inline |
If the last conversion to world or pixel coordinates resulted in an error, report that error.
If the last conversion succeeded, it is undefined what this will return (it might well contain the last error message).
Definition at line 584 of file Coordinate.h.
References error_p.
|
protected |
String casacore::Coordinate::formatQuantity | ( | String & | units, |
Coordinate::formatType | format, | ||
const Quantum< Double > & | worldValue, | ||
uInt | axis, | ||
Bool | isAbsolute = True , |
||
Bool | showAsAbsolute = True , |
||
Int | precision = -1 |
||
) |
|
protected |
Tries to find a canonical unit for input unit (e.g.
GHz -> Hz), and tells you the output name and unit for the Fourier coordinate pairing with the canonical unit
|
protected |
|
virtual |
Provide a common interface to getting formatted representations of coordinate values.
Different derived Coordinate types are formatted
in different ways. For example, an RA/DEC DirectionCoordinate uses an HMS.SS/DMS.SS representation. A Galactic Lat/Long DirectionCoordinate uses floating format in degrees. Other derived Coordinates are formatted with scientific format or floating format. The derived class format functions provide this functionality.
You may specify the format with the format argument and a value from the enum Coordinate::formatType
. If you give it the value Coordinate::DEFAULT
then a sensible default is used.
A mechanism for specifying the precision number of significant digits after decimal point is provided. You can specify the precision directly when calling format if it is unambiguous how the derived Coordinate is going to be formatted. For example, a LinearCoordinate is always formatted with scientific format. However, if you are using these classes polymorphically, you don't want to have to know this and some derived Coordinates may be formatted in multiple ways (such as the DirectionCoordinate examples above). Therefore, the function getPrecision enables you to set default precisions for the different styles of formatting used variously in the base and derived classes. This function chooses the precision from these default values, according to the type of derived Coordinate that your object is and what value for format that you give (refer to the derived classes for details on this).
Some derived classes will format differently depending upon whether you want to format an absolute or offset world value input via absolute (e.g. DirectionCoordinates).
The provided worldValue
must be in the native units of the Coordinate. It may be an absolute (isAbsolute=True
) or relative (isAbsolute=False
) value. You may choose to format the world value as absolute (showAsAbsolute=True
) or relative (showAsAbsolute=False
). axis
specifies which axis of the Coordinate this value belongs to.
units
specifies the units in which the input world value will be formatted. If units
is empty, the native unit for the given axis is used.
Some derived classes will format in units different from the native unit of the Coordinate. The units of the formatted number are returned in units
. If the units
string is provided, the unit must be consistent with the native unit of the coordinate. The input world value will be converted to this unit.
You can also use the Quantum interface. The units of the Quantum can then be anything consistent with the Coordinate.
The default implementation here is to format only with scientific or fixed formats. If precision is negative, a the default precision is used.
Reimplemented in casacore::DirectionCoordinate.
|
static |
Call wcsini on the wcs structure.
|
virtual |
Find the Coordinate for when we Fourier Transform ourselves.
This pointer must be deleted by the caller. Axes specifies which axes of the Coordinate you wish to transform. Shape specifies the shape of the image associated with all the axes of the Coordinate. Currently the output reference pixel is always shape/2.
Reimplemented in casacore::CoordinateSystem, casacore::DirectionCoordinate, casacore::SpectralCoordinate, casacore::TabularCoordinate, and casacore::LinearCoordinate.
Reimplemented in casacore::CoordinateSystem, casacore::TabularCoordinate, casacore::StokesCoordinate, and casacore::QualityCoordinate.
Reimplemented in casacore::CoordinateSystem.
|
private |
Make absolute coordinates relative and vice-versa (with respect to the reference value).
Vectors must be length nPixelAxes()
or nWorldAxes()
or memory access errors will occur
Reimplemented in casacore::CoordinateSystem, casacore::TabularCoordinate, casacore::StokesCoordinate, and casacore::QualityCoordinate.
Batch up a lot of absolute/relative transformations.
Parameters as above for toWorldMany
and toPixelMany
Reimplemented in casacore::CoordinateSystem.
Reimplemented in casacore::CoordinateSystem.
|
virtual |
Make absolute coordinates relative and vice versa with respect to the given reference value.
Add the other functions in this grouping as needed. Vectors must be length nPixelAxes()
or nWorldAxes()
or memory access errors will occur
Reimplemented in casacore::CoordinateSystem, and casacore::DirectionCoordinate.
|
private |
Reimplemented in casacore::CoordinateSystem.
|
pure virtual |
Comparison to fractional tolerance (for floating point values).
Don't compare on specified axes in Coordinate. If the comparison returns False, errorMessage()
contains a message.
Implemented in casacore::CoordinateSystem, casacore::DirectionCoordinate, casacore::SpectralCoordinate, casacore::TabularCoordinate, casacore::StokesCoordinate, casacore::LinearCoordinate, and casacore::QualityCoordinate.
|
pure virtual |
|
pure virtual |
How many world/pixel axes are there in this Coordinate? While the number of world and pixel axes will generally be the same, it is not a requirement.
For example, in CoordinateSystem you could remove a pixel axis and leave the corresponding world axis. Also, if we ever implement a "SlicedCoordinate" class then there would be more world than pixel coordinates (the pixel coordinate would be a pixel number along the slice, whereas the world axes would continue to be RA/DEC).
Implemented in casacore::CoordinateSystem, casacore::DirectionCoordinate, casacore::SpectralCoordinate, casacore::TabularCoordinate, casacore::LinearCoordinate, casacore::StokesCoordinate, and casacore::QualityCoordinate.
|
pure virtual |
|
protected |
Assignment (copy semantics)
|
protected |
Interconvert between wcs PC cards and Matrix xForm format.
|
virtual |
return the result of rotating the coordinate clockwise through the specified angle.
Rotation occurs about the reference pixel. Coordinate must have exactly two pixel axes. The return type is the same as the input type. It is the caller's responsibility to delete the returned pointer when done with it to prevent a memory leak. This method ultimately just changes the input coordinate's linear transform matrix.
|
pure virtual |
Used for persistence.
Derived classes will have similar static restore methods. It will typically only return False if fieldName has already been defined.
Implemented in casacore::CoordinateSystem, casacore::SpectralCoordinate, casacore::DirectionCoordinate, casacore::TabularCoordinate, casacore::StokesCoordinate, casacore::LinearCoordinate, and casacore::QualityCoordinate.
|
protected |
Set error message.
|
static |
Call wcsset on the wcs structure.
|
virtual |
Reimplemented in casacore::CoordinateSystem, casacore::DirectionCoordinate, casacore::StokesCoordinate, and casacore::QualityCoordinate.
Set the requested attribute.
Note that these just change the internal values, they do not cause any recomputation.
Implemented in casacore::CoordinateSystem, casacore::DirectionCoordinate, casacore::SpectralCoordinate, casacore::TabularCoordinate, casacore::LinearCoordinate, casacore::StokesCoordinate, and casacore::QualityCoordinate.
Change the units.
Adjust the increment and reference value by the ratio of the old and new units. This implies that the units must be known Unit strings, and that they must be compatible, e.g. they can't change from time to length.
A default implementation is available which does everything except set the units vector, which must be done in the derived class.
Implemented in casacore::CoordinateSystem, casacore::DirectionCoordinate, casacore::SpectralCoordinate, casacore::TabularCoordinate, casacore::LinearCoordinate, casacore::StokesCoordinate, and casacore::QualityCoordinate.
Set the world min and max ranges, for use in function toMix
, for a lattice of the given shape for this coordinate.
The default implementation here sets the range for pixels dangling 25% off the image. Returns False if fails with a reason in errorMessage()
. setDefaultWorldMixRanges sets the range for each axis to +/-1e99 The ranges remain zero length vectors until you explicitly initialize them.
Reimplemented in casacore::CoordinateSystem, casacore::DirectionCoordinate, casacore::StokesCoordinate, and casacore::QualityCoordinate.
|
pure virtual |
|
static |
Call wcssub on the src/dst pair.
|
protected |
Functions for handling conversion between the current units and the wcs units.
These are called explicitly by the appropriate derived class. convertFrom
|
virtual |
Mixed absolute pixel/world coordinate conversion.
worldIn and worldAxes are vectors of length nWorldAxes
. pixelIn
and pixelAxes
are of length nPixelAxes
. worldAxes(i) = True
specifies you have given a world value in worldIn(i)
to convert to pixel. pixelAxes(i)=True
specifies you have given a pixel value in pixelIn(i)
to convert to world. You cannot specify the same axis via worldAxes
and pixelAxes
. Values in pixelIn
are converted to world and put into worldOut
in the appropriate world axis location. Values in worldIn
are copied to worldOut
. Values in worldIn
are converted to pixel and put into pixelOut
in the appropriate pixel axis location. Values in pixelIn
are copied to pixelOut
. worldMin
and worldMax
specify the range of the world coordinate (in the world axis units of that world axis in the CoordinateSystem) being solved for in a mixed calculation for each world axis. They are only actually needed for DirectionCoordinates and for all other Coordinates the relevant elements can be undefined. If you don't know, use -180 to 180 degrees for longitude, and -90 to 90 for latitude. Removed axes are handled (for example, a removed pixel axis with remaining corresponding world axis will correctly be converted to world using the replacement value). Returns True if the conversion succeeds, otherwise it returns False and errorMessage()
contains an error message. The output vectors are resized.
Reimplemented in casacore::CoordinateSystem, and casacore::DirectionCoordinate.
|
virtual |
Reimplemented in casacore::CoordinateSystem, casacore::DirectionCoordinate, casacore::SpectralCoordinate, and casacore::TabularCoordinate.
|
protected |
|
protected |
|
pure virtual |
Convert an absolute pixel position to an absolute world position or vice versa.
Returns True if the conversion succeeds, otherwise it returns False and method errorMessage contains an error message. The input vector must be of length nPixelAxes
or nWorldAxes
. The output vector is resized appropriately. if useConversionFrame
, if the coordinate has a conversion layer frame (such as can be present in spectral and direction coordinates), it is used. Else, the native frame is used for the conversion.
Implemented in casacore::CoordinateSystem, casacore::DirectionCoordinate, casacore::SpectralCoordinate, casacore::TabularCoordinate, casacore::LinearCoordinate, casacore::StokesCoordinate, and casacore::QualityCoordinate.
|
virtual |
Batch up a lot of transformations.
The first (most rapidly varying) axis of the matrices contain the coordinates. Returns False if any conversion failed and errorMessage()
will hold a message. The failures
array (True for fail, False for success) is the length of the number of conversions and holds an error status for each conversion. The default implementation is provided that works with the "single" version of toWorld
and toPixel
, but for maximum efficiency these should be overridden.
Reimplemented in casacore::CoordinateSystem, casacore::DirectionCoordinate, casacore::SpectralCoordinate, and casacore::TabularCoordinate.
|
protected |
|
protected |
Functions to interconvert pixel<->world via wcs.
These functions are called explicitly by the to{world,Pixel} functions in the appropriate wcs-based derived classes.
|
pure virtual |
List the type of this Coordinate object.
Implemented in casacore::CoordinateSystem, casacore::DirectionCoordinate, casacore::SpectralCoordinate, casacore::TabularCoordinate, casacore::LinearCoordinate, casacore::StokesCoordinate, and casacore::QualityCoordinate.
|
static |
Return the requested attributed.
Implemented in casacore::CoordinateSystem, casacore::DirectionCoordinate, casacore::SpectralCoordinate, casacore::TabularCoordinate, casacore::LinearCoordinate, casacore::StokesCoordinate, and casacore::QualityCoordinate.
Definition at line 265 of file Coordinate.h.
References worldMax_p.
Definition at line 264 of file Coordinate.h.
References worldMin_p.
|
protected |
|
mutableprivate |
Definition at line 570 of file Coordinate.h.
Referenced by errorMessage().
Definition at line 567 of file Coordinate.h.
Referenced by worldMixMax().
toMix ranges.
Should be set by derived class.
Definition at line 567 of file Coordinate.h.
Referenced by worldMixMin().