diff --git a/include/rtkAmsterdamShroudImageFilter.hxx b/include/rtkAmsterdamShroudImageFilter.hxx
index 14a241b50..46f12d0af 100644
--- a/include/rtkAmsterdamShroudImageFilter.hxx
+++ b/include/rtkAmsterdamShroudImageFilter.hxx
@@ -191,8 +191,9 @@ AmsterdamShroudImageFilter<TInputImage>::CropOutsideProjectedBox()
       vnl_vector<double> pCornerVnl = m_Geometry->GetMatrices()[iProj].GetVnlMatrix() * corners[ci].GetVnlVector();
       for (unsigned int i = 0; i < 2; i++)
         pCorner[i] = pCornerVnl[i] / pCornerVnl[2];
-      itk::ContinuousIndex<double, 3> pCornerI;
-      this->GetInput()->TransformPhysicalPointToContinuousIndex(pCorner, pCornerI);
+      using ValueType = typename TInputImage::PointType::ValueType;
+      const itk::ContinuousIndex<double, 3> pCornerI =
+        this->GetInput()->template TransformPhysicalPointToContinuousIndex<ValueType, double>(pCorner);
       if (ci == 0)
       {
         pCornerInf = pCornerI;
diff --git a/include/rtkCudaBackProjectionImageFilter.hxx b/include/rtkCudaBackProjectionImageFilter.hxx
index 6e2cc0de3..a61767d2f 100644
--- a/include/rtkCudaBackProjectionImageFilter.hxx
+++ b/include/rtkCudaBackProjectionImageFilter.hxx
@@ -53,8 +53,9 @@ CudaBackProjectionImageFilter<ImageType>::GPUGenerateData()
   // Rotation center (assumed to be at 0 yet)
   typename ImageType::PointType rotCenterPoint;
   rotCenterPoint.Fill(0.0);
-  itk::ContinuousIndex<double, Dimension> rotCenterIndex;
-  this->GetInput(0)->TransformPhysicalPointToContinuousIndex(rotCenterPoint, rotCenterIndex);
+  using ValueType = typename ImageType::PointType::ValueType;
+  itk::ContinuousIndex<double, Dimension> rotCenterIndex =
+    this->GetInput(0)->template TransformPhysicalPointToContinuousIndex<ValueType, double>(rotCenterPoint);
 
   // Include non-zero index in matrix
   itk::Matrix<double, 4, 4> matrixIdxVol;
diff --git a/include/rtkFDKBackProjectionImageFilter.hxx b/include/rtkFDKBackProjectionImageFilter.hxx
index 7873fbe55..42df52bdd 100644
--- a/include/rtkFDKBackProjectionImageFilter.hxx
+++ b/include/rtkFDKBackProjectionImageFilter.hxx
@@ -82,8 +82,9 @@ FDKBackProjectionImageFilter<TInputImage, TOutputImage>::DynamicThreadedGenerate
   // Rotation center (assumed to be at 0 yet)
   typename TInputImage::PointType rotCenterPoint;
   rotCenterPoint.Fill(0.0);
-  itk::ContinuousIndex<double, Dimension> rotCenterIndex;
-  this->GetInput(0)->TransformPhysicalPointToContinuousIndex(rotCenterPoint, rotCenterIndex);
+  using ValueType = typename TInputImage::PointType::ValueType;
+  const itk::ContinuousIndex<double, Dimension> rotCenterIndex =
+    this->GetInput(0)->template TransformPhysicalPointToContinuousIndex<ValueType, double>(rotCenterPoint);
 
   // Continuous index at which we interpolate
   itk::ContinuousIndex<double, Dimension - 1> pointProj;
diff --git a/include/rtkForwardWarpImageFilter.hxx b/include/rtkForwardWarpImageFilter.hxx
index 99d2b0209..cbbb896be 100644
--- a/include/rtkForwardWarpImageFilter.hxx
+++ b/include/rtkForwardWarpImageFilter.hxx
@@ -48,8 +48,9 @@ ForwardWarpImageFilter<TInputImage, TOutputImage, TDVF>::Protected_EvaluateDispl
 
   const DisplacementFieldType * fieldPtr = this->GetDisplacementField();
 
-  itk::ContinuousIndex<double, TDVF::ImageDimension> index;
-  fieldPtr->TransformPhysicalPointToContinuousIndex(point, index);
+  using ValueType = typename PointType::ValueType;
+  const itk::ContinuousIndex<double, TDVF::ImageDimension> index =
+    fieldPtr->template TransformPhysicalPointToContinuousIndex<ValueType, double>(point);
   unsigned int dim = 0; // index over dimension
   /**
    * Compute base index = closest index below point
@@ -198,8 +199,9 @@ ForwardWarpImageFilter<TInputImage, TOutputImage, TDVF>::GenerateData()
     for (unsigned int j = 0; j < TInputImage::ImageDimension; j++)
       point[j] += displacement[j];
 
-    itk::ContinuousIndex<double, TInputImage::ImageDimension> continuousIndexInOutput;
-    outputPtr->TransformPhysicalPointToContinuousIndex(point, continuousIndexInOutput);
+    using ValueType = typename TOutputImage::PointType::ValueType;
+    itk::ContinuousIndex<double, TInputImage::ImageDimension> continuousIndexInOutput =
+      outputPtr->template TransformPhysicalPointToContinuousIndex<ValueType, double>(point);
 
     // compute the base index in output, ie the closest index below point
     // Check if the baseIndex is in the output's requested region, otherwise skip the splat part
diff --git a/include/rtkZengBackProjectionImageFilter.hxx b/include/rtkZengBackProjectionImageFilter.hxx
index 50290d348..1b38ec2c1 100644
--- a/include/rtkZengBackProjectionImageFilter.hxx
+++ b/include/rtkZengBackProjectionImageFilter.hxx
@@ -191,7 +191,9 @@ ZengBackProjectionImageFilter<TInputImage, TOutputImage>::GenerateOutputInformat
                      static_cast<ContinuousIndexValueType>(sizeVolume[k] - 1) / 2.0;
   }
 
-  this->GetInput(0)->TransformContinuousIndexToPhysicalPoint(centerIndex, m_centerVolume);
+  using ValueType = typename PointType::ValueType;
+  m_centerVolume =
+    this->GetInput(0)->template TransformContinuousIndexToPhysicalPoint<CoordRepType, ValueType>(centerIndex);
 
   PointType centerRotation;
   centerRotation.Fill(0);
diff --git a/include/rtkZengForwardProjectionImageFilter.hxx b/include/rtkZengForwardProjectionImageFilter.hxx
index 90e5521ae..91fed1a16 100644
--- a/include/rtkZengForwardProjectionImageFilter.hxx
+++ b/include/rtkZengForwardProjectionImageFilter.hxx
@@ -192,7 +192,9 @@ ZengForwardProjectionImageFilter<TInputImage, TOutputImage>::GenerateOutputInfor
                      static_cast<ContinuousIndexValueType>(sizeVolume[k] - 1) / 2.0;
   }
 
-  this->GetInput(1)->TransformContinuousIndexToPhysicalPoint(centerIndex, m_centerVolume);
+  using ValueType = typename PointType::ValueType;
+  m_centerVolume =
+    this->GetInput(1)->template TransformContinuousIndexToPhysicalPoint<CoordRepType, ValueType>(centerIndex);
 
   PointType centerRotation;
   centerRotation.Fill(0);