PERF: Use tuned ImageRangeRegion copy in CastImageFilter

Based on performance testing, across converting between Image of
Vector and VectorImage, this loop was the best performing. Key
features of the improved performane:
 - Uses NumericTraits::GetLength over Image::GetNumberOfComponents,
 the former may be constant, while the latter is virtual
 - Uses const InputPixelType & inputPixel
 - OutputPixelType value{ outputIt.Get() } initialized to a reference
 in the output image bufffer and does not perform memory allocation
 for variable length vectors.

Author: blowekamp

Modules/Filtering/ImageFilterBase/include/itkCastImageFilter.hxx

@@ -139,29 +139,28 @@ CastImageFilter<TInputImage, TOutputImage>::DynamicThreadedGenerateDataDispatche
 
   this->CallCopyOutputRegionToInputRegion(inputRegionForThread, outputRegionForThread);
 
-  const unsigned int componentsPerPixel = outputPtr->GetNumberOfComponentsPerPixel();
-
-  // Define the iterators
-  ImageScanlineConstIterator inputIt(inputPtr, inputRegionForThread);
-  ImageScanlineIterator      outputIt(outputPtr, outputRegionForThread);
-
-  while (!inputIt.IsAtEnd())
+  auto inputRange = itk::ImageRegionRange<const TInputImage>(*inputPtr, inputRegionForThread);
+  auto outputRange = itk::ImageRegionRange<TOutputImage>(*outputPtr, outputRegionForThread);
+
+  auto       inputIt = inputRange.begin();
+  auto       outputIt = outputRange.begin();
+  const auto inputEnd = inputRange.end();
+
+  // Note: this loop has been timed for performance with converstions betwee image of vectors and VectorImages and other
+  // combinations. The following is the most efficient usage of iterators. Important considerations are the some times
+  // constant GetLength, the const & to input pixel and free from memory allocations with variable length vectors.
+  const unsigned int componentsPerPixel = itk::NumericTraits<OutputPixelType>::GetLength(*outputIt);
+  while (inputIt != inputEnd)
   {
-    while (!inputIt.IsAtEndOfLine())
+    const InputPixelType & inputPixel = *inputIt;
+    OutputPixelType        outputPixel{ *outputIt };
+    for (unsigned int k = 0; k < componentsPerPixel; ++k)
     {
-      const InputPixelType & inputPixel = inputIt.Get();
-      OutputPixelType        value{ outputIt.Get() };
-      for (unsigned int k = 0; k < componentsPerPixel; ++k)
-      {
-        value[k] = static_cast<typename OutputPixelType::ValueType>(inputPixel[k]);
-      }
-      outputIt.Set(value);
-
-      ++inputIt;
-      ++outputIt;
+      outputPixel[k] = static_cast<typename OutputPixelType::ValueType>(inputPixel[k]);
     }
-    inputIt.NextLine();
-    outputIt.NextLine();
+    *outputIt = outputPixel;
+    ++inputIt;
+    ++outputIt;
   }
 }