1
#include "MultiPoint.hpp"
2
#include "BoundingBox.hpp"
6
MultiPoint::operator Points() const
12
MultiPoint::scale(double factor)
14
for (Points::iterator it = points.begin(); it != points.end(); ++it) {
20
MultiPoint::translate(double x, double y)
22
for (Points::iterator it = points.begin(); it != points.end(); ++it) {
23
(*it).translate(x, y);
28
MultiPoint::rotate(double angle, const Point ¢er)
30
for (Points::iterator it = points.begin(); it != points.end(); ++it) {
31
(*it).rotate(angle, center);
38
std::reverse(this->points.begin(), this->points.end());
42
MultiPoint::first_point() const
44
return this->points.front();
48
MultiPoint::length() const
50
Lines lines = this->lines();
52
for (Lines::iterator it = lines.begin(); it != lines.end(); ++it) {
59
MultiPoint::is_valid() const
61
return this->points.size() >= 2;
65
MultiPoint::find_point(const Point &point) const
67
for (Points::const_iterator it = this->points.begin(); it != this->points.end(); ++it) {
68
if (it->coincides_with(point)) return it - this->points.begin();
70
return -1; // not found
74
MultiPoint::bounding_box(BoundingBox* bb) const
76
*bb = BoundingBox(this->points);
80
MultiPoint::_douglas_peucker(const Points &points, const double tolerance)
85
Line full(points.front(), points.back());
86
for (Points::const_iterator it = points.begin() + 1; it != points.end(); ++it) {
87
double d = it->distance_to(full);
89
index = it - points.begin();
93
if (dmax >= tolerance) {
95
dp0.reserve(index + 1);
96
dp0.insert(dp0.end(), points.begin(), points.begin() + index + 1);
97
Points dp1 = MultiPoint::_douglas_peucker(dp0, tolerance);
98
results.reserve(results.size() + dp1.size() - 1);
99
results.insert(results.end(), dp1.begin(), dp1.end() - 1);
102
dp0.reserve(points.size() - index + 1);
103
dp0.insert(dp0.end(), points.begin() + index, points.end());
104
dp1 = MultiPoint::_douglas_peucker(dp0, tolerance);
105
results.reserve(results.size() + dp1.size());
106
results.insert(results.end(), dp1.begin(), dp1.end());
108
results.push_back(points.front());
109
results.push_back(points.back());
116
MultiPoint::from_SV(SV* poly_sv)
118
AV* poly_av = (AV*)SvRV(poly_sv);
119
const unsigned int num_points = av_len(poly_av)+1;
120
this->points.resize(num_points);
122
for (unsigned int i = 0; i < num_points; i++) {
123
SV** point_sv = av_fetch(poly_av, i, 0);
124
this->points[i].from_SV_check(*point_sv);
129
MultiPoint::from_SV_check(SV* poly_sv)
131
if (sv_isobject(poly_sv) && (SvTYPE(SvRV(poly_sv)) == SVt_PVMG)) {
132
*this = *(MultiPoint*)SvIV((SV*)SvRV( poly_sv ));
134
this->from_SV(poly_sv);
139
MultiPoint::to_AV() {
140
const unsigned int num_points = this->points.size();
142
av_extend(av, num_points-1);
143
for (unsigned int i = 0; i < num_points; i++) {
144
av_store(av, i, perl_to_SV_ref(this->points[i]));
146
return newRV_noinc((SV*)av);
150
MultiPoint::to_SV_pureperl() const {
151
const unsigned int num_points = this->points.size();
153
av_extend(av, num_points-1);
154
for (unsigned int i = 0; i < num_points; i++) {
155
av_store(av, i, this->points[i].to_SV_pureperl());
157
return newRV_noinc((SV*)av);