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
106
107
108
109
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: $RCSfile: itkOtsuMultipleThresholdsCalculatorTest.cxx,v $
Language: C++
Date: $Date: 2005-11-25 15:52:39 $
Version: $Revision: 1.6 $
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.
=========================================================================*/
#include "itkWin32Header.h"
#if defined(_MSC_VER)
#pragma warning ( disable : 4786 )
#endif
#include "itkHistogram.h"
#include "itkOtsuMultipleThresholdsCalculator.h"
int itkOtsuMultipleThresholdsCalculatorTest(int, char*[])
{
typedef float MeasurementType ;
typedef itk::Statistics::Histogram< MeasurementType, 1 > HistogramType ;
HistogramType::Pointer histogram = HistogramType::New() ;
// initialize histogram
HistogramType::SizeType size;
size.Fill(64) ;
HistogramType::MeasurementVectorType lowerBound ;
HistogramType::MeasurementVectorType upperBound ;
lowerBound[0] = 0.0 ;
upperBound[0] = 64.0 ;
histogram->Initialize(size, lowerBound, upperBound ) ;
// create vector of values.
typedef std::vector<MeasurementType> ValuesVectorType;
ValuesVectorType values;
values.push_back(8.0);
values.push_back(16.0);
values.push_back(32.0);
values.push_back(48.0);
MeasurementType range = 2.0;
// create histogram with samples at values +- range.
for (HistogramType::Iterator iter = histogram->Begin(); iter != histogram->End(); ++iter)
{
HistogramType::MeasurementType measurement = iter.GetMeasurementVector()[0];
for (ValuesVectorType::const_iterator viter = values.begin(); viter != values.end(); ++viter)
{
if (measurement > (*viter-range) && measurement < (*viter+range))
{
iter.SetFrequency(1.0);
}
}
}
// Compute numberOfValues - 1 thresholds.
unsigned long numberOfThresholds = values.size() - 1;
typedef itk::OtsuMultipleThresholdsCalculator<HistogramType> OtsuMultipleThresholdCalculatorType;
OtsuMultipleThresholdCalculatorType::Pointer otsuThresholdCalculator = OtsuMultipleThresholdCalculatorType::New();
otsuThresholdCalculator->SetInputHistogram(histogram.GetPointer());
otsuThresholdCalculator->SetNumberOfThresholds(numberOfThresholds);
try
{
otsuThresholdCalculator->Update();
}
catch(itk::ExceptionObject & exp)
{
std::cerr << exp << std::endl;
}
otsuThresholdCalculator->Print (std::cout);
OtsuMultipleThresholdCalculatorType::OutputType otsuThresholds = otsuThresholdCalculator->GetOutput();
bool passed = true;
// Check if thresholds correctly separate values.
for (unsigned long j = 0; j<numberOfThresholds; ++j)
{
if (otsuThresholds[j] < values[j] || otsuThresholds[j] > values[j+1])
{
passed = false;
break;
}
}
if (!passed)
{
std::cout << "Test failed." << std::endl;
std::cout << otsuThresholdCalculator << std::endl;
return EXIT_FAILURE;
}
std::cout << "Test passed." << std::endl;
return EXIT_SUCCESS;
}
|
