ENH: Add GDCM test for 32 bits stored DICOM image

GDCM used to crash the application when it encountered a DICOM image that had 32 bits stored.
This new test checks that GDCM does not crash anymore for such images.

Modules/IO/GDCM/test/CMakeLists.txt

@@ -13,7 +13,8 @@ set(ITKIOGDCMTests
     itkGDCMImageOrientationPatientTest.cxx
     itkGDCMLoadImageSpacingTest.cxx
     itkGDCMLegacyMultiFrameTest.cxx
-    itkGDCMImageIONoPreambleTest.cxx)
+    itkGDCMImageIONoPreambleTest.cxx
+    itkGDCMImageIO32bitsStoredTest.cxx)
 
 createtestdriver(ITKIOGDCM "${ITKIOGDCM-Test_LIBRARIES}" "${ITKIOGDCMTests}")
 
@@ -265,6 +266,14 @@ itk_add_test(
   itkGDCMImageIONoPreambleTest
   DATA{Input/preamble.dcm})
 
+itk_add_test(
+  NAME
+  itkGDCMImageIO32bitsStoredTest
+  COMMAND
+  ITKIOGDCMTestDriver
+  itkGDCMImageIO32bitsStoredTest
+  DATA{Input/image32bitsStoredTest.dcm})
+
 itk_add_test(
   NAME
   itkGDCMImageReadWriteTest_RGB

Modules/IO/GDCM/test/Input/image32bitsStoredTest.dcm.cid

@@ -0,0 +1 @@
+bafybeihwecykoje6v7td77gmoe7m4ogghvxv4mgtkgdbnyo74jfeut7n44

Modules/IO/GDCM/test/itkGDCMImageIO32bitsStoredTest.cxx

@@ -0,0 +1,79 @@
+/*=========================================================================
+ *
+ *  Copyright NumFOCUS
+ *
+ *  Licensed under the Apache License, Version 2.0 (the "License");
+ *  you may not use this file except in compliance with the License.
+ *  You may obtain a copy of the License at
+ *
+ *         https://www.apache.org/licenses/LICENSE-2.0.txt
+ *
+ *  Unless required by applicable law or agreed to in writing, software
+ *  distributed under the License is distributed on an "AS IS" BASIS,
+ *  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+ *  See the License for the specific language governing permissions and
+ *  limitations under the License.
+ *
+ *=========================================================================*/
+
+#include "itkGDCMImageIO.h"
+#include "itkImageFileReader.h"
+#include "itkTestingMacros.h"
+
+// Specific ImageIO test
+
+/** This test verifies itkGDCMImageIO does not crash when it encounters
+ * a file that has 32 bits allocated for a pixel.
+ */
+int
+itkGDCMImageIO32bitsStoredTest(int argc, char * argv[])
+{
+  if (argc < 2)
+  {
+    std::cerr << "Missing parameters." << std::endl;
+    std::cerr << "Usage: " << itkNameOfTestExecutableMacro(argv);
+    std::cerr << " DicomImage" << std::endl;
+    return EXIT_FAILURE;
+  }
+
+  using InputPixelType = short;
+  using InputImageType = itk::Image<InputPixelType, 3>;
+  using ReaderType = itk::ImageFileReader<InputImageType>;
+  using ImageIOType = itk::GDCMImageIO;
+
+  auto dcmImageIO = ImageIOType::New();
+  bool canRead = dcmImageIO->CanReadFile(argv[1]);
+  std::cerr << "GDCM can read file: " << (canRead ? "yes" : "no") << std::endl;
+  if (canRead)
+  {
+    // All we test is that CanReadFile() does not crash, it is fine if it cannot read the file.
+    return EXIT_SUCCESS;
+  }
+
+  auto reader = ReaderType::New();
+  reader->SetFileName(argv[1]);
+  reader->SetImageIO(dcmImageIO);
+
+  try
+  {
+    reader->Update();
+  }
+  catch (const itk::ExceptionObject & e)
+  {
+    std::cerr << "exception in file reader " << std::endl;
+    std::cerr << e << std::endl;
+    return EXIT_FAILURE;
+  }
+
+  InputImageType::SizeType extentSize;
+  extentSize = reader->GetOutput()->GetLargestPossibleRegion().GetSize();
+  std::cout << "Read image dimensions: (" << extentSize[0] << ", " << extentSize[1] << ", " << extentSize[2] << ')'
+            << std::endl;
+  if (extentSize[0] == 0 || extentSize[1] == 0 || extentSize[2] == 0)
+  {
+    std::cerr << "File read but empty " << std::endl;
+    return EXIT_FAILURE;
+  }
+
+  return EXIT_SUCCESS;
+}