diff --git a/Modules/Filtering/Denoising/include/itkPatchBasedDenoisingBaseImageFilter.hxx b/Modules/Filtering/Denoising/include/itkPatchBasedDenoisingBaseImageFilter.hxx
index 7b784ddca8e..bb499f1b3c5 100644
--- a/Modules/Filtering/Denoising/include/itkPatchBasedDenoisingBaseImageFilter.hxx
+++ b/Modules/Filtering/Denoising/include/itkPatchBasedDenoisingBaseImageFilter.hxx
@@ -174,12 +174,7 @@ PatchBasedDenoisingBaseImageFilter<TInputImage, TOutputImage>::GetPatchLengthInV
 {
   const PatchRadiusType diameter = this->GetPatchDiameterInVoxels();
 
-  typename PatchRadiusType::SizeValueType length = 1;
-  for (unsigned int dim = 0; dim < ImageDimension; ++dim)
-  {
-    length *= diameter[dim];
-  }
-  return length;
+  return diameter.CalculateProductOfElements();
 }
 
 template <typename TInputImage, typename TOutputImage>
diff --git a/Modules/Filtering/ImageStatistics/include/itkImagePCADecompositionCalculator.hxx b/Modules/Filtering/ImageStatistics/include/itkImagePCADecompositionCalculator.hxx
index a94b4a54104..34d9b9047b8 100644
--- a/Modules/Filtering/ImageStatistics/include/itkImagePCADecompositionCalculator.hxx
+++ b/Modules/Filtering/ImageStatistics/include/itkImagePCADecompositionCalculator.hxx
@@ -63,11 +63,7 @@ void
 ImagePCADecompositionCalculator<TInputImage, TBasisImage>::CalculateBasisMatrix()
 {
   m_Size = m_BasisImages[0]->GetRequestedRegion().GetSize();
-  m_NumPixels = 1;
-  for (unsigned int i = 0; i < BasisImageDimension; ++i)
-  {
-    m_NumPixels *= m_Size[i];
-  }
+  m_NumPixels = m_Size.CalculateProductOfElements();
 
   m_BasisMatrix = BasisMatrixType(static_cast<unsigned int>(m_BasisImages.size()), m_NumPixels);
 
diff --git a/Modules/Filtering/ImageStatistics/include/itkImagePCAShapeModelEstimator.hxx b/Modules/Filtering/ImageStatistics/include/itkImagePCAShapeModelEstimator.hxx
index 7ae4c1b4982..548ef7ce0cf 100644
--- a/Modules/Filtering/ImageStatistics/include/itkImagePCAShapeModelEstimator.hxx
+++ b/Modules/Filtering/ImageStatistics/include/itkImagePCAShapeModelEstimator.hxx
@@ -271,11 +271,7 @@ ImagePCAShapeModelEstimator<TInputImage, TOutputImage>::CalculateInnerProduct()
   // training data
   m_InputImageSize = (inputImagePointerArray[0])->GetBufferedRegion().GetSize();
 
-  m_NumberOfPixels = 1;
-  for (unsigned int i = 0; i < InputImageDimension; ++i)
-  {
-    m_NumberOfPixels *= m_InputImageSize[i];
-  }
+  m_NumberOfPixels = m_InputImageSize.CalculateProductOfElements();
 
   // Calculate the means
   m_Means.set_size(m_NumberOfPixels);
diff --git a/Modules/IO/ImageBase/test/itkLargeImageWriteReadTest.cxx b/Modules/IO/ImageBase/test/itkLargeImageWriteReadTest.cxx
index dc65cb519b4..20cc46b6114 100644
--- a/Modules/IO/ImageBase/test/itkLargeImageWriteReadTest.cxx
+++ b/Modules/IO/ImageBase/test/itkLargeImageWriteReadTest.cxx
@@ -47,11 +47,7 @@ ActualTest(std::string filename, typename TImageType::SizeType size)
 
     image->SetRegions(region);
 
-    size_t numberOfPixels = 1;
-    for (unsigned int i = 0; i < ImageType::ImageDimension; ++i)
-    {
-      numberOfPixels *= region.GetSize(i);
-    }
+    const size_t numberOfPixels = size.CalculateProductOfElements();
 
     const size_t sizeInMebiBytes = sizeof(PixelType) * numberOfPixels / (1024 * 1024);
 
diff --git a/Modules/IO/Meta/test/itkLargeMetaImageWriteReadTest.cxx b/Modules/IO/Meta/test/itkLargeMetaImageWriteReadTest.cxx
index 7b012884784..ea6b7ae21b5 100644
--- a/Modules/IO/Meta/test/itkLargeMetaImageWriteReadTest.cxx
+++ b/Modules/IO/Meta/test/itkLargeMetaImageWriteReadTest.cxx
@@ -53,11 +53,7 @@ ActualTest(std::string filename, typename TImageType::SizeType size)
 
     image->SetRegions(region);
 
-    size_t numberOfPixels = 1;
-    for (unsigned int i = 0; i < ImageType::ImageDimension; ++i)
-    {
-      numberOfPixels *= region.GetSize(i);
-    }
+    const size_t numberOfPixels = size.CalculateProductOfElements();
 
     const size_t sizeInMebiBytes = sizeof(PixelType) * numberOfPixels / (1024 * 1024);
 
diff --git a/Modules/IO/TIFF/test/itkLargeTIFFImageWriteReadTest.cxx b/Modules/IO/TIFF/test/itkLargeTIFFImageWriteReadTest.cxx
index 7957c52b368..d15ee4acae7 100644
--- a/Modules/IO/TIFF/test/itkLargeTIFFImageWriteReadTest.cxx
+++ b/Modules/IO/TIFF/test/itkLargeTIFFImageWriteReadTest.cxx
@@ -50,11 +50,7 @@ itkLargeTIFFImageWriteReadTestHelper(std::string filename, typename TImage::Size
 
     image->SetRegions(region);
 
-    SizeValueType numberOfPixels = 1;
-    for (unsigned int i = 0; i < ImageType::ImageDimension; ++i)
-    {
-      numberOfPixels *= static_cast<SizeValueType>(region.GetSize(i));
-    }
+    const SizeValueType numberOfPixels = size.CalculateProductOfElements();
 
     constexpr SizeValueType oneMebiByte = static_cast<SizeValueType>(1024) * static_cast<SizeValueType>(1024);