JUTMPosition.hh

Go to the documentation of this file.

#ifndef __JUTM__JUTMPOSITION__
#define __JUTM__JUTMPOSITION__
 
#include <istream>
#include <ostream>
 
#include "JGeometry3D/JPosition3D.hh"
#include "JMath/JMath.hh"
#include "JIO/JSerialisable.hh"
 
 
/**
 * \author mdejong
 */
 
namespace JUTM {}
namespace JPP { using namespace JUTM; }
 
/**
 * Auxiliaries for handling universal transverse mercator coordinate system (UTM).
 */
namespace JUTM {
 
  using JGEOMETRY3D::JVector3D;
  using JGEOMETRY3D::JPosition3D;
  using JIO::JReader;
  using JIO::JWriter;
 
 
  /**
   * Data structure for UTM position.
   * The z-coordinate corresponds to the depth where <tt>z = 0</tt> is at sea level.
   */
  class JUTMPosition :
    public JMATH::JMath<JUTMPosition>
  {
  public:
    /**
     * Default constructor.
     */
    JUTMPosition() :
      east (0.0),
      north(0.0),
      z    (0.0)
    {}
 
 
    /**
     * Constructor.
     *
     * \param  pos              UTM position
     */
    JUTMPosition(const JVector3D& pos) :
      east (pos.getX()),
      north(pos.getY()),
      z    (pos.getZ())
    {}
 
 
    /**
     * Constructor.
     *
     * \param  east             UTM East
     * \param  north            UTM North
     * \param  z                UTM Z
     */
    JUTMPosition(const double east,
                 const double north,
                 const double z) 
    {
      this->east  = east;
      this->north = north;
      this->z     = z;
    }
 
    
    /**
     * Get UTM position.
     *
     * \return                  position
     */
    const JUTMPosition& getUTMPosition() const
    {
      return static_cast<const JUTMPosition&>(*this);
    }
 
 
    /**
     * Get position.
     *
     * \return                  position
     */
    JPosition3D getPosition() const
    {
      return JPosition3D(east, north, z);
    }
 
 
    /**
     * Get UTM east.
     *     
     * \return                  UTM East
     */
    double getUTMEast() const
    {
      return this->east;
    }
 
 
    /**
     * Get UTM north.
     *     
     * \return                  UTM North
     */
    double getUTMNorth() const
    {
      return this->north;
    }
 
 
    /**
     * Get UTM Z.
     *     
     * \return                  UTM Z
     */
    double getUTMZ() const
    {
      return this->z;
    }
 
 
    /**
     * Negate UTM position.
     *
     * \return               this UTM position
     */
    JUTMPosition& negate()
    {
      east  = -east;
      north = -north;
      z     = -z;
 
      return *this;
    }
 
 
    /**
     * Add UTM position.
     *
     * \param  pos           UTM position
     * \return               this UTM position
     */
    JUTMPosition& add(const JUTMPosition& pos)
    {
      east  += pos.getUTMEast();
      north += pos.getUTMNorth();
      z     += pos.getUTMZ();
 
      return *this;
    }
 
 
    /**
     * Subtract UTM position.
     *
     * \param  pos           UTM position
     * \return               this UTM position
     */
    JUTMPosition& sub(const JUTMPosition& pos)
    {
      east  -= pos.getUTMEast();
      north -= pos.getUTMNorth();
      z     -= pos.getUTMZ();
 
      return *this;
    }
 
 
    /**
     * Scale UTM position.
     *
     * \param  factor        multiplication factor
     * \return               this UTM position
     */
    JUTMPosition& mul(const double factor)
    {
      east  *= factor;
      north *= factor;
      z     *= factor;
 
      return *this;
    }
 
 
    /**
     * Scale UTM position.
     *
     * \param  factor        division factor
     * \return               this UTM position
     */
    JUTMPosition& div(const double factor)
    {
      east  /= factor;
      north /= factor;
      z     /= factor;
 
      return *this;
    }
 
 
    /**
     * Read UTM position from input.
     *
     * \param  in               input stream
     * \param  pos              UTM position
     * \return                  input stream
     */
    friend inline std::istream& operator>>(std::istream& in, JUTMPosition& pos)
    {
      return in >> pos.east >> pos.north >> pos.z;
    }
    
 
    /**
     * Write UTM position to output.
     *
     * \param  out              output stream
     * \param  pos              UTM position
     * \return                  output stream
     */
    friend inline std::ostream& operator<<(std::ostream& out, const JUTMPosition& pos)
    {
      return out << pos.east << ' ' << pos.north << ' ' << pos.z;
    }
 
 
 
    /**
     * Read UTM position from input.
     *
     * \param  in               input stream
     * \param  pos              UTM position
     * \return                  input stream
     */
    friend inline JReader& operator>>(JReader& in, JUTMPosition& pos)
    {
      return in >> pos.east >> pos.north >> pos.z;
    }
    
 
    /**
     * Write UTM position to output.
     *
     * \param  out              output stream
     * \param  pos              UTM position
     * \return                  output stream
     */
    friend inline JWriter& operator<<(JWriter& out, const JUTMPosition& pos)
    {
      return out << pos.east << pos.north << pos.z;
    }
 
 
  protected:
    double east;
    double north;
    double z;
  };
}
 
#endif