ITK/Examples/ImageProcessing/MinimumImageFilter

MinimumImageFilter.cxx

<source lang="cpp">

1. include "itkImage.h"
2. include "itkImageFileWriter.h"
3. include "itkMinimumImageFilter.h"
4. include "itkImageFileReader.h"
5. include "itkImageRegionIterator.h"

typedef itk::Image<unsigned char, 2> ImageType;

static void CreateImage1(ImageType* image); static void CreateImage2(ImageType* image);

int main(int, char*[]) {

 ImageType::Pointer image1 = ImageType::New();
 CreateImage1(image1);

 ImageType::Pointer image2 = ImageType::New();
 CreateImage2(image2);

 typedef itk::MinimumImageFilter <ImageType>
   MinimumImageFilterType;

 MinimumImageFilterType::Pointer minimumImageFilter
   = MinimumImageFilterType::New ();
 minimumImageFilter->SetInput(0, image1);
 minimumImageFilter->SetInput(1, image2);
 minimumImageFilter->Update();

 return EXIT_SUCCESS;

}

void CreateImage1(ImageType* image) {

 ImageType::IndexType start;
 start.Fill(0);

 ImageType::SizeType size;
 size.Fill(100);

 ImageType::RegionType region;
 region.SetSize(size);
 region.SetIndex(start);

 image->SetRegions(region);
 image->Allocate();

 itk::ImageRegionIterator<ImageType> imageIterator(image,region);

 while(!imageIterator.IsAtEnd())
   {
   if(imageIterator.GetIndex()[0] < 30)
     {
     imageIterator.Set(255);
     }
   else
     {
     imageIterator.Set(0);
     }

   ++imageIterator;
   }

}

void CreateImage2(ImageType* image) {

 ImageType::IndexType start;
 start.Fill(0);

 ImageType::SizeType size;
 size.Fill(100);

 ImageType::RegionType region;
 region.SetSize(size);
 region.SetIndex(start);

 image->SetRegions(region);
 image->Allocate();

 itk::ImageRegionIterator<ImageType> imageIterator(image,region);

 while(!imageIterator.IsAtEnd())
   {
   if(imageIterator.GetIndex()[0] > 70)
     {
     imageIterator.Set(255);
     }
   else
     {
     imageIterator.Set(0);
     }

   ++imageIterator;
   }

} </source>

CMakeLists.txt

<source lang="cmake"> cmake_minimum_required(VERSION 2.6)

PROJECT(MinimumImageFilter)

FIND_PACKAGE(ITK REQUIRED) INCLUDE(${ITK_USE_FILE})

ADD_EXECUTABLE(MinimumImageFilter MinimumImageFilter.cxx) TARGET_LINK_LIBRARIES(MinimumImageFilter ITKBasicFilters ITKCommon ITKIO) </source>