testCoordinateSystemRegistry.cc

Go to the documentation of this file.

#include <fwk/CentralConfig.h>
#include <fwk/CoordinateSystemRegistry.h>
#include <utl/CoordinateSystem.h>
#include <utl/AugerException.h>
#include <utl/ReferenceEllipsoid.h>
#include <utl/UTMPoint.h>
#include <utl/Point.h>
#include <utl/AugerUnits.h>
#include <cmath>

#include <tst/Verify.h>
#include <cppunit/extensions/HelperMacros.h>

using namespace fwk;
using namespace utl;
using namespace tst;


class CoordinateSystemRegistryTest : public CppUnit::TestFixture {

  CPPUNIT_TEST_SUITE(CoordinateSystemRegistryTest);
  CPPUNIT_TEST(testFactory);
  CPPUNIT_TEST(testPampaAmarilla);
  CPPUNIT_TEST_EXCEPTION(testInvalidCoordinateSystemException,
                         CoordinateSystemRegistry::InvalidCoordinateSystemException);
  CPPUNIT_TEST_EXCEPTION(testInitSequenceException, utl::InitSequenceException);
  CPPUNIT_TEST_SUITE_END();

public:
  void setUp() { }

  void tearDown() { }

  void
  testFactory()
  {
    CentralConfig::GetInstance(BOOTSTRAPFILE);
    CoordinateSystemRegistry::Get("ECEF");
    CoordinateSystemRegistry::Get("Malargue");
    CoordinateSystemRegistry::Get("PampaAmarilla");
    CoordinateSystemRegistry::Get("CLF");
    CoordinateSystemRegistry::Get("XLF");

    CoordinateSystemRegistry::Get("Lamar");
  }

  void
  testPampaAmarilla()
  {
    CentralConfig::GetInstance(BOOTSTRAPFILE);
    const ReferenceEllipsoid wgs84 = ReferenceEllipsoid::GetWGS84();
    const CoordinateSystemPtr pampaAmarillaCS =
      CoordinateSystemRegistry::Get("PampaAmarilla");

    // Verify the origin of the CS
    const Point pampaAmarillaOrigin(0, 0, 0, pampaAmarillaCS);
    const UTMPoint po(pampaAmarillaOrigin, wgs84);

    CPPUNIT_ASSERT(Verify<CloseTo>(po.GetNorthing(), 6099203.68*m));
    CPPUNIT_ASSERT(Verify<CloseTo>(po.GetEasting(), 477256.66*m));
    CPPUNIT_ASSERT(Verify<CloseTo>(po.GetHeight(), 1400*m));
    CPPUNIT_ASSERT(Verify<Equal>(po.GetZone(), 19));
    CPPUNIT_ASSERT(Verify<Equal>(po.GetBand(), 'H'));

    // verify coords in pampaAmarillaCS
    // fd UTM positions as used in CDAS (ints!)
    const int fdUTM[][3] = {
      { 6071872, 459208, 1416 },  // LL
      { 6094570, 498902, 1415 },  // LM
      { 6134066, 480747, 1455 },  // LA
      { 6114142, 445346, 1714 }   // CO
    };
    // corresponding CDAS XYZ result
    const double fdXYZ[][3] = {
      { -18127.757064, -27301.999633, -68.24021981845 },  // LL
      {  21647.377378,  -4690.341619, -23.48079981265 },  // LM
      {   3580.959833,  34875.225069, -41.4046727005  },  // LA
      { -31891.59545,   15028.421417, 216.51114439687 }   // CO
    };
    const double absTolerance[] = { /*X*/ 2*meter, /*Y*/ 2*meter, /*Z*/ 25*centimeter };

    for (int i = 0; i < 4; ++i) {
      const UTMPoint fdUTMPos(fdUTM[i][0]*meter, fdUTM[i][1]*meter, fdUTM[i][2]*meter, 19, 'H', wgs84);
      const Point fdPos(fdUTMPos.GetPoint());
      const double xDiff = fdPos.GetX(pampaAmarillaCS) - fdXYZ[i][0]*meter;
      const double yDiff = fdPos.GetY(pampaAmarillaCS) - fdXYZ[i][1]*meter;
      const double zDiff = fdPos.GetZ(pampaAmarillaCS) - fdXYZ[i][2]*meter;
      std::cout << i << ". diff X: " << xDiff << " Y: " << yDiff << " Z: " << zDiff << std::endl;
      CPPUNIT_ASSERT(Verify<CloseTo>(fdPos.GetX(pampaAmarillaCS), fdXYZ[i][0]*meter, absTolerance[0]));
      CPPUNIT_ASSERT(Verify<CloseTo>(fdPos.GetY(pampaAmarillaCS), fdXYZ[i][1]*meter, absTolerance[1]));
      CPPUNIT_ASSERT(Verify<CloseTo>(fdPos.GetZ(pampaAmarillaCS), fdXYZ[i][2]*meter, absTolerance[2]));
    }
  }

  void
  testInvalidCoordinateSystemException()
  {
    CentralConfig::GetInstance(BOOTSTRAPFILE);
    std::cout << "The coordinate system 'XYZ' should not exist and produce an expected error:" << std::endl;
    CoordinateSystemRegistry::Get("XYZ");
  }

  void
  testInitSequenceException()
  {
    std::cout << "Since initialization was already done, we expect an error that ReferenceEllipsoidRegistry is already initialized." << std::endl;
    CentralConfig::GetInstance(BOOTSTRAPFILE);
    CoordinateSystemRegistry::Init();
  }

};


CPPUNIT_TEST_SUITE_REGISTRATION(CoordinateSystemRegistryTest);


// Configure (x)emacs for this file ...
// Local Variables:
// mode: c++
// End: