ENH: better templates in Transforms

Code/IO/otbPointSetFileReader.txx

@@ -172,7 +172,6 @@ void PointSetFileReader<TOutputPointSet>
       point[1] = p.GetY();
       point[2] = p.GetZ();
 
-
       unsigned long i = output->GetNumberOfPoints();
       output->SetPoint( i, point );
 

Code/Projections/otbCompositeTransform.h

@@ -39,15 +39,15 @@ namespace otb
 * \ingroup Transform
 */
 
-typedef enum
-{
-  PROJDEFAULT=0,
-  PROJIDENTITY=1,
-  PROJMAPFORWARD=2,
-  PROJMAPINVERSE=3,
-  PROJSENSORFORWARD=4,
-  PROJSENSORINVERSE=5
-} ProjectionTypeEnum;
+//typedef enum
+//{
+//  PROJDEFAULT=0,
+//  PROJIDENTITY=1,
+//  PROJMAPFORWARD=2,
+//  PROJMAPINVERSE=3,
+//  PROJSENSORFORWARD=4,
+//  PROJSENSORINVERSE=5
+//} ProjectionTypeEnum;
 
 template <class TFirstTransform,
 class TSecondTransform,

Code/Projections/otbCompositeTransform.txx

@@ -82,7 +82,7 @@ NOutputDimensions>
   SecondTransformOutputPointType outputPoint;
   outputPoint=m_SecondTransform->TransformPoint(geoPoint);
 
-//       otbMsgDevMacro(<< "Converting: " << point1 << " -> " <<  geoPoint<< " -> " << outputPoint);
+//  otbMsgDevMacro(<< std::setprecision(15) << "Converting: " << point1 << " -> " <<  geoPoint<< " -> " << outputPoint);
 
   return outputPoint;
 }

Code/Projections/otbGenericRSTransform.txx

@@ -265,6 +265,8 @@ GenericRSTransform<TScalarType, NInputDimensions, NOutputDimensions>
   outputPoint[0] = (outputPoint[0] - m_OutputOrigin[0]) / m_OutputSpacing[0];
   outputPoint[1] = (outputPoint[1] - m_OutputOrigin[1]) / m_OutputSpacing[1];
 
+//  otbMsgDevMacro("GenericRSTransform: " << point << " -> " << outputPoint);
+
   return outputPoint;
 }
 

Code/Projections/otbInverseSensorModel.txx

@@ -86,7 +86,10 @@ InverseSensorModel< TScalarType, NInputDimensions, NOutputDimensions>
 
   outputPoint[0]=ossimDPoint.x;
   outputPoint[1]=ossimDPoint.y;
-
+  if (OutputPointType::PointDimension == 3)
+    {
+      outputPoint[2] = ossimGPoint.height();
+    }
 //     otbMsgDevMacro(<< "Point in sensor geometry: (" << outputPoint[0] << "," <<  outputPoint[1] << ")");
 
   return outputPoint;

Code/Projections/otbSensorModelBase.txx

@@ -37,7 +37,7 @@ template < class TScalarType,
 unsigned int NInputDimensions,
 unsigned int NOutputDimensions >
 SensorModelBase< TScalarType,NInputDimensions,NOutputDimensions>
-::SensorModelBase()
+::SensorModelBase(): Superclass(OutputSpaceDimension, 0)
 {
   m_Model = NULL;
   m_DEMHandler = DEMHandlerType::New();