1 // Copyright (C) 2007-2008 CEA/DEN, EDF R&D, OPEN CASCADE
3 // Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
6 // This library is free software; you can redistribute it and/or
7 // modify it under the terms of the GNU Lesser General Public
8 // License as published by the Free Software Foundation; either
9 // version 2.1 of the License.
11 // This library is distributed in the hope that it will be useful,
12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 // Lesser General Public License for more details.
16 // You should have received a copy of the GNU Lesser General Public
17 // License along with this library; if not, write to the Free Software
18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
22 #include <Standard_Stream.hxx>
24 #include <GEOMImpl_IMeasureOperations.hxx>
26 #include <GEOMImpl_Types.hxx>
27 #include <GEOMImpl_MeasureDriver.hxx>
28 #include <GEOMImpl_IMeasure.hxx>
30 #include <GEOMAlgo_ShapeInfo.hxx>
31 #include <GEOMAlgo_ShapeInfoFiller.hxx>
33 #include <GEOM_Function.hxx>
34 #include <GEOM_PythonDump.hxx>
36 #include <utilities.h>
38 #include <Utils_ExceptHandlers.hxx>
41 #include <TFunction_DriverTable.hxx>
42 #include <TFunction_Driver.hxx>
43 #include <TFunction_Logbook.hxx>
44 #include <TDF_Tool.hxx>
46 #include <BRep_Tool.hxx>
47 #include <BRepAdaptor_Surface.hxx>
48 #include <BRepBndLib.hxx>
49 #include <BRepCheck.hxx>
50 #include <BRepCheck_Result.hxx>
51 #include <BRepCheck_ListIteratorOfListOfStatus.hxx>
52 #include <BRepExtrema_DistShapeShape.hxx>
53 #include <BRepGProp.hxx>
54 #include <BRepTools.hxx>
56 #include <Bnd_Box.hxx>
58 #include <GProp_GProps.hxx>
59 #include <GProp_PrincipalProps.hxx>
63 #include <TopoDS_Edge.hxx>
64 #include <TopoDS_Face.hxx>
65 #include <TopoDS_Shape.hxx>
66 #include <TopoDS_Vertex.hxx>
67 #include <TopoDS_Iterator.hxx>
68 #include <TopExp_Explorer.hxx>
69 #include <TopTools_MapOfShape.hxx>
70 #include <TopTools_ListOfShape.hxx>
71 #include <TopTools_ListIteratorOfListOfShape.hxx>
73 #include <GeomAbs_SurfaceType.hxx>
74 #include <Geom_Surface.hxx>
75 #include <Geom_Plane.hxx>
76 #include <Geom_SphericalSurface.hxx>
77 #include <Geom_CylindricalSurface.hxx>
78 #include <Geom_ToroidalSurface.hxx>
79 #include <Geom_ConicalSurface.hxx>
80 #include <Geom_SurfaceOfLinearExtrusion.hxx>
81 #include <Geom_SurfaceOfRevolution.hxx>
82 #include <Geom_BezierSurface.hxx>
83 #include <Geom_BSplineSurface.hxx>
84 #include <Geom_RectangularTrimmedSurface.hxx>
85 #include <Geom_OffsetSurface.hxx>
86 #include <Geom_Line.hxx>
91 #include <GeomAPI_ProjectPointOnCurve.hxx>
92 #include <GeomLProp_CLProps.hxx>
93 #include <GeomLProp_SLProps.hxx>
94 #include <ShapeAnalysis.hxx>
95 #include <ShapeAnalysis_Surface.hxx>
97 #include <Standard_Failure.hxx>
98 #include <Standard_ErrorHandler.hxx> // CAREFUL ! position of this file is critic : see Lucien PIGNOLONI / OCC
100 #include <BRepClass3d_SolidClassifier.hxx>
101 #include <BRep_Builder.hxx>
102 #include <GeomAPI_IntSS.hxx>
103 #include <Geom_Circle.hxx>
104 #include <Geom_SphericalSurface.hxx>
105 #include <Geom_ToroidalSurface.hxx>
106 #include <ShapeFix_Shape.hxx>
107 #include <TopoDS_Compound.hxx>
110 //=============================================================================
114 //=============================================================================
115 GEOMImpl_IMeasureOperations::GEOMImpl_IMeasureOperations (GEOM_Engine* theEngine, int theDocID)
116 : GEOM_IOperations(theEngine, theDocID)
118 MESSAGE("GEOMImpl_IMeasureOperations::GEOMImpl_IMeasureOperations");
121 //=============================================================================
125 //=============================================================================
126 GEOMImpl_IMeasureOperations::~GEOMImpl_IMeasureOperations()
128 MESSAGE("GEOMImpl_IMeasureOperations::~GEOMImpl_IMeasureOperations");
131 //=============================================================================
132 /*! Get kind and parameters of the given shape.
134 //=============================================================================
135 GEOMImpl_IMeasureOperations::ShapeKind GEOMImpl_IMeasureOperations::KindOfShape
136 (Handle(GEOM_Object) theShape,
137 Handle(TColStd_HSequenceOfInteger)& theIntegers,
138 Handle(TColStd_HSequenceOfReal)& theDoubles)
141 ShapeKind aKind = SK_NO_SHAPE;
143 if (theIntegers.IsNull()) theIntegers = new TColStd_HSequenceOfInteger;
144 else theIntegers->Clear();
146 if (theDoubles.IsNull()) theDoubles = new TColStd_HSequenceOfReal;
147 else theDoubles->Clear();
149 if (theShape.IsNull())
152 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
153 if (aRefShape.IsNull()) return aKind;
155 TopoDS_Shape aShape = aRefShape->GetValue();
156 if (aShape.IsNull()) return aKind;
159 GEOMAlgo_ShapeInfoFiller aSF;
160 aSF.SetShape(aShape);
162 Standard_Integer iErr = aSF.ErrorStatus();
164 SetErrorCode("Error in GEOMAlgo_ShapeInfoFiller");
167 const GEOMAlgo_ShapeInfo& anInfo = aSF.Info();
169 // Interprete results
170 TopAbs_ShapeEnum aType = anInfo.Type();
173 case TopAbs_COMPOUND:
174 case TopAbs_COMPSOLID:
176 // (+) geompy.kind.COMPOUND nb_solids nb_faces nb_edges nb_vertices
177 // (+) geompy.kind.COMPSOLID nb_solids nb_faces nb_edges nb_vertices
178 // ??? "nb_faces" - all faces or only 'standalone' faces?
179 if (aType == TopAbs_COMPOUND)
182 aKind = SK_COMPSOLID;
184 //theIntegers->Append(anInfo.NbSubShapes(TopAbs_COMPOUND));
185 //theIntegers->Append(anInfo.NbSubShapes(TopAbs_COMPSOLID));
186 theIntegers->Append(anInfo.NbSubShapes(TopAbs_SOLID));
187 theIntegers->Append(anInfo.NbSubShapes(TopAbs_FACE));
188 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
189 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
195 // (+) geompy.kind.SHELL geompy.info.closed nb_faces nb_edges nb_vertices
196 // (+) geompy.kind.SHELL geompy.info.unclosed nb_faces nb_edges nb_vertices
199 theIntegers->Append((int)anInfo.KindOfClosed());
201 theIntegers->Append(anInfo.NbSubShapes(TopAbs_FACE));
202 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
203 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
209 // (+) geompy.kind.WIRE geompy.info.closed nb_edges nb_vertices
210 // (+) geompy.kind.WIRE geompy.info.unclosed nb_edges nb_vertices
213 theIntegers->Append((int)anInfo.KindOfClosed());
215 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
216 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
224 GEOMAlgo_KindOfName aKN = anInfo.KindOfName();
227 case GEOMAlgo_KN_SPHERE:
228 // (+) geompy.kind.SPHERE xc yc zc R
232 gp_Pnt aC = anInfo.Location();
233 theDoubles->Append(aC.X());
234 theDoubles->Append(aC.Y());
235 theDoubles->Append(aC.Z());
237 theDoubles->Append(anInfo.Radius1());
240 case GEOMAlgo_KN_CYLINDER:
241 // (+) geompy.kind.CYLINDER xb yb zb dx dy dz R H
245 gp_Pnt aC = anInfo.Location();
246 theDoubles->Append(aC.X());
247 theDoubles->Append(aC.Y());
248 theDoubles->Append(aC.Z());
250 gp_Ax3 anAx3 = anInfo.Position();
251 gp_Dir aD = anAx3.Direction();
252 theDoubles->Append(aD.X());
253 theDoubles->Append(aD.Y());
254 theDoubles->Append(aD.Z());
256 theDoubles->Append(anInfo.Radius1());
257 theDoubles->Append(anInfo.Height());
260 case GEOMAlgo_KN_BOX:
261 // (+) geompy.kind.BOX xc yc zc ax ay az
265 gp_Pnt aC = anInfo.Location();
266 theDoubles->Append(aC.X());
267 theDoubles->Append(aC.Y());
268 theDoubles->Append(aC.Z());
270 gp_Ax3 anAx3 = anInfo.Position();
271 gp_Dir aD = anAx3.Direction();
272 gp_Dir aX = anAx3.XDirection();
275 if (aD.IsParallel(gp::DZ(), Precision::Angular()) &&
276 aX.IsParallel(gp::DX(), Precision::Angular())) {
277 theDoubles->Append(anInfo.Length()); // ax'
278 theDoubles->Append(anInfo.Width()); // ay'
279 theDoubles->Append(anInfo.Height()); // az'
281 else if (aD.IsParallel(gp::DZ(), Precision::Angular()) &&
282 aX.IsParallel(gp::DY(), Precision::Angular())) {
283 theDoubles->Append(anInfo.Width()); // ay'
284 theDoubles->Append(anInfo.Length()); // ax'
285 theDoubles->Append(anInfo.Height()); // az'
287 else if (aD.IsParallel(gp::DX(), Precision::Angular()) &&
288 aX.IsParallel(gp::DZ(), Precision::Angular())) {
289 theDoubles->Append(anInfo.Height()); // az'
290 theDoubles->Append(anInfo.Width()); // ay'
291 theDoubles->Append(anInfo.Length()); // ax'
293 else if (aD.IsParallel(gp::DX(), Precision::Angular()) &&
294 aX.IsParallel(gp::DY(), Precision::Angular())) {
295 theDoubles->Append(anInfo.Height()); // az'
296 theDoubles->Append(anInfo.Length()); // ax'
297 theDoubles->Append(anInfo.Width()); // ay'
299 else if (aD.IsParallel(gp::DY(), Precision::Angular()) &&
300 aX.IsParallel(gp::DZ(), Precision::Angular())) {
301 theDoubles->Append(anInfo.Width()); // ay'
302 theDoubles->Append(anInfo.Height()); // az'
303 theDoubles->Append(anInfo.Length()); // ax'
305 else if (aD.IsParallel(gp::DY(), Precision::Angular()) &&
306 aX.IsParallel(gp::DX(), Precision::Angular())) {
307 theDoubles->Append(anInfo.Length()); // ax'
308 theDoubles->Append(anInfo.Height()); // az'
309 theDoubles->Append(anInfo.Width()); // ay'
312 // (+) geompy.kind.ROTATED_BOX xo yo zo zx zy zz xx xy xz ax ay az
313 aKind = SK_ROTATED_BOX;
315 // Direction and XDirection
316 theDoubles->Append(aD.X());
317 theDoubles->Append(aD.Y());
318 theDoubles->Append(aD.Z());
320 theDoubles->Append(aX.X());
321 theDoubles->Append(aX.Y());
322 theDoubles->Append(aX.Z());
325 theDoubles->Append(anInfo.Length());
326 theDoubles->Append(anInfo.Width());
327 theDoubles->Append(anInfo.Height());
331 case GEOMAlgo_KN_TORUS:
332 // (+) geompy.kind.TORUS xc yc zc dx dy dz R_1 R_2
336 gp_Pnt aO = anInfo.Location();
337 theDoubles->Append(aO.X());
338 theDoubles->Append(aO.Y());
339 theDoubles->Append(aO.Z());
341 gp_Ax3 anAx3 = anInfo.Position();
342 gp_Dir aD = anAx3.Direction();
343 theDoubles->Append(aD.X());
344 theDoubles->Append(aD.Y());
345 theDoubles->Append(aD.Z());
347 theDoubles->Append(anInfo.Radius1());
348 theDoubles->Append(anInfo.Radius2());
351 case GEOMAlgo_KN_CONE:
352 // (+) geompy.kind.CONE xb yb zb dx dy dz R_1 R_2 H
356 gp_Pnt aO = anInfo.Location();
357 theDoubles->Append(aO.X());
358 theDoubles->Append(aO.Y());
359 theDoubles->Append(aO.Z());
361 gp_Ax3 anAx3 = anInfo.Position();
362 gp_Dir aD = anAx3.Direction();
363 theDoubles->Append(aD.X());
364 theDoubles->Append(aD.Y());
365 theDoubles->Append(aD.Z());
367 theDoubles->Append(anInfo.Radius1());
368 theDoubles->Append(anInfo.Radius2());
369 theDoubles->Append(anInfo.Height());
372 case GEOMAlgo_KN_POLYHEDRON:
373 // (+) geompy.kind.POLYHEDRON nb_faces nb_edges nb_vertices
375 aKind = SK_POLYHEDRON;
377 theIntegers->Append(anInfo.NbSubShapes(TopAbs_FACE));
378 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
379 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
383 // (+) geompy.kind.SOLID nb_faces nb_edges nb_vertices
385 theIntegers->Append(anInfo.NbSubShapes(TopAbs_FACE));
386 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
387 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
397 GEOMAlgo_KindOfName aKN = anInfo.KindOfName();
399 case GEOMAlgo_KN_SPHERE:
400 // (+) geompy.kind.SPHERE2D xc yc zc R
404 gp_Pnt aC = anInfo.Location();
405 theDoubles->Append(aC.X());
406 theDoubles->Append(aC.Y());
407 theDoubles->Append(aC.Z());
409 theDoubles->Append(anInfo.Radius1());
412 case GEOMAlgo_KN_CYLINDER:
413 // (+) geompy.kind.CYLINDER2D xb yb zb dx dy dz R H
415 aKind = SK_CYLINDER2D;
417 gp_Pnt aO = anInfo.Location();
418 theDoubles->Append(aO.X());
419 theDoubles->Append(aO.Y());
420 theDoubles->Append(aO.Z());
422 gp_Ax3 anAx3 = anInfo.Position();
423 gp_Dir aD = anAx3.Direction();
424 theDoubles->Append(aD.X());
425 theDoubles->Append(aD.Y());
426 theDoubles->Append(aD.Z());
428 theDoubles->Append(anInfo.Radius1());
429 theDoubles->Append(anInfo.Height());
432 case GEOMAlgo_KN_TORUS:
433 // (+) geompy.kind.TORUS2D xc yc zc dx dy dz R_1 R_2
437 gp_Pnt aO = anInfo.Location();
438 theDoubles->Append(aO.X());
439 theDoubles->Append(aO.Y());
440 theDoubles->Append(aO.Z());
442 gp_Ax3 anAx3 = anInfo.Position();
443 gp_Dir aD = anAx3.Direction();
444 theDoubles->Append(aD.X());
445 theDoubles->Append(aD.Y());
446 theDoubles->Append(aD.Z());
448 theDoubles->Append(anInfo.Radius1());
449 theDoubles->Append(anInfo.Radius2());
452 case GEOMAlgo_KN_CONE:
453 // (+) geompy.kind.CONE2D xc yc zc dx dy dz R_1 R_2 H
457 gp_Pnt aO = anInfo.Location();
458 theDoubles->Append(aO.X());
459 theDoubles->Append(aO.Y());
460 theDoubles->Append(aO.Z());
462 gp_Ax3 anAx3 = anInfo.Position();
463 gp_Dir aD = anAx3.Direction();
464 theDoubles->Append(aD.X());
465 theDoubles->Append(aD.Y());
466 theDoubles->Append(aD.Z());
468 theDoubles->Append(anInfo.Radius1());
469 theDoubles->Append(anInfo.Radius2());
470 theDoubles->Append(anInfo.Height());
473 case GEOMAlgo_KN_DISKCIRCLE:
474 // (+) geompy.kind.DISK_CIRCLE xc yc zc dx dy dz R
476 aKind = SK_DISK_CIRCLE;
478 gp_Pnt aC = anInfo.Location();
479 theDoubles->Append(aC.X());
480 theDoubles->Append(aC.Y());
481 theDoubles->Append(aC.Z());
483 gp_Ax3 anAx3 = anInfo.Position();
484 gp_Dir aD = anAx3.Direction();
485 theDoubles->Append(aD.X());
486 theDoubles->Append(aD.Y());
487 theDoubles->Append(aD.Z());
489 theDoubles->Append(anInfo.Radius1());
492 case GEOMAlgo_KN_DISKELLIPSE:
493 // (+) geompy.kind.DISK_ELLIPSE xc yc zc dx dy dz R_1 R_2
495 aKind = SK_DISK_ELLIPSE;
497 gp_Pnt aC = anInfo.Location();
498 theDoubles->Append(aC.X());
499 theDoubles->Append(aC.Y());
500 theDoubles->Append(aC.Z());
502 gp_Ax3 anAx3 = anInfo.Position();
503 gp_Dir aD = anAx3.Direction();
504 theDoubles->Append(aD.X());
505 theDoubles->Append(aD.Y());
506 theDoubles->Append(aD.Z());
508 theDoubles->Append(anInfo.Radius1());
509 theDoubles->Append(anInfo.Radius2());
512 case GEOMAlgo_KN_RECTANGLE:
513 case GEOMAlgo_KN_TRIANGLE:
514 case GEOMAlgo_KN_QUADRANGLE:
515 case GEOMAlgo_KN_POLYGON:
516 // (+) geompy.kind.POLYGON xo yo zo dx dy dz nb_edges nb_vertices
520 gp_Pnt aO = anInfo.Location();
521 theDoubles->Append(aO.X());
522 theDoubles->Append(aO.Y());
523 theDoubles->Append(aO.Z());
525 gp_Ax3 anAx3 = anInfo.Position();
526 gp_Dir aD = anAx3.Direction();
527 theDoubles->Append(aD.X());
528 theDoubles->Append(aD.Y());
529 theDoubles->Append(aD.Z());
531 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
532 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
535 case GEOMAlgo_KN_PLANE: // infinite
536 // (+) geompy.kind.PLANE xo yo zo dx dy dz
540 gp_Pnt aC = anInfo.Location();
541 theDoubles->Append(aC.X());
542 theDoubles->Append(aC.Y());
543 theDoubles->Append(aC.Z());
545 gp_Ax3 anAx3 = anInfo.Position();
546 gp_Dir aD = anAx3.Direction();
547 theDoubles->Append(aD.X());
548 theDoubles->Append(aD.Y());
549 theDoubles->Append(aD.Z());
553 if (anInfo.KindOfShape() == GEOMAlgo_KS_PLANE) {
554 // (+) geompy.kind.PLANAR xo yo zo dx dy dz nb_edges nb_vertices
558 gp_Pnt aC = anInfo.Location();
559 theDoubles->Append(aC.X());
560 theDoubles->Append(aC.Y());
561 theDoubles->Append(aC.Z());
563 gp_Ax3 anAx3 = anInfo.Position();
564 gp_Dir aD = anAx3.Direction();
565 theDoubles->Append(aD.X());
566 theDoubles->Append(aD.Y());
567 theDoubles->Append(aD.Z());
569 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
570 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
573 // ??? geompy.kind.FACE nb_edges nb_vertices _surface_type_id_
574 // (+) geompy.kind.FACE nb_edges nb_vertices
576 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
577 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
587 GEOMAlgo_KindOfName aKN = anInfo.KindOfName();
589 case GEOMAlgo_KN_CIRCLE:
591 // (+) geompy.kind.CIRCLE xc yc zc dx dy dz R
594 gp_Pnt aC = anInfo.Location();
595 theDoubles->Append(aC.X());
596 theDoubles->Append(aC.Y());
597 theDoubles->Append(aC.Z());
599 gp_Ax3 anAx3 = anInfo.Position();
600 gp_Dir aD = anAx3.Direction();
601 theDoubles->Append(aD.X());
602 theDoubles->Append(aD.Y());
603 theDoubles->Append(aD.Z());
605 theDoubles->Append(anInfo.Radius1());
608 case GEOMAlgo_KN_ARCCIRCLE:
610 // (+) geompy.kind.ARC_CIRCLE xc yc zc dx dy dz R x1 y1 z1 x2 y2 z2
611 aKind = SK_ARC_CIRCLE;
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());
626 gp_Pnt aP1 = anInfo.Pnt1();
627 theDoubles->Append(aP1.X());
628 theDoubles->Append(aP1.Y());
629 theDoubles->Append(aP1.Z());
631 gp_Pnt aP2 = anInfo.Pnt2();
632 theDoubles->Append(aP2.X());
633 theDoubles->Append(aP2.Y());
634 theDoubles->Append(aP2.Z());
637 case GEOMAlgo_KN_ELLIPSE:
639 // (+) geompy.kind.ELLIPSE xc yc zc dx dy dz R_1 R_2
642 gp_Pnt aC = anInfo.Location();
643 theDoubles->Append(aC.X());
644 theDoubles->Append(aC.Y());
645 theDoubles->Append(aC.Z());
647 gp_Ax3 anAx3 = anInfo.Position();
648 gp_Dir aD = anAx3.Direction();
649 theDoubles->Append(aD.X());
650 theDoubles->Append(aD.Y());
651 theDoubles->Append(aD.Z());
653 theDoubles->Append(anInfo.Radius1());
654 theDoubles->Append(anInfo.Radius2());
657 case GEOMAlgo_KN_ARCELLIPSE:
659 // (+) geompy.kind.ARC_ELLIPSE xc yc zc dx dy dz R_1 R_2 x1 y1 z1 x2 y2 z2
660 aKind = SK_ARC_ELLIPSE;
662 gp_Pnt aC = anInfo.Location();
663 theDoubles->Append(aC.X());
664 theDoubles->Append(aC.Y());
665 theDoubles->Append(aC.Z());
667 gp_Ax3 anAx3 = anInfo.Position();
668 gp_Dir aD = anAx3.Direction();
669 theDoubles->Append(aD.X());
670 theDoubles->Append(aD.Y());
671 theDoubles->Append(aD.Z());
673 theDoubles->Append(anInfo.Radius1());
674 theDoubles->Append(anInfo.Radius2());
676 gp_Pnt aP1 = anInfo.Pnt1();
677 theDoubles->Append(aP1.X());
678 theDoubles->Append(aP1.Y());
679 theDoubles->Append(aP1.Z());
681 gp_Pnt aP2 = anInfo.Pnt2();
682 theDoubles->Append(aP2.X());
683 theDoubles->Append(aP2.Y());
684 theDoubles->Append(aP2.Z());
687 case GEOMAlgo_KN_LINE:
689 // ??? geompy.kind.LINE x1 y1 z1 x2 y2 z2
690 // (+) geompy.kind.LINE x1 y1 z1 dx dy dz
693 gp_Pnt aO = anInfo.Location();
694 theDoubles->Append(aO.X());
695 theDoubles->Append(aO.Y());
696 theDoubles->Append(aO.Z());
698 gp_Dir aD = anInfo.Direction();
699 theDoubles->Append(aD.X());
700 theDoubles->Append(aD.Y());
701 theDoubles->Append(aD.Z());
704 case GEOMAlgo_KN_SEGMENT:
706 // (+) geompy.kind.SEGMENT x1 y1 z1 x2 y2 z2
709 gp_Pnt aP1 = anInfo.Pnt1();
710 theDoubles->Append(aP1.X());
711 theDoubles->Append(aP1.Y());
712 theDoubles->Append(aP1.Z());
714 gp_Pnt aP2 = anInfo.Pnt2();
715 theDoubles->Append(aP2.X());
716 theDoubles->Append(aP2.Y());
717 theDoubles->Append(aP2.Z());
721 // ??? geompy.kind.EDGE nb_vertices _curve_type_id_
722 // (+) geompy.kind.EDGE nb_vertices
723 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
730 // (+) geompy.kind.VERTEX x y z
733 gp_Pnt aP = anInfo.Location();
734 theDoubles->Append(aP.X());
735 theDoubles->Append(aP.Y());
736 theDoubles->Append(aP.Z());
745 //=============================================================================
746 /*! Get LCS, corresponding to the given shape.
747 * Origin of the LCS is situated at the shape's center of mass.
748 * Axes of the LCS are obtained from shape's location or,
749 * if the shape is a planar face, from position of its plane.
751 //=============================================================================
752 gp_Ax3 GEOMImpl_IMeasureOperations::GetPosition (const TopoDS_Shape& theShape)
756 if (theShape.IsNull())
760 aResult.Transform(theShape.Location().Transformation());
761 if (theShape.ShapeType() == TopAbs_FACE) {
762 Handle(Geom_Surface) aGS = BRep_Tool::Surface(TopoDS::Face(theShape));
763 if (!aGS.IsNull() && aGS->IsKind(STANDARD_TYPE(Geom_Plane))) {
764 Handle(Geom_Plane) aGPlane = Handle(Geom_Plane)::DownCast(aGS);
765 gp_Pln aPln = aGPlane->Pln();
766 aResult = aPln.Position();
772 if (theShape.ShapeType() == TopAbs_VERTEX) {
773 aPnt = BRep_Tool::Pnt(TopoDS::Vertex(theShape));
776 GProp_GProps aSystem;
777 if (theShape.ShapeType() == TopAbs_EDGE || theShape.ShapeType() == TopAbs_WIRE)
778 BRepGProp::LinearProperties(theShape, aSystem);
779 else if (theShape.ShapeType() == TopAbs_FACE || theShape.ShapeType() == TopAbs_SHELL)
780 BRepGProp::SurfaceProperties(theShape, aSystem);
782 BRepGProp::VolumeProperties(theShape, aSystem);
784 aPnt = aSystem.CentreOfMass();
787 aResult.SetLocation(aPnt);
792 //=============================================================================
796 //=============================================================================
797 void GEOMImpl_IMeasureOperations::GetPosition
798 (Handle(GEOM_Object) theShape,
799 Standard_Real& Ox, Standard_Real& Oy, Standard_Real& Oz,
800 Standard_Real& Zx, Standard_Real& Zy, Standard_Real& Zz,
801 Standard_Real& Xx, Standard_Real& Xy, Standard_Real& Xz)
805 //Set default values: global CS
806 Ox = Oy = Oz = Zx = Zy = Xy = Xz = 0.;
809 if (theShape.IsNull()) return;
811 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
812 if (aRefShape.IsNull()) return;
814 TopoDS_Shape aShape = aRefShape->GetValue();
815 if (aShape.IsNull()) {
816 SetErrorCode("The Objects has NULL Shape");
821 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
825 gp_Ax3 anAx3 = GetPosition(aShape);
827 gp_Pnt anOri = anAx3.Location();
828 gp_Dir aDirZ = anAx3.Direction();
829 gp_Dir aDirX = anAx3.XDirection();
832 anOri.Coord(Ox, Oy, Oz);
833 aDirZ.Coord(Zx, Zy, Zz);
834 aDirX.Coord(Xx, Xy, Xz);
836 catch (Standard_Failure) {
837 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
838 SetErrorCode(aFail->GetMessageString());
845 //=============================================================================
849 //=============================================================================
850 Handle(GEOM_Object) GEOMImpl_IMeasureOperations::GetCentreOfMass
851 (Handle(GEOM_Object) theShape)
855 if (theShape.IsNull()) return NULL;
857 //Add a new CentreOfMass object
858 Handle(GEOM_Object) aCDG = GetEngine()->AddObject(GetDocID(), GEOM_CDG);
860 //Add a new CentreOfMass function
861 Handle(GEOM_Function) aFunction =
862 aCDG->AddFunction(GEOMImpl_MeasureDriver::GetID(), CDG_MEASURE);
863 if (aFunction.IsNull()) return NULL;
865 //Check if the function is set correctly
866 if (aFunction->GetDriverGUID() != GEOMImpl_MeasureDriver::GetID()) return NULL;
868 GEOMImpl_IMeasure aCI (aFunction);
870 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
871 if (aRefShape.IsNull()) return NULL;
873 aCI.SetBase(aRefShape);
875 //Compute the CentreOfMass value
877 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
880 if (!GetSolver()->ComputeFunction(aFunction)) {
881 SetErrorCode("Measure driver failed to compute centre of mass");
885 catch (Standard_Failure) {
886 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
887 SetErrorCode(aFail->GetMessageString());
891 //Make a Python command
892 GEOM::TPythonDump(aFunction) << aCDG << " = geompy.MakeCDG(" << theShape << ")";
898 //=============================================================================
902 //=============================================================================
903 Handle(GEOM_Object) GEOMImpl_IMeasureOperations::GetVertexByIndex
904 (Handle(GEOM_Object) theShape,
905 Standard_Integer theIndex)
909 if (theShape.IsNull()) return NULL;
911 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
912 if (aRefShape.IsNull()) return NULL;
914 //Add a new Vertex object
915 Handle(GEOM_Object) aVertex = GetEngine()->AddObject(GetDocID(), GEOM_POINT);
918 Handle(GEOM_Function) aFunction =
919 aVertex->AddFunction(GEOMImpl_MeasureDriver::GetID(), VERTEX_BY_INDEX);
920 if (aFunction.IsNull()) return NULL;
922 //Check if the function is set correctly
923 if (aFunction->GetDriverGUID() != GEOMImpl_MeasureDriver::GetID()) return NULL;
925 GEOMImpl_IMeasure aCI (aFunction);
926 aCI.SetBase(aRefShape);
927 aCI.SetIndex(theIndex);
931 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
934 if (!GetSolver()->ComputeFunction(aFunction)) {
935 SetErrorCode("Vertex by index driver failed.");
939 catch (Standard_Failure) {
940 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
941 SetErrorCode(aFail->GetMessageString());
945 //Make a Python command
946 GEOM::TPythonDump(aFunction) << aVertex << " = geompy.GetVertexByIndex(" << theShape << ", " << theIndex << ")";
952 //=============================================================================
956 //=============================================================================
957 Handle(GEOM_Object) GEOMImpl_IMeasureOperations::GetNormal
958 (Handle(GEOM_Object) theFace,
959 Handle(GEOM_Object) theOptionalPoint)
963 if (theFace.IsNull()) return NULL;
965 //Add a new Normale object
966 Handle(GEOM_Object) aNorm = GetEngine()->AddObject(GetDocID(), GEOM_VECTOR);
968 //Add a new Normale function
969 Handle(GEOM_Function) aFunction =
970 aNorm->AddFunction(GEOMImpl_MeasureDriver::GetID(), VECTOR_FACE_NORMALE);
971 if (aFunction.IsNull()) return NULL;
973 //Check if the function is set correctly
974 if (aFunction->GetDriverGUID() != GEOMImpl_MeasureDriver::GetID()) return NULL;
976 GEOMImpl_IMeasure aCI (aFunction);
978 Handle(GEOM_Function) aFace = theFace->GetLastFunction();
979 if (aFace.IsNull()) return NULL;
983 if (!theOptionalPoint.IsNull()) {
984 Handle(GEOM_Function) anOptPnt = theOptionalPoint->GetLastFunction();
985 aCI.SetPoint(anOptPnt);
988 //Compute the Normale value
990 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
993 if (!GetSolver()->ComputeFunction(aFunction)) {
994 SetErrorCode("Measure driver failed to compute normake of face");
998 catch (Standard_Failure) {
999 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1000 SetErrorCode(aFail->GetMessageString());
1004 //Make a Python command
1005 GEOM::TPythonDump pd (aFunction);
1006 pd << aNorm << " = geompy.GetNormal(" << theFace;
1007 if (!theOptionalPoint.IsNull()) {
1008 pd << ", " << theOptionalPoint;
1016 //=============================================================================
1018 * GetBasicProperties
1020 //=============================================================================
1021 void GEOMImpl_IMeasureOperations::GetBasicProperties (Handle(GEOM_Object) theShape,
1022 Standard_Real& theLength,
1023 Standard_Real& theSurfArea,
1024 Standard_Real& theVolume)
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 LProps, SProps;
1042 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
1045 BRepGProp::LinearProperties(aShape, LProps);
1046 theLength = LProps.Mass();
1048 BRepGProp::SurfaceProperties(aShape, SProps);
1049 theSurfArea = SProps.Mass();
1052 if (aShape.ShapeType() < TopAbs_SHELL) {
1053 for (TopExp_Explorer Exp (aShape, TopAbs_SOLID); Exp.More(); Exp.Next()) {
1054 GProp_GProps VProps;
1055 BRepGProp::VolumeProperties(Exp.Current(), VProps);
1056 theVolume += VProps.Mass();
1060 catch (Standard_Failure) {
1061 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1062 SetErrorCode(aFail->GetMessageString());
1069 //=============================================================================
1073 //=============================================================================
1074 void GEOMImpl_IMeasureOperations::GetInertia
1075 (Handle(GEOM_Object) theShape,
1076 Standard_Real& I11, Standard_Real& I12, Standard_Real& I13,
1077 Standard_Real& I21, Standard_Real& I22, Standard_Real& I23,
1078 Standard_Real& I31, Standard_Real& I32, Standard_Real& I33,
1079 Standard_Real& Ix , Standard_Real& Iy , Standard_Real& Iz)
1083 if (theShape.IsNull()) return;
1085 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1086 if (aRefShape.IsNull()) return;
1088 TopoDS_Shape aShape = aRefShape->GetValue();
1089 if (aShape.IsNull()) {
1090 SetErrorCode("The Objects has NULL Shape");
1094 //Compute the parameters
1095 GProp_GProps System;
1098 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
1101 if (aShape.ShapeType() == TopAbs_VERTEX ||
1102 aShape.ShapeType() == TopAbs_EDGE ||
1103 aShape.ShapeType() == TopAbs_WIRE) {
1104 BRepGProp::LinearProperties(aShape, System);
1105 } else if (aShape.ShapeType() == TopAbs_FACE ||
1106 aShape.ShapeType() == TopAbs_SHELL) {
1107 BRepGProp::SurfaceProperties(aShape, System);
1109 BRepGProp::VolumeProperties(aShape, System);
1111 gp_Mat I = System.MatrixOfInertia();
1125 GProp_PrincipalProps Pr = System.PrincipalProperties();
1126 Pr.Moments(Ix,Iy,Iz);
1128 catch (Standard_Failure) {
1129 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1130 SetErrorCode(aFail->GetMessageString());
1137 //=============================================================================
1141 //=============================================================================
1142 void GEOMImpl_IMeasureOperations::GetBoundingBox
1143 (Handle(GEOM_Object) theShape,
1144 Standard_Real& Xmin, Standard_Real& Xmax,
1145 Standard_Real& Ymin, Standard_Real& Ymax,
1146 Standard_Real& Zmin, Standard_Real& Zmax)
1150 if (theShape.IsNull()) return;
1152 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1153 if (aRefShape.IsNull()) return;
1155 TopoDS_Shape aShape = aRefShape->GetValue();
1156 if (aShape.IsNull()) {
1157 SetErrorCode("The Objects has NULL Shape");
1161 //Compute the parameters
1165 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
1168 BRepBndLib::Add(aShape, B);
1169 B.Get(Xmin, Ymin, Zmin, Xmax, Ymax, Zmax);
1171 catch (Standard_Failure) {
1172 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1173 SetErrorCode(aFail->GetMessageString());
1180 //=============================================================================
1184 //=============================================================================
1185 void GEOMImpl_IMeasureOperations::GetTolerance
1186 (Handle(GEOM_Object) theShape,
1187 Standard_Real& FaceMin, Standard_Real& FaceMax,
1188 Standard_Real& EdgeMin, Standard_Real& EdgeMax,
1189 Standard_Real& VertMin, Standard_Real& VertMax)
1193 if (theShape.IsNull()) return;
1195 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1196 if (aRefShape.IsNull()) return;
1198 TopoDS_Shape aShape = aRefShape->GetValue();
1199 if (aShape.IsNull()) {
1200 SetErrorCode("The Objects has NULL Shape");
1204 //Compute the parameters
1206 FaceMin = EdgeMin = VertMin = RealLast();
1207 FaceMax = EdgeMax = VertMax = -RealLast();
1210 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
1213 for (TopExp_Explorer ExF (aShape, TopAbs_FACE); ExF.More(); ExF.Next()) {
1214 TopoDS_Face Face = TopoDS::Face(ExF.Current());
1215 T = BRep_Tool::Tolerance(Face);
1221 for (TopExp_Explorer ExE (aShape, TopAbs_EDGE); ExE.More(); ExE.Next()) {
1222 TopoDS_Edge Edge = TopoDS::Edge(ExE.Current());
1223 T = BRep_Tool::Tolerance(Edge);
1229 for (TopExp_Explorer ExV (aShape, TopAbs_VERTEX); ExV.More(); ExV.Next()) {
1230 TopoDS_Vertex Vertex = TopoDS::Vertex(ExV.Current());
1231 T = BRep_Tool::Tolerance(Vertex);
1238 catch (Standard_Failure) {
1239 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1240 SetErrorCode(aFail->GetMessageString());
1247 //=============================================================================
1251 //=============================================================================
1252 bool GEOMImpl_IMeasureOperations::CheckShape (Handle(GEOM_Object) theShape,
1253 const Standard_Boolean theIsCheckGeom,
1254 TCollection_AsciiString& theDump)
1258 if (theShape.IsNull()) return false;
1260 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1261 if (aRefShape.IsNull()) return false;
1263 TopoDS_Shape aShape = aRefShape->GetValue();
1264 if (aShape.IsNull()) {
1265 SetErrorCode("The Objects has NULL Shape");
1269 //Compute the parameters
1270 bool isValid = false;
1272 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
1275 BRepCheck_Analyzer ana (aShape, theIsCheckGeom);
1276 if (ana.IsValid()) {
1280 StructuralDump(ana, aShape, theDump);
1283 catch (Standard_Failure) {
1284 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1285 SetErrorCode(aFail->GetMessageString());
1293 //=============================================================================
1297 //=============================================================================
1298 TCollection_AsciiString GEOMImpl_IMeasureOperations::WhatIs (Handle(GEOM_Object) theShape)
1302 TCollection_AsciiString Astr;
1304 if (theShape.IsNull()) return Astr;
1306 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1307 if (aRefShape.IsNull()) return Astr;
1309 TopoDS_Shape aShape = aRefShape->GetValue();
1310 if (aShape.IsNull()) {
1311 SetErrorCode("The Objects has NULL Shape");
1315 //Compute the parameters
1316 if (aShape.ShapeType() == TopAbs_EDGE) {
1317 if (BRep_Tool::Degenerated(TopoDS::Edge(aShape))) {
1318 Astr = Astr + " It is a degenerated edge \n";
1322 Astr = Astr + " Number of sub-shapes : \n";
1325 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
1328 int iType, nbTypes [TopAbs_SHAPE];
1329 for (iType = 0; iType < TopAbs_SHAPE; ++iType)
1331 nbTypes[aShape.ShapeType()]++;
1333 TopTools_MapOfShape aMapOfShape;
1334 aMapOfShape.Add(aShape);
1335 TopTools_ListOfShape aListOfShape;
1336 aListOfShape.Append(aShape);
1338 TopTools_ListIteratorOfListOfShape itL (aListOfShape);
1339 for (; itL.More(); itL.Next()) {
1340 TopoDS_Iterator it (itL.Value());
1341 for (; it.More(); it.Next()) {
1342 TopoDS_Shape s = it.Value();
1343 if (aMapOfShape.Add(s)) {
1344 aListOfShape.Append(s);
1345 nbTypes[s.ShapeType()]++;
1350 Astr = Astr + " VERTEX : " + TCollection_AsciiString(nbTypes[TopAbs_VERTEX]) + "\n";
1351 Astr = Astr + " EDGE : " + TCollection_AsciiString(nbTypes[TopAbs_EDGE]) + "\n";
1352 Astr = Astr + " WIRE : " + TCollection_AsciiString(nbTypes[TopAbs_WIRE]) + "\n";
1353 Astr = Astr + " FACE : " + TCollection_AsciiString(nbTypes[TopAbs_FACE]) + "\n";
1354 Astr = Astr + " SHELL : " + TCollection_AsciiString(nbTypes[TopAbs_SHELL]) + "\n";
1355 Astr = Astr + " SOLID : " + TCollection_AsciiString(nbTypes[TopAbs_SOLID]) + "\n";
1356 Astr = Astr + " COMPSOLID : " + TCollection_AsciiString(nbTypes[TopAbs_COMPSOLID]) + "\n";
1357 Astr = Astr + " COMPOUND : " + TCollection_AsciiString(nbTypes[TopAbs_COMPOUND]) + "\n";
1358 Astr = Astr + " SHAPE : " + TCollection_AsciiString(aMapOfShape.Extent());
1360 catch (Standard_Failure) {
1361 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1362 SetErrorCode(aFail->GetMessageString());
1371 //=======================================================================
1372 //function : CheckSingularCase
1373 //purpose : auxilary for GetMinDistance()
1374 // workaround for bugs 19899, 19908 and 19910 from Mantis
1375 //=======================================================================
1376 static double CheckSingularCase(const TopoDS_Shape& aSh1,
1377 const TopoDS_Shape& aSh2,
1378 gp_Pnt& Ptmp1, gp_Pnt& Ptmp2)
1380 bool IsChange1 = false;
1381 double AddDist1 = 0.0;
1382 TopExp_Explorer anExp;
1383 TopoDS_Shape tmpSh1, tmpSh2;
1385 for ( anExp.Init( aSh1, TopAbs_FACE ); anExp.More(); anExp.Next() ) {
1387 tmpSh1 = anExp.Current();
1390 TopoDS_Shape sh = aSh1;
1391 while(sh.ShapeType()==TopAbs_COMPOUND) {
1392 TopoDS_Iterator it(sh);
1395 Handle(Geom_Surface) S = BRep_Tool::Surface(TopoDS::Face(tmpSh1));
1396 if( S->IsKind(STANDARD_TYPE(Geom_SphericalSurface)) ||
1397 S->IsKind(STANDARD_TYPE(Geom_ToroidalSurface)) ) {
1398 if( sh.ShapeType()==TopAbs_SHELL || sh.ShapeType()==TopAbs_FACE ) {
1401 // changes for 0020677: EDF 1219 GEOM: MinDistance gives 0 instead of 20.88
1402 //S->Bounds(U1,U2,V1,V2); changed by
1403 ShapeAnalysis::GetFaceUVBounds(TopoDS::Face(tmpSh1),U1,U2,V1,V2);
1404 // end of changes for 020677 (dmv)
1405 Handle(Geom_RectangularTrimmedSurface) TrS1 =
1406 new Geom_RectangularTrimmedSurface(S,U1,(U1+U2)/2.,V1,V2);
1407 Handle(Geom_RectangularTrimmedSurface) TrS2 =
1408 new Geom_RectangularTrimmedSurface(S,(U1+U2)/2.,U2,V1,V2);
1411 TopoDS_Compound Comp;
1412 B.MakeCompound(Comp);
1413 B.MakeFace(F1,TrS1,1.e-7);
1415 B.MakeFace(F2,TrS2,1.e-7);
1417 Handle(ShapeFix_Shape) sfs = new ShapeFix_Shape;
1419 sfs->SetPrecision(1.e-6);
1420 sfs->SetMaxTolerance(1.0);
1422 tmpSh1 = sfs->Shape();
1426 if( S->IsKind(STANDARD_TYPE(Geom_SphericalSurface)) ) {
1427 Handle(Geom_SphericalSurface) SS = Handle(Geom_SphericalSurface)::DownCast(S);
1428 gp_Pnt PC = SS->Location();
1431 B.MakeVertex(V,PC,1.e-7);
1433 AddDist1 = SS->Radius();
1437 Handle(Geom_ToroidalSurface) TS = Handle(Geom_ToroidalSurface)::DownCast(S);
1438 gp_Ax3 ax3 = TS->Position();
1439 Handle(Geom_Circle) C = new Geom_Circle(ax3.Ax2(),TS->MajorRadius());
1442 B.MakeEdge(E,C,1.e-7);
1444 AddDist1 = TS->MinorRadius();
1454 bool IsChange2 = false;
1455 double AddDist2 = 0.0;
1457 for ( anExp.Init( aSh2, TopAbs_FACE ); anExp.More(); anExp.Next() ) {
1459 tmpSh2 = anExp.Current();
1462 TopoDS_Shape sh = aSh2;
1463 while(sh.ShapeType()==TopAbs_COMPOUND) {
1464 TopoDS_Iterator it(sh);
1467 Handle(Geom_Surface) S = BRep_Tool::Surface(TopoDS::Face(tmpSh2));
1468 if( S->IsKind(STANDARD_TYPE(Geom_SphericalSurface)) ||
1469 S->IsKind(STANDARD_TYPE(Geom_ToroidalSurface)) ) {
1470 if( sh.ShapeType()==TopAbs_SHELL || sh.ShapeType()==TopAbs_FACE ) {
1473 S->Bounds(U1,U2,V1,V2);
1474 Handle(Geom_RectangularTrimmedSurface) TrS1 =
1475 new Geom_RectangularTrimmedSurface(S,U1,(U1+U2)/2.,V1,V2);
1476 Handle(Geom_RectangularTrimmedSurface) TrS2 =
1477 new Geom_RectangularTrimmedSurface(S,(U1+U2)/2.,U2,V1,V2);
1480 TopoDS_Compound Comp;
1481 B.MakeCompound(Comp);
1482 B.MakeFace(F1,TrS1,1.e-7);
1484 B.MakeFace(F2,TrS2,1.e-7);
1486 Handle(ShapeFix_Shape) sfs = new ShapeFix_Shape;
1488 sfs->SetPrecision(1.e-6);
1489 sfs->SetMaxTolerance(1.0);
1491 tmpSh2 = sfs->Shape();
1495 if( S->IsKind(STANDARD_TYPE(Geom_SphericalSurface)) ) {
1496 Handle(Geom_SphericalSurface) SS = Handle(Geom_SphericalSurface)::DownCast(S);
1497 gp_Pnt PC = SS->Location();
1500 B.MakeVertex(V,PC,1.e-7);
1502 AddDist2 = SS->Radius();
1505 else if( S->IsKind(STANDARD_TYPE(Geom_ToroidalSurface)) ) {
1506 Handle(Geom_ToroidalSurface) TS = Handle(Geom_ToroidalSurface)::DownCast(S);
1507 gp_Ax3 ax3 = TS->Position();
1508 Handle(Geom_Circle) C = new Geom_Circle(ax3.Ax2(),TS->MajorRadius());
1511 B.MakeEdge(E,C,1.e-7);
1513 AddDist2 = TS->MinorRadius();
1524 if( !IsChange1 && !IsChange2 )
1527 BRepExtrema_DistShapeShape dst(tmpSh1,tmpSh2);
1529 double MinDist = 1.e9;
1530 gp_Pnt PMin1, PMin2, P1, P2;
1531 for (int i = 1; i <= dst.NbSolution(); i++) {
1532 P1 = dst.PointOnShape1(i);
1533 P2 = dst.PointOnShape2(i);
1534 Standard_Real Dist = P1.Distance(P2);
1535 if (MinDist > Dist) {
1546 gp_Dir aDir(gp_Vec(PMin1,PMin2));
1547 if( MinDist > (AddDist1+AddDist2) ) {
1548 Ptmp1 = gp_Pnt( PMin1.X() + aDir.X()*AddDist1,
1549 PMin1.Y() + aDir.Y()*AddDist1,
1550 PMin1.Z() + aDir.Z()*AddDist1 );
1551 Ptmp2 = gp_Pnt( PMin2.X() - aDir.X()*AddDist2,
1552 PMin2.Y() - aDir.Y()*AddDist2,
1553 PMin2.Z() - aDir.Z()*AddDist2 );
1554 return (MinDist - AddDist1 - AddDist2);
1557 if( AddDist1 > 0 ) {
1558 Ptmp1 = gp_Pnt( PMin1.X() + aDir.X()*AddDist1,
1559 PMin1.Y() + aDir.Y()*AddDist1,
1560 PMin1.Z() + aDir.Z()*AddDist1 );
1564 Ptmp2 = gp_Pnt( PMin2.X() - aDir.X()*AddDist2,
1565 PMin2.Y() - aDir.Y()*AddDist2,
1566 PMin2.Z() - aDir.Z()*AddDist2 );
1571 double res = MinDist - AddDist1 - AddDist2;
1572 if(res<0.) res = 0.0;
1578 static bool CheckSingularCase(const TopoDS_Shape& aSh1,
1579 const TopoDS_Shape& aSh2,
1582 TopExp_Explorer anExp;
1583 TopoDS_Shape tmpSh1, tmpSh2;
1585 for ( anExp.Init( aSh1, TopAbs_FACE ); anExp.More(); anExp.Next() ) {
1587 tmpSh1 = anExp.Current();
1590 Handle(Geom_Surface) S1 = BRep_Tool::Surface(TopoDS::Face(tmpSh1));
1591 if( S1->IsKind(STANDARD_TYPE(Geom_SphericalSurface)) ||
1592 S1->IsKind(STANDARD_TYPE(Geom_ToroidalSurface)) ) {
1594 for ( anExp.Init( aSh2, TopAbs_FACE ); anExp.More(); anExp.Next() ) {
1596 tmpSh2 = anExp.Current();
1597 Handle(Geom_Surface) S2 = BRep_Tool::Surface(TopoDS::Face(tmpSh2));
1598 GeomAPI_IntSS ISS(S1,S2,1.e-7);
1600 for(int i=1; i<=ISS.NbLines(); i++) {
1601 Handle(Geom_Curve) C3d = ISS.Line(i);
1604 B.MakeEdge(E,C3d,1.e-7);
1605 BRepExtrema_DistShapeShape dst(tmpSh2,E);
1607 gp_Pnt PMin1, PMin2, P1, P2;
1608 double MinDist = 1.e9;
1609 for (int i = 1; i <= dst.NbSolution(); i++) {
1610 P1 = dst.PointOnShape1(i);
1611 P2 = dst.PointOnShape2(i);
1612 Standard_Real Dist = P1.Distance(P2);
1613 if (MinDist > Dist) {
1627 for ( anExp.Init( aSh2, TopAbs_FACE ); anExp.More(); anExp.Next() ) {
1629 tmpSh1 = anExp.Current();
1632 Handle(Geom_Surface) S1 = BRep_Tool::Surface(TopoDS::Face(tmpSh1));
1633 if( S1->IsKind(STANDARD_TYPE(Geom_SphericalSurface)) ||
1634 S1->IsKind(STANDARD_TYPE(Geom_ToroidalSurface)) ) {
1636 for ( anExp.Init( aSh1, TopAbs_FACE ); anExp.More(); anExp.Next() ) {
1638 tmpSh2 = anExp.Current();
1639 Handle(Geom_Surface) S2 = BRep_Tool::Surface(TopoDS::Face(tmpSh2));
1640 GeomAPI_IntSS ISS(S1,S2,1.e-7);
1642 for(int i=1; i<=ISS.NbLines(); i++) {
1643 Handle(Geom_Curve) C3d = ISS.Line(i);
1646 B.MakeEdge(E,C3d,1.e-7);
1647 BRepExtrema_DistShapeShape dst(tmpSh2,E);
1650 double MinDist = 1.e9;
1651 for (int i = 1; i <= dst.NbSolution(); i++) {
1652 P1 = dst.PointOnShape1(i);
1653 P2 = dst.PointOnShape2(i);
1654 Standard_Real Dist = P1.Distance(P2);
1655 if (MinDist > Dist) {
1673 //=============================================================================
1677 //=============================================================================
1678 Standard_Real GEOMImpl_IMeasureOperations::GetMinDistance
1679 (Handle(GEOM_Object) theShape1, Handle(GEOM_Object) theShape2,
1680 Standard_Real& X1, Standard_Real& Y1, Standard_Real& Z1,
1681 Standard_Real& X2, Standard_Real& Y2, Standard_Real& Z2)
1684 Standard_Real MinDist = 1.e9;
1686 if (theShape1.IsNull() || theShape2.IsNull()) return MinDist;
1688 Handle(GEOM_Function) aRefShape1 = theShape1->GetLastFunction();
1689 Handle(GEOM_Function) aRefShape2 = theShape2->GetLastFunction();
1690 if (aRefShape1.IsNull() || aRefShape2.IsNull()) return MinDist;
1692 TopoDS_Shape aShape1 = aRefShape1->GetValue();
1693 TopoDS_Shape aShape2 = aRefShape2->GetValue();
1694 if (aShape1.IsNull() || aShape2.IsNull()) {
1695 SetErrorCode("One of Objects has NULL Shape");
1699 //Compute the parameters
1701 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
1705 // Issue 0020231: A min distance bug with torus and vertex.
1706 // Make GetMinDistance() return zero if a sole VERTEX is inside any of SOLIDs
1708 // which of shapes consists of only one vertex?
1709 TopExp_Explorer exp1(aShape1,TopAbs_VERTEX), exp2(aShape2,TopAbs_VERTEX);
1710 TopoDS_Shape V1 = exp1.More() ? exp1.Current() : TopoDS_Shape();
1711 TopoDS_Shape V2 = exp2.More() ? exp2.Current() : TopoDS_Shape();
1712 exp1.Next(); exp2.Next();
1713 if ( exp1.More() ) V1.Nullify();
1714 if ( exp2.More() ) V2.Nullify();
1715 // vertex and container of solids
1716 TopoDS_Shape V = V1.IsNull() ? V2 : V1;
1717 TopoDS_Shape S = V1.IsNull() ? aShape1 : aShape2;
1718 if ( !V.IsNull() ) {
1719 // classify vertex against solids
1720 gp_Pnt p = BRep_Tool::Pnt( TopoDS::Vertex( V ) );
1721 for ( exp1.Init( S, TopAbs_SOLID ); exp1.More(); exp1.Next() ) {
1722 BRepClass3d_SolidClassifier classifier( exp1.Current(), p, 1e-6);
1723 if ( classifier.State() == TopAbs_IN ) {
1724 p.Coord(X1, Y1, Z1);
1725 p.Coord(X2, Y2, Z2);
1731 // End Issue 0020231
1734 // additional workaround for bugs 19899, 19908 and 19910 from Mantis
1735 gp_Pnt Ptmp1, Ptmp2;
1736 double dist = CheckSingularCase(aShape1, aShape2, Ptmp1, Ptmp2);
1738 Ptmp1.Coord(X1, Y1, Z1);
1739 Ptmp2.Coord(X2, Y2, Z2);
1744 BRepExtrema_DistShapeShape dst (aShape1, aShape2);
1746 gp_Pnt PMin1, PMin2, P1, P2;
1748 for (int i = 1; i <= dst.NbSolution(); i++) {
1749 P1 = dst.PointOnShape1(i);
1750 P2 = dst.PointOnShape2(i);
1752 Standard_Real Dist = P1.Distance(P2);
1753 if (MinDist > Dist) {
1760 PMin1.Coord(X1, Y1, Z1);
1761 PMin2.Coord(X2, Y2, Z2);
1764 catch (Standard_Failure) {
1765 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1766 SetErrorCode(aFail->GetMessageString());
1774 //=======================================================================
1776 * Get coordinates of point
1778 //=======================================================================
1779 void GEOMImpl_IMeasureOperations::PointCoordinates (Handle(GEOM_Object) theShape,
1780 Standard_Real& theX, Standard_Real& theY, Standard_Real& theZ)
1784 if (theShape.IsNull())
1787 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1788 if (aRefShape.IsNull())
1791 TopoDS_Shape aShape = aRefShape->GetValue();
1792 if (aShape.IsNull() || aShape.ShapeType() != TopAbs_VERTEX)
1794 SetErrorCode( "Shape must be a vertex" );
1799 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
1802 gp_Pnt aPnt = BRep_Tool::Pnt( TopoDS::Vertex( aShape ) );
1809 catch (Standard_Failure)
1811 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1812 SetErrorCode( aFail->GetMessageString() );
1816 //=======================================================================
1818 * Compute angle (in degrees) between two lines
1820 //=======================================================================
1821 Standard_Real GEOMImpl_IMeasureOperations::GetAngle (Handle(GEOM_Object) theLine1,
1822 Handle(GEOM_Object) theLine2)
1826 Standard_Real anAngle = -1.0;
1828 if (theLine1.IsNull() || theLine2.IsNull())
1831 Handle(GEOM_Function) aRefLine1 = theLine1->GetLastFunction();
1832 Handle(GEOM_Function) aRefLine2 = theLine2->GetLastFunction();
1833 if (aRefLine1.IsNull() || aRefLine2.IsNull())
1836 TopoDS_Shape aLine1 = aRefLine1->GetValue();
1837 TopoDS_Shape aLine2 = aRefLine2->GetValue();
1838 if (aLine1.IsNull() || aLine2.IsNull() ||
1839 aLine1.ShapeType() != TopAbs_EDGE ||
1840 aLine2.ShapeType() != TopAbs_EDGE)
1842 SetErrorCode("Two edges must be given");
1847 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
1850 TopoDS_Edge E1 = TopoDS::Edge(aLine1);
1851 TopoDS_Edge E2 = TopoDS::Edge(aLine2);
1854 Handle(Geom_Curve) C1 = BRep_Tool::Curve(E1,fp,lp);
1855 Handle(Geom_Curve) C2 = BRep_Tool::Curve(E2,fp,lp);
1857 if ( C1.IsNull() || C2.IsNull() ||
1858 !C1->IsKind(STANDARD_TYPE(Geom_Line)) ||
1859 !C2->IsKind(STANDARD_TYPE(Geom_Line)))
1861 SetErrorCode("The edges must be linear");
1865 Handle(Geom_Line) L1 = Handle(Geom_Line)::DownCast(C1);
1866 Handle(Geom_Line) L2 = Handle(Geom_Line)::DownCast(C2);
1868 gp_Lin aLin1 = L1->Lin();
1869 gp_Lin aLin2 = L2->Lin();
1871 anAngle = aLin1.Angle(aLin2);
1872 anAngle /= PI180; // convert radians into degrees
1876 catch (Standard_Failure)
1878 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1879 SetErrorCode(aFail->GetMessageString());
1886 //=============================================================================
1888 * CurveCurvatureByParam
1890 //=============================================================================
1891 Standard_Real GEOMImpl_IMeasureOperations::CurveCurvatureByParam
1892 (Handle(GEOM_Object) theCurve, Standard_Real& theParam)
1895 Standard_Real aRes = -1.0;
1897 if(theCurve.IsNull()) return aRes;
1899 Handle(GEOM_Function) aRefShape = theCurve->GetLastFunction();
1900 if(aRefShape.IsNull()) return aRes;
1902 TopoDS_Shape aShape = aRefShape->GetValue();
1903 if(aShape.IsNull()) {
1904 SetErrorCode("One of Objects has NULL Shape");
1908 Standard_Real aFP, aLP, aP;
1909 Handle(Geom_Curve) aCurve = BRep_Tool::Curve(TopoDS::Edge(aShape), aFP, aLP);
1910 aP = aFP + (aLP - aFP) * theParam;
1912 if(aCurve.IsNull()) return aRes;
1916 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
1919 GeomLProp_CLProps Prop = GeomLProp_CLProps
1920 (aCurve, aP, 2, Precision::Confusion());
1921 aRes = fabs(Prop.Curvature());
1924 catch (Standard_Failure) {
1925 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1926 SetErrorCode(aFail->GetMessageString());
1930 if( aRes > Precision::Confusion() )
1939 //=============================================================================
1941 * CurveCurvatureByPoint
1943 //=============================================================================
1944 Standard_Real GEOMImpl_IMeasureOperations::CurveCurvatureByPoint
1945 (Handle(GEOM_Object) theCurve, Handle(GEOM_Object) thePoint)
1948 Standard_Real aRes = -1.0;
1950 if( theCurve.IsNull() || thePoint.IsNull() ) return aRes;
1952 Handle(GEOM_Function) aRefCurve = theCurve->GetLastFunction();
1953 Handle(GEOM_Function) aRefPoint = thePoint->GetLastFunction();
1954 if( aRefCurve.IsNull() || aRefPoint.IsNull() ) return aRes;
1956 TopoDS_Edge anEdge = TopoDS::Edge(aRefCurve->GetValue());
1957 TopoDS_Vertex aPnt = TopoDS::Vertex(aRefPoint->GetValue());
1958 if( anEdge.IsNull() || aPnt.IsNull() ) {
1959 SetErrorCode("One of Objects has NULL Shape");
1963 Standard_Real aFP, aLP;
1964 Handle(Geom_Curve) aCurve = BRep_Tool::Curve(anEdge, aFP, aLP);
1965 if(aCurve.IsNull()) return aRes;
1966 gp_Pnt aPoint = BRep_Tool::Pnt(aPnt);
1970 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
1973 GeomAPI_ProjectPointOnCurve PPC(aPoint, aCurve, aFP, aLP);
1974 if(PPC.NbPoints()>0) {
1975 GeomLProp_CLProps Prop = GeomLProp_CLProps
1976 (aCurve, PPC.LowerDistanceParameter(), 2, Precision::Confusion());
1977 aRes = fabs(Prop.Curvature());
1981 catch (Standard_Failure) {
1982 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1983 SetErrorCode(aFail->GetMessageString());
1987 if( aRes > Precision::Confusion() )
1996 //=============================================================================
1998 * getSurfaceCurvatures
2000 //=============================================================================
2001 Standard_Real GEOMImpl_IMeasureOperations::getSurfaceCurvatures
2002 (const Handle(Geom_Surface)& aSurf,
2003 Standard_Real theUParam,
2004 Standard_Real theVParam,
2005 Standard_Boolean theNeedMaxCurv)
2008 Standard_Real aRes = 1.0;
2010 if (aSurf.IsNull()) return aRes;
2013 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
2016 GeomLProp_SLProps Prop = GeomLProp_SLProps
2017 (aSurf, theUParam, theVParam, 2, Precision::Confusion());
2018 if(Prop.IsCurvatureDefined()) {
2019 if(Prop.IsUmbilic()) {
2020 //cout<<"is umbilic"<<endl;
2021 aRes = fabs(Prop.MeanCurvature());
2024 //cout<<"is not umbilic"<<endl;
2025 double c1 = fabs(Prop.MaxCurvature());
2026 double c2 = fabs(Prop.MinCurvature());
2035 catch (Standard_Failure) {
2036 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
2037 SetErrorCode(aFail->GetMessageString());
2041 if( fabs(aRes) > Precision::Confusion() )
2050 //=============================================================================
2052 * MaxSurfaceCurvatureByParam
2054 //=============================================================================
2055 Standard_Real GEOMImpl_IMeasureOperations::MaxSurfaceCurvatureByParam
2056 (Handle(GEOM_Object) theSurf,
2057 Standard_Real& theUParam,
2058 Standard_Real& theVParam)
2061 Standard_Real aRes = -1.0;
2063 if (theSurf.IsNull()) return aRes;
2065 Handle(GEOM_Function) aRefShape = theSurf->GetLastFunction();
2066 if(aRefShape.IsNull()) return aRes;
2068 TopoDS_Shape aShape = aRefShape->GetValue();
2069 if(aShape.IsNull()) {
2070 SetErrorCode("One of Objects has NULL Shape");
2074 TopoDS_Face F = TopoDS::Face(aShape);
2075 Handle(Geom_Surface) aSurf = BRep_Tool::Surface(F);
2077 //Compute the parameters
2078 Standard_Real U1,U2,V1,V2;
2079 ShapeAnalysis::GetFaceUVBounds(F,U1,U2,V1,V2);
2080 Standard_Real U = U1 + (U2-U1)*theUParam;
2081 Standard_Real V = V1 + (V2-V1)*theVParam;
2083 return getSurfaceCurvatures(aSurf, U, V, true);
2087 //=============================================================================
2089 * MaxSurfaceCurvatureByPoint
2091 //=============================================================================
2092 Standard_Real GEOMImpl_IMeasureOperations::MaxSurfaceCurvatureByPoint
2093 (Handle(GEOM_Object) theSurf, Handle(GEOM_Object) thePoint)
2096 Standard_Real aRes = -1.0;
2098 if( theSurf.IsNull() || thePoint.IsNull() ) return aRes;
2100 Handle(GEOM_Function) aRefShape = theSurf->GetLastFunction();
2101 Handle(GEOM_Function) aRefPoint = thePoint->GetLastFunction();
2102 if( aRefShape.IsNull() || aRefPoint.IsNull() ) return aRes;
2104 TopoDS_Face aFace = TopoDS::Face(aRefShape->GetValue());
2105 TopoDS_Vertex aPnt = TopoDS::Vertex(aRefPoint->GetValue());
2106 if( aFace.IsNull() || aPnt.IsNull() ) {
2107 SetErrorCode("One of Objects has NULL Shape");
2111 Handle(Geom_Surface) aSurf = BRep_Tool::Surface(aFace);
2112 if(aSurf.IsNull()) return aRes;
2113 gp_Pnt aPoint = BRep_Tool::Pnt(aPnt);
2115 //Compute the parameters
2116 ShapeAnalysis_Surface sas(aSurf);
2117 gp_Pnt2d UV = sas.ValueOfUV(aPoint,Precision::Confusion());
2119 return getSurfaceCurvatures(aSurf, UV.X(), UV.Y(), true);
2123 //=============================================================================
2125 * MinSurfaceCurvatureByParam
2127 //=============================================================================
2128 Standard_Real GEOMImpl_IMeasureOperations::MinSurfaceCurvatureByParam
2129 (Handle(GEOM_Object) theSurf,
2130 Standard_Real& theUParam,
2131 Standard_Real& theVParam)
2134 Standard_Real aRes = -1.0;
2136 if (theSurf.IsNull()) return aRes;
2138 Handle(GEOM_Function) aRefShape = theSurf->GetLastFunction();
2139 if(aRefShape.IsNull()) return aRes;
2141 TopoDS_Shape aShape = aRefShape->GetValue();
2142 if(aShape.IsNull()) {
2143 SetErrorCode("One of Objects has NULL Shape");
2147 TopoDS_Face F = TopoDS::Face(aShape);
2148 Handle(Geom_Surface) aSurf = BRep_Tool::Surface(F);
2150 //Compute the parameters
2151 Standard_Real U1,U2,V1,V2;
2152 ShapeAnalysis::GetFaceUVBounds(F,U1,U2,V1,V2);
2153 Standard_Real U = U1 + (U2-U1)*theUParam;
2154 Standard_Real V = V1 + (V2-V1)*theVParam;
2156 return getSurfaceCurvatures(aSurf, U, V, false);
2160 //=============================================================================
2162 * MinSurfaceCurvatureByPoint
2164 //=============================================================================
2165 Standard_Real GEOMImpl_IMeasureOperations::MinSurfaceCurvatureByPoint
2166 (Handle(GEOM_Object) theSurf, Handle(GEOM_Object) thePoint)
2169 Standard_Real aRes = -1.0;
2171 if( theSurf.IsNull() || thePoint.IsNull() ) return aRes;
2173 Handle(GEOM_Function) aRefShape = theSurf->GetLastFunction();
2174 Handle(GEOM_Function) aRefPoint = thePoint->GetLastFunction();
2175 if( aRefShape.IsNull() || aRefPoint.IsNull() ) return aRes;
2177 TopoDS_Face aFace = TopoDS::Face(aRefShape->GetValue());
2178 TopoDS_Vertex aPnt = TopoDS::Vertex(aRefPoint->GetValue());
2179 if( aFace.IsNull() || aPnt.IsNull() ) {
2180 SetErrorCode("One of Objects has NULL Shape");
2184 Handle(Geom_Surface) aSurf = BRep_Tool::Surface(aFace);
2185 if(aSurf.IsNull()) return aRes;
2186 gp_Pnt aPoint = BRep_Tool::Pnt(aPnt);
2188 //Compute the parameters
2189 ShapeAnalysis_Surface sas(aSurf);
2190 gp_Pnt2d UV = sas.ValueOfUV(aPoint,Precision::Confusion());
2192 return getSurfaceCurvatures(aSurf, UV.X(), UV.Y(), false);
2196 //=======================================================================
2197 //function : StructuralDump
2198 //purpose : Structural (data exchange) style of output.
2199 //=======================================================================
2200 void GEOMImpl_IMeasureOperations::StructuralDump (const BRepCheck_Analyzer& theAna,
2201 const TopoDS_Shape& theShape,
2202 TCollection_AsciiString& theDump)
2206 theDump += " -- The Shape has problems :\n";
2207 theDump += " Check Count\n";
2208 theDump += " ------------------------------------------------\n";
2210 Standard_Integer last_stat = (Standard_Integer)BRepCheck_CheckFail;
2211 Handle(TColStd_HArray1OfInteger) NbProblems =
2212 new TColStd_HArray1OfInteger(1, last_stat);
2213 for (i = 1; i <= last_stat; i++)
2214 NbProblems->SetValue(i,0);
2216 Handle(TopTools_HSequenceOfShape) sl;
2217 sl = new TopTools_HSequenceOfShape();
2218 TopTools_DataMapOfShapeListOfShape theMap;
2220 GetProblemShapes(theAna, theShape, sl, NbProblems, theMap);
2223 Standard_Integer count = 0;
2224 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidPointOnCurve);
2226 theDump += " Invalid Point on Curve ................... ";
2227 theDump += TCollection_AsciiString(count) + "\n";
2229 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidPointOnCurveOnSurface);
2231 theDump += " Invalid Point on CurveOnSurface .......... ";
2232 theDump += TCollection_AsciiString(count) + "\n";
2234 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidPointOnSurface);
2236 theDump += " Invalid Point on Surface ................. ";
2237 theDump += TCollection_AsciiString(count) + "\n";
2239 count = NbProblems->Value((Standard_Integer)BRepCheck_No3DCurve);
2241 theDump += " No 3D Curve .............................. ";
2242 theDump += TCollection_AsciiString(count) + "\n";
2244 count = NbProblems->Value((Standard_Integer)BRepCheck_Multiple3DCurve);
2246 theDump += " Multiple 3D Curve ........................ ";
2247 theDump += TCollection_AsciiString(count) + "\n";
2249 count = NbProblems->Value((Standard_Integer)BRepCheck_Invalid3DCurve);
2251 theDump += " Invalid 3D Curve ......................... ";
2252 theDump += TCollection_AsciiString(count) + "\n";
2254 count = NbProblems->Value((Standard_Integer)BRepCheck_NoCurveOnSurface);
2256 theDump += " No Curve on Surface ...................... ";
2257 theDump += TCollection_AsciiString(count) + "\n";
2259 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidCurveOnSurface);
2261 theDump += " Invalid Curve on Surface ................. ";
2262 theDump += TCollection_AsciiString(count) + "\n";
2264 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidCurveOnClosedSurface);
2266 theDump += " Invalid Curve on closed Surface .......... ";
2267 theDump += TCollection_AsciiString(count) + "\n";
2269 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidSameRangeFlag);
2271 theDump += " Invalid SameRange Flag ................... ";
2272 theDump += TCollection_AsciiString(count) + "\n";
2274 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidSameParameterFlag);
2276 theDump += " Invalid SameParameter Flag ............... ";
2277 theDump += TCollection_AsciiString(count) + "\n";
2279 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidDegeneratedFlag);
2281 theDump += " Invalid Degenerated Flag ................. ";
2282 theDump += TCollection_AsciiString(count) + "\n";
2284 count = NbProblems->Value((Standard_Integer)BRepCheck_FreeEdge);
2286 theDump += " Free Edge ................................ ";
2287 theDump += TCollection_AsciiString(count) + "\n";
2289 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidMultiConnexity);
2291 theDump += " Invalid MultiConnexity ................... ";
2292 theDump += TCollection_AsciiString(count) + "\n";
2294 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidRange);
2296 theDump += " Invalid Range ............................ ";
2297 theDump += TCollection_AsciiString(count) + "\n";
2299 count = NbProblems->Value((Standard_Integer)BRepCheck_EmptyWire);
2301 theDump += " Empty Wire ............................... ";
2302 theDump += TCollection_AsciiString(count) + "\n";
2304 count = NbProblems->Value((Standard_Integer)BRepCheck_RedundantEdge);
2306 theDump += " Redundant Edge ........................... ";
2307 theDump += TCollection_AsciiString(count) + "\n";
2309 count = NbProblems->Value((Standard_Integer)BRepCheck_SelfIntersectingWire);
2311 theDump += " Self Intersecting Wire ................... ";
2312 theDump += TCollection_AsciiString(count) + "\n";
2314 count = NbProblems->Value((Standard_Integer)BRepCheck_NoSurface);
2316 theDump += " No Surface ............................... ";
2317 theDump += TCollection_AsciiString(count) + "\n";
2319 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidWire);
2321 theDump += " Invalid Wire ............................. ";
2322 theDump += TCollection_AsciiString(count) + "\n";
2324 count = NbProblems->Value((Standard_Integer)BRepCheck_RedundantWire);
2326 theDump += " Redundant Wire ........................... ";
2327 theDump += TCollection_AsciiString(count) + "\n";
2329 count = NbProblems->Value((Standard_Integer)BRepCheck_IntersectingWires);
2331 theDump += " Intersecting Wires ....................... ";
2332 theDump += TCollection_AsciiString(count) + "\n";
2334 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidImbricationOfWires);
2336 theDump += " Invalid Imbrication of Wires ............. ";
2337 theDump += TCollection_AsciiString(count) + "\n";
2339 count = NbProblems->Value((Standard_Integer)BRepCheck_EmptyShell);
2341 theDump += " Empty Shell .............................. ";
2342 theDump += TCollection_AsciiString(count) + "\n";
2344 count = NbProblems->Value((Standard_Integer)BRepCheck_RedundantFace);
2346 theDump += " Redundant Face ........................... ";
2347 theDump += TCollection_AsciiString(count) + "\n";
2349 count = NbProblems->Value((Standard_Integer)BRepCheck_UnorientableShape);
2351 theDump += " Unorientable Shape ....................... ";
2352 theDump += TCollection_AsciiString(count) + "\n";
2354 count = NbProblems->Value((Standard_Integer)BRepCheck_NotClosed);
2356 theDump += " Not Closed ............................... ";
2357 theDump += TCollection_AsciiString(count) + "\n";
2359 count = NbProblems->Value((Standard_Integer)BRepCheck_NotConnected);
2361 theDump += " Not Connected ............................ ";
2362 theDump += TCollection_AsciiString(count) + "\n";
2364 count = NbProblems->Value((Standard_Integer)BRepCheck_SubshapeNotInShape);
2366 theDump += " Subshape not in Shape .................... ";
2367 theDump += TCollection_AsciiString(count) + "\n";
2369 count = NbProblems->Value((Standard_Integer)BRepCheck_BadOrientation);
2371 theDump += " Bad Orientation .......................... ";
2372 theDump += TCollection_AsciiString(count) + "\n";
2374 count = NbProblems->Value((Standard_Integer)BRepCheck_BadOrientationOfSubshape);
2376 theDump += " Bad Orientation of Subshape .............. ";
2377 theDump += TCollection_AsciiString(count) + "\n";
2379 count = NbProblems->Value((Standard_Integer)BRepCheck_CheckFail);
2381 theDump += " checkshape failure ....................... ";
2382 theDump += TCollection_AsciiString(count) + "\n";
2385 theDump += " ------------------------------------------------\n";
2386 theDump += "*** Shapes with problems : ";
2387 theDump += TCollection_AsciiString(sl->Length()) + "\n";
2389 Standard_Integer nbv, nbe, nbw, nbf, nbs, nbo;
2390 nbv = nbe = nbw = nbf = nbs = nbo = 0;
2392 for (i = 1; i <= sl->Length(); i++) {
2393 TopoDS_Shape shi = sl->Value(i);
2394 TopAbs_ShapeEnum sti = shi.ShapeType();
2396 case TopAbs_VERTEX : nbv++; break;
2397 case TopAbs_EDGE : nbe++; break;
2398 case TopAbs_WIRE : nbw++; break;
2399 case TopAbs_FACE : nbf++; break;
2400 case TopAbs_SHELL : nbs++; break;
2401 case TopAbs_SOLID : nbo++; break;
2407 theDump += "VERTEX : ";
2408 if (nbv < 10) theDump += " ";
2409 theDump += TCollection_AsciiString(nbv) + "\n";
2412 theDump += "EDGE : ";
2413 if (nbe < 10) theDump += " ";
2414 theDump += TCollection_AsciiString(nbe) + "\n";
2417 theDump += "WIRE : ";
2418 if (nbw < 10) theDump += " ";
2419 theDump += TCollection_AsciiString(nbw) + "\n";
2422 theDump += "FACE : ";
2423 if (nbf < 10) theDump += " ";
2424 theDump += TCollection_AsciiString(nbf) + "\n";
2427 theDump += "SHELL : ";
2428 if (nbs < 10) theDump += " ";
2429 theDump += TCollection_AsciiString(nbs) + "\n";
2432 theDump += "SOLID : ";
2433 if (nbo < 10) theDump += " ";
2434 theDump += TCollection_AsciiString(nbo) + "\n";
2439 //=======================================================================
2440 //function : GetProblemShapes
2441 // purpose : for StructuralDump
2442 //=======================================================================
2443 void GEOMImpl_IMeasureOperations::GetProblemShapes (const BRepCheck_Analyzer& theAna,
2444 const TopoDS_Shape& theShape,
2445 Handle(TopTools_HSequenceOfShape)& sl,
2446 Handle(TColStd_HArray1OfInteger)& NbProblems,
2447 TopTools_DataMapOfShapeListOfShape& theMap)
2449 for (TopoDS_Iterator iter(theShape); iter.More(); iter.Next()) {
2450 GetProblemShapes(theAna, iter.Value(), sl, NbProblems, theMap);
2452 TopAbs_ShapeEnum styp = theShape.ShapeType();
2453 BRepCheck_ListIteratorOfListOfStatus itl;
2454 TopTools_ListOfShape empty;
2455 if (!theMap.IsBound(theShape)) {
2456 theMap.Bind(theShape,empty);
2458 if (!theAna.Result(theShape).IsNull()) {
2459 itl.Initialize(theAna.Result(theShape)->Status());
2460 // !!! May be, we have to print all the problems, not only the first one ?
2461 if (itl.Value() != BRepCheck_NoError) {
2462 sl->Append(theShape);
2463 BRepCheck_Status stat = itl.Value();
2464 NbProblems->SetValue((Standard_Integer)stat,
2465 NbProblems->Value((Standard_Integer)stat) + 1);
2472 GetProblemSub(theAna, theShape, sl, NbProblems, TopAbs_VERTEX, theMap);
2475 GetProblemSub(theAna, theShape, sl, NbProblems, TopAbs_WIRE, theMap);
2476 GetProblemSub(theAna, theShape, sl, NbProblems, TopAbs_EDGE, theMap);
2477 GetProblemSub(theAna, theShape, sl, NbProblems, TopAbs_VERTEX, theMap);
2482 GetProblemSub(theAna, theShape, sl, NbProblems, TopAbs_SHELL, theMap);
2489 //=======================================================================
2490 //function : Contains
2491 //=======================================================================
2492 static Standard_Boolean Contains (const TopTools_ListOfShape& L,
2493 const TopoDS_Shape& S)
2495 TopTools_ListIteratorOfListOfShape it;
2496 for (it.Initialize(L); it.More(); it.Next()) {
2497 if (it.Value().IsSame(S)) {
2498 return Standard_True;
2501 return Standard_False;
2504 //=======================================================================
2505 //function : GetProblemSub
2506 // purpose : for StructuralDump
2507 //=======================================================================
2508 void GEOMImpl_IMeasureOperations::GetProblemSub (const BRepCheck_Analyzer& theAna,
2509 const TopoDS_Shape& theShape,
2510 Handle(TopTools_HSequenceOfShape)& sl,
2511 Handle(TColStd_HArray1OfInteger)& NbProblems,
2512 const TopAbs_ShapeEnum Subtype,
2513 TopTools_DataMapOfShapeListOfShape& theMap)
2515 BRepCheck_ListIteratorOfListOfStatus itl;
2516 TopExp_Explorer exp;
2517 for (exp.Init(theShape, Subtype); exp.More(); exp.Next()) {
2518 const TopoDS_Shape& sub = exp.Current();
2520 const Handle(BRepCheck_Result)& res = theAna.Result(sub);
2521 for (res->InitContextIterator();
2522 res->MoreShapeInContext();
2523 res->NextShapeInContext()) {
2524 if (res->ContextualShape().IsSame(theShape) && !Contains(theMap(sub), theShape)) {
2525 theMap(sub).Append(theShape);
2526 itl.Initialize(res->StatusOnShape());
2528 if (itl.Value() != BRepCheck_NoError) {
2529 Standard_Integer ii = 0;
2531 for (ii = 1; ii <= sl->Length(); ii++)
2532 if (sl->Value(ii).IsSame(sub)) break;
2534 if (ii > sl->Length()) {
2536 BRepCheck_Status stat = itl.Value();
2537 NbProblems->SetValue((Standard_Integer)stat,
2538 NbProblems->Value((Standard_Integer)stat) + 1);
2540 for (ii = 1; ii <= sl->Length(); ii++)
2541 if (sl->Value(ii).IsSame(theShape)) break;
2542 if (ii > sl->Length()) {
2543 sl->Append(theShape);
2544 BRepCheck_Status stat = itl.Value();
2545 NbProblems->SetValue((Standard_Integer)stat,
2546 NbProblems->Value((Standard_Integer)stat) + 1);