Skip to content
Snippets Groups Projects
Commit a3da67be authored by Guillaume Pasero's avatar Guillaume Pasero
Browse files

ENH: use new DEMHandler in test prTvSensorModel

parent 2e9e63e4
Branches
Tags
No related merge requests found
......@@ -23,6 +23,7 @@
#include "otbImageFileReader.h"
#include "otbForwardSensorModel.h"
#include "otbInverseSensorModel.h"
#include "otbDEMHandler.h"
int otbSensorModel(int argc, char* argv[])
{
......@@ -51,6 +52,8 @@ int otbSensorModel(int argc, char* argv[])
file << reader->GetOutput()->GetImageKeywordlist();
file << "\n*** TRANSFORM ***\n";
otb::DEMHandler::Instance()->SetDefaultHeightAboveEllipsoid(16.19688987731934);
typedef otb::ForwardSensorModel<double> ForwardSensorModelType;
ForwardSensorModelType::Pointer forwardSensorModel = ForwardSensorModelType::New();
......@@ -60,7 +63,6 @@ int otbSensorModel(int argc, char* argv[])
std::cout<<"Invalid Model pointer m_Model == NULL!\n The ossim keywordlist is invalid!"<<std::endl;
return EXIT_FAILURE;
}
forwardSensorModel->SetAverageElevation(16.19688987731934);
itk::Point<double, 2> imagePoint;
imagePoint[0] = 10;
......@@ -80,7 +82,6 @@ int otbSensorModel(int argc, char* argv[])
std::cout<<"Invalid Model pointer m_Model == NULL!\n The ossim keywordlist is invalid!"<<std::endl;
return EXIT_FAILURE;
}
inverseSensorModel->SetAverageElevation(16.19688987731934);
itk::Point<double, 2> reversedImagePoint;
reversedImagePoint = inverseSensorModel->TransformPoint(geoPoint);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment