1 // Copyright (C) 2007-2016 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 <Basics_OCCTVersion.hxx>
25 #include <GEOMImpl_IMeasureOperations.hxx>
26 #include <GEOMImpl_IMeasure.hxx>
27 #include <GEOMImpl_MeasureDriver.hxx>
28 #include <GEOMImpl_Types.hxx>
30 #include <GEOMUtils.hxx>
32 #include <GEOMAlgo_AlgoTools.hxx>
33 #include <GEOMAlgo_KindOfName.hxx>
34 #include <GEOMAlgo_ShapeInfoFiller.hxx>
36 #include <GEOM_PythonDump.hxx>
38 #include <utilities.h>
41 #include <Bnd_Box.hxx>
42 #include <BOPAlgo_CheckerSI.hxx>
43 #include <BOPCol_ListOfShape.hxx>
44 #include <BOPDS_DS.hxx>
45 #include <BOPDS_MapOfPassKey.hxx>
46 #include <BRepBndLib.hxx>
47 #include <BRepBuilderAPI_Copy.hxx>
48 #include <BRepCheck_ListIteratorOfListOfStatus.hxx>
49 #include <BRepCheck_Shell.hxx>
50 #include <BRepClass3d_SolidClassifier.hxx>
51 #include <BRepClass_FaceClassifier.hxx>
52 #include <BRepExtrema_DistShapeShape.hxx>
53 #include <BRepExtrema_ShapeProximity.hxx>
54 #if OCC_VERSION_LARGE > 0x06090000
55 #include <BRepExtrema_SelfIntersection.hxx>
56 #include <BRepExtrema_MapOfIntegerPackedMapOfInteger.hxx>
58 #include <BRepGProp.hxx>
59 #include <BRepTools.hxx>
60 #include <BRep_Tool.hxx>
61 #include <Geom_Line.hxx>
62 #include <GeomAPI_ProjectPointOnCurve.hxx>
63 #include <GeomAPI_ProjectPointOnSurf.hxx>
64 #include <GeomLProp_CLProps.hxx>
65 #include <GeomLProp_SLProps.hxx>
66 #include <Geom_Plane.hxx>
67 #include <GProp_GProps.hxx>
68 #include <GProp_PrincipalProps.hxx>
69 #include <ShapeAnalysis.hxx>
70 #include <ShapeAnalysis_Surface.hxx>
72 #include <TopExp_Explorer.hxx>
74 #include <TopoDS_Edge.hxx>
75 #include <TopTools_IndexedMapOfShape.hxx>
76 #include <TopTools_DataMapIteratorOfDataMapOfIntegerListOfShape.hxx>
77 #include <TColStd_MapIteratorOfPackedMapOfInteger.hxx>
78 #include <TopTools_ListIteratorOfListOfShape.hxx>
79 #include <TopTools_ListOfShape.hxx>
80 #include <Standard_ErrorHandler.hxx> // CAREFUL ! position of this file is critic : see Lucien PIGNOLONI / OCC
84 //=============================================================================
88 //=============================================================================
89 GEOMImpl_IMeasureOperations::GEOMImpl_IMeasureOperations (GEOM_Engine* theEngine, int theDocID)
90 : GEOM_IOperations(theEngine, theDocID)
92 MESSAGE("GEOMImpl_IMeasureOperations::GEOMImpl_IMeasureOperations");
95 //=============================================================================
99 //=============================================================================
100 GEOMImpl_IMeasureOperations::~GEOMImpl_IMeasureOperations()
102 MESSAGE("GEOMImpl_IMeasureOperations::~GEOMImpl_IMeasureOperations");
105 //=============================================================================
106 /*! Get kind and parameters of the given shape.
108 //=============================================================================
109 GEOMImpl_IMeasureOperations::ShapeKind GEOMImpl_IMeasureOperations::KindOfShape
110 (Handle(GEOM_Object) theShape,
111 Handle(TColStd_HSequenceOfInteger)& theIntegers,
112 Handle(TColStd_HSequenceOfReal)& theDoubles)
115 ShapeKind aKind = SK_NO_SHAPE;
117 if (theIntegers.IsNull()) theIntegers = new TColStd_HSequenceOfInteger;
118 else theIntegers->Clear();
120 if (theDoubles.IsNull()) theDoubles = new TColStd_HSequenceOfReal;
121 else theDoubles->Clear();
123 if (theShape.IsNull())
126 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
127 if (aRefShape.IsNull()) return aKind;
129 TopoDS_Shape aShape = aRefShape->GetValue();
130 if (aShape.IsNull()) return aKind;
132 int geom_type = theShape->GetType();
134 // check if it's advanced shape
135 if ( geom_type > USER_TYPE ) {
141 GEOMAlgo_ShapeInfoFiller aSF;
142 aSF.SetShape(aShape);
144 Standard_Integer iErr = aSF.ErrorStatus();
146 SetErrorCode("Error in GEOMAlgo_ShapeInfoFiller");
149 const GEOMAlgo_ShapeInfo& anInfo = aSF.Info();
151 // specific processing for some "advandced" objects
152 switch ( geom_type ) {
154 // local coordinate systen
155 // (+) geompy.kind.LCS xc yc zc xx xy xz yx yy yz zx zy zz
157 TopoDS_Face aFace = TopoDS::Face( aShape );
158 Handle(Geom_Plane) aPlane = Handle(Geom_Plane)::DownCast( BRep_Tool::Surface( aFace ) );
159 gp_Pnt aC = aPlane->Pln().Location();
160 gp_Ax3 anAx3 = aPlane->Pln().Position();
162 theDoubles->Append(aC.X());
163 theDoubles->Append(aC.Y());
164 theDoubles->Append(aC.Z());
166 gp_Dir aD = anAx3.XDirection();
167 theDoubles->Append(aD.X());
168 theDoubles->Append(aD.Y());
169 theDoubles->Append(aD.Z());
170 aD = anAx3.YDirection();
171 theDoubles->Append(aD.X());
172 theDoubles->Append(aD.Y());
173 theDoubles->Append(aD.Z());
174 aD = anAx3.Direction();
175 theDoubles->Append(aD.X());
176 theDoubles->Append(aD.Y());
177 theDoubles->Append(aD.Z());
183 // Interprete results
184 TopAbs_ShapeEnum aType = anInfo.Type();
187 case TopAbs_COMPOUND:
188 case TopAbs_COMPSOLID:
190 // (+) geompy.kind.COMPOUND nb_solids nb_faces nb_edges nb_vertices
191 // (+) geompy.kind.COMPSOLID nb_solids nb_faces nb_edges nb_vertices
192 // ??? "nb_faces" - all faces or only 'standalone' faces?
193 if (aType == TopAbs_COMPOUND)
196 aKind = SK_COMPSOLID;
198 //theIntegers->Append(anInfo.NbSubShapes(TopAbs_COMPOUND));
199 //theIntegers->Append(anInfo.NbSubShapes(TopAbs_COMPSOLID));
200 theIntegers->Append(anInfo.NbSubShapes(TopAbs_SOLID));
201 theIntegers->Append(anInfo.NbSubShapes(TopAbs_FACE));
202 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
203 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
209 // (+) geompy.kind.SHELL geompy.info.closed nb_faces nb_edges nb_vertices
210 // (+) geompy.kind.SHELL geompy.info.unclosed nb_faces nb_edges nb_vertices
213 theIntegers->Append((int)anInfo.KindOfClosed());
215 theIntegers->Append(anInfo.NbSubShapes(TopAbs_FACE));
216 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
217 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
223 // (+) geompy.kind.WIRE geompy.info.closed nb_edges nb_vertices
224 // (+) geompy.kind.WIRE geompy.info.unclosed nb_edges nb_vertices
227 theIntegers->Append((int)anInfo.KindOfClosed());
229 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
230 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
238 GEOMAlgo_KindOfName aKN = anInfo.KindOfName();
241 case GEOMAlgo_KN_SPHERE:
242 // (+) geompy.kind.SPHERE xc yc zc R
246 gp_Pnt aC = anInfo.Location();
247 theDoubles->Append(aC.X());
248 theDoubles->Append(aC.Y());
249 theDoubles->Append(aC.Z());
251 theDoubles->Append(anInfo.Radius1());
254 case GEOMAlgo_KN_CYLINDER:
255 // (+) geompy.kind.CYLINDER xb yb zb dx dy dz R H
259 gp_Pnt aC = anInfo.Location();
260 theDoubles->Append(aC.X());
261 theDoubles->Append(aC.Y());
262 theDoubles->Append(aC.Z());
264 gp_Ax3 anAx3 = anInfo.Position();
265 gp_Dir aD = anAx3.Direction();
266 theDoubles->Append(aD.X());
267 theDoubles->Append(aD.Y());
268 theDoubles->Append(aD.Z());
270 theDoubles->Append(anInfo.Radius1());
271 theDoubles->Append(anInfo.Height());
274 case GEOMAlgo_KN_BOX:
275 // (+) geompy.kind.BOX xc yc zc ax ay az
279 gp_Pnt aC = anInfo.Location();
280 theDoubles->Append(aC.X());
281 theDoubles->Append(aC.Y());
282 theDoubles->Append(aC.Z());
284 gp_Ax3 anAx3 = anInfo.Position();
285 gp_Dir aD = anAx3.Direction();
286 gp_Dir aX = anAx3.XDirection();
289 if (aD.IsParallel(gp::DZ(), Precision::Angular()) &&
290 aX.IsParallel(gp::DX(), Precision::Angular())) {
291 theDoubles->Append(anInfo.Length()); // ax'
292 theDoubles->Append(anInfo.Width()); // ay'
293 theDoubles->Append(anInfo.Height()); // az'
295 else if (aD.IsParallel(gp::DZ(), Precision::Angular()) &&
296 aX.IsParallel(gp::DY(), Precision::Angular())) {
297 theDoubles->Append(anInfo.Width()); // ay'
298 theDoubles->Append(anInfo.Length()); // ax'
299 theDoubles->Append(anInfo.Height()); // az'
301 else if (aD.IsParallel(gp::DX(), Precision::Angular()) &&
302 aX.IsParallel(gp::DZ(), Precision::Angular())) {
303 theDoubles->Append(anInfo.Height()); // az'
304 theDoubles->Append(anInfo.Width()); // ay'
305 theDoubles->Append(anInfo.Length()); // ax'
307 else if (aD.IsParallel(gp::DX(), Precision::Angular()) &&
308 aX.IsParallel(gp::DY(), Precision::Angular())) {
309 theDoubles->Append(anInfo.Height()); // az'
310 theDoubles->Append(anInfo.Length()); // ax'
311 theDoubles->Append(anInfo.Width()); // ay'
313 else if (aD.IsParallel(gp::DY(), Precision::Angular()) &&
314 aX.IsParallel(gp::DZ(), Precision::Angular())) {
315 theDoubles->Append(anInfo.Width()); // ay'
316 theDoubles->Append(anInfo.Height()); // az'
317 theDoubles->Append(anInfo.Length()); // ax'
319 else if (aD.IsParallel(gp::DY(), Precision::Angular()) &&
320 aX.IsParallel(gp::DX(), Precision::Angular())) {
321 theDoubles->Append(anInfo.Length()); // ax'
322 theDoubles->Append(anInfo.Height()); // az'
323 theDoubles->Append(anInfo.Width()); // ay'
326 // (+) geompy.kind.ROTATED_BOX xo yo zo zx zy zz xx xy xz ax ay az
327 aKind = SK_ROTATED_BOX;
329 // Direction and XDirection
330 theDoubles->Append(aD.X());
331 theDoubles->Append(aD.Y());
332 theDoubles->Append(aD.Z());
334 theDoubles->Append(aX.X());
335 theDoubles->Append(aX.Y());
336 theDoubles->Append(aX.Z());
339 theDoubles->Append(anInfo.Length());
340 theDoubles->Append(anInfo.Width());
341 theDoubles->Append(anInfo.Height());
345 case GEOMAlgo_KN_TORUS:
346 // (+) geompy.kind.TORUS xc yc zc dx dy dz R_1 R_2
350 gp_Pnt aO = anInfo.Location();
351 theDoubles->Append(aO.X());
352 theDoubles->Append(aO.Y());
353 theDoubles->Append(aO.Z());
355 gp_Ax3 anAx3 = anInfo.Position();
356 gp_Dir aD = anAx3.Direction();
357 theDoubles->Append(aD.X());
358 theDoubles->Append(aD.Y());
359 theDoubles->Append(aD.Z());
361 theDoubles->Append(anInfo.Radius1());
362 theDoubles->Append(anInfo.Radius2());
365 case GEOMAlgo_KN_CONE:
366 // (+) geompy.kind.CONE xb yb zb dx dy dz R_1 R_2 H
370 gp_Pnt aO = anInfo.Location();
371 theDoubles->Append(aO.X());
372 theDoubles->Append(aO.Y());
373 theDoubles->Append(aO.Z());
375 gp_Ax3 anAx3 = anInfo.Position();
376 gp_Dir aD = anAx3.Direction();
377 theDoubles->Append(aD.X());
378 theDoubles->Append(aD.Y());
379 theDoubles->Append(aD.Z());
381 theDoubles->Append(anInfo.Radius1());
382 theDoubles->Append(anInfo.Radius2());
383 theDoubles->Append(anInfo.Height());
386 case GEOMAlgo_KN_POLYHEDRON:
387 // (+) geompy.kind.POLYHEDRON nb_faces nb_edges nb_vertices
389 aKind = SK_POLYHEDRON;
391 theIntegers->Append(anInfo.NbSubShapes(TopAbs_FACE));
392 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
393 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
397 // (+) geompy.kind.SOLID nb_faces nb_edges nb_vertices
399 theIntegers->Append(anInfo.NbSubShapes(TopAbs_FACE));
400 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
401 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
411 GEOMAlgo_KindOfName aKN = anInfo.KindOfName();
413 case GEOMAlgo_KN_SPHERE:
414 // (+) geompy.kind.SPHERE2D xc yc zc R
418 gp_Pnt aC = anInfo.Location();
419 theDoubles->Append(aC.X());
420 theDoubles->Append(aC.Y());
421 theDoubles->Append(aC.Z());
423 theDoubles->Append(anInfo.Radius1());
426 case GEOMAlgo_KN_CYLINDER:
427 // (+) geompy.kind.CYLINDER2D xb yb zb dx dy dz R H
429 aKind = SK_CYLINDER2D;
431 gp_Pnt aO = anInfo.Location();
432 theDoubles->Append(aO.X());
433 theDoubles->Append(aO.Y());
434 theDoubles->Append(aO.Z());
436 gp_Ax3 anAx3 = anInfo.Position();
437 gp_Dir aD = anAx3.Direction();
438 theDoubles->Append(aD.X());
439 theDoubles->Append(aD.Y());
440 theDoubles->Append(aD.Z());
442 theDoubles->Append(anInfo.Radius1());
443 theDoubles->Append(anInfo.Height());
446 case GEOMAlgo_KN_TORUS:
447 // (+) geompy.kind.TORUS2D xc yc zc dx dy dz R_1 R_2
451 gp_Pnt aO = anInfo.Location();
452 theDoubles->Append(aO.X());
453 theDoubles->Append(aO.Y());
454 theDoubles->Append(aO.Z());
456 gp_Ax3 anAx3 = anInfo.Position();
457 gp_Dir aD = anAx3.Direction();
458 theDoubles->Append(aD.X());
459 theDoubles->Append(aD.Y());
460 theDoubles->Append(aD.Z());
462 theDoubles->Append(anInfo.Radius1());
463 theDoubles->Append(anInfo.Radius2());
466 case GEOMAlgo_KN_CONE:
467 // (+) geompy.kind.CONE2D xc yc zc dx dy dz R_1 R_2 H
471 gp_Pnt aO = anInfo.Location();
472 theDoubles->Append(aO.X());
473 theDoubles->Append(aO.Y());
474 theDoubles->Append(aO.Z());
476 gp_Ax3 anAx3 = anInfo.Position();
477 gp_Dir aD = anAx3.Direction();
478 theDoubles->Append(aD.X());
479 theDoubles->Append(aD.Y());
480 theDoubles->Append(aD.Z());
482 theDoubles->Append(anInfo.Radius1());
483 theDoubles->Append(anInfo.Radius2());
484 theDoubles->Append(anInfo.Height());
487 case GEOMAlgo_KN_DISKCIRCLE:
488 // (+) geompy.kind.DISK_CIRCLE xc yc zc dx dy dz R
490 aKind = SK_DISK_CIRCLE;
492 gp_Pnt aC = anInfo.Location();
493 theDoubles->Append(aC.X());
494 theDoubles->Append(aC.Y());
495 theDoubles->Append(aC.Z());
497 gp_Ax3 anAx3 = anInfo.Position();
498 gp_Dir aD = anAx3.Direction();
499 theDoubles->Append(aD.X());
500 theDoubles->Append(aD.Y());
501 theDoubles->Append(aD.Z());
503 theDoubles->Append(anInfo.Radius1());
506 case GEOMAlgo_KN_DISKELLIPSE:
507 // (+) geompy.kind.DISK_ELLIPSE xc yc zc dx dy dz R_1 R_2
509 aKind = SK_DISK_ELLIPSE;
511 gp_Pnt aC = anInfo.Location();
512 theDoubles->Append(aC.X());
513 theDoubles->Append(aC.Y());
514 theDoubles->Append(aC.Z());
516 gp_Ax3 anAx3 = anInfo.Position();
517 gp_Dir aD = anAx3.Direction();
518 theDoubles->Append(aD.X());
519 theDoubles->Append(aD.Y());
520 theDoubles->Append(aD.Z());
522 theDoubles->Append(anInfo.Radius1());
523 theDoubles->Append(anInfo.Radius2());
526 case GEOMAlgo_KN_RECTANGLE:
527 case GEOMAlgo_KN_TRIANGLE:
528 case GEOMAlgo_KN_QUADRANGLE:
529 case GEOMAlgo_KN_POLYGON:
530 // (+) geompy.kind.POLYGON xo yo zo dx dy dz nb_edges nb_vertices
534 gp_Pnt aO = anInfo.Location();
535 theDoubles->Append(aO.X());
536 theDoubles->Append(aO.Y());
537 theDoubles->Append(aO.Z());
539 gp_Ax3 anAx3 = anInfo.Position();
540 gp_Dir aD = anAx3.Direction();
541 theDoubles->Append(aD.X());
542 theDoubles->Append(aD.Y());
543 theDoubles->Append(aD.Z());
545 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
546 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
549 case GEOMAlgo_KN_PLANE: // infinite
550 // (+) geompy.kind.PLANE xo yo zo dx dy dz
554 gp_Pnt aC = anInfo.Location();
555 theDoubles->Append(aC.X());
556 theDoubles->Append(aC.Y());
557 theDoubles->Append(aC.Z());
559 gp_Ax3 anAx3 = anInfo.Position();
560 gp_Dir aD = anAx3.Direction();
561 theDoubles->Append(aD.X());
562 theDoubles->Append(aD.Y());
563 theDoubles->Append(aD.Z());
565 if (anInfo.KindOfBounds() != GEOMAlgo_KB_INFINITE)
567 // (+) geompy.kind.PLANAR xo yo zo dx dy dz nb_edges nb_vertices
570 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
571 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
576 // ??? geompy.kind.FACE nb_edges nb_vertices _surface_type_id_
577 // (+) geompy.kind.FACE nb_edges nb_vertices
578 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
579 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
588 GEOMAlgo_KindOfName aKN = anInfo.KindOfName();
590 case GEOMAlgo_KN_CIRCLE:
592 // (+) geompy.kind.CIRCLE xc yc zc dx dy dz R
595 gp_Pnt aC = anInfo.Location();
596 theDoubles->Append(aC.X());
597 theDoubles->Append(aC.Y());
598 theDoubles->Append(aC.Z());
600 gp_Ax3 anAx3 = anInfo.Position();
601 gp_Dir aD = anAx3.Direction();
602 theDoubles->Append(aD.X());
603 theDoubles->Append(aD.Y());
604 theDoubles->Append(aD.Z());
606 theDoubles->Append(anInfo.Radius1());
609 case GEOMAlgo_KN_ARCCIRCLE:
611 // (+) geompy.kind.ARC_CIRCLE xc yc zc dx dy dz R x1 y1 z1 x2 y2 z2
612 aKind = SK_ARC_CIRCLE;
614 gp_Pnt aC = anInfo.Location();
615 theDoubles->Append(aC.X());
616 theDoubles->Append(aC.Y());
617 theDoubles->Append(aC.Z());
619 gp_Ax3 anAx3 = anInfo.Position();
620 gp_Dir aD = anAx3.Direction();
621 theDoubles->Append(aD.X());
622 theDoubles->Append(aD.Y());
623 theDoubles->Append(aD.Z());
625 theDoubles->Append(anInfo.Radius1());
627 gp_Pnt aP1 = anInfo.Pnt1();
628 theDoubles->Append(aP1.X());
629 theDoubles->Append(aP1.Y());
630 theDoubles->Append(aP1.Z());
632 gp_Pnt aP2 = anInfo.Pnt2();
633 theDoubles->Append(aP2.X());
634 theDoubles->Append(aP2.Y());
635 theDoubles->Append(aP2.Z());
638 case GEOMAlgo_KN_ELLIPSE:
640 // (+) geompy.kind.ELLIPSE xc yc zc dx dy dz R_1 R_2
643 gp_Pnt aC = anInfo.Location();
644 theDoubles->Append(aC.X());
645 theDoubles->Append(aC.Y());
646 theDoubles->Append(aC.Z());
648 gp_Ax3 anAx3 = anInfo.Position();
649 gp_Dir aD = anAx3.Direction();
650 theDoubles->Append(aD.X());
651 theDoubles->Append(aD.Y());
652 theDoubles->Append(aD.Z());
654 theDoubles->Append(anInfo.Radius1());
655 theDoubles->Append(anInfo.Radius2());
658 case GEOMAlgo_KN_ARCELLIPSE:
660 // (+) geompy.kind.ARC_ELLIPSE xc yc zc dx dy dz R_1 R_2 x1 y1 z1 x2 y2 z2
661 aKind = SK_ARC_ELLIPSE;
663 gp_Pnt aC = anInfo.Location();
664 theDoubles->Append(aC.X());
665 theDoubles->Append(aC.Y());
666 theDoubles->Append(aC.Z());
668 gp_Ax3 anAx3 = anInfo.Position();
669 gp_Dir aD = anAx3.Direction();
670 theDoubles->Append(aD.X());
671 theDoubles->Append(aD.Y());
672 theDoubles->Append(aD.Z());
674 theDoubles->Append(anInfo.Radius1());
675 theDoubles->Append(anInfo.Radius2());
677 gp_Pnt aP1 = anInfo.Pnt1();
678 theDoubles->Append(aP1.X());
679 theDoubles->Append(aP1.Y());
680 theDoubles->Append(aP1.Z());
682 gp_Pnt aP2 = anInfo.Pnt2();
683 theDoubles->Append(aP2.X());
684 theDoubles->Append(aP2.Y());
685 theDoubles->Append(aP2.Z());
688 case GEOMAlgo_KN_LINE:
690 // ??? geompy.kind.LINE x1 y1 z1 x2 y2 z2
691 // (+) geompy.kind.LINE x1 y1 z1 dx dy dz
694 gp_Pnt aO = anInfo.Location();
695 theDoubles->Append(aO.X());
696 theDoubles->Append(aO.Y());
697 theDoubles->Append(aO.Z());
699 gp_Dir aD = anInfo.Direction();
700 theDoubles->Append(aD.X());
701 theDoubles->Append(aD.Y());
702 theDoubles->Append(aD.Z());
705 case GEOMAlgo_KN_SEGMENT:
707 // (+) geompy.kind.SEGMENT x1 y1 z1 x2 y2 z2
710 gp_Pnt aP1 = anInfo.Pnt1();
711 theDoubles->Append(aP1.X());
712 theDoubles->Append(aP1.Y());
713 theDoubles->Append(aP1.Z());
715 gp_Pnt aP2 = anInfo.Pnt2();
716 theDoubles->Append(aP2.X());
717 theDoubles->Append(aP2.Y());
718 theDoubles->Append(aP2.Z());
722 // ??? geompy.kind.EDGE nb_vertices _curve_type_id_
723 // (+) geompy.kind.EDGE nb_vertices
724 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
731 // (+) geompy.kind.VERTEX x y z
734 gp_Pnt aP = anInfo.Location();
735 theDoubles->Append(aP.X());
736 theDoubles->Append(aP.Y());
737 theDoubles->Append(aP.Z());
746 //=============================================================================
750 //=============================================================================
751 void GEOMImpl_IMeasureOperations::GetPosition
752 (Handle(GEOM_Object) theShape,
753 Standard_Real& Ox, Standard_Real& Oy, Standard_Real& Oz,
754 Standard_Real& Zx, Standard_Real& Zy, Standard_Real& Zz,
755 Standard_Real& Xx, Standard_Real& Xy, Standard_Real& Xz)
759 //Set default values: global CS
760 Ox = Oy = Oz = Zx = Zy = Xy = Xz = 0.;
763 if (theShape.IsNull()) return;
765 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
766 if (aRefShape.IsNull()) return;
768 TopoDS_Shape aShape = aRefShape->GetValue();
769 if (aShape.IsNull()) {
770 SetErrorCode("The Objects has NULL Shape");
777 gp_Ax3 anAx3 = GEOMUtils::GetPosition(aShape);
779 gp_Pnt anOri = anAx3.Location();
780 gp_Dir aDirZ = anAx3.Direction();
781 gp_Dir aDirX = anAx3.XDirection();
784 anOri.Coord(Ox, Oy, Oz);
785 aDirZ.Coord(Zx, Zy, Zz);
786 aDirX.Coord(Xx, Xy, Xz);
788 catch (Standard_Failure) {
789 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
790 SetErrorCode(aFail->GetMessageString());
797 //=============================================================================
801 //=============================================================================
802 Handle(GEOM_Object) GEOMImpl_IMeasureOperations::GetCentreOfMass
803 (Handle(GEOM_Object) theShape)
807 if (theShape.IsNull()) return NULL;
809 //Add a new CentreOfMass object
810 Handle(GEOM_Object) aCDG = GetEngine()->AddObject(GetDocID(), GEOM_CDG);
812 //Add a new CentreOfMass function
813 Handle(GEOM_Function) aFunction =
814 aCDG->AddFunction(GEOMImpl_MeasureDriver::GetID(), CDG_MEASURE);
815 if (aFunction.IsNull()) return NULL;
817 //Check if the function is set correctly
818 if (aFunction->GetDriverGUID() != GEOMImpl_MeasureDriver::GetID()) return NULL;
820 GEOMImpl_IMeasure aCI (aFunction);
822 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
823 if (aRefShape.IsNull()) return NULL;
825 aCI.SetBase(aRefShape);
827 //Compute the CentreOfMass value
830 if (!GetSolver()->ComputeFunction(aFunction)) {
831 SetErrorCode("Measure driver failed to compute centre of mass");
835 catch (Standard_Failure) {
836 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
837 SetErrorCode(aFail->GetMessageString());
841 //Make a Python command
842 GEOM::TPythonDump(aFunction) << aCDG << " = geompy.MakeCDG(" << theShape << ")";
848 //=============================================================================
852 //=============================================================================
853 Handle(GEOM_Object) GEOMImpl_IMeasureOperations::GetVertexByIndex
854 (Handle(GEOM_Object) theShape,
855 Standard_Integer theIndex)
859 if (theShape.IsNull()) return NULL;
861 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
862 if (aRefShape.IsNull()) return NULL;
864 //Add a new Vertex object
865 Handle(GEOM_Object) aVertex = GetEngine()->AddObject(GetDocID(), GEOM_POINT);
868 Handle(GEOM_Function) aFunction =
869 aVertex->AddFunction(GEOMImpl_MeasureDriver::GetID(), VERTEX_BY_INDEX);
870 if (aFunction.IsNull()) return NULL;
872 //Check if the function is set correctly
873 if (aFunction->GetDriverGUID() != GEOMImpl_MeasureDriver::GetID()) return NULL;
875 GEOMImpl_IMeasure aCI (aFunction);
876 aCI.SetBase(aRefShape);
877 aCI.SetIndex(theIndex);
882 if (!GetSolver()->ComputeFunction(aFunction)) {
883 SetErrorCode("Vertex by index driver failed.");
887 catch (Standard_Failure) {
888 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
889 SetErrorCode(aFail->GetMessageString());
893 //Make a Python command
894 GEOM::TPythonDump(aFunction) << aVertex << " = geompy.GetVertexByIndex(" << theShape << ", " << theIndex << ")";
900 //=============================================================================
904 //=============================================================================
905 Handle(GEOM_Object) GEOMImpl_IMeasureOperations::GetNormal
906 (Handle(GEOM_Object) theFace,
907 Handle(GEOM_Object) theOptionalPoint)
911 if (theFace.IsNull()) return NULL;
913 //Add a new Normale object
914 Handle(GEOM_Object) aNorm = GetEngine()->AddObject(GetDocID(), GEOM_VECTOR);
916 //Add a new Normale function
917 Handle(GEOM_Function) aFunction =
918 aNorm->AddFunction(GEOMImpl_MeasureDriver::GetID(), VECTOR_FACE_NORMALE);
919 if (aFunction.IsNull()) return NULL;
921 //Check if the function is set correctly
922 if (aFunction->GetDriverGUID() != GEOMImpl_MeasureDriver::GetID()) return NULL;
924 GEOMImpl_IMeasure aCI (aFunction);
926 Handle(GEOM_Function) aFace = theFace->GetLastFunction();
927 if (aFace.IsNull()) return NULL;
931 if (!theOptionalPoint.IsNull()) {
932 Handle(GEOM_Function) anOptPnt = theOptionalPoint->GetLastFunction();
933 aCI.SetPoint(anOptPnt);
936 //Compute the Normale value
939 if (!GetSolver()->ComputeFunction(aFunction)) {
940 SetErrorCode("Measure driver failed to compute normake of face");
944 catch (Standard_Failure) {
945 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
946 SetErrorCode(aFail->GetMessageString());
950 //Make a Python command
951 GEOM::TPythonDump pd (aFunction);
952 pd << aNorm << " = geompy.GetNormal(" << theFace;
953 if (!theOptionalPoint.IsNull()) {
954 pd << ", " << theOptionalPoint;
962 //=============================================================================
966 //=============================================================================
967 void GEOMImpl_IMeasureOperations::GetBasicProperties (Handle(GEOM_Object) theShape,
968 Standard_Real& theLength,
969 Standard_Real& theSurfArea,
970 Standard_Real& theVolume)
974 if (theShape.IsNull()) return;
976 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
977 if (aRefShape.IsNull()) return;
979 TopoDS_Shape aShape = aRefShape->GetValue();
980 if (aShape.IsNull()) {
981 SetErrorCode("The Objects has NULL Shape");
985 //Compute the parameters
986 GProp_GProps LProps, SProps;
987 Standard_Real anEps = 1.e-6;
990 BRepGProp::LinearProperties(aShape, LProps);
991 theLength = LProps.Mass();
993 BRepGProp::SurfaceProperties(aShape, SProps, anEps);
994 theSurfArea = SProps.Mass();
997 if (aShape.ShapeType() < TopAbs_SHELL) {
998 for (TopExp_Explorer Exp (aShape, TopAbs_SOLID); Exp.More(); Exp.Next()) {
1000 BRepGProp::VolumeProperties(Exp.Current(), VProps, anEps);
1001 theVolume += VProps.Mass();
1005 catch (Standard_Failure) {
1006 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1007 SetErrorCode(aFail->GetMessageString());
1014 //=============================================================================
1018 //=============================================================================
1019 void GEOMImpl_IMeasureOperations::GetInertia
1020 (Handle(GEOM_Object) theShape,
1021 Standard_Real& I11, Standard_Real& I12, Standard_Real& I13,
1022 Standard_Real& I21, Standard_Real& I22, Standard_Real& I23,
1023 Standard_Real& I31, Standard_Real& I32, Standard_Real& I33,
1024 Standard_Real& Ix , Standard_Real& Iy , Standard_Real& Iz)
1028 if (theShape.IsNull()) return;
1030 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1031 if (aRefShape.IsNull()) return;
1033 TopoDS_Shape aShape = aRefShape->GetValue();
1034 if (aShape.IsNull()) {
1035 SetErrorCode("The Objects has NULL Shape");
1039 //Compute the parameters
1040 GProp_GProps System;
1044 if (aShape.ShapeType() == TopAbs_VERTEX ||
1045 aShape.ShapeType() == TopAbs_EDGE ||
1046 aShape.ShapeType() == TopAbs_WIRE) {
1047 BRepGProp::LinearProperties(aShape, System);
1048 } else if (aShape.ShapeType() == TopAbs_FACE ||
1049 aShape.ShapeType() == TopAbs_SHELL) {
1050 BRepGProp::SurfaceProperties(aShape, System);
1052 BRepGProp::VolumeProperties(aShape, System);
1054 gp_Mat I = System.MatrixOfInertia();
1068 GProp_PrincipalProps Pr = System.PrincipalProperties();
1069 Pr.Moments(Ix,Iy,Iz);
1071 catch (Standard_Failure) {
1072 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1073 SetErrorCode(aFail->GetMessageString());
1080 //=============================================================================
1084 //=============================================================================
1085 void GEOMImpl_IMeasureOperations::GetBoundingBox
1086 (Handle(GEOM_Object) theShape,
1087 const Standard_Boolean precise,
1088 Standard_Real& Xmin, Standard_Real& Xmax,
1089 Standard_Real& Ymin, Standard_Real& Ymax,
1090 Standard_Real& Zmin, Standard_Real& Zmax)
1094 if (theShape.IsNull()) return;
1096 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1097 if (aRefShape.IsNull()) return;
1099 TopoDS_Shape aShape = aRefShape->GetValue();
1100 if (aShape.IsNull()) {
1101 SetErrorCode("The Objects has NULL Shape");
1105 //Compute the parameters
1110 BRepBuilderAPI_Copy aCopyTool (aShape);
1111 if (!aCopyTool.IsDone()) {
1112 SetErrorCode("GetBoundingBox Error: Bad shape detected");
1116 aShape = aCopyTool.Shape();
1118 // remove triangulation to obtain more exact boundaries
1119 BRepTools::Clean(aShape);
1121 BRepBndLib::Add(aShape, B);
1124 if (!GEOMUtils::PreciseBoundingBox(aShape, B)) {
1125 SetErrorCode("GetBoundingBox Error: Bounding box cannot be precised");
1130 B.Get(Xmin, Ymin, Zmin, Xmax, Ymax, Zmax);
1132 catch (Standard_Failure) {
1133 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1134 SetErrorCode(aFail->GetMessageString());
1141 //=============================================================================
1145 //=============================================================================
1146 Handle(GEOM_Object) GEOMImpl_IMeasureOperations::GetBoundingBox
1147 (Handle(GEOM_Object) theShape,
1148 const Standard_Boolean precise)
1152 if (theShape.IsNull()) return NULL;
1154 //Add a new BoundingBox object
1155 Handle(GEOM_Object) aBnd = GetEngine()->AddObject(GetDocID(), GEOM_BOX);
1157 //Add a new BoundingBox function
1158 const int aType = (precise ? BND_BOX_MEASURE_PRECISE : BND_BOX_MEASURE);
1159 Handle(GEOM_Function) aFunction =
1160 aBnd->AddFunction(GEOMImpl_MeasureDriver::GetID(), aType);
1161 if (aFunction.IsNull()) return NULL;
1163 //Check if the function is set correctly
1164 if (aFunction->GetDriverGUID() != GEOMImpl_MeasureDriver::GetID()) return NULL;
1166 GEOMImpl_IMeasure aCI (aFunction);
1168 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1169 if (aRefShape.IsNull()) return NULL;
1171 aCI.SetBase(aRefShape);
1173 //Compute the BoundingBox value
1176 if (!GetSolver()->ComputeFunction(aFunction)) {
1177 SetErrorCode("Measure driver failed to compute a bounding box");
1181 catch (Standard_Failure) {
1182 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1183 SetErrorCode(aFail->GetMessageString());
1187 //Make a Python command
1188 GEOM::TPythonDump aPd(aFunction);
1190 aPd << aBnd << " = geompy.MakeBoundingBox(" << theShape;
1202 //=============================================================================
1206 //=============================================================================
1207 void GEOMImpl_IMeasureOperations::GetTolerance
1208 (Handle(GEOM_Object) theShape,
1209 Standard_Real& FaceMin, Standard_Real& FaceMax,
1210 Standard_Real& EdgeMin, Standard_Real& EdgeMax,
1211 Standard_Real& VertMin, Standard_Real& VertMax)
1215 if (theShape.IsNull()) return;
1217 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1218 if (aRefShape.IsNull()) return;
1220 TopoDS_Shape aShape = aRefShape->GetValue();
1221 if (aShape.IsNull()) {
1222 SetErrorCode("The Objects has NULL Shape");
1226 //Compute the parameters
1228 FaceMin = EdgeMin = VertMin = RealLast();
1229 FaceMax = EdgeMax = VertMax = -RealLast();
1233 for (TopExp_Explorer ExF (aShape, TopAbs_FACE); ExF.More(); ExF.Next()) {
1234 TopoDS_Face Face = TopoDS::Face(ExF.Current());
1235 T = BRep_Tool::Tolerance(Face);
1241 for (TopExp_Explorer ExE (aShape, TopAbs_EDGE); ExE.More(); ExE.Next()) {
1242 TopoDS_Edge Edge = TopoDS::Edge(ExE.Current());
1243 T = BRep_Tool::Tolerance(Edge);
1249 for (TopExp_Explorer ExV (aShape, TopAbs_VERTEX); ExV.More(); ExV.Next()) {
1250 TopoDS_Vertex Vertex = TopoDS::Vertex(ExV.Current());
1251 T = BRep_Tool::Tolerance(Vertex);
1258 catch (Standard_Failure) {
1259 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1260 SetErrorCode(aFail->GetMessageString());
1267 //=============================================================================
1271 //=============================================================================
1272 bool GEOMImpl_IMeasureOperations::CheckShape (Handle(GEOM_Object) theShape,
1273 const Standard_Boolean theIsCheckGeom,
1274 std::list<ShapeError> &theErrors)
1279 if (theShape.IsNull()) return false;
1281 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1282 if (aRefShape.IsNull()) return false;
1284 TopoDS_Shape aShape = aRefShape->GetValue();
1285 if (aShape.IsNull()) {
1286 SetErrorCode("The Objects has NULL Shape");
1290 //Compute the parameters
1291 bool isValid = false;
1294 BRepCheck_Analyzer ana (aShape, theIsCheckGeom);
1295 if (ana.IsValid()) {
1298 FillErrors(ana, aShape, theErrors);
1301 catch (Standard_Failure) {
1302 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1303 SetErrorCode(aFail->GetMessageString());
1311 //=============================================================================
1315 //=============================================================================
1316 TCollection_AsciiString GEOMImpl_IMeasureOperations::PrintShapeErrors
1317 (Handle(GEOM_Object) theShape,
1318 const std::list<ShapeError> &theErrors)
1320 TCollection_AsciiString aDump;
1322 if (theShape.IsNull()) {
1326 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1328 if (aRefShape.IsNull()) {
1332 TopoDS_Shape aShape = aRefShape->GetValue();
1334 if (aShape.IsNull()) {
1335 SetErrorCode("The Objects has NULL Shape");
1339 if (!theErrors.empty()) {
1340 // The shape is not valid. Prepare errors for dump.
1341 TopTools_IndexedMapOfShape anIndices;
1342 std::list<ShapeError>::const_iterator anIter = theErrors.begin();
1343 Standard_Integer nbv, nbe, nbw, nbf, nbs, nbo;
1344 nbv = nbe = nbw = nbf = nbs = nbo = 0;
1346 // Map sub-shapes and their indices
1347 TopExp::MapShapes(aShape, anIndices);
1349 const Standard_Integer aNbSubShapes = anIndices.Extent();
1350 TColStd_MapOfInteger aMapPbInd;
1352 aDump += " -- The Shape has problems :\n";
1353 aDump += " Check Count\n";
1354 aDump += " ------------------------------------------------\n";
1356 for (; anIter != theErrors.end(); anIter++) {
1357 Standard_Integer aNbShapes = anIter->incriminated.size();
1359 switch(anIter->error) {
1360 case BRepCheck_InvalidPointOnCurve:
1361 aDump += " Invalid Point on Curve ................... ";
1363 case BRepCheck_InvalidPointOnCurveOnSurface:
1364 aDump += " Invalid Point on CurveOnSurface .......... ";
1366 case BRepCheck_InvalidPointOnSurface:
1367 aDump += " Invalid Point on Surface ................. ";
1369 case BRepCheck_No3DCurve:
1370 aDump += " No 3D Curve .............................. ";
1372 case BRepCheck_Multiple3DCurve:
1373 aDump += " Multiple 3D Curve ........................ ";
1375 case BRepCheck_Invalid3DCurve:
1376 aDump += " Invalid 3D Curve ......................... ";
1378 case BRepCheck_NoCurveOnSurface:
1379 aDump += " No Curve on Surface ...................... ";
1381 case BRepCheck_InvalidCurveOnSurface:
1382 aDump += " Invalid Curve on Surface ................. ";
1384 case BRepCheck_InvalidCurveOnClosedSurface:
1385 aDump += " Invalid Curve on closed Surface .......... ";
1387 case BRepCheck_InvalidSameRangeFlag:
1388 aDump += " Invalid SameRange Flag ................... ";
1390 case BRepCheck_InvalidSameParameterFlag:
1391 aDump += " Invalid SameParameter Flag ............... ";
1393 case BRepCheck_InvalidDegeneratedFlag:
1394 aDump += " Invalid Degenerated Flag ................. ";
1396 case BRepCheck_FreeEdge:
1397 aDump += " Free Edge ................................ ";
1399 case BRepCheck_InvalidMultiConnexity:
1400 aDump += " Invalid MultiConnexity ................... ";
1402 case BRepCheck_InvalidRange:
1403 aDump += " Invalid Range ............................ ";
1405 case BRepCheck_EmptyWire:
1406 aDump += " Empty Wire ............................... ";
1408 case BRepCheck_RedundantEdge:
1409 aDump += " Redundant Edge ........................... ";
1411 case BRepCheck_SelfIntersectingWire:
1412 aDump += " Self Intersecting Wire ................... ";
1414 case BRepCheck_NoSurface:
1415 aDump += " No Surface ............................... ";
1417 case BRepCheck_InvalidWire:
1418 aDump += " Invalid Wire ............................. ";
1420 case BRepCheck_RedundantWire:
1421 aDump += " Redundant Wire ........................... ";
1423 case BRepCheck_IntersectingWires:
1424 aDump += " Intersecting Wires ....................... ";
1426 case BRepCheck_InvalidImbricationOfWires:
1427 aDump += " Invalid Imbrication of Wires ............. ";
1429 case BRepCheck_EmptyShell:
1430 aDump += " Empty Shell .............................. ";
1432 case BRepCheck_RedundantFace:
1433 aDump += " Redundant Face ........................... ";
1435 case BRepCheck_UnorientableShape:
1436 aDump += " Unorientable Shape ....................... ";
1438 case BRepCheck_NotClosed:
1439 aDump += " Not Closed ............................... ";
1441 case BRepCheck_NotConnected:
1442 aDump += " Not Connected ............................ ";
1444 case BRepCheck_SubshapeNotInShape:
1445 aDump += " Sub-shape not in Shape ................... ";
1447 case BRepCheck_BadOrientation:
1448 aDump += " Bad Orientation .......................... ";
1450 case BRepCheck_BadOrientationOfSubshape:
1451 aDump += " Bad Orientation of Sub-shape ............. ";
1453 case BRepCheck_InvalidToleranceValue:
1454 aDump += " Invalid Tolerance Value .................. ";
1456 case BRepCheck_CheckFail:
1457 aDump += " Check Shape Failure ...................... ";
1463 aDump += TCollection_AsciiString(aNbShapes) + "\n";
1465 // Count types of shape.
1466 std::list<int>::const_iterator aShIter = anIter->incriminated.begin();
1468 for (; aShIter != anIter->incriminated.end(); aShIter++) {
1469 const Standard_Integer anIndex = *aShIter;
1471 if (anIndex > 0 && anIndex <= aNbSubShapes && aMapPbInd.Add(anIndex)) {
1472 const TopoDS_Shape &aSubShape = anIndices.FindKey(anIndex);
1473 const TopAbs_ShapeEnum aType = aSubShape.ShapeType();
1476 case TopAbs_VERTEX : nbv++; break;
1477 case TopAbs_EDGE : nbe++; break;
1478 case TopAbs_WIRE : nbw++; break;
1479 case TopAbs_FACE : nbf++; break;
1480 case TopAbs_SHELL : nbs++; break;
1481 case TopAbs_SOLID : nbo++; break;
1488 const Standard_Integer aNbFaultyShapes = nbv + nbe + nbw + nbf + nbs + nbo;
1489 aDump += " ------------------------------------------------\n";
1490 aDump += "*** Shapes with problems : ";
1491 aDump += TCollection_AsciiString(aNbFaultyShapes) + "\n";
1494 aDump += "VERTEX : ";
1495 if (nbv < 10) aDump += " ";
1496 aDump += TCollection_AsciiString(nbv) + "\n";
1500 if (nbe < 10) aDump += " ";
1501 aDump += TCollection_AsciiString(nbe) + "\n";
1505 if (nbw < 10) aDump += " ";
1506 aDump += TCollection_AsciiString(nbw) + "\n";
1510 if (nbf < 10) aDump += " ";
1511 aDump += TCollection_AsciiString(nbf) + "\n";
1514 aDump += "SHELL : ";
1515 if (nbs < 10) aDump += " ";
1516 aDump += TCollection_AsciiString(nbs) + "\n";
1519 aDump += "SOLID : ";
1520 if (nbo < 10) aDump += " ";
1521 aDump += TCollection_AsciiString(nbo) + "\n";
1528 //=============================================================================
1530 * CheckSelfIntersections
1532 //=============================================================================
1533 bool GEOMImpl_IMeasureOperations::CheckSelfIntersections
1534 (Handle(GEOM_Object) theShape,
1535 const SICheckLevel theCheckLevel,
1536 Handle(TColStd_HSequenceOfInteger)& theIntersections)
1540 if (theIntersections.IsNull())
1541 theIntersections = new TColStd_HSequenceOfInteger;
1543 theIntersections->Clear();
1545 if (theShape.IsNull())
1548 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1549 if (aRefShape.IsNull()) return false;
1551 TopoDS_Shape aShape = aRefShape->GetValue();
1552 if (aShape.IsNull()) return false;
1555 TopoDS_Shape aScopy;
1557 GEOMAlgo_AlgoTools::CopyShape(aShape, aScopy);
1559 // Map sub-shapes and their indices
1560 TopTools_IndexedMapOfShape anIndices;
1561 TopExp::MapShapes(aScopy, anIndices);
1563 BOPCol_ListOfShape aLCS;
1564 aLCS.Append(aScopy);
1566 BOPAlgo_CheckerSI aCSI; // checker of self-interferences
1567 aCSI.SetArguments(aLCS);
1568 aCSI.SetLevelOfCheck(theCheckLevel);
1570 // 1. Launch the checker
1572 Standard_Integer iErr = aCSI.ErrorStatus();
1575 Standard_Integer aNbS, n1, n2;
1576 BOPDS_MapIteratorMapOfPassKey aItMPK;
1578 // 2. Take the shapes from DS
1579 const BOPDS_DS& aDS = aCSI.DS();
1580 aNbS=aDS.NbShapes();
1582 // 3. Get the pairs of interfered shapes
1583 const BOPDS_MapOfPassKey& aMPK=aDS.Interferences();
1584 aItMPK.Initialize(aMPK);
1585 for (; aItMPK.More(); aItMPK.Next()) {
1586 const BOPDS_PassKey& aPK=aItMPK.Value();
1589 if (n1 > aNbS || n2 > aNbS){
1590 return false; // Error
1592 const TopoDS_Shape& aS1 = aDS.Shape(n1);
1593 const TopoDS_Shape& aS2 = aDS.Shape(n2);
1595 theIntersections->Append(anIndices.FindIndex(aS1));
1596 theIntersections->Append(anIndices.FindIndex(aS2));
1603 return theIntersections->IsEmpty();
1606 //=============================================================================
1608 * CheckSelfIntersectionsFast
1610 //=============================================================================
1611 bool GEOMImpl_IMeasureOperations::CheckSelfIntersectionsFast
1612 (Handle(GEOM_Object) theShape,
1613 float theDeflection, double theTolerance,
1614 Handle(TColStd_HSequenceOfInteger)& theIntersections)
1618 if (theIntersections.IsNull())
1619 theIntersections = new TColStd_HSequenceOfInteger;
1621 theIntersections->Clear();
1623 if (theShape.IsNull())
1626 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1627 if (aRefShape.IsNull()) return false;
1629 TopoDS_Shape aShape = aRefShape->GetValue();
1630 if (aShape.IsNull()) return false;
1633 TopoDS_Shape aScopy;
1635 GEOMAlgo_AlgoTools::CopyShape(aShape, aScopy);
1636 GEOMUtils::MeshShape(aScopy, theDeflection);
1638 // Map sub-shapes and their indices
1639 TopTools_IndexedMapOfShape anIndices;
1640 TopExp::MapShapes(aScopy, anIndices);
1642 #if OCC_VERSION_LARGE > 0x06090000
1643 // Checker of fast interferences
1644 BRepExtrema_SelfIntersection aTool(aScopy, (theTolerance <= 0.) ? 0.0 : theTolerance);
1646 // Launch the checker
1649 const BRepExtrema_MapOfIntegerPackedMapOfInteger& intersections = aTool.OverlapElements();
1651 std::set<Standard_Integer> processed;
1653 for (BRepExtrema_MapOfIntegerPackedMapOfInteger::Iterator it(intersections); it.More(); it.Next()) {
1654 Standard_Integer idxLeft = it.Key();
1655 if (processed.count(idxLeft) > 0) continue; // already added
1656 processed.insert(idxLeft);
1657 const TColStd_PackedMapOfInteger& overlaps = it.Value();
1658 for (TColStd_MapIteratorOfPackedMapOfInteger subit(overlaps); subit.More(); subit.Next()) {
1659 Standard_Integer idxRight = subit.Key();
1660 if (processed.count(idxRight) > 0) continue; // already added
1661 const TopoDS_Shape& aS1 = aTool.GetSubShape(idxLeft);
1662 const TopoDS_Shape& aS2 = aTool.GetSubShape(idxRight);
1663 theIntersections->Append(anIndices.FindIndex(aS1));
1664 theIntersections->Append(anIndices.FindIndex(aS2));
1672 return theIntersections->IsEmpty();
1675 //=============================================================================
1679 //=============================================================================
1680 bool GEOMImpl_IMeasureOperations::CheckBOPArguments
1681 (const Handle(GEOM_Object) &theShape)
1685 if (theShape.IsNull()) {
1689 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1691 if (aRefShape.IsNull()) {
1695 TopoDS_Shape aShape = aRefShape->GetValue();
1697 if (aShape.IsNull()) {
1701 //Compute the parameters
1702 bool isValid = GEOMUtils::CheckBOPArguments(aShape);
1709 //=============================================================================
1713 //=============================================================================
1714 bool GEOMImpl_IMeasureOperations::FastIntersect (Handle(GEOM_Object) theShape1, Handle(GEOM_Object) theShape2,
1715 double theTolerance, float theDeflection,
1716 Handle(TColStd_HSequenceOfInteger)& theIntersections1,
1717 Handle(TColStd_HSequenceOfInteger)& theIntersections2)
1720 bool isGood = false;
1722 if (theIntersections1.IsNull())
1723 theIntersections1 = new TColStd_HSequenceOfInteger;
1725 theIntersections1->Clear();
1727 if (theIntersections2.IsNull())
1728 theIntersections2 = new TColStd_HSequenceOfInteger;
1730 theIntersections2->Clear();
1732 if (theShape1.IsNull() || theShape2.IsNull()) {
1733 SetErrorCode("Objects have NULL Shape");
1737 if (theShape1 == theShape2) {
1738 SetErrorCode("Objects are equal");
1741 Handle(GEOM_Function) aRefShape1 = theShape1->GetLastFunction();
1742 Handle(GEOM_Function) aRefShape2 = theShape2->GetLastFunction();
1743 if (aRefShape1.IsNull() || aRefShape2.IsNull()) return isGood;
1745 TopoDS_Shape aShape1 = aRefShape1->GetValue();
1746 TopoDS_Shape aShape2 = aRefShape2->GetValue();
1747 if (aShape1.IsNull() || aShape2.IsNull()) return isGood;
1750 TopoDS_Shape aScopy1, aScopy2;
1751 GEOMAlgo_AlgoTools::CopyShape(aShape1, aScopy1);
1752 GEOMAlgo_AlgoTools::CopyShape(aShape2, aScopy2);
1754 GEOMUtils::MeshShape(aScopy1, theDeflection);
1755 GEOMUtils::MeshShape(aScopy2, theDeflection);
1757 // Map sub-shapes and their indices
1758 TopTools_IndexedMapOfShape anIndices1, anIndices2;
1759 TopExp::MapShapes(aScopy1, anIndices1);
1760 TopExp::MapShapes(aScopy2, anIndices2);
1762 BOPCol_ListOfShape aLCS1, aLCS2;
1763 aLCS1.Append(aScopy1); aLCS2.Append(aScopy2);
1765 BRepExtrema_ShapeProximity aBSP; // checker of fast interferences
1766 aBSP.LoadShape1(aScopy1); aBSP.LoadShape2(aScopy2);
1767 aBSP.SetTolerance((theTolerance <= 0.) ? 0.0 : theTolerance);
1769 // 1. Launch the checker
1772 // 2. Get sets of IDs of overlapped faces
1773 #if OCC_VERSION_LARGE > 0x06090000
1774 for (BRepExtrema_MapOfIntegerPackedMapOfInteger::Iterator anIt1 (aBSP.OverlapSubShapes1()); anIt1.More(); anIt1.Next())
1776 for (BRepExtrema_OverlappedSubShapes::Iterator anIt1 (aBSP.OverlapSubShapes1()); anIt1.More(); anIt1.Next())
1779 const TopoDS_Shape& aS1 = aBSP.GetSubShape1(anIt1.Key());
1780 theIntersections1->Append(anIndices1.FindIndex(aS1));
1783 #if OCC_VERSION_LARGE > 0x06090000
1784 for (BRepExtrema_MapOfIntegerPackedMapOfInteger::Iterator anIt2 (aBSP.OverlapSubShapes2()); anIt2.More(); anIt2.Next())
1786 for (BRepExtrema_OverlappedSubShapes::Iterator anIt2 (aBSP.OverlapSubShapes2()); anIt2.More(); anIt2.Next())
1789 const TopoDS_Shape& aS2 = aBSP.GetSubShape2(anIt2.Key());
1790 theIntersections2->Append(anIndices2.FindIndex(aS2));
1793 isGood = !theIntersections1->IsEmpty() && !theIntersections1->IsEmpty();
1801 //=============================================================================
1805 //=============================================================================
1806 TCollection_AsciiString GEOMImpl_IMeasureOperations::IsGoodForSolid (Handle(GEOM_Object) theShape)
1810 TCollection_AsciiString aRes = "";
1812 if (theShape.IsNull()) {
1813 aRes = "WRN_NULL_OBJECT_OR_SHAPE";
1816 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1817 if (aRefShape.IsNull()) {
1818 aRes = "WRN_NULL_OBJECT_OR_SHAPE";
1821 TopoDS_Shape aShape = aRefShape->GetValue();
1822 if (aShape.IsNull()) {
1823 aRes = "WRN_NULL_OBJECT_OR_SHAPE";
1826 if (aShape.ShapeType() == TopAbs_COMPOUND) {
1827 TopoDS_Iterator It (aShape, Standard_True, Standard_True);
1828 if (It.More()) aShape = It.Value();
1830 if (aShape.ShapeType() == TopAbs_SHELL) {
1831 BRepCheck_Shell chkShell (TopoDS::Shell(aShape));
1832 if (chkShell.Closed() == BRepCheck_NotClosed) {
1833 aRes = "WRN_SHAPE_UNCLOSED";
1837 aRes = "WRN_SHAPE_NOT_SHELL";
1849 //=============================================================================
1853 //=============================================================================
1854 TCollection_AsciiString GEOMImpl_IMeasureOperations::WhatIs (Handle(GEOM_Object) theShape)
1858 TCollection_AsciiString Astr;
1860 if (theShape.IsNull()) return Astr;
1862 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1863 if (aRefShape.IsNull()) return Astr;
1865 TopoDS_Shape aShape = aRefShape->GetValue();
1866 if (aShape.IsNull()) {
1867 SetErrorCode("The Objects has NULL Shape");
1871 //Compute the parameters
1872 if (aShape.ShapeType() == TopAbs_EDGE) {
1873 if (BRep_Tool::Degenerated(TopoDS::Edge(aShape))) {
1874 Astr = Astr + " It is a degenerated edge \n";
1878 Astr = Astr + " Number of sub-shapes : \n";
1882 int iType, nbTypes [TopAbs_SHAPE], nbFlatType [TopAbs_SHAPE];
1883 for (iType = 0; iType < TopAbs_SHAPE; ++iType) {
1885 nbFlatType[iType] = 0;
1887 nbTypes[aShape.ShapeType()]++;
1889 TopTools_MapOfShape aMapOfShape;
1890 aMapOfShape.Add(aShape);
1891 TopTools_ListOfShape aListOfShape;
1892 aListOfShape.Append(aShape);
1894 TopTools_ListIteratorOfListOfShape itL (aListOfShape);
1895 for (; itL.More(); itL.Next()) {
1896 TopoDS_Shape sp = itL.Value();
1897 TopoDS_Iterator it (sp);
1898 for (; it.More(); it.Next()) {
1899 TopoDS_Shape s = it.Value();
1900 if (aMapOfShape.Add(s)) {
1901 aListOfShape.Append(s);
1902 nbTypes[s.ShapeType()]++;
1903 if ((sp.ShapeType() == TopAbs_COMPOUND) || (sp.ShapeType() == TopAbs_COMPSOLID)) {
1904 nbFlatType[s.ShapeType()]++;
1910 Astr = Astr + " VERTEX : " + TCollection_AsciiString(nbTypes[TopAbs_VERTEX]) + "\n";
1911 Astr = Astr + " EDGE : " + TCollection_AsciiString(nbTypes[TopAbs_EDGE]) + "\n";
1912 Astr = Astr + " WIRE : " + TCollection_AsciiString(nbTypes[TopAbs_WIRE]) + "\n";
1913 Astr = Astr + " FACE : " + TCollection_AsciiString(nbTypes[TopAbs_FACE]) + "\n";
1914 Astr = Astr + " SHELL : " + TCollection_AsciiString(nbTypes[TopAbs_SHELL]) + "\n";
1915 Astr = Astr + " SOLID : " + TCollection_AsciiString(nbTypes[TopAbs_SOLID]) + "\n";
1916 Astr = Astr + " COMPSOLID : " + TCollection_AsciiString(nbTypes[TopAbs_COMPSOLID]) + "\n";
1917 Astr = Astr + " COMPOUND : " + TCollection_AsciiString(nbTypes[TopAbs_COMPOUND]) + "\n";
1918 Astr = Astr + " SHAPE : " + TCollection_AsciiString(aMapOfShape.Extent()) + "\n";
1920 if ((aShape.ShapeType() == TopAbs_COMPOUND) || (aShape.ShapeType() == TopAbs_COMPSOLID)){
1921 Astr = Astr + " --------------------- \n Flat content : \n";
1922 if (nbFlatType[TopAbs_VERTEX] > 0)
1923 Astr = Astr + " VERTEX : " + TCollection_AsciiString(nbFlatType[TopAbs_VERTEX]) + "\n";
1924 if (nbFlatType[TopAbs_EDGE] > 0)
1925 Astr = Astr + " EDGE : " + TCollection_AsciiString(nbFlatType[TopAbs_EDGE]) + "\n";
1926 if (nbFlatType[TopAbs_WIRE] > 0)
1927 Astr = Astr + " WIRE : " + TCollection_AsciiString(nbFlatType[TopAbs_WIRE]) + "\n";
1928 if (nbFlatType[TopAbs_FACE] > 0)
1929 Astr = Astr + " FACE : " + TCollection_AsciiString(nbFlatType[TopAbs_FACE]) + "\n";
1930 if (nbFlatType[TopAbs_SHELL] > 0)
1931 Astr = Astr + " SHELL : " + TCollection_AsciiString(nbFlatType[TopAbs_SHELL]) + "\n";
1932 if (nbFlatType[TopAbs_SOLID] > 0)
1933 Astr = Astr + " SOLID : " + TCollection_AsciiString(nbFlatType[TopAbs_SOLID]) + "\n";
1936 catch (Standard_Failure) {
1937 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1938 SetErrorCode(aFail->GetMessageString());
1946 //=============================================================================
1950 //=============================================================================
1952 GEOMImpl_IMeasureOperations::AreCoordsInside(Handle(GEOM_Object) theShape,
1953 const std::vector<double>& coords,
1956 std::vector<bool> isInsideRes;
1957 if (!theShape.IsNull()) {
1958 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1959 if (!aRefShape.IsNull()) {
1960 TopoDS_Shape aShape = aRefShape->GetValue();
1961 if (!aShape.IsNull())
1963 TopTools_IndexedMapOfShape mapShape;
1965 TopExp_Explorer anExp;
1966 for ( anExp.Init( aShape, TopAbs_SOLID ); anExp.More(); anExp.Next() )
1967 mapShape.Add( anExp.Current() );
1968 for ( anExp.Init( aShape, TopAbs_FACE, TopAbs_SOLID ); anExp.More(); anExp.Next() )
1969 mapShape.Add( anExp.Current() );
1970 for ( anExp.Init( aShape, TopAbs_EDGE, TopAbs_FACE ); anExp.More(); anExp.Next() )
1971 mapShape.Add( anExp.Current() );
1972 for ( anExp.Init( aShape, TopAbs_VERTEX, TopAbs_EDGE ); anExp.More(); anExp.Next() )
1973 mapShape.Add( anExp.Current() ); //// ?????????
1975 size_t nb_points = coords.size()/3, nb_points_inside = 0;
1976 isInsideRes.resize( nb_points, false );
1978 for ( int iS = 1; iS <= mapShape.Extent(); ++iS )
1980 if ( nb_points_inside == nb_points )
1982 aShape = mapShape( iS );
1983 switch ( aShape.ShapeType() ) {
1986 BRepClass3d_SolidClassifier SC( TopoDS::Solid( aShape ));
1987 for ( size_t i = 0; i < nb_points; i++)
1989 if ( isInsideRes[ i ]) continue;
1990 gp_Pnt aPnt( coords[3*i], coords[3*i+1], coords[3*i+2] );
1991 SC.Perform( aPnt, tolerance );
1992 isInsideRes[ i ] = (( SC.State() == TopAbs_IN ) || ( SC.State() == TopAbs_ON ));
1993 nb_points_inside += isInsideRes[ i ];
1999 Standard_Real u1,u2,v1,v2;
2000 const TopoDS_Face& face = TopoDS::Face( aShape );
2001 Handle(Geom_Surface) surf = BRep_Tool::Surface( face );
2002 surf->Bounds( u1,u2,v1,v2 );
2003 GeomAPI_ProjectPointOnSurf project;
2004 project.Init(surf, u1,u2, v1,v2, tolerance );
2005 for ( size_t i = 0; i < nb_points; i++)
2007 if ( isInsideRes[ i ]) continue;
2008 gp_Pnt aPnt( coords[3*i], coords[3*i+1], coords[3*i+2] );
2009 project.Perform( aPnt );
2010 if ( project.IsDone() &&
2011 project.NbPoints() > 0 &&
2012 project.LowerDistance() <= tolerance )
2014 Quantity_Parameter u, v;
2015 project.LowerDistanceParameters(u, v);
2016 gp_Pnt2d uv( u, v );
2017 BRepClass_FaceClassifier FC ( face, uv, tolerance );
2018 isInsideRes[ i ] = (( FC.State() == TopAbs_IN ) || ( FC.State() == TopAbs_ON ));
2019 nb_points_inside += isInsideRes[ i ];
2027 const TopoDS_Edge& edge = TopoDS::Edge( aShape );
2028 Handle(Geom_Curve) curve = BRep_Tool::Curve( edge, f, l );
2029 GeomAPI_ProjectPointOnCurve project;
2030 project.Init( curve, f, l );
2031 for ( size_t i = 0; i < nb_points; i++)
2033 if ( isInsideRes[ i ]) continue;
2034 gp_Pnt aPnt( coords[3*i], coords[3*i+1], coords[3*i+2] );
2035 project.Perform( aPnt );
2036 isInsideRes[ i ] = ( project.NbPoints() > 0 &&
2037 project.LowerDistance() <= tolerance );
2038 nb_points_inside += isInsideRes[ i ];
2044 gp_Pnt aVPnt = BRep_Tool::Pnt( TopoDS::Vertex( aShape ));
2045 for ( size_t i = 0; i < nb_points; i++)
2047 if ( isInsideRes[ i ]) continue;
2048 gp_Pnt aPnt( coords[3*i], coords[3*i+1], coords[3*i+2] );
2049 isInsideRes[ i ] = ( aPnt.SquareDistance( aVPnt ) <= tolerance * tolerance );
2050 nb_points_inside += isInsideRes[ i ];
2055 } // switch ( aShape.ShapeType() )
2063 //=============================================================================
2067 //=============================================================================
2069 GEOMImpl_IMeasureOperations::GetMinDistance (Handle(GEOM_Object) theShape1,
2070 Handle(GEOM_Object) theShape2,
2079 Standard_Real MinDist = 1.e9;
2081 if (theShape1.IsNull() || theShape2.IsNull()) return MinDist;
2083 Handle(GEOM_Function) aRefShape1 = theShape1->GetLastFunction();
2084 Handle(GEOM_Function) aRefShape2 = theShape2->GetLastFunction();
2085 if (aRefShape1.IsNull() || aRefShape2.IsNull()) return MinDist;
2087 TopoDS_Shape aShape1 = aRefShape1->GetValue();
2088 TopoDS_Shape aShape2 = aRefShape2->GetValue();
2089 if (aShape1.IsNull() || aShape2.IsNull()) {
2090 SetErrorCode("One of Objects has NULL Shape");
2094 //Compute the parameters
2098 gp_Pnt aPnt1, aPnt2;
2100 MinDist = GEOMUtils::GetMinDistance(aShape1, aShape2, aPnt1, aPnt2);
2102 if (MinDist >= 0.0) {
2103 aPnt1.Coord(X1, Y1, Z1);
2104 aPnt2.Coord(X2, Y2, Z2);
2109 catch (Standard_Failure) {
2110 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
2111 SetErrorCode(aFail->GetMessageString());
2119 //=======================================================================
2121 * Get coordinates of closest points of two shapes
2123 //=======================================================================
2124 Standard_Integer GEOMImpl_IMeasureOperations::ClosestPoints (Handle(GEOM_Object) theShape1,
2125 Handle(GEOM_Object) theShape2,
2126 Handle(TColStd_HSequenceOfReal)& theDoubles)
2129 Standard_Integer nbSolutions = 0;
2131 if (theShape1.IsNull() || theShape2.IsNull()) return nbSolutions;
2133 Handle(GEOM_Function) aRefShape1 = theShape1->GetLastFunction();
2134 Handle(GEOM_Function) aRefShape2 = theShape2->GetLastFunction();
2135 if (aRefShape1.IsNull() || aRefShape2.IsNull()) return nbSolutions;
2137 TopoDS_Shape aShape1 = aRefShape1->GetValue();
2138 TopoDS_Shape aShape2 = aRefShape2->GetValue();
2139 if (aShape1.IsNull() || aShape2.IsNull()) {
2140 SetErrorCode("One of Objects has NULL Shape");
2144 // Compute the extremities
2149 // additional workaround for bugs 19899, 19908 and 19910 from Mantis
2151 double dist = GEOMUtils::GetMinDistanceSingular(aShape1, aShape2, P1, P2);
2155 theDoubles->Append(P1.X());
2156 theDoubles->Append(P1.Y());
2157 theDoubles->Append(P1.Z());
2158 theDoubles->Append(P2.X());
2159 theDoubles->Append(P2.Y());
2160 theDoubles->Append(P2.Z());
2166 BRepExtrema_DistShapeShape dst (aShape1, aShape2);
2168 nbSolutions = dst.NbSolution();
2169 if (theDoubles.IsNull()) theDoubles = new TColStd_HSequenceOfReal;
2172 for (int i = 1; i <= nbSolutions; i++) {
2173 P1 = dst.PointOnShape1(i);
2174 P2 = dst.PointOnShape2(i);
2176 theDoubles->Append(P1.X());
2177 theDoubles->Append(P1.Y());
2178 theDoubles->Append(P1.Z());
2179 theDoubles->Append(P2.X());
2180 theDoubles->Append(P2.Y());
2181 theDoubles->Append(P2.Z());
2185 catch (Standard_Failure) {
2186 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
2187 SetErrorCode(aFail->GetMessageString());
2195 //=======================================================================
2197 * Get coordinates of point
2199 //=======================================================================
2200 void GEOMImpl_IMeasureOperations::PointCoordinates (Handle(GEOM_Object) theShape,
2201 Standard_Real& theX, Standard_Real& theY, Standard_Real& theZ)
2205 if (theShape.IsNull())
2208 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
2209 if (aRefShape.IsNull())
2212 TopoDS_Shape aShape = aRefShape->GetValue();
2213 if (aShape.IsNull() || aShape.ShapeType() != TopAbs_VERTEX)
2215 SetErrorCode( "Shape must be a vertex" );
2221 gp_Pnt aPnt = BRep_Tool::Pnt( TopoDS::Vertex( aShape ) );
2228 catch (Standard_Failure)
2230 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
2231 SetErrorCode( aFail->GetMessageString() );
2235 //=======================================================================
2237 * Compute angle (in degrees) between two lines
2239 //=======================================================================
2240 Standard_Real GEOMImpl_IMeasureOperations::GetAngle (Handle(GEOM_Object) theLine1,
2241 Handle(GEOM_Object) theLine2)
2243 if (theLine1->GetType() == GEOM_VECTOR &&
2244 theLine2->GetType() == GEOM_VECTOR)
2245 return GetAngleBtwVectors(theLine1, theLine2);
2249 Standard_Real anAngle = -1.0;
2251 if (theLine1.IsNull() || theLine2.IsNull())
2254 Handle(GEOM_Function) aRefLine1 = theLine1->GetLastFunction();
2255 Handle(GEOM_Function) aRefLine2 = theLine2->GetLastFunction();
2256 if (aRefLine1.IsNull() || aRefLine2.IsNull())
2259 TopoDS_Shape aLine1 = aRefLine1->GetValue();
2260 TopoDS_Shape aLine2 = aRefLine2->GetValue();
2261 if (aLine1.IsNull() || aLine2.IsNull() ||
2262 aLine1.ShapeType() != TopAbs_EDGE ||
2263 aLine2.ShapeType() != TopAbs_EDGE)
2265 SetErrorCode("Two edges must be given");
2271 TopoDS_Edge E1 = TopoDS::Edge(aLine1);
2272 TopoDS_Edge E2 = TopoDS::Edge(aLine2);
2275 Handle(Geom_Curve) C1 = BRep_Tool::Curve(E1,fp,lp);
2276 Handle(Geom_Curve) C2 = BRep_Tool::Curve(E2,fp,lp);
2278 if ( C1.IsNull() || C2.IsNull() ||
2279 !C1->IsKind(STANDARD_TYPE(Geom_Line)) ||
2280 !C2->IsKind(STANDARD_TYPE(Geom_Line)))
2282 SetErrorCode("The edges must be linear");
2286 Handle(Geom_Line) L1 = Handle(Geom_Line)::DownCast(C1);
2287 Handle(Geom_Line) L2 = Handle(Geom_Line)::DownCast(C2);
2289 gp_Lin aLin1 = L1->Lin();
2290 gp_Lin aLin2 = L2->Lin();
2292 anAngle = aLin1.Angle(aLin2);
2293 anAngle *= 180. / M_PI; // convert radians into degrees
2295 if (anAngle > 90.0) {
2296 anAngle = 180.0 - anAngle;
2301 catch (Standard_Failure)
2303 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
2304 SetErrorCode(aFail->GetMessageString());
2310 //=======================================================================
2312 * Compute angle (in degrees) between two vectors
2314 //=======================================================================
2315 Standard_Real GEOMImpl_IMeasureOperations::GetAngleBtwVectors (Handle(GEOM_Object) theVec1,
2316 Handle(GEOM_Object) theVec2)
2320 Standard_Real anAngle = -1.0;
2322 if (theVec1.IsNull() || theVec2.IsNull())
2325 if (theVec1->GetType() != GEOM_VECTOR || theVec2->GetType() != GEOM_VECTOR) {
2326 SetErrorCode("Two vectors must be given");
2330 Handle(GEOM_Function) aRefVec1 = theVec1->GetLastFunction();
2331 Handle(GEOM_Function) aRefVec2 = theVec2->GetLastFunction();
2332 if (aRefVec1.IsNull() || aRefVec2.IsNull())
2335 TopoDS_Shape aVec1 = aRefVec1->GetValue();
2336 TopoDS_Shape aVec2 = aRefVec2->GetValue();
2337 if (aVec1.IsNull() || aVec2.IsNull() ||
2338 aVec1.ShapeType() != TopAbs_EDGE ||
2339 aVec2.ShapeType() != TopAbs_EDGE)
2341 SetErrorCode("Two edges must be given");
2347 TopoDS_Edge aE1 = TopoDS::Edge(aVec1);
2348 TopoDS_Edge aE2 = TopoDS::Edge(aVec2);
2350 TopoDS_Vertex aP11, aP12, aP21, aP22;
2351 TopExp::Vertices(aE1, aP11, aP12, Standard_True);
2352 TopExp::Vertices(aE2, aP21, aP22, Standard_True);
2353 if (aP11.IsNull() || aP12.IsNull() || aP21.IsNull() || aP22.IsNull()) {
2354 SetErrorCode("Bad edge given");
2358 gp_Vec aV1 (BRep_Tool::Pnt(aP11), BRep_Tool::Pnt(aP12));
2359 gp_Vec aV2 (BRep_Tool::Pnt(aP21), BRep_Tool::Pnt(aP22)) ;
2361 anAngle = aV1.Angle(aV2);
2362 anAngle *= 180. / M_PI; // convert radians into degrees
2366 catch (Standard_Failure)
2368 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
2369 SetErrorCode(aFail->GetMessageString());
2376 //=============================================================================
2378 * CurveCurvatureByParam
2380 //=============================================================================
2381 Standard_Real GEOMImpl_IMeasureOperations::CurveCurvatureByParam
2382 (Handle(GEOM_Object) theCurve, Standard_Real& theParam)
2385 Standard_Real aRes = -1.0;
2387 if(theCurve.IsNull()) return aRes;
2389 Handle(GEOM_Function) aRefShape = theCurve->GetLastFunction();
2390 if(aRefShape.IsNull()) return aRes;
2392 TopoDS_Shape aShape = aRefShape->GetValue();
2393 if(aShape.IsNull()) {
2394 SetErrorCode("One of Objects has NULL Shape");
2398 Standard_Real aFP, aLP, aP;
2399 Handle(Geom_Curve) aCurve = BRep_Tool::Curve(TopoDS::Edge(aShape), aFP, aLP);
2400 aP = aFP + (aLP - aFP) * theParam;
2402 if(aCurve.IsNull()) return aRes;
2407 GeomLProp_CLProps Prop = GeomLProp_CLProps
2408 (aCurve, aP, 2, Precision::Confusion());
2409 aRes = fabs(Prop.Curvature());
2412 catch (Standard_Failure) {
2413 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
2414 SetErrorCode(aFail->GetMessageString());
2418 if( aRes > Precision::Confusion() )
2427 //=============================================================================
2429 * CurveCurvatureByPoint
2431 //=============================================================================
2432 Standard_Real GEOMImpl_IMeasureOperations::CurveCurvatureByPoint
2433 (Handle(GEOM_Object) theCurve, Handle(GEOM_Object) thePoint)
2436 Standard_Real aRes = -1.0;
2438 if( theCurve.IsNull() || thePoint.IsNull() ) return aRes;
2440 Handle(GEOM_Function) aRefCurve = theCurve->GetLastFunction();
2441 Handle(GEOM_Function) aRefPoint = thePoint->GetLastFunction();
2442 if( aRefCurve.IsNull() || aRefPoint.IsNull() ) return aRes;
2444 TopoDS_Edge anEdge = TopoDS::Edge(aRefCurve->GetValue());
2445 TopoDS_Vertex aPnt = TopoDS::Vertex(aRefPoint->GetValue());
2446 if( anEdge.IsNull() || aPnt.IsNull() ) {
2447 SetErrorCode("One of Objects has NULL Shape");
2451 Standard_Real aFP, aLP;
2452 Handle(Geom_Curve) aCurve = BRep_Tool::Curve(anEdge, aFP, aLP);
2453 if(aCurve.IsNull()) return aRes;
2454 gp_Pnt aPoint = BRep_Tool::Pnt(aPnt);
2459 GeomAPI_ProjectPointOnCurve PPCurve(aPoint, aCurve, aFP, aLP);
2460 if(PPCurve.NbPoints()>0) {
2461 GeomLProp_CLProps Prop = GeomLProp_CLProps
2462 (aCurve, PPCurve.LowerDistanceParameter(), 2, Precision::Confusion());
2463 aRes = fabs(Prop.Curvature());
2467 catch (Standard_Failure) {
2468 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
2469 SetErrorCode(aFail->GetMessageString());
2473 if( aRes > Precision::Confusion() )
2482 //=============================================================================
2484 * getSurfaceCurvatures
2486 //=============================================================================
2487 Standard_Real GEOMImpl_IMeasureOperations::getSurfaceCurvatures
2488 (const Handle(Geom_Surface)& aSurf,
2489 Standard_Real theUParam,
2490 Standard_Real theVParam,
2491 Standard_Boolean theNeedMaxCurv)
2494 Standard_Real aRes = 1.0;
2496 if (aSurf.IsNull()) return aRes;
2500 GeomLProp_SLProps Prop = GeomLProp_SLProps
2501 (aSurf, theUParam, theVParam, 2, Precision::Confusion());
2502 if(Prop.IsCurvatureDefined()) {
2503 if(Prop.IsUmbilic()) {
2504 //cout<<"is umbilic"<<endl;
2505 aRes = fabs(Prop.MeanCurvature());
2508 //cout<<"is not umbilic"<<endl;
2509 double c1 = fabs(Prop.MaxCurvature());
2510 double c2 = fabs(Prop.MinCurvature());
2519 catch (Standard_Failure) {
2520 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
2521 SetErrorCode(aFail->GetMessageString());
2525 if( fabs(aRes) > Precision::Confusion() )
2534 //=============================================================================
2536 * MaxSurfaceCurvatureByParam
2538 //=============================================================================
2539 Standard_Real GEOMImpl_IMeasureOperations::MaxSurfaceCurvatureByParam
2540 (Handle(GEOM_Object) theSurf,
2541 Standard_Real& theUParam,
2542 Standard_Real& theVParam)
2545 Standard_Real aRes = -1.0;
2547 if (theSurf.IsNull()) return aRes;
2549 Handle(GEOM_Function) aRefShape = theSurf->GetLastFunction();
2550 if(aRefShape.IsNull()) return aRes;
2552 TopoDS_Shape aShape = aRefShape->GetValue();
2553 if(aShape.IsNull()) {
2554 SetErrorCode("One of Objects has NULL Shape");
2558 TopoDS_Face F = TopoDS::Face(aShape);
2559 Handle(Geom_Surface) aSurf = BRep_Tool::Surface(F);
2561 //Compute the parameters
2562 Standard_Real U1,U2,V1,V2;
2563 ShapeAnalysis::GetFaceUVBounds(F,U1,U2,V1,V2);
2564 Standard_Real U = U1 + (U2-U1)*theUParam;
2565 Standard_Real V = V1 + (V2-V1)*theVParam;
2567 return getSurfaceCurvatures(aSurf, U, V, true);
2571 //=============================================================================
2573 * MaxSurfaceCurvatureByPoint
2575 //=============================================================================
2576 Standard_Real GEOMImpl_IMeasureOperations::MaxSurfaceCurvatureByPoint
2577 (Handle(GEOM_Object) theSurf, Handle(GEOM_Object) thePoint)
2580 Standard_Real aRes = -1.0;
2582 if( theSurf.IsNull() || thePoint.IsNull() ) return aRes;
2584 Handle(GEOM_Function) aRefShape = theSurf->GetLastFunction();
2585 Handle(GEOM_Function) aRefPoint = thePoint->GetLastFunction();
2586 if( aRefShape.IsNull() || aRefPoint.IsNull() ) return aRes;
2588 TopoDS_Face aFace = TopoDS::Face(aRefShape->GetValue());
2589 TopoDS_Vertex aPnt = TopoDS::Vertex(aRefPoint->GetValue());
2590 if( aFace.IsNull() || aPnt.IsNull() ) {
2591 SetErrorCode("One of Objects has NULL Shape");
2595 Handle(Geom_Surface) aSurf = BRep_Tool::Surface(aFace);
2596 if(aSurf.IsNull()) return aRes;
2597 gp_Pnt aPoint = BRep_Tool::Pnt(aPnt);
2599 //Compute the parameters
2600 ShapeAnalysis_Surface sas(aSurf);
2601 gp_Pnt2d UV = sas.ValueOfUV(aPoint,Precision::Confusion());
2603 return getSurfaceCurvatures(aSurf, UV.X(), UV.Y(), true);
2607 //=============================================================================
2609 * MinSurfaceCurvatureByParam
2611 //=============================================================================
2612 Standard_Real GEOMImpl_IMeasureOperations::MinSurfaceCurvatureByParam
2613 (Handle(GEOM_Object) theSurf,
2614 Standard_Real& theUParam,
2615 Standard_Real& theVParam)
2618 Standard_Real aRes = -1.0;
2620 if (theSurf.IsNull()) return aRes;
2622 Handle(GEOM_Function) aRefShape = theSurf->GetLastFunction();
2623 if(aRefShape.IsNull()) return aRes;
2625 TopoDS_Shape aShape = aRefShape->GetValue();
2626 if(aShape.IsNull()) {
2627 SetErrorCode("One of Objects has NULL Shape");
2631 TopoDS_Face F = TopoDS::Face(aShape);
2632 Handle(Geom_Surface) aSurf = BRep_Tool::Surface(F);
2634 //Compute the parameters
2635 Standard_Real U1,U2,V1,V2;
2636 ShapeAnalysis::GetFaceUVBounds(F,U1,U2,V1,V2);
2637 Standard_Real U = U1 + (U2-U1)*theUParam;
2638 Standard_Real V = V1 + (V2-V1)*theVParam;
2640 return getSurfaceCurvatures(aSurf, U, V, false);
2644 //=============================================================================
2646 * MinSurfaceCurvatureByPoint
2648 //=============================================================================
2649 Standard_Real GEOMImpl_IMeasureOperations::MinSurfaceCurvatureByPoint
2650 (Handle(GEOM_Object) theSurf, Handle(GEOM_Object) thePoint)
2653 Standard_Real aRes = -1.0;
2655 if( theSurf.IsNull() || thePoint.IsNull() ) return aRes;
2657 Handle(GEOM_Function) aRefShape = theSurf->GetLastFunction();
2658 Handle(GEOM_Function) aRefPoint = thePoint->GetLastFunction();
2659 if( aRefShape.IsNull() || aRefPoint.IsNull() ) return aRes;
2661 TopoDS_Face aFace = TopoDS::Face(aRefShape->GetValue());
2662 TopoDS_Vertex aPnt = TopoDS::Vertex(aRefPoint->GetValue());
2663 if( aFace.IsNull() || aPnt.IsNull() ) {
2664 SetErrorCode("One of Objects has NULL Shape");
2668 Handle(Geom_Surface) aSurf = BRep_Tool::Surface(aFace);
2669 if(aSurf.IsNull()) return aRes;
2670 gp_Pnt aPoint = BRep_Tool::Pnt(aPnt);
2672 //Compute the parameters
2673 ShapeAnalysis_Surface sas(aSurf);
2674 gp_Pnt2d UV = sas.ValueOfUV(aPoint,Precision::Confusion());
2676 return getSurfaceCurvatures(aSurf, UV.X(), UV.Y(), false);
2679 //=======================================================================
2680 //function : FillErrorsSub
2681 //purpose : Fill the errors list of subshapes on shape.
2682 //=======================================================================
2683 void GEOMImpl_IMeasureOperations::FillErrorsSub
2684 (const BRepCheck_Analyzer &theAna,
2685 const TopoDS_Shape &theShape,
2686 const TopAbs_ShapeEnum theSubType,
2687 TopTools_DataMapOfIntegerListOfShape &theMapErrors) const
2689 TopExp_Explorer anExp(theShape, theSubType);
2690 TopTools_MapOfShape aMapSubShapes;
2692 for (; anExp.More(); anExp.Next()) {
2693 const TopoDS_Shape &aSubShape = anExp.Current();
2695 if (aMapSubShapes.Add(aSubShape)) {
2696 const Handle(BRepCheck_Result) &aRes = theAna.Result(aSubShape);
2698 for (aRes->InitContextIterator();
2699 aRes->MoreShapeInContext();
2700 aRes->NextShapeInContext()) {
2701 if (aRes->ContextualShape().IsSame(theShape)) {
2702 BRepCheck_ListIteratorOfListOfStatus itl(aRes->StatusOnShape());
2704 if (itl.Value() != BRepCheck_NoError) {
2705 // Add all errors for theShape and its sub-shape.
2706 for (;itl.More(); itl.Next()) {
2707 const Standard_Integer aStat = (Standard_Integer)itl.Value();
2709 if (!theMapErrors.IsBound(aStat)) {
2710 TopTools_ListOfShape anEmpty;
2712 theMapErrors.Bind(aStat, anEmpty);
2715 TopTools_ListOfShape &theShapes = theMapErrors.ChangeFind(aStat);
2717 theShapes.Append(aSubShape);
2718 theShapes.Append(theShape);
2729 //=======================================================================
2730 //function : FillErrors
2731 //purpose : Fill the errors list.
2732 //=======================================================================
2733 void GEOMImpl_IMeasureOperations::FillErrors
2734 (const BRepCheck_Analyzer &theAna,
2735 const TopoDS_Shape &theShape,
2736 TopTools_DataMapOfIntegerListOfShape &theMapErrors,
2737 TopTools_MapOfShape &theMapShapes) const
2739 if (theMapShapes.Add(theShape)) {
2740 // Fill errors of child shapes.
2741 for (TopoDS_Iterator iter(theShape); iter.More(); iter.Next()) {
2742 FillErrors(theAna, iter.Value(), theMapErrors, theMapShapes);
2745 // Fill errors of theShape.
2746 const Handle(BRepCheck_Result) &aRes = theAna.Result(theShape);
2748 if (!aRes.IsNull()) {
2749 BRepCheck_ListIteratorOfListOfStatus itl(aRes->Status());
2751 if (itl.Value() != BRepCheck_NoError) {
2752 // Add all errors for theShape.
2753 for (;itl.More(); itl.Next()) {
2754 const Standard_Integer aStat = (Standard_Integer)itl.Value();
2756 if (!theMapErrors.IsBound(aStat)) {
2757 TopTools_ListOfShape anEmpty;
2759 theMapErrors.Bind(aStat, anEmpty);
2762 theMapErrors.ChangeFind(aStat).Append(theShape);
2767 // Add errors of subshapes on theShape.
2768 const TopAbs_ShapeEnum aType = theShape.ShapeType();
2772 FillErrorsSub(theAna, theShape, TopAbs_VERTEX, theMapErrors);
2775 FillErrorsSub(theAna, theShape, TopAbs_WIRE, theMapErrors);
2776 FillErrorsSub(theAna, theShape, TopAbs_EDGE, theMapErrors);
2777 FillErrorsSub(theAna, theShape, TopAbs_VERTEX, theMapErrors);
2780 FillErrorsSub(theAna, theShape, TopAbs_SHELL, theMapErrors);
2788 //=======================================================================
2789 //function : FillErrors
2790 //purpose : Fill the errors list.
2791 //=======================================================================
2792 void GEOMImpl_IMeasureOperations::FillErrors
2793 (const BRepCheck_Analyzer &theAna,
2794 const TopoDS_Shape &theShape,
2795 std::list<ShapeError> &theErrors) const
2797 // Fill the errors map.
2798 TopTools_DataMapOfIntegerListOfShape aMapErrors;
2799 TopTools_MapOfShape aMapShapes;
2801 FillErrors(theAna, theShape, aMapErrors, aMapShapes);
2803 // Map sub-shapes and their indices
2804 TopTools_IndexedMapOfShape anIndices;
2806 TopExp::MapShapes(theShape, anIndices);
2808 TopTools_DataMapIteratorOfDataMapOfIntegerListOfShape aMapIter(aMapErrors);
2810 for (; aMapIter.More(); aMapIter.Next()) {
2813 anError.error = (BRepCheck_Status)aMapIter.Key();
2815 TopTools_ListIteratorOfListOfShape aListIter(aMapIter.Value());
2816 TopTools_MapOfShape aMapUnique;
2818 for (; aListIter.More(); aListIter.Next()) {
2819 const TopoDS_Shape &aShape = aListIter.Value();
2821 if (aMapUnique.Add(aShape)) {
2822 const Standard_Integer anIndex = anIndices.FindIndex(aShape);
2824 anError.incriminated.push_back(anIndex);
2828 if (!anError.incriminated.empty()) {
2829 theErrors.push_back(anError);