1 // Copyright (C) 2007-2011 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.
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
22 #include <Standard_Stream.hxx>
24 #include <GEOMImpl_IMeasureOperations.hxx>
26 #include <GEOMImpl_Types.hxx>
27 #include <GEOMImpl_MeasureDriver.hxx>
28 #include <GEOMImpl_IMeasure.hxx>
29 #include <GEOMImpl_IShapesOperations.hxx>
31 #include <GEOMAlgo_ShapeInfo.hxx>
32 #include <GEOMAlgo_ShapeInfoFiller.hxx>
34 #include <GEOM_Function.hxx>
35 #include <GEOM_PythonDump.hxx>
37 #include <NMTTools_CheckerSI.hxx>
39 #include <NMTDS_Tools.hxx>
40 #include <NMTDS_InterfPool.hxx>
41 #include <NMTDS_PInterfPool.hxx>
42 //#include <NMTDS_PassKeyBoolean.hxx>
43 #include <NMTDS_PairBoolean.hxx>
44 #include <NMTDS_ShapesDataStructure.hxx>
45 //#include <NMTDS_ListIteratorOfListOfPassKeyBoolean.hxx>
46 #include <NMTDS_ListIteratorOfListOfPairBoolean.hxx>
48 #include <Basics_OCCTVersion.hxx>
50 #include <utilities.h>
52 #include <Utils_ExceptHandlers.hxx>
55 #include <TFunction_DriverTable.hxx>
56 #include <TFunction_Driver.hxx>
57 #include <TFunction_Logbook.hxx>
58 #include <TDF_Tool.hxx>
60 #include <BRep_Tool.hxx>
61 #include <BRepAdaptor_Surface.hxx>
62 #include <BRepBndLib.hxx>
63 #include <BRepCheck.hxx>
64 #include <BRepCheck_ListIteratorOfListOfStatus.hxx>
65 #include <BRepCheck_Result.hxx>
66 #include <BRepCheck_Shell.hxx>
67 #include <BRepExtrema_DistShapeShape.hxx>
68 #include <BRepGProp.hxx>
69 #include <BRepTools.hxx>
71 #include <Bnd_Box.hxx>
73 #include <GProp_GProps.hxx>
74 #include <GProp_PrincipalProps.hxx>
79 #include <TopoDS_Edge.hxx>
80 #include <TopoDS_Face.hxx>
81 #include <TopoDS_Shape.hxx>
82 #include <TopoDS_Vertex.hxx>
83 #include <TopoDS_Iterator.hxx>
84 #include <TopExp_Explorer.hxx>
85 #include <TopTools_MapOfShape.hxx>
86 #include <TopTools_ListOfShape.hxx>
87 #include <TopTools_ListIteratorOfListOfShape.hxx>
89 #include <GeomAbs_SurfaceType.hxx>
90 #include <Geom_Surface.hxx>
91 #include <Geom_Plane.hxx>
92 #include <Geom_SphericalSurface.hxx>
93 #include <Geom_CylindricalSurface.hxx>
94 #include <Geom_ToroidalSurface.hxx>
95 #include <Geom_ConicalSurface.hxx>
96 #include <Geom_SurfaceOfLinearExtrusion.hxx>
97 #include <Geom_SurfaceOfRevolution.hxx>
98 #include <Geom_BezierSurface.hxx>
99 #include <Geom_BSplineSurface.hxx>
100 #include <Geom_RectangularTrimmedSurface.hxx>
101 #include <Geom_OffsetSurface.hxx>
102 #include <Geom_Line.hxx>
104 #include <gp_Pln.hxx>
105 #include <gp_Lin.hxx>
107 #include <GeomAPI_ProjectPointOnCurve.hxx>
108 #include <GeomLProp_CLProps.hxx>
109 #include <GeomLProp_SLProps.hxx>
110 #include <ShapeAnalysis.hxx>
111 #include <ShapeAnalysis_Surface.hxx>
113 #include <Standard_Failure.hxx>
114 #include <Standard_ErrorHandler.hxx> // CAREFUL ! position of this file is critic : see Lucien PIGNOLONI / OCC
116 #include <BRepClass3d_SolidClassifier.hxx>
117 #include <BRep_Builder.hxx>
118 #include <GeomAPI_IntSS.hxx>
119 #include <Geom_Circle.hxx>
120 #include <Geom_SphericalSurface.hxx>
121 #include <Geom_ToroidalSurface.hxx>
122 #include <ShapeFix_Shape.hxx>
123 #include <TopoDS_Compound.hxx>
126 //=============================================================================
130 //=============================================================================
131 GEOMImpl_IMeasureOperations::GEOMImpl_IMeasureOperations (GEOM_Engine* theEngine, int theDocID)
132 : GEOM_IOperations(theEngine, theDocID)
134 MESSAGE("GEOMImpl_IMeasureOperations::GEOMImpl_IMeasureOperations");
137 //=============================================================================
141 //=============================================================================
142 GEOMImpl_IMeasureOperations::~GEOMImpl_IMeasureOperations()
144 MESSAGE("GEOMImpl_IMeasureOperations::~GEOMImpl_IMeasureOperations");
147 //=============================================================================
148 /*! Get kind and parameters of the given shape.
150 //=============================================================================
151 GEOMImpl_IMeasureOperations::ShapeKind GEOMImpl_IMeasureOperations::KindOfShape
152 (Handle(GEOM_Object) theShape,
153 Handle(TColStd_HSequenceOfInteger)& theIntegers,
154 Handle(TColStd_HSequenceOfReal)& theDoubles)
157 ShapeKind aKind = SK_NO_SHAPE;
159 if (theIntegers.IsNull()) theIntegers = new TColStd_HSequenceOfInteger;
160 else theIntegers->Clear();
162 if (theDoubles.IsNull()) theDoubles = new TColStd_HSequenceOfReal;
163 else theDoubles->Clear();
165 if (theShape.IsNull())
168 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
169 if (aRefShape.IsNull()) return aKind;
171 TopoDS_Shape aShape = aRefShape->GetValue();
172 if (aShape.IsNull()) return aKind;
174 int geom_type = theShape->GetType();
176 // check if it's advanced shape
177 if ( geom_type > ADVANCED_BASE ) {
183 GEOMAlgo_ShapeInfoFiller aSF;
184 aSF.SetShape(aShape);
186 Standard_Integer iErr = aSF.ErrorStatus();
188 SetErrorCode("Error in GEOMAlgo_ShapeInfoFiller");
191 const GEOMAlgo_ShapeInfo& anInfo = aSF.Info();
193 // Interprete results
194 TopAbs_ShapeEnum aType = anInfo.Type();
197 case TopAbs_COMPOUND:
198 case TopAbs_COMPSOLID:
200 // (+) geompy.kind.COMPOUND nb_solids nb_faces nb_edges nb_vertices
201 // (+) geompy.kind.COMPSOLID nb_solids nb_faces nb_edges nb_vertices
202 // ??? "nb_faces" - all faces or only 'standalone' faces?
203 if (aType == TopAbs_COMPOUND)
206 aKind = SK_COMPSOLID;
208 //theIntegers->Append(anInfo.NbSubShapes(TopAbs_COMPOUND));
209 //theIntegers->Append(anInfo.NbSubShapes(TopAbs_COMPSOLID));
210 theIntegers->Append(anInfo.NbSubShapes(TopAbs_SOLID));
211 theIntegers->Append(anInfo.NbSubShapes(TopAbs_FACE));
212 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
213 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
219 // (+) geompy.kind.SHELL geompy.info.closed nb_faces nb_edges nb_vertices
220 // (+) geompy.kind.SHELL geompy.info.unclosed nb_faces nb_edges nb_vertices
223 theIntegers->Append((int)anInfo.KindOfClosed());
225 theIntegers->Append(anInfo.NbSubShapes(TopAbs_FACE));
226 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
227 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
233 // (+) geompy.kind.WIRE geompy.info.closed nb_edges nb_vertices
234 // (+) geompy.kind.WIRE geompy.info.unclosed nb_edges nb_vertices
237 theIntegers->Append((int)anInfo.KindOfClosed());
239 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
240 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
248 GEOMAlgo_KindOfName aKN = anInfo.KindOfName();
251 case GEOMAlgo_KN_SPHERE:
252 // (+) geompy.kind.SPHERE xc yc zc R
256 gp_Pnt aC = anInfo.Location();
257 theDoubles->Append(aC.X());
258 theDoubles->Append(aC.Y());
259 theDoubles->Append(aC.Z());
261 theDoubles->Append(anInfo.Radius1());
264 case GEOMAlgo_KN_CYLINDER:
265 // (+) geompy.kind.CYLINDER xb yb zb dx dy dz R H
269 gp_Pnt aC = anInfo.Location();
270 theDoubles->Append(aC.X());
271 theDoubles->Append(aC.Y());
272 theDoubles->Append(aC.Z());
274 gp_Ax3 anAx3 = anInfo.Position();
275 gp_Dir aD = anAx3.Direction();
276 theDoubles->Append(aD.X());
277 theDoubles->Append(aD.Y());
278 theDoubles->Append(aD.Z());
280 theDoubles->Append(anInfo.Radius1());
281 theDoubles->Append(anInfo.Height());
284 case GEOMAlgo_KN_BOX:
285 // (+) geompy.kind.BOX xc yc zc ax ay az
289 gp_Pnt aC = anInfo.Location();
290 theDoubles->Append(aC.X());
291 theDoubles->Append(aC.Y());
292 theDoubles->Append(aC.Z());
294 gp_Ax3 anAx3 = anInfo.Position();
295 gp_Dir aD = anAx3.Direction();
296 gp_Dir aX = anAx3.XDirection();
299 if (aD.IsParallel(gp::DZ(), Precision::Angular()) &&
300 aX.IsParallel(gp::DX(), Precision::Angular())) {
301 theDoubles->Append(anInfo.Length()); // ax'
302 theDoubles->Append(anInfo.Width()); // ay'
303 theDoubles->Append(anInfo.Height()); // az'
305 else if (aD.IsParallel(gp::DZ(), Precision::Angular()) &&
306 aX.IsParallel(gp::DY(), Precision::Angular())) {
307 theDoubles->Append(anInfo.Width()); // ay'
308 theDoubles->Append(anInfo.Length()); // ax'
309 theDoubles->Append(anInfo.Height()); // az'
311 else if (aD.IsParallel(gp::DX(), Precision::Angular()) &&
312 aX.IsParallel(gp::DZ(), Precision::Angular())) {
313 theDoubles->Append(anInfo.Height()); // az'
314 theDoubles->Append(anInfo.Width()); // ay'
315 theDoubles->Append(anInfo.Length()); // ax'
317 else if (aD.IsParallel(gp::DX(), Precision::Angular()) &&
318 aX.IsParallel(gp::DY(), Precision::Angular())) {
319 theDoubles->Append(anInfo.Height()); // az'
320 theDoubles->Append(anInfo.Length()); // ax'
321 theDoubles->Append(anInfo.Width()); // ay'
323 else if (aD.IsParallel(gp::DY(), Precision::Angular()) &&
324 aX.IsParallel(gp::DZ(), Precision::Angular())) {
325 theDoubles->Append(anInfo.Width()); // ay'
326 theDoubles->Append(anInfo.Height()); // az'
327 theDoubles->Append(anInfo.Length()); // ax'
329 else if (aD.IsParallel(gp::DY(), Precision::Angular()) &&
330 aX.IsParallel(gp::DX(), Precision::Angular())) {
331 theDoubles->Append(anInfo.Length()); // ax'
332 theDoubles->Append(anInfo.Height()); // az'
333 theDoubles->Append(anInfo.Width()); // ay'
336 // (+) geompy.kind.ROTATED_BOX xo yo zo zx zy zz xx xy xz ax ay az
337 aKind = SK_ROTATED_BOX;
339 // Direction and XDirection
340 theDoubles->Append(aD.X());
341 theDoubles->Append(aD.Y());
342 theDoubles->Append(aD.Z());
344 theDoubles->Append(aX.X());
345 theDoubles->Append(aX.Y());
346 theDoubles->Append(aX.Z());
349 theDoubles->Append(anInfo.Length());
350 theDoubles->Append(anInfo.Width());
351 theDoubles->Append(anInfo.Height());
355 case GEOMAlgo_KN_TORUS:
356 // (+) geompy.kind.TORUS xc yc zc dx dy dz R_1 R_2
360 gp_Pnt aO = anInfo.Location();
361 theDoubles->Append(aO.X());
362 theDoubles->Append(aO.Y());
363 theDoubles->Append(aO.Z());
365 gp_Ax3 anAx3 = anInfo.Position();
366 gp_Dir aD = anAx3.Direction();
367 theDoubles->Append(aD.X());
368 theDoubles->Append(aD.Y());
369 theDoubles->Append(aD.Z());
371 theDoubles->Append(anInfo.Radius1());
372 theDoubles->Append(anInfo.Radius2());
375 case GEOMAlgo_KN_CONE:
376 // (+) geompy.kind.CONE xb yb zb dx dy dz R_1 R_2 H
380 gp_Pnt aO = anInfo.Location();
381 theDoubles->Append(aO.X());
382 theDoubles->Append(aO.Y());
383 theDoubles->Append(aO.Z());
385 gp_Ax3 anAx3 = anInfo.Position();
386 gp_Dir aD = anAx3.Direction();
387 theDoubles->Append(aD.X());
388 theDoubles->Append(aD.Y());
389 theDoubles->Append(aD.Z());
391 theDoubles->Append(anInfo.Radius1());
392 theDoubles->Append(anInfo.Radius2());
393 theDoubles->Append(anInfo.Height());
396 case GEOMAlgo_KN_POLYHEDRON:
397 // (+) geompy.kind.POLYHEDRON nb_faces nb_edges nb_vertices
399 aKind = SK_POLYHEDRON;
401 theIntegers->Append(anInfo.NbSubShapes(TopAbs_FACE));
402 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
403 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
407 // (+) geompy.kind.SOLID nb_faces nb_edges nb_vertices
409 theIntegers->Append(anInfo.NbSubShapes(TopAbs_FACE));
410 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
411 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
421 GEOMAlgo_KindOfName aKN = anInfo.KindOfName();
423 case GEOMAlgo_KN_SPHERE:
424 // (+) geompy.kind.SPHERE2D xc yc zc R
428 gp_Pnt aC = anInfo.Location();
429 theDoubles->Append(aC.X());
430 theDoubles->Append(aC.Y());
431 theDoubles->Append(aC.Z());
433 theDoubles->Append(anInfo.Radius1());
436 case GEOMAlgo_KN_CYLINDER:
437 // (+) geompy.kind.CYLINDER2D xb yb zb dx dy dz R H
439 aKind = SK_CYLINDER2D;
441 gp_Pnt aO = anInfo.Location();
442 theDoubles->Append(aO.X());
443 theDoubles->Append(aO.Y());
444 theDoubles->Append(aO.Z());
446 gp_Ax3 anAx3 = anInfo.Position();
447 gp_Dir aD = anAx3.Direction();
448 theDoubles->Append(aD.X());
449 theDoubles->Append(aD.Y());
450 theDoubles->Append(aD.Z());
452 theDoubles->Append(anInfo.Radius1());
453 theDoubles->Append(anInfo.Height());
456 case GEOMAlgo_KN_TORUS:
457 // (+) geompy.kind.TORUS2D xc yc zc dx dy dz R_1 R_2
461 gp_Pnt aO = anInfo.Location();
462 theDoubles->Append(aO.X());
463 theDoubles->Append(aO.Y());
464 theDoubles->Append(aO.Z());
466 gp_Ax3 anAx3 = anInfo.Position();
467 gp_Dir aD = anAx3.Direction();
468 theDoubles->Append(aD.X());
469 theDoubles->Append(aD.Y());
470 theDoubles->Append(aD.Z());
472 theDoubles->Append(anInfo.Radius1());
473 theDoubles->Append(anInfo.Radius2());
476 case GEOMAlgo_KN_CONE:
477 // (+) geompy.kind.CONE2D xc yc zc dx dy dz R_1 R_2 H
481 gp_Pnt aO = anInfo.Location();
482 theDoubles->Append(aO.X());
483 theDoubles->Append(aO.Y());
484 theDoubles->Append(aO.Z());
486 gp_Ax3 anAx3 = anInfo.Position();
487 gp_Dir aD = anAx3.Direction();
488 theDoubles->Append(aD.X());
489 theDoubles->Append(aD.Y());
490 theDoubles->Append(aD.Z());
492 theDoubles->Append(anInfo.Radius1());
493 theDoubles->Append(anInfo.Radius2());
494 theDoubles->Append(anInfo.Height());
497 case GEOMAlgo_KN_DISKCIRCLE:
498 // (+) geompy.kind.DISK_CIRCLE xc yc zc dx dy dz R
500 aKind = SK_DISK_CIRCLE;
502 gp_Pnt aC = anInfo.Location();
503 theDoubles->Append(aC.X());
504 theDoubles->Append(aC.Y());
505 theDoubles->Append(aC.Z());
507 gp_Ax3 anAx3 = anInfo.Position();
508 gp_Dir aD = anAx3.Direction();
509 theDoubles->Append(aD.X());
510 theDoubles->Append(aD.Y());
511 theDoubles->Append(aD.Z());
513 theDoubles->Append(anInfo.Radius1());
516 case GEOMAlgo_KN_DISKELLIPSE:
517 // (+) geompy.kind.DISK_ELLIPSE xc yc zc dx dy dz R_1 R_2
519 aKind = SK_DISK_ELLIPSE;
521 gp_Pnt aC = anInfo.Location();
522 theDoubles->Append(aC.X());
523 theDoubles->Append(aC.Y());
524 theDoubles->Append(aC.Z());
526 gp_Ax3 anAx3 = anInfo.Position();
527 gp_Dir aD = anAx3.Direction();
528 theDoubles->Append(aD.X());
529 theDoubles->Append(aD.Y());
530 theDoubles->Append(aD.Z());
532 theDoubles->Append(anInfo.Radius1());
533 theDoubles->Append(anInfo.Radius2());
536 case GEOMAlgo_KN_RECTANGLE:
537 case GEOMAlgo_KN_TRIANGLE:
538 case GEOMAlgo_KN_QUADRANGLE:
539 case GEOMAlgo_KN_POLYGON:
540 // (+) geompy.kind.POLYGON xo yo zo dx dy dz nb_edges nb_vertices
544 gp_Pnt aO = anInfo.Location();
545 theDoubles->Append(aO.X());
546 theDoubles->Append(aO.Y());
547 theDoubles->Append(aO.Z());
549 gp_Ax3 anAx3 = anInfo.Position();
550 gp_Dir aD = anAx3.Direction();
551 theDoubles->Append(aD.X());
552 theDoubles->Append(aD.Y());
553 theDoubles->Append(aD.Z());
555 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
556 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
559 case GEOMAlgo_KN_PLANE: // infinite
560 // (+) geompy.kind.PLANE xo yo zo dx dy dz
564 gp_Pnt aC = anInfo.Location();
565 theDoubles->Append(aC.X());
566 theDoubles->Append(aC.Y());
567 theDoubles->Append(aC.Z());
569 gp_Ax3 anAx3 = anInfo.Position();
570 gp_Dir aD = anAx3.Direction();
571 theDoubles->Append(aD.X());
572 theDoubles->Append(aD.Y());
573 theDoubles->Append(aD.Z());
577 if (anInfo.KindOfShape() == GEOMAlgo_KS_PLANE) {
578 // (+) geompy.kind.PLANAR xo yo zo dx dy dz nb_edges nb_vertices
582 gp_Pnt aC = anInfo.Location();
583 theDoubles->Append(aC.X());
584 theDoubles->Append(aC.Y());
585 theDoubles->Append(aC.Z());
587 gp_Ax3 anAx3 = anInfo.Position();
588 gp_Dir aD = anAx3.Direction();
589 theDoubles->Append(aD.X());
590 theDoubles->Append(aD.Y());
591 theDoubles->Append(aD.Z());
593 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
594 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
597 // ??? geompy.kind.FACE nb_edges nb_vertices _surface_type_id_
598 // (+) geompy.kind.FACE nb_edges nb_vertices
600 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
601 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
611 GEOMAlgo_KindOfName aKN = anInfo.KindOfName();
613 case GEOMAlgo_KN_CIRCLE:
615 // (+) geompy.kind.CIRCLE xc yc zc dx dy dz R
618 gp_Pnt aC = anInfo.Location();
619 theDoubles->Append(aC.X());
620 theDoubles->Append(aC.Y());
621 theDoubles->Append(aC.Z());
623 gp_Ax3 anAx3 = anInfo.Position();
624 gp_Dir aD = anAx3.Direction();
625 theDoubles->Append(aD.X());
626 theDoubles->Append(aD.Y());
627 theDoubles->Append(aD.Z());
629 theDoubles->Append(anInfo.Radius1());
632 case GEOMAlgo_KN_ARCCIRCLE:
634 // (+) geompy.kind.ARC_CIRCLE xc yc zc dx dy dz R x1 y1 z1 x2 y2 z2
635 aKind = SK_ARC_CIRCLE;
637 gp_Pnt aC = anInfo.Location();
638 theDoubles->Append(aC.X());
639 theDoubles->Append(aC.Y());
640 theDoubles->Append(aC.Z());
642 gp_Ax3 anAx3 = anInfo.Position();
643 gp_Dir aD = anAx3.Direction();
644 theDoubles->Append(aD.X());
645 theDoubles->Append(aD.Y());
646 theDoubles->Append(aD.Z());
648 theDoubles->Append(anInfo.Radius1());
650 gp_Pnt aP1 = anInfo.Pnt1();
651 theDoubles->Append(aP1.X());
652 theDoubles->Append(aP1.Y());
653 theDoubles->Append(aP1.Z());
655 gp_Pnt aP2 = anInfo.Pnt2();
656 theDoubles->Append(aP2.X());
657 theDoubles->Append(aP2.Y());
658 theDoubles->Append(aP2.Z());
661 case GEOMAlgo_KN_ELLIPSE:
663 // (+) geompy.kind.ELLIPSE xc yc zc dx dy dz R_1 R_2
666 gp_Pnt aC = anInfo.Location();
667 theDoubles->Append(aC.X());
668 theDoubles->Append(aC.Y());
669 theDoubles->Append(aC.Z());
671 gp_Ax3 anAx3 = anInfo.Position();
672 gp_Dir aD = anAx3.Direction();
673 theDoubles->Append(aD.X());
674 theDoubles->Append(aD.Y());
675 theDoubles->Append(aD.Z());
677 theDoubles->Append(anInfo.Radius1());
678 theDoubles->Append(anInfo.Radius2());
681 case GEOMAlgo_KN_ARCELLIPSE:
683 // (+) geompy.kind.ARC_ELLIPSE xc yc zc dx dy dz R_1 R_2 x1 y1 z1 x2 y2 z2
684 aKind = SK_ARC_ELLIPSE;
686 gp_Pnt aC = anInfo.Location();
687 theDoubles->Append(aC.X());
688 theDoubles->Append(aC.Y());
689 theDoubles->Append(aC.Z());
691 gp_Ax3 anAx3 = anInfo.Position();
692 gp_Dir aD = anAx3.Direction();
693 theDoubles->Append(aD.X());
694 theDoubles->Append(aD.Y());
695 theDoubles->Append(aD.Z());
697 theDoubles->Append(anInfo.Radius1());
698 theDoubles->Append(anInfo.Radius2());
700 gp_Pnt aP1 = anInfo.Pnt1();
701 theDoubles->Append(aP1.X());
702 theDoubles->Append(aP1.Y());
703 theDoubles->Append(aP1.Z());
705 gp_Pnt aP2 = anInfo.Pnt2();
706 theDoubles->Append(aP2.X());
707 theDoubles->Append(aP2.Y());
708 theDoubles->Append(aP2.Z());
711 case GEOMAlgo_KN_LINE:
713 // ??? geompy.kind.LINE x1 y1 z1 x2 y2 z2
714 // (+) geompy.kind.LINE x1 y1 z1 dx dy dz
717 gp_Pnt aO = anInfo.Location();
718 theDoubles->Append(aO.X());
719 theDoubles->Append(aO.Y());
720 theDoubles->Append(aO.Z());
722 gp_Dir aD = anInfo.Direction();
723 theDoubles->Append(aD.X());
724 theDoubles->Append(aD.Y());
725 theDoubles->Append(aD.Z());
728 case GEOMAlgo_KN_SEGMENT:
730 // (+) geompy.kind.SEGMENT x1 y1 z1 x2 y2 z2
733 gp_Pnt aP1 = anInfo.Pnt1();
734 theDoubles->Append(aP1.X());
735 theDoubles->Append(aP1.Y());
736 theDoubles->Append(aP1.Z());
738 gp_Pnt aP2 = anInfo.Pnt2();
739 theDoubles->Append(aP2.X());
740 theDoubles->Append(aP2.Y());
741 theDoubles->Append(aP2.Z());
745 // ??? geompy.kind.EDGE nb_vertices _curve_type_id_
746 // (+) geompy.kind.EDGE nb_vertices
747 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
754 // (+) geompy.kind.VERTEX x y z
757 gp_Pnt aP = anInfo.Location();
758 theDoubles->Append(aP.X());
759 theDoubles->Append(aP.Y());
760 theDoubles->Append(aP.Z());
769 //=============================================================================
770 /*! Get LCS, corresponding to the given shape.
771 * Origin of the LCS is situated at the shape's center of mass.
772 * Axes of the LCS are obtained from shape's location or,
773 * if the shape is a planar face, from position of its plane.
775 //=============================================================================
776 gp_Ax3 GEOMImpl_IMeasureOperations::GetPosition (const TopoDS_Shape& theShape)
780 if (theShape.IsNull())
784 aResult.Transform(theShape.Location().Transformation());
785 if (theShape.ShapeType() == TopAbs_FACE) {
786 Handle(Geom_Surface) aGS = BRep_Tool::Surface(TopoDS::Face(theShape));
787 if (!aGS.IsNull() && aGS->IsKind(STANDARD_TYPE(Geom_Plane))) {
788 Handle(Geom_Plane) aGPlane = Handle(Geom_Plane)::DownCast(aGS);
789 gp_Pln aPln = aGPlane->Pln();
790 aResult = aPln.Position();
791 // In case of reverse orinetation of the face invert the plane normal
792 // (the face's normal does not mathc the plane's normal in this case)
793 if(theShape.Orientation() == TopAbs_REVERSED)
795 gp_Dir Vx = aResult.XDirection();
796 gp_Dir N = aResult.Direction().Mirrored(Vx);
797 gp_Pnt P = aResult.Location();
798 aResult = gp_Ax3(P, N, Vx);
806 TopAbs_ShapeEnum aShType = theShape.ShapeType();
808 if (aShType == TopAbs_VERTEX) {
809 aPnt = BRep_Tool::Pnt(TopoDS::Vertex(theShape));
812 if (aShType == TopAbs_COMPOUND) {
813 aShType = GEOMImpl_IShapesOperations::GetTypeOfSimplePart(theShape);
816 GProp_GProps aSystem;
817 if (aShType == TopAbs_EDGE || aShType == TopAbs_WIRE)
818 BRepGProp::LinearProperties(theShape, aSystem);
819 else if (aShType == TopAbs_FACE || aShType == TopAbs_SHELL)
820 BRepGProp::SurfaceProperties(theShape, aSystem);
822 BRepGProp::VolumeProperties(theShape, aSystem);
824 aPnt = aSystem.CentreOfMass();
827 aResult.SetLocation(aPnt);
832 //=============================================================================
836 //=============================================================================
837 void GEOMImpl_IMeasureOperations::GetPosition
838 (Handle(GEOM_Object) theShape,
839 Standard_Real& Ox, Standard_Real& Oy, Standard_Real& Oz,
840 Standard_Real& Zx, Standard_Real& Zy, Standard_Real& Zz,
841 Standard_Real& Xx, Standard_Real& Xy, Standard_Real& Xz)
845 //Set default values: global CS
846 Ox = Oy = Oz = Zx = Zy = Xy = Xz = 0.;
849 if (theShape.IsNull()) return;
851 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
852 if (aRefShape.IsNull()) return;
854 TopoDS_Shape aShape = aRefShape->GetValue();
855 if (aShape.IsNull()) {
856 SetErrorCode("The Objects has NULL Shape");
861 #if OCC_VERSION_LARGE > 0x06010000
865 gp_Ax3 anAx3 = GetPosition(aShape);
867 gp_Pnt anOri = anAx3.Location();
868 gp_Dir aDirZ = anAx3.Direction();
869 gp_Dir aDirX = anAx3.XDirection();
872 anOri.Coord(Ox, Oy, Oz);
873 aDirZ.Coord(Zx, Zy, Zz);
874 aDirX.Coord(Xx, Xy, Xz);
876 catch (Standard_Failure) {
877 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
878 SetErrorCode(aFail->GetMessageString());
885 //=============================================================================
889 //=============================================================================
890 Handle(GEOM_Object) GEOMImpl_IMeasureOperations::GetCentreOfMass
891 (Handle(GEOM_Object) theShape)
895 if (theShape.IsNull()) return NULL;
897 //Add a new CentreOfMass object
898 Handle(GEOM_Object) aCDG = GetEngine()->AddObject(GetDocID(), GEOM_CDG);
900 //Add a new CentreOfMass function
901 Handle(GEOM_Function) aFunction =
902 aCDG->AddFunction(GEOMImpl_MeasureDriver::GetID(), CDG_MEASURE);
903 if (aFunction.IsNull()) return NULL;
905 //Check if the function is set correctly
906 if (aFunction->GetDriverGUID() != GEOMImpl_MeasureDriver::GetID()) return NULL;
908 GEOMImpl_IMeasure aCI (aFunction);
910 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
911 if (aRefShape.IsNull()) return NULL;
913 aCI.SetBase(aRefShape);
915 //Compute the CentreOfMass value
917 #if OCC_VERSION_LARGE > 0x06010000
920 if (!GetSolver()->ComputeFunction(aFunction)) {
921 SetErrorCode("Measure driver failed to compute centre of mass");
925 catch (Standard_Failure) {
926 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
927 SetErrorCode(aFail->GetMessageString());
931 //Make a Python command
932 GEOM::TPythonDump(aFunction) << aCDG << " = geompy.MakeCDG(" << theShape << ")";
938 //=============================================================================
942 //=============================================================================
943 Handle(GEOM_Object) GEOMImpl_IMeasureOperations::GetVertexByIndex
944 (Handle(GEOM_Object) theShape,
945 Standard_Integer theIndex)
949 if (theShape.IsNull()) return NULL;
951 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
952 if (aRefShape.IsNull()) return NULL;
954 //Add a new Vertex object
955 Handle(GEOM_Object) aVertex = GetEngine()->AddObject(GetDocID(), GEOM_POINT);
958 Handle(GEOM_Function) aFunction =
959 aVertex->AddFunction(GEOMImpl_MeasureDriver::GetID(), VERTEX_BY_INDEX);
960 if (aFunction.IsNull()) return NULL;
962 //Check if the function is set correctly
963 if (aFunction->GetDriverGUID() != GEOMImpl_MeasureDriver::GetID()) return NULL;
965 GEOMImpl_IMeasure aCI (aFunction);
966 aCI.SetBase(aRefShape);
967 aCI.SetIndex(theIndex);
971 #if OCC_VERSION_LARGE > 0x06010000
974 if (!GetSolver()->ComputeFunction(aFunction)) {
975 SetErrorCode("Vertex by index driver failed.");
979 catch (Standard_Failure) {
980 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
981 SetErrorCode(aFail->GetMessageString());
985 //Make a Python command
986 GEOM::TPythonDump(aFunction) << aVertex << " = geompy.GetVertexByIndex(" << theShape << ", " << theIndex << ")";
992 //=============================================================================
996 //=============================================================================
997 Handle(GEOM_Object) GEOMImpl_IMeasureOperations::GetNormal
998 (Handle(GEOM_Object) theFace,
999 Handle(GEOM_Object) theOptionalPoint)
1003 if (theFace.IsNull()) return NULL;
1005 //Add a new Normale object
1006 Handle(GEOM_Object) aNorm = GetEngine()->AddObject(GetDocID(), GEOM_VECTOR);
1008 //Add a new Normale function
1009 Handle(GEOM_Function) aFunction =
1010 aNorm->AddFunction(GEOMImpl_MeasureDriver::GetID(), VECTOR_FACE_NORMALE);
1011 if (aFunction.IsNull()) return NULL;
1013 //Check if the function is set correctly
1014 if (aFunction->GetDriverGUID() != GEOMImpl_MeasureDriver::GetID()) return NULL;
1016 GEOMImpl_IMeasure aCI (aFunction);
1018 Handle(GEOM_Function) aFace = theFace->GetLastFunction();
1019 if (aFace.IsNull()) return NULL;
1023 if (!theOptionalPoint.IsNull()) {
1024 Handle(GEOM_Function) anOptPnt = theOptionalPoint->GetLastFunction();
1025 aCI.SetPoint(anOptPnt);
1028 //Compute the Normale value
1030 #if OCC_VERSION_LARGE > 0x06010000
1033 if (!GetSolver()->ComputeFunction(aFunction)) {
1034 SetErrorCode("Measure driver failed to compute normake of face");
1038 catch (Standard_Failure) {
1039 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1040 SetErrorCode(aFail->GetMessageString());
1044 //Make a Python command
1045 GEOM::TPythonDump pd (aFunction);
1046 pd << aNorm << " = geompy.GetNormal(" << theFace;
1047 if (!theOptionalPoint.IsNull()) {
1048 pd << ", " << theOptionalPoint;
1056 //=============================================================================
1058 * GetBasicProperties
1060 //=============================================================================
1061 void GEOMImpl_IMeasureOperations::GetBasicProperties (Handle(GEOM_Object) theShape,
1062 Standard_Real& theLength,
1063 Standard_Real& theSurfArea,
1064 Standard_Real& theVolume)
1068 if (theShape.IsNull()) return;
1070 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1071 if (aRefShape.IsNull()) return;
1073 TopoDS_Shape aShape = aRefShape->GetValue();
1074 if (aShape.IsNull()) {
1075 SetErrorCode("The Objects has NULL Shape");
1079 //Compute the parameters
1080 GProp_GProps LProps, SProps;
1082 #if OCC_VERSION_LARGE > 0x06010000
1085 BRepGProp::LinearProperties(aShape, LProps);
1086 theLength = LProps.Mass();
1088 BRepGProp::SurfaceProperties(aShape, SProps);
1089 theSurfArea = SProps.Mass();
1092 if (aShape.ShapeType() < TopAbs_SHELL) {
1093 for (TopExp_Explorer Exp (aShape, TopAbs_SOLID); Exp.More(); Exp.Next()) {
1094 GProp_GProps VProps;
1095 BRepGProp::VolumeProperties(Exp.Current(), VProps);
1096 theVolume += VProps.Mass();
1100 catch (Standard_Failure) {
1101 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1102 SetErrorCode(aFail->GetMessageString());
1109 //=============================================================================
1113 //=============================================================================
1114 void GEOMImpl_IMeasureOperations::GetInertia
1115 (Handle(GEOM_Object) theShape,
1116 Standard_Real& I11, Standard_Real& I12, Standard_Real& I13,
1117 Standard_Real& I21, Standard_Real& I22, Standard_Real& I23,
1118 Standard_Real& I31, Standard_Real& I32, Standard_Real& I33,
1119 Standard_Real& Ix , Standard_Real& Iy , Standard_Real& Iz)
1123 if (theShape.IsNull()) return;
1125 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1126 if (aRefShape.IsNull()) return;
1128 TopoDS_Shape aShape = aRefShape->GetValue();
1129 if (aShape.IsNull()) {
1130 SetErrorCode("The Objects has NULL Shape");
1134 //Compute the parameters
1135 GProp_GProps System;
1138 #if OCC_VERSION_LARGE > 0x06010000
1141 if (aShape.ShapeType() == TopAbs_VERTEX ||
1142 aShape.ShapeType() == TopAbs_EDGE ||
1143 aShape.ShapeType() == TopAbs_WIRE) {
1144 BRepGProp::LinearProperties(aShape, System);
1145 } else if (aShape.ShapeType() == TopAbs_FACE ||
1146 aShape.ShapeType() == TopAbs_SHELL) {
1147 BRepGProp::SurfaceProperties(aShape, System);
1149 BRepGProp::VolumeProperties(aShape, System);
1151 gp_Mat I = System.MatrixOfInertia();
1165 GProp_PrincipalProps Pr = System.PrincipalProperties();
1166 Pr.Moments(Ix,Iy,Iz);
1168 catch (Standard_Failure) {
1169 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1170 SetErrorCode(aFail->GetMessageString());
1177 //=============================================================================
1181 //=============================================================================
1182 void GEOMImpl_IMeasureOperations::GetBoundingBox
1183 (Handle(GEOM_Object) theShape,
1184 Standard_Real& Xmin, Standard_Real& Xmax,
1185 Standard_Real& Ymin, Standard_Real& Ymax,
1186 Standard_Real& Zmin, Standard_Real& Zmax)
1190 if (theShape.IsNull()) return;
1192 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1193 if (aRefShape.IsNull()) return;
1195 TopoDS_Shape aShape = aRefShape->GetValue();
1196 if (aShape.IsNull()) {
1197 SetErrorCode("The Objects has NULL Shape");
1201 //Compute the parameters
1205 #if OCC_VERSION_LARGE > 0x06010000
1208 BRepBndLib::Add(aShape, B);
1209 B.Get(Xmin, Ymin, Zmin, Xmax, Ymax, Zmax);
1211 catch (Standard_Failure) {
1212 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1213 SetErrorCode(aFail->GetMessageString());
1220 //=============================================================================
1224 //=============================================================================
1225 void GEOMImpl_IMeasureOperations::GetTolerance
1226 (Handle(GEOM_Object) theShape,
1227 Standard_Real& FaceMin, Standard_Real& FaceMax,
1228 Standard_Real& EdgeMin, Standard_Real& EdgeMax,
1229 Standard_Real& VertMin, Standard_Real& VertMax)
1233 if (theShape.IsNull()) return;
1235 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1236 if (aRefShape.IsNull()) return;
1238 TopoDS_Shape aShape = aRefShape->GetValue();
1239 if (aShape.IsNull()) {
1240 SetErrorCode("The Objects has NULL Shape");
1244 //Compute the parameters
1246 FaceMin = EdgeMin = VertMin = RealLast();
1247 FaceMax = EdgeMax = VertMax = -RealLast();
1250 #if OCC_VERSION_LARGE > 0x06010000
1253 for (TopExp_Explorer ExF (aShape, TopAbs_FACE); ExF.More(); ExF.Next()) {
1254 TopoDS_Face Face = TopoDS::Face(ExF.Current());
1255 T = BRep_Tool::Tolerance(Face);
1261 for (TopExp_Explorer ExE (aShape, TopAbs_EDGE); ExE.More(); ExE.Next()) {
1262 TopoDS_Edge Edge = TopoDS::Edge(ExE.Current());
1263 T = BRep_Tool::Tolerance(Edge);
1269 for (TopExp_Explorer ExV (aShape, TopAbs_VERTEX); ExV.More(); ExV.Next()) {
1270 TopoDS_Vertex Vertex = TopoDS::Vertex(ExV.Current());
1271 T = BRep_Tool::Tolerance(Vertex);
1278 catch (Standard_Failure) {
1279 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1280 SetErrorCode(aFail->GetMessageString());
1287 //=============================================================================
1291 //=============================================================================
1292 bool GEOMImpl_IMeasureOperations::CheckShape (Handle(GEOM_Object) theShape,
1293 const Standard_Boolean theIsCheckGeom,
1294 TCollection_AsciiString& theDump)
1298 if (theShape.IsNull()) return false;
1300 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1301 if (aRefShape.IsNull()) return false;
1303 TopoDS_Shape aShape = aRefShape->GetValue();
1304 if (aShape.IsNull()) {
1305 SetErrorCode("The Objects has NULL Shape");
1309 //Compute the parameters
1310 bool isValid = false;
1312 #if OCC_VERSION_LARGE > 0x06010000
1315 BRepCheck_Analyzer ana (aShape, theIsCheckGeom);
1316 if (ana.IsValid()) {
1320 StructuralDump(ana, aShape, theDump);
1323 catch (Standard_Failure) {
1324 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1325 SetErrorCode(aFail->GetMessageString());
1333 //=============================================================================
1335 * CheckSelfIntersections
1337 //=============================================================================
1338 bool GEOMImpl_IMeasureOperations::CheckSelfIntersections
1339 (Handle(GEOM_Object) theShape,
1340 Handle(TColStd_HSequenceOfInteger)& theIntersections)
1343 bool isGood = false;
1345 if (theIntersections.IsNull())
1346 theIntersections = new TColStd_HSequenceOfInteger;
1348 theIntersections->Clear();
1350 if (theShape.IsNull())
1353 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1354 if (aRefShape.IsNull()) return isGood;
1356 TopoDS_Shape aShape = aRefShape->GetValue();
1357 if (aShape.IsNull()) return isGood;
1361 TopoDS_Compound aCS;
1362 TopoDS_Shape aScopy;
1363 NMTDS_Tools::CopyShape(aShape, aScopy);
1365 // Map sub-shapes and their indices
1366 TopTools_IndexedMapOfShape anIndices;
1367 TopExp::MapShapes(aScopy, anIndices);
1369 aBB.MakeCompound(aCS);
1370 aBB.Add(aCS, aScopy);
1372 NMTTools_CheckerSI aCSI; // checker of self-interferences
1373 aCSI.SetCompositeShape(aCS);
1375 // 1. Launch the checker
1377 Standard_Integer iErr = aCSI.StopStatus();
1379 return false; // Error
1384 // 2. Take the shapes from DS
1385 const NMTDS_ShapesDataStructure& aDS = *(aCSI.DS());
1386 Standard_Integer aNbS = aDS.NumberOfShapesOfTheObject();
1388 // 3. Get the pairs of interfered shapes
1389 NMTDS_PInterfPool pIP = aCSI.IP();
1390 //const NMTDS_ListOfPassKeyBoolean& aLPKB = pIP->Get();
1391 const NMTDS_ListOfPairBoolean& aLPKB = pIP->Get();
1393 Standard_Integer n1, n2;
1394 //NMTDS_ListIteratorOfListOfPassKeyBoolean aIt;
1395 NMTDS_ListIteratorOfListOfPairBoolean aIt;
1397 aIt.Initialize(aLPKB);
1398 for (; aIt.More(); aIt.Next()) {
1399 //const NMTDS_PassKeyBoolean& aPKB = aIt.Value();
1400 const NMTDS_PairBoolean& aPKB = aIt.Value();
1403 if (n1 > aNbS || n2 > aNbS)
1404 return false; // Error
1406 const TopoDS_Shape& aS1 = aDS.Shape(n1);
1407 const TopoDS_Shape& aS2 = aDS.Shape(n2);
1409 theIntersections->Append(anIndices.FindIndex(aS1));
1410 theIntersections->Append(anIndices.FindIndex(aS2));
1418 //=============================================================================
1422 //=============================================================================
1423 TCollection_AsciiString GEOMImpl_IMeasureOperations::IsGoodForSolid (Handle(GEOM_Object) theShape)
1427 TCollection_AsciiString aRes = "";
1429 if (theShape.IsNull()) {
1430 aRes = "WRN_NULL_OBJECT_OR_SHAPE";
1433 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1434 if (aRefShape.IsNull()) {
1435 aRes = "WRN_NULL_OBJECT_OR_SHAPE";
1438 TopoDS_Shape aShape = aRefShape->GetValue();
1439 if (aShape.IsNull()) {
1440 aRes = "WRN_NULL_OBJECT_OR_SHAPE";
1443 if (aShape.ShapeType() == TopAbs_COMPOUND) {
1444 TopoDS_Iterator It (aShape, Standard_True, Standard_True);
1445 if (It.More()) aShape = It.Value();
1447 if (aShape.ShapeType() == TopAbs_SHELL) {
1448 BRepCheck_Shell chkShell (TopoDS::Shell(aShape));
1449 if (chkShell.Closed() == BRepCheck_NotClosed) {
1450 aRes = "WRN_SHAPE_UNCLOSED";
1454 aRes = "WRN_SHAPE_NOT_SHELL";
1466 //=============================================================================
1470 //=============================================================================
1471 TCollection_AsciiString GEOMImpl_IMeasureOperations::WhatIs (Handle(GEOM_Object) theShape)
1475 TCollection_AsciiString Astr;
1477 if (theShape.IsNull()) return Astr;
1479 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1480 if (aRefShape.IsNull()) return Astr;
1482 TopoDS_Shape aShape = aRefShape->GetValue();
1483 if (aShape.IsNull()) {
1484 SetErrorCode("The Objects has NULL Shape");
1488 //Compute the parameters
1489 if (aShape.ShapeType() == TopAbs_EDGE) {
1490 if (BRep_Tool::Degenerated(TopoDS::Edge(aShape))) {
1491 Astr = Astr + " It is a degenerated edge \n";
1495 Astr = Astr + " Number of sub-shapes : \n";
1498 #if OCC_VERSION_LARGE > 0x06010000
1501 int iType, nbTypes [TopAbs_SHAPE];
1502 for (iType = 0; iType < TopAbs_SHAPE; ++iType)
1504 nbTypes[aShape.ShapeType()]++;
1506 TopTools_MapOfShape aMapOfShape;
1507 aMapOfShape.Add(aShape);
1508 TopTools_ListOfShape aListOfShape;
1509 aListOfShape.Append(aShape);
1511 TopTools_ListIteratorOfListOfShape itL (aListOfShape);
1512 for (; itL.More(); itL.Next()) {
1513 TopoDS_Iterator it (itL.Value());
1514 for (; it.More(); it.Next()) {
1515 TopoDS_Shape s = it.Value();
1516 if (aMapOfShape.Add(s)) {
1517 aListOfShape.Append(s);
1518 nbTypes[s.ShapeType()]++;
1523 Astr = Astr + " VERTEX : " + TCollection_AsciiString(nbTypes[TopAbs_VERTEX]) + "\n";
1524 Astr = Astr + " EDGE : " + TCollection_AsciiString(nbTypes[TopAbs_EDGE]) + "\n";
1525 Astr = Astr + " WIRE : " + TCollection_AsciiString(nbTypes[TopAbs_WIRE]) + "\n";
1526 Astr = Astr + " FACE : " + TCollection_AsciiString(nbTypes[TopAbs_FACE]) + "\n";
1527 Astr = Astr + " SHELL : " + TCollection_AsciiString(nbTypes[TopAbs_SHELL]) + "\n";
1528 Astr = Astr + " SOLID : " + TCollection_AsciiString(nbTypes[TopAbs_SOLID]) + "\n";
1529 Astr = Astr + " COMPSOLID : " + TCollection_AsciiString(nbTypes[TopAbs_COMPSOLID]) + "\n";
1530 Astr = Astr + " COMPOUND : " + TCollection_AsciiString(nbTypes[TopAbs_COMPOUND]) + "\n";
1531 Astr = Astr + " SHAPE : " + TCollection_AsciiString(aMapOfShape.Extent());
1533 catch (Standard_Failure) {
1534 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1535 SetErrorCode(aFail->GetMessageString());
1544 //=======================================================================
1545 //function : CheckSingularCase
1546 //purpose : auxilary for GetMinDistance()
1547 // workaround for bugs 19899, 19908 and 19910 from Mantis
1548 //=======================================================================
1549 static double CheckSingularCase(const TopoDS_Shape& aSh1,
1550 const TopoDS_Shape& aSh2,
1551 gp_Pnt& Ptmp1, gp_Pnt& Ptmp2)
1553 bool IsChange1 = false;
1554 double AddDist1 = 0.0;
1555 TopExp_Explorer anExp;
1556 TopoDS_Shape tmpSh1, tmpSh2;
1558 for ( anExp.Init( aSh1, TopAbs_FACE ); anExp.More(); anExp.Next() ) {
1560 tmpSh1 = anExp.Current();
1563 TopoDS_Shape sh = aSh1;
1564 while(sh.ShapeType()==TopAbs_COMPOUND) {
1565 TopoDS_Iterator it(sh);
1568 Handle(Geom_Surface) S = BRep_Tool::Surface(TopoDS::Face(tmpSh1));
1569 if( S->IsKind(STANDARD_TYPE(Geom_SphericalSurface)) ||
1570 S->IsKind(STANDARD_TYPE(Geom_ToroidalSurface)) ) {
1571 if( sh.ShapeType()==TopAbs_SHELL || sh.ShapeType()==TopAbs_FACE ) {
1574 // changes for 0020677: EDF 1219 GEOM: MinDistance gives 0 instead of 20.88
1575 //S->Bounds(U1,U2,V1,V2); changed by
1576 ShapeAnalysis::GetFaceUVBounds(TopoDS::Face(tmpSh1),U1,U2,V1,V2);
1577 // end of changes for 020677 (dmv)
1578 Handle(Geom_RectangularTrimmedSurface) TrS1 =
1579 new Geom_RectangularTrimmedSurface(S,U1,(U1+U2)/2.,V1,V2);
1580 Handle(Geom_RectangularTrimmedSurface) TrS2 =
1581 new Geom_RectangularTrimmedSurface(S,(U1+U2)/2.,U2,V1,V2);
1584 TopoDS_Compound Comp;
1585 B.MakeCompound(Comp);
1586 B.MakeFace(F1,TrS1,1.e-7);
1588 B.MakeFace(F2,TrS2,1.e-7);
1590 Handle(ShapeFix_Shape) sfs = new ShapeFix_Shape;
1592 sfs->SetPrecision(1.e-6);
1593 sfs->SetMaxTolerance(1.0);
1595 tmpSh1 = sfs->Shape();
1599 if( S->IsKind(STANDARD_TYPE(Geom_SphericalSurface)) ) {
1600 Handle(Geom_SphericalSurface) SS = Handle(Geom_SphericalSurface)::DownCast(S);
1601 gp_Pnt PC = SS->Location();
1604 B.MakeVertex(V,PC,1.e-7);
1606 AddDist1 = SS->Radius();
1610 Handle(Geom_ToroidalSurface) TS = Handle(Geom_ToroidalSurface)::DownCast(S);
1611 gp_Ax3 ax3 = TS->Position();
1612 Handle(Geom_Circle) C = new Geom_Circle(ax3.Ax2(),TS->MajorRadius());
1615 B.MakeEdge(E,C,1.e-7);
1617 AddDist1 = TS->MinorRadius();
1627 bool IsChange2 = false;
1628 double AddDist2 = 0.0;
1630 for ( anExp.Init( aSh2, TopAbs_FACE ); anExp.More(); anExp.Next() ) {
1632 tmpSh2 = anExp.Current();
1635 TopoDS_Shape sh = aSh2;
1636 while(sh.ShapeType()==TopAbs_COMPOUND) {
1637 TopoDS_Iterator it(sh);
1640 Handle(Geom_Surface) S = BRep_Tool::Surface(TopoDS::Face(tmpSh2));
1641 if( S->IsKind(STANDARD_TYPE(Geom_SphericalSurface)) ||
1642 S->IsKind(STANDARD_TYPE(Geom_ToroidalSurface)) ) {
1643 if( sh.ShapeType()==TopAbs_SHELL || sh.ShapeType()==TopAbs_FACE ) {
1646 //S->Bounds(U1,U2,V1,V2);
1647 ShapeAnalysis::GetFaceUVBounds(TopoDS::Face(tmpSh2),U1,U2,V1,V2);
1648 Handle(Geom_RectangularTrimmedSurface) TrS1 =
1649 new Geom_RectangularTrimmedSurface(S,U1,(U1+U2)/2.,V1,V2);
1650 Handle(Geom_RectangularTrimmedSurface) TrS2 =
1651 new Geom_RectangularTrimmedSurface(S,(U1+U2)/2.,U2,V1,V2);
1654 TopoDS_Compound Comp;
1655 B.MakeCompound(Comp);
1656 B.MakeFace(F1,TrS1,1.e-7);
1658 B.MakeFace(F2,TrS2,1.e-7);
1660 Handle(ShapeFix_Shape) sfs = new ShapeFix_Shape;
1662 sfs->SetPrecision(1.e-6);
1663 sfs->SetMaxTolerance(1.0);
1665 tmpSh2 = sfs->Shape();
1669 if( S->IsKind(STANDARD_TYPE(Geom_SphericalSurface)) ) {
1670 Handle(Geom_SphericalSurface) SS = Handle(Geom_SphericalSurface)::DownCast(S);
1671 gp_Pnt PC = SS->Location();
1674 B.MakeVertex(V,PC,1.e-7);
1676 AddDist2 = SS->Radius();
1679 else if( S->IsKind(STANDARD_TYPE(Geom_ToroidalSurface)) ) {
1680 Handle(Geom_ToroidalSurface) TS = Handle(Geom_ToroidalSurface)::DownCast(S);
1681 gp_Ax3 ax3 = TS->Position();
1682 Handle(Geom_Circle) C = new Geom_Circle(ax3.Ax2(),TS->MajorRadius());
1685 B.MakeEdge(E,C,1.e-7);
1687 AddDist2 = TS->MinorRadius();
1698 if( !IsChange1 && !IsChange2 )
1701 BRepExtrema_DistShapeShape dst(tmpSh1,tmpSh2);
1703 double MinDist = 1.e9;
1704 gp_Pnt PMin1, PMin2, P1, P2;
1705 for (int i = 1; i <= dst.NbSolution(); i++) {
1706 P1 = dst.PointOnShape1(i);
1707 P2 = dst.PointOnShape2(i);
1708 Standard_Real Dist = P1.Distance(P2);
1709 if (MinDist > Dist) {
1720 gp_Dir aDir(gp_Vec(PMin1,PMin2));
1721 if( MinDist > (AddDist1+AddDist2) ) {
1722 Ptmp1 = gp_Pnt( PMin1.X() + aDir.X()*AddDist1,
1723 PMin1.Y() + aDir.Y()*AddDist1,
1724 PMin1.Z() + aDir.Z()*AddDist1 );
1725 Ptmp2 = gp_Pnt( PMin2.X() - aDir.X()*AddDist2,
1726 PMin2.Y() - aDir.Y()*AddDist2,
1727 PMin2.Z() - aDir.Z()*AddDist2 );
1728 return (MinDist - AddDist1 - AddDist2);
1731 if( AddDist1 > 0 ) {
1732 Ptmp1 = gp_Pnt( PMin1.X() + aDir.X()*AddDist1,
1733 PMin1.Y() + aDir.Y()*AddDist1,
1734 PMin1.Z() + aDir.Z()*AddDist1 );
1738 Ptmp2 = gp_Pnt( PMin2.X() - aDir.X()*AddDist2,
1739 PMin2.Y() - aDir.Y()*AddDist2,
1740 PMin2.Z() - aDir.Z()*AddDist2 );
1745 double res = MinDist - AddDist1 - AddDist2;
1746 if(res<0.) res = 0.0;
1752 static bool CheckSingularCase(const TopoDS_Shape& aSh1,
1753 const TopoDS_Shape& aSh2,
1756 TopExp_Explorer anExp;
1757 TopoDS_Shape tmpSh1, tmpSh2;
1759 for ( anExp.Init( aSh1, TopAbs_FACE ); anExp.More(); anExp.Next() ) {
1761 tmpSh1 = anExp.Current();
1764 Handle(Geom_Surface) S1 = BRep_Tool::Surface(TopoDS::Face(tmpSh1));
1765 if( S1->IsKind(STANDARD_TYPE(Geom_SphericalSurface)) ||
1766 S1->IsKind(STANDARD_TYPE(Geom_ToroidalSurface)) ) {
1768 for ( anExp.Init( aSh2, TopAbs_FACE ); anExp.More(); anExp.Next() ) {
1770 tmpSh2 = anExp.Current();
1771 Handle(Geom_Surface) S2 = BRep_Tool::Surface(TopoDS::Face(tmpSh2));
1772 GeomAPI_IntSS ISS(S1,S2,1.e-7);
1774 for(int i=1; i<=ISS.NbLines(); i++) {
1775 Handle(Geom_Curve) C3d = ISS.Line(i);
1778 B.MakeEdge(E,C3d,1.e-7);
1779 BRepExtrema_DistShapeShape dst(tmpSh2,E);
1781 gp_Pnt PMin1, PMin2, P1, P2;
1782 double MinDist = 1.e9;
1783 for (int i = 1; i <= dst.NbSolution(); i++) {
1784 P1 = dst.PointOnShape1(i);
1785 P2 = dst.PointOnShape2(i);
1786 Standard_Real Dist = P1.Distance(P2);
1787 if (MinDist > Dist) {
1801 for ( anExp.Init( aSh2, TopAbs_FACE ); anExp.More(); anExp.Next() ) {
1803 tmpSh1 = anExp.Current();
1806 Handle(Geom_Surface) S1 = BRep_Tool::Surface(TopoDS::Face(tmpSh1));
1807 if( S1->IsKind(STANDARD_TYPE(Geom_SphericalSurface)) ||
1808 S1->IsKind(STANDARD_TYPE(Geom_ToroidalSurface)) ) {
1810 for ( anExp.Init( aSh1, TopAbs_FACE ); anExp.More(); anExp.Next() ) {
1812 tmpSh2 = anExp.Current();
1813 Handle(Geom_Surface) S2 = BRep_Tool::Surface(TopoDS::Face(tmpSh2));
1814 GeomAPI_IntSS ISS(S1,S2,1.e-7);
1816 for(int i=1; i<=ISS.NbLines(); i++) {
1817 Handle(Geom_Curve) C3d = ISS.Line(i);
1820 B.MakeEdge(E,C3d,1.e-7);
1821 BRepExtrema_DistShapeShape dst(tmpSh2,E);
1824 double MinDist = 1.e9;
1825 for (int i = 1; i <= dst.NbSolution(); i++) {
1826 P1 = dst.PointOnShape1(i);
1827 P2 = dst.PointOnShape2(i);
1828 Standard_Real Dist = P1.Distance(P2);
1829 if (MinDist > Dist) {
1847 //=============================================================================
1851 //=============================================================================
1852 std::vector<bool> GEOMImpl_IMeasureOperations::AreCoordsInside(Handle(GEOM_Object) theShape,
1853 const std::vector<double>& coords,
1856 std::vector<bool> res;
1857 if (!theShape.IsNull()) {
1858 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1859 if (!aRefShape.IsNull()) {
1860 TopoDS_Shape aShape = aRefShape->GetValue();
1861 if (!aShape.IsNull()) {
1862 BRepClass3d_SolidClassifier SC(aShape);
1863 unsigned int nb_points = coords.size()/3;
1864 for (int i = 0; i < nb_points; i++) {
1865 double x = coords[3*i];
1866 double y = coords[3*i+1];
1867 double z = coords[3*i+2];
1868 gp_Pnt aPnt(x, y, z);
1869 SC.Perform(aPnt, tolerance);
1870 res.push_back( ( SC.State() == TopAbs_IN ) || ( SC.State() == TopAbs_ON ) );
1878 //=============================================================================
1882 //=============================================================================
1883 Standard_Real GEOMImpl_IMeasureOperations::GetMinDistance
1884 (Handle(GEOM_Object) theShape1, Handle(GEOM_Object) theShape2,
1885 Standard_Real& X1, Standard_Real& Y1, Standard_Real& Z1,
1886 Standard_Real& X2, Standard_Real& Y2, Standard_Real& Z2)
1889 Standard_Real MinDist = 1.e9;
1891 if (theShape1.IsNull() || theShape2.IsNull()) return MinDist;
1893 Handle(GEOM_Function) aRefShape1 = theShape1->GetLastFunction();
1894 Handle(GEOM_Function) aRefShape2 = theShape2->GetLastFunction();
1895 if (aRefShape1.IsNull() || aRefShape2.IsNull()) return MinDist;
1897 TopoDS_Shape aShape1 = aRefShape1->GetValue();
1898 TopoDS_Shape aShape2 = aRefShape2->GetValue();
1899 if (aShape1.IsNull() || aShape2.IsNull()) {
1900 SetErrorCode("One of Objects has NULL Shape");
1904 //Compute the parameters
1906 #if OCC_VERSION_LARGE > 0x06010000
1910 // Issue 0020231: A min distance bug with torus and vertex.
1911 // Make GetMinDistance() return zero if a sole VERTEX is inside any of SOLIDs
1913 // which of shapes consists of only one vertex?
1914 TopExp_Explorer exp1(aShape1,TopAbs_VERTEX), exp2(aShape2,TopAbs_VERTEX);
1915 TopoDS_Shape V1 = exp1.More() ? exp1.Current() : TopoDS_Shape();
1916 TopoDS_Shape V2 = exp2.More() ? exp2.Current() : TopoDS_Shape();
1917 exp1.Next(); exp2.Next();
1918 if ( exp1.More() ) V1.Nullify();
1919 if ( exp2.More() ) V2.Nullify();
1920 // vertex and container of solids
1921 TopoDS_Shape V = V1.IsNull() ? V2 : V1;
1922 TopoDS_Shape S = V1.IsNull() ? aShape1 : aShape2;
1923 if ( !V.IsNull() ) {
1924 // classify vertex against solids
1925 gp_Pnt p = BRep_Tool::Pnt( TopoDS::Vertex( V ) );
1926 for ( exp1.Init( S, TopAbs_SOLID ); exp1.More(); exp1.Next() ) {
1927 BRepClass3d_SolidClassifier classifier( exp1.Current(), p, 1e-6);
1928 if ( classifier.State() == TopAbs_IN ) {
1929 p.Coord(X1, Y1, Z1);
1930 p.Coord(X2, Y2, Z2);
1936 // End Issue 0020231
1939 // additional workaround for bugs 19899, 19908 and 19910 from Mantis
1940 gp_Pnt Ptmp1, Ptmp2;
1941 double dist = CheckSingularCase(aShape1, aShape2, Ptmp1, Ptmp2);
1943 Ptmp1.Coord(X1, Y1, Z1);
1944 Ptmp2.Coord(X2, Y2, Z2);
1949 BRepExtrema_DistShapeShape dst (aShape1, aShape2);
1951 gp_Pnt PMin1, PMin2, P1, P2;
1953 for (int i = 1; i <= dst.NbSolution(); i++) {
1954 P1 = dst.PointOnShape1(i);
1955 P2 = dst.PointOnShape2(i);
1957 Standard_Real Dist = P1.Distance(P2);
1958 if (MinDist > Dist) {
1965 PMin1.Coord(X1, Y1, Z1);
1966 PMin2.Coord(X2, Y2, Z2);
1969 catch (Standard_Failure) {
1970 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1971 SetErrorCode(aFail->GetMessageString());
1979 //=======================================================================
1981 * Get coordinates of point
1983 //=======================================================================
1984 void GEOMImpl_IMeasureOperations::PointCoordinates (Handle(GEOM_Object) theShape,
1985 Standard_Real& theX, Standard_Real& theY, Standard_Real& theZ)
1989 if (theShape.IsNull())
1992 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1993 if (aRefShape.IsNull())
1996 TopoDS_Shape aShape = aRefShape->GetValue();
1997 if (aShape.IsNull() || aShape.ShapeType() != TopAbs_VERTEX)
1999 SetErrorCode( "Shape must be a vertex" );
2004 #if OCC_VERSION_LARGE > 0x06010000
2007 gp_Pnt aPnt = BRep_Tool::Pnt( TopoDS::Vertex( aShape ) );
2014 catch (Standard_Failure)
2016 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
2017 SetErrorCode( aFail->GetMessageString() );
2021 //=======================================================================
2023 * Compute angle (in degrees) between two lines
2025 //=======================================================================
2026 Standard_Real GEOMImpl_IMeasureOperations::GetAngle (Handle(GEOM_Object) theLine1,
2027 Handle(GEOM_Object) theLine2)
2029 if (theLine1->GetType() == GEOM_VECTOR &&
2030 theLine2->GetType() == GEOM_VECTOR)
2031 return GetAngleBtwVectors(theLine1, theLine2);
2035 Standard_Real anAngle = -1.0;
2037 if (theLine1.IsNull() || theLine2.IsNull())
2040 Handle(GEOM_Function) aRefLine1 = theLine1->GetLastFunction();
2041 Handle(GEOM_Function) aRefLine2 = theLine2->GetLastFunction();
2042 if (aRefLine1.IsNull() || aRefLine2.IsNull())
2045 TopoDS_Shape aLine1 = aRefLine1->GetValue();
2046 TopoDS_Shape aLine2 = aRefLine2->GetValue();
2047 if (aLine1.IsNull() || aLine2.IsNull() ||
2048 aLine1.ShapeType() != TopAbs_EDGE ||
2049 aLine2.ShapeType() != TopAbs_EDGE)
2051 SetErrorCode("Two edges must be given");
2056 #if OCC_VERSION_LARGE > 0x06010000
2059 TopoDS_Edge E1 = TopoDS::Edge(aLine1);
2060 TopoDS_Edge E2 = TopoDS::Edge(aLine2);
2063 Handle(Geom_Curve) C1 = BRep_Tool::Curve(E1,fp,lp);
2064 Handle(Geom_Curve) C2 = BRep_Tool::Curve(E2,fp,lp);
2066 if ( C1.IsNull() || C2.IsNull() ||
2067 !C1->IsKind(STANDARD_TYPE(Geom_Line)) ||
2068 !C2->IsKind(STANDARD_TYPE(Geom_Line)))
2070 SetErrorCode("The edges must be linear");
2074 Handle(Geom_Line) L1 = Handle(Geom_Line)::DownCast(C1);
2075 Handle(Geom_Line) L2 = Handle(Geom_Line)::DownCast(C2);
2077 gp_Lin aLin1 = L1->Lin();
2078 gp_Lin aLin2 = L2->Lin();
2080 anAngle = aLin1.Angle(aLin2);
2081 anAngle *= 180. / M_PI; // convert radians into degrees
2083 if (anAngle > 90.0) {
2084 anAngle = 180.0 - anAngle;
2089 catch (Standard_Failure)
2091 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
2092 SetErrorCode(aFail->GetMessageString());
2098 //=======================================================================
2100 * Compute angle (in degrees) between two vectors
2102 //=======================================================================
2103 Standard_Real GEOMImpl_IMeasureOperations::GetAngleBtwVectors (Handle(GEOM_Object) theVec1,
2104 Handle(GEOM_Object) theVec2)
2108 Standard_Real anAngle = -1.0;
2110 if (theVec1.IsNull() || theVec2.IsNull())
2113 Handle(GEOM_Function) aRefVec1 = theVec1->GetLastFunction();
2114 Handle(GEOM_Function) aRefVec2 = theVec2->GetLastFunction();
2115 if (aRefVec1.IsNull() || aRefVec2.IsNull())
2118 TopoDS_Shape aVec1 = aRefVec1->GetValue();
2119 TopoDS_Shape aVec2 = aRefVec2->GetValue();
2120 if (aVec1.IsNull() || aVec2.IsNull() ||
2121 aVec1.ShapeType() != TopAbs_EDGE ||
2122 aVec2.ShapeType() != TopAbs_EDGE)
2124 SetErrorCode("Two edges must be given");
2129 #if OCC_VERSION_LARGE > 0x06010000
2132 TopoDS_Edge aE1 = TopoDS::Edge(aVec1);
2133 TopoDS_Edge aE2 = TopoDS::Edge(aVec2);
2135 TopoDS_Vertex aP11, aP12, aP21, aP22;
2136 TopExp::Vertices(aE1, aP11, aP12, Standard_True);
2137 TopExp::Vertices(aE2, aP21, aP22, Standard_True);
2138 if (aP11.IsNull() || aP12.IsNull() || aP21.IsNull() || aP22.IsNull()) {
2139 SetErrorCode("Bad edge given");
2143 gp_Vec aV1 (BRep_Tool::Pnt(aP11), BRep_Tool::Pnt(aP12));
2144 gp_Vec aV2 (BRep_Tool::Pnt(aP21), BRep_Tool::Pnt(aP22)) ;
2146 anAngle = aV1.Angle(aV2);
2147 anAngle *= 180. / M_PI; // convert radians into degrees
2151 catch (Standard_Failure)
2153 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
2154 SetErrorCode(aFail->GetMessageString());
2161 //=============================================================================
2163 * CurveCurvatureByParam
2165 //=============================================================================
2166 Standard_Real GEOMImpl_IMeasureOperations::CurveCurvatureByParam
2167 (Handle(GEOM_Object) theCurve, Standard_Real& theParam)
2170 Standard_Real aRes = -1.0;
2172 if(theCurve.IsNull()) return aRes;
2174 Handle(GEOM_Function) aRefShape = theCurve->GetLastFunction();
2175 if(aRefShape.IsNull()) return aRes;
2177 TopoDS_Shape aShape = aRefShape->GetValue();
2178 if(aShape.IsNull()) {
2179 SetErrorCode("One of Objects has NULL Shape");
2183 Standard_Real aFP, aLP, aP;
2184 Handle(Geom_Curve) aCurve = BRep_Tool::Curve(TopoDS::Edge(aShape), aFP, aLP);
2185 aP = aFP + (aLP - aFP) * theParam;
2187 if(aCurve.IsNull()) return aRes;
2191 #if OCC_VERSION_LARGE > 0x06010000
2194 GeomLProp_CLProps Prop = GeomLProp_CLProps
2195 (aCurve, aP, 2, Precision::Confusion());
2196 aRes = fabs(Prop.Curvature());
2199 catch (Standard_Failure) {
2200 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
2201 SetErrorCode(aFail->GetMessageString());
2205 if( aRes > Precision::Confusion() )
2214 //=============================================================================
2216 * CurveCurvatureByPoint
2218 //=============================================================================
2219 Standard_Real GEOMImpl_IMeasureOperations::CurveCurvatureByPoint
2220 (Handle(GEOM_Object) theCurve, Handle(GEOM_Object) thePoint)
2223 Standard_Real aRes = -1.0;
2225 if( theCurve.IsNull() || thePoint.IsNull() ) return aRes;
2227 Handle(GEOM_Function) aRefCurve = theCurve->GetLastFunction();
2228 Handle(GEOM_Function) aRefPoint = thePoint->GetLastFunction();
2229 if( aRefCurve.IsNull() || aRefPoint.IsNull() ) return aRes;
2231 TopoDS_Edge anEdge = TopoDS::Edge(aRefCurve->GetValue());
2232 TopoDS_Vertex aPnt = TopoDS::Vertex(aRefPoint->GetValue());
2233 if( anEdge.IsNull() || aPnt.IsNull() ) {
2234 SetErrorCode("One of Objects has NULL Shape");
2238 Standard_Real aFP, aLP;
2239 Handle(Geom_Curve) aCurve = BRep_Tool::Curve(anEdge, aFP, aLP);
2240 if(aCurve.IsNull()) return aRes;
2241 gp_Pnt aPoint = BRep_Tool::Pnt(aPnt);
2245 #if OCC_VERSION_LARGE > 0x06010000
2248 GeomAPI_ProjectPointOnCurve PPCurve(aPoint, aCurve, aFP, aLP);
2249 if(PPCurve.NbPoints()>0) {
2250 GeomLProp_CLProps Prop = GeomLProp_CLProps
2251 (aCurve, PPCurve.LowerDistanceParameter(), 2, Precision::Confusion());
2252 aRes = fabs(Prop.Curvature());
2256 catch (Standard_Failure) {
2257 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
2258 SetErrorCode(aFail->GetMessageString());
2262 if( aRes > Precision::Confusion() )
2271 //=============================================================================
2273 * getSurfaceCurvatures
2275 //=============================================================================
2276 Standard_Real GEOMImpl_IMeasureOperations::getSurfaceCurvatures
2277 (const Handle(Geom_Surface)& aSurf,
2278 Standard_Real theUParam,
2279 Standard_Real theVParam,
2280 Standard_Boolean theNeedMaxCurv)
2283 Standard_Real aRes = 1.0;
2285 if (aSurf.IsNull()) return aRes;
2288 #if OCC_VERSION_LARGE > 0x06010000
2291 GeomLProp_SLProps Prop = GeomLProp_SLProps
2292 (aSurf, theUParam, theVParam, 2, Precision::Confusion());
2293 if(Prop.IsCurvatureDefined()) {
2294 if(Prop.IsUmbilic()) {
2295 //cout<<"is umbilic"<<endl;
2296 aRes = fabs(Prop.MeanCurvature());
2299 //cout<<"is not umbilic"<<endl;
2300 double c1 = fabs(Prop.MaxCurvature());
2301 double c2 = fabs(Prop.MinCurvature());
2310 catch (Standard_Failure) {
2311 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
2312 SetErrorCode(aFail->GetMessageString());
2316 if( fabs(aRes) > Precision::Confusion() )
2325 //=============================================================================
2327 * MaxSurfaceCurvatureByParam
2329 //=============================================================================
2330 Standard_Real GEOMImpl_IMeasureOperations::MaxSurfaceCurvatureByParam
2331 (Handle(GEOM_Object) theSurf,
2332 Standard_Real& theUParam,
2333 Standard_Real& theVParam)
2336 Standard_Real aRes = -1.0;
2338 if (theSurf.IsNull()) return aRes;
2340 Handle(GEOM_Function) aRefShape = theSurf->GetLastFunction();
2341 if(aRefShape.IsNull()) return aRes;
2343 TopoDS_Shape aShape = aRefShape->GetValue();
2344 if(aShape.IsNull()) {
2345 SetErrorCode("One of Objects has NULL Shape");
2349 TopoDS_Face F = TopoDS::Face(aShape);
2350 Handle(Geom_Surface) aSurf = BRep_Tool::Surface(F);
2352 //Compute the parameters
2353 Standard_Real U1,U2,V1,V2;
2354 ShapeAnalysis::GetFaceUVBounds(F,U1,U2,V1,V2);
2355 Standard_Real U = U1 + (U2-U1)*theUParam;
2356 Standard_Real V = V1 + (V2-V1)*theVParam;
2358 return getSurfaceCurvatures(aSurf, U, V, true);
2362 //=============================================================================
2364 * MaxSurfaceCurvatureByPoint
2366 //=============================================================================
2367 Standard_Real GEOMImpl_IMeasureOperations::MaxSurfaceCurvatureByPoint
2368 (Handle(GEOM_Object) theSurf, Handle(GEOM_Object) thePoint)
2371 Standard_Real aRes = -1.0;
2373 if( theSurf.IsNull() || thePoint.IsNull() ) return aRes;
2375 Handle(GEOM_Function) aRefShape = theSurf->GetLastFunction();
2376 Handle(GEOM_Function) aRefPoint = thePoint->GetLastFunction();
2377 if( aRefShape.IsNull() || aRefPoint.IsNull() ) return aRes;
2379 TopoDS_Face aFace = TopoDS::Face(aRefShape->GetValue());
2380 TopoDS_Vertex aPnt = TopoDS::Vertex(aRefPoint->GetValue());
2381 if( aFace.IsNull() || aPnt.IsNull() ) {
2382 SetErrorCode("One of Objects has NULL Shape");
2386 Handle(Geom_Surface) aSurf = BRep_Tool::Surface(aFace);
2387 if(aSurf.IsNull()) return aRes;
2388 gp_Pnt aPoint = BRep_Tool::Pnt(aPnt);
2390 //Compute the parameters
2391 ShapeAnalysis_Surface sas(aSurf);
2392 gp_Pnt2d UV = sas.ValueOfUV(aPoint,Precision::Confusion());
2394 return getSurfaceCurvatures(aSurf, UV.X(), UV.Y(), true);
2398 //=============================================================================
2400 * MinSurfaceCurvatureByParam
2402 //=============================================================================
2403 Standard_Real GEOMImpl_IMeasureOperations::MinSurfaceCurvatureByParam
2404 (Handle(GEOM_Object) theSurf,
2405 Standard_Real& theUParam,
2406 Standard_Real& theVParam)
2409 Standard_Real aRes = -1.0;
2411 if (theSurf.IsNull()) return aRes;
2413 Handle(GEOM_Function) aRefShape = theSurf->GetLastFunction();
2414 if(aRefShape.IsNull()) return aRes;
2416 TopoDS_Shape aShape = aRefShape->GetValue();
2417 if(aShape.IsNull()) {
2418 SetErrorCode("One of Objects has NULL Shape");
2422 TopoDS_Face F = TopoDS::Face(aShape);
2423 Handle(Geom_Surface) aSurf = BRep_Tool::Surface(F);
2425 //Compute the parameters
2426 Standard_Real U1,U2,V1,V2;
2427 ShapeAnalysis::GetFaceUVBounds(F,U1,U2,V1,V2);
2428 Standard_Real U = U1 + (U2-U1)*theUParam;
2429 Standard_Real V = V1 + (V2-V1)*theVParam;
2431 return getSurfaceCurvatures(aSurf, U, V, false);
2435 //=============================================================================
2437 * MinSurfaceCurvatureByPoint
2439 //=============================================================================
2440 Standard_Real GEOMImpl_IMeasureOperations::MinSurfaceCurvatureByPoint
2441 (Handle(GEOM_Object) theSurf, Handle(GEOM_Object) thePoint)
2444 Standard_Real aRes = -1.0;
2446 if( theSurf.IsNull() || thePoint.IsNull() ) return aRes;
2448 Handle(GEOM_Function) aRefShape = theSurf->GetLastFunction();
2449 Handle(GEOM_Function) aRefPoint = thePoint->GetLastFunction();
2450 if( aRefShape.IsNull() || aRefPoint.IsNull() ) return aRes;
2452 TopoDS_Face aFace = TopoDS::Face(aRefShape->GetValue());
2453 TopoDS_Vertex aPnt = TopoDS::Vertex(aRefPoint->GetValue());
2454 if( aFace.IsNull() || aPnt.IsNull() ) {
2455 SetErrorCode("One of Objects has NULL Shape");
2459 Handle(Geom_Surface) aSurf = BRep_Tool::Surface(aFace);
2460 if(aSurf.IsNull()) return aRes;
2461 gp_Pnt aPoint = BRep_Tool::Pnt(aPnt);
2463 //Compute the parameters
2464 ShapeAnalysis_Surface sas(aSurf);
2465 gp_Pnt2d UV = sas.ValueOfUV(aPoint,Precision::Confusion());
2467 return getSurfaceCurvatures(aSurf, UV.X(), UV.Y(), false);
2471 //=======================================================================
2472 //function : StructuralDump
2473 //purpose : Structural (data exchange) style of output.
2474 //=======================================================================
2475 void GEOMImpl_IMeasureOperations::StructuralDump (const BRepCheck_Analyzer& theAna,
2476 const TopoDS_Shape& theShape,
2477 TCollection_AsciiString& theDump)
2481 theDump += " -- The Shape has problems :\n";
2482 theDump += " Check Count\n";
2483 theDump += " ------------------------------------------------\n";
2485 Standard_Integer last_stat = (Standard_Integer)BRepCheck_CheckFail;
2486 Handle(TColStd_HArray1OfInteger) NbProblems =
2487 new TColStd_HArray1OfInteger(1, last_stat);
2488 for (i = 1; i <= last_stat; i++)
2489 NbProblems->SetValue(i,0);
2491 Handle(TopTools_HSequenceOfShape) sl;
2492 sl = new TopTools_HSequenceOfShape();
2493 TopTools_DataMapOfShapeListOfShape theMap;
2495 GetProblemShapes(theAna, theShape, sl, NbProblems, theMap);
2498 Standard_Integer count = 0;
2499 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidPointOnCurve);
2501 theDump += " Invalid Point on Curve ................... ";
2502 theDump += TCollection_AsciiString(count) + "\n";
2504 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidPointOnCurveOnSurface);
2506 theDump += " Invalid Point on CurveOnSurface .......... ";
2507 theDump += TCollection_AsciiString(count) + "\n";
2509 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidPointOnSurface);
2511 theDump += " Invalid Point on Surface ................. ";
2512 theDump += TCollection_AsciiString(count) + "\n";
2514 count = NbProblems->Value((Standard_Integer)BRepCheck_No3DCurve);
2516 theDump += " No 3D Curve .............................. ";
2517 theDump += TCollection_AsciiString(count) + "\n";
2519 count = NbProblems->Value((Standard_Integer)BRepCheck_Multiple3DCurve);
2521 theDump += " Multiple 3D Curve ........................ ";
2522 theDump += TCollection_AsciiString(count) + "\n";
2524 count = NbProblems->Value((Standard_Integer)BRepCheck_Invalid3DCurve);
2526 theDump += " Invalid 3D Curve ......................... ";
2527 theDump += TCollection_AsciiString(count) + "\n";
2529 count = NbProblems->Value((Standard_Integer)BRepCheck_NoCurveOnSurface);
2531 theDump += " No Curve on Surface ...................... ";
2532 theDump += TCollection_AsciiString(count) + "\n";
2534 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidCurveOnSurface);
2536 theDump += " Invalid Curve on Surface ................. ";
2537 theDump += TCollection_AsciiString(count) + "\n";
2539 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidCurveOnClosedSurface);
2541 theDump += " Invalid Curve on closed Surface .......... ";
2542 theDump += TCollection_AsciiString(count) + "\n";
2544 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidSameRangeFlag);
2546 theDump += " Invalid SameRange Flag ................... ";
2547 theDump += TCollection_AsciiString(count) + "\n";
2549 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidSameParameterFlag);
2551 theDump += " Invalid SameParameter Flag ............... ";
2552 theDump += TCollection_AsciiString(count) + "\n";
2554 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidDegeneratedFlag);
2556 theDump += " Invalid Degenerated Flag ................. ";
2557 theDump += TCollection_AsciiString(count) + "\n";
2559 count = NbProblems->Value((Standard_Integer)BRepCheck_FreeEdge);
2561 theDump += " Free Edge ................................ ";
2562 theDump += TCollection_AsciiString(count) + "\n";
2564 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidMultiConnexity);
2566 theDump += " Invalid MultiConnexity ................... ";
2567 theDump += TCollection_AsciiString(count) + "\n";
2569 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidRange);
2571 theDump += " Invalid Range ............................ ";
2572 theDump += TCollection_AsciiString(count) + "\n";
2574 count = NbProblems->Value((Standard_Integer)BRepCheck_EmptyWire);
2576 theDump += " Empty Wire ............................... ";
2577 theDump += TCollection_AsciiString(count) + "\n";
2579 count = NbProblems->Value((Standard_Integer)BRepCheck_RedundantEdge);
2581 theDump += " Redundant Edge ........................... ";
2582 theDump += TCollection_AsciiString(count) + "\n";
2584 count = NbProblems->Value((Standard_Integer)BRepCheck_SelfIntersectingWire);
2586 theDump += " Self Intersecting Wire ................... ";
2587 theDump += TCollection_AsciiString(count) + "\n";
2589 count = NbProblems->Value((Standard_Integer)BRepCheck_NoSurface);
2591 theDump += " No Surface ............................... ";
2592 theDump += TCollection_AsciiString(count) + "\n";
2594 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidWire);
2596 theDump += " Invalid Wire ............................. ";
2597 theDump += TCollection_AsciiString(count) + "\n";
2599 count = NbProblems->Value((Standard_Integer)BRepCheck_RedundantWire);
2601 theDump += " Redundant Wire ........................... ";
2602 theDump += TCollection_AsciiString(count) + "\n";
2604 count = NbProblems->Value((Standard_Integer)BRepCheck_IntersectingWires);
2606 theDump += " Intersecting Wires ....................... ";
2607 theDump += TCollection_AsciiString(count) + "\n";
2609 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidImbricationOfWires);
2611 theDump += " Invalid Imbrication of Wires ............. ";
2612 theDump += TCollection_AsciiString(count) + "\n";
2614 count = NbProblems->Value((Standard_Integer)BRepCheck_EmptyShell);
2616 theDump += " Empty Shell .............................. ";
2617 theDump += TCollection_AsciiString(count) + "\n";
2619 count = NbProblems->Value((Standard_Integer)BRepCheck_RedundantFace);
2621 theDump += " Redundant Face ........................... ";
2622 theDump += TCollection_AsciiString(count) + "\n";
2624 count = NbProblems->Value((Standard_Integer)BRepCheck_UnorientableShape);
2626 theDump += " Unorientable Shape ....................... ";
2627 theDump += TCollection_AsciiString(count) + "\n";
2629 count = NbProblems->Value((Standard_Integer)BRepCheck_NotClosed);
2631 theDump += " Not Closed ............................... ";
2632 theDump += TCollection_AsciiString(count) + "\n";
2634 count = NbProblems->Value((Standard_Integer)BRepCheck_NotConnected);
2636 theDump += " Not Connected ............................ ";
2637 theDump += TCollection_AsciiString(count) + "\n";
2639 count = NbProblems->Value((Standard_Integer)BRepCheck_SubshapeNotInShape);
2641 theDump += " Sub-shape not in Shape .................... ";
2642 theDump += TCollection_AsciiString(count) + "\n";
2644 count = NbProblems->Value((Standard_Integer)BRepCheck_BadOrientation);
2646 theDump += " Bad Orientation .......................... ";
2647 theDump += TCollection_AsciiString(count) + "\n";
2649 count = NbProblems->Value((Standard_Integer)BRepCheck_BadOrientationOfSubshape);
2651 theDump += " Bad Orientation of Sub-shape .............. ";
2652 theDump += TCollection_AsciiString(count) + "\n";
2654 count = NbProblems->Value((Standard_Integer)BRepCheck_CheckFail);
2656 theDump += " checkshape failure ....................... ";
2657 theDump += TCollection_AsciiString(count) + "\n";
2660 theDump += " ------------------------------------------------\n";
2661 theDump += "*** Shapes with problems : ";
2662 theDump += TCollection_AsciiString(sl->Length()) + "\n";
2664 Standard_Integer nbv, nbe, nbw, nbf, nbs, nbo;
2665 nbv = nbe = nbw = nbf = nbs = nbo = 0;
2667 for (i = 1; i <= sl->Length(); i++) {
2668 TopoDS_Shape shi = sl->Value(i);
2669 TopAbs_ShapeEnum sti = shi.ShapeType();
2671 case TopAbs_VERTEX : nbv++; break;
2672 case TopAbs_EDGE : nbe++; break;
2673 case TopAbs_WIRE : nbw++; break;
2674 case TopAbs_FACE : nbf++; break;
2675 case TopAbs_SHELL : nbs++; break;
2676 case TopAbs_SOLID : nbo++; break;
2682 theDump += "VERTEX : ";
2683 if (nbv < 10) theDump += " ";
2684 theDump += TCollection_AsciiString(nbv) + "\n";
2687 theDump += "EDGE : ";
2688 if (nbe < 10) theDump += " ";
2689 theDump += TCollection_AsciiString(nbe) + "\n";
2692 theDump += "WIRE : ";
2693 if (nbw < 10) theDump += " ";
2694 theDump += TCollection_AsciiString(nbw) + "\n";
2697 theDump += "FACE : ";
2698 if (nbf < 10) theDump += " ";
2699 theDump += TCollection_AsciiString(nbf) + "\n";
2702 theDump += "SHELL : ";
2703 if (nbs < 10) theDump += " ";
2704 theDump += TCollection_AsciiString(nbs) + "\n";
2707 theDump += "SOLID : ";
2708 if (nbo < 10) theDump += " ";
2709 theDump += TCollection_AsciiString(nbo) + "\n";
2714 //=======================================================================
2715 //function : GetProblemShapes
2716 // purpose : for StructuralDump
2717 //=======================================================================
2718 void GEOMImpl_IMeasureOperations::GetProblemShapes (const BRepCheck_Analyzer& theAna,
2719 const TopoDS_Shape& theShape,
2720 Handle(TopTools_HSequenceOfShape)& sl,
2721 Handle(TColStd_HArray1OfInteger)& NbProblems,
2722 TopTools_DataMapOfShapeListOfShape& theMap)
2724 for (TopoDS_Iterator iter(theShape); iter.More(); iter.Next()) {
2725 GetProblemShapes(theAna, iter.Value(), sl, NbProblems, theMap);
2727 TopAbs_ShapeEnum styp = theShape.ShapeType();
2728 BRepCheck_ListIteratorOfListOfStatus itl;
2729 TopTools_ListOfShape empty;
2730 if (!theMap.IsBound(theShape)) {
2731 theMap.Bind(theShape,empty);
2733 if (!theAna.Result(theShape).IsNull()) {
2734 itl.Initialize(theAna.Result(theShape)->Status());
2735 // !!! May be, we have to print all the problems, not only the first one ?
2736 if (itl.Value() != BRepCheck_NoError) {
2737 sl->Append(theShape);
2738 BRepCheck_Status stat = itl.Value();
2739 NbProblems->SetValue((Standard_Integer)stat,
2740 NbProblems->Value((Standard_Integer)stat) + 1);
2747 GetProblemSub(theAna, theShape, sl, NbProblems, TopAbs_VERTEX, theMap);
2750 GetProblemSub(theAna, theShape, sl, NbProblems, TopAbs_WIRE, theMap);
2751 GetProblemSub(theAna, theShape, sl, NbProblems, TopAbs_EDGE, theMap);
2752 GetProblemSub(theAna, theShape, sl, NbProblems, TopAbs_VERTEX, theMap);
2757 GetProblemSub(theAna, theShape, sl, NbProblems, TopAbs_SHELL, theMap);
2764 //=======================================================================
2765 //function : Contains
2766 //=======================================================================
2767 static Standard_Boolean Contains (const TopTools_ListOfShape& L,
2768 const TopoDS_Shape& S)
2770 TopTools_ListIteratorOfListOfShape it;
2771 for (it.Initialize(L); it.More(); it.Next()) {
2772 if (it.Value().IsSame(S)) {
2773 return Standard_True;
2776 return Standard_False;
2779 //=======================================================================
2780 //function : GetProblemSub
2781 // purpose : for StructuralDump
2782 //=======================================================================
2783 void GEOMImpl_IMeasureOperations::GetProblemSub (const BRepCheck_Analyzer& theAna,
2784 const TopoDS_Shape& theShape,
2785 Handle(TopTools_HSequenceOfShape)& sl,
2786 Handle(TColStd_HArray1OfInteger)& NbProblems,
2787 const TopAbs_ShapeEnum Subtype,
2788 TopTools_DataMapOfShapeListOfShape& theMap)
2790 BRepCheck_ListIteratorOfListOfStatus itl;
2791 TopExp_Explorer exp;
2792 for (exp.Init(theShape, Subtype); exp.More(); exp.Next()) {
2793 const TopoDS_Shape& sub = exp.Current();
2795 const Handle(BRepCheck_Result)& res = theAna.Result(sub);
2796 for (res->InitContextIterator();
2797 res->MoreShapeInContext();
2798 res->NextShapeInContext()) {
2799 if (res->ContextualShape().IsSame(theShape) && !Contains(theMap(sub), theShape)) {
2800 theMap(sub).Append(theShape);
2801 itl.Initialize(res->StatusOnShape());
2803 if (itl.Value() != BRepCheck_NoError) {
2804 Standard_Integer ii = 0;
2806 for (ii = 1; ii <= sl->Length(); ii++)
2807 if (sl->Value(ii).IsSame(sub)) break;
2809 if (ii > sl->Length()) {
2811 BRepCheck_Status stat = itl.Value();
2812 NbProblems->SetValue((Standard_Integer)stat,
2813 NbProblems->Value((Standard_Integer)stat) + 1);
2815 for (ii = 1; ii <= sl->Length(); ii++)
2816 if (sl->Value(ii).IsSame(theShape)) break;
2817 if (ii > sl->Length()) {
2818 sl->Append(theShape);
2819 BRepCheck_Status stat = itl.Value();
2820 NbProblems->SetValue((Standard_Integer)stat,
2821 NbProblems->Value((Standard_Integer)stat) + 1);