ITK  4.9.0
Insight Segmentation and Registration Toolkit
WikiExamples/Segmentation/OtsuMultipleThresholdsImageFilter.cxx
#include "itkImage.h"
#include "itksys/SystemTools.hxx"
#include <sstream>
#include "QuickView.h"
typedef unsigned char PixelType;
typedef itk::Image<PixelType, 2> ImageType;
typedef itk::RGBPixel<unsigned char> RGBPixelType;
typedef itk::Image<RGBPixelType, 2> RGBImageType;
static void CreateImage(ImageType::Pointer image);
int main(int argc, char *argv[])
{
unsigned int numberOfThresholds = 1;
ImageType::Pointer image;
if( argc < 2 )
{
image = ImageType::New();
CreateImage(image.GetPointer());
}
else
{
typedef itk::ImageFileReader<ImageType> ReaderType;
ReaderType::Pointer reader =
ReaderType::New();
reader->SetFileName(argv[1]);
reader->Update();
if (argc > 2)
{
numberOfThresholds = atoi(argv[2]);
}
image = reader->GetOutput();
}
FilterType;
FilterType::Pointer otsuFilter
= FilterType::New();
otsuFilter->SetInput(image);
otsuFilter->SetNumberOfThresholds(numberOfThresholds);
otsuFilter->Update(); // To compute threshold
RGBFilterType::Pointer colormapImageFilter = RGBFilterType::New();
colormapImageFilter->SetInput(otsuFilter->GetOutput());
QuickView viewer;
viewer.AddImage(
image.GetPointer(),true,
argc > 1 ? itksys::SystemTools::GetFilenameName(argv[1]) : "Generated image");
std::stringstream desc;
desc << "Otsu Threshold\n"
<< otsuFilter->GetNumberOfThresholds()
<< " thresholds: ";
FilterType::ThresholdVectorType thresholds = otsuFilter->GetThresholds();
for (unsigned int i = 0; i < thresholds.size(); i++)
{
desc << itk::NumericTraits<FilterType::InputPixelType>::PrintType(thresholds[i]) << " ";
}
viewer.AddRGBImage(
colormapImageFilter->GetOutput(),
true,
desc.str());
viewer.Visualize();
return EXIT_SUCCESS;
}
void CreateImage(ImageType::Pointer image)
{
// Create an 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();
// Make the whole image white
itk::ImageRegionIterator<ImageType> iterator(image,image->GetLargestPossibleRegion());
/*
//Create a square
while(!iterator.IsAtEnd())
{
iterator.Set(255);
++iterator;
}
*/
}