X-Git-Url: http://git.salome-platform.org/gitweb/?a=blobdiff_plain;f=src%2FGeomAPI%2FGeomAPI_Edge.cpp;h=ede977f30df08670d940ebbf007986f89286d43f;hb=6e421e939851e0de46554ae45a3ca0e1f67cd91d;hp=f9572daa79949e4da668fc374466dc078f598de2;hpb=9aad36a6b453e73b05356af4662a66a32cb71f1e;p=modules%2Fshaper.git diff --git a/src/GeomAPI/GeomAPI_Edge.cpp b/src/GeomAPI/GeomAPI_Edge.cpp index f9572daa7..ede977f30 100644 --- a/src/GeomAPI/GeomAPI_Edge.cpp +++ b/src/GeomAPI/GeomAPI_Edge.cpp @@ -1,28 +1,63 @@ -// Copyright (C) 2014-20xx CEA/DEN, EDF R&D - -// File: GeomAPI_Edge.cpp -// Created: 24 Jul 2014 -// Author: Artem ZHIDKOV +// Copyright (C) 2014-2019 CEA/DEN, EDF R&D +// +// This library is free software; you can redistribute it and/or +// modify it under the terms of the GNU Lesser General Public +// License as published by the Free Software Foundation; either +// version 2.1 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 +// Lesser General Public License for more details. +// +// You should have received a copy of the GNU Lesser General Public +// License along with this library; if not, write to the Free Software +// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA +// +// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com +// #include +#include #include #include #include #include +#include +#include + +#include #include #include #include +#include #include +#include +#include #include #include #include +#include +#include +#include +#include #include #include +#include +#include +#include + +#include GeomAPI_Edge::GeomAPI_Edge() - : GeomAPI_Shape() { + TopoDS_Edge* anEdge = new TopoDS_Edge; + + BRep_Builder aBuilder; + aBuilder.MakeEdge(*anEdge); + + setImpl(anEdge); } GeomAPI_Edge::GeomAPI_Edge(const std::shared_ptr& theShape) @@ -37,18 +72,36 @@ bool GeomAPI_Edge::isLine() const const TopoDS_Shape& aShape = const_cast(this)->impl(); double aFirst, aLast; Handle(Geom_Curve) aCurve = BRep_Tool::Curve((const TopoDS_Edge&)aShape, aFirst, aLast); + if (aCurve.IsNull()) // degenerative edge + return false; if (aCurve->IsKind(STANDARD_TYPE(Geom_Line))) return true; return false; } +/// extracts a circle curve from the arbitrary curve, returns null is it is different type +static Handle(Geom_Circle) circ(const Handle(Geom_Curve) theCurve) +{ + Handle(Geom_Circle) aResult = Handle(Geom_Circle)::DownCast(theCurve); + if (!aResult.IsNull()) + return aResult; + // check this may be a trimmed curve that contains circle inside + Handle(Geom_TrimmedCurve) aTrimmed = Handle(Geom_TrimmedCurve)::DownCast(theCurve); + while(!aTrimmed.IsNull()) { + aResult = Handle(Geom_Circle)::DownCast(aTrimmed->BasisCurve()); + if (!aResult.IsNull()) + return aResult; + aTrimmed = Handle(Geom_TrimmedCurve)::DownCast(aTrimmed->BasisCurve()); + } + return aResult; // null, not circle +} + bool GeomAPI_Edge::isCircle() const { const TopoDS_Shape& aShape = const_cast(this)->impl(); double aFirst, aLast; Handle(Geom_Curve) aCurve = BRep_Tool::Curve((const TopoDS_Edge&)aShape, aFirst, aLast); - if (aCurve->IsKind(STANDARD_TYPE(Geom_Circle))) - { + if (!circ(aCurve).IsNull()) { // Check the difference of first and last parameters to be equal to the curve period if (Abs(aLast - aFirst - aCurve->Period()) < Precision::PConfusion()) return true; @@ -61,8 +114,7 @@ bool GeomAPI_Edge::isArc() const const TopoDS_Shape& aShape = const_cast(this)->impl(); double aFirst, aLast; Handle(Geom_Curve) aCurve = BRep_Tool::Curve((const TopoDS_Edge&)aShape, aFirst, aLast); - if (aCurve->IsKind(STANDARD_TYPE(Geom_Circle))) - { + if (!circ(aCurve).IsNull()) { // Check the difference of first and last parameters is not equal the curve period if (Abs(aLast - aFirst - aCurve->Period()) >= Precision::PConfusion()) return true; @@ -70,6 +122,18 @@ bool GeomAPI_Edge::isArc() const return false; } +bool GeomAPI_Edge::isEllipse() const +{ + const TopoDS_Shape& aShape = const_cast(this)->impl(); + double aFirst, aLast; + Handle(Geom_Curve) aCurve = BRep_Tool::Curve((const TopoDS_Edge&)aShape, aFirst, aLast); + if (aCurve.IsNull()) // degenerative edge + return false; + if (aCurve->IsKind(STANDARD_TYPE(Geom_Ellipse))) + return true; + return false; +} + std::shared_ptr GeomAPI_Edge::firstPoint() { const TopoDS_Shape& aShape = const_cast(this)->impl(); @@ -90,25 +154,40 @@ std::shared_ptr GeomAPI_Edge::lastPoint() return std::shared_ptr(new GeomAPI_Pnt(aPoint.X(), aPoint.Y(), aPoint.Z())); } -std::shared_ptr GeomAPI_Edge::circle() +std::shared_ptr GeomAPI_Edge::circle() const { const TopoDS_Shape& aShape = const_cast(this)->impl(); double aFirst, aLast; Handle(Geom_Curve) aCurve = BRep_Tool::Curve((const TopoDS_Edge&)aShape, aFirst, aLast); - if (aCurve) { - Handle(Geom_Circle) aCirc = Handle(Geom_Circle)::DownCast(aCurve); - if (aCirc) { - gp_Pnt aLoc = aCirc->Location(); - std::shared_ptr aCenter(new GeomAPI_Pnt(aLoc.X(), aLoc.Y(), aLoc.Z())); - gp_Dir anAxis = aCirc->Axis().Direction(); - std::shared_ptr aDir(new GeomAPI_Dir(anAxis.X(), anAxis.Y(), anAxis.Z())); - return std::shared_ptr(new GeomAPI_Circ(aCenter, aDir, aCirc->Radius())); - } + Handle(Geom_Circle) aCirc = circ(aCurve); + if (!aCirc.IsNull()) { + gp_Pnt aLoc = aCirc->Location(); + std::shared_ptr aCenter(new GeomAPI_Pnt(aLoc.X(), aLoc.Y(), aLoc.Z())); + gp_Dir anAxis = aCirc->Axis().Direction(); + std::shared_ptr aDir(new GeomAPI_Dir(anAxis.X(), anAxis.Y(), anAxis.Z())); + return std::shared_ptr(new GeomAPI_Circ(aCenter, aDir, aCirc->Radius())); } return std::shared_ptr(); // not circle } -std::shared_ptr GeomAPI_Edge::line() +std::shared_ptr GeomAPI_Edge::ellipse() const +{ + const TopoDS_Shape& aShape = const_cast(this)->impl(); + double aFirst, aLast; + Handle(Geom_Curve) aCurve = BRep_Tool::Curve((const TopoDS_Edge&)aShape, aFirst, aLast); + if (!aCurve.IsNull()) { + Handle(Geom_Ellipse) aElips = Handle(Geom_Ellipse)::DownCast(aCurve); + if (!aElips.IsNull()) { + gp_Elips aGpElips = aElips->Elips(); + std::shared_ptr aEllipse(new GeomAPI_Ellipse()); + aEllipse->setImpl(new gp_Elips(aGpElips)); + return aEllipse; + } + } + return std::shared_ptr(); // not ellipse +} + +std::shared_ptr GeomAPI_Edge::line() const { const TopoDS_Shape& aShape = const_cast(this)->impl(); double aFirst, aLast; @@ -131,9 +210,14 @@ std::shared_ptr GeomAPI_Edge::line() bool GeomAPI_Edge::isEqual(const std::shared_ptr theEdge) const { + if (!theEdge.get() || ! theEdge->isEdge()) + return false; const TopoDS_Shape& aMyShape = const_cast(this)->impl(); const TopoDS_Shape& aInShape = theEdge->impl(); + if (aMyShape.IsNull() || aInShape.IsNull()) + return false; + if (aMyShape.ShapeType() != aInShape.ShapeType()) return false; @@ -142,25 +226,173 @@ bool GeomAPI_Edge::isEqual(const std::shared_ptr theEdge) const double aInStart, aInEnd; Handle(Geom_Curve) aInCurve = BRep_Tool::Curve(TopoDS::Edge(aInShape), aInStart, aInEnd); + // Check that end point parameters are the same + if ((aMyStart != aInStart) || (aMyEnd != aInEnd)) + return false; + // Check that curves a the same type GeomAdaptor_Curve aMyAdaptor(aMyCurve); GeomAdaptor_Curve aInAdaptor(aInCurve); if (aMyAdaptor.GetType() != aInAdaptor.GetType()) return false; - // Check that end point parameters are the same - if ((aMyStart != aInStart) || (aMyEnd != aInEnd)) - return false; - // Check that end points are equal gp_Pnt aMyPnt1 = aMyAdaptor.Value(aMyStart); gp_Pnt aMyPnt2 = aMyAdaptor.Value(aMyEnd); gp_Pnt aInPnt1 = aInAdaptor.Value(aInStart); gp_Pnt aInPnt2 = aInAdaptor.Value(aInEnd); - if ((!aMyPnt1.IsEqual(aInPnt1, Precision::Confusion())) || + if ((!aMyPnt1.IsEqual(aInPnt1, Precision::Confusion())) || (!aMyPnt2.IsEqual(aInPnt2, Precision::Confusion()))) return false; return true; } + +// LCOV_EXCL_START +void GeomAPI_Edge::getRange(double& theFirst, double& theLast) const +{ + const TopoDS_Shape& aShape = const_cast(this)->impl(); + Handle(Geom_Curve) aCurve = BRep_Tool::Curve((const TopoDS_Edge&)aShape, theFirst, theLast); +} + +bool GeomAPI_Edge::isInPlane(std::shared_ptr thePlane) const +{ + double aFirst, aLast; + const TopoDS_Shape& aShape = const_cast(this)->impl(); + Handle(Geom_Curve) aCurve = BRep_Tool::Curve((const TopoDS_Edge&)aShape, aFirst, aLast); + if (aCurve.IsNull()) + return false; + + double A, B, C, D; + thePlane->coefficients(A, B, C, D); + gp_Pln aPlane(A, B, C, D); + + bool inPlane = false; + if (aCurve->IsKind(STANDARD_TYPE(Geom_Line))) { + // check start and end points on the plane + gp_Pnt aFirstPnt = aCurve->Value(aFirst); + gp_Pnt aLastPnt = aCurve->Value(aLast); + inPlane = aPlane.SquareDistance(aFirstPnt) < Precision::SquareConfusion() && + aPlane.SquareDistance(aLastPnt) < Precision::SquareConfusion(); + } else { + Handle(Geom_Circle) aCirc = circ(aCurve); + if (!aCirc.IsNull()) { + gp_Pnt aCenter = aCirc->Location(); + Standard_Real aDot = aPlane.Axis().Direction().Dot(aCirc->Axis().Direction()); + inPlane = aPlane.SquareDistance(aCenter) < Precision::SquareConfusion() && + Abs(Abs(aDot) - 1.0) < Precision::Confusion(); + } else { + // three points checking + gp_Pnt aFirstPnt = aCurve->Value(aFirst); + gp_Pnt aMidPnt = aCurve->Value((aFirst + aLast) / 2.); + gp_Pnt aLastPnt = aCurve->Value(aLast); + inPlane = aPlane.SquareDistance(aFirstPnt) < Precision::SquareConfusion() && + aPlane.SquareDistance(aMidPnt) < Precision::SquareConfusion() && + aPlane.SquareDistance(aLastPnt) < Precision::SquareConfusion(); + } + } + return inPlane; +} +// LCOV_EXCL_STOP + +void GeomAPI_Edge::intersectWithPlane(const std::shared_ptr thePlane, + std::list>& theResult) const +{ + double aFirst, aLast; + const TopoDS_Shape& aShape = const_cast(this)->impl(); + Handle(Geom_Curve) aCurve = BRep_Tool::Curve((const TopoDS_Edge&)aShape, aFirst, aLast); + if (!aCurve.IsNull()) { + double A, B, C, D; + thePlane->coefficients(A, B, C, D); + gp_Pln aPln(A, B, C, D); + + Handle(Geom_Plane) aPlane = new Geom_Plane(aPln); + GeomAPI_IntCS aIntersect; + aIntersect.Perform(aCurve, aPlane); + if (aIntersect.IsDone() && (aIntersect.NbPoints() > 0)) { + gp_Pnt aPnt; + for (int i = 1; i <= aIntersect.NbPoints(); i++) { + // check the parameter of intersection in the edge range + aIntersect.Parameters(i, A, B, C); + if (aCurve->IsPeriodic()) + C = ElCLib::InPeriod(C, aFirst, aFirst + aCurve->Period()); + if (C < aFirst - Precision::PConfusion() || C > aLast + Precision::PConfusion()) + continue; + + // obtain intersection point + aPnt = aIntersect.Point(i); + std::shared_ptr aPntPtr(new GeomAPI_Pnt(aPnt.X(), aPnt.Y(), aPnt.Z())); + theResult.push_back(aPntPtr); + } + } + } +} + +double GeomAPI_Edge::length() const +{ + const TopoDS_Edge& anEdge = TopoDS::Edge(impl()); + BRepAdaptor_Curve aBRepAdaptor = BRepAdaptor_Curve(anEdge); + Adaptor3d_Curve* anAdaptor3d = &aBRepAdaptor; + return GCPnts_AbscissaPoint::Length(*anAdaptor3d); +} + +bool GeomAPI_Edge::isClosed() const +{ + const TopoDS_Shape& aShape = const_cast(this)->impl(); + if (aShape.IsNull()) + return false; + double aFirst, aLast; + Handle(Geom_Curve) aCurve = BRep_Tool::Curve((const TopoDS_Edge&)aShape, aFirst, aLast); + if (aCurve.IsNull() || !aCurve->IsPeriodic()) + return false; + aLast += aLast > aFirst ? -aCurve->Period() : aCurve->Period();; + + return fabs(aFirst - aLast) < 1.e-9; +} + +bool GeomAPI_Edge::isDegenerated() const +{ + const TopoDS_Shape& aShape = const_cast(this)->impl(); + if (aShape.IsNull() || aShape.ShapeType() != TopAbs_EDGE) + return false; + return BRep_Tool::Degenerated(TopoDS::Edge(aShape)); +} + +void GeomAPI_Edge::setFirstPointTolerance(const double theTolerance) +{ + TopoDS_Edge anEdge = impl(); + TopoDS_Vertex aVFirst, aVLast; + TopExp::Vertices(anEdge, aVFirst, aVLast); + BRep_Builder().UpdateVertex(aVFirst, theTolerance); +} + +void GeomAPI_Edge::setLastPointTolerance(const double theTolerance) +{ + TopoDS_Edge anEdge = impl(); + TopoDS_Vertex aVFirst, aVLast; + TopExp::Vertices(anEdge, aVFirst, aVLast); + BRep_Builder().UpdateVertex(aVLast, theTolerance); +} + +GeomPointPtr GeomAPI_Edge::middlePoint() const +{ + GeomPointPtr aMiddlePoint; + + const TopoDS_Edge& anEdge = impl(); + if (anEdge.IsNull()) + return aMiddlePoint; + double aFirst, aLast; + Handle(Geom_Curve) aCurve = BRep_Tool::Curve(anEdge, aFirst, aLast); + if (aCurve.IsNull()) + return aMiddlePoint; + + static const int NB_POINTS = 3; + GeomAdaptor_Curve aCurveAdaptor(aCurve, aFirst, aLast); + GCPnts_UniformAbscissa anAlgo(aCurveAdaptor, NB_POINTS); + if (anAlgo.IsDone()) { + gp_Pnt aPnt = aCurveAdaptor.Value(anAlgo.Parameter(2)); + aMiddlePoint = GeomPointPtr(new GeomAPI_Pnt(aPnt.X(), aPnt.Y(), aPnt.Z())); + } + return aMiddlePoint; +}