testReferenceEllipsoid.cc

Go to the documentation of this file.

#include <utl/ReferenceEllipsoid.h>
#include <utl/Point.h>
#include <utl/AugerUnits.h>

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

#include <cppunit/extensions/HelperMacros.h>

using namespace utl;
using namespace tst;


class testReferenceEllipsoid : public CppUnit::TestFixture {

  CPPUNIT_TEST_SUITE(testReferenceEllipsoid);
  CPPUNIT_TEST(testEllipsoidCreation);
  CPPUNIT_TEST_EXCEPTION(testGetEllipsoid0,
                         ReferenceEllipsoid::InvalidEllipsoidException);
  CPPUNIT_TEST_EXCEPTION(testGetEllipsoid999, AugerException);
  CPPUNIT_TEST(testEllipsoidComparison);
  CPPUNIT_TEST(testLatLongConversion);
  CPPUNIT_TEST_SUITE_END();

public:
  void setUp() { }

  void tearDown() { }

  void
  testEllipsoidCreation()
  {
    const ReferenceEllipsoid r = ReferenceEllipsoid::Get(ReferenceEllipsoid::eWGS84);
    CPPUNIT_ASSERT(Verify<CloseTo>(r.GetEquatorialRadius(), 6378137.));
    CPPUNIT_ASSERT(Verify<CloseTo>(r.GetEccentricity2(), 0.00669437999013));
    CPPUNIT_ASSERT(Verify<CloseTo>(r.GetPolarRadius(), 6356752.3142));
    CPPUNIT_ASSERT(Verify<CloseTo>(1/r.GetFlattening(), 298.257223563));
  }

  void
  testGetEllipsoid0()
  {
    typedef ReferenceEllipsoid::EllipsoidID EllipsoidID;
    // catch specific exception
    tst::Expected();
    ReferenceEllipsoid::Get(EllipsoidID(0));
    // ReferenceEllipsoid::InvalidEllipsoidException;
  }

  void
  testGetEllipsoid999()
  {
    typedef ReferenceEllipsoid::EllipsoidID EllipsoidID;
    // catch generic auger exception
    tst::Expected();
    ReferenceEllipsoid::Get(EllipsoidID(999));
    // AugerException
  }

  void
  testEllipsoidComparison()
  {
    const ReferenceEllipsoid wgs84(ReferenceEllipsoid::Get(ReferenceEllipsoid::eWGS84));
    // same parameters as in ReferenceEllipsoid code
    const ReferenceEllipsoid e(6378137., 0.00669437999013,
                               CoordinateSystem::GetRootCoordinateSystem());
    CPPUNIT_ASSERT(wgs84 == e);

    const ReferenceEllipsoid wgs84a(ReferenceEllipsoid::GetWGS84());

    CPPUNIT_ASSERT(wgs84 == wgs84a);
  }

  void
  testLatLongConversion()
  {
    const ReferenceEllipsoid wgs84(ReferenceEllipsoid::Get(ReferenceEllipsoid::eWGS84));
    const CoordinateSystemPtr ECEF = wgs84.GetECEF();

    const Point GmEq1 = Point(wgs84.GetEquatorialRadius(), 0, 0, ECEF);
    CPPUNIT_ASSERT(Verify<CloseTo>(
      wgs84.PointToLatitudeLongitudeHeight(GmEq1), Triple(0,0,0)));

    const Point GmEq2 = wgs84.LatitudeLongitudeHeightToPoint(0, 0, 0);
    CPPUNIT_ASSERT(Verify<CloseTo>(
      GmEq2.GetCoordinates(ECEF), Triple(wgs84.GetEquatorialRadius(),0,0)));

    const Point EastEq1 = Point(0, wgs84.GetEquatorialRadius(), 0, ECEF);
    CPPUNIT_ASSERT(Verify<CloseTo>(
      wgs84.PointToLatitudeLongitudeHeight(EastEq1), Triple(0,90*deg,0)));

    const Point EastEq2 = wgs84.LatitudeLongitudeHeightToPoint(0, 90*deg, 0);
    CPPUNIT_ASSERT(Verify<CloseTo>(
      EastEq2.GetCoordinates(ECEF), Triple(0, wgs84.GetEquatorialRadius(), 0)));

    const Point NPole1 = Point(0, 0, wgs84.GetPolarRadius(), ECEF);
    CPPUNIT_ASSERT(Verify<CloseTo>(
      wgs84.PointToLatitudeLongitudeHeight(NPole1), Triple(90*deg,0,0)));

    const Point NPole2 = wgs84.LatitudeLongitudeHeightToPoint(90*deg, 0, 0);
    CPPUNIT_ASSERT(Verify<CloseTo>(
      NPole2.GetCoordinates(ECEF), Triple(0,0,wgs84.GetPolarRadius())));
  }

};


CPPUNIT_TEST_SUITE_REGISTRATION(testReferenceEllipsoid);


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