1 // Copyright (C) 2007-2015 CEA/DEN, EDF R&D, OPEN CASCADE
3 // Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
6 // This library is free software; you can redistribute it and/or
7 // modify it under the terms of the GNU Lesser General Public
8 // License as published by the Free Software Foundation; either
9 // version 2.1 of the License, or (at your option) any later version.
11 // This library is distributed in the hope that it will be useful,
12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 // Lesser General Public License for more details.
16 // You should have received a copy of the GNU Lesser General Public
17 // License along with this library; if not, write to the Free Software
18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
23 #include <Standard_Stream.hxx>
25 #include <GEOMImpl_ProjectionDriver.hxx>
27 #include <GEOMImpl_IMirror.hxx>
28 #include <GEOMImpl_IProjection.hxx>
29 #include <GEOMImpl_IProjOnCyl.hxx>
30 #include <GEOMImpl_Types.hxx>
31 #include <GEOM_Function.hxx>
32 #include <GEOMUtils.hxx>
33 #include <GEOMUtils_HTrsfCurve2d.hxx>
35 #include <Approx_Curve2d.hxx>
36 #include <Bnd_Box2d.hxx>
37 #include <BndLib_Add2dCurve.hxx>
38 #include <BRep_Tool.hxx>
39 #include <BRepAdaptor_Curve2d.hxx>
40 #include <BRepBuilderAPI_Transform.hxx>
41 #include <BRepBuilderAPI_MakeEdge.hxx>
42 #include <BRepBuilderAPI_MakeFace.hxx>
43 #include <BRepBuilderAPI_MakeVertex.hxx>
44 #include <BRepBuilderAPI_MakeWire.hxx>
45 #include <BRepClass_FaceClassifier.hxx>
46 #include <BRepExtrema_DistShapeShape.hxx>
47 #include <BRepLib.hxx>
48 #include <BRepOffsetAPI_NormalProjection.hxx>
49 #include <BRepTools.hxx>
50 #include <BRepTools_WireExplorer.hxx>
55 #include <TopoDS_Shape.hxx>
56 #include <TopoDS_Edge.hxx>
57 #include <TopoDS_Face.hxx>
58 #include <TopoDS_Vertex.hxx>
59 #include <TopoDS_Wire.hxx>
60 #include <TopTools_ListIteratorOfListOfShape.hxx>
62 #include <GeomAPI_ProjectPointOnSurf.hxx>
63 #include <Geom_Curve.hxx>
64 #include <Geom_CylindricalSurface.hxx>
65 #include <Geom_Plane.hxx>
66 #include <Geom2d_TrimmedCurve.hxx>
68 #include <gp_Trsf.hxx>
72 //=======================================================================
75 //=======================================================================
76 const Standard_GUID& GEOMImpl_ProjectionDriver::GetID()
78 static Standard_GUID aProjectionDriver ("FF1BBB70-5D14-4df2-980B-3A668264EA16");
79 return aProjectionDriver;
83 //=======================================================================
84 //function : GEOMImpl_ProjectionDriver
86 //=======================================================================
88 GEOMImpl_ProjectionDriver::GEOMImpl_ProjectionDriver()
92 //=======================================================================
95 //=======================================================================
96 Standard_Integer GEOMImpl_ProjectionDriver::Execute(TFunction_Logbook& log) const
98 if (Label().IsNull()) return 0;
99 Handle(GEOM_Function) aFunction = GEOM_Function::GetFunction(Label());
101 if (aFunction.IsNull()) return 0;
103 Standard_Integer aType = aFunction->GetType();
105 if (aType == PROJECTION_COPY) {
110 GEOMImpl_IMirror TI (aFunction);
112 Handle(GEOM_Function) anOriginalFunction = TI.GetOriginal();
113 if (anOriginalFunction.IsNull()) return 0;
115 TopoDS_Shape anOriginal = anOriginalFunction->GetValue();
116 if (anOriginal.IsNull()) return 0;
118 // Source shape (point, edge or wire)
119 if (anOriginal.ShapeType() != TopAbs_VERTEX &&
120 anOriginal.ShapeType() != TopAbs_EDGE &&
121 anOriginal.ShapeType() != TopAbs_WIRE) {
122 Standard_ConstructionError::Raise
123 ("Projection aborted : the source shape is neither a vertex, nor an edge or a wire");
127 Handle(GEOM_Function) aTargetFunction = TI.GetPlane();
128 if (aTargetFunction.IsNull()) return 0;
129 TopoDS_Shape aFaceShape = aTargetFunction->GetValue();
130 //if (aFaceShape.IsNull() || aFaceShape.ShapeType() != TopAbs_FACE) {
131 // Standard_ConstructionError::Raise
132 // ("Projection aborted : the target shape is not a face");
135 Standard_Real tol = 1.e-4;
137 if (anOriginal.ShapeType() == TopAbs_VERTEX) {
138 if (aFaceShape.IsNull() || aFaceShape.ShapeType() != TopAbs_FACE) {
139 Standard_ConstructionError::Raise
140 ("Projection aborted : the target shape is not a face");
142 TopoDS_Face aFace = TopoDS::Face(aFaceShape);
143 Handle(Geom_Surface) surface = BRep_Tool::Surface(aFace);
144 double U1, U2, V1, V2;
145 //surface->Bounds(U1, U2, V1, V2);
146 BRepTools::UVBounds(aFace, U1, U2, V1, V2);
149 GeomAPI_ProjectPointOnSurf proj;
150 proj.Init(surface, U1, U2, V1, V2, tol);
152 gp_Pnt aPnt = BRep_Tool::Pnt(TopoDS::Vertex(anOriginal));
154 if (!proj.IsDone()) {
155 Standard_ConstructionError::Raise
156 ("Projection aborted : the algorithm failed");
158 int nbPoints = proj.NbPoints();
160 Standard_ConstructionError::Raise("No solution found");
163 Quantity_Parameter U, V;
164 proj.LowerDistanceParameters(U, V);
165 gp_Pnt2d aProjPnt (U, V);
168 BRepClass_FaceClassifier aClsf (aFace, aProjPnt, tol);
169 if (aClsf.State() != TopAbs_IN && aClsf.State() != TopAbs_ON) {
171 double minDist = RealLast();
172 for (int i = 1; i <= nbPoints; i++) {
173 Quantity_Parameter Ui, Vi;
174 proj.Parameters(i, Ui, Vi);
175 aProjPnt = gp_Pnt2d(Ui, Vi);
176 aClsf.Perform(aFace, aProjPnt, tol);
177 if (aClsf.State() == TopAbs_IN || aClsf.State() == TopAbs_ON) {
179 double dist = proj.Distance(i);
180 if (dist < minDist) {
188 Standard_ConstructionError::Raise("No solution found");
192 gp_Pnt surfPnt = surface->Value(U, V);
194 aShape = BRepBuilderAPI_MakeVertex(surfPnt).Shape();
197 //see BRepTest_BasicCommands.cxx for example of BRepOffsetAPI_NormalProjection
198 BRepOffsetAPI_NormalProjection OrtProj (aFaceShape);
199 OrtProj.Add(anOriginal);
201 // Compute maximal tolerance of projection.
202 TopExp_Explorer anExp(anOriginal,TopAbs_VERTEX);
203 Standard_Real aMaxTol = Precision::Confusion();
205 for(; anExp.More(); anExp.Next()) {
206 const TopoDS_Vertex aVtx = TopoDS::Vertex(anExp.Current());
207 const Standard_Real aCurTol = BRep_Tool::Tolerance(aVtx);
209 if (aMaxTol < aCurTol) {
214 Standard_Real tol2d = Pow(aMaxTol, 2./3);
215 GeomAbs_Shape Continuity = GeomAbs_C2;
216 Standard_Integer MaxDeg = 14;
217 Standard_Integer MaxSeg = 16;
219 OrtProj.SetParams(aMaxTol, tol2d, Continuity, MaxDeg, MaxSeg);
223 } catch (Standard_Failure) {
224 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
225 TCollection_AsciiString aMsg (aFail->GetMessageString());
227 aMsg = "Projection aborted : possibly the source shape intersects the cylinder's axis";
228 Standard_ConstructionError::Raise(aMsg.ToCString());
230 if (!OrtProj.IsDone()) {
231 Standard_ConstructionError::Raise
232 ("Projection aborted : BRepOffsetAPI_NormalProjection failed");
235 aShape = OrtProj.Shape();
237 // check that the result shape is an empty compound
238 // (IPAL22905: TC650: Projection on face dialog problems)
239 if( !aShape.IsNull() && aShape.ShapeType() == TopAbs_COMPOUND )
241 TopoDS_Iterator anIter( aShape );
243 Standard_ConstructionError::Raise("Projection aborted : empty compound produced");
247 if (aShape.IsNull()) return 0;
249 aFunction->SetValue(aShape);
250 log.SetTouched(Label());
251 } else if (aType == PROJECTION_ON_WIRE) {
252 // Perform projection of point on a wire or an edge.
253 GEOMImpl_IProjection aProj (aFunction);
254 Handle(GEOM_Function) aPointFunction = aProj.GetPoint();
255 Handle(GEOM_Function) aShapeFunction = aProj.GetShape();
257 if (aPointFunction.IsNull() || aShapeFunction.IsNull()) {
261 TopoDS_Shape aPoint = aPointFunction->GetValue();
262 TopoDS_Shape aShape = aShapeFunction->GetValue();
264 if (aPoint.IsNull() || aShape.IsNull()) {
268 // Check shape types.
269 if (aPoint.ShapeType() != TopAbs_VERTEX) {
270 Standard_ConstructionError::Raise
271 ("Projection aborted : the point is not a vertex");
274 if (aShape.ShapeType() != TopAbs_EDGE &&
275 aShape.ShapeType() != TopAbs_WIRE) {
276 Standard_ConstructionError::Raise
277 ("Projection aborted : the shape is neither an edge nor a wire");
280 // Perform projection.
281 BRepExtrema_DistShapeShape aDistShSh(aPoint, aShape, Extrema_ExtFlag_MIN);
283 if (aDistShSh.IsDone() == Standard_False) {
284 Standard_ConstructionError::Raise("Projection not done");
287 Standard_Boolean hasValidSolution = Standard_False;
288 Standard_Integer aNbSolutions = aDistShSh.NbSolution();
291 Standard_Real aTolConf = BRep_Tool::Tolerance(TopoDS::Vertex(aPoint));
292 Standard_Real aTolAng = 1.e-4;
294 for (i = 1; i <= aNbSolutions; i++) {
295 Standard_Boolean isValid = Standard_False;
296 BRepExtrema_SupportType aSupportType = aDistShSh.SupportTypeShape2(i);
297 TopoDS_Shape aSupportShape = aDistShSh.SupportOnShape2(i);
299 if (aSupportType == BRepExtrema_IsOnEdge) {
300 // Minimal distance inside edge is really a projection.
301 isValid = Standard_True;
302 aDistShSh.ParOnEdgeS2(i, aParam);
303 } else if (aSupportType == BRepExtrema_IsVertex) {
304 TopExp_Explorer anExp(aShape, TopAbs_EDGE);
306 if (aDistShSh.Value() <= aTolConf) {
307 // The point lies on the shape. This means this point
308 // is really a projection.
309 for (; anExp.More() && !isValid; anExp.Next()) {
310 TopoDS_Edge aCurEdge = TopoDS::Edge(anExp.Current());
312 if (aCurEdge.IsNull() == Standard_False) {
313 TopoDS_Vertex aVtx[2];
315 TopExp::Vertices(aCurEdge, aVtx[0], aVtx[1]);
317 for (int j = 0; j < 2; j++) {
318 if (aSupportShape.IsSame(aVtx[j])) {
319 // The current edge is a projection edge.
320 isValid = Standard_True;
321 aSupportShape = aCurEdge;
322 aParam = BRep_Tool::Parameter(aVtx[j], aCurEdge);
329 // Minimal distance to vertex is not always a real projection.
330 gp_Pnt aPnt = BRep_Tool::Pnt(TopoDS::Vertex(aPoint));
331 gp_Pnt aPrjPnt = BRep_Tool::Pnt(TopoDS::Vertex(aSupportShape));
332 gp_Vec aDProjP(aPrjPnt, aPnt);
334 for (; anExp.More() && !isValid; anExp.Next()) {
335 TopoDS_Edge aCurEdge = TopoDS::Edge(anExp.Current());
337 if (aCurEdge.IsNull() == Standard_False) {
338 TopoDS_Vertex aVtx[2];
340 TopExp::Vertices(aCurEdge, aVtx[0], aVtx[1]);
342 for (int j = 0; j < 2; j++) {
343 if (aSupportShape.IsSame(aVtx[j])) {
344 // Check if the point is a projection to the current edge.
345 Standard_Real anEdgePars[2];
346 Handle(Geom_Curve) aCurve =
347 BRep_Tool::Curve(aCurEdge, anEdgePars[0], anEdgePars[1]);
351 aParam = BRep_Tool::Parameter(aVtx[j], aCurEdge);
352 aCurve->D1(aParam, aVal, aD1);
354 if (Abs(aD1.Dot(aDProjP)) <= aTolAng) {
355 // The current edge is a projection edge.
356 isValid = Standard_True;
357 aSupportShape = aCurEdge;
369 if (hasValidSolution) {
370 Standard_ConstructionError::Raise
371 ("Projection aborted : multiple solutions");
374 // Store the valid solution.
375 hasValidSolution = Standard_True;
377 // Normalize parameter.
378 TopoDS_Edge aSupportEdge = TopoDS::Edge(aSupportShape);
379 Standard_Real aF, aL;
381 BRep_Tool::Range(aSupportEdge, aF, aL);
383 if (Abs(aL - aF) <= aTolConf) {
384 Standard_ConstructionError::Raise
385 ("Projection aborted : degenerated projection edge");
388 aParam = (aParam - aF)/(aL - aF);
391 // Compute edge index.
392 TopExp_Explorer anExp(aShape, TopAbs_EDGE);
395 for (; anExp.More(); anExp.Next(), anIndex++) {
396 if (aSupportShape.IsSame(anExp.Current())) {
397 aProj.SetIndex(anIndex);
403 Standard_ConstructionError::Raise
404 ("Projection aborted : Can't define edge index");
407 // Construct a projection vertex.
408 const gp_Pnt &aPntProj = aDistShSh.PointOnShape2(i);
409 TopoDS_Shape aProj = BRepBuilderAPI_MakeVertex(aPntProj).Shape();
411 aFunction->SetValue(aProj);
415 if (!hasValidSolution) {
416 Standard_ConstructionError::Raise("Projection aborted : no projection");
418 } else if (aType == PROJECTION_ON_CYLINDER) {
419 GEOMImpl_IProjOnCyl aProj (aFunction);
420 Handle(GEOM_Function) aShapeFunction = aProj.GetShape();
422 if (aShapeFunction.IsNull()) {
426 TopoDS_Shape aShape = aShapeFunction->GetValue();
428 if (aShape.IsNull()) {
433 const TopAbs_ShapeEnum aType = aShape.ShapeType();
434 const Standard_Real aRadius = aProj.GetRadius();
435 const Standard_Real aStartAngle = aProj.GetStartAngle();
436 const Standard_Real aLengthAngle = aProj.GetAngleLength();
437 const Standard_Real aRotationAngle = aProj.GetAngleRotation();
439 if (aType != TopAbs_WIRE && aType != TopAbs_FACE) {
443 if (aRadius <= Precision::Confusion()) {
447 TopoDS_Shape aProjShape = projectOnCylinder
448 (aShape, aRadius, aStartAngle, aLengthAngle, aRotationAngle);
450 if (aProjShape.IsNull()) {
454 aFunction->SetValue(aProjShape);
460 //================================================================================
462 * \brief Returns a name of creation operation and names and values of creation parameters
464 //================================================================================
466 bool GEOMImpl_ProjectionDriver::
467 GetCreationInformation(std::string& theOperationName,
468 std::vector<GEOM_Param>& theParams)
470 if (Label().IsNull()) return 0;
471 Handle(GEOM_Function) function = GEOM_Function::GetFunction(Label());
473 Standard_Integer aType = function->GetType();
475 theOperationName = "PROJECTION";
478 case PROJECTION_COPY:
480 GEOMImpl_IMirror aCI( function );
482 AddParam( theParams, "Source object", aCI.GetOriginal() );
483 AddParam( theParams, "Target face", aCI.GetPlane() );
486 case PROJECTION_ON_WIRE:
488 GEOMImpl_IProjection aProj (function);
490 AddParam(theParams, "Point", aProj.GetPoint());
491 AddParam(theParams, "Shape", aProj.GetShape());
495 case PROJECTION_ON_CYLINDER:
497 theOperationName = "PROJ_ON_CYL";
499 GEOMImpl_IProjOnCyl aProj (function);
500 const Standard_Real aLengthAngle = aProj.GetAngleLength();
502 AddParam(theParams, "Shape", aProj.GetShape());
503 AddParam(theParams, "Radius", aProj.GetRadius());
504 AddParam(theParams, "Start angle", aProj.GetStartAngle());
506 if (aLengthAngle >= 0.) {
507 AddParam(theParams, "Length angle", aLengthAngle);
510 AddParam(theParams, "Rotation angle", aProj.GetAngleRotation());
521 //================================================================================
523 * \brief Performs projection of a planar wire or a face on a cylinder.
525 //================================================================================
527 TopoDS_Shape GEOMImpl_ProjectionDriver::projectOnCylinder
528 (const TopoDS_Shape &theShape,
529 const Standard_Real theRadius,
530 const Standard_Real theStartAngle,
531 const Standard_Real theAngleLength,
532 const Standard_Real theAngleRotation) const
534 TopoDS_Shape aResult;
537 const TopAbs_ShapeEnum aType = theShape.ShapeType();
540 if (aType == TopAbs_WIRE) {
541 // Try to create a planar face.
542 TopoDS_Wire aWire = TopoDS::Wire(theShape);
543 BRepBuilderAPI_MakeFace aMkFace(aWire, Standard_True);
545 if (aMkFace.IsDone()) {
546 aFace = aMkFace.Face();
548 // Check if the wire is a straight line.
549 TopExp_Explorer anEExp(aWire, TopAbs_EDGE);
552 for (; anEExp.More(); anEExp.Next()) {
553 anEdge = TopoDS::Edge(anEExp.Current());
555 if (!BRep_Tool::Degenerated(anEdge)) {
561 // Not degenerated edge found. Try to create a plane.
562 Standard_Real aPar[2];
563 Handle(Geom_Curve) aCurve = BRep_Tool::Curve(anEdge, aPar[0], aPar[1]);
564 gp_Pnt aP0 = aCurve->Value(aPar[0]);
565 gp_Pnt aP1 = aCurve->Value(0.5*(aPar[1] + aPar[0]));
566 gp_Vec aX(aP1.XYZ().Subtracted(aP0.XYZ()));
567 Standard_Real aTolConf = Precision::Confusion();
569 if (aX.Magnitude() > aTolConf) {
572 // Get the plane normal ortogonal to Z axis.
573 gp_Vec aZ(0., 0., 1.);
574 gp_Vec aN = aX.Crossed(aZ);
576 if (aN.Magnitude() <= aTolConf) {
577 // aX is parallel to aZ. Get the plane normal ortogonal to Y axis.
578 gp_Vec aY(0., 1., 0.);
583 if (aN.Magnitude() > aTolConf) {
584 gp_Ax3 anAxis(aP0, gp_Dir(aN), gp_Dir(aX));
585 Handle(Geom_Plane) aPlane = new Geom_Plane(anAxis);
586 BRepBuilderAPI_MakeFace aMkFace(aPlane, aWire);
588 if (aMkFace.IsDone()) {
589 aFace = aMkFace.Face();
595 } else if (aType == TopAbs_FACE) {
596 aFace = TopoDS::Face(theShape);
599 if (aFace.IsNull()) {
603 // Compute 2d translation transformation.
604 TopoDS_Wire anOuterWire = BRepTools::OuterWire(aFace);
605 BRepTools_WireExplorer aOWExp(anOuterWire, aFace);
607 if (!aOWExp.More()) {
612 // Rotate 2D presentation of face.
613 TopoDS_Vertex aFirstVertex = aOWExp.CurrentVertex();
614 TopoDS_Edge aFirstEdge = aOWExp.Current();
615 gp_Pnt aPnt = BRep_Tool::Pnt(aFirstVertex);
616 BRepAdaptor_Curve2d anAdaptorCurve(aFirstEdge, aFace);
617 Standard_Real aParam =
618 BRep_Tool::Parameter(aFirstVertex, aFirstEdge, aFace);
621 gp_Vec2d aVecU0(1., 0);
623 anAdaptorCurve.D1(aParam, aPntUV, aVecUV);
625 if (aVecUV.Magnitude() <= gp::Resolution()) {
629 if (aFirstEdge.Orientation() == TopAbs_REVERSED) {
633 const Standard_Real anAngle = aVecUV.Angle(aVecU0) + theAngleRotation;
634 const Standard_Boolean isToRotate = Abs(anAngle) > Precision::Angular();
637 Standard_Real aPar[2];
640 aRotTrsf.SetRotation(aPntUV, anAngle);
643 for (; aOWExp.More(); aOWExp.Next()) {
644 TopoDS_Edge anEdge = aOWExp.Current();
645 Handle(Geom2d_Curve) aCurve =
646 BRep_Tool::CurveOnSurface(anEdge, aFace, aPar[0], aPar[1]);
648 if (aCurve.IsNull()) {
653 aCurve = Handle(Geom2d_Curve)::DownCast(aCurve->Transformed(aRotTrsf));
656 BndLib_Add2dCurve::Add(aCurve, aPar[0], aPar[1], 0., aUVBox);
662 aUVBox.Get(aU[0], aV[0], aU[1], aV[1]);
664 // Compute anisotropic transformation from a face's 2d space
665 // to cylinder's 2d space.
666 GEOMUtils::Trsf2d aTrsf2d
667 (1./theRadius, 0., theStartAngle - aU[0]/theRadius,
668 0., 1., aPnt.Z() - aPntUV.Y());
670 // Compute scaling trsf.
671 const Standard_Boolean isToScale = theAngleLength >= Precision::Angular();
672 gp_Trsf2d aScaleTrsf;
675 // Perform 2d scaling.
676 gp_Pnt2d aMidPnt(0.5*(aU[1] + aU[0]), 0.5*(aV[1] + aV[0]));
677 const Standard_Real aScaleFactor = theAngleLength*theRadius/(aU[1] - aU[0]);
679 aTrsf2d.TransformD0(aMidPnt);
681 aScaleTrsf.SetScale(aMidPnt, aScaleFactor);
684 // Get 2d presentation of a face.
685 Handle(Geom_Surface) aCylinder =
686 new Geom_CylindricalSurface(gp_Ax3(), theRadius);
687 GeomAdaptor_Surface aGACyl(aCylinder);
688 TopExp_Explorer anExp(aFace, TopAbs_WIRE);
689 TopTools_ListOfShape aWires;
690 Standard_Real aUResol = aGACyl.UResolution(Precision::Confusion());
691 Standard_Real aVResol = aGACyl.VResolution(Precision::Confusion());
693 for (; anExp.More(); anExp.Next()) {
694 TopoDS_Wire aWire = TopoDS::Wire(anExp.Current());
695 BRepTools_WireExplorer aWExp(aWire, aFace);
696 BRepBuilderAPI_MakeWire aMkWire;
698 for (; aWExp.More(); aWExp.Next()) {
699 TopoDS_Edge anEdge = aWExp.Current();
700 Handle(Geom2d_Curve) aCurve =
701 BRep_Tool::CurveOnSurface(anEdge, aFace, aPar[0], aPar[1]);
703 if (aCurve.IsNull()) {
708 aCurve = Handle(Geom2d_Curve)::DownCast(aCurve->Transformed(aRotTrsf));
711 // Transform the curve to cylinder's parametric space.
712 GEOMUtils::Handle(HTrsfCurve2d) aTrsfCurve =
713 new GEOMUtils::HTrsfCurve2d(aCurve, aPar[0], aPar[1], aTrsf2d);
714 Approx_Curve2d aConv (aTrsfCurve, aPar[0], aPar[1],
715 aUResol, aVResol, GeomAbs_C1,
718 if (!aConv.IsDone() && !aConv.HasResult()) {
722 Handle(Geom2d_Curve) aCylCurve = aConv.Curve();
725 aCylCurve->Transform(aScaleTrsf);
728 // Create edge and add it to the wire.
729 BRepBuilderAPI_MakeEdge aMkEdge(aCylCurve, aCylinder);
731 if (!aMkEdge.IsDone()) {
735 aMkWire.Add(aMkEdge.Edge());
737 if (!aMkWire.IsDone()) {
742 if (aWire.IsSame(anOuterWire)) {
743 // Make the outer wire first.
744 aWires.Prepend(aMkWire.Wire());
746 aWires.Append(aMkWire.Wire());
751 if (aWires.IsEmpty()) {
755 TopTools_ListIteratorOfListOfShape aWIter(aWires);
756 TopoDS_Wire aWire = TopoDS::Wire(aWIter.Value());
757 BRepBuilderAPI_MakeFace aMkFace(aCylinder, aWire);
759 if (!aMkFace.IsDone()) {
763 for (aWIter.Next(); aWIter.More(); aWIter.Next()) {
764 aWire = TopoDS::Wire(aWIter.Value());
767 if (!aMkFace.IsDone()) {
773 TopoDS_Face aCylFace = aMkFace.Face();
774 TopoDS_Shape aResShape;
776 BRepLib::BuildCurves3d(aCylFace);
779 if (aType == TopAbs_WIRE) {
780 TopExp_Explorer aResExp(aCylFace, TopAbs_WIRE);
782 if (aResExp.More()) {
783 aResShape = aResExp.Current();
786 aResShape = aCylFace;
789 if (aResShape.IsNull() == Standard_False) {
790 if (!GEOMUtils::CheckShape(aResShape, true)) {
791 if (!GEOMUtils::FixShapeTolerance(aResShape)) {
802 IMPLEMENT_STANDARD_HANDLE (GEOMImpl_ProjectionDriver,GEOM_BaseDriver);
803 IMPLEMENT_STANDARD_RTTIEXT (GEOMImpl_ProjectionDriver,GEOM_BaseDriver);