X-Git-Url: http://git.salome-platform.org/gitweb/?a=blobdiff_plain;f=src%2FGeomAPI%2FGeomAPI_Edge.cpp;h=03868875ba53de38688ae1abddfed3e399a942e7;hb=1f38083fe812d086f053a46910c2cec3c8167868;hp=b5d0de4b673c35ad38c0d532e3b8cbf8df7d87b3;hpb=1c3738ae81b02ba62136ac03a53a81a532b95141;p=modules%2Fshaper.git diff --git a/src/GeomAPI/GeomAPI_Edge.cpp b/src/GeomAPI/GeomAPI_Edge.cpp index b5d0de4b6..03868875b 100644 --- a/src/GeomAPI/GeomAPI_Edge.cpp +++ b/src/GeomAPI/GeomAPI_Edge.cpp @@ -1,4 +1,4 @@ -// Copyright (C) 2014-2017 CEA/DEN, EDF R&D +// 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 @@ -12,10 +12,9 @@ // // 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 +// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA // -// See http://www.salome-platform.org/ or -// email : webmaster.salome@opencascade.com +// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com // #include @@ -26,6 +25,7 @@ #include #include #include +#include #include @@ -35,11 +35,16 @@ #include #include #include +#include +#include #include #include #include +#include #include #include +#include +#include #include #include #include @@ -66,6 +71,48 @@ GeomAPI_Edge::GeomAPI_Edge(const std::shared_ptr& theShape) } } +void GeomAPI_Edge::vertices(std::shared_ptr& theStartVertex, + std::shared_ptr& theEndVertex) const +{ + const TopoDS_Edge& anEdge = impl(); + TopoDS_Vertex aStart, aEnd; + TopExp::Vertices(anEdge, aStart, aEnd); + theStartVertex.reset(new GeomAPI_Vertex); + theStartVertex->setImpl(new TopoDS_Vertex(aStart)); + theEndVertex.reset(new GeomAPI_Vertex); + theEndVertex->setImpl(new TopoDS_Vertex(aEnd)); +} + +static Handle(Geom_Curve) baseCurve(const TopoDS_Edge& theEdge) +{ + double aFirst, aLast; + Handle(Geom_Curve) aCurve = BRep_Tool::Curve(theEdge, aFirst, aLast); + while (aCurve->IsKind(STANDARD_TYPE(Geom_TrimmedCurve))) { + Handle(Geom_TrimmedCurve) tc = Handle(Geom_TrimmedCurve)::DownCast(aCurve); + aCurve = tc->BasisCurve(); + } + return aCurve; +} + +bool GeomAPI_Edge::isSameGeometry(const std::shared_ptr theShape) const +{ + if (!theShape->isEdge()) + return false; + if (isSame(theShape)) + return true; + + TopoDS_Edge anOwnEdge = TopoDS::Edge(impl()); + TopoDS_Edge anOtherEdge = TopoDS::Edge(theShape->impl()); + + Handle(Geom_Curve) anOwnCurve = baseCurve(anOwnEdge); + Handle(Geom_Curve) anOtherCurve = baseCurve(anOtherEdge); + GeomAPI_ExtremaCurveCurve anExtrema(anOwnCurve, anOtherCurve); + + bool isSame = anExtrema.Extrema().IsParallel() && + anExtrema.TotalLowerDistance() < Precision::Confusion(); + return isSame; +} + bool GeomAPI_Edge::isLine() const { const TopoDS_Shape& aShape = const_cast(this)->impl(); @@ -78,15 +125,29 @@ bool GeomAPI_Edge::isLine() const 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.IsNull()) // degenerative edge - return false; - 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; @@ -99,10 +160,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.IsNull()) // degenerative edge - return false; - 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; @@ -122,6 +180,18 @@ bool GeomAPI_Edge::isEllipse() const return false; } +bool GeomAPI_Edge::isBSpline() 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; + while (aCurve->IsKind(STANDARD_TYPE(Geom_TrimmedCurve))) + aCurve = Handle(Geom_TrimmedCurve)::DownCast(aCurve)->BasisCurve(); + return aCurve->IsKind(STANDARD_TYPE(Geom_BSplineCurve)); +} + std::shared_ptr GeomAPI_Edge::firstPoint() { const TopoDS_Shape& aShape = const_cast(this)->impl(); @@ -147,15 +217,13 @@ 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.IsNull()) { - Handle(Geom_Circle) aCirc = Handle(Geom_Circle)::DownCast(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())); - } + 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 } @@ -174,7 +242,7 @@ std::shared_ptr GeomAPI_Edge::ellipse() const return aEllipse; } } - return std::shared_ptr(); // not elipse + return std::shared_ptr(); // not ellipse } std::shared_ptr GeomAPI_Edge::line() const @@ -216,16 +284,16 @@ 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); @@ -239,12 +307,26 @@ bool GeomAPI_Edge::isEqual(const std::shared_ptr theEdge) const return true; } +void GeomAPI_Edge::setRange(const double& theFirst, const double& theLast) +{ + TopoDS_Edge anEdge = impl(); + double aFirst, aLast; + Handle(Geom_Curve) aCurve = BRep_Tool::Curve(anEdge, aFirst, aLast); + double aTolerance = BRep_Tool::Tolerance(anEdge); + if (aCurve->DynamicType() == STANDARD_TYPE(Geom_TrimmedCurve)) { + aCurve = Handle(Geom_TrimmedCurve)::DownCast(aCurve)->BasisCurve(); + BRep_Builder().UpdateEdge(anEdge, aCurve, aTolerance); + } + BRep_Builder().Range(anEdge, theFirst, theLast); +} + 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); } +// LCOV_EXCL_START bool GeomAPI_Edge::isInPlane(std::shared_ptr thePlane) const { double aFirst, aLast; @@ -264,24 +346,26 @@ bool GeomAPI_Edge::isInPlane(std::shared_ptr thePlane) const gp_Pnt aLastPnt = aCurve->Value(aLast); inPlane = aPlane.SquareDistance(aFirstPnt) < Precision::SquareConfusion() && aPlane.SquareDistance(aLastPnt) < Precision::SquareConfusion(); - } else if (aCurve->IsKind(STANDARD_TYPE(Geom_Circle))) { - // check the center on the plane and normals are collinear - Handle(Geom_Circle) aCirc = Handle(Geom_Circle)::DownCast(aCurve); - 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(); + 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 @@ -293,8 +377,9 @@ void GeomAPI_Edge::intersectWithPlane(const std::shared_ptr thePlan 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); + + // intersect the plane with the curve GeomAPI_IntCS aIntersect; aIntersect.Perform(aCurve, aPlane); if (aIntersect.IsDone() && (aIntersect.NbPoints() > 0)) { @@ -313,6 +398,19 @@ void GeomAPI_Edge::intersectWithPlane(const std::shared_ptr thePlan theResult.push_back(aPntPtr); } } + else { + // find minimal distance between the plane and the curve + GeomAPI_ExtremaCurveSurface anExtrema(aCurve, aPlane); + double aTolerance = BRep_Tool::Tolerance(TopoDS::Edge(aShape)); + if (anExtrema.NbExtrema() > 0 && + anExtrema.LowerDistance() < aTolerance) { + // distance is lower than tolerance => tangent case + gp_Pnt aPntC, aPntS; + anExtrema.NearestPoints(aPntC, aPntS); + std::shared_ptr aPntPtr(new GeomAPI_Pnt(aPntS.X(), aPntS.Y(), aPntS.Z())); + theResult.push_back(aPntPtr); + } + } } } @@ -361,3 +459,25 @@ void GeomAPI_Edge::setLastPointTolerance(const double theTolerance) 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; +}