itkMorphologicalContourInterpolationTest.cxx

/*=========================================================================
 *
 *  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 "itkImageFileReader.h"
#include "itkImageFileWriter.h"
#include "itkMorphologicalContourInterpolator.h"
#include <cstdlib>
#include <string>

template< typename ImageType >
void
doTest( std::string inFilename, std::string outFilename, bool UseDistanceTransform, bool ball, int axis, int label )
{
  using ReaderType = itk::ImageFileReader< ImageType >;
  typename ReaderType::Pointer reader = ReaderType::New();
  reader->SetFileName( inFilename );
  reader->Update();

  typename ImageType::Pointer test = reader->GetOutput();

  // region for partial coverage
  typename ImageType::RegionType reg = test->GetLargestPossibleRegion();
  // for (int i = 0; i < ImageType::ImageDimension; i++)
  //  {
  //  reg.GetModifiableIndex()[i] += (reg.GetSize(i) - 1) / 4;
  //  reg.SetSize(i, (reg.GetSize(i) + 1) / 2);
  //  }

  using mciType = itk::MorphologicalContourInterpolator< ImageType >;
  typename mciType::Pointer mci = mciType::New();
  mci->SetInput( test );
  mci->SetUseDistanceTransform( UseDistanceTransform );
  mci->SetUseBallStructuringElement( ball );
  mci->SetAxis( axis );
  mci->SetLabel( label );

  using WriterType = itk::ImageFileWriter< ImageType >;
  typename WriterType::Pointer writer = WriterType::New();
  writer->SetFileName( outFilename );
  writer->SetInput( mci->GetOutput() );
  writer->SetUseCompression( true );
  writer->Update();
}

int
itkMorphologicalContourInterpolationTest( int argc, char* argv[] )
{
  if ( argc < 3 )
    {
      std::cerr << "Usage: " << argv[0];
      std::cerr << " inputImage outputImage [algorithm] [axis] [label]\n";
      std::cerr << " algorithms:\n";
      std::cerr << "  B = repeated dilations with ball structuring element";
      std::cerr << "  C = repeated dilations with cross structuring element";
      std::cerr << "  T = distance transform (not geodesic!)";
      std::cerr << " defaults: algo B, axis -1 (all axes), label 0 (all labels)";
      std::cerr << std::endl;
      return EXIT_FAILURE;
    }
  const char* inputImageFileName = argv[1];
  const char* outputImageFileName = argv[2];
  bool        dt = false; // DistanceTransform
  bool        ball = true;
  int         axis = -1, label = 0;
  if ( argc >= 4 )
    {
      char algo = toupper( argv[3][0] );
      if ( algo == 'T' )
        {
          dt = true;
        }
      else if ( algo == 'C' )
        {
          ball = false;
        }
      // else B
    }
  if ( argc >= 5 )
    {
      axis = strtol( argv[4], nullptr, 10 );
    }
  if ( argc >= 6 )
    {
      label = strtol( argv[5], nullptr, 10 );
    }

  using ScalarPixelType = itk::CommonEnums::IOComponent;
  itk::ImageIOBase::Pointer imageIO = itk::ImageIOFactory::CreateImageIO( inputImageFileName, itk::CommonEnums::IOFileMode::ReadMode );
  if ( !imageIO )
    {
      std::cerr << "Could not CreateImageIO for: " << inputImageFileName << std::endl;
      return EXIT_FAILURE;
    }
  imageIO->SetFileName( inputImageFileName );
  imageIO->ReadImageInformation();
  const ScalarPixelType pixelType = imageIO->GetComponentType();
  const size_t          numDimensions = imageIO->GetNumberOfDimensions();

  try
    {
      // unused cases are not instantiated because they greatly increase compile time
      if ( numDimensions == 3 && ( pixelType == itk::CommonEnums::IOComponent::SHORT || pixelType == itk::CommonEnums::IOComponent::USHORT ) )
        {
          doTest< itk::Image< short, 3 > >( inputImageFileName, outputImageFileName, dt, ball, axis, label );
          return EXIT_SUCCESS;
        }
      if ( numDimensions == 4 && pixelType == itk::CommonEnums::IOComponent::UCHAR )
        {
          doTest< itk::Image< unsigned char, 4 > >( inputImageFileName, outputImageFileName, dt, ball, axis, label );
          return EXIT_SUCCESS;
        }

      std::cerr << "Unsupported image type:\n  Dimensions: " << numDimensions;
      std::cerr << "\n  Pixel type:" << imageIO->GetComponentTypeAsString( pixelType ) << std::endl;
      return EXIT_FAILURE;
    }
  catch ( itk::ExceptionObject& error )
    {
      std::cerr << "Error: " << error << std::endl;
      return EXIT_FAILURE;
    }
  return EXIT_SUCCESS;
}