STYLE: Use itk::ReadImage and itk::WriteImage in src/Filtering examples

Author: josetascon

src/Filtering/AnisotropicSmoothing/ComputeCurvatureAnisotropicDiffusion/Code.cxx

@@ -48,13 +48,11 @@ main(int argc, char * argv[])
   const InputPixelType timeStep = std::stod(argv[4]);
   const InputPixelType conductance = std::stod(argv[5]);
 
-  using ReaderType = itk::ImageFileReader<InputImageType>;
-  ReaderType::Pointer reader = ReaderType::New();
-  reader->SetFileName(inputFileName);
+  const auto input = itk::ReadImage<InputImageType>(inputFileName);
 
   using FilterType = itk::CurvatureAnisotropicDiffusionImageFilter<InputImageType, InputImageType>;
   FilterType::Pointer filter = FilterType::New();
-  filter->SetInput(reader->GetOutput());
+  filter->SetInput(input);
   filter->SetNumberOfIterations(numberOfIterations);
   filter->SetTimeStep(timeStep);
   filter->SetConductanceParameter(conductance);
@@ -65,14 +63,9 @@ main(int argc, char * argv[])
   rescaler->SetOutputMinimum(itk::NumericTraits<OutputPixelType>::min());
   rescaler->SetOutputMaximum(itk::NumericTraits<OutputPixelType>::max());
 
-  using WriterType = itk::ImageFileWriter<OutputImageType>;
-  WriterType::Pointer writer = WriterType::New();
-  writer->SetFileName(outputFileName);
-  writer->SetInput(rescaler->GetOutput());
-
   try
   {
-    writer->Update();
+    itk::WriteImage(rescaler->GetOutput(), outputFileName);
   }
   catch (itk::ExceptionObject & error)
   {

src/Filtering/AnisotropicSmoothing/ComputeCurvatureFlow/Code.cxx

@@ -47,13 +47,11 @@ main(int argc, char * argv[])
   const int            numberOfIterations = std::stoi(argv[3]);
   const InputPixelType timeStep = std::stod(argv[4]);
 
-  using ReaderType = itk::ImageFileReader<InputImageType>;
-  ReaderType::Pointer reader = ReaderType::New();
-  reader->SetFileName(inputFileName);
+  const auto input = itk::ReadImage<InputImageType>(inputFileName);
 
   using FilterType = itk::CurvatureFlowImageFilter<InputImageType, InputImageType>;
   FilterType::Pointer filter = FilterType::New();
-  filter->SetInput(reader->GetOutput());
+  filter->SetInput(input);
   filter->SetNumberOfIterations(numberOfIterations);
   filter->SetTimeStep(timeStep);
 
@@ -63,14 +61,9 @@ main(int argc, char * argv[])
   rescaler->SetOutputMinimum(itk::NumericTraits<OutputPixelType>::min());
   rescaler->SetOutputMaximum(itk::NumericTraits<OutputPixelType>::max());
 
-  using WriterType = itk::ImageFileWriter<OutputImageType>;
-  WriterType::Pointer writer = WriterType::New();
-  writer->SetFileName(outputFileName);
-  writer->SetInput(rescaler->GetOutput());
-
   try
   {
-    writer->Update();
+    itk::WriteImage(rescaler->GetOutput(), outputFileName);
   }
   catch (itk::ExceptionObject & error)
   {

src/Filtering/AnisotropicSmoothing/ComputeGradientAnisotropicDiffusion/Code.cxx

@@ -48,13 +48,11 @@ main(int argc, char * argv[])
   const InputPixelType timeStep = std::stod(argv[4]);
   const InputPixelType conductance = std::stod(argv[5]);
 
-  using ReaderType = itk::ImageFileReader<InputImageType>;
-  ReaderType::Pointer reader = ReaderType::New();
-  reader->SetFileName(inputFileName);
+  const auto input = itk::ReadImage<InputImageType>(inputFileName);
 
   using FilterType = itk::GradientAnisotropicDiffusionImageFilter<InputImageType, InputImageType>;
   FilterType::Pointer filter = FilterType::New();
-  filter->SetInput(reader->GetOutput());
+  filter->SetInput(input);
   filter->SetNumberOfIterations(numberOfIterations);
   filter->SetTimeStep(timeStep);
   filter->SetConductanceParameter(conductance);
@@ -65,14 +63,9 @@ main(int argc, char * argv[])
   rescaler->SetOutputMinimum(itk::NumericTraits<OutputPixelType>::min());
   rescaler->SetOutputMaximum(itk::NumericTraits<OutputPixelType>::max());
 
-  using WriterType = itk::ImageFileWriter<OutputImageType>;
-  WriterType::Pointer writer = WriterType::New();
-  writer->SetFileName(outputFileName);
-  writer->SetInput(rescaler->GetOutput());
-
   try
   {
-    writer->Update();
+    itk::WriteImage(rescaler->GetOutput(), outputFileName);
   }
   catch (itk::ExceptionObject & error)
   {

src/Filtering/AnisotropicSmoothing/ComputePeronaMalikAnisotropicDiffusion/Code.cxx

@@ -40,27 +40,20 @@ main(int argc, char * argv[])
   using InputPixelType = unsigned char;
   using InputImageType = itk::Image<InputPixelType, Dimension>;
 
-  using ReaderType = itk::ImageFileReader<InputImageType>;
-  ReaderType::Pointer reader = ReaderType::New();
-  reader->SetFileName(argv[1]);
+  const auto input = itk::ReadImage<InputImageType>(argv[1]);
 
   using OutputPixelType = float;
   using OutputImageType = itk::Image<OutputPixelType, Dimension>;
   using FilterType = itk::GradientAnisotropicDiffusionImageFilter<InputImageType, OutputImageType>;
   FilterType::Pointer filter = FilterType::New();
-  filter->SetInput(reader->GetOutput());
+  filter->SetInput(input);
   filter->SetNumberOfIterations(std::stoi(argv[3]));
   filter->SetTimeStep(0.125);
   filter->SetConductanceParameter(std::stod(argv[4]));
 
-  using WriterType = itk::ImageFileWriter<OutputImageType>;
-  WriterType::Pointer writer = WriterType::New();
-  writer->SetFileName(argv[2]);
-  writer->SetInput(filter->GetOutput());
-
   try
   {
-    writer->Update();
+    itk::WriteImage(filter->GetOutput(), argv[2]);
   }
   catch (itk::ExceptionObject & error)
   {

src/Filtering/AnisotropicSmoothing/SmoothImageWhilePreservingEdges/Code.cxx

@@ -51,14 +51,12 @@ main(int argc, char * argv[])
   using RGBImageType = itk::Image<itk::RGBPixel<float>, 2>;
 
   // 1) Read the RGB image
-  using ReaderType = itk::ImageFileReader<RGBImageType>;
-  ReaderType::Pointer reader = ReaderType::New();
-  reader->SetFileName(inputFileName);
+  const auto input = itk::ReadImage<RGBImageType>(inputFileName);
 
   // 2) Cast to Vector image for processing
   using AdaptorInputType = itk::RGBToVectorImageAdaptor<RGBImageType>;
   AdaptorInputType::Pointer adaptInput = AdaptorInputType::New();
-  adaptInput->SetImage(reader->GetOutput());
+  adaptInput->SetImage(input);
   using CastInputType = itk::CastImageFilter<AdaptorInputType, FloatImageType>;
   CastInputType::Pointer castInput = CastInputType::New();
   castInput->SetInput(adaptInput);
@@ -90,7 +88,7 @@ main(int argc, char * argv[])
   // 5) Display the input and smoothed images
 #ifdef ENABLE_QUICKVIEW
   QuickView viewer;
-  viewer.AddRGBImage(reader->GetOutput(), true, itksys::SystemTools::GetFilenameName(inputFileName));
+  viewer.AddRGBImage(input, true, itksys::SystemTools::GetFilenameName(inputFileName));
 
   std::stringstream desc;
   desc << "VectorGradientAnisotropicDiffusionImageFilter\niterations: " << filter->GetNumberOfIterations()

src/Filtering/AnisotropicSmoothing/SmoothImageWhilePreservingEdges2/Code.cxx

@@ -50,14 +50,12 @@ main(int argc, char * argv[])
   using RGBImageType = itk::Image<itk::RGBPixel<float>, 2>;
 
   // 1) Read the RGB image
-  using ReaderType = itk::ImageFileReader<RGBImageType>;
-  ReaderType::Pointer reader = ReaderType::New();
-  reader->SetFileName(inputFileName);
+  const auto input = itk::ReadImage<RGBImageType>(inputFileName);
 
   // 2) Cast to Vector image for processing
   using AdaptorInputType = itk::RGBToVectorImageAdaptor<RGBImageType>;
   AdaptorInputType::Pointer adaptInput = AdaptorInputType::New();
-  adaptInput->SetImage(reader->GetOutput());
+  adaptInput->SetImage(input);
   using CastInputType = itk::CastImageFilter<AdaptorInputType, FloatImageType>;
   CastInputType::Pointer castInput = CastInputType::New();
   castInput->SetInput(adaptInput);
@@ -89,7 +87,7 @@ main(int argc, char * argv[])
   // 5) Display the input and smoothed images
 #ifdef ENABLE_QUICKVIEW
   QuickView viewer;
-  viewer.AddRGBImage(reader->GetOutput(), true, itksys::SystemTools::GetFilenameName(inputFileName));
+  viewer.AddRGBImage(input, true, itksys::SystemTools::GetFilenameName(inputFileName));
 
   std::stringstream desc;
   desc << "VectorCurvatureAnisotropicDiffusionImageFilter\niterations: " << filter->GetNumberOfIterations()

src/Filtering/AntiAlias/SmoothBinaryImageBeforeSurfaceExtraction/Code.cxx

@@ -52,29 +52,21 @@ main(int argc, char * argv[])
   }
 
   constexpr unsigned int Dimension = 2;
-
   using PixelType = float;
-
   using ImageType = itk::Image<PixelType, Dimension>;
 
-  using ReaderType = itk::ImageFileReader<ImageType>;
-  ReaderType::Pointer reader = ReaderType::New();
-  reader->SetFileName(inputFileName);
+  const auto input = itk::ReadImage<ImageType>(inputFileName);
 
   using FilterType = itk::AntiAliasBinaryImageFilter<ImageType, ImageType>;
   FilterType::Pointer filter = FilterType::New();
-  filter->SetInput(reader->GetOutput());
+  filter->SetInput(input);
   filter->SetMaximumRMSError(maximumRMSError);
   filter->SetNumberOfIterations(numberOfIterations);
   filter->SetNumberOfLayers(numberOfLayers);
 
-  using WriterType = itk::ImageFileWriter<ImageType>;
-  WriterType::Pointer writer = WriterType::New();
-  writer->SetFileName(outputFileName);
-  writer->SetInput(filter->GetOutput());
   try
   {
-    writer->Update();
+    itk::WriteImage(filter->GetOutput(), outputFileName);
   }
   catch (itk::ExceptionObject & error)
   {

src/Filtering/BinaryMathematicalMorphology/ClosingBinaryImage/Code.cxx

@@ -45,11 +45,7 @@ main(int argc, char * argv[])
   }
   else
   {
-    using ReaderType = itk::ImageFileReader<ImageType>;
-    ReaderType::Pointer reader = ReaderType::New();
-    reader->SetFileName(argv[1]);
-    reader->Update();
-    image = reader->GetOutput();
+    image = itk::ReadImage<ImageType>(argv[1]);
   }
 
   unsigned int radius = 5;

src/Filtering/BinaryMathematicalMorphology/DilateABinaryImage/Code.cxx

@@ -38,11 +38,9 @@ main(int argc, char * argv[])
 
   using PixelType = unsigned char;
   constexpr unsigned int Dimension = 2;
-
   using ImageType = itk::Image<PixelType, Dimension>;
-  using ReaderType = itk::ImageFileReader<ImageType>;
-  ReaderType::Pointer reader = ReaderType::New();
-  reader->SetFileName(inputImage);
+
+  const auto input = itk::ReadImage<ImageType>(inputImage);
 
   using StructuringElementType = itk::FlatStructuringElement<Dimension>;
   StructuringElementType::RadiusType radius;
@@ -52,18 +50,13 @@ main(int argc, char * argv[])
   using BinaryDilateImageFilterType = itk::BinaryDilateImageFilter<ImageType, ImageType, StructuringElementType>;
 
   BinaryDilateImageFilterType::Pointer dilateFilter = BinaryDilateImageFilterType::New();
-  dilateFilter->SetInput(reader->GetOutput());
+  dilateFilter->SetInput(input);
   dilateFilter->SetKernel(structuringElement);
   dilateFilter->SetForegroundValue(255); // Value to dilate
 
-  using WriterType = itk::ImageFileWriter<ImageType>;
-  WriterType::Pointer writer = WriterType::New();
-  writer->SetInput(dilateFilter->GetOutput());
-  writer->SetFileName(outputImage);
-
   try
   {
-    writer->Update();
+    itk::WriteImage(dilateFilter->GetOutput(), outputImage);
   }
   catch (itk::ExceptionObject & error)
   {

src/Filtering/BinaryMathematicalMorphology/ErodeABinaryImage/Code.cxx

@@ -38,11 +38,9 @@ main(int argc, char * argv[])
 
   using PixelType = unsigned char;
   constexpr unsigned int Dimension = 2;
-
   using ImageType = itk::Image<PixelType, Dimension>;
-  using ReaderType = itk::ImageFileReader<ImageType>;
-  ReaderType::Pointer reader = ReaderType::New();
-  reader->SetFileName(inputImage);
+
+  const auto input = itk::ReadImage<ImageType>(inputImage);
 
   using StructuringElementType = itk::FlatStructuringElement<Dimension>;
   StructuringElementType::RadiusType radius;
@@ -52,19 +50,14 @@ main(int argc, char * argv[])
   using BinaryErodeImageFilterType = itk::BinaryErodeImageFilter<ImageType, ImageType, StructuringElementType>;
 
   BinaryErodeImageFilterType::Pointer erodeFilter = BinaryErodeImageFilterType::New();
-  erodeFilter->SetInput(reader->GetOutput());
+  erodeFilter->SetInput(input);
   erodeFilter->SetKernel(structuringElement);
   erodeFilter->SetForegroundValue(255); // Intensity value to erode
   erodeFilter->SetBackgroundValue(0);   // Replacement value for eroded voxels
 
-  using WriterType = itk::ImageFileWriter<ImageType>;
-  WriterType::Pointer writer = WriterType::New();
-  writer->SetInput(erodeFilter->GetOutput());
-  writer->SetFileName(outputImage);
-
   try
   {
-    writer->Update();
+    itk::WriteImage(erodeFilter->GetOutput(), outputImage);
   }
   catch (itk::ExceptionObject & e)
   {