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
|
// MyReader.cpp
#include "MyReader.hpp"
#include <pdal/util/ProgramArgs.hpp>
namespace pdal
{
static PluginInfo const s_info
{
"readers.myreader",
"My Awesome Reader",
"http://link/to/documentation"
};
CREATE_SHARED_STAGE(MyReader, s_info)
std::string MyReader::getName() const { return s_info.name; }
void MyReader::addArgs(ProgramArgs& args)
{
args.add("z_scale", "Z Scaling", m_scale_z, 1.0);
}
void MyReader::addDimensions(PointLayoutPtr layout)
{
layout->registerDim(Dimension::Id::X);
layout->registerDim(Dimension::Id::Y);
layout->registerDim(Dimension::Id::Z);
layout->registerOrAssignDim("MyData", Dimension::Type::Unsigned64);
}
void MyReader::ready(PointTableRef)
{
SpatialReference ref("EPSG:4385");
setSpatialReference(ref);
}
template <typename T>
T convert(const StringList& s, const std::string& name, size_t fieldno)
{
T output;
bool bConverted = Utils::fromString(s[fieldno], output);
if (!bConverted)
{
std::stringstream oss;
oss << "Unable to convert " << name << ", " << s[fieldno] <<
", to double";
throw pdal_error(oss.str());
}
return output;
}
point_count_t MyReader::read(PointViewPtr view, point_count_t count)
{
PointLayoutPtr layout = view->layout();
PointId nextId = view->size();
PointId idx = m_index;
point_count_t numRead = 0;
m_stream.reset(new ILeStream(m_filename));
size_t HEADERSIZE(1);
size_t skip_lines((std::max)(HEADERSIZE, (size_t)m_index));
size_t line_no(1);
for (std::string line; std::getline(*m_stream->stream(), line); line_no++)
{
if (line_no <= skip_lines)
{
continue;
}
// MyReader format: X::Y::Z::Data
StringList s = Utils::split2(line, ':');
unsigned long u64(0);
if (s.size() != 4)
{
std::stringstream oss;
oss << "Unable to split proper number of fields. Expected 4, got "
<< s.size();
throw pdal_error(oss.str());
}
std::string name("X");
view->setField(Dimension::Id::X, nextId, convert<double>(s, name, 0));
name = "Y";
view->setField(Dimension::Id::Y, nextId, convert<double>(s, name, 1));
name = "Z";
double z = convert<double>(s, name, 2) * m_scale_z;
view->setField(Dimension::Id::Z, nextId, z);
name = "MyData";
view->setField(layout->findProprietaryDim(name),
nextId,
convert<unsigned int>(s, name, 3));
nextId++;
if (m_cb)
m_cb(*view, nextId);
}
m_index = nextId;
numRead = nextId;
return numRead;
}
void MyReader::done(PointTableRef)
{
m_stream.reset();
}
} //namespace pdal
|
