2010-01-05 22:09:01 +03:00
|
|
|
/*
|
2011-01-21 04:21:09 +03:00
|
|
|
* OpenSCAD (www.openscad.org)
|
|
|
|
* Copyright (C) 2009-2011 Clifford Wolf <clifford@clifford.at> and
|
|
|
|
* Marius Kintel <marius@kintel.net>
|
2010-01-05 22:09:01 +03:00
|
|
|
*
|
|
|
|
* 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; either version 2 of the License, or
|
|
|
|
* (at your option) any later version.
|
|
|
|
*
|
2010-02-01 12:34:18 +03:00
|
|
|
* As a special exception, you have permission to link this program
|
|
|
|
* with the CGAL library and distribute executables, as long as you
|
|
|
|
* follow the requirements of the GNU GPL in regard to all of the
|
|
|
|
* software in the executable aside from CGAL.
|
|
|
|
*
|
2010-01-05 22:09:01 +03:00
|
|
|
* 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.
|
|
|
|
*
|
|
|
|
* You should have received a copy of the GNU General Public License
|
|
|
|
* along with this program; if not, write to the Free Software
|
|
|
|
* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
|
|
|
|
*
|
|
|
|
*/
|
|
|
|
|
2010-01-30 07:26:05 +03:00
|
|
|
#include "dxfdata.h"
|
|
|
|
#include "grid.h"
|
2011-09-01 06:09:06 +04:00
|
|
|
#include "CGAL_Nef_polyhedron.h"
|
2010-01-30 07:26:05 +03:00
|
|
|
#include "cgal.h"
|
2012-10-16 05:17:12 +04:00
|
|
|
#include "cgalutils.h"
|
2013-03-11 06:35:30 +04:00
|
|
|
#include <boost/variant.hpp>
|
|
|
|
#include "polyset.h"
|
|
|
|
#include "dxftess.h"
|
|
|
|
#include "CGALEvaluator.h"
|
|
|
|
#include "Tree.h"
|
2010-01-05 22:09:01 +03:00
|
|
|
|
2010-02-01 13:01:05 +03:00
|
|
|
#ifdef ENABLE_CGAL
|
|
|
|
|
2011-09-01 06:09:06 +04:00
|
|
|
DxfData *CGAL_Nef_polyhedron::convertToDxfData() const
|
2010-01-05 22:09:01 +03:00
|
|
|
{
|
2011-09-01 06:09:06 +04:00
|
|
|
assert(this->dim == 2);
|
|
|
|
DxfData *dxfdata = new DxfData();
|
2010-01-14 11:14:52 +03:00
|
|
|
Grid2d<int> grid(GRID_COARSE);
|
2010-01-05 22:09:01 +03:00
|
|
|
|
|
|
|
typedef CGAL_Nef_polyhedron2::Explorer Explorer;
|
|
|
|
typedef Explorer::Face_const_iterator fci_t;
|
|
|
|
typedef Explorer::Halfedge_around_face_const_circulator heafcc_t;
|
2011-09-01 06:09:06 +04:00
|
|
|
Explorer E = this->p2->explorer();
|
2010-01-05 22:09:01 +03:00
|
|
|
|
2011-01-18 01:50:33 +03:00
|
|
|
for (fci_t fit = E.faces_begin(), facesend = E.faces_end(); fit != facesend; ++fit)
|
2010-01-05 22:09:01 +03:00
|
|
|
{
|
|
|
|
heafcc_t fcirc(E.halfedge(fit)), fend(fcirc);
|
2010-01-14 11:14:52 +03:00
|
|
|
int first_point = -1, last_point = -1;
|
2010-01-05 22:09:01 +03:00
|
|
|
CGAL_For_all(fcirc, fend) {
|
|
|
|
if (E.is_standard(E.target(fcirc))) {
|
|
|
|
Explorer::Point ep = E.point(E.target(fcirc));
|
|
|
|
double x = to_double(ep.x()), y = to_double(ep.y());
|
|
|
|
int this_point = -1;
|
|
|
|
if (grid.has(x, y)) {
|
|
|
|
this_point = grid.align(x, y);
|
|
|
|
} else {
|
2011-09-01 06:09:06 +04:00
|
|
|
this_point = grid.align(x, y) = dxfdata->points.size();
|
2011-09-04 00:44:41 +04:00
|
|
|
dxfdata->points.push_back(Vector2d(x, y));
|
2010-01-05 22:09:01 +03:00
|
|
|
}
|
|
|
|
if (first_point < 0) {
|
2011-09-04 00:44:41 +04:00
|
|
|
dxfdata->paths.push_back(DxfData::Path());
|
2010-01-05 22:09:01 +03:00
|
|
|
first_point = this_point;
|
|
|
|
}
|
2010-01-14 11:14:52 +03:00
|
|
|
if (this_point != last_point) {
|
2011-09-04 00:44:41 +04:00
|
|
|
dxfdata->paths.back().indices.push_back(this_point);
|
2010-01-14 11:14:52 +03:00
|
|
|
last_point = this_point;
|
|
|
|
}
|
2010-01-05 22:09:01 +03:00
|
|
|
}
|
|
|
|
}
|
|
|
|
if (first_point >= 0) {
|
2011-09-04 00:44:41 +04:00
|
|
|
dxfdata->paths.back().is_closed = 1;
|
|
|
|
dxfdata->paths.back().indices.push_back(first_point);
|
2010-01-05 22:09:01 +03:00
|
|
|
}
|
|
|
|
}
|
2010-01-07 23:49:27 +03:00
|
|
|
|
2011-09-01 06:09:06 +04:00
|
|
|
dxfdata->fixup_path_direction();
|
|
|
|
return dxfdata;
|
2010-01-05 22:09:01 +03:00
|
|
|
}
|
|
|
|
|
2012-10-16 05:17:12 +04:00
|
|
|
std::string CGAL_Nef_polyhedron::dump() const
|
|
|
|
{
|
|
|
|
if (this->dim==2)
|
2012-10-28 17:56:23 +04:00
|
|
|
return OpenSCAD::dump_svg( *this->p2 );
|
2012-10-16 05:17:12 +04:00
|
|
|
else if (this->dim==3)
|
2012-10-28 17:56:23 +04:00
|
|
|
return OpenSCAD::dump_svg( *this->p3 );
|
2012-10-16 05:17:12 +04:00
|
|
|
else
|
2012-10-21 21:10:06 +04:00
|
|
|
return std::string("Nef Polyhedron with dimension != 2 or 3");
|
2012-02-11 06:15:16 +04:00
|
|
|
}
|
|
|
|
|
2013-03-10 06:28:43 +04:00
|
|
|
|
2013-03-11 06:35:30 +04:00
|
|
|
void CGAL_Nef_polyhedron::transform( const Transform3d &matrix )
|
2013-03-10 06:28:43 +04:00
|
|
|
{
|
2013-03-11 06:35:30 +04:00
|
|
|
if (!this->isNull()) {
|
|
|
|
if (this->dim == 2) {
|
|
|
|
// Unfortunately CGAL provides no transform method for CGAL_Nef_polyhedron2
|
|
|
|
// objects. So we convert in to our internal 2d data format, transform it,
|
|
|
|
// tesselate it and create a new CGAL_Nef_polyhedron2 from it.. What a hack!
|
|
|
|
Eigen::Matrix2f testmat;
|
|
|
|
testmat << matrix(0,0), matrix(0,1), matrix(1,0), matrix(1,1);
|
|
|
|
if (testmat.determinant() == 0) {
|
|
|
|
PRINT("Warning: Scaling a 2D object with 0 - removing object");
|
|
|
|
this->reset();
|
2013-03-11 19:59:52 +04:00
|
|
|
return;
|
2013-03-10 06:28:43 +04:00
|
|
|
}
|
2013-03-11 06:35:30 +04:00
|
|
|
else {
|
|
|
|
CGAL_Aff_transformation2 t(
|
|
|
|
matrix(0,0), matrix(0,1), matrix(0,3),
|
|
|
|
matrix(1,0), matrix(1,1), matrix(1,3), matrix(3,3));
|
|
|
|
|
|
|
|
DxfData *dd = this->convertToDxfData();
|
|
|
|
for (size_t i=0; i < dd->points.size(); i++) {
|
|
|
|
CGAL_Kernel2::Point_2 p = CGAL_Kernel2::Point_2(dd->points[i][0], dd->points[i][1]);
|
|
|
|
p = t.transform(p);
|
|
|
|
dd->points[i][0] = to_double(p.x());
|
|
|
|
dd->points[i][1] = to_double(p.y());
|
|
|
|
}
|
2013-03-10 06:28:43 +04:00
|
|
|
|
2013-03-11 06:35:30 +04:00
|
|
|
PolySet ps;
|
|
|
|
ps.is2d = true;
|
|
|
|
dxf_tesselate(&ps, *dd, 0, true, false, 0);
|
2013-03-10 06:28:43 +04:00
|
|
|
|
2013-03-11 06:35:30 +04:00
|
|
|
Tree nulltree;
|
|
|
|
CGALEvaluator tmpeval(nulltree);
|
|
|
|
CGAL_Nef_polyhedron N = tmpeval.evaluateCGALMesh(ps);
|
2013-03-12 04:40:02 +04:00
|
|
|
if ( N.p2 ) this->p2.reset( new CGAL_Nef_polyhedron2( *N.p2 ) );
|
2013-03-11 06:35:30 +04:00
|
|
|
delete dd;
|
|
|
|
}
|
2013-03-10 06:28:43 +04:00
|
|
|
}
|
2013-03-11 06:35:30 +04:00
|
|
|
else if (this->dim == 3) {
|
|
|
|
if (matrix.matrix().determinant() == 0) {
|
|
|
|
PRINT("Warning: Scaling a 3D object with 0 - removing object");
|
|
|
|
this->reset();
|
|
|
|
}
|
|
|
|
else {
|
|
|
|
CGAL_Aff_transformation t(
|
|
|
|
matrix(0,0), matrix(0,1), matrix(0,2), matrix(0,3),
|
|
|
|
matrix(1,0), matrix(1,1), matrix(1,2), matrix(1,3),
|
|
|
|
matrix(2,0), matrix(2,1), matrix(2,2), matrix(2,3), matrix(3,3));
|
|
|
|
this->p3->transform(t);
|
|
|
|
}
|
|
|
|
}
|
|
|
|
}
|
2013-03-10 06:28:43 +04:00
|
|
|
}
|
|
|
|
|
2010-02-01 13:01:05 +03:00
|
|
|
#endif // ENABLE_CGAL
|