Skip to content
Snippets Groups Projects
Commit 1d8821db authored by Julien Michel's avatar Julien Michel
Browse files

Correction tests en erreurs (suite).

parent c148a74f
No related branches found
No related tags found
No related merge requests found
...@@ -123,8 +123,13 @@ ADD_TEST(bfTuInverseLogPolarTransformNew ${BASICFILTERS_TESTS} ...@@ -123,8 +123,13 @@ ADD_TEST(bfTuInverseLogPolarTransformNew ${BASICFILTERS_TESTS}
otbInverseLogPolarTransformNew) otbInverseLogPolarTransformNew)
ADD_TEST(bfTvInverseLogPolarTransform ${BASICFILTERS_TESTS} ADD_TEST(bfTvInverseLogPolarTransform ${BASICFILTERS_TESTS}
--compare-ascii ${TOL}
${BASELINE_FILES}/bfInverseLogPolarTransformResults.txt
${TEMP}/bfInverseLogPolarTransformResults.txt
otbInverseLogPolarTransform otbInverseLogPolarTransform
1 1 6 0 0 10 20 45 60 3.14 3.14 1 0 0 1 1 1
${TEMP}/bfInverseLogPolarTransformResults.txt
6 0 0 10 20 45 60 3.14 3.14 1 0 0 1
) )
ADD_TEST(bfTvInverseLogPolarTransformResample ${BASICFILTERS_TESTS} ADD_TEST(bfTvInverseLogPolarTransformResample ${BASICFILTERS_TESTS}
......
...@@ -17,6 +17,7 @@ ...@@ -17,6 +17,7 @@
=========================================================================*/ =========================================================================*/
#include "otbInverseLogPolarTransform.h" #include "otbInverseLogPolarTransform.h"
#include "otbMacro.h" #include "otbMacro.h"
#include <fstream>
int otbInverseLogPolarTransform(int argc, char* argv[]) int otbInverseLogPolarTransform(int argc, char* argv[])
{ {
...@@ -24,20 +25,25 @@ try ...@@ -24,20 +25,25 @@ try
{ {
double radialStep = atof(argv[1]); double radialStep = atof(argv[1]);
double angularStep = atof(argv[2]); 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 double PrecisionType;
typedef otb::InverseLogPolarTransform<PrecisionType> InverseLogPolarTransformType; typedef otb::InverseLogPolarTransform<PrecisionType> InverseLogPolarTransformType;
typedef itk::Point<PrecisionType,2> PointType; typedef itk::Point<PrecisionType,2> PointType;
typedef std::vector<PointType> PointsVectorType; typedef std::vector<PointType> PointsVectorType;
std::ofstream file;
file.open(outputFilename);
file << "input points retrieval : "<<std::endl;
// input points retrieval // input points retrieval
PointsVectorType vect; PointsVectorType vect;
for(unsigned int i=0;i<nbPoints;++i) for(unsigned int i=0;i<nbPoints;++i)
{ {
PointType p; PointType p;
p[0]=atof(argv[4+2*i]); p[0]=atof(argv[5+2*i]);
p[1]=atof(argv[5+2*i]); p[1]=atof(argv[6+2*i]);
std::cout<<"Adding point "<<p<<"."<<std::endl; file <<"Adding point "<<p<<"."<<std::endl;
vect.push_back(p); vect.push_back(p);
} }
// Instantiation // Instantiation
...@@ -48,6 +54,8 @@ try ...@@ -48,6 +54,8 @@ try
params[2]=radialStep; params[2]=radialStep;
params[3]=angularStep; params[3]=angularStep;
transform->SetParameters(params); transform->SetParameters(params);
file << "Transform calculation ... :" <<std::endl;
for(PointsVectorType::iterator it=vect.begin();it!=vect.end();++it) for(PointsVectorType::iterator it=vect.begin();it!=vect.end();++it)
{ {
...@@ -74,10 +82,9 @@ try ...@@ -74,10 +82,9 @@ try
pprime[1]=0; pprime[1]=0;
} }
std::cout<<"Original Point: "<<(*it)<<", Reference point: "<<pprime<<", Transformed point: "<<p<<std::endl; file <<"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.close();
} }
catch( itk::ExceptionObject & err ) catch( itk::ExceptionObject & err )
{ {
......
...@@ -73,8 +73,6 @@ try ...@@ -73,8 +73,6 @@ try
pprime[1]=exp(logRho) * sin(theta); pprime[1]=exp(logRho) * sin(theta);
file <<"Original Point: "<<(*it)<<", Reference point: "<<pprime<<", Transformed point: "<<p<<std::endl<<std::endl; 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(); file.close();
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment