Skip to content
Snippets Groups Projects
Commit b3111d5b authored by Emmanuel Christophe's avatar Emmanuel Christophe
Browse files

COMP: remove switch case

parent 40b19cfd
No related branches found
No related tags found
No related merge requests found
...@@ -148,55 +148,46 @@ GenericMapProjection<Transform, TScalarType, NInputDimensions, NOutputDimensions ...@@ -148,55 +148,46 @@ GenericMapProjection<Transform, TScalarType, NInputDimensions, NOutputDimensions
{ {
OutputPointType outputPoint; OutputPointType outputPoint;
// otbMsgDevMacro(<< "DirectionOfMapping: " << DirectionOfMapping); // otbMsgDevMacro(<< "DirectionOfMapping: " <<
switch (DirectionOfMapping) // DirectionOfMapping);
if (DirectionOfMapping == INVERSE)
{ {
case INVERSE:
{
// otbMsgDevMacro(<< "Cartographic coordinates: (" << point[0] << "," << point[1] << ")"); // otbMsgDevMacro(<< "Cartographic coordinates: (" << point[0] << "," << point[1] << ")");
//from "itk::point" to "ossim::ossimDpt" //from "itk::point" to "ossim::ossimDpt"
ossimDpt ossimDPoint(point[0], point[1]); ossimDpt ossimDPoint(point[0], point[1]);
//map projection //map projection
ossimGpt ossimGPoint; ossimGpt ossimGPoint;
ossimGPoint = this->GetMapProjection()->inverse(ossimDPoint); ossimGPoint = this->GetMapProjection()->inverse(ossimDPoint);
ossimGPoint.changeDatum(ossimDatumFactory::instance()->wgs84()); ossimGPoint.changeDatum(ossimDatumFactory::instance()->wgs84());
// otbGenericMsgDebugMacro(<< "Inverse : " << std::endl << m_MapProjection->print(std::cout)); // otbGenericMsgDebugMacro(<< "Inverse : " << std::endl << m_MapProjection->print(std::cout));
outputPoint[0] = ossimGPoint.lon; outputPoint[0] = ossimGPoint.lon;
outputPoint[1] = ossimGPoint.lat; outputPoint[1] = ossimGPoint.lat;
// otbMsgDevMacro(<< "Geographic coordinates (lon, lat) : (" << outputPoint[0] << "," << outputPoint[1] << ")"); // otbMsgDevMacro(<< "Geographic coordinates (lon, lat) : (" << outputPoint[0] << "," << outputPoint[1] << ")");
if ((InputPointType::PointDimension == 3) && (OutputPointType::PointDimension == 3)) if ((InputPointType::PointDimension == 3) && (OutputPointType::PointDimension == 3))
{
outputPoint[2] = point[2];
}
break;
}
case FORWARD:
{ {
outputPoint[2] = point[2];
}
}
if (DirectionOfMapping == FORWARD)
{
// otbMsgDevMacro(<< "Geographic coordinates (lon, lat) : (" << point[1] << "," << point[0] << ")"); // otbMsgDevMacro(<< "Geographic coordinates (lon, lat) : (" << point[1] << "," << point[0] << ")");
//from "itk::point" to "ossim::ossimGpt" //from "itk::point" to "ossim::ossimGpt"
ossimGpt ossimGPoint(point[1], point[0]); ossimGpt ossimGPoint(point[1], point[0]);
//map projection //map projection
ossimDpt ossimDPoint; ossimDpt ossimDPoint;
ossimDPoint = this->GetMapProjection()->forward(ossimGPoint); ossimDPoint = this->GetMapProjection()->forward(ossimGPoint);
// otbGenericMsgDebugMacro(<< "Forward : ========================= " << std::endl << m_MapProjection->print(std::cout)); // otbGenericMsgDebugMacro(<< "Forward : ========================= " << std::endl << m_MapProjection->print(std::cout));
outputPoint[0] = ossimDPoint.x; outputPoint[0] = ossimDPoint.x;
outputPoint[1] = ossimDPoint.y; outputPoint[1] = ossimDPoint.y;
// otbMsgDevMacro(<< "Cartographic coordinates: (" << outputPoint[0] << "," << outputPoint[1] << ")"); // otbMsgDevMacro(<< "Cartographic coordinates: (" << outputPoint[0] << "," << outputPoint[1] << ")");
if ((InputPointType::PointDimension == 3) && (OutputPointType::PointDimension == 3)) if ((InputPointType::PointDimension == 3) && (OutputPointType::PointDimension == 3))
{
outputPoint[2] = point[2];
}
break;
}
default:
{ {
itkExceptionMacro(<< "Model is INVERSE or FORWARD only !!"); outputPoint[2] = point[2];
break;
} }
} }
......
...@@ -117,48 +117,37 @@ MapProjection<TOssimMapProjection, Transform, TScalarType, NInputDimensions, NOu ...@@ -117,48 +117,37 @@ MapProjection<TOssimMapProjection, Transform, TScalarType, NInputDimensions, NOu
{ {
OutputPointType outputPoint; OutputPointType outputPoint;
switch (DirectionOfMapping) if (DirectionOfMapping == INVERSE)
{ {
case INVERSE:
{
// otbMsgDevMacro(<< "Cartographic coordinates: (" << point[0] << "," << point[1] << ")"); // otbMsgDevMacro(<< "Cartographic coordinates: (" << point[0] << "," << point[1] << ")");
//from "itk::point" to "ossim::ossimDpt" //from "itk::point" to "ossim::ossimDpt"
ossimDpt ossimDPoint(point[0], point[1]); ossimDpt ossimDPoint(point[0], point[1]);
//map projection //map projection
ossimGpt ossimGPoint; ossimGpt ossimGPoint;
ossimGPoint = m_MapProjection->inverse(ossimDPoint); ossimGPoint = m_MapProjection->inverse(ossimDPoint);
ossimGPoint.changeDatum(ossimDatumFactory::instance()->wgs84()); ossimGPoint.changeDatum(ossimDatumFactory::instance()->wgs84());
// otbGenericMsgDebugMacro(<< "Inverse : " << std::endl << m_MapProjection->print(std::cout)); // otbGenericMsgDebugMacro(<< "Inverse : " << std::endl << m_MapProjection->print(std::cout));
outputPoint[0] = ossimGPoint.lon; outputPoint[0] = ossimGPoint.lon;
outputPoint[1] = ossimGPoint.lat; outputPoint[1] = ossimGPoint.lat;
// otbMsgDevMacro(<< "Geographic coordinates (lon, lat) : (" << outputPoint[0] << "," << outputPoint[1] << ")"); // otbMsgDevMacro(<< "Geographic coordinates (lon, lat) : (" << outputPoint[0] << "," << outputPoint[1] << ")");
break; }
} if (DirectionOfMapping == FORWARD)
case FORWARD: {
{
// otbMsgDevMacro(<< "Geographic coordinates (lon, lat) : (" << point[1] << "," << point[0] << ")"); // otbMsgDevMacro(<< "Geographic coordinates (lon, lat) : (" << point[1] << "," << point[0] << ")");
//from "itk::point" to "ossim::ossimGpt" //from "itk::point" to "ossim::ossimGpt"
ossimGpt ossimGPoint(point[1], point[0]); ossimGpt ossimGPoint(point[1], point[0]);
//map projection //map projection
ossimDpt ossimDPoint; ossimDpt ossimDPoint;
ossimDPoint = m_MapProjection->forward(ossimGPoint); ossimDPoint = m_MapProjection->forward(ossimGPoint);
// otbGenericMsgDebugMacro(<< "Forward : ========================= " << std::endl << m_MapProjection->print(std::cout)); // otbGenericMsgDebugMacro(<< "Forward : ========================= " << std::endl << m_MapProjection->print(std::cout));
outputPoint[0] = ossimDPoint.x; outputPoint[0] = ossimDPoint.x;
outputPoint[1] = ossimDPoint.y; outputPoint[1] = ossimDPoint.y;
// otbMsgDevMacro(<< "Cartographic coordinates: (" << outputPoint[0] << "," << outputPoint[1] << ")"); // otbMsgDevMacro(<< "Cartographic coordinates: (" << outputPoint[0] << "," << outputPoint[1] << ")");
break;
}
default:
{
itkExceptionMacro(<< "Model is INVERSE or FORWARD only !!");
break;
}
} }
return outputPoint; return outputPoint;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment