+ PLate surface approximation
This commit is contained in:
parent
13e2b24ddf
commit
c8e241c9f3
|
@ -445,9 +445,6 @@ makeCompound(PyObject *self, PyObject *args)
|
|||
|
||||
static PyObject * makeFilledFace(PyObject *self, PyObject *args)
|
||||
{
|
||||
// http://opencascade.blogspot.com/2010/03/surface-modeling-part6.html
|
||||
// TODO: GeomPlate_BuildPlateSurface
|
||||
// TODO: GeomPlate_MakeApprox
|
||||
// TODO: BRepFeat_SplitShape
|
||||
PyObject *obj;
|
||||
PyObject *surf=0;
|
||||
|
|
|
@ -112,6 +112,7 @@
|
|||
#include "ConePy.h"
|
||||
#include "CylinderPy.h"
|
||||
#include "OffsetSurfacePy.h"
|
||||
#include "PlateSurfacePy.h"
|
||||
#include "PlanePy.h"
|
||||
#include "RectangularTrimmedSurfacePy.h"
|
||||
#include "SpherePy.h"
|
||||
|
@ -3328,8 +3329,7 @@ void GeomPlateSurface::Restore(Base::XMLReader &/*reader*/)
|
|||
|
||||
PyObject *GeomPlateSurface::getPyObject(void)
|
||||
{
|
||||
throw Base::NotImplementedError("GeomPlateSurface::getPyObject");
|
||||
// return new PlateSurfacePy(static_cast<GeomPlateSurface*>(this->clone()));
|
||||
return new PlateSurfacePy(static_cast<GeomPlateSurface*>(this->clone()));
|
||||
}
|
||||
|
||||
// -------------------------------------------------
|
||||
|
|
|
@ -22,12 +22,27 @@
|
|||
|
||||
#include "PreCompiled.h"
|
||||
#ifndef _PreComp_
|
||||
# include <cassert>
|
||||
# include <gp_Pln.hxx>
|
||||
# include <gp_Lin.hxx>
|
||||
# include <Adaptor3d_HCurveOnSurface.hxx>
|
||||
# include <Geom_Plane.hxx>
|
||||
# include <GeomAdaptor_HCurve.hxx>
|
||||
# include <GeomAPI_IntSS.hxx>
|
||||
# include <Geom_Line.hxx>
|
||||
# include <Geom_Point.hxx>
|
||||
# include <GeomPlate_BuildPlateSurface.hxx>
|
||||
# include <GeomPlate_CurveConstraint.hxx>
|
||||
# include <GeomPlate_MakeApprox.hxx>
|
||||
# include <GeomPlate_PlateG0Criterion.hxx>
|
||||
# include <GeomPlate_PointConstraint.hxx>
|
||||
# include <Precision.hxx>
|
||||
# include <Standard_Mutex.hxx>
|
||||
# include <Standard_TypeMismatch.hxx>
|
||||
# include <TColStd_ListOfTransient.hxx>
|
||||
# include <TColStd_ListIteratorOfListOfTransient.hxx>
|
||||
# include <TColgp_SequenceOfXY.hxx>
|
||||
# include <TColgp_SequenceOfXYZ.hxx>
|
||||
#endif
|
||||
|
||||
#include <Base/Vector3D.h>
|
||||
|
@ -81,3 +96,95 @@ bool Part::intersect(const gp_Pln& pln1, const gp_Pln& pln2, gp_Lin& lin)
|
|||
|
||||
return found;
|
||||
}
|
||||
|
||||
/*! The objects in \a theBoundaries must be of the type Adaptor3d_HCurveOnSurface or
|
||||
GeomAdaptor_HCurve or Geom_Point indicating type of a constraint. Otherwise an exception
|
||||
Standard_TypeMismatch is thrown.
|
||||
|
||||
If the \a theBoundaries list is empty then Standard_ConstructionError is thrown.
|
||||
|
||||
If the algorithm fails it returns a null surface.
|
||||
\see http://opencascade.blogspot.com/2010/03/surface-modeling-part6.html
|
||||
*/
|
||||
Handle_Geom_Surface
|
||||
Part::Tools::makeSurface(const TColStd_ListOfTransient &theBoundaries,
|
||||
const Standard_Real theTol,
|
||||
const Standard_Integer theNbPnts,
|
||||
const Standard_Integer theNbIter,
|
||||
const Standard_Integer theMaxDeg)
|
||||
{
|
||||
//constants for algorithm
|
||||
const Standard_Integer aNbIter = theNbIter; //number of algorithm iterations
|
||||
const Standard_Integer aNbPnts = theNbPnts; //sample points per each constraint
|
||||
const Standard_Integer aDeg = 3; //requested surface degree ?
|
||||
const Standard_Integer aMaxDeg = theMaxDeg;
|
||||
const Standard_Integer aMaxSeg = 10000;
|
||||
const Standard_Real aTol3d = 1.e-04;
|
||||
const Standard_Real aTol2d = 1.e-05;
|
||||
const Standard_Real anAngTol = 1.e-02; //angular
|
||||
const Standard_Real aCurvTol = 1.e-01; //curvature
|
||||
|
||||
Handle(Geom_Surface) aRes;
|
||||
GeomPlate_BuildPlateSurface aPlateBuilder (aDeg, aNbPnts, aNbIter, aTol2d, aTol3d, anAngTol, aCurvTol);
|
||||
|
||||
TColStd_ListIteratorOfListOfTransient anIt (theBoundaries);
|
||||
if (anIt.More()) {
|
||||
int i = 1;
|
||||
for (; anIt.More(); anIt.Next(), i++) {
|
||||
const Handle(Standard_Transient)& aCur = anIt.Value();
|
||||
if (aCur.IsNull()) {
|
||||
assert (0);
|
||||
Standard_ConstructionError::Raise ("Tools::makeSurface()");
|
||||
}
|
||||
else if (aCur->IsKind (STANDARD_TYPE (Adaptor3d_HCurveOnSurface))) {
|
||||
//G1 constraint
|
||||
const Handle(Adaptor3d_HCurveOnSurface)& aHCOS = Handle(Adaptor3d_HCurveOnSurface)::DownCast (aCur);
|
||||
Handle (GeomPlate_CurveConstraint) aConst = new GeomPlate_CurveConstraint (aHCOS, 1 /*GeomAbs_G1*/,aNbPnts, aTol3d, anAngTol, aCurvTol);
|
||||
aPlateBuilder.Add (aConst);
|
||||
}
|
||||
else if (aCur->IsKind (STANDARD_TYPE (GeomAdaptor_HCurve))) {
|
||||
//G0 constraint
|
||||
const Handle(GeomAdaptor_HCurve)& aHC = Handle(GeomAdaptor_HCurve)::DownCast (aCur);
|
||||
Handle (GeomPlate_CurveConstraint) aConst = new GeomPlate_CurveConstraint (aHC, 0 /*GeomAbs_G0*/, aNbPnts, aTol3d);
|
||||
aPlateBuilder.Add (aConst);
|
||||
}
|
||||
else if (aCur->IsKind (STANDARD_TYPE (Geom_Point))) {
|
||||
//Point constraint
|
||||
const Handle(Geom_Point)& aGP = Handle(Geom_Point)::DownCast (aCur);
|
||||
Handle(GeomPlate_PointConstraint) aConst = new GeomPlate_PointConstraint(aGP->Pnt(),0);
|
||||
aPlateBuilder.Add(aConst);
|
||||
}
|
||||
else {
|
||||
Standard_TypeMismatch::Raise ("Tools::makeSurface()");
|
||||
}
|
||||
}
|
||||
}
|
||||
else {
|
||||
Standard_ConstructionError::Raise ("Tools::makeSurface()");
|
||||
}
|
||||
|
||||
//construct
|
||||
aPlateBuilder.Perform();
|
||||
|
||||
if (!aPlateBuilder.IsDone()) {
|
||||
return aRes;
|
||||
}
|
||||
|
||||
const Handle(GeomPlate_Surface)& aPlate = aPlateBuilder.Surface();
|
||||
//approximation (see BRepFill_Filling - when no initial surface was given)
|
||||
Standard_Real aDMax = aPlateBuilder.G0Error();
|
||||
TColgp_SequenceOfXY aS2d;
|
||||
TColgp_SequenceOfXYZ aS3d;
|
||||
aPlateBuilder.Disc2dContour (4, aS2d);
|
||||
aPlateBuilder.Disc3dContour (4, 0, aS3d);
|
||||
Standard_Real aMax = Max (aTol3d, 10. * aDMax);
|
||||
GeomPlate_PlateG0Criterion aCriterion (aS2d, aS3d, aMax);
|
||||
{
|
||||
//data races in AdvApp2Var used by GeomApprox_Surface, use global mutex
|
||||
//Standard_Mutex::Sentry aSentry (theBSMutex);
|
||||
GeomPlate_MakeApprox aMakeApprox (aPlate, aCriterion, aTol3d, aMaxSeg, aMaxDeg);
|
||||
aRes = aMakeApprox.Surface();
|
||||
}
|
||||
|
||||
return aRes;
|
||||
}
|
||||
|
|
|
@ -28,9 +28,11 @@
|
|||
#include <gp_Vec.hxx>
|
||||
#include <gp_Dir.hxx>
|
||||
#include <gp_XYZ.hxx>
|
||||
#include <Handle_Geom_Surface.hxx>
|
||||
|
||||
class gp_Lin;
|
||||
class gp_Pln;
|
||||
class TColStd_ListOfTransient;
|
||||
|
||||
namespace Base {
|
||||
// Specialization for gp_Pnt
|
||||
|
@ -93,6 +95,17 @@ bool intersect(const gp_Pln& pln1, const gp_Pln& pln2, gp_Lin& lin);
|
|||
PartExport
|
||||
bool tangentialArc(const gp_Pnt& p0, const gp_Vec& v0, const gp_Pnt& p1, gp_Pnt& c, gp_Dir& a);
|
||||
|
||||
class PartExport Tools
|
||||
{
|
||||
public:
|
||||
Handle_Geom_Surface makeSurface (const TColStd_ListOfTransient& theBoundaries,
|
||||
const Standard_Real theTol,
|
||||
const Standard_Integer theNbPnts,
|
||||
const Standard_Integer theNbIter,
|
||||
const Standard_Integer theMaxDeg);
|
||||
|
||||
};
|
||||
|
||||
} //namespace Part
|
||||
|
||||
|
||||
|
|
Loading…
Reference in New Issue
Block a user