969 lines
32 KiB
C++
969 lines
32 KiB
C++
/***************************************************************************
|
|
* Copyright (c) 2015 Yorik van Havre (yorik@uncreated.net) *
|
|
* *
|
|
* This file is part of the FreeCAD CAx development system. *
|
|
* *
|
|
* This library is free software; you can redistribute it and/or *
|
|
* modify it under the terms of the GNU Library General Public *
|
|
* License as published by the Free Software Foundation; either *
|
|
* version 2 of the License, or (at your option) any later version. *
|
|
* *
|
|
* This library 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 Library General Public License for more details. *
|
|
* *
|
|
* You should have received a copy of the GNU Library General Public *
|
|
* License along with this library; see the file COPYING.LIB. If not, *
|
|
* write to the Free Software Foundation, Inc., 59 Temple Place, *
|
|
* Suite 330, Boston, MA 02111-1307, USA *
|
|
* *
|
|
***************************************************************************/
|
|
|
|
#include "PreCompiled.h"
|
|
|
|
#ifndef _PreComp_
|
|
#endif
|
|
|
|
#include "ImpExpDxf.h"
|
|
|
|
#include <Approx_Curve3d.hxx>
|
|
#include <BRepAdaptor_Curve.hxx>
|
|
#include <BRep_Builder.hxx>
|
|
#include <BRepBuilderAPI_MakeEdge.hxx>
|
|
#include <BRepBuilderAPI_MakeVertex.hxx>
|
|
#include <GCPnts_UniformAbscissa.hxx>
|
|
#include <GeomAPI_PointsToBSpline.hxx>
|
|
#include <GeomAPI_Interpolate.hxx>
|
|
#include <GeomAdaptor_Curve.hxx>
|
|
#include <Geom_BSplineCurve.hxx>
|
|
#include <Geom_BezierCurve.hxx>
|
|
#include <gp_Circ.hxx>
|
|
#include <gp_Elips.hxx>
|
|
#include <gp_Lin.hxx>
|
|
#include <gp_Ax1.hxx>
|
|
#include <gp_Ax2.hxx>
|
|
#include <gp_Pnt.hxx>
|
|
#include <gp_Dir.hxx>
|
|
#include <gp_Vec.hxx>
|
|
#include <Standard_Version.hxx>
|
|
#include <TopoDS.hxx>
|
|
#include <TopoDS_Edge.hxx>
|
|
#include <TopoDS_Vertex.hxx>
|
|
#include <TopoDS_Compound.hxx>
|
|
#include <TopoDS_Shape.hxx>
|
|
#include <TopExp_Explorer.hxx>
|
|
#include <TopoDS_Edge.hxx>
|
|
#include <TColgp_Array1OfPnt.hxx>
|
|
#if OCC_VERSION_HEX < 0x070600
|
|
#include <BRepAdaptor_HCurve.hxx>
|
|
#endif
|
|
|
|
#include <Base/Console.h>
|
|
#include <Base/Parameter.h>
|
|
#include <Base/Matrix.h>
|
|
#include <Base/Vector3D.h>
|
|
#include <Base/Interpreter.h>
|
|
#include <App/Application.h>
|
|
#include <App/Document.h>
|
|
#include <App/Annotation.h>
|
|
#include <Mod/Part/App/PartFeature.h>
|
|
|
|
using namespace Import;
|
|
|
|
#if OCC_VERSION_HEX >= 0x070600
|
|
using BRepAdaptor_HCurve = BRepAdaptor_Curve;
|
|
#endif
|
|
|
|
|
|
//******************************************************************************
|
|
// reading
|
|
ImpExpDxfRead::ImpExpDxfRead(std::string filepath, App::Document *pcDoc) : CDxfRead(filepath.c_str())
|
|
{
|
|
document = pcDoc;
|
|
setOptionSource("User parameter:BaseApp/Preferences/Mod/Draft");
|
|
setOptions();
|
|
}
|
|
|
|
void ImpExpDxfRead::setOptions()
|
|
{
|
|
ParameterGrp::handle hGrp = App::GetApplication().GetParameterGroupByPath(getOptionSource().c_str());
|
|
optionGroupLayers = hGrp->GetBool("groupLayers",false);
|
|
optionImportAnnotations = hGrp->GetBool("dxftext",false);
|
|
optionScaling = hGrp->GetFloat("dxfScaling",1.0);
|
|
}
|
|
|
|
gp_Pnt ImpExpDxfRead::makePoint(const double* p)
|
|
{
|
|
double sp1(p[0]);
|
|
double sp2(p[1]);
|
|
double sp3(p[2]);
|
|
if (optionScaling != 1.0) {
|
|
sp1 = sp1 * optionScaling;
|
|
sp2 = sp2 * optionScaling;
|
|
sp3 = sp3 * optionScaling;
|
|
}
|
|
return gp_Pnt(sp1,sp2,sp3);
|
|
}
|
|
|
|
void ImpExpDxfRead::OnReadLine(const double* s, const double* e, bool /*hidden*/)
|
|
{
|
|
gp_Pnt p0 = makePoint(s);
|
|
gp_Pnt p1 = makePoint(e);
|
|
if (p0.IsEqual(p1,0.00000001))
|
|
return;
|
|
BRepBuilderAPI_MakeEdge makeEdge(p0, p1);
|
|
TopoDS_Edge edge = makeEdge.Edge();
|
|
AddObject(new Part::TopoShape(edge));
|
|
}
|
|
|
|
|
|
void ImpExpDxfRead::OnReadPoint(const double* s)
|
|
{
|
|
BRepBuilderAPI_MakeVertex makeVertex(makePoint(s));
|
|
TopoDS_Vertex vertex = makeVertex.Vertex();
|
|
AddObject(new Part::TopoShape(vertex));
|
|
}
|
|
|
|
|
|
void ImpExpDxfRead::OnReadArc(const double* s, const double* e, const double* c, bool dir, bool /*hidden*/)
|
|
{
|
|
gp_Pnt p0 = makePoint(s);
|
|
gp_Pnt p1 = makePoint(e);
|
|
gp_Dir up(0, 0, 1);
|
|
if (!dir)
|
|
up = -up;
|
|
gp_Pnt pc = makePoint(c);
|
|
gp_Circ circle(gp_Ax2(pc, up), p0.Distance(pc));
|
|
if (circle.Radius() > 0) {
|
|
BRepBuilderAPI_MakeEdge makeEdge(circle, p0, p1);
|
|
TopoDS_Edge edge = makeEdge.Edge();
|
|
AddObject(new Part::TopoShape(edge));
|
|
}
|
|
else {
|
|
Base::Console().Warning("ImpExpDxf - ignore degenerate arc of circle\n");
|
|
}
|
|
}
|
|
|
|
|
|
void ImpExpDxfRead::OnReadCircle(const double* s, const double* c, bool dir, bool /*hidden*/)
|
|
{
|
|
gp_Pnt p0 = makePoint(s);
|
|
gp_Dir up(0, 0, 1);
|
|
if (!dir)
|
|
up = -up;
|
|
gp_Pnt pc = makePoint(c);
|
|
gp_Circ circle(gp_Ax2(pc, up), p0.Distance(pc));
|
|
if (circle.Radius() > 0) {
|
|
BRepBuilderAPI_MakeEdge makeEdge(circle);
|
|
TopoDS_Edge edge = makeEdge.Edge();
|
|
AddObject(new Part::TopoShape(edge));
|
|
}
|
|
else {
|
|
Base::Console().Warning("ImpExpDxf - ignore degenerate circle\n");
|
|
}
|
|
}
|
|
|
|
|
|
Handle(Geom_BSplineCurve) getSplineFromPolesAndKnots(struct SplineData& sd)
|
|
{
|
|
std::size_t numPoles = sd.control_points;
|
|
if (sd.controlx.size() > numPoles || sd.controly.size() > numPoles ||
|
|
sd.controlz.size() > numPoles || sd.weight.size() > numPoles) {
|
|
return nullptr;
|
|
}
|
|
|
|
// handle the poles
|
|
TColgp_Array1OfPnt occpoles(1, sd.control_points);
|
|
int index = 1;
|
|
for (auto x : sd.controlx) {
|
|
occpoles(index++).SetX(x);
|
|
}
|
|
|
|
index = 1;
|
|
for (auto y : sd.controly) {
|
|
occpoles(index++).SetY(y);
|
|
}
|
|
|
|
index = 1;
|
|
for (auto z : sd.controlz) {
|
|
occpoles(index++).SetZ(z);
|
|
}
|
|
|
|
// handle knots and mults
|
|
std::set<double> unique;
|
|
unique.insert(sd.knot.begin(), sd.knot.end());
|
|
|
|
int numKnots = int(unique.size());
|
|
TColStd_Array1OfInteger occmults(1, numKnots);
|
|
TColStd_Array1OfReal occknots(1, numKnots);
|
|
index = 1;
|
|
for (auto k : unique) {
|
|
size_t m = std::count(sd.knot.begin(), sd.knot.end(), k);
|
|
occknots(index) = k;
|
|
occmults(index) = m;
|
|
index++;
|
|
}
|
|
|
|
// handle weights
|
|
TColStd_Array1OfReal occweights(1, sd.control_points);
|
|
if (sd.weight.size() == std::size_t(sd.control_points)) {
|
|
index = 1;
|
|
for (auto w : sd.weight) {
|
|
occweights(index++) = w;
|
|
}
|
|
}
|
|
else {
|
|
// non-rational
|
|
for (int i=occweights.Lower(); i<=occweights.Upper(); i++) {
|
|
occweights(i) = 1.0;
|
|
}
|
|
}
|
|
|
|
Standard_Boolean periodic = sd.flag == 2;
|
|
Handle(Geom_BSplineCurve) geom = new Geom_BSplineCurve(occpoles, occweights, occknots, occmults, sd.degree, periodic);
|
|
return geom;
|
|
}
|
|
|
|
Handle(Geom_BSplineCurve) getInterpolationSpline(struct SplineData& sd)
|
|
{
|
|
std::size_t numPoints = sd.fit_points;
|
|
if (sd.fitx.size() > numPoints || sd.fity.size() > numPoints || sd.fitz.size() > numPoints) {
|
|
return nullptr;
|
|
}
|
|
|
|
// handle the poles
|
|
Handle(TColgp_HArray1OfPnt) fitpoints = new TColgp_HArray1OfPnt(1, sd.fit_points);
|
|
int index = 1;
|
|
for (auto x : sd.fitx) {
|
|
fitpoints->ChangeValue(index++).SetX(x);
|
|
}
|
|
|
|
index = 1;
|
|
for (auto y : sd.fity) {
|
|
fitpoints->ChangeValue(index++).SetY(y);
|
|
}
|
|
|
|
index = 1;
|
|
for (auto z : sd.fitz) {
|
|
fitpoints->ChangeValue(index++).SetZ(z);
|
|
}
|
|
|
|
Standard_Boolean periodic = sd.flag == 2;
|
|
GeomAPI_Interpolate interp(fitpoints, periodic, Precision::Confusion());
|
|
interp.Perform();
|
|
return interp.Curve();
|
|
}
|
|
|
|
void ImpExpDxfRead::OnReadSpline(struct SplineData& sd)
|
|
{
|
|
// https://documentation.help/AutoCAD-DXF/WS1a9193826455f5ff18cb41610ec0a2e719-79e1.htm
|
|
// Flags:
|
|
// 1: Closed, 2: Periodic, 4: Rational, 8: Planar, 16: Linear
|
|
|
|
try {
|
|
Handle(Geom_BSplineCurve) geom;
|
|
if (sd.control_points > 0)
|
|
geom = getSplineFromPolesAndKnots(sd);
|
|
else if (sd.fit_points > 0)
|
|
geom = getInterpolationSpline(sd);
|
|
|
|
if (geom.IsNull())
|
|
throw Standard_Failure();
|
|
|
|
BRepBuilderAPI_MakeEdge makeEdge(geom);
|
|
TopoDS_Edge edge = makeEdge.Edge();
|
|
AddObject(new Part::TopoShape(edge));
|
|
}
|
|
catch (const Standard_Failure&) {
|
|
Base::Console().Warning("ImpExpDxf - failed to create bspline\n");
|
|
}
|
|
}
|
|
|
|
|
|
void ImpExpDxfRead::OnReadEllipse(const double* c, double major_radius, double minor_radius, double rotation, double /*start_angle*/, double /*end_angle*/, bool dir)
|
|
{
|
|
gp_Dir up(0, 0, 1);
|
|
if(!dir)
|
|
up = -up;
|
|
gp_Pnt pc = makePoint(c);
|
|
gp_Elips ellipse(gp_Ax2(pc, up), major_radius * optionScaling, minor_radius * optionScaling);
|
|
ellipse.Rotate(gp_Ax1(pc,up),rotation);
|
|
if (ellipse.MinorRadius() > 0) {
|
|
BRepBuilderAPI_MakeEdge makeEdge(ellipse);
|
|
TopoDS_Edge edge = makeEdge.Edge();
|
|
AddObject(new Part::TopoShape(edge));
|
|
}
|
|
else {
|
|
Base::Console().Warning("ImpExpDxf - ignore degenerate ellipse\n");
|
|
}
|
|
}
|
|
|
|
|
|
void ImpExpDxfRead::OnReadText(const double *point, const double /*height*/, const char* text)
|
|
{
|
|
if (optionImportAnnotations) {
|
|
Base::Vector3d pt(point[0] * optionScaling, point[1] * optionScaling, point[2] * optionScaling);
|
|
if(LayerName().substr(0, 6) != "BLOCKS") {
|
|
App::Annotation *pcFeature = (App::Annotation *)document->addObject("App::Annotation", "Text");
|
|
pcFeature->LabelText.setValue(Deformat(text));
|
|
pcFeature->Position.setValue(pt);
|
|
}
|
|
//else std::cout << "skipped text in block: " << LayerName() << std::endl;
|
|
}
|
|
}
|
|
|
|
|
|
void ImpExpDxfRead::OnReadInsert(const double* point, const double* scale, const char* name, double rotation)
|
|
{
|
|
//std::cout << "Inserting block " << name << " rotation " << rotation << " pos " << point[0] << "," << point[1] << "," << point[2] << " scale " << scale[0] << "," << scale[1] << "," << scale[2] << std::endl;
|
|
std::string prefix = "BLOCKS ";
|
|
prefix += name;
|
|
prefix += " ";
|
|
auto checkScale = [=](double v) {
|
|
return v != 0.0 ? v : 1.0;
|
|
};
|
|
for(std::map<std::string,std::vector<Part::TopoShape*> > ::const_iterator i = layers.begin(); i != layers.end(); ++i) {
|
|
std::string k = i->first;
|
|
if(k.substr(0, prefix.size()) == prefix) {
|
|
BRep_Builder builder;
|
|
TopoDS_Compound comp;
|
|
builder.MakeCompound(comp);
|
|
std::vector<Part::TopoShape*> v = i->second;
|
|
for(std::vector<Part::TopoShape*>::const_iterator j = v.begin(); j != v.end(); ++j) {
|
|
const TopoDS_Shape& sh = (*j)->getShape();
|
|
if (!sh.IsNull())
|
|
builder.Add(comp, sh);
|
|
}
|
|
if (!comp.IsNull()) {
|
|
Part::TopoShape* pcomp = new Part::TopoShape(comp);
|
|
Base::Matrix4D mat;
|
|
mat.scale(checkScale(scale[0]),checkScale(scale[1]),checkScale(scale[2]));
|
|
mat.rotZ(rotation);
|
|
mat.move(point[0]*optionScaling,point[1]*optionScaling,point[2]*optionScaling);
|
|
pcomp->transformShape(mat,true);
|
|
AddObject(pcomp);
|
|
}
|
|
}
|
|
}
|
|
}
|
|
|
|
|
|
void ImpExpDxfRead::OnReadDimension(const double* s, const double* e, const double* point, double /*rotation*/)
|
|
{
|
|
if (optionImportAnnotations) {
|
|
Base::Interpreter().runString("import Draft");
|
|
Base::Interpreter().runStringArg("p1=FreeCAD.Vector(%f,%f,%f)",s[0]*optionScaling,s[1]*optionScaling,s[2]*optionScaling);
|
|
Base::Interpreter().runStringArg("p2=FreeCAD.Vector(%f,%f,%f)",e[0]*optionScaling,e[1]*optionScaling,e[2]*optionScaling);
|
|
Base::Interpreter().runStringArg("p3=FreeCAD.Vector(%f,%f,%f)",point[0]*optionScaling,point[1]*optionScaling,point[2]*optionScaling);
|
|
Base::Interpreter().runString("Draft.makeDimension(p1,p2,p3)");
|
|
}
|
|
}
|
|
|
|
|
|
void ImpExpDxfRead::AddObject(Part::TopoShape *shape)
|
|
{
|
|
//std::cout << "layer:" << LayerName() << std::endl;
|
|
std::vector <Part::TopoShape*> vec;
|
|
if (layers.count(LayerName()))
|
|
vec = layers[LayerName()];
|
|
vec.push_back(shape);
|
|
layers[LayerName()] = vec;
|
|
if (!optionGroupLayers) {
|
|
if(LayerName().substr(0, 6) != "BLOCKS") {
|
|
Part::Feature *pcFeature = (Part::Feature *)document->addObject("Part::Feature", "Shape");
|
|
pcFeature->Shape.setValue(shape->getShape());
|
|
}
|
|
}
|
|
}
|
|
|
|
|
|
std::string ImpExpDxfRead::Deformat(const char* text)
|
|
{
|
|
// this function removes DXF formatting from texts
|
|
std::stringstream ss;
|
|
bool escape = false; // turned on when finding an escape character
|
|
bool longescape = false; // turned on for certain escape codes that expect additional chars
|
|
for(unsigned int i = 0; i<strlen(text); i++) {
|
|
if (text[i] == '\\')
|
|
escape = true;
|
|
else if (escape) {
|
|
if (longescape) {
|
|
if (text[i] == ';') {
|
|
escape = false;
|
|
longescape = false;
|
|
}
|
|
} else {
|
|
if ( (text[i] == 'H') || (text[i] == 'h') ||
|
|
(text[i] == 'Q') || (text[i] == 'q') ||
|
|
(text[i] == 'W') || (text[i] == 'w') ||
|
|
(text[i] == 'F') || (text[i] == 'f') ||
|
|
(text[i] == 'A') || (text[i] == 'a') ||
|
|
(text[i] == 'C') || (text[i] == 'c') ||
|
|
(text[i] == 'T') || (text[i] == 't') )
|
|
longescape = true;
|
|
else {
|
|
if ( (text[i] == 'P') || (text[i] == 'p') )
|
|
ss << "\n";
|
|
escape = false;
|
|
}
|
|
}
|
|
}
|
|
else if ( (text[i] != '{') && (text[i] != '}') ) {
|
|
ss << text[i];
|
|
}
|
|
}
|
|
return ss.str();
|
|
}
|
|
|
|
|
|
void ImpExpDxfRead::AddGraphics() const
|
|
{
|
|
if (optionGroupLayers) {
|
|
for(std::map<std::string,std::vector<Part::TopoShape*> > ::const_iterator i = layers.begin(); i != layers.end(); ++i) {
|
|
BRep_Builder builder;
|
|
TopoDS_Compound comp;
|
|
builder.MakeCompound(comp);
|
|
std::string k = i->first;
|
|
if (k == "0") // FreeCAD doesn't like an object name being '0'...
|
|
k = "LAYER_0";
|
|
std::vector<Part::TopoShape*> v = i->second;
|
|
if(k.substr(0, 6) != "BLOCKS") {
|
|
for(std::vector<Part::TopoShape*>::const_iterator j = v.begin(); j != v.end(); ++j) {
|
|
const TopoDS_Shape& sh = (*j)->getShape();
|
|
if (!sh.IsNull())
|
|
builder.Add(comp, sh);
|
|
}
|
|
if (!comp.IsNull()) {
|
|
Part::Feature *pcFeature = (Part::Feature *)document->addObject("Part::Feature", k.c_str());
|
|
pcFeature->Shape.setValue(comp);
|
|
}
|
|
}
|
|
}
|
|
}
|
|
}
|
|
|
|
//******************************************************************************
|
|
// writing
|
|
|
|
void gPntToTuple(double* result, gp_Pnt& p)
|
|
{
|
|
result[0] = p.X();
|
|
result[1] = p.Y();
|
|
result[2] = p.Z();
|
|
}
|
|
|
|
point3D gPntTopoint3D(gp_Pnt& p)
|
|
{
|
|
point3D result;
|
|
result.x = p.X();
|
|
result.y = p.Y();
|
|
result.z = p.Z();
|
|
return result;
|
|
}
|
|
|
|
ImpExpDxfWrite::ImpExpDxfWrite(std::string filepath) :
|
|
CDxfWrite(filepath.c_str())
|
|
{
|
|
setOptionSource("User parameter:BaseApp/Preferences/Mod/Import");
|
|
setOptions();
|
|
}
|
|
|
|
ImpExpDxfWrite::~ImpExpDxfWrite()
|
|
{
|
|
}
|
|
|
|
void ImpExpDxfWrite::setOptions()
|
|
{
|
|
ParameterGrp::handle hGrp = App::GetApplication().GetParameterGroupByPath(getOptionSource().c_str());
|
|
optionMaxLength = hGrp->GetFloat("maxsegmentlength",5.0);
|
|
optionExpPoints = hGrp->GetBool("ExportPoints",false);
|
|
m_version = hGrp->GetInt("DxfVersionOut",14);
|
|
optionPolyLine = hGrp->GetBool("DiscretizeEllipses",false);
|
|
m_polyOverride = hGrp->GetBool("DiscretizeEllipses",false);
|
|
setDataDir(App::Application::getResourceDir() + "Mod/Import/DxfPlate/");
|
|
}
|
|
|
|
void ImpExpDxfWrite::exportShape(const TopoDS_Shape input)
|
|
{
|
|
//export Edges
|
|
TopExp_Explorer edges(input, TopAbs_EDGE);
|
|
for (int i = 1 ; edges.More(); edges.Next(),i++) {
|
|
const TopoDS_Edge& edge = TopoDS::Edge(edges.Current());
|
|
BRepAdaptor_Curve adapt(edge);
|
|
if (adapt.GetType() == GeomAbs_Circle) {
|
|
double f = adapt.FirstParameter();
|
|
double l = adapt.LastParameter();
|
|
gp_Pnt s = adapt.Value(f);
|
|
gp_Pnt e = adapt.Value(l);
|
|
if (fabs(l-f) > 1.0 && s.SquareDistance(e) < 0.001) {
|
|
exportCircle(adapt);
|
|
} else {
|
|
exportArc(adapt);
|
|
}
|
|
} else if (adapt.GetType() == GeomAbs_Ellipse) {
|
|
double f = adapt.FirstParameter();
|
|
double l = adapt.LastParameter();
|
|
gp_Pnt s = adapt.Value(f);
|
|
gp_Pnt e = adapt.Value(l);
|
|
if (fabs(l-f) > 1.0 && s.SquareDistance(e) < 0.001) {
|
|
if (m_polyOverride) {
|
|
if (m_version >= 14) {
|
|
exportLWPoly(adapt);
|
|
} else { //m_version < 14
|
|
exportPolyline(adapt);
|
|
}
|
|
} else if (optionPolyLine) {
|
|
if (m_version >= 14) {
|
|
exportLWPoly(adapt);
|
|
} else { //m_version < 14
|
|
exportPolyline(adapt);
|
|
}
|
|
} else { //no overrides, do what's right!
|
|
if (m_version < 14) {
|
|
exportPolyline(adapt);
|
|
} else {
|
|
exportEllipse(adapt);
|
|
}
|
|
}
|
|
} else { // it's an arc
|
|
if (m_polyOverride) {
|
|
if (m_version >= 14) {
|
|
exportLWPoly(adapt);
|
|
} else { //m_version < 14
|
|
exportPolyline(adapt);
|
|
}
|
|
} else if (optionPolyLine) {
|
|
if (m_version >= 14) {
|
|
exportLWPoly(adapt);
|
|
} else { //m_version < 14
|
|
exportPolyline(adapt);
|
|
}
|
|
} else { //no overrides, do what's right!
|
|
if (m_version < 14) {
|
|
exportPolyline(adapt);
|
|
} else {
|
|
exportEllipseArc(adapt);
|
|
}
|
|
}
|
|
}
|
|
} else if (adapt.GetType() == GeomAbs_BSplineCurve) {
|
|
if (m_polyOverride) {
|
|
if (m_version >= 14) {
|
|
exportLWPoly(adapt);
|
|
} else { //m_version < 14
|
|
exportPolyline(adapt);
|
|
}
|
|
} else if (optionPolyLine) {
|
|
if (m_version >= 14) {
|
|
exportLWPoly(adapt);
|
|
} else { //m_version < 14
|
|
exportPolyline(adapt);
|
|
}
|
|
} else { //no overrides, do what's right!
|
|
if (m_version < 14) {
|
|
exportPolyline(adapt);
|
|
} else {
|
|
exportBSpline(adapt);
|
|
}
|
|
}
|
|
} else if (adapt.GetType() == GeomAbs_BezierCurve) {
|
|
exportBCurve(adapt);
|
|
} else if (adapt.GetType() == GeomAbs_Line) {
|
|
exportLine(adapt);
|
|
} else {
|
|
Base::Console().Warning("ImpExpDxf - unknown curve type: %d\n",adapt.GetType());
|
|
}
|
|
}
|
|
|
|
if (optionExpPoints) {
|
|
TopExp_Explorer verts(input, TopAbs_VERTEX);
|
|
std::vector<gp_Pnt> duplicates;
|
|
for (int i = 1 ; verts.More(); verts.Next(),i++) {
|
|
const TopoDS_Vertex& v = TopoDS::Vertex(verts.Current());
|
|
gp_Pnt p = BRep_Tool::Pnt(v);
|
|
duplicates.push_back(p);
|
|
}
|
|
|
|
std::sort(duplicates.begin(),duplicates.end(),ImpExpDxfWrite::gp_PntCompare);
|
|
auto newEnd = std::unique(duplicates.begin(),duplicates.end(),ImpExpDxfWrite::gp_PntEqual);
|
|
std::vector<gp_Pnt> uniquePts(duplicates.begin(),newEnd);
|
|
for (auto& p: uniquePts) {
|
|
double point[3] = {0,0,0};
|
|
gPntToTuple(point, p);
|
|
writePoint(point);
|
|
}
|
|
}
|
|
}
|
|
|
|
bool ImpExpDxfWrite::gp_PntEqual(gp_Pnt p1, gp_Pnt p2)
|
|
{
|
|
bool result = false;
|
|
if (p1.IsEqual(p2,Precision::Confusion())) {
|
|
result = true;
|
|
}
|
|
return result;
|
|
}
|
|
|
|
//is p1 "less than" p2?
|
|
bool ImpExpDxfWrite::gp_PntCompare(gp_Pnt p1, gp_Pnt p2)
|
|
{
|
|
bool result = false;
|
|
if (!(p1.IsEqual(p2,Precision::Confusion()))) { //ie v1 != v2
|
|
if (!(fabs(p1.X() - p2.X()) < Precision::Confusion())) { // x1 != x2
|
|
result = p1.X() < p2.X();
|
|
} else if (!(fabs(p1.Y() - p2.Y()) < Precision::Confusion())) { // y1 != y2
|
|
result = p1.Y() < p2.Y();
|
|
} else {
|
|
result = p1.Z() < p2.Z();
|
|
}
|
|
}
|
|
return result;
|
|
}
|
|
|
|
|
|
void ImpExpDxfWrite::exportCircle(BRepAdaptor_Curve& c)
|
|
{
|
|
gp_Circ circ = c.Circle();
|
|
gp_Pnt p = circ.Location();
|
|
double center[3] = {0,0,0};
|
|
gPntToTuple(center, p);
|
|
|
|
double radius = circ.Radius();
|
|
|
|
writeCircle(center, radius);
|
|
}
|
|
|
|
void ImpExpDxfWrite::exportEllipse(BRepAdaptor_Curve& c)
|
|
{
|
|
gp_Elips ellp = c.Ellipse();
|
|
gp_Pnt p = ellp.Location();
|
|
double center[3] = {0,0,0};
|
|
gPntToTuple(center, p);
|
|
|
|
double major = ellp.MajorRadius();
|
|
double minor = ellp.MinorRadius();
|
|
|
|
gp_Dir xaxis = ellp.XAxis().Direction(); //direction of major axis
|
|
//rotation appears to be the clockwise(?) angle between major & +Y??
|
|
double rotation = xaxis.AngleWithRef(gp_Dir(0, 1, 0), gp_Dir(0, 0, 1));
|
|
|
|
//2*M_PI = 6.28319 is invalid(doesn't display in LibreCAD), but 2PI = 6.28318 is valid!
|
|
//writeEllipse(center, major, minor, rotation, 0.0, 2 * M_PI, true );
|
|
writeEllipse(center, major, minor, rotation, 0.0, 6.28318, true );
|
|
}
|
|
|
|
void ImpExpDxfWrite::exportArc(BRepAdaptor_Curve& c)
|
|
{
|
|
gp_Circ circ = c.Circle();
|
|
gp_Pnt p = circ.Location();
|
|
double center[3] = {0,0,0};
|
|
gPntToTuple(center, p);
|
|
|
|
double f = c.FirstParameter();
|
|
double l = c.LastParameter();
|
|
gp_Pnt s = c.Value(f);
|
|
double start[3];
|
|
gPntToTuple(start, s);
|
|
gp_Pnt m = c.Value((l+f)/2.0);
|
|
gp_Pnt e = c.Value(l);
|
|
double end[3] = {0,0,0};
|
|
gPntToTuple(end, e);
|
|
|
|
gp_Vec v1(m,s);
|
|
gp_Vec v2(m,e);
|
|
gp_Vec v3(0,0,1);
|
|
double a = v3.DotCross(v1,v2);
|
|
|
|
bool dir = (a < 0) ? true: false;
|
|
writeArc(start, end, center, dir );
|
|
}
|
|
|
|
void ImpExpDxfWrite::exportEllipseArc(BRepAdaptor_Curve& c)
|
|
{
|
|
gp_Elips ellp = c.Ellipse();
|
|
gp_Pnt p = ellp.Location();
|
|
double center[3] = {0,0,0};
|
|
gPntToTuple(center, p);
|
|
|
|
double major = ellp.MajorRadius();
|
|
double minor = ellp.MinorRadius();
|
|
|
|
gp_Dir xaxis = ellp.XAxis().Direction(); //direction of major axis
|
|
//rotation appears to be the clockwise angle between major & +Y??
|
|
double rotation = xaxis.AngleWithRef(gp_Dir(0, 1, 0), gp_Dir(0, 0, 1));
|
|
|
|
double f = c.FirstParameter();
|
|
double l = c.LastParameter();
|
|
gp_Pnt s = c.Value(f);
|
|
gp_Pnt m = c.Value((l+f)/2.0);
|
|
gp_Pnt e = c.Value(l);
|
|
|
|
gp_Vec v1(m,s);
|
|
gp_Vec v2(m,e);
|
|
gp_Vec v3(0,0,1);
|
|
double a = v3.DotCross(v1,v2); // a = v3 dot (v1 cross v2)
|
|
// relates to "handedness" of 3 vectors
|
|
// a > 0 ==> v2 is CCW from v1 (righthanded)?
|
|
// a < 0 ==> v2 is CW from v1 (lefthanded)?
|
|
|
|
double startAngle = fmod(f,2.0*M_PI); //revolutions
|
|
double endAngle = fmod(l,2.0*M_PI);
|
|
bool endIsCW = (a < 0) ? true: false; //if !endIsCW swap(start,end)
|
|
//not sure if this is a hack or not. seems to make valid arcs.
|
|
if (!endIsCW) {
|
|
startAngle = -startAngle;
|
|
endAngle = -endAngle;
|
|
}
|
|
|
|
writeEllipse(center, major, minor, rotation, startAngle, endAngle, endIsCW);
|
|
}
|
|
|
|
void ImpExpDxfWrite::exportBSpline(BRepAdaptor_Curve& c)
|
|
{
|
|
SplineDataOut sd;
|
|
Handle(Geom_BSplineCurve) spline;
|
|
double f,l;
|
|
gp_Pnt s,ePt;
|
|
|
|
Standard_Real tol3D = 0.001;
|
|
Standard_Integer maxDegree = 3, maxSegment = 200;
|
|
Handle(BRepAdaptor_HCurve) hCurve = new BRepAdaptor_HCurve(c);
|
|
Approx_Curve3d approx(hCurve, tol3D, GeomAbs_C0, maxSegment, maxDegree);
|
|
if (approx.IsDone() && approx.HasResult()) {
|
|
spline = approx.Curve();
|
|
} else {
|
|
if (approx.HasResult()) { //result, but not within tolerance
|
|
spline = approx.Curve();
|
|
Base::Console().Message("DxfWrite::exportBSpline - result not within tolerance\n");
|
|
} else {
|
|
f = c.FirstParameter();
|
|
l = c.LastParameter();
|
|
s = c.Value(f);
|
|
ePt = c.Value(l);
|
|
Base::Console().Message("DxfWrite::exportBSpline - no result- from:(%.3f,%.3f) to:(%.3f,%.3f) poles: %d\n",
|
|
s.X(),s.Y(),ePt.X(),ePt.Y(),spline->NbPoles());
|
|
TColgp_Array1OfPnt controlPoints(0,1);
|
|
controlPoints.SetValue(0,s);
|
|
controlPoints.SetValue(1,ePt);
|
|
spline = GeomAPI_PointsToBSpline(controlPoints,1).Curve();
|
|
}
|
|
}
|
|
//WF? norm of surface containing curve??
|
|
sd.norm.x = 0.0;
|
|
sd.norm.y = 0.0;
|
|
sd.norm.z = 1.0;
|
|
|
|
sd.flag = spline->IsClosed();
|
|
sd.flag += spline->IsPeriodic()*2;
|
|
sd.flag += spline->IsRational()*4;
|
|
sd.flag += 8; //planar spline
|
|
|
|
sd.degree = spline->Degree();
|
|
sd.control_points = spline->NbPoles();
|
|
sd.knots = spline->NbKnots();
|
|
gp_Pnt p;
|
|
spline->D0(spline->FirstParameter(),p);
|
|
sd.starttan = gPntTopoint3D(p);
|
|
spline->D0(spline->LastParameter(),p);
|
|
sd.endtan = gPntTopoint3D(p);
|
|
|
|
//next bit is from DrawingExport.cpp (Dan Falk?).
|
|
Standard_Integer m = 0;
|
|
if (spline->IsPeriodic()) {
|
|
m = spline->NbPoles() + 2*spline->Degree() - spline->Multiplicity(1) + 2;
|
|
}
|
|
else {
|
|
for (int i=1; i<= spline->NbKnots(); i++)
|
|
m += spline->Multiplicity(i);
|
|
}
|
|
TColStd_Array1OfReal knotsequence(1,m);
|
|
spline->KnotSequence(knotsequence);
|
|
for (int i = knotsequence.Lower() ; i <= knotsequence.Upper(); i++) {
|
|
sd.knot.push_back(knotsequence(i));
|
|
}
|
|
sd.knots = knotsequence.Length();
|
|
|
|
TColgp_Array1OfPnt poles(1,spline->NbPoles());
|
|
spline->Poles(poles);
|
|
for (int i = poles.Lower(); i <= poles.Upper(); i++) {
|
|
sd.control.push_back(gPntTopoint3D(poles(i)));
|
|
}
|
|
//OCC doesn't have separate lists for control points and fit points.
|
|
|
|
writeSpline(sd);
|
|
}
|
|
|
|
void ImpExpDxfWrite::exportBCurve(BRepAdaptor_Curve& c)
|
|
{
|
|
(void) c;
|
|
Base::Console().Message("BCurve dxf export not yet supported\n");
|
|
}
|
|
|
|
void ImpExpDxfWrite::exportLine(BRepAdaptor_Curve& c)
|
|
{
|
|
double f = c.FirstParameter();
|
|
double l = c.LastParameter();
|
|
gp_Pnt s = c.Value(f);
|
|
double start[3] = {0,0,0};
|
|
gPntToTuple(start, s);
|
|
gp_Pnt e = c.Value(l);
|
|
double end[3] = {0,0,0};
|
|
gPntToTuple(end, e);
|
|
writeLine(start, end);
|
|
}
|
|
|
|
void ImpExpDxfWrite::exportLWPoly(BRepAdaptor_Curve& c)
|
|
{
|
|
LWPolyDataOut pd;
|
|
pd.Flag = c.IsClosed();
|
|
pd.Elev = 0.0;
|
|
pd.Thick = 0.0;
|
|
pd.Extr.x = 0.0;
|
|
pd.Extr.y = 0.0;
|
|
pd.Extr.z = 1.0;
|
|
pd.nVert = 0;
|
|
|
|
GCPnts_UniformAbscissa discretizer;
|
|
discretizer.Initialize (c, optionMaxLength);
|
|
std::vector<point3D> points;
|
|
if (discretizer.IsDone () && discretizer.NbPoints () > 0) {
|
|
int nbPoints = discretizer.NbPoints ();
|
|
for (int i=1; i<=nbPoints; i++) {
|
|
gp_Pnt p = c.Value (discretizer.Parameter (i));
|
|
pd.Verts.push_back(gPntTopoint3D(p));
|
|
}
|
|
pd.nVert = discretizer.NbPoints ();
|
|
writeLWPolyLine(pd);
|
|
}
|
|
}
|
|
|
|
void ImpExpDxfWrite::exportPolyline(BRepAdaptor_Curve& c)
|
|
{
|
|
LWPolyDataOut pd;
|
|
pd.Flag = c.IsClosed();
|
|
pd.Elev = 0.0;
|
|
pd.Thick = 0.0;
|
|
pd.Extr.x = 0.0;
|
|
pd.Extr.y = 0.0;
|
|
pd.Extr.z = 1.0;
|
|
pd.nVert = 0;
|
|
|
|
GCPnts_UniformAbscissa discretizer;
|
|
discretizer.Initialize (c, optionMaxLength);
|
|
std::vector<point3D> points;
|
|
if (discretizer.IsDone () && discretizer.NbPoints () > 0) {
|
|
int nbPoints = discretizer.NbPoints ();
|
|
for (int i=1; i<=nbPoints; i++) {
|
|
gp_Pnt p = c.Value (discretizer.Parameter (i));
|
|
pd.Verts.push_back(gPntTopoint3D(p));
|
|
}
|
|
pd.nVert = discretizer.NbPoints ();
|
|
writePolyline(pd);
|
|
}
|
|
}
|
|
|
|
void ImpExpDxfWrite::exportText(const char* text, Base::Vector3d position1, Base::Vector3d position2, double size, int just)
|
|
{
|
|
double location1[3] = {0,0,0};
|
|
location1[0] = position1.x;
|
|
location1[1] = position1.y;
|
|
location1[2] = position1.z;
|
|
double location2[3] = {0,0,0};
|
|
location2[0] = position2.x;
|
|
location2[1] = position2.y;
|
|
location2[2] = position2.z;
|
|
|
|
writeText(text, location1, location2, size, just);
|
|
}
|
|
|
|
void ImpExpDxfWrite::exportLinearDim(Base::Vector3d textLocn, Base::Vector3d lineLocn,
|
|
Base::Vector3d extLine1Start, Base::Vector3d extLine2Start,
|
|
char* dimText, int type)
|
|
{
|
|
double text[3] = {0,0,0};
|
|
text[0] = textLocn.x;
|
|
text[1] = textLocn.y;
|
|
text[2] = textLocn.z;
|
|
double line[3] = {0,0,0};
|
|
line[0] = lineLocn.x;
|
|
line[1] = lineLocn.y;
|
|
line[2] = lineLocn.z;
|
|
double ext1[3] = {0,0,0};
|
|
ext1[0] = extLine1Start.x;
|
|
ext1[1] = extLine1Start.y;
|
|
ext1[2] = extLine1Start.z;
|
|
double ext2[3] = {0,0,0};
|
|
ext2[0] = extLine2Start.x;
|
|
ext2[1] = extLine2Start.y;
|
|
ext2[2] = extLine2Start.z;
|
|
writeLinearDim(text, line, ext1,ext2,dimText, type);
|
|
}
|
|
|
|
void ImpExpDxfWrite::exportAngularDim(Base::Vector3d textLocn, Base::Vector3d lineLocn,
|
|
Base::Vector3d extLine1End, Base::Vector3d extLine2End,
|
|
Base::Vector3d apexPoint,
|
|
char* dimText)
|
|
{
|
|
double text[3] = {0,0,0};
|
|
text[0] = textLocn.x;
|
|
text[1] = textLocn.y;
|
|
text[2] = textLocn.z;
|
|
double line[3] = {0,0,0};
|
|
line[0] = lineLocn.x;
|
|
line[1] = lineLocn.y;
|
|
line[2] = lineLocn.z;
|
|
double ext1[3] = {0,0,0};
|
|
ext1[0] = extLine1End.x;
|
|
ext1[1] = extLine1End.y;
|
|
ext1[2] = extLine1End.z;
|
|
double ext2[3] = {0,0,0};
|
|
ext2[0] = extLine2End.x;
|
|
ext2[1] = extLine2End.y;
|
|
ext2[2] = extLine2End.z;
|
|
double apex[3] = {0,0,0};
|
|
apex[0] = apexPoint.x;
|
|
apex[1] = apexPoint.y;
|
|
apex[2] = apexPoint.z;
|
|
writeAngularDim(text, line, apex, ext1, apex, ext2, dimText);
|
|
}
|
|
|
|
void ImpExpDxfWrite::exportRadialDim(Base::Vector3d centerPoint, Base::Vector3d textLocn,
|
|
Base::Vector3d arcPoint,
|
|
char* dimText)
|
|
{
|
|
double center[3] = {0,0,0};
|
|
center[0] = centerPoint.x;
|
|
center[1] = centerPoint.y;
|
|
center[2] = centerPoint.z;
|
|
double text[3] = {0,0,0};
|
|
text[0] = textLocn.x;
|
|
text[1] = textLocn.y;
|
|
text[2] = textLocn.z;
|
|
double arc[3] = {0,0,0};
|
|
arc[0] = arcPoint.x;
|
|
arc[1] = arcPoint.y;
|
|
arc[2] = arcPoint.z;
|
|
writeRadialDim(center, text, arc, dimText);
|
|
|
|
}
|
|
|
|
void ImpExpDxfWrite::exportDiametricDim(Base::Vector3d textLocn,
|
|
Base::Vector3d arcPoint1, Base::Vector3d arcPoint2,
|
|
char* dimText)
|
|
{
|
|
double text[3] = {0,0,0};
|
|
text[0] = textLocn.x;
|
|
text[1] = textLocn.y;
|
|
text[2] = textLocn.z;
|
|
double arc1[3] = {0,0,0};
|
|
arc1[0] = arcPoint1.x;
|
|
arc1[1] = arcPoint1.y;
|
|
arc1[2] = arcPoint1.z;
|
|
double arc2[3] = {0,0,0};
|
|
arc2[0] = arcPoint2.x;
|
|
arc2[1] = arcPoint2.y;
|
|
arc2[2] = arcPoint2.z;
|
|
writeDiametricDim(text, arc1, arc2, dimText);
|
|
}
|