testUTMPoint.cc

Go to the documentation of this file.

#include "Geometry/testUTMPointData.h"
#include <utl/UTMPoint.h>
#include <utl/Point.h>
#include <utl/ReferenceEllipsoid.h>

#include <tst/Verify.h>
#include <utl/Triple.h>

#include <cppunit/extensions/HelperMacros.h>

#include <iostream>

using namespace utl;
using namespace tst;

#define TEST_CLOSETO3(n1_, n2_, e1_, e2_, h1_, h2_) \
  CPPUNIT_ASSERT(Verify<CloseTo>(n1_, n2_)); \
  CPPUNIT_ASSERT(Verify<CloseTo>(e1_, e2_)); \
  CPPUNIT_ASSERT(Verify<CloseTo>(h1_, h2_))


class testUTMPoint : public CppUnit::TestFixture {

  CPPUNIT_TEST_SUITE(testUTMPoint);
  CPPUNIT_TEST(testConstructors);
  CPPUNIT_TEST(testAssignement);
  CPPUNIT_TEST(testLatLongConversions);
  CPPUNIT_TEST(testUTMConversions);
  CPPUNIT_TEST_SUITE_END();

private:
  void
  oneLatLongConversion(const double la, const double lo, const double n, const double e,
                       const double h, const int z, const char b)
  {
    const UTMPoint p(la, lo, h, ReferenceEllipsoid::eWGS84);
    CPPUNIT_ASSERT(Verify<CloseTo>(
      p.GetGeodeticCoordinates(), Triple(la, lo, h)));
    CPPUNIT_ASSERT(p.GetZone() == z);
    CPPUNIT_ASSERT(p.GetBand() == b);
    /*TEST_CLOSETO3(p.GetNorthing(), n,
                  p.GetEasting(), e,
                  p.GetHeight(), h);*/
    CPPUNIT_ASSERT(Verify<CloseTo>(
      p.GetCoordinates3(), Triple(n, e, h)));

    const UTMPoint q(n, e, h, z, b, ReferenceEllipsoid::eWGS84);
    CPPUNIT_ASSERT(Verify<CloseTo>(
      q.GetGeodeticCoordinates(), Triple(la, lo, h)));
    CPPUNIT_ASSERT(q.GetZone() == z);
    CPPUNIT_ASSERT(q.GetBand() == b);
    /*TEST_CLOSETO3(q.GetNorthing(), n,
                  q.GetEasting(), e,
                  q.GetHeight(), h);*/
    CPPUNIT_ASSERT(Verify<CloseTo>(
      q.GetCoordinates3(), Triple(n, e, h)));
  }

  void
  oneUTMConversion(const double /*la*/, const double /*lo*/, const double n, const double e,
                   const double h, const int z, const char b)
  {
    const UTMPoint up(n, e, h, z, b, ReferenceEllipsoid::eWGS84);
    const Point cp(up.GetPoint());
    const UTMPoint uq(cp, ReferenceEllipsoid::eWGS84);

    CPPUNIT_ASSERT(Verify<CloseTo>(
      up.GetCoordinates3(), uq.GetCoordinates3()));
    CPPUNIT_ASSERT(up.GetZone() == uq.GetZone());
    CPPUNIT_ASSERT(up.GetBand() == uq.GetBand());
  }

public:
  void setUp() { }

  void tearDown() { }

  void dummy(UTMPoint /*p*/) { };

  void
  testConstructors()
  {
    const ReferenceEllipsoid e = ReferenceEllipsoid::Get(ReferenceEllipsoid::eWGS84);
    const UTMPoint p1(10000, 500000, 0, 10, 'I', e);
    CPPUNIT_ASSERT(p1.GetNorthing() == 10000);
    CPPUNIT_ASSERT(p1.GetEasting() == 500000);
    CPPUNIT_ASSERT(p1.GetHeight() == 0);
    CPPUNIT_ASSERT(p1.GetZone() == 10);
    CPPUNIT_ASSERT(p1.GetBand() == 'I');
    CPPUNIT_ASSERT(p1.GetEllipsoid() == e);
    CPPUNIT_ASSERT(p1.GetCoordinates() ==
                   boost::make_tuple(10000, 500000, 0, 10., 'I'));
    CPPUNIT_ASSERT(p1.GetCoordinates4() ==
                   boost::make_tuple(10000, 500000, 0, 10));
    CPPUNIT_ASSERT(p1.GetCoordinates3() ==
                   Triple(10000, 500000, 0));

    const UTMPoint p2(9345711., 821935., 1432., 7, 'L', ReferenceEllipsoid::eWGS84);
    CPPUNIT_ASSERT(p2.GetCoordinates() ==
                   boost::make_tuple(9345711., 821935., 1432., 7, 'L'));
    CPPUNIT_ASSERT(p2.GetEllipsoid() == e);

    const UTMPoint p3(p2, e);
    CPPUNIT_ASSERT(p2 == p3);
  }

  void
  testAssignement()
  {
    const ReferenceEllipsoid e = ReferenceEllipsoid::Get(ReferenceEllipsoid::eWGS84);
    const UTMPoint p0(10000, 500000, 0, 10, 'I', e);
    UTMPoint p1(0, 0, 100, 11, 'J', e);

    CPPUNIT_ASSERT(p1 != p0);

    p1 = p0;
    CPPUNIT_ASSERT(p1 == p0);
    CPPUNIT_ASSERT(p1.GetNorthing() == 10000);
    CPPUNIT_ASSERT(p1.GetEasting() == 500000);
    CPPUNIT_ASSERT(p1.GetHeight() == 0);
    CPPUNIT_ASSERT(p1.GetZone() == 10);
    CPPUNIT_ASSERT(p1.GetBand() == 'I');
    CPPUNIT_ASSERT(p1.GetEllipsoid() == e);
    CPPUNIT_ASSERT(p1.GetCoordinates() ==
                   boost::make_tuple(10000, 500000, 0, 10., 'I'));
    CPPUNIT_ASSERT(p1.GetCoordinates4() ==
                   boost::make_tuple(10000, 500000, 0, 10));
    CPPUNIT_ASSERT(p1.GetCoordinates3() ==
                   Triple(10000, 500000, 0));

  }

  void
  testLatLongConversions()
  {
    for (LatLongNorthEast* p = llneData;
         p < llneData + sizeof(llneData)/sizeof(LatLongNorthEast);
         ++p) {
      oneLatLongConversion(p->fLat, p->fLong, p->fNorth, p->fEast,
                           testHeight, testZone, testBand);
    }
  }

  void
  testUTMConversions()
  {
    // Note: in all these tests, be aware that rounding errors near
    // zone or band boundaries can place the point you convert from
    // lat/long to UTM into the neighbouring zone or band. These
    // errors cancel when you convert back. The resulting UTM
    // representation is valid, since the UTM zones have some
    // overlap.
    const CoordinateSystemPtr cs = CoordinateSystem::GetRootCoordinateSystem();
    const ReferenceEllipsoid& e = ReferenceEllipsoid::Get(ReferenceEllipsoid::eWGS84);

    const UTMPoint up1(0*deg, 0*deg, 0*m, ReferenceEllipsoid::eWGS84);
    CPPUNIT_ASSERT(up1.GetZone() == 31);
    CPPUNIT_ASSERT(up1.GetBand() == 'N');

    /*TEST_CLOSETO3(up1.GetNorthing(), 0*m,
                  up1.GetEasting(), 166021.44*m,
                  up1.GetHeight(), 0*m);*/
    CPPUNIT_ASSERT(Verify<CloseTo>(
      up1.GetCoordinates3(), Triple(0*m, 166021.44*m, 0*m)));
    CPPUNIT_ASSERT(Verify<CloseTo>(
      up1.GetGeodeticCoordinates(), Triple(0.*deg, 0.*deg, 0.*m)));
    const Point p1(up1.GetPoint());
    CPPUNIT_ASSERT(Verify<CloseTo>(
      p1.GetCoordinates(cs), Triple(e.GetEquatorialRadius(), 0, 0)));

    // made longitude slightly larger than 90 deg to avoid rounding problems
    const UTMPoint up2(0*deg, 90.0000000001*deg, 0*m, ReferenceEllipsoid::eWGS84);
    CPPUNIT_ASSERT(up2.GetZone() == 46);
    CPPUNIT_ASSERT(up2.GetBand() == 'N');
    CPPUNIT_ASSERT(Verify<CloseTo>(
      up2.GetCoordinates3(), Triple(0*m, 166021.44*m, 0*m)));
    CPPUNIT_ASSERT(Verify<CloseTo>(
      up2.GetGeodeticCoordinates(), Triple(0.*deg, 90.*deg, 0.*m)));
    const Point p2(up2.GetPoint());
    CPPUNIT_ASSERT(Verify<CloseTo>(
      p2.GetCoordinates(cs), Triple(0., e.GetEquatorialRadius(), 0.)));

    const UTMPoint up3(1*deg, 90*deg, 0*m, ReferenceEllipsoid::eWGS84);
    CPPUNIT_ASSERT(up3.GetZone() == 45);
    CPPUNIT_ASSERT(up3.GetBand() == 'N');
    CPPUNIT_ASSERT(Verify<CloseTo>(
      up3.GetGeodeticCoordinates(), Triple(1.*deg, 90.*deg, 0.*m)));

    const UTMPoint up4(-1*deg, 90*deg, 0*m, ReferenceEllipsoid::eWGS84);
    CPPUNIT_ASSERT(up4.GetZone() == 45);
    CPPUNIT_ASSERT(up4.GetBand() == 'M');
    CPPUNIT_ASSERT(Verify<CloseTo>(
      up4.GetGeodeticCoordinates(), Triple(-1.*deg, 90.*deg, 0.*m)));

    const UTMPoint up5(0*deg, 91*deg, 0*m, ReferenceEllipsoid::eWGS84);
    CPPUNIT_ASSERT(up5.GetZone() == 46);
    CPPUNIT_ASSERT(up5.GetBand() == 'N');
    CPPUNIT_ASSERT(Verify<CloseTo>(
      up5.GetGeodeticCoordinates(), Triple(0.*deg, 91.*deg, 0.*m)));

    for (LatLongNorthEast* p = llneData;
         p < llneData + sizeof(llneData)/sizeof(LatLongNorthEast);
         ++p) {
      oneUTMConversion(p->fLat, p->fLong, p->fNorth, p->fEast,
                             testHeight, testZone, testBand);
    }
  }

};


CPPUNIT_TEST_SUITE_REGISTRATION(testUTMPoint);


// Configure (x)emacs for this file ...
// Local Variables:
// mode: c++
// compile-command: "make -C .. -k testGeometry && ../testGeometry "
// End: