implement converting paths...
Stefan Schuermans authored 11 years ago
|
1) /* drawing (DXF) to G-code (NGC) converter
2) * Copyright 2013 Stefan Schuermans <stefan@schuermans.info>
3) * Copyleft: CC-BY-SA http://creativecommons.org/licenses/by-sa/3.0/
4) */
5)
6) #include <boost/shared_ptr.hpp>
7) #include <iostream>
8) #include <vector>
9)
|
implement converting paths...
Stefan Schuermans authored 11 years ago
|
15)
16) #include "layer.h"
17) #include "path.h"
18) #include "point.h"
19) #include "polygons.h"
20)
21) /**
22) * @brief clear all polygons
23) */
24) void Polygons::clear()
25) {
26) mPolys.clear();
27) }
28)
29) /**
30) * @brief add polygons from layer
31) * @param[in] layer layer to obtain the polygons from
32) * @param[in] eqDist maximum distance of two points to be considered equal
33) * @return if the layer could be converted to polygons and imported
34) */
35) bool Polygons::addLayer(const Layer &layer, double eqDist)
36) {
37) // convert all paths to simple polygons
38) CgPolyVec simples;
39) Layer::Paths::const_iterator path;
40) for (path = layer.mPaths.begin(); path != layer.mPaths.end(); ++path) {
41) // skip empty paths
42) if (path->mPoints.empty())
43) continue;
44) // check that path is closed
45) if (!path->mPoints.front().equals(path->mPoints.back(), eqDist)) {
46) std::cerr << "path not closed" << std::endl;
47) return false;
48) }
49) // create simple polygon from path
50) // - do not add last point to polygon, as it is the same as the first one
51) CgPoly simple;
52) Path::Points::const_iterator pt;
53) for (pt = path->mPoints.begin(); pt + 1 != path->mPoints.end(); ++pt)
54) simple.push_back(CgPoint(pt->mX, pt->mY));
55) // check that polygon is simple
56) if (!simple.is_simple()) {
57) std::cerr << "path is not simple (maybe self-itersecting?)" << std::endl;
58) return false;
59) }
60) // ensure orientation is clockwise (must be the same for all polygons)
61) if (simple.is_clockwise_oriented())
62) simple.reverse_orientation();
63) // collect polygons
64) simples.push_back(simple);
65) } // for path
66)
67) // check that no partly overlaps are present
68) // deterine which polygon is contained in which one
69) std::vector<ssize_t> contained_in; // idx: contained poly, val: containing poly
70) ssize_t ia, ib;
71) CgPolyVec::const_iterator a, b;
72) for (a = simples.begin(); a != simples.end(); ++a)
73) contained_in.push_back(-1); // not contained in any polygon
74) for (a = simples.begin(), ia = 0; a != simples.end(); ++a, ++ia) {
75) b = a;
76) ib = ia;
77) for (++b, ++ib; b != simples.end(); ++b, ++ib) {
78) // A and B intersect each other
79) if (CGAL::do_intersect(*a, *b)) {
80) // compute differences A - B and B - A
81) CgPolyHolesVec a_b, b_a;
82) CGAL::difference(*a, *b, std::back_inserter(a_b));
83) CGAL::difference(*b, *a, std::back_inserter(b_a));
84) if (a_b.empty()) {
85) // A and B are identical
86) if (b_a.empty()) {
87) std::cerr << "duplicate polygons found - meaning is undefined"
88) << std::endl;
89) return false;
90) }
91) // B contains A (A is hole in B)
92) else {
93) if (contained_in.at(ia) >= 0) {
94) std::cerr << "polygon contained in multiple polygons"
95) " - meaning is undefined" << std::endl;
96) return false;
97) }
98) contained_in.at(ia) = ib;
99) }
100) } else {
101) // A contains B (B is hole in A)
102) if (b_a.empty()) {
103) if (contained_in.at(ib) >= 0) {
104) std::cerr << "polygon contained in multiple polygons"
105) " - meaning is undefined" << std::endl;
106) return false;
107) }
108) contained_in.at(ib) = ia;
109) }
110) // A and B overlap partly
111) else {
112) std::cerr << "polygons overlap partly - meaning is undefined"
113) << std::endl;
114) return false;
115) }
116) }
117) }
118) } // for b
119) } // for a
120)
121) // add found polygons to member variable
122) for (ia = 0; (size_t)ia < contained_in.size(); ++ia) {
123) // outer polygon
124) if (contained_in.at(ia) < 0) {
125) CgPolyHoles poly(simples.at(ia));
126) // find and add holes
127) for (ib = 0; (size_t)ib < contained_in.size(); ++ib)
128) if (contained_in.at(ib) == ia)
129) poly.add_hole(simples.at(ib));
130) // add polygon to member variable
131) mPolys.push_back(poly);
132) }
133) }
134)
135) return true;
136) }
137)
138) /**
139) * @brief load polygons from layer
140) * @param[in] layer layer to obtain the polygons from
141) * @param[in] eqDist maximum distance of two points to be considered equal
142) * @return if the layer could be converted to polygons and imported
143) */
144) bool Polygons::loadLayer(const Layer &layer, double eqDist)
145) {
146) clear();
147) return addLayer(layer, eqDist);
148) }
149)
|
implement inside cutting (h...
Stefan Schuermans authored 11 years ago
|
150) /**
151) * @brief create inner offset polygons
152) * @param[in] offset offset, > 0.0
153) * @param[out] offsetPolys offset polygons (!= *this)
154) */
155) void Polygons::createInnerOffset(double offset, Polygons &offsetPolys) const
156) {
157) // clear output polygons
158) offsetPolys.mPolys.clear();
159) // leave if tool diameter is invalid
160) if (offset <= 0.0)
161) return;
162)
163) // process all polygons
164) CgPolyHolesVec::const_iterator poly;
165) for (poly = mPolys.begin(); poly != mPolys.end(); ++poly) {
166)
167) // create inner offset polygons
168) CgPolyHolesPtrVec offPolys =
169) CGAL::create_interior_skeleton_and_offset_polygons_with_holes_2(
170) offset, *poly);
171)
172) // add offset polygons to output polygons
173) CgPolyHolesPtrVec::const_iterator offPoly;
174) for (offPoly = offPolys.begin(); offPoly != offPolys.end(); ++offPoly)
175) offsetPolys.mPolys.push_back(**offPoly);
176)
177) } // for poly
178) }
179)
|