Merge pull request InsightSoftwareConsortium#5541 from N-Dekker/Use-CalculateProductOfElements

Use CalculateProductOfElements in Denoising, ImageStatistics, and IO tests

Author: thewtex

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>

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);
 

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);

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);
 

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);
 

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);