1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: $RCSfile: itkScalarImageToHistogramGeneratorTest.cxx,v $
Language: C++
Date: $Date: 2005-02-08 03:18:41 $
Version: $Revision: 1.2 $
Copyright (c) Insight Software Consortium. All rights reserved.
See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notices for more information.
=========================================================================*/
#if defined(_MSC_VER)
#pragma warning ( disable : 4786 )
#endif
#include "itkScalarImageToHistogramGenerator.h"
#include "itkImage.h"
#include "itkImageRegionIteratorWithIndex.h"
int itkScalarImageToHistogramGeneratorTest( int , char * [] )
{
typedef unsigned char PixelType;
const unsigned int Dimension = 2;
typedef itk::Image< PixelType, Dimension > ImageType;
ImageType::Pointer image = ImageType::New();
ImageType::RegionType region;
ImageType::SizeType size;
ImageType::IndexType start;
size[0] = 255;
size[1] = 255;
start[0] = 0;
start[1] = 0;
region.SetIndex( start );
region.SetSize( size );
image->SetRegions( region );
image->Allocate();
// Now fill up the image will all the combinations of
// values from 0-255.
itk::ImageRegionIteratorWithIndex< ImageType > it( image, region );
it.GoToBegin();
while( !it.IsAtEnd() )
{
it.Set( it.GetIndex()[0] );
++it;
}
typedef itk::Statistics::ScalarImageToHistogramGenerator<
ImageType
> HistogramGeneratorType;
HistogramGeneratorType::Pointer histogramGenerator = HistogramGeneratorType::New();
histogramGenerator->SetInput( image );
histogramGenerator->SetNumberOfBins( 255 );
histogramGenerator->SetMarginalScale( 10.0 );
histogramGenerator->Compute();
typedef HistogramGeneratorType::HistogramType HistogramType;
const HistogramType * histogram = histogramGenerator->GetOutput();
const unsigned int histogramSize = histogram->Size();
std::cout << "Histogram size " << histogramSize << std::endl;
unsigned int channel = 0; // the only channel
for( unsigned int bin=0; bin < histogramSize; bin++ )
{
if( histogram->GetFrequency( bin, channel ) != 255 )
{
std::cerr << "Error in bin= " << bin << " channel = " << channel << std::endl;
std::cerr << "Frequency was= " << histogram->GetFrequency( bin, channel ) << " Instead of the expected " << 255 << std::endl;
return EXIT_FAILURE;
}
}
return EXIT_SUCCESS;
}
|
