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
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: $RCSfile: itkSimplexMeshToTriangleMeshFilter.h,v $
Language: C++
Date: $Date: 2005-08-14 21:21:28 $
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.
=========================================================================*/
#ifndef __SimplexMeshToTriangleMeshFilter_h
#define __SimplexMeshToTriangleMeshFilter_h
#include <itkMesh.h>
#include <itkLineCell.h>
#include <itkPolygonCell.h>
#include <itkVertexCell.h>
#include <itkMapContainer.h>
#include "itkSimplexMesh.h"
#include "itkMeshToMeshFilter.h"
#include "itkVectorContainer.h"
#include "itkMapContainer.h"
#include "itkAutomaticTopologyMeshSource.h"
namespace itk
{
/** \class SimplexMeshToTriangleMeshFilter
* \brief This filter converts a 2-simplex mesh into a triangle mesh
*
* Convert a simplex mesh into a triangle mesh. Therefore the center of each
* simplex cell is computed. These centers are taken as the points for the
* triangle mesh then the points are connected.
*
*
* \author Thomas Boettger. Division Medical and Biological Informatics, German Cancer Research Center, Heidelberg.
*
*/
template <class TInputMesh, class TOutputMesh>
class SimplexMeshToTriangleMeshFilter : public MeshToMeshFilter<TInputMesh, TOutputMesh>
{
public:
/** Standard "Self" typedef. */
typedef SimplexMeshToTriangleMeshFilter Self;
/** Standard "Superclass" typedef. */
typedef MeshToMeshFilter<TInputMesh, TOutputMesh> Superclass;
/** Smart pointer typedef support */
typedef SmartPointer<Self> Pointer;
typedef SmartPointer<const Self> ConstPointer;
/** Method of creation through the object factory. */
itkNewMacro(Self);
/** Run-time type information (and related methods). */
itkTypeMacro(SimplexMeshToTriangleMeshFilter,MeshToMeshFilter);
typedef TInputMesh InputMeshType;
typedef typename InputMeshType::Pointer InputMeshPointer;
typedef typename InputMeshType::PointType InputPointType;
typedef typename InputMeshType::PixelType InputPixelType;
typedef typename InputMeshType::MeshTraits::CellTraits InputCellTraitsType;
typedef typename InputMeshType::PointsContainer InputPointsContainer;
typedef typename InputPointsContainer::Pointer InputPointsContainerPointer;
typedef typename InputPointsContainer::Iterator InputPointsContainerIterator;
typedef typename InputMeshType::NeighborListType InputNeighbors;
typedef typename InputMeshType::NeighborListType::iterator InputNeighborsIterator;
typedef itk::AutomaticTopologyMeshSource<TOutputMesh> AutoMeshSourceType;
typedef typename InputMeshType::CellType SimplexCellType;
typedef itk::PolygonCell<SimplexCellType> SimplexPolygonType;
// stores the center for each simplex mesh cell, key is the point id
typedef itk::MapContainer<unsigned long, InputPointType> PointMapType;
typedef typename PointMapType::Pointer PointMapPointer;
/**
* This class provides methods for visiting
* each simplex cell of a simplex mesh
* It computes the center of each visited cell.
*/
class SimplexCellVisitor
{
public:
/**
* default constructor
*/
SimplexCellVisitor()
{
m_CenterMap = PointMapType::New();
}
/**
* \brief visits all polygon cells and compute the cell centers
*/
void Visit(unsigned long cellId, SimplexPolygonType * poly)
{
typedef typename SimplexPolygonType::PointIdIterator PointIdIterator;
PointIdIterator it = poly->PointIdsBegin();
InputPointType center,p;
center.Fill(0);
while ( it != poly->PointIdsEnd() )
{
this->m_Mesh->GetPoint(*it, &p);
center += p.GetVectorFromOrigin();
it++;
}
center[0] /= poly->GetNumberOfPoints();
center[1] /= poly->GetNumberOfPoints();
center[2] /= poly->GetNumberOfPoints();
m_CenterMap->InsertElement(cellId, center);
//std::cout << "cellId: " << cellId << " center = " << center << std::endl;
}
PointMapPointer GetCenterMap()
{
return m_CenterMap;
}
void SetMesh(InputMeshPointer mesh)
{
this->m_Mesh = mesh;
}
protected:
InputMeshPointer m_Mesh;
PointMapPointer m_CenterMap;
};
typedef itk::CellInterfaceVisitorImplementation<InputPixelType,
InputCellTraitsType,
SimplexPolygonType,
SimplexCellVisitor>
SimplexVisitorInterfaceType;
typedef typename SimplexVisitorInterfaceType::Pointer SimplexVisitorInterfacePointer;
typedef typename SimplexCellType::MultiVisitor CellMultiVisitorType;
typedef typename CellMultiVisitorType::Pointer CellMultiVisitorPointer;
protected:
SimplexMeshToTriangleMeshFilter();
~SimplexMeshToTriangleMeshFilter();
SimplexMeshToTriangleMeshFilter(const Self&) {}
void operator=(const Self&) {}
void PrintSelf(std::ostream& os, Indent indent) const;
virtual void GenerateData();
void Initialize();
/** creates dual triangles for all simplex cells */
void CreateTriangles();
/** part of algorithm */
unsigned long FindCellId(unsigned long id1, unsigned long id2, unsigned long id3);
/** attribute stores the result of the simplex cell visitor */
PointMapPointer m_Centers;
};
} //end of namespace
#ifndef ITK_MANUAL_INSTANTIATION
#include "itkSimplexMeshToTriangleMeshFilter.txx"
#endif
#endif //__SimplexMeshToTriangleMeshFilter_h
|
