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
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: PointSetToAdaptor.cxx
Language: C++
Date: $Date$
Version: $Revision$
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
// Software Guide : BeginLatex
//
// We will describe how to use \doxygen{PointSet} as a \code{Sample}
// using an adaptor in this example.
//
// \index{itk::Sample!PointSetToListAdaptor}
//
// \subdoxygen{Statistics}{PointSetToListAdaptor} class requires the type
// of input \doxygen{PointSet} object. The \doxygen{PointSet} class is an
// associative data container. Each point in a \code{PointSet} object can have
// its associated data value (optional). For the statistics subsystem, current
// implementation of \code{PointSetToListAdaptor} takes only the point part
// into consideration. In other words, the measurement vectors from a
// \code{PointSetToListAdaptor} object are points from the \code{PointSet}
// object that is plugged-into the adaptor object.
//
// To use, an \doxygen{PointSetToListAdaptor} object, we include the
// header file for the class.
//
// Software Guide : EndLatex
// Software Guide : BeginCodeSnippet
#include "itkPointSetToListAdaptor.h"
// Software Guide : EndCodeSnippet
// Software Guide : BeginLatex
// Since, we are using an adaptor, we also include the header file for
// the \doxygen{PointSet} class.
// Software Guide :EndLatex
// Software Guide : BeginCodeSnippet
#include "itkPointSet.h"
// Software Guide : EndCodeSnippet
int main()
{
// Software Guide : BeginLatex
//
// We assume you already know how to create an \doxygen{PointSet} object. The
// following code snippet will create a 2D image of float pixels filled
// with random values.
//
// Software Guide :EndLatex
// Software Guide : BeginCodeSnippet
typedef itk::PointSet<float,2> FloatPointSet2DType ;
itk::RandomPointSetSource<FloatPointSet2DType>::Pointer random ;
random = itk::RandomPointSetSource<FloatPointSet2DType>::New() ;
random->SetMin(0.0) ;
random->SetMax(1000.0) ;
unsigned long size[2] = {20, 20} ;
random->SetSize(size) ;
float spacing[2] = {0.7, 2.1} ;
random->SetSpacing( spacing ) ;
float origin[2] = {15, 400} ;
random->SetOrigin( origin ) ;
// Software Guide : EndCodeSnippet
// Software Guide : BeginLatex
// We now have an \doxygen{PointSet} object and need to cast it to an
// \doxygen{PointSet} object with array type (anything derived from
// the \doxygen{FixedArray} class) pixels.
//
// Since, the \doxygen{PointSet} object's pixel type is \code{float},
// We will use single element \code{float} \doxygen{FixedArray}
// as our measurement vector type. And that will also be our pixel
// type for the cast filter.
// Software Guide : EndLatex
// Software Guide : BeginCodeSnippet
typedef itk::FixedArray< float, 1 > MeasurementVectorType ;
typedef itk::PointSet< MeasurementVectorType, 2 > ArrayPointSetType ;
typedef itk::ScalarToArrayCastPointSetFilter< FloatPointSet2DType,
ArrayPointSetType > CasterType ;
CasterType::Pointer caster = CasterType::New() ;
caster->SetInput( random->GetOutput() ) ;
caster->Update() ;
// Software Guide : EndCodeSnippet
// Software Guide : BeginLatex
// Up to now, we spend most of time to prepare an \doxygen{PointSet} object
// suitable for the adaptor. Actually, the hard part of this example is
// done. Now, we must define an adaptor with the image type and
// instantiate an object.
// Software Guide : EndLatex
// Software Guide : BeginCodeSnippet
typedef itk::Statistics::PointSetToListAdaptor< ArrayPointSetType > SampleType ;
SampleType::Pointer sample = SampleType::New() ;
// Software Guide : EndCodeSnippet
// Software Guide : BeginLatex
// The final thing we have to is to plug-in the image object to the adaptor.
// After that, we can use the common methods and iterator interfaces
// shown in \ref{sec:SampleInterface}.
// Software Guide : EndLatex
// Software Guide : BeginCodeSnippet
sample->SetPointSet( caster->GetOutput() ) ;
// Software Guide : EndCodeSnippet
return 0 ;
}
|
