ITK  4.13.0
Insight Segmentation and Registration Toolkit
WikiExamples/Statistics/ImageToHistogramFilterRGB.cxx
#include "itkRGBPixel.h"
#include "itkImage.h"
typedef itk::RGBPixel<unsigned char> RGBPixelType;
typedef itk::Image< RGBPixelType, 2> ImageType;
void CreateImage(ImageType::Pointer image);
int main(int, char *[])
{
const unsigned int MeasurementVectorSize = 3; // RGB
const unsigned int binsPerDimension = 4;
ImageType::Pointer image = ImageType::New();
CreateImage(image);
typedef itk::Statistics::ImageToHistogramFilter< ImageType > ImageToHistogramFilterType;
ImageToHistogramFilterType::HistogramType::MeasurementVectorType lowerBound(MeasurementVectorSize);
lowerBound.Fill(0);
ImageToHistogramFilterType::HistogramType::MeasurementVectorType upperBound(MeasurementVectorSize);
upperBound.Fill(255) ;
size.Fill(binsPerDimension);
ImageToHistogramFilterType::Pointer imageToHistogramFilter = ImageToHistogramFilterType::New();
imageToHistogramFilter->SetInput(image);
imageToHistogramFilter->SetHistogramBinMinimum(lowerBound);
imageToHistogramFilter->SetHistogramBinMaximum(upperBound);
imageToHistogramFilter->SetHistogramSize(size);
imageToHistogramFilter->SetAutoMinimumMaximum(false);
imageToHistogramFilter->Update();
ImageToHistogramFilterType::HistogramType* histogram = imageToHistogramFilter->GetOutput();
unsigned int totalBins = 1;
for(unsigned int i = 0; i < size.GetNumberOfElements(); ++i)
{
totalBins *= histogram->GetSize()[i];
}
std::cout << "totalBins: " << totalBins << std::endl;
std::cout << "Frequency = ";
for(unsigned int i = 0; i < totalBins; ++i)
{
std::cout << histogram->GetFrequency(i) << " ";
}
std::cout << std::endl;
return EXIT_SUCCESS;
}
void CreateImage(ImageType::Pointer image)
{
// Create a black image with a red square and a green square.
// This should produce a histogram with very strong spikes.
size.Fill(10);
start.Fill(0);
itk::ImageRegion<2> region(start, size);
image->SetRegions(region);
image->Allocate();
RGBPixelType zeroPixel;
zeroPixel.SetRed(0);
zeroPixel.SetGreen(0);
zeroPixel.SetBlue(0);
image->FillBuffer(zeroPixel);
RGBPixelType pixel;
pixel.SetRed(110);
pixel.SetGreen(110);
pixel.SetBlue(110);
itk::ImageRandomIteratorWithIndex<ImageType> imageIterator(image, image->GetLargestPossibleRegion());
imageIterator.SetNumberOfSamples(10);
while(!imageIterator.IsAtEnd())
{
std::cout << "Setting pixel " << imageIterator.GetIndex() << " to " << pixel << std::endl;
imageIterator.Set(pixel);
++imageIterator;
}
}