X-Git-Url: http://git.salome-platform.org/gitweb/?a=blobdiff_plain;f=src%2FGeomAlgoAPI%2FGeomAlgoAPI_PointBuilder.cpp;h=7b46a50ff1ed6737d00aa27467175433abe0bdc6;hb=6e421e939851e0de46554ae45a3ca0e1f67cd91d;hp=086cc419b429c43205fea01aba229792954b2955;hpb=9e869ede4d8c56262bb20534543c2bf56cd6a91b;p=modules%2Fshaper.git diff --git a/src/GeomAlgoAPI/GeomAlgoAPI_PointBuilder.cpp b/src/GeomAlgoAPI/GeomAlgoAPI_PointBuilder.cpp index 086cc419b..7b46a50ff 100644 --- a/src/GeomAlgoAPI/GeomAlgoAPI_PointBuilder.cpp +++ b/src/GeomAlgoAPI/GeomAlgoAPI_PointBuilder.cpp @@ -1,20 +1,196 @@ -// File: GeomAlgoAPI_PointBuilder.cpp -// Created: 02 Jun 2014 -// Author: Mikhail PONIKAROV +// 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 "GeomAlgoAPI_PointBuilder.h" + +#include +#include +#include +#include #include #include +#include + +#include #include +#include +#include +#include +#include +#include +#include +#include +#include +#include #include -std::shared_ptr GeomAlgoAPI_PointBuilder::point( - std::shared_ptr thePoint) +//================================================================================================== +std::shared_ptr + GeomAlgoAPI_PointBuilder::vertex(const std::shared_ptr thePoint) { const gp_Pnt& aPnt = thePoint->impl(); BRepBuilderAPI_MakeVertex aMaker(aPnt); TopoDS_Vertex aVertex = aMaker.Vertex(); - std::shared_ptr aRes(new GeomAPI_Shape); + std::shared_ptr aRes(new GeomAPI_Vertex); + aRes->setImpl(new TopoDS_Shape(aVertex)); + return aRes; +} + +//================================================================================================== +std::shared_ptr GeomAlgoAPI_PointBuilder::vertex(const double theX, + const double theY, + const double theZ) +{ + const gp_Pnt aPnt(theX, theY, theZ); + BRepBuilderAPI_MakeVertex aMaker(aPnt); + TopoDS_Vertex aVertex = aMaker.Vertex(); + std::shared_ptr aRes(new GeomAPI_Vertex); aRes->setImpl(new TopoDS_Shape(aVertex)); return aRes; } + +//================================================================================================== +std::shared_ptr + GeomAlgoAPI_PointBuilder::point(const std::shared_ptr theVertex) +{ + TopoDS_Shape aShape = theVertex->impl(); + if ((!aShape.IsNull()) && (aShape.ShapeType() == TopAbs_VERTEX)) { + TopoDS_Vertex aVertex = TopoDS::Vertex(aShape); + gp_Pnt aPoint = BRep_Tool::Pnt(aVertex); + std::shared_ptr aPnt(new GeomAPI_Pnt(aPoint.X(), aPoint.Y(), aPoint.Z())); + return aPnt; + } + return std::shared_ptr(); +} + +//================================================================================================== +std::shared_ptr GeomAlgoAPI_PointBuilder::vertexOnEdge( + const std::shared_ptr theEdge, + const double theValue, + const bool theIsPercent, + const bool theIsReverse) +{ + std::shared_ptr aVertex; + + if(!theEdge.get()) { + return aVertex; + } + + double aValue = theValue; + if(theIsPercent) { + aValue = theEdge->length() / 100.0 * aValue; + } + + const TopoDS_Edge& anEdge = TopoDS::Edge(theEdge->impl()); + Standard_Real aUFirst, aULast; + Handle(Geom_Curve) anEdgeCurve = BRep_Tool::Curve(anEdge, aUFirst, aULast); + + if(!anEdgeCurve.IsNull() ) { + Handle(Geom_Curve) aReOrientedCurve = anEdgeCurve; + + if(theIsReverse) { + aReOrientedCurve = anEdgeCurve->Reversed(); + aUFirst = anEdgeCurve->ReversedParameter(aULast); + } + + // Get the point by length + GeomAdaptor_Curve anAdapCurve = GeomAdaptor_Curve(aReOrientedCurve); + GCPnts_AbscissaPoint anAbsPnt(anAdapCurve, aValue, aUFirst); + Standard_Real aParam = anAbsPnt.Parameter(); + gp_Pnt aPnt = anAdapCurve.Value(aParam); + BRepBuilderAPI_MakeVertex aMkVertex(aPnt); + const TopoDS_Vertex& aShape = aMkVertex.Vertex(); + aVertex.reset(new GeomAPI_Vertex()); + aVertex->setImpl(new TopoDS_Vertex(aShape)); + } + + return aVertex; +} + +//================================================================================================== +std::shared_ptr GeomAlgoAPI_PointBuilder::vertexByProjection( + const std::shared_ptr theVertex, + const std::shared_ptr theEdge) +{ + std::shared_ptr aVertex; + + if (!theVertex.get() || !theEdge.get()) { + return aVertex; + } + + std::shared_ptr aProjPnt = theVertex->point(); + gp_Pnt aPnt(aProjPnt->x(), aProjPnt->y(), aProjPnt->z()); + + TopoDS_Edge anEdge = TopoDS::Edge(theEdge->impl()); + double aFirstOnCurve, aLastOnCurve; + Handle(Geom_Curve) aCurve = BRep_Tool::Curve(anEdge, aFirstOnCurve, aLastOnCurve); + + if (aCurve.IsNull()) { + return aVertex; + } + + GeomAPI_ProjectPointOnCurve aProjection(aPnt, aCurve); + + if (aProjection.NbPoints() != 0) { + gp_Pnt aNearestPoint = aProjection.NearestPoint(); + aVertex.reset(new GeomAPI_Vertex(aNearestPoint.X(), aNearestPoint.Y(), aNearestPoint.Z())); + } + + return aVertex; +} + +//================================================================================================== +std::shared_ptr GeomAlgoAPI_PointBuilder::vertexByProjection( + const std::shared_ptr theVertex, + const std::shared_ptr theFace) +{ + std::shared_ptr aVertex; + + if (theVertex.get() && theFace.get() && theFace->isPlanar()) { + std::shared_ptr aProjPnt = theVertex->point(); + std::shared_ptr aProjPln = theFace->getPlane(); + + std::shared_ptr aPnt = aProjPln->project(aProjPnt); + + if (aPnt.get()) + aVertex.reset(new GeomAPI_Vertex(aPnt->x(), aPnt->y(), aPnt->z())); + } + + return aVertex; +} + +//================================================================================================== +std::shared_ptr GeomAlgoAPI_PointBuilder::vertexByIntersection( + const std::shared_ptr theEdge1, + const std::shared_ptr theEdge2) +{ + std::shared_ptr aVertex; + + if (theEdge1.get() && theEdge2.get() && theEdge1->isLine() && theEdge2->isLine()) { + std::shared_ptr aLin1 = theEdge1->line(); + std::shared_ptr aLin2 = theEdge2->line(); + + std::shared_ptr aPnt = aLin1->intersect(aLin2); + + if (aPnt.get()) + aVertex.reset(new GeomAPI_Vertex(aPnt->x(), aPnt->y(), aPnt->z())); + } + + return aVertex; +}