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
|
//##########################################################################
//# #
//# CLOUDCOMPARE #
//# #
//# This program is free software; you can redistribute it and/or modify #
//# it under the terms of the GNU General Public License as published by #
//# the Free Software Foundation; version 2 or later of the License. #
//# #
//# This program is distributed in the hope that it will be useful, #
//# but WITHOUT ANY WARRANTY; without even the implied warranty of #
//# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the #
//# GNU General Public License for more details. #
//# #
//# COPYRIGHT: EDF R&D / TELECOM ParisTech (ENST-TSI) #
//# #
//##########################################################################
#include "HeightProfileFilter.h"
//qCC_db
#include <ccPolyline.h>
//Qt
#include <QFile>
#include <QTextStream>
bool HeightProfileFilter::canSave(CC_CLASS_ENUM type, bool& multiple, bool& exclusive) const
{
if (type == CC_TYPES::POLY_LINE)
{
multiple = false;
exclusive = true;
return true;
}
return false;
}
CC_FILE_ERROR HeightProfileFilter::saveToFile(ccHObject* entity, const QString& filename, const SaveParameters& parameters)
{
if (!entity || filename.isEmpty())
{
return CC_FERR_BAD_ARGUMENT;
}
//get the polyline
if (!entity->isA(CC_TYPES::POLY_LINE))
{
return CC_FERR_BAD_ENTITY_TYPE;
}
ccPolyline* poly = static_cast<ccPolyline*>(entity);
unsigned vertCount = poly->size();
if (vertCount == 0)
{
//invalid size
ccLog::Warning(QString("[Height profile] Polyline '%1' is empty").arg(poly->getName()));
return CC_FERR_NO_SAVE;
}
//open ASCII file for writing
QFile file(filename);
if (!file.open(QIODevice::WriteOnly | QIODevice::Text))
{
return CC_FERR_WRITING;
}
QTextStream outFile(&file);
outFile.setRealNumberNotation(QTextStream::FixedNotation);
outFile.setRealNumberPrecision(sizeof(PointCoordinateType) == 4 && !poly->isShifted() ? 8 : 12);
outFile << "Curvilinear abscissa; Z" << endl;
//curvilinear abscissa
double s = 0;
const CCVector3* lastP = 0;
for (unsigned j = 0; j < vertCount; ++j)
{
const CCVector3* P = poly->getPoint(j);
//update the curvilinear abscissa
if (lastP)
{
s += (*P - *lastP).normd();
}
lastP = P;
//convert to 'local' coordinate system
CCVector3d Pg = poly->toGlobal3d(*P);
outFile << s << "; " << Pg.z << endl;
}
file.close();
return CC_FERR_NO_ERROR;
}
|
