1 // Copyright (C) 2007-2014 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 <GEOMImpl_IMeasureOperations.hxx>
24 #include <GEOMImpl_IMeasure.hxx>
25 #include <GEOMImpl_MeasureDriver.hxx>
26 #include <GEOMImpl_Types.hxx>
28 #include <GEOMUtils.hxx>
30 #include <GEOMAlgo_AlgoTools.hxx>
31 #include <GEOMAlgo_KindOfName.hxx>
32 #include <GEOMAlgo_ShapeInfoFiller.hxx>
34 #include <GEOM_PythonDump.hxx>
36 #include <utilities.h>
39 #include <Bnd_Box.hxx>
40 #include <BOPAlgo_CheckerSI.hxx>
41 #include <BOPCol_ListOfShape.hxx>
42 #include <BOPDS_DS.hxx>
43 #include <BOPDS_MapOfPassKey.hxx>
44 #include <BRepBndLib.hxx>
45 #include <BRepBuilderAPI_Copy.hxx>
46 #include <BRepCheck_ListIteratorOfListOfStatus.hxx>
47 #include <BRepCheck_Shell.hxx>
48 #include <BRepClass3d_SolidClassifier.hxx>
49 #include <BRepClass_FaceClassifier.hxx>
50 #include <BRepExtrema_DistShapeShape.hxx>
51 #include <BRepGProp.hxx>
52 #include <BRepTools.hxx>
53 #include <BRep_Tool.hxx>
54 #include <Geom_Line.hxx>
55 #include <GeomAPI_ProjectPointOnCurve.hxx>
56 #include <GeomAPI_ProjectPointOnSurf.hxx>
57 #include <GeomLProp_CLProps.hxx>
58 #include <GeomLProp_SLProps.hxx>
59 #include <GProp_GProps.hxx>
60 #include <GProp_PrincipalProps.hxx>
61 #include <ShapeAnalysis.hxx>
62 #include <ShapeAnalysis_Surface.hxx>
64 #include <TopExp_Explorer.hxx>
66 #include <TopoDS_Edge.hxx>
67 #include <TopTools_IndexedMapOfShape.hxx>
68 #include <TopTools_DataMapIteratorOfDataMapOfIntegerListOfShape.hxx>
69 #include <TopTools_ListIteratorOfListOfShape.hxx>
70 #include <TopTools_ListOfShape.hxx>
71 #include <Standard_ErrorHandler.hxx> // CAREFUL ! position of this file is critic : see Lucien PIGNOLONI / OCC
73 //=============================================================================
77 //=============================================================================
78 GEOMImpl_IMeasureOperations::GEOMImpl_IMeasureOperations (GEOM_Engine* theEngine, int theDocID)
79 : GEOM_IOperations(theEngine, theDocID)
81 MESSAGE("GEOMImpl_IMeasureOperations::GEOMImpl_IMeasureOperations");
84 //=============================================================================
88 //=============================================================================
89 GEOMImpl_IMeasureOperations::~GEOMImpl_IMeasureOperations()
91 MESSAGE("GEOMImpl_IMeasureOperations::~GEOMImpl_IMeasureOperations");
94 //=============================================================================
95 /*! Get kind and parameters of the given shape.
97 //=============================================================================
98 GEOMImpl_IMeasureOperations::ShapeKind GEOMImpl_IMeasureOperations::KindOfShape
99 (Handle(GEOM_Object) theShape,
100 Handle(TColStd_HSequenceOfInteger)& theIntegers,
101 Handle(TColStd_HSequenceOfReal)& theDoubles)
104 ShapeKind aKind = SK_NO_SHAPE;
106 if (theIntegers.IsNull()) theIntegers = new TColStd_HSequenceOfInteger;
107 else theIntegers->Clear();
109 if (theDoubles.IsNull()) theDoubles = new TColStd_HSequenceOfReal;
110 else theDoubles->Clear();
112 if (theShape.IsNull())
115 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
116 if (aRefShape.IsNull()) return aKind;
118 TopoDS_Shape aShape = aRefShape->GetValue();
119 if (aShape.IsNull()) return aKind;
121 int geom_type = theShape->GetType();
123 // check if it's advanced shape
124 if ( geom_type > ADVANCED_BASE ) {
130 GEOMAlgo_ShapeInfoFiller aSF;
131 aSF.SetShape(aShape);
133 Standard_Integer iErr = aSF.ErrorStatus();
135 SetErrorCode("Error in GEOMAlgo_ShapeInfoFiller");
138 const GEOMAlgo_ShapeInfo& anInfo = aSF.Info();
140 // Interprete results
141 TopAbs_ShapeEnum aType = anInfo.Type();
144 case TopAbs_COMPOUND:
145 case TopAbs_COMPSOLID:
147 // (+) geompy.kind.COMPOUND nb_solids nb_faces nb_edges nb_vertices
148 // (+) geompy.kind.COMPSOLID nb_solids nb_faces nb_edges nb_vertices
149 // ??? "nb_faces" - all faces or only 'standalone' faces?
150 if (aType == TopAbs_COMPOUND)
153 aKind = SK_COMPSOLID;
155 //theIntegers->Append(anInfo.NbSubShapes(TopAbs_COMPOUND));
156 //theIntegers->Append(anInfo.NbSubShapes(TopAbs_COMPSOLID));
157 theIntegers->Append(anInfo.NbSubShapes(TopAbs_SOLID));
158 theIntegers->Append(anInfo.NbSubShapes(TopAbs_FACE));
159 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
160 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
166 // (+) geompy.kind.SHELL geompy.info.closed nb_faces nb_edges nb_vertices
167 // (+) geompy.kind.SHELL geompy.info.unclosed nb_faces nb_edges nb_vertices
170 theIntegers->Append((int)anInfo.KindOfClosed());
172 theIntegers->Append(anInfo.NbSubShapes(TopAbs_FACE));
173 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
174 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
180 // (+) geompy.kind.WIRE geompy.info.closed nb_edges nb_vertices
181 // (+) geompy.kind.WIRE geompy.info.unclosed nb_edges nb_vertices
184 theIntegers->Append((int)anInfo.KindOfClosed());
186 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
187 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
195 GEOMAlgo_KindOfName aKN = anInfo.KindOfName();
198 case GEOMAlgo_KN_SPHERE:
199 // (+) geompy.kind.SPHERE xc yc zc R
203 gp_Pnt aC = anInfo.Location();
204 theDoubles->Append(aC.X());
205 theDoubles->Append(aC.Y());
206 theDoubles->Append(aC.Z());
208 theDoubles->Append(anInfo.Radius1());
211 case GEOMAlgo_KN_CYLINDER:
212 // (+) geompy.kind.CYLINDER xb yb zb dx dy dz R H
216 gp_Pnt aC = anInfo.Location();
217 theDoubles->Append(aC.X());
218 theDoubles->Append(aC.Y());
219 theDoubles->Append(aC.Z());
221 gp_Ax3 anAx3 = anInfo.Position();
222 gp_Dir aD = anAx3.Direction();
223 theDoubles->Append(aD.X());
224 theDoubles->Append(aD.Y());
225 theDoubles->Append(aD.Z());
227 theDoubles->Append(anInfo.Radius1());
228 theDoubles->Append(anInfo.Height());
231 case GEOMAlgo_KN_BOX:
232 // (+) geompy.kind.BOX xc yc zc ax ay az
236 gp_Pnt aC = anInfo.Location();
237 theDoubles->Append(aC.X());
238 theDoubles->Append(aC.Y());
239 theDoubles->Append(aC.Z());
241 gp_Ax3 anAx3 = anInfo.Position();
242 gp_Dir aD = anAx3.Direction();
243 gp_Dir aX = anAx3.XDirection();
246 if (aD.IsParallel(gp::DZ(), Precision::Angular()) &&
247 aX.IsParallel(gp::DX(), Precision::Angular())) {
248 theDoubles->Append(anInfo.Length()); // ax'
249 theDoubles->Append(anInfo.Width()); // ay'
250 theDoubles->Append(anInfo.Height()); // az'
252 else if (aD.IsParallel(gp::DZ(), Precision::Angular()) &&
253 aX.IsParallel(gp::DY(), Precision::Angular())) {
254 theDoubles->Append(anInfo.Width()); // ay'
255 theDoubles->Append(anInfo.Length()); // ax'
256 theDoubles->Append(anInfo.Height()); // az'
258 else if (aD.IsParallel(gp::DX(), Precision::Angular()) &&
259 aX.IsParallel(gp::DZ(), Precision::Angular())) {
260 theDoubles->Append(anInfo.Height()); // az'
261 theDoubles->Append(anInfo.Width()); // ay'
262 theDoubles->Append(anInfo.Length()); // ax'
264 else if (aD.IsParallel(gp::DX(), Precision::Angular()) &&
265 aX.IsParallel(gp::DY(), Precision::Angular())) {
266 theDoubles->Append(anInfo.Height()); // az'
267 theDoubles->Append(anInfo.Length()); // ax'
268 theDoubles->Append(anInfo.Width()); // ay'
270 else if (aD.IsParallel(gp::DY(), Precision::Angular()) &&
271 aX.IsParallel(gp::DZ(), Precision::Angular())) {
272 theDoubles->Append(anInfo.Width()); // ay'
273 theDoubles->Append(anInfo.Height()); // az'
274 theDoubles->Append(anInfo.Length()); // ax'
276 else if (aD.IsParallel(gp::DY(), Precision::Angular()) &&
277 aX.IsParallel(gp::DX(), Precision::Angular())) {
278 theDoubles->Append(anInfo.Length()); // ax'
279 theDoubles->Append(anInfo.Height()); // az'
280 theDoubles->Append(anInfo.Width()); // ay'
283 // (+) geompy.kind.ROTATED_BOX xo yo zo zx zy zz xx xy xz ax ay az
284 aKind = SK_ROTATED_BOX;
286 // Direction and XDirection
287 theDoubles->Append(aD.X());
288 theDoubles->Append(aD.Y());
289 theDoubles->Append(aD.Z());
291 theDoubles->Append(aX.X());
292 theDoubles->Append(aX.Y());
293 theDoubles->Append(aX.Z());
296 theDoubles->Append(anInfo.Length());
297 theDoubles->Append(anInfo.Width());
298 theDoubles->Append(anInfo.Height());
302 case GEOMAlgo_KN_TORUS:
303 // (+) geompy.kind.TORUS xc yc zc dx dy dz R_1 R_2
307 gp_Pnt aO = anInfo.Location();
308 theDoubles->Append(aO.X());
309 theDoubles->Append(aO.Y());
310 theDoubles->Append(aO.Z());
312 gp_Ax3 anAx3 = anInfo.Position();
313 gp_Dir aD = anAx3.Direction();
314 theDoubles->Append(aD.X());
315 theDoubles->Append(aD.Y());
316 theDoubles->Append(aD.Z());
318 theDoubles->Append(anInfo.Radius1());
319 theDoubles->Append(anInfo.Radius2());
322 case GEOMAlgo_KN_CONE:
323 // (+) geompy.kind.CONE xb yb zb dx dy dz R_1 R_2 H
327 gp_Pnt aO = anInfo.Location();
328 theDoubles->Append(aO.X());
329 theDoubles->Append(aO.Y());
330 theDoubles->Append(aO.Z());
332 gp_Ax3 anAx3 = anInfo.Position();
333 gp_Dir aD = anAx3.Direction();
334 theDoubles->Append(aD.X());
335 theDoubles->Append(aD.Y());
336 theDoubles->Append(aD.Z());
338 theDoubles->Append(anInfo.Radius1());
339 theDoubles->Append(anInfo.Radius2());
340 theDoubles->Append(anInfo.Height());
343 case GEOMAlgo_KN_POLYHEDRON:
344 // (+) geompy.kind.POLYHEDRON nb_faces nb_edges nb_vertices
346 aKind = SK_POLYHEDRON;
348 theIntegers->Append(anInfo.NbSubShapes(TopAbs_FACE));
349 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
350 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
354 // (+) geompy.kind.SOLID nb_faces nb_edges nb_vertices
356 theIntegers->Append(anInfo.NbSubShapes(TopAbs_FACE));
357 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
358 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
368 GEOMAlgo_KindOfName aKN = anInfo.KindOfName();
370 case GEOMAlgo_KN_SPHERE:
371 // (+) geompy.kind.SPHERE2D xc yc zc R
375 gp_Pnt aC = anInfo.Location();
376 theDoubles->Append(aC.X());
377 theDoubles->Append(aC.Y());
378 theDoubles->Append(aC.Z());
380 theDoubles->Append(anInfo.Radius1());
383 case GEOMAlgo_KN_CYLINDER:
384 // (+) geompy.kind.CYLINDER2D xb yb zb dx dy dz R H
386 aKind = SK_CYLINDER2D;
388 gp_Pnt aO = anInfo.Location();
389 theDoubles->Append(aO.X());
390 theDoubles->Append(aO.Y());
391 theDoubles->Append(aO.Z());
393 gp_Ax3 anAx3 = anInfo.Position();
394 gp_Dir aD = anAx3.Direction();
395 theDoubles->Append(aD.X());
396 theDoubles->Append(aD.Y());
397 theDoubles->Append(aD.Z());
399 theDoubles->Append(anInfo.Radius1());
400 theDoubles->Append(anInfo.Height());
403 case GEOMAlgo_KN_TORUS:
404 // (+) geompy.kind.TORUS2D xc yc zc dx dy dz R_1 R_2
408 gp_Pnt aO = anInfo.Location();
409 theDoubles->Append(aO.X());
410 theDoubles->Append(aO.Y());
411 theDoubles->Append(aO.Z());
413 gp_Ax3 anAx3 = anInfo.Position();
414 gp_Dir aD = anAx3.Direction();
415 theDoubles->Append(aD.X());
416 theDoubles->Append(aD.Y());
417 theDoubles->Append(aD.Z());
419 theDoubles->Append(anInfo.Radius1());
420 theDoubles->Append(anInfo.Radius2());
423 case GEOMAlgo_KN_CONE:
424 // (+) geompy.kind.CONE2D xc yc zc dx dy dz R_1 R_2 H
428 gp_Pnt aO = anInfo.Location();
429 theDoubles->Append(aO.X());
430 theDoubles->Append(aO.Y());
431 theDoubles->Append(aO.Z());
433 gp_Ax3 anAx3 = anInfo.Position();
434 gp_Dir aD = anAx3.Direction();
435 theDoubles->Append(aD.X());
436 theDoubles->Append(aD.Y());
437 theDoubles->Append(aD.Z());
439 theDoubles->Append(anInfo.Radius1());
440 theDoubles->Append(anInfo.Radius2());
441 theDoubles->Append(anInfo.Height());
444 case GEOMAlgo_KN_DISKCIRCLE:
445 // (+) geompy.kind.DISK_CIRCLE xc yc zc dx dy dz R
447 aKind = SK_DISK_CIRCLE;
449 gp_Pnt aC = anInfo.Location();
450 theDoubles->Append(aC.X());
451 theDoubles->Append(aC.Y());
452 theDoubles->Append(aC.Z());
454 gp_Ax3 anAx3 = anInfo.Position();
455 gp_Dir aD = anAx3.Direction();
456 theDoubles->Append(aD.X());
457 theDoubles->Append(aD.Y());
458 theDoubles->Append(aD.Z());
460 theDoubles->Append(anInfo.Radius1());
463 case GEOMAlgo_KN_DISKELLIPSE:
464 // (+) geompy.kind.DISK_ELLIPSE xc yc zc dx dy dz R_1 R_2
466 aKind = SK_DISK_ELLIPSE;
468 gp_Pnt aC = anInfo.Location();
469 theDoubles->Append(aC.X());
470 theDoubles->Append(aC.Y());
471 theDoubles->Append(aC.Z());
473 gp_Ax3 anAx3 = anInfo.Position();
474 gp_Dir aD = anAx3.Direction();
475 theDoubles->Append(aD.X());
476 theDoubles->Append(aD.Y());
477 theDoubles->Append(aD.Z());
479 theDoubles->Append(anInfo.Radius1());
480 theDoubles->Append(anInfo.Radius2());
483 case GEOMAlgo_KN_RECTANGLE:
484 case GEOMAlgo_KN_TRIANGLE:
485 case GEOMAlgo_KN_QUADRANGLE:
486 case GEOMAlgo_KN_POLYGON:
487 // (+) geompy.kind.POLYGON xo yo zo dx dy dz nb_edges nb_vertices
491 gp_Pnt aO = anInfo.Location();
492 theDoubles->Append(aO.X());
493 theDoubles->Append(aO.Y());
494 theDoubles->Append(aO.Z());
496 gp_Ax3 anAx3 = anInfo.Position();
497 gp_Dir aD = anAx3.Direction();
498 theDoubles->Append(aD.X());
499 theDoubles->Append(aD.Y());
500 theDoubles->Append(aD.Z());
502 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
503 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
506 case GEOMAlgo_KN_PLANE: // infinite
507 // (+) geompy.kind.PLANE xo yo zo dx dy dz
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());
524 if (anInfo.KindOfShape() == GEOMAlgo_KS_PLANE) {
525 // (+) geompy.kind.PLANAR xo yo zo dx dy dz nb_edges nb_vertices
529 gp_Pnt aC = anInfo.Location();
530 theDoubles->Append(aC.X());
531 theDoubles->Append(aC.Y());
532 theDoubles->Append(aC.Z());
534 gp_Ax3 anAx3 = anInfo.Position();
535 gp_Dir aD = anAx3.Direction();
536 theDoubles->Append(aD.X());
537 theDoubles->Append(aD.Y());
538 theDoubles->Append(aD.Z());
540 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
541 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
544 // ??? geompy.kind.FACE nb_edges nb_vertices _surface_type_id_
545 // (+) geompy.kind.FACE nb_edges nb_vertices
547 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
548 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
558 GEOMAlgo_KindOfName aKN = anInfo.KindOfName();
560 case GEOMAlgo_KN_CIRCLE:
562 // (+) geompy.kind.CIRCLE xc yc zc dx dy dz R
565 gp_Pnt aC = anInfo.Location();
566 theDoubles->Append(aC.X());
567 theDoubles->Append(aC.Y());
568 theDoubles->Append(aC.Z());
570 gp_Ax3 anAx3 = anInfo.Position();
571 gp_Dir aD = anAx3.Direction();
572 theDoubles->Append(aD.X());
573 theDoubles->Append(aD.Y());
574 theDoubles->Append(aD.Z());
576 theDoubles->Append(anInfo.Radius1());
579 case GEOMAlgo_KN_ARCCIRCLE:
581 // (+) geompy.kind.ARC_CIRCLE xc yc zc dx dy dz R x1 y1 z1 x2 y2 z2
582 aKind = SK_ARC_CIRCLE;
584 gp_Pnt aC = anInfo.Location();
585 theDoubles->Append(aC.X());
586 theDoubles->Append(aC.Y());
587 theDoubles->Append(aC.Z());
589 gp_Ax3 anAx3 = anInfo.Position();
590 gp_Dir aD = anAx3.Direction();
591 theDoubles->Append(aD.X());
592 theDoubles->Append(aD.Y());
593 theDoubles->Append(aD.Z());
595 theDoubles->Append(anInfo.Radius1());
597 gp_Pnt aP1 = anInfo.Pnt1();
598 theDoubles->Append(aP1.X());
599 theDoubles->Append(aP1.Y());
600 theDoubles->Append(aP1.Z());
602 gp_Pnt aP2 = anInfo.Pnt2();
603 theDoubles->Append(aP2.X());
604 theDoubles->Append(aP2.Y());
605 theDoubles->Append(aP2.Z());
608 case GEOMAlgo_KN_ELLIPSE:
610 // (+) geompy.kind.ELLIPSE xc yc zc dx dy dz R_1 R_2
613 gp_Pnt aC = anInfo.Location();
614 theDoubles->Append(aC.X());
615 theDoubles->Append(aC.Y());
616 theDoubles->Append(aC.Z());
618 gp_Ax3 anAx3 = anInfo.Position();
619 gp_Dir aD = anAx3.Direction();
620 theDoubles->Append(aD.X());
621 theDoubles->Append(aD.Y());
622 theDoubles->Append(aD.Z());
624 theDoubles->Append(anInfo.Radius1());
625 theDoubles->Append(anInfo.Radius2());
628 case GEOMAlgo_KN_ARCELLIPSE:
630 // (+) geompy.kind.ARC_ELLIPSE xc yc zc dx dy dz R_1 R_2 x1 y1 z1 x2 y2 z2
631 aKind = SK_ARC_ELLIPSE;
633 gp_Pnt aC = anInfo.Location();
634 theDoubles->Append(aC.X());
635 theDoubles->Append(aC.Y());
636 theDoubles->Append(aC.Z());
638 gp_Ax3 anAx3 = anInfo.Position();
639 gp_Dir aD = anAx3.Direction();
640 theDoubles->Append(aD.X());
641 theDoubles->Append(aD.Y());
642 theDoubles->Append(aD.Z());
644 theDoubles->Append(anInfo.Radius1());
645 theDoubles->Append(anInfo.Radius2());
647 gp_Pnt aP1 = anInfo.Pnt1();
648 theDoubles->Append(aP1.X());
649 theDoubles->Append(aP1.Y());
650 theDoubles->Append(aP1.Z());
652 gp_Pnt aP2 = anInfo.Pnt2();
653 theDoubles->Append(aP2.X());
654 theDoubles->Append(aP2.Y());
655 theDoubles->Append(aP2.Z());
658 case GEOMAlgo_KN_LINE:
660 // ??? geompy.kind.LINE x1 y1 z1 x2 y2 z2
661 // (+) geompy.kind.LINE x1 y1 z1 dx dy dz
664 gp_Pnt aO = anInfo.Location();
665 theDoubles->Append(aO.X());
666 theDoubles->Append(aO.Y());
667 theDoubles->Append(aO.Z());
669 gp_Dir aD = anInfo.Direction();
670 theDoubles->Append(aD.X());
671 theDoubles->Append(aD.Y());
672 theDoubles->Append(aD.Z());
675 case GEOMAlgo_KN_SEGMENT:
677 // (+) geompy.kind.SEGMENT x1 y1 z1 x2 y2 z2
680 gp_Pnt aP1 = anInfo.Pnt1();
681 theDoubles->Append(aP1.X());
682 theDoubles->Append(aP1.Y());
683 theDoubles->Append(aP1.Z());
685 gp_Pnt aP2 = anInfo.Pnt2();
686 theDoubles->Append(aP2.X());
687 theDoubles->Append(aP2.Y());
688 theDoubles->Append(aP2.Z());
692 // ??? geompy.kind.EDGE nb_vertices _curve_type_id_
693 // (+) geompy.kind.EDGE nb_vertices
694 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
701 // (+) geompy.kind.VERTEX x y z
704 gp_Pnt aP = anInfo.Location();
705 theDoubles->Append(aP.X());
706 theDoubles->Append(aP.Y());
707 theDoubles->Append(aP.Z());
716 //=============================================================================
720 //=============================================================================
721 void GEOMImpl_IMeasureOperations::GetPosition
722 (Handle(GEOM_Object) theShape,
723 Standard_Real& Ox, Standard_Real& Oy, Standard_Real& Oz,
724 Standard_Real& Zx, Standard_Real& Zy, Standard_Real& Zz,
725 Standard_Real& Xx, Standard_Real& Xy, Standard_Real& Xz)
729 //Set default values: global CS
730 Ox = Oy = Oz = Zx = Zy = Xy = Xz = 0.;
733 if (theShape.IsNull()) return;
735 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
736 if (aRefShape.IsNull()) return;
738 TopoDS_Shape aShape = aRefShape->GetValue();
739 if (aShape.IsNull()) {
740 SetErrorCode("The Objects has NULL Shape");
745 #if OCC_VERSION_LARGE > 0x06010000
749 gp_Ax3 anAx3 = GEOMUtils::GetPosition(aShape);
751 gp_Pnt anOri = anAx3.Location();
752 gp_Dir aDirZ = anAx3.Direction();
753 gp_Dir aDirX = anAx3.XDirection();
756 anOri.Coord(Ox, Oy, Oz);
757 aDirZ.Coord(Zx, Zy, Zz);
758 aDirX.Coord(Xx, Xy, Xz);
760 catch (Standard_Failure) {
761 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
762 SetErrorCode(aFail->GetMessageString());
769 //=============================================================================
773 //=============================================================================
774 Handle(GEOM_Object) GEOMImpl_IMeasureOperations::GetCentreOfMass
775 (Handle(GEOM_Object) theShape)
779 if (theShape.IsNull()) return NULL;
781 //Add a new CentreOfMass object
782 Handle(GEOM_Object) aCDG = GetEngine()->AddObject(GetDocID(), GEOM_CDG);
784 //Add a new CentreOfMass function
785 Handle(GEOM_Function) aFunction =
786 aCDG->AddFunction(GEOMImpl_MeasureDriver::GetID(), CDG_MEASURE);
787 if (aFunction.IsNull()) return NULL;
789 //Check if the function is set correctly
790 if (aFunction->GetDriverGUID() != GEOMImpl_MeasureDriver::GetID()) return NULL;
792 GEOMImpl_IMeasure aCI (aFunction);
794 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
795 if (aRefShape.IsNull()) return NULL;
797 aCI.SetBase(aRefShape);
799 //Compute the CentreOfMass value
801 #if OCC_VERSION_LARGE > 0x06010000
804 if (!GetSolver()->ComputeFunction(aFunction)) {
805 SetErrorCode("Measure driver failed to compute centre of mass");
809 catch (Standard_Failure) {
810 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
811 SetErrorCode(aFail->GetMessageString());
815 //Make a Python command
816 GEOM::TPythonDump(aFunction) << aCDG << " = geompy.MakeCDG(" << theShape << ")";
822 //=============================================================================
826 //=============================================================================
827 Handle(GEOM_Object) GEOMImpl_IMeasureOperations::GetVertexByIndex
828 (Handle(GEOM_Object) theShape,
829 Standard_Integer theIndex)
833 if (theShape.IsNull()) return NULL;
835 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
836 if (aRefShape.IsNull()) return NULL;
838 //Add a new Vertex object
839 Handle(GEOM_Object) aVertex = GetEngine()->AddObject(GetDocID(), GEOM_POINT);
842 Handle(GEOM_Function) aFunction =
843 aVertex->AddFunction(GEOMImpl_MeasureDriver::GetID(), VERTEX_BY_INDEX);
844 if (aFunction.IsNull()) return NULL;
846 //Check if the function is set correctly
847 if (aFunction->GetDriverGUID() != GEOMImpl_MeasureDriver::GetID()) return NULL;
849 GEOMImpl_IMeasure aCI (aFunction);
850 aCI.SetBase(aRefShape);
851 aCI.SetIndex(theIndex);
855 #if OCC_VERSION_LARGE > 0x06010000
858 if (!GetSolver()->ComputeFunction(aFunction)) {
859 SetErrorCode("Vertex by index driver failed.");
863 catch (Standard_Failure) {
864 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
865 SetErrorCode(aFail->GetMessageString());
869 //Make a Python command
870 GEOM::TPythonDump(aFunction) << aVertex << " = geompy.GetVertexByIndex(" << theShape << ", " << theIndex << ")";
876 //=============================================================================
880 //=============================================================================
881 Handle(GEOM_Object) GEOMImpl_IMeasureOperations::GetNormal
882 (Handle(GEOM_Object) theFace,
883 Handle(GEOM_Object) theOptionalPoint)
887 if (theFace.IsNull()) return NULL;
889 //Add a new Normale object
890 Handle(GEOM_Object) aNorm = GetEngine()->AddObject(GetDocID(), GEOM_VECTOR);
892 //Add a new Normale function
893 Handle(GEOM_Function) aFunction =
894 aNorm->AddFunction(GEOMImpl_MeasureDriver::GetID(), VECTOR_FACE_NORMALE);
895 if (aFunction.IsNull()) return NULL;
897 //Check if the function is set correctly
898 if (aFunction->GetDriverGUID() != GEOMImpl_MeasureDriver::GetID()) return NULL;
900 GEOMImpl_IMeasure aCI (aFunction);
902 Handle(GEOM_Function) aFace = theFace->GetLastFunction();
903 if (aFace.IsNull()) return NULL;
907 if (!theOptionalPoint.IsNull()) {
908 Handle(GEOM_Function) anOptPnt = theOptionalPoint->GetLastFunction();
909 aCI.SetPoint(anOptPnt);
912 //Compute the Normale value
914 #if OCC_VERSION_LARGE > 0x06010000
917 if (!GetSolver()->ComputeFunction(aFunction)) {
918 SetErrorCode("Measure driver failed to compute normake of face");
922 catch (Standard_Failure) {
923 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
924 SetErrorCode(aFail->GetMessageString());
928 //Make a Python command
929 GEOM::TPythonDump pd (aFunction);
930 pd << aNorm << " = geompy.GetNormal(" << theFace;
931 if (!theOptionalPoint.IsNull()) {
932 pd << ", " << theOptionalPoint;
940 //=============================================================================
944 //=============================================================================
945 void GEOMImpl_IMeasureOperations::GetBasicProperties (Handle(GEOM_Object) theShape,
946 Standard_Real& theLength,
947 Standard_Real& theSurfArea,
948 Standard_Real& theVolume)
952 if (theShape.IsNull()) return;
954 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
955 if (aRefShape.IsNull()) return;
957 TopoDS_Shape aShape = aRefShape->GetValue();
958 if (aShape.IsNull()) {
959 SetErrorCode("The Objects has NULL Shape");
963 //Compute the parameters
964 GProp_GProps LProps, SProps;
966 #if OCC_VERSION_LARGE > 0x06010000
969 BRepGProp::LinearProperties(aShape, LProps);
970 theLength = LProps.Mass();
972 BRepGProp::SurfaceProperties(aShape, SProps);
973 theSurfArea = SProps.Mass();
976 if (aShape.ShapeType() < TopAbs_SHELL) {
977 for (TopExp_Explorer Exp (aShape, TopAbs_SOLID); Exp.More(); Exp.Next()) {
979 BRepGProp::VolumeProperties(Exp.Current(), VProps);
980 theVolume += VProps.Mass();
984 catch (Standard_Failure) {
985 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
986 SetErrorCode(aFail->GetMessageString());
993 //=============================================================================
997 //=============================================================================
998 void GEOMImpl_IMeasureOperations::GetInertia
999 (Handle(GEOM_Object) theShape,
1000 Standard_Real& I11, Standard_Real& I12, Standard_Real& I13,
1001 Standard_Real& I21, Standard_Real& I22, Standard_Real& I23,
1002 Standard_Real& I31, Standard_Real& I32, Standard_Real& I33,
1003 Standard_Real& Ix , Standard_Real& Iy , Standard_Real& Iz)
1007 if (theShape.IsNull()) return;
1009 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1010 if (aRefShape.IsNull()) return;
1012 TopoDS_Shape aShape = aRefShape->GetValue();
1013 if (aShape.IsNull()) {
1014 SetErrorCode("The Objects has NULL Shape");
1018 //Compute the parameters
1019 GProp_GProps System;
1022 #if OCC_VERSION_LARGE > 0x06010000
1025 if (aShape.ShapeType() == TopAbs_VERTEX ||
1026 aShape.ShapeType() == TopAbs_EDGE ||
1027 aShape.ShapeType() == TopAbs_WIRE) {
1028 BRepGProp::LinearProperties(aShape, System);
1029 } else if (aShape.ShapeType() == TopAbs_FACE ||
1030 aShape.ShapeType() == TopAbs_SHELL) {
1031 BRepGProp::SurfaceProperties(aShape, System);
1033 BRepGProp::VolumeProperties(aShape, System);
1035 gp_Mat I = System.MatrixOfInertia();
1049 GProp_PrincipalProps Pr = System.PrincipalProperties();
1050 Pr.Moments(Ix,Iy,Iz);
1052 catch (Standard_Failure) {
1053 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1054 SetErrorCode(aFail->GetMessageString());
1061 //=============================================================================
1065 //=============================================================================
1066 void GEOMImpl_IMeasureOperations::GetBoundingBox
1067 (Handle(GEOM_Object) theShape,
1068 const Standard_Boolean precise,
1069 Standard_Real& Xmin, Standard_Real& Xmax,
1070 Standard_Real& Ymin, Standard_Real& Ymax,
1071 Standard_Real& Zmin, Standard_Real& Zmax)
1075 if (theShape.IsNull()) return;
1077 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1078 if (aRefShape.IsNull()) return;
1080 TopoDS_Shape aShape = aRefShape->GetValue();
1081 if (aShape.IsNull()) {
1082 SetErrorCode("The Objects has NULL Shape");
1086 //Compute the parameters
1090 #if OCC_VERSION_LARGE > 0x06010000
1093 BRepBuilderAPI_Copy aCopyTool (aShape);
1094 if (!aCopyTool.IsDone()) {
1095 SetErrorCode("GetBoundingBox Error: Bad shape detected");
1099 aShape = aCopyTool.Shape();
1101 // remove triangulation to obtain more exact boundaries
1102 BRepTools::Clean(aShape);
1104 BRepBndLib::Add(aShape, B);
1107 if (!GEOMUtils::PreciseBoundingBox(aShape, B)) {
1108 SetErrorCode("GetBoundingBox Error: Bounding box cannot be precised");
1113 B.Get(Xmin, Ymin, Zmin, Xmax, Ymax, Zmax);
1115 catch (Standard_Failure) {
1116 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1117 SetErrorCode(aFail->GetMessageString());
1124 //=============================================================================
1128 //=============================================================================
1129 Handle(GEOM_Object) GEOMImpl_IMeasureOperations::GetBoundingBox
1130 (Handle(GEOM_Object) theShape,
1131 const Standard_Boolean precise)
1135 if (theShape.IsNull()) return NULL;
1137 //Add a new BoundingBox object
1138 Handle(GEOM_Object) aBnd = GetEngine()->AddObject(GetDocID(), GEOM_BOX);
1140 //Add a new BoundingBox function
1141 const int aType = (precise ? BND_BOX_MEASURE_PRECISE : BND_BOX_MEASURE);
1142 Handle(GEOM_Function) aFunction =
1143 aBnd->AddFunction(GEOMImpl_MeasureDriver::GetID(), aType);
1144 if (aFunction.IsNull()) return NULL;
1146 //Check if the function is set correctly
1147 if (aFunction->GetDriverGUID() != GEOMImpl_MeasureDriver::GetID()) return NULL;
1149 GEOMImpl_IMeasure aCI (aFunction);
1151 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1152 if (aRefShape.IsNull()) return NULL;
1154 aCI.SetBase(aRefShape);
1156 //Compute the BoundingBox value
1158 #if OCC_VERSION_LARGE > 0x06010000
1161 if (!GetSolver()->ComputeFunction(aFunction)) {
1162 SetErrorCode("Measure driver failed to compute a bounding box");
1166 catch (Standard_Failure) {
1167 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1168 SetErrorCode(aFail->GetMessageString());
1172 //Make a Python command
1173 GEOM::TPythonDump aPd(aFunction);
1175 aPd << aBnd << " = geompy.MakeBoundingBox(" << theShape;
1187 //=============================================================================
1191 //=============================================================================
1192 void GEOMImpl_IMeasureOperations::GetTolerance
1193 (Handle(GEOM_Object) theShape,
1194 Standard_Real& FaceMin, Standard_Real& FaceMax,
1195 Standard_Real& EdgeMin, Standard_Real& EdgeMax,
1196 Standard_Real& VertMin, Standard_Real& VertMax)
1200 if (theShape.IsNull()) return;
1202 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1203 if (aRefShape.IsNull()) return;
1205 TopoDS_Shape aShape = aRefShape->GetValue();
1206 if (aShape.IsNull()) {
1207 SetErrorCode("The Objects has NULL Shape");
1211 //Compute the parameters
1213 FaceMin = EdgeMin = VertMin = RealLast();
1214 FaceMax = EdgeMax = VertMax = -RealLast();
1217 #if OCC_VERSION_LARGE > 0x06010000
1220 for (TopExp_Explorer ExF (aShape, TopAbs_FACE); ExF.More(); ExF.Next()) {
1221 TopoDS_Face Face = TopoDS::Face(ExF.Current());
1222 T = BRep_Tool::Tolerance(Face);
1228 for (TopExp_Explorer ExE (aShape, TopAbs_EDGE); ExE.More(); ExE.Next()) {
1229 TopoDS_Edge Edge = TopoDS::Edge(ExE.Current());
1230 T = BRep_Tool::Tolerance(Edge);
1236 for (TopExp_Explorer ExV (aShape, TopAbs_VERTEX); ExV.More(); ExV.Next()) {
1237 TopoDS_Vertex Vertex = TopoDS::Vertex(ExV.Current());
1238 T = BRep_Tool::Tolerance(Vertex);
1245 catch (Standard_Failure) {
1246 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1247 SetErrorCode(aFail->GetMessageString());
1254 //=============================================================================
1258 //=============================================================================
1259 bool GEOMImpl_IMeasureOperations::CheckShape (Handle(GEOM_Object) theShape,
1260 const Standard_Boolean theIsCheckGeom,
1261 std::list<ShapeError> &theErrors)
1266 if (theShape.IsNull()) return false;
1268 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1269 if (aRefShape.IsNull()) return false;
1271 TopoDS_Shape aShape = aRefShape->GetValue();
1272 if (aShape.IsNull()) {
1273 SetErrorCode("The Objects has NULL Shape");
1277 //Compute the parameters
1278 bool isValid = false;
1280 #if OCC_VERSION_LARGE > 0x06010000
1283 BRepCheck_Analyzer ana (aShape, theIsCheckGeom);
1284 if (ana.IsValid()) {
1287 FillErrors(ana, aShape, theErrors);
1290 catch (Standard_Failure) {
1291 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1292 SetErrorCode(aFail->GetMessageString());
1300 //=============================================================================
1304 //=============================================================================
1305 TCollection_AsciiString GEOMImpl_IMeasureOperations::PrintShapeErrors
1306 (Handle(GEOM_Object) theShape,
1307 const std::list<ShapeError> &theErrors)
1309 TCollection_AsciiString aDump;
1311 if (theShape.IsNull()) {
1315 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1317 if (aRefShape.IsNull()) {
1321 TopoDS_Shape aShape = aRefShape->GetValue();
1323 if (aShape.IsNull()) {
1324 SetErrorCode("The Objects has NULL Shape");
1328 if (!theErrors.empty()) {
1329 // The shape is not valid. Prepare errors for dump.
1330 TopTools_IndexedMapOfShape anIndices;
1331 std::list<ShapeError>::const_iterator anIter = theErrors.begin();
1332 Standard_Integer nbv, nbe, nbw, nbf, nbs, nbo;
1333 nbv = nbe = nbw = nbf = nbs = nbo = 0;
1335 // Map sub-shapes and their indices
1336 TopExp::MapShapes(aShape, anIndices);
1338 const Standard_Integer aNbSubShapes = anIndices.Extent();
1339 TColStd_MapOfInteger aMapPbInd;
1341 aDump += " -- The Shape has problems :\n";
1342 aDump += " Check Count\n";
1343 aDump += " ------------------------------------------------\n";
1345 for (; anIter != theErrors.end(); anIter++) {
1346 Standard_Integer aNbShapes = anIter->incriminated.size();
1348 switch(anIter->error) {
1349 case BRepCheck_InvalidPointOnCurve:
1350 aDump += " Invalid Point on Curve ................... ";
1352 case BRepCheck_InvalidPointOnCurveOnSurface:
1353 aDump += " Invalid Point on CurveOnSurface .......... ";
1355 case BRepCheck_InvalidPointOnSurface:
1356 aDump += " Invalid Point on Surface ................. ";
1358 case BRepCheck_No3DCurve:
1359 aDump += " No 3D Curve .............................. ";
1361 case BRepCheck_Multiple3DCurve:
1362 aDump += " Multiple 3D Curve ........................ ";
1364 case BRepCheck_Invalid3DCurve:
1365 aDump += " Invalid 3D Curve ......................... ";
1367 case BRepCheck_NoCurveOnSurface:
1368 aDump += " No Curve on Surface ...................... ";
1370 case BRepCheck_InvalidCurveOnSurface:
1371 aDump += " Invalid Curve on Surface ................. ";
1373 case BRepCheck_InvalidCurveOnClosedSurface:
1374 aDump += " Invalid Curve on closed Surface .......... ";
1376 case BRepCheck_InvalidSameRangeFlag:
1377 aDump += " Invalid SameRange Flag ................... ";
1379 case BRepCheck_InvalidSameParameterFlag:
1380 aDump += " Invalid SameParameter Flag ............... ";
1382 case BRepCheck_InvalidDegeneratedFlag:
1383 aDump += " Invalid Degenerated Flag ................. ";
1385 case BRepCheck_FreeEdge:
1386 aDump += " Free Edge ................................ ";
1388 case BRepCheck_InvalidMultiConnexity:
1389 aDump += " Invalid MultiConnexity ................... ";
1391 case BRepCheck_InvalidRange:
1392 aDump += " Invalid Range ............................ ";
1394 case BRepCheck_EmptyWire:
1395 aDump += " Empty Wire ............................... ";
1397 case BRepCheck_RedundantEdge:
1398 aDump += " Redundant Edge ........................... ";
1400 case BRepCheck_SelfIntersectingWire:
1401 aDump += " Self Intersecting Wire ................... ";
1403 case BRepCheck_NoSurface:
1404 aDump += " No Surface ............................... ";
1406 case BRepCheck_InvalidWire:
1407 aDump += " Invalid Wire ............................. ";
1409 case BRepCheck_RedundantWire:
1410 aDump += " Redundant Wire ........................... ";
1412 case BRepCheck_IntersectingWires:
1413 aDump += " Intersecting Wires ....................... ";
1415 case BRepCheck_InvalidImbricationOfWires:
1416 aDump += " Invalid Imbrication of Wires ............. ";
1418 case BRepCheck_EmptyShell:
1419 aDump += " Empty Shell .............................. ";
1421 case BRepCheck_RedundantFace:
1422 aDump += " Redundant Face ........................... ";
1424 case BRepCheck_UnorientableShape:
1425 aDump += " Unorientable Shape ....................... ";
1427 case BRepCheck_NotClosed:
1428 aDump += " Not Closed ............................... ";
1430 case BRepCheck_NotConnected:
1431 aDump += " Not Connected ............................ ";
1433 case BRepCheck_SubshapeNotInShape:
1434 aDump += " Sub-shape not in Shape ................... ";
1436 case BRepCheck_BadOrientation:
1437 aDump += " Bad Orientation .......................... ";
1439 case BRepCheck_BadOrientationOfSubshape:
1440 aDump += " Bad Orientation of Sub-shape ............. ";
1442 case BRepCheck_InvalidToleranceValue:
1443 aDump += " Invalid Tolerance Value .................. ";
1445 case BRepCheck_CheckFail:
1446 aDump += " Check Shape Failure ...................... ";
1452 aDump += TCollection_AsciiString(aNbShapes) + "\n";
1454 // Count types of shape.
1455 std::list<int>::const_iterator aShIter = anIter->incriminated.begin();
1457 for (; aShIter != anIter->incriminated.end(); aShIter++) {
1458 const Standard_Integer anIndex = *aShIter;
1460 if (anIndex > 0 && anIndex <= aNbSubShapes && aMapPbInd.Add(anIndex)) {
1461 const TopoDS_Shape &aSubShape = anIndices.FindKey(anIndex);
1462 const TopAbs_ShapeEnum aType = aSubShape.ShapeType();
1465 case TopAbs_VERTEX : nbv++; break;
1466 case TopAbs_EDGE : nbe++; break;
1467 case TopAbs_WIRE : nbw++; break;
1468 case TopAbs_FACE : nbf++; break;
1469 case TopAbs_SHELL : nbs++; break;
1470 case TopAbs_SOLID : nbo++; break;
1477 const Standard_Integer aNbFaultyShapes = nbv + nbe + nbw + nbf + nbs + nbo;
1478 aDump += " ------------------------------------------------\n";
1479 aDump += "*** Shapes with problems : ";
1480 aDump += TCollection_AsciiString(aNbFaultyShapes) + "\n";
1483 aDump += "VERTEX : ";
1484 if (nbv < 10) aDump += " ";
1485 aDump += TCollection_AsciiString(nbv) + "\n";
1489 if (nbe < 10) aDump += " ";
1490 aDump += TCollection_AsciiString(nbe) + "\n";
1494 if (nbw < 10) aDump += " ";
1495 aDump += TCollection_AsciiString(nbw) + "\n";
1499 if (nbf < 10) aDump += " ";
1500 aDump += TCollection_AsciiString(nbf) + "\n";
1503 aDump += "SHELL : ";
1504 if (nbs < 10) aDump += " ";
1505 aDump += TCollection_AsciiString(nbs) + "\n";
1508 aDump += "SOLID : ";
1509 if (nbo < 10) aDump += " ";
1510 aDump += TCollection_AsciiString(nbo) + "\n";
1517 //=============================================================================
1519 * CheckSelfIntersections
1521 //=============================================================================
1522 bool GEOMImpl_IMeasureOperations::CheckSelfIntersections
1523 (Handle(GEOM_Object) theShape,
1524 Handle(TColStd_HSequenceOfInteger)& theIntersections)
1527 bool isGood = false;
1529 if (theIntersections.IsNull())
1530 theIntersections = new TColStd_HSequenceOfInteger;
1532 theIntersections->Clear();
1534 if (theShape.IsNull())
1537 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1538 if (aRefShape.IsNull()) return isGood;
1540 TopoDS_Shape aShape = aRefShape->GetValue();
1541 if (aShape.IsNull()) return isGood;
1544 TopoDS_Shape aScopy;
1546 GEOMAlgo_AlgoTools::CopyShape(aShape, aScopy);
1548 // Map sub-shapes and their indices
1549 TopTools_IndexedMapOfShape anIndices;
1550 TopExp::MapShapes(aScopy, anIndices);
1552 BOPCol_ListOfShape aLCS;
1553 aLCS.Append(aScopy);
1555 BOPAlgo_CheckerSI aCSI; // checker of self-interferences
1556 aCSI.SetArguments(aLCS);
1558 // 1. Launch the checker
1560 Standard_Integer iErr = aCSI.ErrorStatus();
1564 Standard_Integer aNbS, n1, n2;
1565 BOPDS_MapIteratorMapOfPassKey aItMPK;
1567 // 2. Take the shapes from DS
1568 const BOPDS_DS& aDS = aCSI.DS();
1569 aNbS=aDS.NbShapes();
1571 // 3. Get the pairs of interfered shapes
1572 const BOPDS_MapOfPassKey& aMPK=aDS.Interferences();
1573 aItMPK.Initialize(aMPK);
1574 for (; aItMPK.More(); aItMPK.Next()) {
1575 const BOPDS_PassKey& aPK=aItMPK.Value();
1578 if (n1 > aNbS || n2 > aNbS){
1579 return false; // Error
1581 const TopoDS_Shape& aS1 = aDS.Shape(n1);
1582 const TopoDS_Shape& aS2 = aDS.Shape(n2);
1584 theIntersections->Append(anIndices.FindIndex(aS1));
1585 theIntersections->Append(anIndices.FindIndex(aS2));
1596 //=============================================================================
1600 //=============================================================================
1601 TCollection_AsciiString GEOMImpl_IMeasureOperations::IsGoodForSolid (Handle(GEOM_Object) theShape)
1605 TCollection_AsciiString aRes = "";
1607 if (theShape.IsNull()) {
1608 aRes = "WRN_NULL_OBJECT_OR_SHAPE";
1611 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1612 if (aRefShape.IsNull()) {
1613 aRes = "WRN_NULL_OBJECT_OR_SHAPE";
1616 TopoDS_Shape aShape = aRefShape->GetValue();
1617 if (aShape.IsNull()) {
1618 aRes = "WRN_NULL_OBJECT_OR_SHAPE";
1621 if (aShape.ShapeType() == TopAbs_COMPOUND) {
1622 TopoDS_Iterator It (aShape, Standard_True, Standard_True);
1623 if (It.More()) aShape = It.Value();
1625 if (aShape.ShapeType() == TopAbs_SHELL) {
1626 BRepCheck_Shell chkShell (TopoDS::Shell(aShape));
1627 if (chkShell.Closed() == BRepCheck_NotClosed) {
1628 aRes = "WRN_SHAPE_UNCLOSED";
1632 aRes = "WRN_SHAPE_NOT_SHELL";
1644 //=============================================================================
1648 //=============================================================================
1649 TCollection_AsciiString GEOMImpl_IMeasureOperations::WhatIs (Handle(GEOM_Object) theShape)
1653 TCollection_AsciiString Astr;
1655 if (theShape.IsNull()) return Astr;
1657 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1658 if (aRefShape.IsNull()) return Astr;
1660 TopoDS_Shape aShape = aRefShape->GetValue();
1661 if (aShape.IsNull()) {
1662 SetErrorCode("The Objects has NULL Shape");
1666 //Compute the parameters
1667 if (aShape.ShapeType() == TopAbs_EDGE) {
1668 if (BRep_Tool::Degenerated(TopoDS::Edge(aShape))) {
1669 Astr = Astr + " It is a degenerated edge \n";
1673 Astr = Astr + " Number of sub-shapes : \n";
1676 #if OCC_VERSION_LARGE > 0x06010000
1679 int iType, nbTypes [TopAbs_SHAPE];
1680 for (iType = 0; iType < TopAbs_SHAPE; ++iType)
1682 nbTypes[aShape.ShapeType()]++;
1684 TopTools_MapOfShape aMapOfShape;
1685 aMapOfShape.Add(aShape);
1686 TopTools_ListOfShape aListOfShape;
1687 aListOfShape.Append(aShape);
1689 TopTools_ListIteratorOfListOfShape itL (aListOfShape);
1690 for (; itL.More(); itL.Next()) {
1691 TopoDS_Iterator it (itL.Value());
1692 for (; it.More(); it.Next()) {
1693 TopoDS_Shape s = it.Value();
1694 if (aMapOfShape.Add(s)) {
1695 aListOfShape.Append(s);
1696 nbTypes[s.ShapeType()]++;
1701 Astr = Astr + " VERTEX : " + TCollection_AsciiString(nbTypes[TopAbs_VERTEX]) + "\n";
1702 Astr = Astr + " EDGE : " + TCollection_AsciiString(nbTypes[TopAbs_EDGE]) + "\n";
1703 Astr = Astr + " WIRE : " + TCollection_AsciiString(nbTypes[TopAbs_WIRE]) + "\n";
1704 Astr = Astr + " FACE : " + TCollection_AsciiString(nbTypes[TopAbs_FACE]) + "\n";
1705 Astr = Astr + " SHELL : " + TCollection_AsciiString(nbTypes[TopAbs_SHELL]) + "\n";
1706 Astr = Astr + " SOLID : " + TCollection_AsciiString(nbTypes[TopAbs_SOLID]) + "\n";
1707 Astr = Astr + " COMPSOLID : " + TCollection_AsciiString(nbTypes[TopAbs_COMPSOLID]) + "\n";
1708 Astr = Astr + " COMPOUND : " + TCollection_AsciiString(nbTypes[TopAbs_COMPOUND]) + "\n";
1709 Astr = Astr + " SHAPE : " + TCollection_AsciiString(aMapOfShape.Extent());
1711 catch (Standard_Failure) {
1712 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1713 SetErrorCode(aFail->GetMessageString());
1721 //=============================================================================
1725 //=============================================================================
1727 GEOMImpl_IMeasureOperations::AreCoordsInside(Handle(GEOM_Object) theShape,
1728 const std::vector<double>& coords,
1731 std::vector<bool> isInsideRes;
1732 if (!theShape.IsNull()) {
1733 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1734 if (!aRefShape.IsNull()) {
1735 TopoDS_Shape aShape = aRefShape->GetValue();
1736 if (!aShape.IsNull())
1738 TopTools_IndexedMapOfShape mapShape;
1740 TopExp_Explorer anExp;
1741 for ( anExp.Init( aShape, TopAbs_SOLID ); anExp.More(); anExp.Next() )
1742 mapShape.Add( anExp.Current() );
1743 for ( anExp.Init( aShape, TopAbs_FACE, TopAbs_SOLID ); anExp.More(); anExp.Next() )
1744 mapShape.Add( anExp.Current() );
1745 for ( anExp.Init( aShape, TopAbs_EDGE, TopAbs_FACE ); anExp.More(); anExp.Next() )
1746 mapShape.Add( anExp.Current() );
1747 for ( anExp.Init( aShape, TopAbs_VERTEX, TopAbs_EDGE ); anExp.More(); anExp.Next() )
1748 mapShape.Add( anExp.Current() ); //// ?????????
1750 size_t nb_points = coords.size()/3, nb_points_inside = 0;
1751 isInsideRes.resize( nb_points, false );
1753 for ( int iS = 1; iS <= mapShape.Extent(); ++iS )
1755 if ( nb_points_inside == nb_points )
1757 aShape = mapShape( iS );
1758 switch ( aShape.ShapeType() ) {
1761 BRepClass3d_SolidClassifier SC( TopoDS::Solid( aShape ));
1762 for ( size_t i = 0; i < nb_points; i++)
1764 if ( isInsideRes[ i ]) continue;
1765 gp_Pnt aPnt( coords[3*i], coords[3*i+1], coords[3*i+2] );
1766 SC.Perform( aPnt, tolerance );
1767 isInsideRes[ i ] = (( SC.State() == TopAbs_IN ) || ( SC.State() == TopAbs_ON ));
1768 nb_points_inside += isInsideRes[ i ];
1774 Standard_Real u1,u2,v1,v2;
1775 const TopoDS_Face& face = TopoDS::Face( aShape );
1776 Handle(Geom_Surface) surf = BRep_Tool::Surface( face );
1777 surf->Bounds( u1,u2,v1,v2 );
1778 GeomAPI_ProjectPointOnSurf project;
1779 project.Init(surf, u1,u2, v1,v2, tolerance );
1780 for ( size_t i = 0; i < nb_points; i++)
1782 if ( isInsideRes[ i ]) continue;
1783 gp_Pnt aPnt( coords[3*i], coords[3*i+1], coords[3*i+2] );
1784 project.Perform( aPnt );
1785 if ( project.IsDone() &&
1786 project.NbPoints() > 0 &&
1787 project.LowerDistance() <= tolerance )
1789 Quantity_Parameter u, v;
1790 project.LowerDistanceParameters(u, v);
1791 gp_Pnt2d uv( u, v );
1792 BRepClass_FaceClassifier FC ( face, uv, tolerance );
1793 isInsideRes[ i ] = (( FC.State() == TopAbs_IN ) || ( FC.State() == TopAbs_ON ));
1794 nb_points_inside += isInsideRes[ i ];
1802 const TopoDS_Edge& edge = TopoDS::Edge( aShape );
1803 Handle(Geom_Curve) curve = BRep_Tool::Curve( edge, f, l );
1804 GeomAPI_ProjectPointOnCurve project;
1805 project.Init( curve, f, l );
1806 for ( size_t i = 0; i < nb_points; i++)
1808 if ( isInsideRes[ i ]) continue;
1809 gp_Pnt aPnt( coords[3*i], coords[3*i+1], coords[3*i+2] );
1810 project.Perform( aPnt );
1811 isInsideRes[ i ] = ( project.NbPoints() > 0 &&
1812 project.LowerDistance() <= tolerance );
1813 nb_points_inside += isInsideRes[ i ];
1819 gp_Pnt aVPnt = BRep_Tool::Pnt( TopoDS::Vertex( aShape ));
1820 for ( size_t i = 0; i < nb_points; i++)
1822 if ( isInsideRes[ i ]) continue;
1823 gp_Pnt aPnt( coords[3*i], coords[3*i+1], coords[3*i+2] );
1824 isInsideRes[ i ] = ( aPnt.SquareDistance( aVPnt ) <= tolerance * tolerance );
1825 nb_points_inside += isInsideRes[ i ];
1830 } // switch ( aShape.ShapeType() )
1838 //=============================================================================
1842 //=============================================================================
1844 GEOMImpl_IMeasureOperations::GetMinDistance (Handle(GEOM_Object) theShape1,
1845 Handle(GEOM_Object) theShape2,
1854 Standard_Real MinDist = 1.e9;
1856 if (theShape1.IsNull() || theShape2.IsNull()) return MinDist;
1858 Handle(GEOM_Function) aRefShape1 = theShape1->GetLastFunction();
1859 Handle(GEOM_Function) aRefShape2 = theShape2->GetLastFunction();
1860 if (aRefShape1.IsNull() || aRefShape2.IsNull()) return MinDist;
1862 TopoDS_Shape aShape1 = aRefShape1->GetValue();
1863 TopoDS_Shape aShape2 = aRefShape2->GetValue();
1864 if (aShape1.IsNull() || aShape2.IsNull()) {
1865 SetErrorCode("One of Objects has NULL Shape");
1869 //Compute the parameters
1871 #if OCC_VERSION_LARGE > 0x06010000
1875 gp_Pnt aPnt1, aPnt2;
1877 MinDist = GEOMUtils::GetMinDistance(aShape1, aShape2, aPnt1, aPnt2);
1879 if (MinDist >= 0.0) {
1880 aPnt1.Coord(X1, Y1, Z1);
1881 aPnt2.Coord(X2, Y2, Z2);
1886 catch (Standard_Failure) {
1887 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1888 SetErrorCode(aFail->GetMessageString());
1896 //=======================================================================
1898 * Get coordinates of closest points of two shapes
1900 //=======================================================================
1901 Standard_Integer GEOMImpl_IMeasureOperations::ClosestPoints (Handle(GEOM_Object) theShape1,
1902 Handle(GEOM_Object) theShape2,
1903 Handle(TColStd_HSequenceOfReal)& theDoubles)
1906 Standard_Integer nbSolutions = 0;
1908 if (theShape1.IsNull() || theShape2.IsNull()) return nbSolutions;
1910 Handle(GEOM_Function) aRefShape1 = theShape1->GetLastFunction();
1911 Handle(GEOM_Function) aRefShape2 = theShape2->GetLastFunction();
1912 if (aRefShape1.IsNull() || aRefShape2.IsNull()) return nbSolutions;
1914 TopoDS_Shape aShape1 = aRefShape1->GetValue();
1915 TopoDS_Shape aShape2 = aRefShape2->GetValue();
1916 if (aShape1.IsNull() || aShape2.IsNull()) {
1917 SetErrorCode("One of Objects has NULL Shape");
1921 // Compute the extremities
1923 #if OCC_VERSION_LARGE > 0x06010000
1928 // additional workaround for bugs 19899, 19908 and 19910 from Mantis
1930 double dist = GEOMUtils::GetMinDistanceSingular(aShape1, aShape2, P1, P2);
1934 theDoubles->Append(P1.X());
1935 theDoubles->Append(P1.Y());
1936 theDoubles->Append(P1.Z());
1937 theDoubles->Append(P2.X());
1938 theDoubles->Append(P2.Y());
1939 theDoubles->Append(P2.Z());
1945 BRepExtrema_DistShapeShape dst (aShape1, aShape2);
1947 nbSolutions = dst.NbSolution();
1948 if (theDoubles.IsNull()) theDoubles = new TColStd_HSequenceOfReal;
1951 for (int i = 1; i <= nbSolutions; i++) {
1952 P1 = dst.PointOnShape1(i);
1953 P2 = dst.PointOnShape2(i);
1955 theDoubles->Append(P1.X());
1956 theDoubles->Append(P1.Y());
1957 theDoubles->Append(P1.Z());
1958 theDoubles->Append(P2.X());
1959 theDoubles->Append(P2.Y());
1960 theDoubles->Append(P2.Z());
1964 catch (Standard_Failure) {
1965 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1966 SetErrorCode(aFail->GetMessageString());
1974 //=======================================================================
1976 * Get coordinates of point
1978 //=======================================================================
1979 void GEOMImpl_IMeasureOperations::PointCoordinates (Handle(GEOM_Object) theShape,
1980 Standard_Real& theX, Standard_Real& theY, Standard_Real& theZ)
1984 if (theShape.IsNull())
1987 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1988 if (aRefShape.IsNull())
1991 TopoDS_Shape aShape = aRefShape->GetValue();
1992 if (aShape.IsNull() || aShape.ShapeType() != TopAbs_VERTEX)
1994 SetErrorCode( "Shape must be a vertex" );
1999 #if OCC_VERSION_LARGE > 0x06010000
2002 gp_Pnt aPnt = BRep_Tool::Pnt( TopoDS::Vertex( aShape ) );
2009 catch (Standard_Failure)
2011 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
2012 SetErrorCode( aFail->GetMessageString() );
2016 //=======================================================================
2018 * Compute angle (in degrees) between two lines
2020 //=======================================================================
2021 Standard_Real GEOMImpl_IMeasureOperations::GetAngle (Handle(GEOM_Object) theLine1,
2022 Handle(GEOM_Object) theLine2)
2024 if (theLine1->GetType() == GEOM_VECTOR &&
2025 theLine2->GetType() == GEOM_VECTOR)
2026 return GetAngleBtwVectors(theLine1, theLine2);
2030 Standard_Real anAngle = -1.0;
2032 if (theLine1.IsNull() || theLine2.IsNull())
2035 Handle(GEOM_Function) aRefLine1 = theLine1->GetLastFunction();
2036 Handle(GEOM_Function) aRefLine2 = theLine2->GetLastFunction();
2037 if (aRefLine1.IsNull() || aRefLine2.IsNull())
2040 TopoDS_Shape aLine1 = aRefLine1->GetValue();
2041 TopoDS_Shape aLine2 = aRefLine2->GetValue();
2042 if (aLine1.IsNull() || aLine2.IsNull() ||
2043 aLine1.ShapeType() != TopAbs_EDGE ||
2044 aLine2.ShapeType() != TopAbs_EDGE)
2046 SetErrorCode("Two edges must be given");
2051 #if OCC_VERSION_LARGE > 0x06010000
2054 TopoDS_Edge E1 = TopoDS::Edge(aLine1);
2055 TopoDS_Edge E2 = TopoDS::Edge(aLine2);
2058 Handle(Geom_Curve) C1 = BRep_Tool::Curve(E1,fp,lp);
2059 Handle(Geom_Curve) C2 = BRep_Tool::Curve(E2,fp,lp);
2061 if ( C1.IsNull() || C2.IsNull() ||
2062 !C1->IsKind(STANDARD_TYPE(Geom_Line)) ||
2063 !C2->IsKind(STANDARD_TYPE(Geom_Line)))
2065 SetErrorCode("The edges must be linear");
2069 Handle(Geom_Line) L1 = Handle(Geom_Line)::DownCast(C1);
2070 Handle(Geom_Line) L2 = Handle(Geom_Line)::DownCast(C2);
2072 gp_Lin aLin1 = L1->Lin();
2073 gp_Lin aLin2 = L2->Lin();
2075 anAngle = aLin1.Angle(aLin2);
2076 anAngle *= 180. / M_PI; // convert radians into degrees
2078 if (anAngle > 90.0) {
2079 anAngle = 180.0 - anAngle;
2084 catch (Standard_Failure)
2086 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
2087 SetErrorCode(aFail->GetMessageString());
2093 //=======================================================================
2095 * Compute angle (in degrees) between two vectors
2097 //=======================================================================
2098 Standard_Real GEOMImpl_IMeasureOperations::GetAngleBtwVectors (Handle(GEOM_Object) theVec1,
2099 Handle(GEOM_Object) theVec2)
2103 Standard_Real anAngle = -1.0;
2105 if (theVec1.IsNull() || theVec2.IsNull())
2108 if (theVec1->GetType() != GEOM_VECTOR || theVec2->GetType() != GEOM_VECTOR) {
2109 SetErrorCode("Two vectors must be given");
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);
2470 //=======================================================================
2471 //function : FillErrorsSub
2472 //purpose : Fill the errors list of subshapes on shape.
2473 //=======================================================================
2474 void GEOMImpl_IMeasureOperations::FillErrorsSub
2475 (const BRepCheck_Analyzer &theAna,
2476 const TopoDS_Shape &theShape,
2477 const TopAbs_ShapeEnum theSubType,
2478 TopTools_DataMapOfIntegerListOfShape &theMapErrors) const
2480 TopExp_Explorer anExp(theShape, theSubType);
2481 TopTools_MapOfShape aMapSubShapes;
2483 for (; anExp.More(); anExp.Next()) {
2484 const TopoDS_Shape &aSubShape = anExp.Current();
2486 if (aMapSubShapes.Add(aSubShape)) {
2487 const Handle(BRepCheck_Result) &aRes = theAna.Result(aSubShape);
2489 for (aRes->InitContextIterator();
2490 aRes->MoreShapeInContext();
2491 aRes->NextShapeInContext()) {
2492 if (aRes->ContextualShape().IsSame(theShape)) {
2493 BRepCheck_ListIteratorOfListOfStatus itl(aRes->StatusOnShape());
2495 if (itl.Value() != BRepCheck_NoError) {
2496 // Add all errors for theShape and its sub-shape.
2497 for (;itl.More(); itl.Next()) {
2498 const Standard_Integer aStat = (Standard_Integer)itl.Value();
2500 if (!theMapErrors.IsBound(aStat)) {
2501 TopTools_ListOfShape anEmpty;
2503 theMapErrors.Bind(aStat, anEmpty);
2506 TopTools_ListOfShape &theShapes = theMapErrors.ChangeFind(aStat);
2508 theShapes.Append(aSubShape);
2509 theShapes.Append(theShape);
2520 //=======================================================================
2521 //function : FillErrors
2522 //purpose : Fill the errors list.
2523 //=======================================================================
2524 void GEOMImpl_IMeasureOperations::FillErrors
2525 (const BRepCheck_Analyzer &theAna,
2526 const TopoDS_Shape &theShape,
2527 TopTools_DataMapOfIntegerListOfShape &theMapErrors,
2528 TopTools_MapOfShape &theMapShapes) const
2530 if (theMapShapes.Add(theShape)) {
2531 // Fill errors of child shapes.
2532 for (TopoDS_Iterator iter(theShape); iter.More(); iter.Next()) {
2533 FillErrors(theAna, iter.Value(), theMapErrors, theMapShapes);
2536 // Fill errors of theShape.
2537 const Handle(BRepCheck_Result) &aRes = theAna.Result(theShape);
2539 if (!aRes.IsNull()) {
2540 BRepCheck_ListIteratorOfListOfStatus itl(aRes->Status());
2542 if (itl.Value() != BRepCheck_NoError) {
2543 // Add all errors for theShape.
2544 for (;itl.More(); itl.Next()) {
2545 const Standard_Integer aStat = (Standard_Integer)itl.Value();
2547 if (!theMapErrors.IsBound(aStat)) {
2548 TopTools_ListOfShape anEmpty;
2550 theMapErrors.Bind(aStat, anEmpty);
2553 theMapErrors.ChangeFind(aStat).Append(theShape);
2558 // Add errors of subshapes on theShape.
2559 const TopAbs_ShapeEnum aType = theShape.ShapeType();
2563 FillErrorsSub(theAna, theShape, TopAbs_VERTEX, theMapErrors);
2566 FillErrorsSub(theAna, theShape, TopAbs_WIRE, theMapErrors);
2567 FillErrorsSub(theAna, theShape, TopAbs_EDGE, theMapErrors);
2568 FillErrorsSub(theAna, theShape, TopAbs_VERTEX, theMapErrors);
2571 FillErrorsSub(theAna, theShape, TopAbs_SHELL, theMapErrors);
2579 //=======================================================================
2580 //function : FillErrors
2581 //purpose : Fill the errors list.
2582 //=======================================================================
2583 void GEOMImpl_IMeasureOperations::FillErrors
2584 (const BRepCheck_Analyzer &theAna,
2585 const TopoDS_Shape &theShape,
2586 std::list<ShapeError> &theErrors) const
2588 // Fill the errors map.
2589 TopTools_DataMapOfIntegerListOfShape aMapErrors;
2590 TopTools_MapOfShape aMapShapes;
2592 FillErrors(theAna, theShape, aMapErrors, aMapShapes);
2594 // Map sub-shapes and their indices
2595 TopTools_IndexedMapOfShape anIndices;
2597 TopExp::MapShapes(theShape, anIndices);
2599 TopTools_DataMapIteratorOfDataMapOfIntegerListOfShape aMapIter(aMapErrors);
2601 for (; aMapIter.More(); aMapIter.Next()) {
2604 anError.error = (BRepCheck_Status)aMapIter.Key();
2606 TopTools_ListIteratorOfListOfShape aListIter(aMapIter.Value());
2607 TopTools_MapOfShape aMapUnique;
2609 for (; aListIter.More(); aListIter.Next()) {
2610 const TopoDS_Shape &aShape = aListIter.Value();
2612 if (aMapUnique.Add(aShape)) {
2613 const Standard_Integer anIndex = anIndices.FindIndex(aShape);
2615 anError.incriminated.push_back(anIndex);
2619 if (!anError.incriminated.empty()) {
2620 theErrors.push_back(anError);