blob: 99371d562c0f5dcf6ea42b1ef7fe6860c89486ef (
plain)
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
|
/* This file is part of the Gudhi Library - https://gudhi.inria.fr/ - which is released under MIT.
* See file LICENSE or go to https://gudhi.inria.fr/licensing/ for full license details.
* Author(s): Vincent Rouvreau
*
* Copyright (C) 2015 Inria
*
* Modification(s):
* - YYYY/MM Author: Description of the modification
*/
#ifndef POINTS_OFF_IO_H_
#define POINTS_OFF_IO_H_
#include <gudhi/Off_reader.h>
#include <string>
#include <vector>
#include <fstream>
#include <map>
namespace Gudhi {
/**
* \brief OFF file visitor implementation according to Off_reader in order to read points from an OFF file.
*/
template<typename Point_d>
class Points_off_visitor_reader {
private:
std::vector<Point_d> point_cloud;
public:
/** \brief Off_reader visitor init implementation.
*
* The init parameters are set from OFF file header.
* Dimension value is required in order to construct a vector of points.
*
* @param[in] dim space dimension of vertices.
* @param[in] num_vertices number of vertices in the OFF file (not used).
* @param[in] num_faces number of faces in the OFF file (not used).
* @param[in] num_edges number of edges in the OFF file (not used).
*/
void init(int dim, int num_vertices, int num_faces, int num_edges) {
#ifdef DEBUG_TRACES
std::cout << "Points_off_visitor_reader::init - dim=" << dim << " - num_vertices=" <<
num_vertices << " - num_faces=" << num_faces << " - num_edges=" << num_edges << std::endl;
#endif // DEBUG_TRACES
if (num_faces > 0) {
std::cerr << "Points_off_visitor_reader::init faces are not taken into account from OFF file for Points.\n";
}
if (num_edges > 0) {
std::cerr << "Points_off_visitor_reader::init edges are not taken into account from OFF file for Points.\n";
}
}
/** @brief Off_reader visitor point implementation.
*
* The point function is called on each vertex line from OFF file.
* This function inserts the vertex in the vector of points.
*
* @param[in] point vector of vertex coordinates.
*
* @details
* Point_d must have a constructor with the following form:
*
* @code template<class InputIterator > Point_d::Point_d(InputIterator first, InputIterator last) @endcode
*
*/
void point(const std::vector<double>& point) {
#ifdef DEBUG_TRACES
std::cout << "Points_off_visitor_reader::point ";
for (auto coordinate : point) {
std::cout << coordinate << " | ";
}
std::cout << std::endl;
#endif // DEBUG_TRACES
// Fill the point cloud
point_cloud.push_back(Point_d(point.begin(), point.end()));
}
// Off_reader visitor maximal_face implementation - Only points are read
void maximal_face(const std::vector<int>& face) { }
// Off_reader visitor done implementation - Only points are read
void done() { }
/** \brief Point cloud getter.
*
* @return point_cloud.
*/
const std::vector<Point_d>& get_point_cloud() const {
return point_cloud;
}
};
/**
* \brief OFF file reader implementation in order to read points from an OFF file.
*
* This class is using the Points_off_visitor_reader to visit the OFF file according to Off_reader.
*
* Point_d must have a constructor with the following form:
*
* \code template<class InputIterator > Point_d::Point_d(int d, InputIterator first, InputIterator last) \endcode
*
* where d is the point dimension.
*
* \section pointoffioexample Example
*
* This example loads points from an OFF file and builds a vector of points (vector of double).
* Then, it is asked to display the points.
*
* \include common/example_vector_double_points_off_reader.cpp
*
* When launching:
*
* \code $> ./vector_double_off_reader ../../data/points/alphacomplexdoc.off
* \endcode
*
* the program outputs a file ../../data/points/alphacomplexdoc.off.txt:
*
* \include common/vectordoubleoffreader_result.txt
*/
template<typename Point_d>
class Points_off_reader {
public:
/** \brief Reads the OFF file and constructs a vector of points from the points
* that are in the OFF file.
*
* @param[in] name_file OFF file to read.
*
* \post Check with is_valid() function to see if read operation was successful.
*/
Points_off_reader(const std::string& name_file)
: valid_(false) {
std::ifstream stream(name_file);
if (stream.is_open()) {
Off_reader off_reader(stream);
Points_off_visitor_reader<Point_d> off_visitor;
valid_ = off_reader.read(off_visitor);
if (valid_) {
point_cloud = off_visitor.get_point_cloud();
}
} else {
std::cerr << "Points_off_reader::Points_off_reader could not open file " << name_file << "\n";
}
}
/** \brief Returns if the OFF file read operation was successful or not.
*
* @return OFF file read status.
*/
bool is_valid() const {
return valid_;
}
/** \brief Point cloud getter.
*
* @return point_cloud.
*/
const std::vector<Point_d>& get_point_cloud() const {
return point_cloud;
}
private:
/** \brief point_cloud.*/
std::vector<Point_d> point_cloud;
/** \brief OFF file read status.*/
bool valid_;
};
} // namespace Gudhi
#endif // POINTS_OFF_IO_H_
|
