1 // Copyright (C) 2007-2010 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>
64 #include <TopoDS_Edge.hxx>
65 #include <TopoDS_Face.hxx>
66 #include <TopoDS_Shape.hxx>
67 #include <TopoDS_Vertex.hxx>
68 #include <TopoDS_Iterator.hxx>
69 #include <TopExp_Explorer.hxx>
70 #include <TopTools_MapOfShape.hxx>
71 #include <TopTools_ListOfShape.hxx>
72 #include <TopTools_ListIteratorOfListOfShape.hxx>
74 #include <GeomAbs_SurfaceType.hxx>
75 #include <Geom_Surface.hxx>
76 #include <Geom_Plane.hxx>
77 #include <Geom_SphericalSurface.hxx>
78 #include <Geom_CylindricalSurface.hxx>
79 #include <Geom_ToroidalSurface.hxx>
80 #include <Geom_ConicalSurface.hxx>
81 #include <Geom_SurfaceOfLinearExtrusion.hxx>
82 #include <Geom_SurfaceOfRevolution.hxx>
83 #include <Geom_BezierSurface.hxx>
84 #include <Geom_BSplineSurface.hxx>
85 #include <Geom_RectangularTrimmedSurface.hxx>
86 #include <Geom_OffsetSurface.hxx>
87 #include <Geom_Line.hxx>
92 #include <GeomAPI_ProjectPointOnCurve.hxx>
93 #include <GeomLProp_CLProps.hxx>
94 #include <GeomLProp_SLProps.hxx>
95 #include <ShapeAnalysis.hxx>
96 #include <ShapeAnalysis_Surface.hxx>
98 #include <Standard_Failure.hxx>
99 #include <Standard_ErrorHandler.hxx> // CAREFUL ! position of this file is critic : see Lucien PIGNOLONI / OCC
101 #include <BRepClass3d_SolidClassifier.hxx>
102 #include <BRep_Builder.hxx>
103 #include <GeomAPI_IntSS.hxx>
104 #include <Geom_Circle.hxx>
105 #include <Geom_SphericalSurface.hxx>
106 #include <Geom_ToroidalSurface.hxx>
107 #include <ShapeFix_Shape.hxx>
108 #include <TopoDS_Compound.hxx>
111 //=============================================================================
115 //=============================================================================
116 GEOMImpl_IMeasureOperations::GEOMImpl_IMeasureOperations (GEOM_Engine* theEngine, int theDocID)
117 : GEOM_IOperations(theEngine, theDocID)
119 MESSAGE("GEOMImpl_IMeasureOperations::GEOMImpl_IMeasureOperations");
122 //=============================================================================
126 //=============================================================================
127 GEOMImpl_IMeasureOperations::~GEOMImpl_IMeasureOperations()
129 MESSAGE("GEOMImpl_IMeasureOperations::~GEOMImpl_IMeasureOperations");
132 //=============================================================================
133 /*! Get kind and parameters of the given shape.
135 //=============================================================================
136 GEOMImpl_IMeasureOperations::ShapeKind GEOMImpl_IMeasureOperations::KindOfShape
137 (Handle(GEOM_Object) theShape,
138 Handle(TColStd_HSequenceOfInteger)& theIntegers,
139 Handle(TColStd_HSequenceOfReal)& theDoubles)
142 ShapeKind aKind = SK_NO_SHAPE;
144 if (theIntegers.IsNull()) theIntegers = new TColStd_HSequenceOfInteger;
145 else theIntegers->Clear();
147 if (theDoubles.IsNull()) theDoubles = new TColStd_HSequenceOfReal;
148 else theDoubles->Clear();
150 if (theShape.IsNull())
153 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
154 if (aRefShape.IsNull()) return aKind;
156 TopoDS_Shape aShape = aRefShape->GetValue();
157 if (aShape.IsNull()) return aKind;
159 int geom_type = theShape->GetType();
161 // check if it's advanced shape
162 if ( geom_type > ADVANCED_BASE ) {
168 GEOMAlgo_ShapeInfoFiller aSF;
169 aSF.SetShape(aShape);
171 Standard_Integer iErr = aSF.ErrorStatus();
173 SetErrorCode("Error in GEOMAlgo_ShapeInfoFiller");
176 const GEOMAlgo_ShapeInfo& anInfo = aSF.Info();
178 // Interprete results
179 TopAbs_ShapeEnum aType = anInfo.Type();
182 case TopAbs_COMPOUND:
183 case TopAbs_COMPSOLID:
185 // (+) geompy.kind.COMPOUND nb_solids nb_faces nb_edges nb_vertices
186 // (+) geompy.kind.COMPSOLID nb_solids nb_faces nb_edges nb_vertices
187 // ??? "nb_faces" - all faces or only 'standalone' faces?
188 if (aType == TopAbs_COMPOUND)
191 aKind = SK_COMPSOLID;
193 //theIntegers->Append(anInfo.NbSubShapes(TopAbs_COMPOUND));
194 //theIntegers->Append(anInfo.NbSubShapes(TopAbs_COMPSOLID));
195 theIntegers->Append(anInfo.NbSubShapes(TopAbs_SOLID));
196 theIntegers->Append(anInfo.NbSubShapes(TopAbs_FACE));
197 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
198 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
204 // (+) geompy.kind.SHELL geompy.info.closed nb_faces nb_edges nb_vertices
205 // (+) geompy.kind.SHELL geompy.info.unclosed nb_faces nb_edges nb_vertices
208 theIntegers->Append((int)anInfo.KindOfClosed());
210 theIntegers->Append(anInfo.NbSubShapes(TopAbs_FACE));
211 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
212 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
218 // (+) geompy.kind.WIRE geompy.info.closed nb_edges nb_vertices
219 // (+) geompy.kind.WIRE geompy.info.unclosed nb_edges nb_vertices
222 theIntegers->Append((int)anInfo.KindOfClosed());
224 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
225 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
233 GEOMAlgo_KindOfName aKN = anInfo.KindOfName();
236 case GEOMAlgo_KN_SPHERE:
237 // (+) geompy.kind.SPHERE xc yc zc R
241 gp_Pnt aC = anInfo.Location();
242 theDoubles->Append(aC.X());
243 theDoubles->Append(aC.Y());
244 theDoubles->Append(aC.Z());
246 theDoubles->Append(anInfo.Radius1());
249 case GEOMAlgo_KN_CYLINDER:
250 // (+) geompy.kind.CYLINDER xb yb zb dx dy dz R H
254 gp_Pnt aC = anInfo.Location();
255 theDoubles->Append(aC.X());
256 theDoubles->Append(aC.Y());
257 theDoubles->Append(aC.Z());
259 gp_Ax3 anAx3 = anInfo.Position();
260 gp_Dir aD = anAx3.Direction();
261 theDoubles->Append(aD.X());
262 theDoubles->Append(aD.Y());
263 theDoubles->Append(aD.Z());
265 theDoubles->Append(anInfo.Radius1());
266 theDoubles->Append(anInfo.Height());
269 case GEOMAlgo_KN_BOX:
270 // (+) geompy.kind.BOX xc yc zc ax ay az
274 gp_Pnt aC = anInfo.Location();
275 theDoubles->Append(aC.X());
276 theDoubles->Append(aC.Y());
277 theDoubles->Append(aC.Z());
279 gp_Ax3 anAx3 = anInfo.Position();
280 gp_Dir aD = anAx3.Direction();
281 gp_Dir aX = anAx3.XDirection();
284 if (aD.IsParallel(gp::DZ(), Precision::Angular()) &&
285 aX.IsParallel(gp::DX(), Precision::Angular())) {
286 theDoubles->Append(anInfo.Length()); // ax'
287 theDoubles->Append(anInfo.Width()); // ay'
288 theDoubles->Append(anInfo.Height()); // az'
290 else if (aD.IsParallel(gp::DZ(), Precision::Angular()) &&
291 aX.IsParallel(gp::DY(), Precision::Angular())) {
292 theDoubles->Append(anInfo.Width()); // ay'
293 theDoubles->Append(anInfo.Length()); // ax'
294 theDoubles->Append(anInfo.Height()); // az'
296 else if (aD.IsParallel(gp::DX(), Precision::Angular()) &&
297 aX.IsParallel(gp::DZ(), Precision::Angular())) {
298 theDoubles->Append(anInfo.Height()); // az'
299 theDoubles->Append(anInfo.Width()); // ay'
300 theDoubles->Append(anInfo.Length()); // ax'
302 else if (aD.IsParallel(gp::DX(), Precision::Angular()) &&
303 aX.IsParallel(gp::DY(), Precision::Angular())) {
304 theDoubles->Append(anInfo.Height()); // az'
305 theDoubles->Append(anInfo.Length()); // ax'
306 theDoubles->Append(anInfo.Width()); // ay'
308 else if (aD.IsParallel(gp::DY(), Precision::Angular()) &&
309 aX.IsParallel(gp::DZ(), Precision::Angular())) {
310 theDoubles->Append(anInfo.Width()); // ay'
311 theDoubles->Append(anInfo.Height()); // az'
312 theDoubles->Append(anInfo.Length()); // ax'
314 else if (aD.IsParallel(gp::DY(), Precision::Angular()) &&
315 aX.IsParallel(gp::DX(), Precision::Angular())) {
316 theDoubles->Append(anInfo.Length()); // ax'
317 theDoubles->Append(anInfo.Height()); // az'
318 theDoubles->Append(anInfo.Width()); // ay'
321 // (+) geompy.kind.ROTATED_BOX xo yo zo zx zy zz xx xy xz ax ay az
322 aKind = SK_ROTATED_BOX;
324 // Direction and XDirection
325 theDoubles->Append(aD.X());
326 theDoubles->Append(aD.Y());
327 theDoubles->Append(aD.Z());
329 theDoubles->Append(aX.X());
330 theDoubles->Append(aX.Y());
331 theDoubles->Append(aX.Z());
334 theDoubles->Append(anInfo.Length());
335 theDoubles->Append(anInfo.Width());
336 theDoubles->Append(anInfo.Height());
340 case GEOMAlgo_KN_TORUS:
341 // (+) geompy.kind.TORUS xc yc zc dx dy dz R_1 R_2
345 gp_Pnt aO = anInfo.Location();
346 theDoubles->Append(aO.X());
347 theDoubles->Append(aO.Y());
348 theDoubles->Append(aO.Z());
350 gp_Ax3 anAx3 = anInfo.Position();
351 gp_Dir aD = anAx3.Direction();
352 theDoubles->Append(aD.X());
353 theDoubles->Append(aD.Y());
354 theDoubles->Append(aD.Z());
356 theDoubles->Append(anInfo.Radius1());
357 theDoubles->Append(anInfo.Radius2());
360 case GEOMAlgo_KN_CONE:
361 // (+) geompy.kind.CONE xb yb zb dx dy dz R_1 R_2 H
365 gp_Pnt aO = anInfo.Location();
366 theDoubles->Append(aO.X());
367 theDoubles->Append(aO.Y());
368 theDoubles->Append(aO.Z());
370 gp_Ax3 anAx3 = anInfo.Position();
371 gp_Dir aD = anAx3.Direction();
372 theDoubles->Append(aD.X());
373 theDoubles->Append(aD.Y());
374 theDoubles->Append(aD.Z());
376 theDoubles->Append(anInfo.Radius1());
377 theDoubles->Append(anInfo.Radius2());
378 theDoubles->Append(anInfo.Height());
381 case GEOMAlgo_KN_POLYHEDRON:
382 // (+) geompy.kind.POLYHEDRON nb_faces nb_edges nb_vertices
384 aKind = SK_POLYHEDRON;
386 theIntegers->Append(anInfo.NbSubShapes(TopAbs_FACE));
387 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
388 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
392 // (+) geompy.kind.SOLID nb_faces nb_edges nb_vertices
394 theIntegers->Append(anInfo.NbSubShapes(TopAbs_FACE));
395 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
396 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
406 GEOMAlgo_KindOfName aKN = anInfo.KindOfName();
408 case GEOMAlgo_KN_SPHERE:
409 // (+) geompy.kind.SPHERE2D xc yc zc R
413 gp_Pnt aC = anInfo.Location();
414 theDoubles->Append(aC.X());
415 theDoubles->Append(aC.Y());
416 theDoubles->Append(aC.Z());
418 theDoubles->Append(anInfo.Radius1());
421 case GEOMAlgo_KN_CYLINDER:
422 // (+) geompy.kind.CYLINDER2D xb yb zb dx dy dz R H
424 aKind = SK_CYLINDER2D;
426 gp_Pnt aO = anInfo.Location();
427 theDoubles->Append(aO.X());
428 theDoubles->Append(aO.Y());
429 theDoubles->Append(aO.Z());
431 gp_Ax3 anAx3 = anInfo.Position();
432 gp_Dir aD = anAx3.Direction();
433 theDoubles->Append(aD.X());
434 theDoubles->Append(aD.Y());
435 theDoubles->Append(aD.Z());
437 theDoubles->Append(anInfo.Radius1());
438 theDoubles->Append(anInfo.Height());
441 case GEOMAlgo_KN_TORUS:
442 // (+) geompy.kind.TORUS2D xc yc zc dx dy dz R_1 R_2
446 gp_Pnt aO = anInfo.Location();
447 theDoubles->Append(aO.X());
448 theDoubles->Append(aO.Y());
449 theDoubles->Append(aO.Z());
451 gp_Ax3 anAx3 = anInfo.Position();
452 gp_Dir aD = anAx3.Direction();
453 theDoubles->Append(aD.X());
454 theDoubles->Append(aD.Y());
455 theDoubles->Append(aD.Z());
457 theDoubles->Append(anInfo.Radius1());
458 theDoubles->Append(anInfo.Radius2());
461 case GEOMAlgo_KN_CONE:
462 // (+) geompy.kind.CONE2D xc yc zc dx dy dz R_1 R_2 H
466 gp_Pnt aO = anInfo.Location();
467 theDoubles->Append(aO.X());
468 theDoubles->Append(aO.Y());
469 theDoubles->Append(aO.Z());
471 gp_Ax3 anAx3 = anInfo.Position();
472 gp_Dir aD = anAx3.Direction();
473 theDoubles->Append(aD.X());
474 theDoubles->Append(aD.Y());
475 theDoubles->Append(aD.Z());
477 theDoubles->Append(anInfo.Radius1());
478 theDoubles->Append(anInfo.Radius2());
479 theDoubles->Append(anInfo.Height());
482 case GEOMAlgo_KN_DISKCIRCLE:
483 // (+) geompy.kind.DISK_CIRCLE xc yc zc dx dy dz R
485 aKind = SK_DISK_CIRCLE;
487 gp_Pnt aC = anInfo.Location();
488 theDoubles->Append(aC.X());
489 theDoubles->Append(aC.Y());
490 theDoubles->Append(aC.Z());
492 gp_Ax3 anAx3 = anInfo.Position();
493 gp_Dir aD = anAx3.Direction();
494 theDoubles->Append(aD.X());
495 theDoubles->Append(aD.Y());
496 theDoubles->Append(aD.Z());
498 theDoubles->Append(anInfo.Radius1());
501 case GEOMAlgo_KN_DISKELLIPSE:
502 // (+) geompy.kind.DISK_ELLIPSE xc yc zc dx dy dz R_1 R_2
504 aKind = SK_DISK_ELLIPSE;
506 gp_Pnt aC = anInfo.Location();
507 theDoubles->Append(aC.X());
508 theDoubles->Append(aC.Y());
509 theDoubles->Append(aC.Z());
511 gp_Ax3 anAx3 = anInfo.Position();
512 gp_Dir aD = anAx3.Direction();
513 theDoubles->Append(aD.X());
514 theDoubles->Append(aD.Y());
515 theDoubles->Append(aD.Z());
517 theDoubles->Append(anInfo.Radius1());
518 theDoubles->Append(anInfo.Radius2());
521 case GEOMAlgo_KN_RECTANGLE:
522 case GEOMAlgo_KN_TRIANGLE:
523 case GEOMAlgo_KN_QUADRANGLE:
524 case GEOMAlgo_KN_POLYGON:
525 // (+) geompy.kind.POLYGON xo yo zo dx dy dz nb_edges nb_vertices
529 gp_Pnt aO = anInfo.Location();
530 theDoubles->Append(aO.X());
531 theDoubles->Append(aO.Y());
532 theDoubles->Append(aO.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 case GEOMAlgo_KN_PLANE: // infinite
545 // (+) geompy.kind.PLANE xo yo zo dx dy dz
549 gp_Pnt aC = anInfo.Location();
550 theDoubles->Append(aC.X());
551 theDoubles->Append(aC.Y());
552 theDoubles->Append(aC.Z());
554 gp_Ax3 anAx3 = anInfo.Position();
555 gp_Dir aD = anAx3.Direction();
556 theDoubles->Append(aD.X());
557 theDoubles->Append(aD.Y());
558 theDoubles->Append(aD.Z());
562 if (anInfo.KindOfShape() == GEOMAlgo_KS_PLANE) {
563 // (+) geompy.kind.PLANAR xo yo zo dx dy dz nb_edges nb_vertices
567 gp_Pnt aC = anInfo.Location();
568 theDoubles->Append(aC.X());
569 theDoubles->Append(aC.Y());
570 theDoubles->Append(aC.Z());
572 gp_Ax3 anAx3 = anInfo.Position();
573 gp_Dir aD = anAx3.Direction();
574 theDoubles->Append(aD.X());
575 theDoubles->Append(aD.Y());
576 theDoubles->Append(aD.Z());
578 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
579 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
582 // ??? geompy.kind.FACE nb_edges nb_vertices _surface_type_id_
583 // (+) geompy.kind.FACE nb_edges nb_vertices
585 theIntegers->Append(anInfo.NbSubShapes(TopAbs_EDGE));
586 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
596 GEOMAlgo_KindOfName aKN = anInfo.KindOfName();
598 case GEOMAlgo_KN_CIRCLE:
600 // (+) geompy.kind.CIRCLE xc yc zc dx dy dz R
603 gp_Pnt aC = anInfo.Location();
604 theDoubles->Append(aC.X());
605 theDoubles->Append(aC.Y());
606 theDoubles->Append(aC.Z());
608 gp_Ax3 anAx3 = anInfo.Position();
609 gp_Dir aD = anAx3.Direction();
610 theDoubles->Append(aD.X());
611 theDoubles->Append(aD.Y());
612 theDoubles->Append(aD.Z());
614 theDoubles->Append(anInfo.Radius1());
617 case GEOMAlgo_KN_ARCCIRCLE:
619 // (+) geompy.kind.ARC_CIRCLE xc yc zc dx dy dz R x1 y1 z1 x2 y2 z2
620 aKind = SK_ARC_CIRCLE;
622 gp_Pnt aC = anInfo.Location();
623 theDoubles->Append(aC.X());
624 theDoubles->Append(aC.Y());
625 theDoubles->Append(aC.Z());
627 gp_Ax3 anAx3 = anInfo.Position();
628 gp_Dir aD = anAx3.Direction();
629 theDoubles->Append(aD.X());
630 theDoubles->Append(aD.Y());
631 theDoubles->Append(aD.Z());
633 theDoubles->Append(anInfo.Radius1());
635 gp_Pnt aP1 = anInfo.Pnt1();
636 theDoubles->Append(aP1.X());
637 theDoubles->Append(aP1.Y());
638 theDoubles->Append(aP1.Z());
640 gp_Pnt aP2 = anInfo.Pnt2();
641 theDoubles->Append(aP2.X());
642 theDoubles->Append(aP2.Y());
643 theDoubles->Append(aP2.Z());
646 case GEOMAlgo_KN_ELLIPSE:
648 // (+) geompy.kind.ELLIPSE xc yc zc dx dy dz R_1 R_2
651 gp_Pnt aC = anInfo.Location();
652 theDoubles->Append(aC.X());
653 theDoubles->Append(aC.Y());
654 theDoubles->Append(aC.Z());
656 gp_Ax3 anAx3 = anInfo.Position();
657 gp_Dir aD = anAx3.Direction();
658 theDoubles->Append(aD.X());
659 theDoubles->Append(aD.Y());
660 theDoubles->Append(aD.Z());
662 theDoubles->Append(anInfo.Radius1());
663 theDoubles->Append(anInfo.Radius2());
666 case GEOMAlgo_KN_ARCELLIPSE:
668 // (+) geompy.kind.ARC_ELLIPSE xc yc zc dx dy dz R_1 R_2 x1 y1 z1 x2 y2 z2
669 aKind = SK_ARC_ELLIPSE;
671 gp_Pnt aC = anInfo.Location();
672 theDoubles->Append(aC.X());
673 theDoubles->Append(aC.Y());
674 theDoubles->Append(aC.Z());
676 gp_Ax3 anAx3 = anInfo.Position();
677 gp_Dir aD = anAx3.Direction();
678 theDoubles->Append(aD.X());
679 theDoubles->Append(aD.Y());
680 theDoubles->Append(aD.Z());
682 theDoubles->Append(anInfo.Radius1());
683 theDoubles->Append(anInfo.Radius2());
685 gp_Pnt aP1 = anInfo.Pnt1();
686 theDoubles->Append(aP1.X());
687 theDoubles->Append(aP1.Y());
688 theDoubles->Append(aP1.Z());
690 gp_Pnt aP2 = anInfo.Pnt2();
691 theDoubles->Append(aP2.X());
692 theDoubles->Append(aP2.Y());
693 theDoubles->Append(aP2.Z());
696 case GEOMAlgo_KN_LINE:
698 // ??? geompy.kind.LINE x1 y1 z1 x2 y2 z2
699 // (+) geompy.kind.LINE x1 y1 z1 dx dy dz
702 gp_Pnt aO = anInfo.Location();
703 theDoubles->Append(aO.X());
704 theDoubles->Append(aO.Y());
705 theDoubles->Append(aO.Z());
707 gp_Dir aD = anInfo.Direction();
708 theDoubles->Append(aD.X());
709 theDoubles->Append(aD.Y());
710 theDoubles->Append(aD.Z());
713 case GEOMAlgo_KN_SEGMENT:
715 // (+) geompy.kind.SEGMENT x1 y1 z1 x2 y2 z2
718 gp_Pnt aP1 = anInfo.Pnt1();
719 theDoubles->Append(aP1.X());
720 theDoubles->Append(aP1.Y());
721 theDoubles->Append(aP1.Z());
723 gp_Pnt aP2 = anInfo.Pnt2();
724 theDoubles->Append(aP2.X());
725 theDoubles->Append(aP2.Y());
726 theDoubles->Append(aP2.Z());
730 // ??? geompy.kind.EDGE nb_vertices _curve_type_id_
731 // (+) geompy.kind.EDGE nb_vertices
732 theIntegers->Append(anInfo.NbSubShapes(TopAbs_VERTEX));
739 // (+) geompy.kind.VERTEX x y z
742 gp_Pnt aP = anInfo.Location();
743 theDoubles->Append(aP.X());
744 theDoubles->Append(aP.Y());
745 theDoubles->Append(aP.Z());
754 //=============================================================================
755 /*! Get LCS, corresponding to the given shape.
756 * Origin of the LCS is situated at the shape's center of mass.
757 * Axes of the LCS are obtained from shape's location or,
758 * if the shape is a planar face, from position of its plane.
760 //=============================================================================
761 gp_Ax3 GEOMImpl_IMeasureOperations::GetPosition (const TopoDS_Shape& theShape)
765 if (theShape.IsNull())
769 aResult.Transform(theShape.Location().Transformation());
770 if (theShape.ShapeType() == TopAbs_FACE) {
771 Handle(Geom_Surface) aGS = BRep_Tool::Surface(TopoDS::Face(theShape));
772 if (!aGS.IsNull() && aGS->IsKind(STANDARD_TYPE(Geom_Plane))) {
773 Handle(Geom_Plane) aGPlane = Handle(Geom_Plane)::DownCast(aGS);
774 gp_Pln aPln = aGPlane->Pln();
775 aResult = aPln.Position();
781 if (theShape.ShapeType() == TopAbs_VERTEX) {
782 aPnt = BRep_Tool::Pnt(TopoDS::Vertex(theShape));
785 GProp_GProps aSystem;
786 if (theShape.ShapeType() == TopAbs_EDGE || theShape.ShapeType() == TopAbs_WIRE)
787 BRepGProp::LinearProperties(theShape, aSystem);
788 else if (theShape.ShapeType() == TopAbs_FACE || theShape.ShapeType() == TopAbs_SHELL)
789 BRepGProp::SurfaceProperties(theShape, aSystem);
791 BRepGProp::VolumeProperties(theShape, aSystem);
793 aPnt = aSystem.CentreOfMass();
796 aResult.SetLocation(aPnt);
801 //=============================================================================
805 //=============================================================================
806 void GEOMImpl_IMeasureOperations::GetPosition
807 (Handle(GEOM_Object) theShape,
808 Standard_Real& Ox, Standard_Real& Oy, Standard_Real& Oz,
809 Standard_Real& Zx, Standard_Real& Zy, Standard_Real& Zz,
810 Standard_Real& Xx, Standard_Real& Xy, Standard_Real& Xz)
814 //Set default values: global CS
815 Ox = Oy = Oz = Zx = Zy = Xy = Xz = 0.;
818 if (theShape.IsNull()) return;
820 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
821 if (aRefShape.IsNull()) return;
823 TopoDS_Shape aShape = aRefShape->GetValue();
824 if (aShape.IsNull()) {
825 SetErrorCode("The Objects has NULL Shape");
830 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
834 gp_Ax3 anAx3 = GetPosition(aShape);
836 gp_Pnt anOri = anAx3.Location();
837 gp_Dir aDirZ = anAx3.Direction();
838 gp_Dir aDirX = anAx3.XDirection();
841 anOri.Coord(Ox, Oy, Oz);
842 aDirZ.Coord(Zx, Zy, Zz);
843 aDirX.Coord(Xx, Xy, Xz);
845 catch (Standard_Failure) {
846 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
847 SetErrorCode(aFail->GetMessageString());
854 //=============================================================================
858 //=============================================================================
859 Handle(GEOM_Object) GEOMImpl_IMeasureOperations::GetCentreOfMass
860 (Handle(GEOM_Object) theShape)
864 if (theShape.IsNull()) return NULL;
866 //Add a new CentreOfMass object
867 Handle(GEOM_Object) aCDG = GetEngine()->AddObject(GetDocID(), GEOM_CDG);
869 //Add a new CentreOfMass function
870 Handle(GEOM_Function) aFunction =
871 aCDG->AddFunction(GEOMImpl_MeasureDriver::GetID(), CDG_MEASURE);
872 if (aFunction.IsNull()) return NULL;
874 //Check if the function is set correctly
875 if (aFunction->GetDriverGUID() != GEOMImpl_MeasureDriver::GetID()) return NULL;
877 GEOMImpl_IMeasure aCI (aFunction);
879 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
880 if (aRefShape.IsNull()) return NULL;
882 aCI.SetBase(aRefShape);
884 //Compute the CentreOfMass value
886 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
889 if (!GetSolver()->ComputeFunction(aFunction)) {
890 SetErrorCode("Measure driver failed to compute centre of mass");
894 catch (Standard_Failure) {
895 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
896 SetErrorCode(aFail->GetMessageString());
900 //Make a Python command
901 GEOM::TPythonDump(aFunction) << aCDG << " = geompy.MakeCDG(" << theShape << ")";
907 //=============================================================================
911 //=============================================================================
912 Handle(GEOM_Object) GEOMImpl_IMeasureOperations::GetVertexByIndex
913 (Handle(GEOM_Object) theShape,
914 Standard_Integer theIndex)
918 if (theShape.IsNull()) return NULL;
920 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
921 if (aRefShape.IsNull()) return NULL;
923 //Add a new Vertex object
924 Handle(GEOM_Object) aVertex = GetEngine()->AddObject(GetDocID(), GEOM_POINT);
927 Handle(GEOM_Function) aFunction =
928 aVertex->AddFunction(GEOMImpl_MeasureDriver::GetID(), VERTEX_BY_INDEX);
929 if (aFunction.IsNull()) return NULL;
931 //Check if the function is set correctly
932 if (aFunction->GetDriverGUID() != GEOMImpl_MeasureDriver::GetID()) return NULL;
934 GEOMImpl_IMeasure aCI (aFunction);
935 aCI.SetBase(aRefShape);
936 aCI.SetIndex(theIndex);
940 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
943 if (!GetSolver()->ComputeFunction(aFunction)) {
944 SetErrorCode("Vertex by index driver failed.");
948 catch (Standard_Failure) {
949 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
950 SetErrorCode(aFail->GetMessageString());
954 //Make a Python command
955 GEOM::TPythonDump(aFunction) << aVertex << " = geompy.GetVertexByIndex(" << theShape << ", " << theIndex << ")";
961 //=============================================================================
965 //=============================================================================
966 Handle(GEOM_Object) GEOMImpl_IMeasureOperations::GetNormal
967 (Handle(GEOM_Object) theFace,
968 Handle(GEOM_Object) theOptionalPoint)
972 if (theFace.IsNull()) return NULL;
974 //Add a new Normale object
975 Handle(GEOM_Object) aNorm = GetEngine()->AddObject(GetDocID(), GEOM_VECTOR);
977 //Add a new Normale function
978 Handle(GEOM_Function) aFunction =
979 aNorm->AddFunction(GEOMImpl_MeasureDriver::GetID(), VECTOR_FACE_NORMALE);
980 if (aFunction.IsNull()) return NULL;
982 //Check if the function is set correctly
983 if (aFunction->GetDriverGUID() != GEOMImpl_MeasureDriver::GetID()) return NULL;
985 GEOMImpl_IMeasure aCI (aFunction);
987 Handle(GEOM_Function) aFace = theFace->GetLastFunction();
988 if (aFace.IsNull()) return NULL;
992 if (!theOptionalPoint.IsNull()) {
993 Handle(GEOM_Function) anOptPnt = theOptionalPoint->GetLastFunction();
994 aCI.SetPoint(anOptPnt);
997 //Compute the Normale value
999 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
1002 if (!GetSolver()->ComputeFunction(aFunction)) {
1003 SetErrorCode("Measure driver failed to compute normake of face");
1007 catch (Standard_Failure) {
1008 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1009 SetErrorCode(aFail->GetMessageString());
1013 //Make a Python command
1014 GEOM::TPythonDump pd (aFunction);
1015 pd << aNorm << " = geompy.GetNormal(" << theFace;
1016 if (!theOptionalPoint.IsNull()) {
1017 pd << ", " << theOptionalPoint;
1025 //=============================================================================
1027 * GetBasicProperties
1029 //=============================================================================
1030 void GEOMImpl_IMeasureOperations::GetBasicProperties (Handle(GEOM_Object) theShape,
1031 Standard_Real& theLength,
1032 Standard_Real& theSurfArea,
1033 Standard_Real& theVolume)
1037 if (theShape.IsNull()) return;
1039 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1040 if (aRefShape.IsNull()) return;
1042 TopoDS_Shape aShape = aRefShape->GetValue();
1043 if (aShape.IsNull()) {
1044 SetErrorCode("The Objects has NULL Shape");
1048 //Compute the parameters
1049 GProp_GProps LProps, SProps;
1051 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
1054 BRepGProp::LinearProperties(aShape, LProps);
1055 theLength = LProps.Mass();
1057 BRepGProp::SurfaceProperties(aShape, SProps);
1058 theSurfArea = SProps.Mass();
1061 if (aShape.ShapeType() < TopAbs_SHELL) {
1062 for (TopExp_Explorer Exp (aShape, TopAbs_SOLID); Exp.More(); Exp.Next()) {
1063 GProp_GProps VProps;
1064 BRepGProp::VolumeProperties(Exp.Current(), VProps);
1065 theVolume += VProps.Mass();
1069 catch (Standard_Failure) {
1070 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1071 SetErrorCode(aFail->GetMessageString());
1078 //=============================================================================
1082 //=============================================================================
1083 void GEOMImpl_IMeasureOperations::GetInertia
1084 (Handle(GEOM_Object) theShape,
1085 Standard_Real& I11, Standard_Real& I12, Standard_Real& I13,
1086 Standard_Real& I21, Standard_Real& I22, Standard_Real& I23,
1087 Standard_Real& I31, Standard_Real& I32, Standard_Real& I33,
1088 Standard_Real& Ix , Standard_Real& Iy , Standard_Real& Iz)
1092 if (theShape.IsNull()) return;
1094 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1095 if (aRefShape.IsNull()) return;
1097 TopoDS_Shape aShape = aRefShape->GetValue();
1098 if (aShape.IsNull()) {
1099 SetErrorCode("The Objects has NULL Shape");
1103 //Compute the parameters
1104 GProp_GProps System;
1107 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
1110 if (aShape.ShapeType() == TopAbs_VERTEX ||
1111 aShape.ShapeType() == TopAbs_EDGE ||
1112 aShape.ShapeType() == TopAbs_WIRE) {
1113 BRepGProp::LinearProperties(aShape, System);
1114 } else if (aShape.ShapeType() == TopAbs_FACE ||
1115 aShape.ShapeType() == TopAbs_SHELL) {
1116 BRepGProp::SurfaceProperties(aShape, System);
1118 BRepGProp::VolumeProperties(aShape, System);
1120 gp_Mat I = System.MatrixOfInertia();
1134 GProp_PrincipalProps Pr = System.PrincipalProperties();
1135 Pr.Moments(Ix,Iy,Iz);
1137 catch (Standard_Failure) {
1138 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1139 SetErrorCode(aFail->GetMessageString());
1146 //=============================================================================
1150 //=============================================================================
1151 void GEOMImpl_IMeasureOperations::GetBoundingBox
1152 (Handle(GEOM_Object) theShape,
1153 Standard_Real& Xmin, Standard_Real& Xmax,
1154 Standard_Real& Ymin, Standard_Real& Ymax,
1155 Standard_Real& Zmin, Standard_Real& Zmax)
1159 if (theShape.IsNull()) return;
1161 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1162 if (aRefShape.IsNull()) return;
1164 TopoDS_Shape aShape = aRefShape->GetValue();
1165 if (aShape.IsNull()) {
1166 SetErrorCode("The Objects has NULL Shape");
1170 //Compute the parameters
1174 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
1177 BRepBndLib::Add(aShape, B);
1178 B.Get(Xmin, Ymin, Zmin, Xmax, Ymax, Zmax);
1180 catch (Standard_Failure) {
1181 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1182 SetErrorCode(aFail->GetMessageString());
1189 //=============================================================================
1193 //=============================================================================
1194 void GEOMImpl_IMeasureOperations::GetTolerance
1195 (Handle(GEOM_Object) theShape,
1196 Standard_Real& FaceMin, Standard_Real& FaceMax,
1197 Standard_Real& EdgeMin, Standard_Real& EdgeMax,
1198 Standard_Real& VertMin, Standard_Real& VertMax)
1202 if (theShape.IsNull()) return;
1204 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1205 if (aRefShape.IsNull()) return;
1207 TopoDS_Shape aShape = aRefShape->GetValue();
1208 if (aShape.IsNull()) {
1209 SetErrorCode("The Objects has NULL Shape");
1213 //Compute the parameters
1215 FaceMin = EdgeMin = VertMin = RealLast();
1216 FaceMax = EdgeMax = VertMax = -RealLast();
1219 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
1222 for (TopExp_Explorer ExF (aShape, TopAbs_FACE); ExF.More(); ExF.Next()) {
1223 TopoDS_Face Face = TopoDS::Face(ExF.Current());
1224 T = BRep_Tool::Tolerance(Face);
1230 for (TopExp_Explorer ExE (aShape, TopAbs_EDGE); ExE.More(); ExE.Next()) {
1231 TopoDS_Edge Edge = TopoDS::Edge(ExE.Current());
1232 T = BRep_Tool::Tolerance(Edge);
1238 for (TopExp_Explorer ExV (aShape, TopAbs_VERTEX); ExV.More(); ExV.Next()) {
1239 TopoDS_Vertex Vertex = TopoDS::Vertex(ExV.Current());
1240 T = BRep_Tool::Tolerance(Vertex);
1247 catch (Standard_Failure) {
1248 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1249 SetErrorCode(aFail->GetMessageString());
1256 //=============================================================================
1260 //=============================================================================
1261 bool GEOMImpl_IMeasureOperations::CheckShape (Handle(GEOM_Object) theShape,
1262 const Standard_Boolean theIsCheckGeom,
1263 TCollection_AsciiString& theDump)
1267 if (theShape.IsNull()) return false;
1269 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1270 if (aRefShape.IsNull()) return false;
1272 TopoDS_Shape aShape = aRefShape->GetValue();
1273 if (aShape.IsNull()) {
1274 SetErrorCode("The Objects has NULL Shape");
1278 //Compute the parameters
1279 bool isValid = false;
1281 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
1284 BRepCheck_Analyzer ana (aShape, theIsCheckGeom);
1285 if (ana.IsValid()) {
1289 StructuralDump(ana, aShape, theDump);
1292 catch (Standard_Failure) {
1293 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1294 SetErrorCode(aFail->GetMessageString());
1302 //=============================================================================
1306 //=============================================================================
1307 TCollection_AsciiString GEOMImpl_IMeasureOperations::IsGoodForSolid (Handle(GEOM_Object) theShape)
1311 TCollection_AsciiString aRes = "";
1313 if (theShape.IsNull()) {
1314 aRes = "WRN_NULL_OBJECT_OR_SHAPE";
1317 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1318 if (aRefShape.IsNull()) {
1319 aRes = "WRN_NULL_OBJECT_OR_SHAPE";
1322 TopoDS_Shape aShape = aRefShape->GetValue();
1323 if (aShape.IsNull()) {
1324 aRes = "WRN_NULL_OBJECT_OR_SHAPE";
1327 if (aShape.ShapeType() == TopAbs_COMPOUND) {
1328 TopoDS_Iterator It (aShape, Standard_True, Standard_True);
1329 if (It.More()) aShape = It.Value();
1331 if (aShape.ShapeType() == TopAbs_SHELL) {
1332 if (!aShape.Closed()) {
1333 aRes = "WRN_SHAPE_UNCLOSED";
1337 aRes = "WRN_SHAPE_NOT_SHELL";
1349 //=============================================================================
1353 //=============================================================================
1354 TCollection_AsciiString GEOMImpl_IMeasureOperations::WhatIs (Handle(GEOM_Object) theShape)
1358 TCollection_AsciiString Astr;
1360 if (theShape.IsNull()) return Astr;
1362 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1363 if (aRefShape.IsNull()) return Astr;
1365 TopoDS_Shape aShape = aRefShape->GetValue();
1366 if (aShape.IsNull()) {
1367 SetErrorCode("The Objects has NULL Shape");
1371 //Compute the parameters
1372 if (aShape.ShapeType() == TopAbs_EDGE) {
1373 if (BRep_Tool::Degenerated(TopoDS::Edge(aShape))) {
1374 Astr = Astr + " It is a degenerated edge \n";
1378 Astr = Astr + " Number of sub-shapes : \n";
1381 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
1384 int iType, nbTypes [TopAbs_SHAPE];
1385 for (iType = 0; iType < TopAbs_SHAPE; ++iType)
1387 nbTypes[aShape.ShapeType()]++;
1389 TopTools_MapOfShape aMapOfShape;
1390 aMapOfShape.Add(aShape);
1391 TopTools_ListOfShape aListOfShape;
1392 aListOfShape.Append(aShape);
1394 TopTools_ListIteratorOfListOfShape itL (aListOfShape);
1395 for (; itL.More(); itL.Next()) {
1396 TopoDS_Iterator it (itL.Value());
1397 for (; it.More(); it.Next()) {
1398 TopoDS_Shape s = it.Value();
1399 if (aMapOfShape.Add(s)) {
1400 aListOfShape.Append(s);
1401 nbTypes[s.ShapeType()]++;
1406 Astr = Astr + " VERTEX : " + TCollection_AsciiString(nbTypes[TopAbs_VERTEX]) + "\n";
1407 Astr = Astr + " EDGE : " + TCollection_AsciiString(nbTypes[TopAbs_EDGE]) + "\n";
1408 Astr = Astr + " WIRE : " + TCollection_AsciiString(nbTypes[TopAbs_WIRE]) + "\n";
1409 Astr = Astr + " FACE : " + TCollection_AsciiString(nbTypes[TopAbs_FACE]) + "\n";
1410 Astr = Astr + " SHELL : " + TCollection_AsciiString(nbTypes[TopAbs_SHELL]) + "\n";
1411 Astr = Astr + " SOLID : " + TCollection_AsciiString(nbTypes[TopAbs_SOLID]) + "\n";
1412 Astr = Astr + " COMPSOLID : " + TCollection_AsciiString(nbTypes[TopAbs_COMPSOLID]) + "\n";
1413 Astr = Astr + " COMPOUND : " + TCollection_AsciiString(nbTypes[TopAbs_COMPOUND]) + "\n";
1414 Astr = Astr + " SHAPE : " + TCollection_AsciiString(aMapOfShape.Extent());
1416 catch (Standard_Failure) {
1417 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1418 SetErrorCode(aFail->GetMessageString());
1427 //=======================================================================
1428 //function : CheckSingularCase
1429 //purpose : auxilary for GetMinDistance()
1430 // workaround for bugs 19899, 19908 and 19910 from Mantis
1431 //=======================================================================
1432 static double CheckSingularCase(const TopoDS_Shape& aSh1,
1433 const TopoDS_Shape& aSh2,
1434 gp_Pnt& Ptmp1, gp_Pnt& Ptmp2)
1436 bool IsChange1 = false;
1437 double AddDist1 = 0.0;
1438 TopExp_Explorer anExp;
1439 TopoDS_Shape tmpSh1, tmpSh2;
1441 for ( anExp.Init( aSh1, TopAbs_FACE ); anExp.More(); anExp.Next() ) {
1443 tmpSh1 = anExp.Current();
1446 TopoDS_Shape sh = aSh1;
1447 while(sh.ShapeType()==TopAbs_COMPOUND) {
1448 TopoDS_Iterator it(sh);
1451 Handle(Geom_Surface) S = BRep_Tool::Surface(TopoDS::Face(tmpSh1));
1452 if( S->IsKind(STANDARD_TYPE(Geom_SphericalSurface)) ||
1453 S->IsKind(STANDARD_TYPE(Geom_ToroidalSurface)) ) {
1454 if( sh.ShapeType()==TopAbs_SHELL || sh.ShapeType()==TopAbs_FACE ) {
1457 // changes for 0020677: EDF 1219 GEOM: MinDistance gives 0 instead of 20.88
1458 //S->Bounds(U1,U2,V1,V2); changed by
1459 ShapeAnalysis::GetFaceUVBounds(TopoDS::Face(tmpSh1),U1,U2,V1,V2);
1460 // end of changes for 020677 (dmv)
1461 Handle(Geom_RectangularTrimmedSurface) TrS1 =
1462 new Geom_RectangularTrimmedSurface(S,U1,(U1+U2)/2.,V1,V2);
1463 Handle(Geom_RectangularTrimmedSurface) TrS2 =
1464 new Geom_RectangularTrimmedSurface(S,(U1+U2)/2.,U2,V1,V2);
1467 TopoDS_Compound Comp;
1468 B.MakeCompound(Comp);
1469 B.MakeFace(F1,TrS1,1.e-7);
1471 B.MakeFace(F2,TrS2,1.e-7);
1473 Handle(ShapeFix_Shape) sfs = new ShapeFix_Shape;
1475 sfs->SetPrecision(1.e-6);
1476 sfs->SetMaxTolerance(1.0);
1478 tmpSh1 = sfs->Shape();
1482 if( S->IsKind(STANDARD_TYPE(Geom_SphericalSurface)) ) {
1483 Handle(Geom_SphericalSurface) SS = Handle(Geom_SphericalSurface)::DownCast(S);
1484 gp_Pnt PC = SS->Location();
1487 B.MakeVertex(V,PC,1.e-7);
1489 AddDist1 = SS->Radius();
1493 Handle(Geom_ToroidalSurface) TS = Handle(Geom_ToroidalSurface)::DownCast(S);
1494 gp_Ax3 ax3 = TS->Position();
1495 Handle(Geom_Circle) C = new Geom_Circle(ax3.Ax2(),TS->MajorRadius());
1498 B.MakeEdge(E,C,1.e-7);
1500 AddDist1 = TS->MinorRadius();
1510 bool IsChange2 = false;
1511 double AddDist2 = 0.0;
1513 for ( anExp.Init( aSh2, TopAbs_FACE ); anExp.More(); anExp.Next() ) {
1515 tmpSh2 = anExp.Current();
1518 TopoDS_Shape sh = aSh2;
1519 while(sh.ShapeType()==TopAbs_COMPOUND) {
1520 TopoDS_Iterator it(sh);
1523 Handle(Geom_Surface) S = BRep_Tool::Surface(TopoDS::Face(tmpSh2));
1524 if( S->IsKind(STANDARD_TYPE(Geom_SphericalSurface)) ||
1525 S->IsKind(STANDARD_TYPE(Geom_ToroidalSurface)) ) {
1526 if( sh.ShapeType()==TopAbs_SHELL || sh.ShapeType()==TopAbs_FACE ) {
1529 //S->Bounds(U1,U2,V1,V2);
1530 ShapeAnalysis::GetFaceUVBounds(TopoDS::Face(tmpSh2),U1,U2,V1,V2);
1531 Handle(Geom_RectangularTrimmedSurface) TrS1 =
1532 new Geom_RectangularTrimmedSurface(S,U1,(U1+U2)/2.,V1,V2);
1533 Handle(Geom_RectangularTrimmedSurface) TrS2 =
1534 new Geom_RectangularTrimmedSurface(S,(U1+U2)/2.,U2,V1,V2);
1537 TopoDS_Compound Comp;
1538 B.MakeCompound(Comp);
1539 B.MakeFace(F1,TrS1,1.e-7);
1541 B.MakeFace(F2,TrS2,1.e-7);
1543 Handle(ShapeFix_Shape) sfs = new ShapeFix_Shape;
1545 sfs->SetPrecision(1.e-6);
1546 sfs->SetMaxTolerance(1.0);
1548 tmpSh2 = sfs->Shape();
1552 if( S->IsKind(STANDARD_TYPE(Geom_SphericalSurface)) ) {
1553 Handle(Geom_SphericalSurface) SS = Handle(Geom_SphericalSurface)::DownCast(S);
1554 gp_Pnt PC = SS->Location();
1557 B.MakeVertex(V,PC,1.e-7);
1559 AddDist2 = SS->Radius();
1562 else if( S->IsKind(STANDARD_TYPE(Geom_ToroidalSurface)) ) {
1563 Handle(Geom_ToroidalSurface) TS = Handle(Geom_ToroidalSurface)::DownCast(S);
1564 gp_Ax3 ax3 = TS->Position();
1565 Handle(Geom_Circle) C = new Geom_Circle(ax3.Ax2(),TS->MajorRadius());
1568 B.MakeEdge(E,C,1.e-7);
1570 AddDist2 = TS->MinorRadius();
1581 if( !IsChange1 && !IsChange2 )
1584 BRepExtrema_DistShapeShape dst(tmpSh1,tmpSh2);
1586 double MinDist = 1.e9;
1587 gp_Pnt PMin1, PMin2, P1, P2;
1588 for (int i = 1; i <= dst.NbSolution(); i++) {
1589 P1 = dst.PointOnShape1(i);
1590 P2 = dst.PointOnShape2(i);
1591 Standard_Real Dist = P1.Distance(P2);
1592 if (MinDist > Dist) {
1603 gp_Dir aDir(gp_Vec(PMin1,PMin2));
1604 if( MinDist > (AddDist1+AddDist2) ) {
1605 Ptmp1 = gp_Pnt( PMin1.X() + aDir.X()*AddDist1,
1606 PMin1.Y() + aDir.Y()*AddDist1,
1607 PMin1.Z() + aDir.Z()*AddDist1 );
1608 Ptmp2 = gp_Pnt( PMin2.X() - aDir.X()*AddDist2,
1609 PMin2.Y() - aDir.Y()*AddDist2,
1610 PMin2.Z() - aDir.Z()*AddDist2 );
1611 return (MinDist - AddDist1 - AddDist2);
1614 if( AddDist1 > 0 ) {
1615 Ptmp1 = gp_Pnt( PMin1.X() + aDir.X()*AddDist1,
1616 PMin1.Y() + aDir.Y()*AddDist1,
1617 PMin1.Z() + aDir.Z()*AddDist1 );
1621 Ptmp2 = gp_Pnt( PMin2.X() - aDir.X()*AddDist2,
1622 PMin2.Y() - aDir.Y()*AddDist2,
1623 PMin2.Z() - aDir.Z()*AddDist2 );
1628 double res = MinDist - AddDist1 - AddDist2;
1629 if(res<0.) res = 0.0;
1635 static bool CheckSingularCase(const TopoDS_Shape& aSh1,
1636 const TopoDS_Shape& aSh2,
1639 TopExp_Explorer anExp;
1640 TopoDS_Shape tmpSh1, tmpSh2;
1642 for ( anExp.Init( aSh1, TopAbs_FACE ); anExp.More(); anExp.Next() ) {
1644 tmpSh1 = anExp.Current();
1647 Handle(Geom_Surface) S1 = BRep_Tool::Surface(TopoDS::Face(tmpSh1));
1648 if( S1->IsKind(STANDARD_TYPE(Geom_SphericalSurface)) ||
1649 S1->IsKind(STANDARD_TYPE(Geom_ToroidalSurface)) ) {
1651 for ( anExp.Init( aSh2, TopAbs_FACE ); anExp.More(); anExp.Next() ) {
1653 tmpSh2 = anExp.Current();
1654 Handle(Geom_Surface) S2 = BRep_Tool::Surface(TopoDS::Face(tmpSh2));
1655 GeomAPI_IntSS ISS(S1,S2,1.e-7);
1657 for(int i=1; i<=ISS.NbLines(); i++) {
1658 Handle(Geom_Curve) C3d = ISS.Line(i);
1661 B.MakeEdge(E,C3d,1.e-7);
1662 BRepExtrema_DistShapeShape dst(tmpSh2,E);
1664 gp_Pnt PMin1, PMin2, P1, P2;
1665 double MinDist = 1.e9;
1666 for (int i = 1; i <= dst.NbSolution(); i++) {
1667 P1 = dst.PointOnShape1(i);
1668 P2 = dst.PointOnShape2(i);
1669 Standard_Real Dist = P1.Distance(P2);
1670 if (MinDist > Dist) {
1684 for ( anExp.Init( aSh2, TopAbs_FACE ); anExp.More(); anExp.Next() ) {
1686 tmpSh1 = anExp.Current();
1689 Handle(Geom_Surface) S1 = BRep_Tool::Surface(TopoDS::Face(tmpSh1));
1690 if( S1->IsKind(STANDARD_TYPE(Geom_SphericalSurface)) ||
1691 S1->IsKind(STANDARD_TYPE(Geom_ToroidalSurface)) ) {
1693 for ( anExp.Init( aSh1, TopAbs_FACE ); anExp.More(); anExp.Next() ) {
1695 tmpSh2 = anExp.Current();
1696 Handle(Geom_Surface) S2 = BRep_Tool::Surface(TopoDS::Face(tmpSh2));
1697 GeomAPI_IntSS ISS(S1,S2,1.e-7);
1699 for(int i=1; i<=ISS.NbLines(); i++) {
1700 Handle(Geom_Curve) C3d = ISS.Line(i);
1703 B.MakeEdge(E,C3d,1.e-7);
1704 BRepExtrema_DistShapeShape dst(tmpSh2,E);
1707 double MinDist = 1.e9;
1708 for (int i = 1; i <= dst.NbSolution(); i++) {
1709 P1 = dst.PointOnShape1(i);
1710 P2 = dst.PointOnShape2(i);
1711 Standard_Real Dist = P1.Distance(P2);
1712 if (MinDist > Dist) {
1730 //=============================================================================
1734 //=============================================================================
1735 std::vector<bool> GEOMImpl_IMeasureOperations::AreCoordsInside(Handle(GEOM_Object) theShape,
1736 const std::vector<double>& coords,
1739 std::vector<bool> res;
1740 if (!theShape.IsNull()) {
1741 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1742 if (!aRefShape.IsNull()) {
1743 TopoDS_Shape aShape = aRefShape->GetValue();
1744 if (!aShape.IsNull()) {
1745 BRepClass3d_SolidClassifier SC(aShape);
1746 unsigned int nb_points = coords.size()/3;
1747 for (int i = 0; i < nb_points; i++) {
1748 double x = coords[3*i];
1749 double y = coords[3*i+1];
1750 double z = coords[3*i+2];
1751 gp_Pnt aPnt(x, y, z);
1752 SC.Perform(aPnt, tolerance);
1753 res.push_back( ( SC.State() == TopAbs_IN ) || ( SC.State() == TopAbs_ON ) );
1761 //=============================================================================
1765 //=============================================================================
1766 Standard_Real GEOMImpl_IMeasureOperations::GetMinDistance
1767 (Handle(GEOM_Object) theShape1, Handle(GEOM_Object) theShape2,
1768 Standard_Real& X1, Standard_Real& Y1, Standard_Real& Z1,
1769 Standard_Real& X2, Standard_Real& Y2, Standard_Real& Z2)
1772 Standard_Real MinDist = 1.e9;
1774 if (theShape1.IsNull() || theShape2.IsNull()) return MinDist;
1776 Handle(GEOM_Function) aRefShape1 = theShape1->GetLastFunction();
1777 Handle(GEOM_Function) aRefShape2 = theShape2->GetLastFunction();
1778 if (aRefShape1.IsNull() || aRefShape2.IsNull()) return MinDist;
1780 TopoDS_Shape aShape1 = aRefShape1->GetValue();
1781 TopoDS_Shape aShape2 = aRefShape2->GetValue();
1782 if (aShape1.IsNull() || aShape2.IsNull()) {
1783 SetErrorCode("One of Objects has NULL Shape");
1787 //Compute the parameters
1789 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
1793 // Issue 0020231: A min distance bug with torus and vertex.
1794 // Make GetMinDistance() return zero if a sole VERTEX is inside any of SOLIDs
1796 // which of shapes consists of only one vertex?
1797 TopExp_Explorer exp1(aShape1,TopAbs_VERTEX), exp2(aShape2,TopAbs_VERTEX);
1798 TopoDS_Shape V1 = exp1.More() ? exp1.Current() : TopoDS_Shape();
1799 TopoDS_Shape V2 = exp2.More() ? exp2.Current() : TopoDS_Shape();
1800 exp1.Next(); exp2.Next();
1801 if ( exp1.More() ) V1.Nullify();
1802 if ( exp2.More() ) V2.Nullify();
1803 // vertex and container of solids
1804 TopoDS_Shape V = V1.IsNull() ? V2 : V1;
1805 TopoDS_Shape S = V1.IsNull() ? aShape1 : aShape2;
1806 if ( !V.IsNull() ) {
1807 // classify vertex against solids
1808 gp_Pnt p = BRep_Tool::Pnt( TopoDS::Vertex( V ) );
1809 for ( exp1.Init( S, TopAbs_SOLID ); exp1.More(); exp1.Next() ) {
1810 BRepClass3d_SolidClassifier classifier( exp1.Current(), p, 1e-6);
1811 if ( classifier.State() == TopAbs_IN ) {
1812 p.Coord(X1, Y1, Z1);
1813 p.Coord(X2, Y2, Z2);
1819 // End Issue 0020231
1822 // additional workaround for bugs 19899, 19908 and 19910 from Mantis
1823 gp_Pnt Ptmp1, Ptmp2;
1824 double dist = CheckSingularCase(aShape1, aShape2, Ptmp1, Ptmp2);
1826 Ptmp1.Coord(X1, Y1, Z1);
1827 Ptmp2.Coord(X2, Y2, Z2);
1832 BRepExtrema_DistShapeShape dst (aShape1, aShape2);
1834 gp_Pnt PMin1, PMin2, P1, P2;
1836 for (int i = 1; i <= dst.NbSolution(); i++) {
1837 P1 = dst.PointOnShape1(i);
1838 P2 = dst.PointOnShape2(i);
1840 Standard_Real Dist = P1.Distance(P2);
1841 if (MinDist > Dist) {
1848 PMin1.Coord(X1, Y1, Z1);
1849 PMin2.Coord(X2, Y2, Z2);
1852 catch (Standard_Failure) {
1853 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1854 SetErrorCode(aFail->GetMessageString());
1862 //=======================================================================
1864 * Get coordinates of point
1866 //=======================================================================
1867 void GEOMImpl_IMeasureOperations::PointCoordinates (Handle(GEOM_Object) theShape,
1868 Standard_Real& theX, Standard_Real& theY, Standard_Real& theZ)
1872 if (theShape.IsNull())
1875 Handle(GEOM_Function) aRefShape = theShape->GetLastFunction();
1876 if (aRefShape.IsNull())
1879 TopoDS_Shape aShape = aRefShape->GetValue();
1880 if (aShape.IsNull() || aShape.ShapeType() != TopAbs_VERTEX)
1882 SetErrorCode( "Shape must be a vertex" );
1887 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
1890 gp_Pnt aPnt = BRep_Tool::Pnt( TopoDS::Vertex( aShape ) );
1897 catch (Standard_Failure)
1899 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1900 SetErrorCode( aFail->GetMessageString() );
1904 //=======================================================================
1906 * Compute angle (in degrees) between two lines
1908 //=======================================================================
1909 Standard_Real GEOMImpl_IMeasureOperations::GetAngle (Handle(GEOM_Object) theLine1,
1910 Handle(GEOM_Object) theLine2)
1914 Standard_Real anAngle = -1.0;
1916 if (theLine1.IsNull() || theLine2.IsNull())
1919 Handle(GEOM_Function) aRefLine1 = theLine1->GetLastFunction();
1920 Handle(GEOM_Function) aRefLine2 = theLine2->GetLastFunction();
1921 if (aRefLine1.IsNull() || aRefLine2.IsNull())
1924 TopoDS_Shape aLine1 = aRefLine1->GetValue();
1925 TopoDS_Shape aLine2 = aRefLine2->GetValue();
1926 if (aLine1.IsNull() || aLine2.IsNull() ||
1927 aLine1.ShapeType() != TopAbs_EDGE ||
1928 aLine2.ShapeType() != TopAbs_EDGE)
1930 SetErrorCode("Two edges must be given");
1935 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
1938 TopoDS_Edge E1 = TopoDS::Edge(aLine1);
1939 TopoDS_Edge E2 = TopoDS::Edge(aLine2);
1942 Handle(Geom_Curve) C1 = BRep_Tool::Curve(E1,fp,lp);
1943 Handle(Geom_Curve) C2 = BRep_Tool::Curve(E2,fp,lp);
1945 if ( C1.IsNull() || C2.IsNull() ||
1946 !C1->IsKind(STANDARD_TYPE(Geom_Line)) ||
1947 !C2->IsKind(STANDARD_TYPE(Geom_Line)))
1949 SetErrorCode("The edges must be linear");
1953 Handle(Geom_Line) L1 = Handle(Geom_Line)::DownCast(C1);
1954 Handle(Geom_Line) L2 = Handle(Geom_Line)::DownCast(C2);
1956 gp_Lin aLin1 = L1->Lin();
1957 gp_Lin aLin2 = L2->Lin();
1959 anAngle = aLin1.Angle(aLin2);
1960 anAngle /= PI180; // convert radians into degrees
1962 if (anAngle > 90.0) {
1963 anAngle = 180.0 - anAngle;
1968 catch (Standard_Failure)
1970 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
1971 SetErrorCode(aFail->GetMessageString());
1977 //=======================================================================
1979 * Compute angle (in degrees) between two vectors
1981 //=======================================================================
1982 Standard_Real GEOMImpl_IMeasureOperations::GetAngleBtwVectors (Handle(GEOM_Object) theVec1,
1983 Handle(GEOM_Object) theVec2)
1987 Standard_Real anAngle = -1.0;
1989 if (theVec1.IsNull() || theVec2.IsNull())
1992 Handle(GEOM_Function) aRefVec1 = theVec1->GetLastFunction();
1993 Handle(GEOM_Function) aRefVec2 = theVec2->GetLastFunction();
1994 if (aRefVec1.IsNull() || aRefVec2.IsNull())
1997 TopoDS_Shape aVec1 = aRefVec1->GetValue();
1998 TopoDS_Shape aVec2 = aRefVec2->GetValue();
1999 if (aVec1.IsNull() || aVec2.IsNull() ||
2000 aVec1.ShapeType() != TopAbs_EDGE ||
2001 aVec2.ShapeType() != TopAbs_EDGE)
2003 SetErrorCode("Two edges must be given");
2008 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
2011 TopoDS_Edge aE1 = TopoDS::Edge(aVec1);
2012 TopoDS_Edge aE2 = TopoDS::Edge(aVec2);
2014 TopoDS_Vertex aP11, aP12, aP21, aP22;
2015 TopExp::Vertices(aE1, aP11, aP12, Standard_True);
2016 TopExp::Vertices(aE2, aP21, aP22, Standard_True);
2017 if (aP11.IsNull() || aP12.IsNull() || aP21.IsNull() || aP22.IsNull()) {
2018 SetErrorCode("Bad edge given");
2022 gp_Vec aV1 (BRep_Tool::Pnt(aP11), BRep_Tool::Pnt(aP12));
2023 gp_Vec aV2 (BRep_Tool::Pnt(aP21), BRep_Tool::Pnt(aP22)) ;
2025 anAngle = aV1.Angle(aV2);
2026 anAngle /= PI180; // convert radians into degrees
2030 catch (Standard_Failure)
2032 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
2033 SetErrorCode(aFail->GetMessageString());
2040 //=============================================================================
2042 * CurveCurvatureByParam
2044 //=============================================================================
2045 Standard_Real GEOMImpl_IMeasureOperations::CurveCurvatureByParam
2046 (Handle(GEOM_Object) theCurve, Standard_Real& theParam)
2049 Standard_Real aRes = -1.0;
2051 if(theCurve.IsNull()) return aRes;
2053 Handle(GEOM_Function) aRefShape = theCurve->GetLastFunction();
2054 if(aRefShape.IsNull()) return aRes;
2056 TopoDS_Shape aShape = aRefShape->GetValue();
2057 if(aShape.IsNull()) {
2058 SetErrorCode("One of Objects has NULL Shape");
2062 Standard_Real aFP, aLP, aP;
2063 Handle(Geom_Curve) aCurve = BRep_Tool::Curve(TopoDS::Edge(aShape), aFP, aLP);
2064 aP = aFP + (aLP - aFP) * theParam;
2066 if(aCurve.IsNull()) return aRes;
2070 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
2073 GeomLProp_CLProps Prop = GeomLProp_CLProps
2074 (aCurve, aP, 2, Precision::Confusion());
2075 aRes = fabs(Prop.Curvature());
2078 catch (Standard_Failure) {
2079 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
2080 SetErrorCode(aFail->GetMessageString());
2084 if( aRes > Precision::Confusion() )
2093 //=============================================================================
2095 * CurveCurvatureByPoint
2097 //=============================================================================
2098 Standard_Real GEOMImpl_IMeasureOperations::CurveCurvatureByPoint
2099 (Handle(GEOM_Object) theCurve, Handle(GEOM_Object) thePoint)
2102 Standard_Real aRes = -1.0;
2104 if( theCurve.IsNull() || thePoint.IsNull() ) return aRes;
2106 Handle(GEOM_Function) aRefCurve = theCurve->GetLastFunction();
2107 Handle(GEOM_Function) aRefPoint = thePoint->GetLastFunction();
2108 if( aRefCurve.IsNull() || aRefPoint.IsNull() ) return aRes;
2110 TopoDS_Edge anEdge = TopoDS::Edge(aRefCurve->GetValue());
2111 TopoDS_Vertex aPnt = TopoDS::Vertex(aRefPoint->GetValue());
2112 if( anEdge.IsNull() || aPnt.IsNull() ) {
2113 SetErrorCode("One of Objects has NULL Shape");
2117 Standard_Real aFP, aLP;
2118 Handle(Geom_Curve) aCurve = BRep_Tool::Curve(anEdge, aFP, aLP);
2119 if(aCurve.IsNull()) return aRes;
2120 gp_Pnt aPoint = BRep_Tool::Pnt(aPnt);
2124 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
2127 GeomAPI_ProjectPointOnCurve PPCurve(aPoint, aCurve, aFP, aLP);
2128 if(PPCurve.NbPoints()>0) {
2129 GeomLProp_CLProps Prop = GeomLProp_CLProps
2130 (aCurve, PPCurve.LowerDistanceParameter(), 2, Precision::Confusion());
2131 aRes = fabs(Prop.Curvature());
2135 catch (Standard_Failure) {
2136 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
2137 SetErrorCode(aFail->GetMessageString());
2141 if( aRes > Precision::Confusion() )
2150 //=============================================================================
2152 * getSurfaceCurvatures
2154 //=============================================================================
2155 Standard_Real GEOMImpl_IMeasureOperations::getSurfaceCurvatures
2156 (const Handle(Geom_Surface)& aSurf,
2157 Standard_Real theUParam,
2158 Standard_Real theVParam,
2159 Standard_Boolean theNeedMaxCurv)
2162 Standard_Real aRes = 1.0;
2164 if (aSurf.IsNull()) return aRes;
2167 #if (OCC_VERSION_MAJOR << 16 | OCC_VERSION_MINOR << 8 | OCC_VERSION_MAINTENANCE) > 0x060100
2170 GeomLProp_SLProps Prop = GeomLProp_SLProps
2171 (aSurf, theUParam, theVParam, 2, Precision::Confusion());
2172 if(Prop.IsCurvatureDefined()) {
2173 if(Prop.IsUmbilic()) {
2174 //cout<<"is umbilic"<<endl;
2175 aRes = fabs(Prop.MeanCurvature());
2178 //cout<<"is not umbilic"<<endl;
2179 double c1 = fabs(Prop.MaxCurvature());
2180 double c2 = fabs(Prop.MinCurvature());
2189 catch (Standard_Failure) {
2190 Handle(Standard_Failure) aFail = Standard_Failure::Caught();
2191 SetErrorCode(aFail->GetMessageString());
2195 if( fabs(aRes) > Precision::Confusion() )
2204 //=============================================================================
2206 * MaxSurfaceCurvatureByParam
2208 //=============================================================================
2209 Standard_Real GEOMImpl_IMeasureOperations::MaxSurfaceCurvatureByParam
2210 (Handle(GEOM_Object) theSurf,
2211 Standard_Real& theUParam,
2212 Standard_Real& theVParam)
2215 Standard_Real aRes = -1.0;
2217 if (theSurf.IsNull()) return aRes;
2219 Handle(GEOM_Function) aRefShape = theSurf->GetLastFunction();
2220 if(aRefShape.IsNull()) return aRes;
2222 TopoDS_Shape aShape = aRefShape->GetValue();
2223 if(aShape.IsNull()) {
2224 SetErrorCode("One of Objects has NULL Shape");
2228 TopoDS_Face F = TopoDS::Face(aShape);
2229 Handle(Geom_Surface) aSurf = BRep_Tool::Surface(F);
2231 //Compute the parameters
2232 Standard_Real U1,U2,V1,V2;
2233 ShapeAnalysis::GetFaceUVBounds(F,U1,U2,V1,V2);
2234 Standard_Real U = U1 + (U2-U1)*theUParam;
2235 Standard_Real V = V1 + (V2-V1)*theVParam;
2237 return getSurfaceCurvatures(aSurf, U, V, true);
2241 //=============================================================================
2243 * MaxSurfaceCurvatureByPoint
2245 //=============================================================================
2246 Standard_Real GEOMImpl_IMeasureOperations::MaxSurfaceCurvatureByPoint
2247 (Handle(GEOM_Object) theSurf, Handle(GEOM_Object) thePoint)
2250 Standard_Real aRes = -1.0;
2252 if( theSurf.IsNull() || thePoint.IsNull() ) return aRes;
2254 Handle(GEOM_Function) aRefShape = theSurf->GetLastFunction();
2255 Handle(GEOM_Function) aRefPoint = thePoint->GetLastFunction();
2256 if( aRefShape.IsNull() || aRefPoint.IsNull() ) return aRes;
2258 TopoDS_Face aFace = TopoDS::Face(aRefShape->GetValue());
2259 TopoDS_Vertex aPnt = TopoDS::Vertex(aRefPoint->GetValue());
2260 if( aFace.IsNull() || aPnt.IsNull() ) {
2261 SetErrorCode("One of Objects has NULL Shape");
2265 Handle(Geom_Surface) aSurf = BRep_Tool::Surface(aFace);
2266 if(aSurf.IsNull()) return aRes;
2267 gp_Pnt aPoint = BRep_Tool::Pnt(aPnt);
2269 //Compute the parameters
2270 ShapeAnalysis_Surface sas(aSurf);
2271 gp_Pnt2d UV = sas.ValueOfUV(aPoint,Precision::Confusion());
2273 return getSurfaceCurvatures(aSurf, UV.X(), UV.Y(), true);
2277 //=============================================================================
2279 * MinSurfaceCurvatureByParam
2281 //=============================================================================
2282 Standard_Real GEOMImpl_IMeasureOperations::MinSurfaceCurvatureByParam
2283 (Handle(GEOM_Object) theSurf,
2284 Standard_Real& theUParam,
2285 Standard_Real& theVParam)
2288 Standard_Real aRes = -1.0;
2290 if (theSurf.IsNull()) return aRes;
2292 Handle(GEOM_Function) aRefShape = theSurf->GetLastFunction();
2293 if(aRefShape.IsNull()) return aRes;
2295 TopoDS_Shape aShape = aRefShape->GetValue();
2296 if(aShape.IsNull()) {
2297 SetErrorCode("One of Objects has NULL Shape");
2301 TopoDS_Face F = TopoDS::Face(aShape);
2302 Handle(Geom_Surface) aSurf = BRep_Tool::Surface(F);
2304 //Compute the parameters
2305 Standard_Real U1,U2,V1,V2;
2306 ShapeAnalysis::GetFaceUVBounds(F,U1,U2,V1,V2);
2307 Standard_Real U = U1 + (U2-U1)*theUParam;
2308 Standard_Real V = V1 + (V2-V1)*theVParam;
2310 return getSurfaceCurvatures(aSurf, U, V, false);
2314 //=============================================================================
2316 * MinSurfaceCurvatureByPoint
2318 //=============================================================================
2319 Standard_Real GEOMImpl_IMeasureOperations::MinSurfaceCurvatureByPoint
2320 (Handle(GEOM_Object) theSurf, Handle(GEOM_Object) thePoint)
2323 Standard_Real aRes = -1.0;
2325 if( theSurf.IsNull() || thePoint.IsNull() ) return aRes;
2327 Handle(GEOM_Function) aRefShape = theSurf->GetLastFunction();
2328 Handle(GEOM_Function) aRefPoint = thePoint->GetLastFunction();
2329 if( aRefShape.IsNull() || aRefPoint.IsNull() ) return aRes;
2331 TopoDS_Face aFace = TopoDS::Face(aRefShape->GetValue());
2332 TopoDS_Vertex aPnt = TopoDS::Vertex(aRefPoint->GetValue());
2333 if( aFace.IsNull() || aPnt.IsNull() ) {
2334 SetErrorCode("One of Objects has NULL Shape");
2338 Handle(Geom_Surface) aSurf = BRep_Tool::Surface(aFace);
2339 if(aSurf.IsNull()) return aRes;
2340 gp_Pnt aPoint = BRep_Tool::Pnt(aPnt);
2342 //Compute the parameters
2343 ShapeAnalysis_Surface sas(aSurf);
2344 gp_Pnt2d UV = sas.ValueOfUV(aPoint,Precision::Confusion());
2346 return getSurfaceCurvatures(aSurf, UV.X(), UV.Y(), false);
2350 //=======================================================================
2351 //function : StructuralDump
2352 //purpose : Structural (data exchange) style of output.
2353 //=======================================================================
2354 void GEOMImpl_IMeasureOperations::StructuralDump (const BRepCheck_Analyzer& theAna,
2355 const TopoDS_Shape& theShape,
2356 TCollection_AsciiString& theDump)
2360 theDump += " -- The Shape has problems :\n";
2361 theDump += " Check Count\n";
2362 theDump += " ------------------------------------------------\n";
2364 Standard_Integer last_stat = (Standard_Integer)BRepCheck_CheckFail;
2365 Handle(TColStd_HArray1OfInteger) NbProblems =
2366 new TColStd_HArray1OfInteger(1, last_stat);
2367 for (i = 1; i <= last_stat; i++)
2368 NbProblems->SetValue(i,0);
2370 Handle(TopTools_HSequenceOfShape) sl;
2371 sl = new TopTools_HSequenceOfShape();
2372 TopTools_DataMapOfShapeListOfShape theMap;
2374 GetProblemShapes(theAna, theShape, sl, NbProblems, theMap);
2377 Standard_Integer count = 0;
2378 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidPointOnCurve);
2380 theDump += " Invalid Point on Curve ................... ";
2381 theDump += TCollection_AsciiString(count) + "\n";
2383 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidPointOnCurveOnSurface);
2385 theDump += " Invalid Point on CurveOnSurface .......... ";
2386 theDump += TCollection_AsciiString(count) + "\n";
2388 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidPointOnSurface);
2390 theDump += " Invalid Point on Surface ................. ";
2391 theDump += TCollection_AsciiString(count) + "\n";
2393 count = NbProblems->Value((Standard_Integer)BRepCheck_No3DCurve);
2395 theDump += " No 3D Curve .............................. ";
2396 theDump += TCollection_AsciiString(count) + "\n";
2398 count = NbProblems->Value((Standard_Integer)BRepCheck_Multiple3DCurve);
2400 theDump += " Multiple 3D Curve ........................ ";
2401 theDump += TCollection_AsciiString(count) + "\n";
2403 count = NbProblems->Value((Standard_Integer)BRepCheck_Invalid3DCurve);
2405 theDump += " Invalid 3D Curve ......................... ";
2406 theDump += TCollection_AsciiString(count) + "\n";
2408 count = NbProblems->Value((Standard_Integer)BRepCheck_NoCurveOnSurface);
2410 theDump += " No Curve on Surface ...................... ";
2411 theDump += TCollection_AsciiString(count) + "\n";
2413 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidCurveOnSurface);
2415 theDump += " Invalid Curve on Surface ................. ";
2416 theDump += TCollection_AsciiString(count) + "\n";
2418 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidCurveOnClosedSurface);
2420 theDump += " Invalid Curve on closed Surface .......... ";
2421 theDump += TCollection_AsciiString(count) + "\n";
2423 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidSameRangeFlag);
2425 theDump += " Invalid SameRange Flag ................... ";
2426 theDump += TCollection_AsciiString(count) + "\n";
2428 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidSameParameterFlag);
2430 theDump += " Invalid SameParameter Flag ............... ";
2431 theDump += TCollection_AsciiString(count) + "\n";
2433 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidDegeneratedFlag);
2435 theDump += " Invalid Degenerated Flag ................. ";
2436 theDump += TCollection_AsciiString(count) + "\n";
2438 count = NbProblems->Value((Standard_Integer)BRepCheck_FreeEdge);
2440 theDump += " Free Edge ................................ ";
2441 theDump += TCollection_AsciiString(count) + "\n";
2443 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidMultiConnexity);
2445 theDump += " Invalid MultiConnexity ................... ";
2446 theDump += TCollection_AsciiString(count) + "\n";
2448 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidRange);
2450 theDump += " Invalid Range ............................ ";
2451 theDump += TCollection_AsciiString(count) + "\n";
2453 count = NbProblems->Value((Standard_Integer)BRepCheck_EmptyWire);
2455 theDump += " Empty Wire ............................... ";
2456 theDump += TCollection_AsciiString(count) + "\n";
2458 count = NbProblems->Value((Standard_Integer)BRepCheck_RedundantEdge);
2460 theDump += " Redundant Edge ........................... ";
2461 theDump += TCollection_AsciiString(count) + "\n";
2463 count = NbProblems->Value((Standard_Integer)BRepCheck_SelfIntersectingWire);
2465 theDump += " Self Intersecting Wire ................... ";
2466 theDump += TCollection_AsciiString(count) + "\n";
2468 count = NbProblems->Value((Standard_Integer)BRepCheck_NoSurface);
2470 theDump += " No Surface ............................... ";
2471 theDump += TCollection_AsciiString(count) + "\n";
2473 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidWire);
2475 theDump += " Invalid Wire ............................. ";
2476 theDump += TCollection_AsciiString(count) + "\n";
2478 count = NbProblems->Value((Standard_Integer)BRepCheck_RedundantWire);
2480 theDump += " Redundant Wire ........................... ";
2481 theDump += TCollection_AsciiString(count) + "\n";
2483 count = NbProblems->Value((Standard_Integer)BRepCheck_IntersectingWires);
2485 theDump += " Intersecting Wires ....................... ";
2486 theDump += TCollection_AsciiString(count) + "\n";
2488 count = NbProblems->Value((Standard_Integer)BRepCheck_InvalidImbricationOfWires);
2490 theDump += " Invalid Imbrication of Wires ............. ";
2491 theDump += TCollection_AsciiString(count) + "\n";
2493 count = NbProblems->Value((Standard_Integer)BRepCheck_EmptyShell);
2495 theDump += " Empty Shell .............................. ";
2496 theDump += TCollection_AsciiString(count) + "\n";
2498 count = NbProblems->Value((Standard_Integer)BRepCheck_RedundantFace);
2500 theDump += " Redundant Face ........................... ";
2501 theDump += TCollection_AsciiString(count) + "\n";
2503 count = NbProblems->Value((Standard_Integer)BRepCheck_UnorientableShape);
2505 theDump += " Unorientable Shape ....................... ";
2506 theDump += TCollection_AsciiString(count) + "\n";
2508 count = NbProblems->Value((Standard_Integer)BRepCheck_NotClosed);
2510 theDump += " Not Closed ............................... ";
2511 theDump += TCollection_AsciiString(count) + "\n";
2513 count = NbProblems->Value((Standard_Integer)BRepCheck_NotConnected);
2515 theDump += " Not Connected ............................ ";
2516 theDump += TCollection_AsciiString(count) + "\n";
2518 count = NbProblems->Value((Standard_Integer)BRepCheck_SubshapeNotInShape);
2520 theDump += " Subshape not in Shape .................... ";
2521 theDump += TCollection_AsciiString(count) + "\n";
2523 count = NbProblems->Value((Standard_Integer)BRepCheck_BadOrientation);
2525 theDump += " Bad Orientation .......................... ";
2526 theDump += TCollection_AsciiString(count) + "\n";
2528 count = NbProblems->Value((Standard_Integer)BRepCheck_BadOrientationOfSubshape);
2530 theDump += " Bad Orientation of Subshape .............. ";
2531 theDump += TCollection_AsciiString(count) + "\n";
2533 count = NbProblems->Value((Standard_Integer)BRepCheck_CheckFail);
2535 theDump += " checkshape failure ....................... ";
2536 theDump += TCollection_AsciiString(count) + "\n";
2539 theDump += " ------------------------------------------------\n";
2540 theDump += "*** Shapes with problems : ";
2541 theDump += TCollection_AsciiString(sl->Length()) + "\n";
2543 Standard_Integer nbv, nbe, nbw, nbf, nbs, nbo;
2544 nbv = nbe = nbw = nbf = nbs = nbo = 0;
2546 for (i = 1; i <= sl->Length(); i++) {
2547 TopoDS_Shape shi = sl->Value(i);
2548 TopAbs_ShapeEnum sti = shi.ShapeType();
2550 case TopAbs_VERTEX : nbv++; break;
2551 case TopAbs_EDGE : nbe++; break;
2552 case TopAbs_WIRE : nbw++; break;
2553 case TopAbs_FACE : nbf++; break;
2554 case TopAbs_SHELL : nbs++; break;
2555 case TopAbs_SOLID : nbo++; break;
2561 theDump += "VERTEX : ";
2562 if (nbv < 10) theDump += " ";
2563 theDump += TCollection_AsciiString(nbv) + "\n";
2566 theDump += "EDGE : ";
2567 if (nbe < 10) theDump += " ";
2568 theDump += TCollection_AsciiString(nbe) + "\n";
2571 theDump += "WIRE : ";
2572 if (nbw < 10) theDump += " ";
2573 theDump += TCollection_AsciiString(nbw) + "\n";
2576 theDump += "FACE : ";
2577 if (nbf < 10) theDump += " ";
2578 theDump += TCollection_AsciiString(nbf) + "\n";
2581 theDump += "SHELL : ";
2582 if (nbs < 10) theDump += " ";
2583 theDump += TCollection_AsciiString(nbs) + "\n";
2586 theDump += "SOLID : ";
2587 if (nbo < 10) theDump += " ";
2588 theDump += TCollection_AsciiString(nbo) + "\n";
2593 //=======================================================================
2594 //function : GetProblemShapes
2595 // purpose : for StructuralDump
2596 //=======================================================================
2597 void GEOMImpl_IMeasureOperations::GetProblemShapes (const BRepCheck_Analyzer& theAna,
2598 const TopoDS_Shape& theShape,
2599 Handle(TopTools_HSequenceOfShape)& sl,
2600 Handle(TColStd_HArray1OfInteger)& NbProblems,
2601 TopTools_DataMapOfShapeListOfShape& theMap)
2603 for (TopoDS_Iterator iter(theShape); iter.More(); iter.Next()) {
2604 GetProblemShapes(theAna, iter.Value(), sl, NbProblems, theMap);
2606 TopAbs_ShapeEnum styp = theShape.ShapeType();
2607 BRepCheck_ListIteratorOfListOfStatus itl;
2608 TopTools_ListOfShape empty;
2609 if (!theMap.IsBound(theShape)) {
2610 theMap.Bind(theShape,empty);
2612 if (!theAna.Result(theShape).IsNull()) {
2613 itl.Initialize(theAna.Result(theShape)->Status());
2614 // !!! May be, we have to print all the problems, not only the first one ?
2615 if (itl.Value() != BRepCheck_NoError) {
2616 sl->Append(theShape);
2617 BRepCheck_Status stat = itl.Value();
2618 NbProblems->SetValue((Standard_Integer)stat,
2619 NbProblems->Value((Standard_Integer)stat) + 1);
2626 GetProblemSub(theAna, theShape, sl, NbProblems, TopAbs_VERTEX, theMap);
2629 GetProblemSub(theAna, theShape, sl, NbProblems, TopAbs_WIRE, theMap);
2630 GetProblemSub(theAna, theShape, sl, NbProblems, TopAbs_EDGE, theMap);
2631 GetProblemSub(theAna, theShape, sl, NbProblems, TopAbs_VERTEX, theMap);
2636 GetProblemSub(theAna, theShape, sl, NbProblems, TopAbs_SHELL, theMap);
2643 //=======================================================================
2644 //function : Contains
2645 //=======================================================================
2646 static Standard_Boolean Contains (const TopTools_ListOfShape& L,
2647 const TopoDS_Shape& S)
2649 TopTools_ListIteratorOfListOfShape it;
2650 for (it.Initialize(L); it.More(); it.Next()) {
2651 if (it.Value().IsSame(S)) {
2652 return Standard_True;
2655 return Standard_False;
2658 //=======================================================================
2659 //function : GetProblemSub
2660 // purpose : for StructuralDump
2661 //=======================================================================
2662 void GEOMImpl_IMeasureOperations::GetProblemSub (const BRepCheck_Analyzer& theAna,
2663 const TopoDS_Shape& theShape,
2664 Handle(TopTools_HSequenceOfShape)& sl,
2665 Handle(TColStd_HArray1OfInteger)& NbProblems,
2666 const TopAbs_ShapeEnum Subtype,
2667 TopTools_DataMapOfShapeListOfShape& theMap)
2669 BRepCheck_ListIteratorOfListOfStatus itl;
2670 TopExp_Explorer exp;
2671 for (exp.Init(theShape, Subtype); exp.More(); exp.Next()) {
2672 const TopoDS_Shape& sub = exp.Current();
2674 const Handle(BRepCheck_Result)& res = theAna.Result(sub);
2675 for (res->InitContextIterator();
2676 res->MoreShapeInContext();
2677 res->NextShapeInContext()) {
2678 if (res->ContextualShape().IsSame(theShape) && !Contains(theMap(sub), theShape)) {
2679 theMap(sub).Append(theShape);
2680 itl.Initialize(res->StatusOnShape());
2682 if (itl.Value() != BRepCheck_NoError) {
2683 Standard_Integer ii = 0;
2685 for (ii = 1; ii <= sl->Length(); ii++)
2686 if (sl->Value(ii).IsSame(sub)) break;
2688 if (ii > sl->Length()) {
2690 BRepCheck_Status stat = itl.Value();
2691 NbProblems->SetValue((Standard_Integer)stat,
2692 NbProblems->Value((Standard_Integer)stat) + 1);
2694 for (ii = 1; ii <= sl->Length(); ii++)
2695 if (sl->Value(ii).IsSame(theShape)) break;
2696 if (ii > sl->Length()) {
2697 sl->Append(theShape);
2698 BRepCheck_Status stat = itl.Value();
2699 NbProblems->SetValue((Standard_Integer)stat,
2700 NbProblems->Value((Standard_Integer)stat) + 1);