Merge pull request InsightSoftwareConsortium#5992 from N-Dekker/Remove-GoToBegin-on-ImageIterators

hjmjohnson · web-flow · commit 2b7e6d4f1d54 · 2026-03-31T13:38:19.000-05:00
Remove GoToBegin() calls from Image Iterators in hxx files
diff --git a/Modules/Filtering/BiasCorrection/include/itkMRASlabIdentifier.hxx b/Modules/Filtering/BiasCorrection/include/itkMRASlabIdentifier.hxx
@@ -78,7 +78,6 @@ MRASlabIdentifier<TInputImage>::GenerateSlabRegions()
     region.SetIndex(index);
 
     ImageRegionConstIterator iter(m_Image, region);
-    iter.GoToBegin();
 
     std::priority_queue<ImagePixelType> mins;
     for (unsigned int i = 0; i < m_NumberOfSamples; ++i)
diff --git a/Modules/Filtering/BiasCorrection/include/itkMRIBiasFieldCorrectionFilter.hxx b/Modules/Filtering/BiasCorrection/include/itkMRIBiasFieldCorrectionFilter.hxx
@@ -85,7 +85,6 @@ MRIBiasEnergyFunction<TImage, TImageMask, TBiasField>::GetValue(const Parameters
     typename TBiasField::SimpleForwardIterator bIter(m_BiasField);
 
     bIter.Begin();
-    iIter.GoToBegin();
     total = 0.0;
 
     // Fastest for full sampling
@@ -102,7 +101,6 @@ MRIBiasEnergyFunction<TImage, TImageMask, TBiasField>::GetValue(const Parameters
     else
     {
       ImageRegionIterator mIter(m_Mask, m_Region);
-      mIter.GoToBegin();
       while (!bIter.IsAtEnd())
       {
         if (mIter.Get() > 0.0)
@@ -547,7 +545,6 @@ MRIBiasFieldCorrectionFilter<TInputImage, TOutputImage, TMaskImage>::CorrectImag
   {
     itkDebugMacro("Output mask is being used");
     ImageRegionIterator mIter(m_OutputMask, region);
-    mIter.GoToBegin();
     while (!bIter.IsAtEnd())
     {
       const double inputPixel = iIter.Get();
diff --git a/Modules/Filtering/DistanceMap/include/itkSignedMaurerDistanceMapImageFilter.hxx b/Modules/Filtering/DistanceMap/include/itkSignedMaurerDistanceMapImageFilter.hxx
@@ -253,9 +253,6 @@ SignedMaurerDistanceMapImageFilter<TInputImage, TOutputImage>::ThreadedGenerateD
     ImageRegionIterator      Ot(outputPtr, outputRegion);
     ImageRegionConstIterator It(m_InputCache, outputRegion);
 
-    Ot.GoToBegin();
-    It.GoToBegin();
-
     ProgressReporter progress2(this,
                                threadId,
                                outputRegionForThread.GetNumberOfPixels(),
diff --git a/Modules/Filtering/ImageGrid/include/itkBSplineScatteredDataPointSetToImageFilter.hxx b/Modules/Filtering/ImageGrid/include/itkBSplineScatteredDataPointSetToImageFilter.hxx
@@ -658,8 +658,6 @@ BSplineScatteredDataPointSetToImageFilter<TInputPointSet, TOutputImage>::AfterTh
 
       ItD.GoToBegin();
       ItO.GoToBegin();
-      Itd.GoToBegin();
-      Ito.GoToBegin();
       while (!ItD.IsAtEnd())
       {
         ItD.Set(ItD.Get() + Itd.Get());
diff --git a/Modules/Filtering/ImageIntensity/include/itkDiscreteGaussianDerivativeImageFunction.hxx b/Modules/Filtering/ImageIntensity/include/itkDiscreteGaussianDerivativeImageFunction.hxx
@@ -136,8 +136,7 @@ DiscreteGaussianDerivativeImageFunction<TInputImage, TOutput>::RecomputeGaussian
 
   // Copy kernel image to neighborhood. Do not copy boundaries.
   ImageRegionConstIterator it(kernelImage, kernelRegion);
-  it.GoToBegin();
-  unsigned int idx = 0;
+  unsigned int             idx = 0;
 
   while (!it.IsAtEnd())
   {
diff --git a/Modules/Filtering/ImageIntensity/include/itkDiscreteGradientMagnitudeGaussianImageFunction.hxx b/Modules/Filtering/ImageIntensity/include/itkDiscreteGradientMagnitudeGaussianImageFunction.hxx
@@ -158,8 +158,7 @@ DiscreteGradientMagnitudeGaussianImageFunction<TInputImage, TOutput>::RecomputeG
 
     // Copy kernel image to neighborhood. Do not copy boundaries.
     ImageRegionConstIterator it(kernelImage, kernelRegion);
-    it.GoToBegin();
-    unsigned int idx = 0;
+    unsigned int             idx = 0;
 
     while (!it.IsAtEnd())
     {
diff --git a/Modules/Filtering/ImageIntensity/include/itkDiscreteHessianGaussianImageFunction.hxx b/Modules/Filtering/ImageIntensity/include/itkDiscreteHessianGaussianImageFunction.hxx
@@ -164,8 +164,7 @@ DiscreteHessianGaussianImageFunction<TInputImage, TOutput>::RecomputeGaussianKer
 
       // Copy kernel image to neighborhood. Do not copy boundaries.
       ImageRegionConstIterator it(kernelImage, kernelRegion);
-      it.GoToBegin();
-      unsigned int idx = 0;
+      unsigned int             idx = 0;
       while (!it.IsAtEnd())
       {
         m_KernelArray[kernelidx][idx] = it.Get();
diff --git a/Modules/Filtering/ImageIntensity/include/itkHistogramMatchingImageFilter.hxx b/Modules/Filtering/ImageIntensity/include/itkHistogramMatchingImageFilter.hxx
@@ -437,8 +437,6 @@ HistogramMatchingImageFilter<TInputImage, TOutputImage, THistogramMeasurement>::
   {
     // put each image pixel into the histogram
     ImageRegionConstIterator iter(image, image->GetBufferedRegion());
-
-    iter.GoToBegin();
     while (!iter.IsAtEnd())
     {
       const InputPixelType & value = iter.Value();
diff --git a/Modules/Filtering/ImageStatistics/include/itkAccumulateImageFilter.hxx b/Modules/Filtering/ImageStatistics/include/itkAccumulateImageFilter.hxx
@@ -181,8 +181,7 @@ AccumulateImageFilter<TInputImage, TOutputImage>::GenerateData()
     }
     const typename TInputImage::RegionType AccumulatedRegion(AccumulatedIndex, AccumulatedSize);
     ImageRegionConstIterator               inputIter(inputImage, AccumulatedRegion);
-    inputIter.GoToBegin();
-    AccumulateType Value{};
+    AccumulateType                         Value{};
     while (!inputIter.IsAtEnd())
     {
       Value += static_cast<AccumulateType>(inputIter.Get());
diff --git a/Modules/Filtering/ImageStatistics/include/itkImagePCAShapeModelEstimator.hxx b/Modules/Filtering/ImageStatistics/include/itkImagePCAShapeModelEstimator.hxx
@@ -149,7 +149,6 @@ ImagePCAShapeModelEstimator<TInputImage, TOutputImage>::GenerateData()
     ImageRegionIterator outIterJ(this->GetOutput(j), region);
 
     unsigned int idx = 0;
-    outIterJ.GoToBegin();
     while (!outIterJ.IsAtEnd())
     {
       outIterJ.Set(static_cast<typename OutputImageType::PixelType>(m_OneEigenVector[idx]));
@@ -166,8 +165,6 @@ ImagePCAShapeModelEstimator<TInputImage, TOutputImage>::GenerateData()
   {
     region = this->GetOutput(j)->GetRequestedRegion();
     ImageRegionIterator outIterJ(this->GetOutput(j), region);
-
-    outIterJ.GoToBegin();
     while (!outIterJ.IsAtEnd())
     {
       outIterJ.Set(0);
diff --git a/Modules/Filtering/LabelMap/include/itkLabelMapToBinaryImageFilter.hxx b/Modules/Filtering/LabelMap/include/itkLabelMapToBinaryImageFilter.hxx
@@ -101,9 +101,6 @@ LabelMapToBinaryImageFilter<TInputImage, TOutputImage>::DynamicThreadedGenerateD
     ImageRegionConstIterator bgIt(this->GetBackgroundImage(), outputRegionForThread);
     ImageRegionIterator      oIt(output, outputRegionForThread);
 
-    bgIt.GoToBegin();
-    oIt.GoToBegin();
-
     while (!oIt.IsAtEnd())
     {
       const OutputImagePixelType & bg = bgIt.Get();
@@ -123,7 +120,6 @@ LabelMapToBinaryImageFilter<TInputImage, TOutputImage>::DynamicThreadedGenerateD
   {
     // fill the background with the background value
     ImageRegionIterator oIt(output, outputRegionForThread);
-    oIt.GoToBegin();
 
     while (!oIt.IsAtEnd())
     {
diff --git a/Modules/Filtering/Smoothing/include/itkBinomialBlurImageFilter.hxx b/Modules/Filtering/Smoothing/include/itkBinomialBlurImageFilter.hxx
@@ -184,8 +184,6 @@ BinomialBlurImageFilter<TInputImage, TOutputImage>::GenerateData()
       //----------------------Reverse pass----------------------
       ImageRegionReverseIterator tempReverseIt(tempPtr, tempPtr->GetRequestedRegion());
 
-      tempReverseIt.GoToBegin();
-
       while (!tempReverseIt.IsAtEnd())
       {
         // determine the index of the output pixel
diff --git a/Modules/Nonunit/Review/include/itkAttributeMorphologyBaseImageFilter.hxx b/Modules/Nonunit/Review/include/itkAttributeMorphologyBaseImageFilter.hxx
@@ -173,7 +173,6 @@ AttributeMorphologyBaseImageFilter<TInputImage, TOutputImage, TAttribute, TFunct
   // resolving phase
   // copy pixels back
   ImageRegionIterator ORegIt(output, output->GetRequestedRegion());
-  ORegIt.GoToBegin();
 
   // fill Raw - worry about iteration details later.
   // We aren't filling m_Parent, as suggested in the paper, because it
diff --git a/Modules/Nonunit/Review/include/itkMultiphaseDenseFiniteDifferenceImageFilter.hxx b/Modules/Nonunit/Review/include/itkMultiphaseDenseFiniteDifferenceImageFilter.hxx
@@ -64,9 +64,6 @@ MultiphaseDenseFiniteDifferenceImageFilter<TInputImage, TFeatureImage, TOutputIm
     // Fill the output pointer
     auto p = static_cast<OutputPixelType>(this->m_Lookup[i]);
 
-    in.GoToBegin();
-    out.GoToBegin();
-
     while (!out.IsAtEnd())
     {
       if (in.Get() < 0)
@@ -132,7 +129,6 @@ typename MultiphaseDenseFiniteDifferenceImageFilter<TInputImage, TFeatureImage,
       ImageRegionIterator      nU(m_UpdateBuffers[i], *fIt);
 
       nD.GoToBegin();
-      nU.GoToBegin();
 
       while (!nD.IsAtEnd())
       {
@@ -206,9 +202,6 @@ MultiphaseDenseFiniteDifferenceImageFilter<TInputImage, TFeatureImage, TOutputIm
     ImageRegionIterator u(m_UpdateBuffers[i], region);
     ImageRegionIterator o(this->m_LevelSet[i], region);
 
-    u.GoToBegin();
-    o.GoToBegin();
-
     while (!u.IsAtEnd())
     {
       InputPixelType val = static_cast<InputPixelType>(dt) * u.Get();
@@ -240,7 +233,6 @@ MultiphaseDenseFiniteDifferenceImageFilter<TInputImage, TFeatureImage, TOutputIm
       rms_change_accumulator = 0;
 
       o.GoToBegin();
-      it.GoToBegin();
 
       while (!o.IsAtEnd())
       {
diff --git a/Modules/Nonunit/Review/include/itkMultiphaseSparseFiniteDifferenceImageFilter.hxx b/Modules/Nonunit/Review/include/itkMultiphaseSparseFiniteDifferenceImageFilter.hxx
@@ -138,9 +138,6 @@ MultiphaseSparseFiniteDifferenceImageFilter<TInputImage, TFeatureImage, TOutputI
     ImageRegionIterator lIt(input, region);
     ImageRegionIterator tIt(tempImage, region);
 
-    lIt.GoToBegin();
-    tIt.GoToBegin();
-
     while (!lIt.IsAtEnd())
     {
       tIt.Set(lIt.Get());
@@ -160,7 +157,6 @@ MultiphaseSparseFiniteDifferenceImageFilter<TInputImage, TFeatureImage, TOutputI
     ImageRegionIterator zIt(zeroCrossingFilter->GetOutput(), region);
 
     lIt.GoToBegin();
-    zIt.GoToBegin();
 
     while (!lIt.IsAtEnd())
     {
@@ -1015,8 +1011,6 @@ MultiphaseSparseFiniteDifferenceImageFilter<TInputImage, TFeatureImage, TOutputI
       // For each region, set the pixel in m_StatusImage to
       // m_StatusBoundaryPixel
       ImageRegionIterator statusIt(sparsePtr->m_StatusImage, *fit);
-
-      statusIt.GoToBegin();
       while (!statusIt.IsAtEnd())
       {
         statusIt.Set(m_StatusBoundaryPixel);
@@ -1124,9 +1118,6 @@ MultiphaseSparseFiniteDifferenceImageFilter<TInputImage, TFeatureImage, TOutputI
 
     ImageRegionIterator outputIt(this->m_LevelSet[fId], this->m_LevelSet[fId]->GetRequestedRegion());
 
-    outputIt.GoToBegin();
-    statusIt.GoToBegin();
-
     while (!outputIt.IsAtEnd())
     {
       if (statusIt.Get() == m_StatusNull || statusIt.Get() == m_StatusBoundaryPixel)
@@ -1333,9 +1324,6 @@ MultiphaseSparseFiniteDifferenceImageFilter<TInputImage, TFeatureImage, TOutputI
 
     auto p = static_cast<OutputPixelType>(this->m_Lookup[fId]);
 
-    inIt.GoToBegin();
-    outIt.GoToBegin();
-
     while (!outIt.IsAtEnd())
     {
       if (inIt.Get() < 0)
diff --git a/Modules/Numerics/Statistics/include/itkImageClassifierFilter.hxx b/Modules/Numerics/Statistics/include/itkImageClassifierFilter.hxx
@@ -149,9 +149,6 @@ ImageClassifierFilter<TSample, TInputImage, TOutputImage>::GenerateData()
   ImageRegionConstIterator inpItr(inputImage, inputImage->GetBufferedRegion());
   ImageRegionIterator      outItr(outputImage, outputImage->GetBufferedRegion());
 
-  inpItr.GoToBegin();
-  outItr.GoToBegin();
-
   while (!inpItr.IsAtEnd())
   {
     MeasurementVectorType measurements;
diff --git a/Modules/Numerics/Statistics/include/itkImageToHistogramFilter.hxx b/Modules/Numerics/Statistics/include/itkImageToHistogramFilter.hxx
@@ -199,8 +199,7 @@ ImageToHistogramFilter<TImage>::ThreadedComputeMinimumAndMaximum(const RegionTyp
   HistogramMeasurementVectorType min(nbOfComponents);
   HistogramMeasurementVectorType max(nbOfComponents);
 
-  ImageRegionConstIterator inputIt(this->GetInput(), inputRegionForThread);
-  inputIt.GoToBegin();
+  ImageRegionConstIterator       inputIt(this->GetInput(), inputRegionForThread);
   HistogramMeasurementVectorType m(nbOfComponents);
 
   min.Fill(NumericTraits<ValueType>::max());
@@ -238,8 +237,7 @@ ImageToHistogramFilter<TImage>::ThreadedStreamedGenerateData(const RegionType &
   histogram->Initialize(outputHistogram->GetSize(), m_Minimum, m_Maximum);
 
 
-  ImageRegionConstIterator inputIt(this->GetInput(), inputRegionForThread);
-  inputIt.GoToBegin();
+  ImageRegionConstIterator       inputIt(this->GetInput(), inputRegionForThread);
   HistogramMeasurementVectorType m(nbOfComponents);
 
   typename HistogramType::IndexType index;
diff --git a/Modules/Numerics/Statistics/include/itkImageToListSampleFilter.hxx b/Modules/Numerics/Statistics/include/itkImageToListSampleFilter.hxx
@@ -125,12 +125,10 @@ ImageToListSampleFilter<TImage, TMaskImage>::GenerateData()
   output->Clear();
 
   ImageRegionConstIterator it(input, input->GetBufferedRegion());
-  it.GoToBegin();
 
   if (maskImage) // mask specified
   {
     ImageRegionConstIterator mit(maskImage, maskImage->GetBufferedRegion());
-    mit.GoToBegin();
     while (!it.IsAtEnd())
     {
       if (mit.Get() == this->m_MaskValue)
diff --git a/Modules/Numerics/Statistics/include/itkMaskedImageToHistogramFilter.hxx b/Modules/Numerics/Statistics/include/itkMaskedImageToHistogramFilter.hxx
@@ -41,10 +41,8 @@ MaskedImageToHistogramFilter<TImage, TMaskImage>::ThreadedComputeMinimumAndMaxim
 
   const MaskPixelType maskValue = this->GetMaskValue();
 
-  ImageRegionConstIterator inputIt(this->GetInput(), inputRegionForThread);
-  ImageRegionConstIterator maskIt(this->GetMaskImage(), inputRegionForThread);
-  inputIt.GoToBegin();
-  maskIt.GoToBegin();
+  ImageRegionConstIterator       inputIt(this->GetInput(), inputRegionForThread);
+  ImageRegionConstIterator       maskIt(this->GetMaskImage(), inputRegionForThread);
   HistogramMeasurementVectorType m(nbOfComponents);
 
   min.Fill(NumericTraits<ValueType>::max());
@@ -84,10 +82,8 @@ MaskedImageToHistogramFilter<TImage, TMaskImage>::ThreadedStreamedGenerateData(c
   histogram->SetMeasurementVectorSize(nbOfComponents);
   histogram->Initialize(outputHistogram->GetSize(), this->m_Minimum, this->m_Maximum);
 
-  ImageRegionConstIterator inputIt(this->GetInput(), inputRegionForThread);
-  ImageRegionConstIterator maskIt(this->GetMaskImage(), inputRegionForThread);
-  inputIt.GoToBegin();
-  maskIt.GoToBegin();
+  ImageRegionConstIterator       inputIt(this->GetInput(), inputRegionForThread);
+  ImageRegionConstIterator       maskIt(this->GetMaskImage(), inputRegionForThread);
   HistogramMeasurementVectorType m(nbOfComponents);
   const MaskPixelType            maskValue = this->GetMaskValue();
 
diff --git a/Modules/Registration/Common/include/itkHistogramImageToImageMetric.hxx b/Modules/Registration/Common/include/itkHistogramImageToImageMetric.hxx
@@ -89,9 +89,8 @@ HistogramImageToImageMetric<TFixedImage, TMovingImage>::Initialize()
     // Calculate min and max image values in fixed image.
     const FixedImageConstPointerType pFixedImage = this->m_FixedImage;
     ImageRegionConstIterator         fiIt(pFixedImage, pFixedImage->GetBufferedRegion());
-    fiIt.GoToBegin();
-    FixedImagePixelType minFixed = fiIt.Value();
-    FixedImagePixelType maxFixed = fiIt.Value();
+    FixedImagePixelType              minFixed = fiIt.Value();
+    FixedImagePixelType              maxFixed = fiIt.Value();
     ++fiIt;
     while (!fiIt.IsAtEnd())
     {
@@ -112,9 +111,8 @@ HistogramImageToImageMetric<TFixedImage, TMovingImage>::Initialize()
     // Calculate min and max image values in moving image.
     const MovingImageConstPointerType pMovingImage = this->m_MovingImage;
     ImageRegionConstIterator          miIt(pMovingImage, pMovingImage->GetBufferedRegion());
-    miIt.GoToBegin();
-    MovingImagePixelType minMoving = miIt.Value();
-    MovingImagePixelType maxMoving = miIt.Value();
+    MovingImagePixelType              minMoving = miIt.Value();
+    MovingImagePixelType              maxMoving = miIt.Value();
     ++miIt;
     while (!miIt.IsAtEnd())
     {
diff --git a/Modules/Registration/Metricsv4/include/itkJointHistogramMutualInformationComputeJointPDFThreaderBase.hxx b/Modules/Registration/Metricsv4/include/itkJointHistogramMutualInformationComputeJointPDFThreaderBase.hxx
@@ -125,7 +125,6 @@ JointHistogramMutualInformationComputeJointPDFThreaderBase<TDomainPartitioner,
   }
 
   ImageRegionIterator jointPDFIt(this->m_Associate->m_JointPDF, this->m_Associate->m_JointPDF->GetBufferedRegion());
-  jointPDFIt.GoToBegin();
   std::vector<ImageRegionConstIterator<JointHistogramType>> jointHistogramPerThreadIts;
   for (ThreadIdType i = 0; i < numberOfWorkUnitsUsed; ++i)
   {
diff --git a/Modules/Segmentation/Classifiers/include/itkBayesianClassifierImageFilter.hxx b/Modules/Segmentation/Classifiers/include/itkBayesianClassifierImageFilter.hxx
@@ -285,8 +285,6 @@ BayesianClassifierImageFilter<TInputVectorImage, TLabelsType, TPosteriorsPrecisi
 
       itrPosteriorImage.GoToBegin();
       ImageRegionIterator it(extractedComponentImage, extractedComponentImage->GetBufferedRegion());
-
-      it.GoToBegin();
       while (!itrPosteriorImage.IsAtEnd())
       {
         it.Set(itrPosteriorImage.Get()[componentToExtract]);
@@ -301,7 +299,6 @@ BayesianClassifierImageFilter<TInputVectorImage, TLabelsType, TPosteriorsPrecisi
       itrPosteriorImage.GoToBegin();
 
       ImageRegionIterator sit(m_SmoothingFilter->GetOutput(), m_SmoothingFilter->GetOutput()->GetBufferedRegion());
-      sit.GoToBegin();
       while (!itrPosteriorImage.IsAtEnd())
       {
         PosteriorsPixelType posteriorPixel = itrPosteriorImage.Get();
diff --git a/Modules/Segmentation/Classifiers/include/itkBayesianClassifierInitializationImageFilter.hxx b/Modules/Segmentation/Classifiers/include/itkBayesianClassifierInitializationImageFilter.hxx
@@ -137,9 +137,6 @@ BayesianClassifierInitializationImageFilter<TInputImage, TProbabilityPrecisionTy
   const typename InputImageType::RegionType imageRegion = inputImage->GetLargestPossibleRegion();
   InputImageIteratorType                    itrInputImage(inputImage, imageRegion);
 
-  itrInputImage.GoToBegin();
-  itrKMeansImage.GoToBegin();
-
   // find sumsOfSquares, sums, and classCount by indexing using the kmeans
   // output labelmap
   while (!itrInputImage.IsAtEnd())
diff --git a/Modules/Segmentation/Classifiers/include/itkScalarImageKmeansImageFilter.hxx b/Modules/Segmentation/Classifiers/include/itkScalarImageKmeansImageFilter.hxx
@@ -161,7 +161,6 @@ ScalarImageKmeansImageFilter<TInputImage, TOutputImage>::GenerateData()
   }
 
   ImageRegionIterator pixel(outputPtr, region);
-  pixel.GoToBegin();
 
   using ClassifierOutputType = typename ClassifierType::MembershipSampleType;
   const ClassifierOutputType * membershipSample = classifier->GetOutput();
diff --git a/Modules/Segmentation/ConnectedComponents/include/itkConnectedComponentFunctorImageFilter.hxx b/Modules/Segmentation/ConnectedComponents/include/itkConnectedComponentFunctorImageFilter.hxx
@@ -98,8 +98,6 @@ ConnectedComponentFunctorImageFilter<TInputImage, TOutputImage, TFunctor, TMaskI
   if (const typename TMaskImage::ConstPointer mask = this->GetMaskImage())
   {
     ImageRegionConstIterator mit(mask, output->GetRequestedRegion());
-
-    mit.GoToBegin();
     oit.GoToBegin();
     while (!mit.IsAtEnd())
     {
@@ -120,7 +118,6 @@ ConnectedComponentFunctorImageFilter<TInputImage, TOutputImage, TFunctor, TMaskI
   // current pixel
   init.GoToBegin();
   onit.GoToBegin();
-  it.GoToBegin();
   oit.GoToBegin();
   OutputPixelType maxLabel{};
   // create an equivalency table
diff --git a/Modules/Segmentation/ConnectedComponents/include/itkHardConnectedComponentImageFilter.hxx b/Modules/Segmentation/ConnectedComponents/include/itkHardConnectedComponentImageFilter.hxx
@@ -49,8 +49,6 @@ HardConnectedComponentImageFilter<TInputImage, TOutputImage>::GenerateData()
   ImageRegionIteratorWithIndex ot(output, output->GetRequestedRegion());
 
   ProgressReporter progress(this, 0, output->GetRequestedRegion().GetNumberOfPixels());
-  it.GoToBegin();
-  ot.GoToBegin();
   for (; !it.IsAtEnd(); ++it, ++ot)
   {
     if (Math::NotExactlyEquals(it.Get(), typename ImageRegionConstIterator<TInputImage>::PixelType{}))
diff --git a/Modules/Segmentation/LevelSets/include/itkExtensionVelocitiesImageFilter.hxx b/Modules/Segmentation/LevelSets/include/itkExtensionVelocitiesImageFilter.hxx
diff --git a/Modules/Segmentation/Watersheds/include/itkIsolatedWatershedImageFilter.hxx b/Modules/Segmentation/Watersheds/include/itkIsolatedWatershedImageFilter.hxx
diff --git a/Modules/Segmentation/Watersheds/include/itkWatershedEquivalenceRelabeler.hxx b/Modules/Segmentation/Watersheds/include/itkWatershedEquivalenceRelabeler.hxx
diff --git a/Modules/Segmentation/Watersheds/include/itkWatershedRelabeler.hxx b/Modules/Segmentation/Watersheds/include/itkWatershedRelabeler.hxx
diff --git a/Modules/Segmentation/Watersheds/include/itkWatershedSegmenter.hxx b/Modules/Segmentation/Watersheds/include/itkWatershedSegmenter.hxx

Original file line number	Diff line number	Diff line change
`@@ -85,7 +85,6 @@ MRIBiasEnergyFunction<TImage, TImageMask, TBiasField>::GetValue(const Parameters`
`85`	`85`	`typename TBiasField::SimpleForwardIterator bIter(m_BiasField);`
`86`	`86`
`87`	`87`	`bIter.Begin();`
`88`		`- iIter.GoToBegin();`
`89`	`88`	`total = 0.0;`
`90`	`89`
`91`	`90`	`// Fastest for full sampling`
`@@ -102,7 +101,6 @@ MRIBiasEnergyFunction<TImage, TImageMask, TBiasField>::GetValue(const Parameters`
`102`	`101`	`else`
`103`	`102`	`{`
`104`	`103`	`ImageRegionIterator mIter(m_Mask, m_Region);`
`105`		`- mIter.GoToBegin();`
`106`	`104`	`while (!bIter.IsAtEnd())`
`107`	`105`	`{`
`108`	`106`	`if (mIter.Get() > 0.0)`
`@@ -547,7 +545,6 @@ MRIBiasFieldCorrectionFilter<TInputImage, TOutputImage, TMaskImage>::CorrectImag`
`547`	`545`	`{`
`548`	`546`	`itkDebugMacro("Output mask is being used");`
`549`	`547`	`ImageRegionIterator mIter(m_OutputMask, region);`
`550`		`- mIter.GoToBegin();`
`551`	`548`	`while (!bIter.IsAtEnd())`
`552`	`549`	`{`
`553`	`550`	`const double inputPixel = iIter.Get();`
Original file line number	Diff line number	Diff line change
`@@ -658,8 +658,6 @@ BSplineScatteredDataPointSetToImageFilter<TInputPointSet, TOutputImage>::AfterTh`
`658`	`658`
`659`	`659`	`ItD.GoToBegin();`
`660`	`660`	`ItO.GoToBegin();`
`661`		`- Itd.GoToBegin();`
`662`		`- Ito.GoToBegin();`
`663`	`661`	`while (!ItD.IsAtEnd())`
`664`	`662`	`{`
`665`	`663`	`ItD.Set(ItD.Get() + Itd.Get());`
Original file line number	Diff line number	Diff line change
`@@ -136,8 +136,7 @@ DiscreteGaussianDerivativeImageFunction<TInputImage, TOutput>::RecomputeGaussian`
`136`	`136`
`137`	`137`	`// Copy kernel image to neighborhood. Do not copy boundaries.`
`138`	`138`	`ImageRegionConstIterator it(kernelImage, kernelRegion);`
`139`		`- it.GoToBegin();`
`140`		`- unsigned int idx = 0;`
	`139`	`+ unsigned int idx = 0;`
`141`	`140`
`142`	`141`	`while (!it.IsAtEnd())`
`143`	`142`	`{`
Original file line number	Diff line number	Diff line change
`@@ -158,8 +158,7 @@ DiscreteGradientMagnitudeGaussianImageFunction<TInputImage, TOutput>::RecomputeG`
`158`	`158`
`159`	`159`	`// Copy kernel image to neighborhood. Do not copy boundaries.`
`160`	`160`	`ImageRegionConstIterator it(kernelImage, kernelRegion);`
`161`		`- it.GoToBegin();`
`162`		`- unsigned int idx = 0;`
	`161`	`+ unsigned int idx = 0;`
`163`	`162`
`164`	`163`	`while (!it.IsAtEnd())`
`165`	`164`	`{`
Original file line number	Diff line number	Diff line change
`@@ -164,8 +164,7 @@ DiscreteHessianGaussianImageFunction<TInputImage, TOutput>::RecomputeGaussianKer`
`164`	`164`
`165`	`165`	`// Copy kernel image to neighborhood. Do not copy boundaries.`
`166`	`166`	`ImageRegionConstIterator it(kernelImage, kernelRegion);`
`167`		`- it.GoToBegin();`
`168`		`- unsigned int idx = 0;`
	`167`	`+ unsigned int idx = 0;`
`169`	`168`	`while (!it.IsAtEnd())`
`170`	`169`	`{`
`171`	`170`	`m_KernelArray[kernelidx][idx] = it.Get();`
Original file line number	Diff line number	Diff line change
`@@ -437,8 +437,6 @@ HistogramMatchingImageFilter<TInputImage, TOutputImage, THistogramMeasurement>::`
`437`	`437`	`{`
`438`	`438`	`// put each image pixel into the histogram`
`439`	`439`	`ImageRegionConstIterator iter(image, image->GetBufferedRegion());`
`440`		`-`
`441`		`- iter.GoToBegin();`
`442`	`440`	`while (!iter.IsAtEnd())`
`443`	`441`	`{`
`444`	`442`	`const InputPixelType & value = iter.Value();`
Original file line number	Diff line number	Diff line change
`@@ -181,8 +181,7 @@ AccumulateImageFilter<TInputImage, TOutputImage>::GenerateData()`
`181`	`181`	`}`
`182`	`182`	`const typename TInputImage::RegionType AccumulatedRegion(AccumulatedIndex, AccumulatedSize);`
`183`	`183`	`ImageRegionConstIterator inputIter(inputImage, AccumulatedRegion);`
`184`		`- inputIter.GoToBegin();`
`185`		`- AccumulateType Value{};`
	`184`	`+ AccumulateType Value{};`
`186`	`185`	`while (!inputIter.IsAtEnd())`
`187`	`186`	`{`
`188`	`187`	`Value += static_cast<AccumulateType>(inputIter.Get());`
Original file line number	Diff line number	Diff line change
`@@ -149,7 +149,6 @@ ImagePCAShapeModelEstimator<TInputImage, TOutputImage>::GenerateData()`
`149`	`149`	`ImageRegionIterator outIterJ(this->GetOutput(j), region);`
`150`	`150`
`151`	`151`	`unsigned int idx = 0;`
`152`		`- outIterJ.GoToBegin();`
`153`	`152`	`while (!outIterJ.IsAtEnd())`
`154`	`153`	`{`
`155`	`154`	`outIterJ.Set(static_cast<typename OutputImageType::PixelType>(m_OneEigenVector[idx]));`
`@@ -166,8 +165,6 @@ ImagePCAShapeModelEstimator<TInputImage, TOutputImage>::GenerateData()`
`166`	`165`	`{`
`167`	`166`	`region = this->GetOutput(j)->GetRequestedRegion();`
`168`	`167`	`ImageRegionIterator outIterJ(this->GetOutput(j), region);`
`169`		`-`
`170`		`- outIterJ.GoToBegin();`
`171`	`168`	`while (!outIterJ.IsAtEnd())`
`172`	`169`	`{`
`173`	`170`	`outIterJ.Set(0);`