From 4afe3ea3a8399d133e605ab5c9cbefdaa16c2d07 Mon Sep 17 00:00:00 2001
From: Emmanuel Christophe <emmanuel.christophe@orfeo-toolbox.org>
Date: Sat, 26 Mar 2011 20:41:44 -0700
Subject: [PATCH] TEST: properly exercise the transforms
---
Testing/Code/Projections/otbMapProjection.cxx | 64 ++++++++++++++++---
1 file changed, 56 insertions(+), 8 deletions(-)
diff --git a/Testing/Code/Projections/otbMapProjection.cxx b/Testing/Code/Projections/otbMapProjection.cxx
index 2e21b5aacb..b4f29e7c51 100644
--- a/Testing/Code/Projections/otbMapProjection.cxx
+++ b/Testing/Code/Projections/otbMapProjection.cxx
@@ -24,32 +24,80 @@ int otbMapProjection(int argc, char* argv[])
{
const char * outFileName = argv[1];
std::ofstream file;
+
+ itk::Point<double, 2> point;
+ point[0] = 1.44;
+ point[1] = 43.6050;
+ itk::Point<double, 2> point2;
file.open(outFileName);
+ file << std::setprecision(15);
/** Test the output of the projection in Wkt format*/
+ otb::UtmForwardProjection::Pointer lUtmProjection2 = otb::UtmForwardProjection::New();
+ point2 = lUtmProjection2->TransformPoint(point);
+ file << lUtmProjection2->GetWkt() << std::endl << std::endl;
+ file << "Forward projection: " << std::endl;
+ file << point << " -> ";
+ file << point2;
+ file << std::endl << std::endl;
+
+
otb::UtmInverseProjection::Pointer lUtmProjection = otb::UtmInverseProjection::New();
file << lUtmProjection->GetWkt() << std::endl << std::endl;
+ file << "Inverse projection: " << std::endl;
+ file << point2 << " -> ";
+ file << lUtmProjection->TransformPoint(point2);
+ file << std::endl << std::endl;
- otb::UtmForwardProjection::Pointer lUtmProjection2 = otb::UtmForwardProjection::New();
- file << lUtmProjection2->GetWkt() << std::endl << std::endl;
+
+ otb::Lambert2EtenduForwardProjection::Pointer lLambert2Etendu2 = otb::Lambert2EtenduForwardProjection::New();
+ point2 = lLambert2Etendu2->TransformPoint(point);
+ file << lLambert2Etendu2->GetWkt() << std::endl << std::endl;
+ file << "Forward projection: " << std::endl;
+ file << point << " -> ";
+ file << point2;
+ file << std::endl << std::endl;
otb::Lambert2EtenduInverseProjection::Pointer lLambert2Etendu = otb::Lambert2EtenduInverseProjection::New();
file << lLambert2Etendu->GetWkt() << std::endl << std::endl;
+ file << "Inverse projection: " << std::endl;
+ file << point2 << " -> ";
+ file << lLambert2Etendu->TransformPoint(point2);
+ file << std::endl << std::endl;
- otb::Lambert2EtenduForwardProjection::Pointer lLambert2Etendu2 = otb::Lambert2EtenduForwardProjection::New();
- file << lLambert2Etendu2->GetWkt() << std::endl << std::endl;
+
+ otb::Lambert93ForwardProjection::Pointer lLambert93_2 = otb::Lambert93ForwardProjection::New();
+ point2 = lLambert93_2->TransformPoint(point);
+ file << lLambert93_2->GetWkt() << std::endl << std::endl;
+ file << "Forward projection: " << std::endl;
+ file << point << " -> ";
+ file << point2;
+ file << std::endl << std::endl;
otb::Lambert93InverseProjection::Pointer lLambert93 = otb::Lambert93InverseProjection::New();
file << lLambert93->GetWkt() << std::endl << std::endl;
+ file << "Inverse projection: " << std::endl;
+ file << point2 << " -> ";
+ file << lLambert93->TransformPoint(point2);
+ file << std::endl << std::endl;
- otb::Lambert93ForwardProjection::Pointer lLambert93_2 = otb::Lambert93ForwardProjection::New();
- file << lLambert93_2->GetWkt() << std::endl << std::endl;
- otb::MercatorInverseProjection::Pointer lMercatorProjection = otb::MercatorInverseProjection::New();
- file << lMercatorProjection->GetWkt() << std::endl << std::endl;
otb::MercatorForwardProjection::Pointer lMercatorProjection2 = otb::MercatorForwardProjection::New();
+ point2 = lMercatorProjection2->TransformPoint(point);
file << lMercatorProjection2->GetWkt() << std::endl << std::endl;
+ file << "Forward projection: " << std::endl;
+ file << point << " -> ";
+ file << point2;
+ file << std::endl << std::endl;
+
+
+ otb::MercatorInverseProjection::Pointer lMercatorProjection = otb::MercatorInverseProjection::New();
+ file << lMercatorProjection->GetWkt() << std::endl << std::endl;
+ file << "Inverse projection: " << std::endl;
+ file << point2 << " -> ";
+ file << lMercatorProjection->TransformPoint(point2);
+ file << std::endl << std::endl;
// Moving to GenericMapProjection, change the access
--
GitLab