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
|
#include "vtkF3DExampleReader.h"
//----------------------------------------------------------------------------
vtkStandardNewMacro(vtkF3DExampleReader);
//----------------------------------------------------------------------------
vtkF3DExampleReader::vtkF3DExampleReader()
{
// All the readers have no input, this have to be specified to the VTK pipeline
this->SetNumberOfInputPorts(0);
}
//----------------------------------------------------------------------------
vtkF3DExampleReader::~vtkF3DExampleReader() = default;
//----------------------------------------------------------------------------
void vtkF3DExampleReader::PrintSelf(ostream& os, vtkIndent indent)
{
this->Superclass::PrintSelf(os, indent);
os << indent << "FileName: " << (this->FileName.empty() ? "(none)" : this->FileName) << "\n";
}
//----------------------------------------------------------------------------
int vtkF3DExampleReader::RequestData(
vtkInformation*, vtkInformationVector**, vtkInformationVector* outputVector)
{
// vtkPolyData is the VTK class of surfacic mesh
// Other data structure exists for volumetric mesh like vtkImageData or vtkUnstructuredGrid
vtkPolyData* output = vtkPolyData::GetData(outputVector);
std::ifstream f(this->FileName);
if (f.is_open())
{
std::vector<std::string> tokens;
std::copy(std::istream_iterator<std::string>(f), std::istream_iterator<std::string>(),
std::back_inserter(tokens));
if (tokens.size() % 3 != 0)
{
vtkErrorMacro("Found " << tokens.size() << " numbers, which is not a multiple of 3.");
return 0;
}
// Since this example is only for reading point clouds, we just have to fill a vtkPoints
// In case you need polygons too, you would have to fill a vtkCellArray too
// and give it to vtkPolyData::SetPolys
vtkNew<vtkPoints> points;
for (size_t i = 0; i < tokens.size(); i += 3)
{
// This is a simple call to add a new 3D point given 3 floating point values
points->InsertNextPoint(std::atof(tokens[i].c_str()), std::atof(tokens[i + 1].c_str()),
std::atof(tokens[i + 2].c_str()));
}
// Attach the points to the polydata
output->SetPoints(points);
f.close();
}
else
{
vtkErrorMacro("File " << this->FileName << " cannot be found.");
return 0;
}
return 1;
}
|
