AnnRegion.h

Go to the documentation of this file.

//# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
//# FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Library General Public
//# License for more details.
//#
//# You should have received a copy of the GNU Library General Public License
//# along with this library; if not, write to the Free Software Foundation,
//# Inc., 675 Massachusetts Ave, Cambridge, MA 02139, USA.
//#
//# Correspondence concerning AIPS++ should be addressed as follows:
//#        Internet email: aips2-request@nrao.edu.
//#        Postal address: AIPS++ Project Office
//#                        National Radio Astronomy Observatory
//#                        520 Edgemont Road
//#                        Charlottesville, VA 22903-2475 USA
//#

#ifndef REGIONS_ANNREGION_H
#define REGIONS_ANNREGION_H

#include <coordinates/Coordinates/CoordinateSystem.h>
#include <imageanalysis/Annotations/AnnotationBase.h>
#include <images/Regions/WCBox.h>
#include <images/Regions/ImageRegion.h>
#include <measures/Measures/MDirection.h>
#include <measures/Measures/MFrequency.h>

namespace casa {

// <summary>
// This class represents a annotation referring to a region specified in an ascii region file as proposed
// in CAS-2285
// </summary>
// <author>Dave Mehringer</author>
// <use visibility=export>
// <reviewed reviewer="" date="yyyy/mm/dd" tests="" demos="">
// </reviewed>
// <prerequisite>

// </prerequisite>

// <etymology>
// Holds the specification of an annotation containing a region as specified in ASCII format.
// </etymology>

// <synopsis>
// This class is meant to be a container for all parameters necessary to specify a region per
// the format proposal attached to CAS-2285 (https://bugs.nrao.edu/browse/CAS-2285).

// Conversions of frequency from one reference frame to another are done here. The position of
// the reference pixel in the supplied coordinate system is used when converting frequencies.
// </synopsis>

class AnnRegion: public AnnotationBase {

public:

        virtual ~AnnRegion();

        void setAnnotationOnly(const Bool isAnnotationOnly);

        // is this region an annotation only? ie just for graphical rendering?
        Bool isAnnotationOnly() const;

        virtual TableRecord asRecord() const;

        virtual ImageRegion asImageRegion() const;

        virtual WCRegion* getRegion() const;

        // returns True unless overridden.
        virtual Bool isRegion() const;

        void setDifference(const Bool difference);

        Bool isDifference() const;

        // get the pixel range included in the spectral selection.
        // If there is no spectral axis, a zero length vector is returned. Otherwise,
        // a vector of two values is returned. The zeroth value will always be less
        // than or equal to the first.
        vector<Double> getSpectralPixelRange() const;



protected:

        // only to be called by subclasses

        // beginFreq and endFreq can both be 0, in which case
        // all the spectral range is used if a spectral axis exists in
        // <src>csys</csys>. If one of <src>beginFreq</src> or <src>endFreq</src>
        // is given, the other must be given. Frequency units can either conform
        // to Hz, m/s, or pix. If <src>beginFreq</src> and <src>endFreq</src> are not
        // specifed or if they are specified in pixel units,
        // <src>freqRefFrame</src>, <src>dopplerString</src>, and
        // <src>restfreq</src> are ignored. If provided, <src>beginFreq</src>
        // and <src>endFreq</src> must conform to the same units.
        AnnRegion(
                const Type shape,
                const String& dirRefFrameString,
                const CoordinateSystem& csys,
                const IPosition& imShape,
                const Quantity& beginFreq,
                const Quantity& endFreq,
                const String& freqRefFrame,
                const String& dopplerString,
                const Quantity& restfreq,
                const Vector<Stokes::StokesTypes> stokes,
                const Bool annotationOnly
        );

        // use if all coordinate values will be specified in
        // the same frames as the input coordinate system. frequencies
        // and the annotationOnly flag can be set after
        // construction. By default, all frequencies and all polarizations
        // are used, and the annotationOnly flag is False
        AnnRegion(
                const Type shape,
                const CoordinateSystem& csys,
                const IPosition& imShape,
                const Vector<Stokes::StokesTypes>& stokes
        );

        // implicitly defined copy constructor is fine
        // AnnRegion(AnnRegion& other);

        // copy constructor
        AnnRegion(const AnnRegion& other);

        // assignment operator
        AnnRegion& operator= (const AnnRegion& rhs);

        // extend the direction plane region over spectral and/or polarization
        // coordinates
        void _extend();

        void _toRecord(const ImageRegion& region);

        // convert a length in pixels to an angle.
        Quantity _lengthToAngle(
                const Quantity& quantity, const uInt pixelAxis
        ) const;

        virtual void _printPrefix(ostream& os) const;

        // subclasses must call this at construction to set their base region
        // defined in the direction plane
        void _setDirectionRegion(const ImageRegion& region);

        Bool _setFrequencyLimits(
                const Quantity& beginFreq,
                const Quantity& endFreq,
                const String& freqRefFrame,
                const String& dopplerString,
                const Quantity& restfreq
        );


private:

        Bool _isAnnotationOnly;
        Bool _isDifference, _constructing;
        ImageRegion _imageRegion, _directionRegion;
        IPosition _imShape;
        vector<Double> _spectralPixelRange;

        static const String _class;


        WCBox _makeExtensionBox(
                const Vector<Quantity>& freqRange,
                const Vector<Stokes::StokesTypes>& stokesRange,
                const IPosition& pixelAxes
        ) const;

        void _init();

};

// Just need a identifable expection class, compiler can generate implementation implicitly
class ToLCRegionConversionError : public AipsError {
public:
        ToLCRegionConversionError(String msg) : AipsError(msg) {}
};

}



#endif