VcstMapDataSet.h

Go to the documentation of this file.

/**************************************************************************
 * 
 * NAME: VcstMapDataSet
 *
 **************************************************************************/
 
#ifndef _VcstMapDataSet_H_
#define _VcstMapDataSet_H_
 
#include <algorithm>
#include <VcstCmnConsts.h>
#include <mds.h>
#include <novas.h>
 
//------------------------------------------------------------------------------
// Constants for instantiating N. & S. Mds.  The values for the lat and lon are
// just in the opposite hemisphere that is being used, so that the MDS
// encapsulates the entire hemisphere.
//------------------------------------------------------------------------------
const double NORTH_UPPER_LEFT_LAT = (DEG2RAD * -20.747206768662e0);
const double NORTH_UPPER_LEFT_LON = (DEG2RAD * 145.0e0);
const double NORTH_LOWER_RIGHT_LAT = NORTH_UPPER_LEFT_LAT;
const double NORTH_LOWER_RIGHT_LON = (NORTH_UPPER_LEFT_LON - PI);
const double SOUTH_UPPER_LEFT_LAT = -NORTH_UPPER_LEFT_LAT;
const double SOUTH_UPPER_LEFT_LON = (DEG2RAD * (-125.0e0));
//const double SOUTH_LOWER_RIGHT_LAT  =  SOUTH_UPPER_LEFT_LAT;
const double SOUTH_LOWER_RIGHT_LON = (SOUTH_UPPER_LEFT_LON + PI);
 
//------------------------------------------------------------------------------
// Constant used to allow for small float-to-double conversion errors
// when checking values against PI.
//------------------------------------------------------------------------------
const double PI_EPSILON = 1.0e-7;
 
enum Hemisphere {
    NORTHERN, SOUTHERN, EQUATOR
};
 
enum GridType {
    GRIDTYPE_NONE = 0,
    MERCATOR_COAX = 11,
    POLAR_STEREO_NORTH = 21,
    POLAR_STEREO_SOUTH = 25,
    LAMBERT_TANGENT_NORTH = 31,
    LAMBERT_TANGENT_SOUTH = 35,
    LAMBERT_SECANT_NORTH = 41,
    LAMBERT_SECANT_SOUTH = 45,
    CYLINDRICAL_EQUIDISTANT = 51,
    POLAR_AZMUTHAL_NORTH = 61,
    POLAR_AZMUTHAL_SOUTH = 65
};
 
class VcstMapDataSet {
public:
 
    VcstMapDataSet(Hemisphere theHemisphere);
 
    virtual ~VcstMapDataSet();
 
    void fillMds(mds_type * const anMds) const;
 
    /* inline */
    double getUpperLeftLatitude() const;
 
    /* inline */
    double getUpperLeftLongitude() const;
 
    /* inline */
    double getLowerRightLatitude() const;
 
    /* inline */
    double getLowerRightLongitude() const;
 
    /* inline */
    double getBiggestLatitude() const;
 
    /* inline */
    double getBiggestLongitude() const;
 
    /* inline */
    double getSmallestLatitude() const;
 
    /* inline */
    double getSmallestLongitude() const;
 
    /* inline */
    Hemisphere getHemisphere() const;
 
    /* inline */
    int getGridType() const;
 
    virtual int latLonToGrid(
            double rlat,
            double rlon,
            double* row,
            double* column) = 0;
 
    virtual int gridToLatLon(
            double row,
            double column,
            double* rlat,
            double* rlon) = 0;
 
 
protected:
 
    GridType gridType_;
 
    Hemisphere hemisphere_;
 
    short wedgeRotation_;
 
    int mdsNum_;
 
    double stanLat1_;
 
    double stanLat2_;
 
    double baseLongitude_;
 
    double gridIncConstant_;
 
    double maxRow_;
 
    double maxCol_;
 
    double upperLeftLatitude_;
 
    double upperLeftLongitude_;
 
    double lowerRightLatitude_;
 
    double lowerRightLongitude_;
 
    double upperLeftX_;
 
    double upperLeftY_;
 
    double lowerRightX_;
 
    double lowerRightY_;
 
    double smallestLatitude_;
 
    double biggestLatitude_;
 
    double smallestLongitude_[2];
 
    double biggestLongitude_[2];
 
    short numberOfRanges_;
 
    int checkLat(double &rlat);
 
    int checkLon(double &rlon);
 
private:
 
};
 
//-----------------------------------------------------------------------------
 
inline double
VcstMapDataSet::getUpperLeftLatitude() const {
    return upperLeftLatitude_;
}
 
//-----------------------------------------------------------------------------
 
inline double
VcstMapDataSet::getUpperLeftLongitude() const {
    return upperLeftLongitude_;
}
 
//-----------------------------------------------------------------------------
 
inline double
VcstMapDataSet::getLowerRightLatitude() const {
    return lowerRightLatitude_;
}
 
//-----------------------------------------------------------------------------
 
inline double
VcstMapDataSet::getLowerRightLongitude() const {
    return lowerRightLongitude_;
}
 
//-----------------------------------------------------------------------------
 
inline double
VcstMapDataSet::getBiggestLatitude() const {
    return biggestLatitude_;
}
 
//-----------------------------------------------------------------------------
 
inline double
VcstMapDataSet::getSmallestLatitude() const {
    return smallestLatitude_;
}
 
//-----------------------------------------------------------------------------
 
inline double
VcstMapDataSet::getBiggestLongitude() const {
    if (numberOfRanges_ == 1) {
        return biggestLongitude_[0];
    } else {
        return std::max(biggestLongitude_[0], biggestLongitude_[1]);
    }
}
 
//-----------------------------------------------------------------------------
 
inline double
VcstMapDataSet::getSmallestLongitude() const {
    if (numberOfRanges_ == 1) {
        return smallestLongitude_[0];
    } else {
        return std::min(smallestLongitude_[0], smallestLongitude_[1]);
    }
}
 
//-----------------------------------------------------------------------------
 
inline Hemisphere
VcstMapDataSet::getHemisphere() const {
    return hemisphere_;
}
 
//-----------------------------------------------------------------------------
 
inline int
VcstMapDataSet::getGridType() const {
    return gridType_;
}
 
 
#endif