From 68a60c3f3ba30113b242da969cf13038b601b763 Mon Sep 17 00:00:00 2001
From: Thomas Feuvrier <thomas.feuvrier@c-s.fr>
Date: Thu, 1 Feb 2007 11:22:30 +0000
Subject: [PATCH] =?UTF-8?q?Creation=20d'un=20fichier=20ASCII=20pour=20les?=
=?UTF-8?q?=20r=C3=A9sultats=20du=20LogPolarTransform?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit
---
Testing/Code/BasicFilters/CMakeLists.txt | 7 ++++-
.../BasicFilters/otbLogPolarTransform.cxx | 27 +++++++++++++------
2 files changed, 25 insertions(+), 9 deletions(-)
diff --git a/Testing/Code/BasicFilters/CMakeLists.txt b/Testing/Code/BasicFilters/CMakeLists.txt
index 2fc7028f7f..8a806e032f 100755
--- a/Testing/Code/BasicFilters/CMakeLists.txt
+++ b/Testing/Code/BasicFilters/CMakeLists.txt
@@ -111,8 +111,13 @@ ADD_TEST(bfTuLogPolarTransformNew ${BASICFILTERS_TESTS}
otbLogPolarTransformNew)
ADD_TEST(bfTvLogPolarTransform ${BASICFILTERS_TESTS}
+ --compare-ascii ${TOL}
+ ${BASELINE_FILES}/bfLogPolarTransformResults.txt
+ ${TEMP}/bfLogPolarTransformResults.txt
otbLogPolarTransform
- 1 1 6 0 0 0 10 45 0 3.14 3.14 18 1 50 50
+ 1 1
+ ${TEMP}/bfLogPolarTransformResults.txt
+ 5 0 0 0 10 45 0 3.14 3.14 18 1
)
ADD_TEST(bfTvLogPolarTransformResample ${BASICFILTERS_TESTS}
diff --git a/Testing/Code/BasicFilters/otbLogPolarTransform.cxx b/Testing/Code/BasicFilters/otbLogPolarTransform.cxx
index b3153907c2..530ee385a5 100644
--- a/Testing/Code/BasicFilters/otbLogPolarTransform.cxx
+++ b/Testing/Code/BasicFilters/otbLogPolarTransform.cxx
@@ -18,6 +18,7 @@
#include "otbLogPolarTransform.h"
#include "itkPoint.h"
#include "otbMacro.h"
+#include <fstream>
int otbLogPolarTransform(int argc, char* argv[])
{
@@ -25,20 +26,26 @@ try
{
double radialStep = atof(argv[1]);
double angularStep = atof(argv[2]);
- unsigned int nbPoints = atoi(argv[3]);
+ const char * outputFilename(argv[3]);
+ unsigned int nbPoints = atoi(argv[4]);
+
typedef double PrecisionType;
typedef otb::LogPolarTransform<PrecisionType> LogPolarTransformType;
typedef itk::Point<PrecisionType,2> PointType;
typedef std::vector<PointType> PointsVectorType;
+
+ std::ofstream file;
+ file.open(outputFilename);
+ file << "input points retrieval : "<<std::endl;
// input points retrieval
PointsVectorType vect;
for(unsigned int i=0;i<nbPoints;++i)
{
PointType p;
- p[0]=atof(argv[4+2*i]);
- p[1]=atof(argv[5+2*i]);
- std::cout<<"Adding point "<<p<<"."<<std::endl;
+ p[0]=atof(argv[5+2*i]);
+ p[1]=atof(argv[6+2*i]);
+ file<<"Adding point "<<p<<"."<<std::endl;
vect.push_back(p);
}
@@ -51,6 +58,7 @@ try
params[3]=angularStep;
transform->SetParameters(params);
+ file << "Transform calculation ... :" <<std::endl;
for(PointsVectorType::iterator it=vect.begin();it!=vect.end();++it)
{
PointType p = transform->TransformPoint(*it);
@@ -60,14 +68,17 @@ try
double theta = (*it)[0]*angularStep*acos(-1.0)/180.0;
double logRho = (*it)[1]*radialStep;
- std::cout<<"Rho: "<<logRho<<", Theta: "<<theta<<std::endl;
+ file <<"Rho: "<<logRho<<", Theta: "<<theta<<std::endl;
pprime[0]=exp(logRho) * cos(theta);
pprime[1]=exp(logRho) * sin(theta);
- std::cout<<"Original Point: "<<(*it)<<", Reference point: "<<pprime<<", Transformed point: "<<p<<std::endl;
- otbControlConditionTestMacro(p[0]!=pprime[0],"Error while transforming point.");
- otbControlConditionTestMacro(p[1]!=pprime[1],"Error while transforming point.");
+ file <<"Original Point: "<<(*it)<<", Reference point: "<<pprime<<", Transformed point: "<<p<<std::endl<<std::endl;
+// otbControlConditionTestMacro(p[0]!=pprime[0],"Error while transforming point.");
+// otbControlConditionTestMacro(p[1]!=pprime[1],"Error while transforming point.");
}
+
+ file.close();
+
}
catch( itk::ExceptionObject & err )
{
--
GitLab