projects / occt.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
0031642: Visualization - crash in Graphic3d_Structure::SetVisual() on redisplaying...
[occt.git] / src / SWDRAW / SWDRAW_ShapeAnalysis.cxx
CommitLineData
b311480e 1// Created on: 1999-03-09
2// Created by: data exchange team
3// Copyright (c) 1999-1999 Matra Datavision
973c2be1 4// Copyright (c) 1999-2014 OPEN CASCADE SAS
b311480e 5//
973c2be1 6// This file is part of Open CASCADE Technology software library.
b311480e 7//
d5f74e42 8// This library is free software; you can redistribute it and/or modify it under
9// the terms of the GNU Lesser General Public License version 2.1 as published
973c2be1 10// by the Free Software Foundation, with special exception defined in the file
11// OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
12// distribution for complete text of the license and disclaimer of any warranty.
b311480e 13//
973c2be1 14// Alternatively, this file may be used under the terms of Open CASCADE
15// commercial license or contractual agreement.
b311480e 16
7fd59977 17// sln 19.11.2001. Bug 2: Correction of output of 'statshape' draw function.
18
42cf5bc1 19#include <BRep_Builder.hxx>
20#include <BRep_Tool.hxx>
21#include <BRepBuilderAPI.hxx>
22#include <BRepBuilderAPI_MakeFace.hxx>
23#include <BRepGProp.hxx>
24#include <BRepTopAdaptor_FClass2d.hxx>
7fd59977 25#include <DBRep.hxx>
42cf5bc1 26#include <Draw.hxx>
7fd59977 27#include <DrawTrSurf.hxx>
42cf5bc1 28#include <Geom2d_Curve.hxx>
7fd59977 29#include <Geom_Curve.hxx>
30#include <Geom_Surface.hxx>
31#include <GeomAPI_ProjectPointOnSurf.hxx>
7fd59977 32#include <gp_Pnt.hxx>
33#include <gp_Pnt2d.hxx>
7fd59977 34#include <GProp_GProps.hxx>
42cf5bc1 35#include <Precision.hxx>
36#include <ShapeAnalysis.hxx>
37#include <ShapeAnalysis_Curve.hxx>
38#include <ShapeAnalysis_Edge.hxx>
39#include <ShapeAnalysis_FreeBoundData.hxx>
7fd59977 40#include <ShapeAnalysis_FreeBounds.hxx>
41#include <ShapeAnalysis_FreeBoundsProperties.hxx>
42cf5bc1 42#include <ShapeAnalysis_ShapeContents.hxx>
43#include <ShapeAnalysis_ShapeTolerance.hxx>
44#include <ShapeAnalysis_Surface.hxx>
ef57920a 45#include <ShapeAnalysis_Wire.hxx>
42cf5bc1 46#include <ShapeExtend_WireData.hxx>
47#include <ShapeFix_FreeBounds.hxx>
48#include <SWDRAW.hxx>
49#include <SWDRAW_ShapeAnalysis.hxx>
50#include <TCollection_AsciiString.hxx>
51#include <TColStd_MapOfInteger.hxx>
52#include <TopExp.hxx>
53#include <TopExp_Explorer.hxx>
54#include <TopoDS.hxx>
55#include <TopoDS_Compound.hxx>
56#include <TopoDS_Edge.hxx>
57#include <TopoDS_Face.hxx>
58#include <TopoDS_Iterator.hxx>
59#include <TopoDS_Shape.hxx>
60#include <TopoDS_Vertex.hxx>
61#include <TopoDS_Wire.hxx>
62#include <TopTools_HSequenceOfShape.hxx>
f6b08ecf 63#include <Adaptor3d_CurveOnSurface.hxx>
64#include <BRepAdaptor_HCurve2d.hxx>
65#include <BRepAdaptor_HSurface.hxx>
7fd59977 66
42cf5bc1 67#include <stdio.h>
7fd59977 68static Standard_Integer tolerance
69 (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
70{
71 if (argc < 2) { di<<"myshape : analyses a shape\n"
72 << "myshape [mode] valmin : sub-shapes over valmin\n"
73 <<"myshape [mode] valmin valmax : between valmin and valmax\n"
74 <<"myshape [mode] 0 valmax : below valmax\n"
586db386 75 <<" more : all shapes(D) v vertices e edges f faces c combined(faces)\n";
7fd59977 76 return (argc < 2 ? 0 : 1 /* Error */); }
77 Standard_CString arg1 = argv[1];
78 TopoDS_Shape Shape = DBRep::Get(arg1);
79 if (Shape.IsNull()) { di<<"Shape unknown : "<<arg1<<"\n"; return 1 /* Error */; }
80
81 ShapeAnalysis_ShapeTolerance sat;
82 if (argc == 2) {
83 sat.InitTolerance();
84 sat.AddTolerance (Shape);
85 di<<"Tolerance MAX="<<sat.GlobalTolerance(1)<<" AVG="<<sat.GlobalTolerance(0)<<" MIN="<<sat.GlobalTolerance(-1)<<"\n";
86 sat.InitTolerance();
87 sat.AddTolerance (Shape,TopAbs_FACE);
88 di<<"FACE : MAX="<<sat.GlobalTolerance(1)<<" AVG="<<sat.GlobalTolerance(0)<<" MIN="<<sat.GlobalTolerance(-1)<<"\n";
89 sat.InitTolerance();
90 sat.AddTolerance (Shape,TopAbs_EDGE);
91 di<<"EDGE : MAX="<<sat.GlobalTolerance(1)<<" AVG="<<sat.GlobalTolerance(0)<<" MIN="<<sat.GlobalTolerance(-1)<<"\n";
92 sat.InitTolerance();
93 sat.AddTolerance (Shape,TopAbs_VERTEX);
94 di<<"VERTEX : MAX="<<sat.GlobalTolerance(1)<<" AVG="<<sat.GlobalTolerance(0)<<" MIN="<<sat.GlobalTolerance(-1)<<"\n";
95 } else {
96// Analyse plus fouillee
97 Standard_Integer nextarg = 2;
98 char opt = argv[2][0];
99 TopAbs_ShapeEnum type = TopAbs_SHAPE;
100 if (opt == 'a') nextarg = 3;
101 else if (opt == 'c') { nextarg = 3; type = TopAbs_SHELL; }
102 else if (opt == 'e') { nextarg = 3; type = TopAbs_EDGE; }
103 else if (opt == 'f') { nextarg = 3; type = TopAbs_FACE; }
104 else if (opt == 'v') { nextarg = 3; type = TopAbs_VERTEX; }
105 if (nextarg < 3) opt = 'a';
106
107 Standard_Real tol1= 0., tol2= 0.;
91322f44 108 if (nextarg < argc) tol1 = Draw::Atof (argv[nextarg]);
109 if (nextarg < argc-1) tol2 = Draw::Atof (argv[nextarg+1]);
110// Standard_Real tol = Draw::Atof (argv[2]);
7fd59977 111 Handle(TopTools_HSequenceOfShape) list = sat.InTolerance (Shape,tol1,tol2,type);
112 Standard_Integer i, nb = list->Length();
113 switch (type) {
114 case TopAbs_VERTEX : di<<"Analysing Vertices gives "; break;
115 case TopAbs_EDGE : di<<"Analysing Edges gives "; break;
116 case TopAbs_FACE : di<<"Analysing Faces gives "; break;
117 case TopAbs_SHELL : di<<"Analysing Shells,Faces+content gives "; break;
118 default : di<<"Analysing all sub-shapes gives ";
119 }
120 if (tol1 == 0) di<<nb<<" Shapes below tol="<<tol2;
121 else if (tol2 == 0) di<<nb<<" Shapes over tol="<<tol1;
122 else di<<nb<<" Shapes between tol1="<<tol1<<" and tol2="<<tol2;
123 if (nb == 1) di<<" , named tol_1";
124 if (nb > 1) di<<" , named tol_1 to tol_"<<nb;
125 di<<"\n";
126 char nomsh[30];
127 for (i = 1; i <= nb; i ++) {
91322f44 128 Sprintf (nomsh,"tol_%d",i);
7fd59977 129 DBRep::Set (nomsh,list->Value(i));
130 }
131 }
132 return 0;
133}
134
135static Standard_Integer projface
136 (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
137{
586db386 138 if (argc < 4) { di<<"Give FACE name and X Y [Z]\n"; return 1 /* Error */; }
7fd59977 139 Standard_CString arg1 = argv[1];
140 TopoDS_Shape Shape = DBRep::Get(arg1);
141 if (Shape.IsNull()) { di<<"Shape unknown : "<<arg1<<"\n"; return 1 /* Error */; }
586db386 142 if (Shape.ShapeType() != TopAbs_FACE) { di<<"Not a face\n"; return 1 /* Error */; }
7fd59977 143 TopoDS_Face F = TopoDS::Face (Shape);
144 Handle(Geom_Surface) thesurf = BRep_Tool::Surface (F); // pas locface
145// On y va
146 Standard_Real X,Y,Z,U,V;
91322f44 147 X = U = Draw::Atof (argv[2]);
148 Y = V = Draw::Atof (argv[3]);
7fd59977 149 if (argc > 4) {
91322f44 150 Z = Draw::Atof (argv[4]);
7fd59977 151 gp_Pnt P3D (X,Y,Z);
152 di<<" Point 3D X = "<<X<<" Y = "<<Y<<" Z = "<<Z<<"\n";
153 Standard_Real uf, ul, vf, vl;
154 thesurf->Bounds(uf, ul, vf, vl);
155
156 if (Precision::IsInfinite(uf)) uf = -1000;
157 if (Precision::IsInfinite(ul)) ul = 1000;
158 if (Precision::IsInfinite(vf)) vf = -1000;
159 if (Precision::IsInfinite(vl)) vl = 1000;
160 Standard_Real du = Abs(ul-uf)/10; Standard_Real dv = Abs(vl-vf)/10;
161
162 GeomAPI_ProjectPointOnSurf proj(P3D, thesurf, uf-du, ul+du, vf-dv, vl+dv);
163 Standard_Integer sol, nPSurf = proj.NbPoints();
586db386 164 di<<" Found "<<nPSurf<<" Points\n";
7fd59977 165 for (sol = 1; sol <= nPSurf; sol ++) {
166 di<<"n0 "<<sol<<" Distance "<<proj.Distance(sol);
167 proj.Parameters(sol, U,V);
168 di<<" U = "<<U<<" V = "<<V<<"\n";
169// reprojection
170 P3D = thesurf->Value (U,V);
171 di<<" => reproj X = "<<P3D.X()<<" Y = "<<P3D.Y()<<" Z = "<<P3D.Z()<<"\n";
172 }
173// Que donne ShapeTool ?
174 P3D.SetCoord (X,Y,Z);
175 Handle(ShapeAnalysis_Surface) su = new ShapeAnalysis_Surface(thesurf);
176 gp_Pnt2d suval = su->ValueOfUV (P3D,BRep_Tool::Tolerance(F));
177 suval.Coord(U,V);
178 di<<"** ShapeAnalysis_Surface gives U = "<<U<<" V = "<<V<<"\n";
179 P3D = thesurf->Value(U,V);
180 di<<" => reproj X = "<<P3D.X()<<" Y = "<<P3D.Y()<<" Z = "<<P3D.Z()<<"\n";
181
182 } else {
183 di<<" Point UV U = "<<U<<" V = "<<V<<"\n";
184 gp_Pnt P3D = thesurf->Value(U,V);
185 di<<" => proj X = "<<P3D.X()<<" Y = "<<P3D.Y()<<" Z = "<<P3D.Z()<<"\n";
186 }
187 return 0;
188}
189
190static Standard_Integer projcurve
191 (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
192{
193// admet une EDGE ou une CURVE
586db386 194 if (argc < 5) { di<<"Give : EDGE X Y Z or CURVE3D X Y Z or CURVE3D first last X Y Z\n"; return 1 /* Error */; }
7fd59977 195 const char * arg1 = argv[1];
196 Standard_Integer i0 = 0;
197 TopoDS_Shape Shape = DBRep::Get(arg1);
198
199 Handle(Geom_Curve) C; Standard_Real cf,cl;
200 if (!Shape.IsNull()) {
201 if (Shape.ShapeType() == TopAbs_EDGE) {
202 TopoDS_Edge E = TopoDS::Edge(Shape);
203 C = BRep_Tool::Curve (E,cf,cl);
204 di<<"Edge "<<arg1<<" Params from "<<cf<<" to "<<cl<<"\n";
205 } else {
586db386 206 di<<"SHAPE "<<arg1<<" Not an EDGE\n"; return 1 /* Error */;
7fd59977 207 }
208 } else {
209 const char* a1 = arg1;
210 C = DrawTrSurf::GetCurve (a1);
211 if (C.IsNull())
586db386 212 { di<<arg1<<" neither EDGE nor CURVE 3D\n"; return 1 /* Error */; }
7fd59977 213 cf = C->FirstParameter(); cl = C->LastParameter();
91322f44 214 if (argc >= 7) { cf = Draw::Atof (argv[2]); cl = Draw::Atof (argv[3]); i0 = 2; }
7fd59977 215 di<<"Curve 3D "<<arg1<<" Params from "<<cf<<" to "<<cl<<"\n";
216 }
217
218 Standard_Real X,Y,Z;
91322f44 219 X = Draw::Atof (argv[2+i0]);
220 Y = Draw::Atof (argv[3+i0]);
221 Z = Draw::Atof (argv[4+i0]);
7fd59977 222 di<<"Precision (BRepBuilderAPI) : "<<BRepBuilderAPI::Precision()
223 <<" Projection : "<<X<<" "<<Y<<" "<<Z<<"\n";
224
225 gp_Pnt P3D (X,Y,Z);
226 gp_Pnt res; Standard_Real param,dist;
227
228 dist = ShapeAnalysis_Curve().Project (C,P3D,BRepBuilderAPI::Precision(),res,param, cf,cl);
229 res.Coord(X,Y,Z);
586db386 230 di<<"Result : "<<X<<" "<<Y<<" "<<Z<<"\nParam = "<<param<<" Gap = "<<dist<<"\n";
7fd59977 231 return 0;
232}
f6b08ecf 233static Standard_Integer projpcurve
234(Draw_Interpretor& di, Standard_Integer argc, const char** argv)
235{
236 // admet une EDGE ou une CURVE
237 if (argc < 7)
238 {
239 di << "Give : projpcurve edge face tol X Y Z [start_param]\n";
240 return 1;
241 }
242
243 TopoDS_Edge aEdge = TopoDS::Edge(DBRep::Get(argv[1]));
244 if (aEdge.IsNull())
245 {
246 di << "SHAPE " << argv[1] << " is not an EDGE\n"; return 1 /* Error */;
247 }
248 TopoDS_Face aFace = TopoDS::Face(DBRep::Get(argv[2]));
249 if (aFace.IsNull())
250 {
251 di << "SHAPE " << argv[2] << " is not a FACE\n"; return 1 /* Error */;
252 }
253
254 Standard_Real aTol = Draw::Atof(argv[3]);
255
256
257 Standard_Real X = Draw::Atof(argv[4]);
258 Standard_Real Y = Draw::Atof(argv[5]);
259 Standard_Real Z = Draw::Atof(argv[6]);
260
261 gp_Pnt aP3D(X, Y, Z);
262
263 Standard_Boolean IsStartPoint = Standard_False;
264 Standard_Real startpar = 0.;
265 if (argc > 7)
266 {
267 startpar = Draw::Atof(argv[7]);
268 IsStartPoint = Standard_True;
269 }
270
271 Adaptor3d_CurveOnSurface aCOnS =
272 Adaptor3d_CurveOnSurface(new BRepAdaptor_HCurve2d(BRepAdaptor_Curve2d(aEdge, aFace)),
273 new BRepAdaptor_HSurface(BRepAdaptor_Surface(aFace, Standard_False)));
274
275 gp_Pnt aPnt;
276 Standard_Real aParam;
277 ShapeAnalysis_Curve aTool;
278 Standard_Real aDist = RealLast();
279 if (IsStartPoint)
280 {
281 aDist = aTool.NextProject(startpar, aCOnS, aP3D, aTol, aPnt, aParam);
282 }
283 else
284 {
285 aDist = aTool.Project(aCOnS, aP3D, aTol, aPnt, aParam, Standard_False);
286 }
287
288 di << "Point:" << "\n" << aPnt.X() << " " << aPnt.Y() << " " << aPnt.Z() << "\n";
289 di << "Param: " << aParam << "\n";
290 di << "Dist: " << aDist << "\n";
291 return 0;
292}
7fd59977 293
294static Standard_Integer anaface
295 (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
296{
586db386 297 if (argc < 2) { di<<"donner un nom de face\n"; return 1 /* Error */; }
7fd59977 298 Standard_CString arg1 = argv[1];
299 TopoDS_Shape Shape = DBRep::Get(arg1);
586db386 300 if (Shape.IsNull()) { di<<arg1<<" inconnu\n"; return 1 /* Error */; }
301// if (Shape.ShapeType() != TopAbs_FACE) { di<<"Pas une FACE\n"; return 1 /* Error */; }
7fd59977 302
303// On regarde les wires, sont-ils bien fermes
304 Standard_Integer nbw = 0, nbe = 0;
305 Handle(Geom_Surface) surface;
306 TopoDS_Face Face;
307 if (Shape.ShapeType() == TopAbs_FACE) {
308 Face = TopoDS::Face(Shape);
309 surface = BRep_Tool::Surface (Face); // pas locface
310// TopLoc_Location locface;
311 di<<"Face, surface type = "<<surface->DynamicType()->Name()<<" Tol="<<BRep_Tool::Tolerance(Face);
312 if (Face.Orientation() == TopAbs_REVERSED) {
313 di<<" (REV)";
314// Face.Orientation(TopAbs_FORWARD);
315 }
316 di<<"\n";
317 }
586db386 318 else { di<<"Analyse Wires, Ignore Faces\n"; }
7fd59977 319 Standard_Boolean iasurf = !surface.IsNull();
320//:sw if (!Face.IsNull()) STW.SetFace (Face);
321//:sw else STW.SetPrecision (BRepBuilderAPI::Precision());
322
323 for (TopExp_Explorer exp(Shape,TopAbs_WIRE); exp.More(); exp.Next()) {
324 TopoDS_Wire Wire = TopoDS::Wire (exp.Current()); nbw ++; nbe = 0;
325 di<<"WIRE "<<nbw;
586db386 326 if (Wire.Orientation() == TopAbs_FORWARD) di<<" (FWD)\n";
327 else di<<" (REV)\n";
7fd59977 328 gp_Pnt fin,debut;
329 gp_Pnt2d finuv,debuv; gp_XY baseuv;
330 TopoDS_Vertex fv,lv;
331 Standard_Real maxp3d = 0., maxuv = 0., maxvtx = 0.;
332 Standard_Real df3d= 0.,dl3d= 0.,dvtx= 0.,duv= 0.;
333 Standard_Boolean iaw2d = iasurf;
334 Standard_Real u1= 0.,u2= 0.,v1= 0.,v2= 0.,umin= 0.,umax= 0.,vmin= 0.,vmax= 0., totcross= 0.;
335
336//:sw STW.ClearWire(); STW.Load (Wire);
586db386 337// if (STW.Reorder()) di<<"Wire reordered to explore\n";
7fd59977 338// for (TopExp_Explorer exe(Wire.Oriented(TopAbs_FORWARD),TopAbs_EDGE); exe.More(); exe.Next()) {
339// TopoDS_Edge Edge = TopoDS::Edge (exe.Current()); nbe ++;
340//:sw for (nbe = 1; nbe <= STW.NbEdges(); nbe ++) {
341//:sw TopoDS_Edge Edge = STW.Edge(nbe);
342 TopoDS_Iterator it (Wire); //:sw: ShapeTool_Wire -> TopoDS_Iterator (it is the same)
343 for (; it.More(); it.Next()) {
344 TopoDS_Edge Edge = TopoDS::Edge (it.Value());
345
346 di<<"Wire "<<nbw<<", Edge "<<nbe;
347 if (Edge.Orientation() == TopAbs_FORWARD) di<<" (FWD";
348 else di<<" (REV";
586db386 349 di<<" , Tol= "<<BRep_Tool::Tolerance (Edge)<<" )\n";
7fd59977 350 Standard_Real f3d,l3d,f2d,l2d;
351 Handle(Geom_Curve) curve3d = BRep_Tool::Curve (Edge,f3d,l3d);
352 Handle(Geom2d_Curve) curve2d;
353 if (iasurf) curve2d = BRep_Tool::CurveOnSurface(Edge,Face,f2d,l2d);
354 Standard_Boolean ia2d = !curve2d.IsNull(); iaw2d |= ia2d;
586db386 355 if (!ia2d) di<<"-- No PCurve\n";
356 if (curve3d.IsNull()) di<<"-- no Curve 3d\n";
7fd59977 357
358// On va tacher de calculer les positions et les comparer
359 gp_Pnt2d fuv,luv;
360 if (Edge.Orientation() == TopAbs_FORWARD)
c48e2889 361 {
362 TopExp::Vertices (Edge, fv, lv);
363 if (ia2d) BRep_Tool::UVPoints (Edge, Face, fuv, luv);
364 }
7fd59977 365 else
c48e2889 366 {
367 TopExp::Vertices (Edge, lv, fv);
368 if (ia2d) BRep_Tool::UVPoints (Edge, Face, luv, fuv);
369 }
7fd59977 370 gp_Pnt fp = BRep_Tool::Pnt (fv);
371 gp_Pnt lp = BRep_Tool::Pnt (lv);
372 gp_Pnt fxyz, lxyz;
c48e2889 373 if (ia2d)
374 {
7fd59977 375 surface->D0 (fuv.X(),fuv.Y(),fxyz);
376 surface->D0 (luv.X(),luv.Y(),lxyz);
377 df3d = fp.Distance (fxyz); maxp3d = Max (maxp3d,df3d);
378 dl3d = lp.Distance (lxyz); maxp3d = Max (maxp3d,dl3d);
c48e2889 379 if (nbe > 1)
380 {
381 duv = finuv.Distance (fuv);
382 maxuv = Max (maxuv, duv);
383 }
7fd59977 384// et les min-max
385 u1 = Min (fuv.X(),luv.X()); u2 = Max (fuv.X(),luv.X());
386 v1 = Min (fuv.Y(),luv.Y()); v2 = Max (fuv.Y(),luv.Y());
387 if (nbe == 1) { umin = u1; umax = u2; vmin = v1; vmax = v2; }
388 else { umin = Min(umin,u1); umax = Max(umax,u2); vmin = Min(vmin,v1); vmax = Max(vmax,v2); }
389// et la classification directe
390 if (nbe == 1) { baseuv = fuv.XY(); totcross = 0; }
391 else {
392 gp_XY buv1 = fuv.XY() - baseuv;
393 gp_XY buv2 = luv.XY() - baseuv;
394 totcross += buv2.Crossed(buv1);
395 }
396 }
397
398// Resultats ...
399 if (nbe == 1) { debut = fxyz; debuv = fuv; }
400 else {
401 duv = finuv.Distance (fuv); maxuv = Max (maxuv, duv);
402 dvtx = fin.Distance (fxyz); maxvtx = Max (maxvtx,dvtx);
403 di<<" Fin("<<nbe-1<<")-Debut("<<nbe<<"): DISTANCE="<<dvtx;
404 if (ia2d) di<<" DeltaUV="<<duv;
405 di<<" Tol(Fin)="<<BRep_Tool::Tolerance(lv)<<"\n";
406 }
407 fin = lxyz; finuv = luv;
408
409 di<<"-- Deb : VTX="<<fp.X()<<" , "<<fp.Y()<<" , "<<fp.Z()
410 <<"\n XYZ="<<fxyz.X()<<" , "<<fxyz.Y()<<" , "<<fxyz.Z()
411 <<"\n UV="<<fuv.X()<<" , "<<fuv.Y()<<" -- D.UV/3D="<<df3d<<"\n";
412 di<<"-- Fin : VTX="<<lp.X()<<" , "<<lp.Y()<<" , "<<lp.Z()
413 <<"\n XYZ="<<lxyz.X()<<" , "<<lxyz.Y()<<" , "<<lxyz.Z()
414 <<"\n UV="<<luv.X()<<" , "<<luv.Y()<<" -- D.UV/3D="<<dl3d<<"\n";
415 }
416 dvtx = fin.Distance (debut); maxvtx = Max (maxvtx,dvtx);
417 if (iaw2d) { duv = finuv.Distance (debuv); maxuv = Max (maxuv, duv); }
418 di<<" Fin("<<nbe<<")-Debut(1): DISTANCE="<<dvtx;
419 if (iaw2d) di <<" DeltaUV="<<duv;
420 di<<" Tol(Fin)="<<BRep_Tool::Tolerance(lv)<<"\n";
421
422 di<<" Wire "<<nbw<<" Max : Dist.Vertex="<<maxvtx;
423 if (iaw2d) di<<" Ecart UV/3D="<<maxp3d<<" DeltaUV="<<maxuv;
424 di<<"\n";
425// Min Max
426 if (iaw2d) {
427 di<<"TotCross="<<totcross;
428 di<<" UMin-Max:"<<umin<<" , "<<umax<<" VMin-Max:"<<vmin<<" , "<<vmax<<"\n";
429 Standard_Real difu = umax - umin, difv = vmax - vmin;
430 GProp_GProps G;
431 BRepGProp::SurfaceProperties(Face,G);
586db386 432 if (G.Mass() > 0) di<<"GProps:Mass Out\n";
433 else di<<"GProps:Mass In\n";
7fd59977 434/// return (G.Mass() > 0);
435 BRepTopAdaptor_FClass2d fcl (Face,BRep_Tool::Tolerance(Face));
586db386 436 if (fcl.PerformInfinitePoint () == TopAbs_OUT) di<<"Classifier Infinite : Out\n";
437 else di<<"Classifier Infinite : In\n";
7fd59977 438 gp_Pnt2d pcl;
439 pcl.SetCoord(umin-difu,vmin-difv);
586db386 440 if (fcl.Perform (pcl) == TopAbs_OUT) di<<"Classifier UMin-VMin : Out\n";
7fd59977 441 pcl.SetCoord(umin-difu,vmax+difv);
586db386 442 if (fcl.Perform (pcl) == TopAbs_OUT) di<<"Classifier UMin-VMax : Out\n";
7fd59977 443 pcl.SetCoord(umax+difu,vmin-difv);
586db386 444 if (fcl.Perform (pcl) == TopAbs_OUT) di<<"Classifier UMax-VMin : Out\n";
7fd59977 445 pcl.SetCoord(umax+difu,vmax+difv);
586db386 446 if (fcl.Perform (pcl) == TopAbs_OUT) di<<"Classifier UMax-VMax : Out\n";
7fd59977 447 }
448 }
586db386 449 if (ShapeAnalysis::IsOuterBound (Face)) di<<"ShapeAnalysis: Outer Bound\n";
450 else di<<"ShapeAnalysis: Not Outer Bound\n";
451 di<<" Total "<<nbw<<" Wire(s)\n";
7fd59977 452 return 0;
453}
454
455static Standard_Integer XSHAPE_statshape(Draw_Interpretor& di, Standard_Integer argc, const char** argv)
456{
586db386 457 if (argc < 2) { di<<"ShapeName; + options : prefix for particular cases : what to get\n"; return 1 /* Error */; }
7fd59977 458 Standard_CString arg1 = argv[1];
459 Standard_CString arg2 = NULL;
460 Standard_CString arg3 = NULL;
461 if (argc > 2) arg2 = argv[2];
462 if (argc > 3) arg3 = argv[3];
463 TopoDS_Shape Shape = DBRep::Get(arg1);
464 if (Shape.IsNull()) { di<<"Shape unknown : "<<arg1<<"\n"; return 1 /* Error */; }
465 char nompart[100];
466 Standard_Integer nb;
467
468 ShapeAnalysis_ShapeContents analyzer;
469 if(argc > 3) {
470 analyzer.ModifyBigSplineMode()=(strstr("bigspl",arg3)!=NULL);
471 analyzer.ModifyIndirectMode()=(strstr("indsur",arg3)!=NULL);
472 analyzer.ModifyOffestSurfaceMode()=(strstr("ofsur",arg3)!=NULL);
473 analyzer.ModifyTrimmed3dMode()=(strstr("trc3d",arg3)!=NULL);
474 analyzer.ModifyOffsetCurveMode()=(strstr("ofcur",arg3)!=NULL);
475 analyzer.ModifyTrimmed2dMode()=(strstr("trc2d",arg3)!=NULL);
476 }
477 analyzer.Perform(Shape);
478
586db386 479 di<<"Count Item\n----- ----\n";
7fd59977 480 nb = analyzer.NbEdges();
586db386 481 if(nb>0) di<<nb<<" EDGE (Oriented)\n";
7fd59977 482 nb = analyzer.NbSharedEdges();
586db386 483 if(nb>0) di<<nb<<" EDGE (Shared)\n";
7fd59977 484 nb = analyzer.NbFreeEdges();
586db386 485 if(nb>0) di<<nb<<" EDGE (Free)\n";
7fd59977 486 nb = analyzer.NbFaces();
586db386 487 if(nb>0) di<<nb<<" FACE\n";
7fd59977 488 nb = analyzer.NbFreeFaces();
586db386 489 if(nb>0) di<<nb<<" FACE (Free)\n";
7fd59977 490 nb = analyzer.NbFreeWires();
586db386 491 if(nb>0) di<<nb<<" WIRE (Free)\n";
7fd59977 492 nb = analyzer.NbShells();
586db386 493 if(nb>0) di<<nb<<" SHELL\n";
7fd59977 494 nb = analyzer.NbSolids();
586db386 495 if(nb>0) di<<nb<<" SOLID\n";
7fd59977 496 nb = analyzer.NbVertices();
586db386 497 if(nb>0) di<<nb<<" VERTEX (Oriented)\n";
7fd59977 498 nb = analyzer.NbSharedVertices();
586db386 499 if(nb>0) di<<nb<<" VERTEX (Shared)\n";
7fd59977 500 nb = analyzer.NbWires();
586db386 501 if(nb>0) di<<nb<<" WIRE\n";
7fd59977 502 nb = analyzer.NbFaceWithSevWires();
586db386 503 if(nb>0) di<<nb<<" Face with more than one wire\n";
7fd59977 504 nb = analyzer.NbNoPCurve();
586db386 505 if(nb>0) di<<nb<<" No pcurve\n";
7fd59977 506 nb = analyzer.NbSolidsWithVoids();
586db386 507 if(nb>0) di<<nb<<" SOLID with voids\n";
7fd59977 508 nb = analyzer.NbWireWitnSeam();
586db386 509 if(nb>0) di<<nb<<" Wire(s) with one seam edge\n";
7fd59977 510 nb = analyzer.NbWireWithSevSeams();
586db386 511 if(nb>0) di<<nb<<" Wire(s) with several seam edges\n";
7fd59977 512 nb = analyzer.NbBigSplines();
586db386 513 if(nb>0) di<<nb<<" bigspl : BSpline > 8192 poles\n";
7fd59977 514 nb = analyzer.NbBezierSurf();
586db386 515 if(nb>0) di<<nb<<" bezsur : BezierSurface\n";
7fd59977 516 nb = analyzer.NbBSplibeSurf();
586db386 517 if(nb>0) di<<nb<<" bspsur : BSplineSurface\n";
7fd59977 518 nb = analyzer.NbC0Curves();
586db386 519 if(nb>0) di<<nb<<" c0curv : Curve Only C0\n";
7fd59977 520 nb = analyzer.NbC0Surfaces();
586db386 521 if(nb>0) di<<nb<<" c0surf : Surface Only C0\n";
7fd59977 522 nb = analyzer.NbIndirectSurf();
586db386 523 if(nb>0) di<<nb<<" indsur : Indirect Surface\n";
7fd59977 524 nb = analyzer.NbOffsetCurves();
586db386 525 if(nb>0) di<<nb<<" ofcur : Offset Curve(s)\n";
7fd59977 526 nb = analyzer.NbOffsetSurf();
586db386 527 if(nb>0) di<<nb<<" ofsur : Offset Surface\n";
7fd59977 528 nb = analyzer.NbTrimmedCurve2d();
586db386 529 if(nb>0) di<<nb<<" trc2d : Trimmed Curve2d\n";
7fd59977 530 nb = analyzer.NbTrimmedCurve3d();
586db386 531 if(nb>0) di<<nb<<" trc3d : Trimmed Curve3d\n";
7fd59977 532 nb = analyzer.NbTrimSurf();
586db386 533 if(nb>0) di<<nb<<" trimsu : RectangularTrimmedSurface\n";
7fd59977 534
535 if(arg3 ==NULL) return 0;
536
537 Handle(TopTools_HSequenceOfShape) sec;
538 if(analyzer.ModifyBigSplineMode()) {
539 sec = analyzer.BigSplineSec();
540 for(Standard_Integer i = 1; i <= sec->Length(); i++) {
91322f44 541 Sprintf(nompart,"%s_bigspl_%d",arg2,i);
7fd59977 542 DBRep::Set (nompart,sec->Value(i));
543 }
544 }
545 if(analyzer.ModifyIndirectMode()) {
546 sec = analyzer.IndirectSec();
547 for(Standard_Integer i = 1; i <= sec->Length(); i++) {
91322f44 548 Sprintf(nompart,"%s_indsur_%d",arg2,i);
7fd59977 549 DBRep::Set (nompart,sec->Value(i));
550 }
551 }
552 if(analyzer.ModifyOffestSurfaceMode()) {
553 sec = analyzer.OffsetSurfaceSec();
554 for(Standard_Integer i = 1; i <= sec->Length(); i++) {
91322f44 555 Sprintf(nompart,"%s_ofsur_%d",arg2,i);
7fd59977 556 DBRep::Set (nompart,sec->Value(i));
557 }
558 }
559 if(analyzer.ModifyTrimmed3dMode()) {
560 sec = analyzer.Trimmed3dSec();
561 for(Standard_Integer i = 1; i <= sec->Length(); i++) {
91322f44 562 Sprintf(nompart,"%s_trc3d_%d",arg2,i);
7fd59977 563 DBRep::Set (nompart,sec->Value(i));
564 }
565 }
566 if(analyzer.ModifyOffsetCurveMode()) {
567 sec = analyzer.OffsetCurveSec();
568 for(Standard_Integer i = 1; i <= sec->Length(); i++) {
91322f44 569 Sprintf(nompart,"%s_ofcur_%d",arg2,i);
7fd59977 570 DBRep::Set (nompart,sec->Value(i));
571 }
572 }
573 if(analyzer.ModifyTrimmed2dMode()) {
574 sec = analyzer.Trimmed2dSec();
575 for(Standard_Integer i = 1; i <= sec->Length(); i++) {
91322f44 576 Sprintf(nompart,"%s_trc2d_%d",arg2,i);
7fd59977 577 DBRep::Set (nompart,sec->Value(i));
578 }
579 }
580 return 0;
581}
582
583static Standard_Integer XSHAPE_comptoledge
584 (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
585{
586 if ( argc <2 ) {
587 di << "Use:\n"
588 "> comptol shape [nbpoints=371] [prefix]\n\n"
589 "Computes real tolerance of edges in the shape \n"
590 "as maximal deviation of 3d curve and pcurves.\n"
591 "Deviation is computed by nbpoints sample points (default is 371).\n"
592 "Gives the max, min and average value on all edges in the shape\n"
593 "If prefix is defined, edges with maximal real tolerance and\n"
586db386 594 "relation (and corresponding faces) will be saved with such names\n\n";
7fd59977 595 return 0;
596 }
597
598 TopoDS_Shape shape = DBRep::Get(argv[1]);
586db386 599 if (shape.IsNull()) { di << "Shape name \"" << argv[1] << "\" is invalid\n"; return 1; }
7fd59977 600 Standard_Integer nbpnts = 371;
601 Standard_CString prefix = 0;
602 if ( argc >2 ) {
603 if ( IsDigit(argv[2][0]) ) {
91322f44 604 nbpnts = Draw::Atoi(argv[2]);
7fd59977 605 if ( nbpnts <2 ) nbpnts = 2;
606 if ( argc >3 ) prefix = argv[3];
607 }
608 else prefix = argv[2];
609 }
610
611 Standard_Integer num = 0;
612 Standard_Real max=0, min=0, ave=0.;
613 Standard_Real relmax=0, relmin=0, relave=0.;
614 ShapeAnalysis_Edge sae;
615 TopoDS_Edge edmax, edmaxrel;
616 for ( TopExp_Explorer exp(shape,TopAbs_EDGE); exp.More(); exp.Next() ) {
617 Standard_Real tol;
618 TopoDS_Edge edge = TopoDS::Edge(exp.Current());
619 sae.CheckSameParameter ( edge, tol, nbpnts );
620 Standard_Real t = BRep_Tool::Tolerance(edge);
621 Standard_Real rel = tol / ( t > Precision::Confusion() ? t : Precision::Confusion() );
622 ave += tol;
623 relave += rel;
624 if ( ! num ) { max = min = tol; relmax = relmin = rel; edmax = edmaxrel = edge; }
625 else {
626 if ( max < tol ) { max = tol; edmax = edge; }
627 if ( min > tol ) min = tol;
628 if ( relmax < rel ) { relmax = rel; edmaxrel = edge; }
629 if ( relmin > rel ) relmin = rel;
630 }
631 num++;
632 }
586db386 633 if ( ! num ) { di << "No edges found in the shape\n"; return 1; }
7fd59977 634
635 di << "Edges tolerance computed by " << nbpnts << " points: \n"
636 "MAX=" << max << " AVG=" << ave/num << " MIN=" << min << "\n";
637 di << "Relation real tolerance / tolerance set in edge\n"
638 "MAX=" << relmax << " AVG=" << relave/num << " MIN=" << relmin << "\n";
639 if ( prefix && prefix[0] ) {
640 char name[21];
91322f44 641 Sprintf ( name, "%.10s_edge_tol", prefix );
7fd59977 642 DBRep::Set (name,edmax);
643 di << "Edge with max tolerance saved to " << name;
644 if ( edmax.IsSame ( edmaxrel ) ) di << "\n";
645 else {
91322f44 646 Sprintf ( name, "%.10s_edge_rel", prefix );
7fd59977 647 DBRep::Set (name,edmaxrel);
648 di << "; edge with max relation saved to " << name << "\n";
649 }
650 Standard_Integer num1=0;
651 for ( TopExp_Explorer fac(shape,TopAbs_FACE); fac.More(); fac.Next() ) {
652 TopoDS_Face face = TopoDS::Face ( fac.Current() );
653 for ( TopExp_Explorer ed(face,TopAbs_EDGE); ed.More(); ed.Next() ) {
654 TopoDS_Edge edge = TopoDS::Edge ( ed.Current() );
655 if ( edge.IsSame ( edmax ) || edge.IsSame ( edmaxrel ) ) {
656 if ( ! num1 ) di << "Concerned faces saved to shapes ";
91322f44 657 Sprintf ( name, "%.10s_%d", prefix, num1+1 );
7fd59977 658 DBRep::Set (name,face);
04232180 659 //std::cout << ( num1 ? ", " : "" ) << name;
7fd59977 660 if (num1 == 0) {
661 di << "" << name;
662 } else {
663 di << ", " << name;
664 }
665 num1++;
666 break;
667 }
668 }
669 }
670 if ( num1 >0 ) di << "\n";
671 }
672 return 0;
673}
674
675
676//=======================================================================
677//function : freebounds
678//purpose :
679//=======================================================================
680static Standard_Integer freebounds (Draw_Interpretor& di,
681 Standard_Integer n, const char** a)
682{
683 if ((n < 3) || (n > 5)) return 1;
684 TopoDS_Shape shape = DBRep::Get(a[1]);
685 if (shape.IsNull()) return 1;
91322f44 686 Standard_Real toler = Draw::Atof (a[2]);
7fd59977 687 Standard_Boolean splitclosed = Standard_False, splitopen = Standard_False;
dde68833 688 if ( n > 3) splitclosed = Draw::Atoi (a[3]) != 0;
689 if ( n > 4) splitopen = Draw::Atoi (a[4]) != 0;
7fd59977 690
691 ShapeAnalysis_FreeBounds F;
692 if (toler <= 0)
693 F = ShapeAnalysis_FreeBounds (shape, splitclosed, splitopen);
694 else
695 F = ShapeAnalysis_FreeBounds (shape, toler, splitclosed, splitopen);
696
697 char name[100];
698 TopoDS_Shape wires = F.GetClosedWires();
91322f44 699 Sprintf (name, "%s_c", a[1]);
7fd59977 700 DBRep::Set (name, wires);
586db386 701 di << name << " - closed wires\n";
7fd59977 702 wires = F.GetOpenWires();
91322f44 703 Sprintf (name, "%s_o", a[1]);
7fd59977 704 DBRep::Set (name, wires);
586db386 705 di << name << " - open wires\n";
7fd59977 706
707 return 0;
708}
709
710
711//=======================================================================
712//function : PrintProps
713//purpose : auxilary for FreeBoundsProps
714//=======================================================================
715static void PrintProps(Standard_Integer i,
716 const Handle(ShapeAnalysis_FreeBoundData)& fbd,
717 Draw_Interpretor& di)
718{
719 char str[100];
720
721 Standard_Real area = fbd->Area();
722 Standard_Real perimeter = fbd->Perimeter();
723 Standard_Real ratio = fbd->Ratio();
724 Standard_Real width = fbd->Width();
725 Standard_Integer notch = fbd->NbNotches();
91322f44 726 Sprintf(str," %d\t%12.5f\t%12.5f\t%12.5f\t%12.5f\t%d", i, area, perimeter, ratio, width, notch);
7fd59977 727 di<<str<<"\n";
728}
729
730
731//=======================================================================
732//function : FreeBoundsProps
733//purpose :
734//=======================================================================
735static Standard_Integer FreeBoundsProps(Draw_Interpretor& di,
736 Standard_Integer n, const char** a)
737{
738 if ( (n < 2)||(n > 5) ) {
586db386 739 di<<"Usage : freeprops shapename [tolerance [splitclosed [splitopen]]]\n";
7fd59977 740 return 1;
741 }
742
743 TopoDS_Shape source = DBRep::Get(a[1]);
744 if (source.IsNull()) {
745 di<<"Error : unknown shape "<<a[1]<<"\n";
746 return 1;
747 }
748
749 Standard_Real toler = 0.;
750 Standard_Boolean splitclosed = Standard_False, splitopen = Standard_False;
91322f44 751 if (n > 2) toler = Draw::Atof(a[2]);
dde68833 752 if (n > 3) splitclosed = Draw::Atoi(a[3]) != 0;
753 if (n > 4) splitopen = Draw::Atoi(a[4]) != 0;
7fd59977 754 ShapeAnalysis_FreeBoundsProperties analyzer;
755 if (toler > 0)
756 analyzer.Init(source, toler, splitclosed, splitopen);
757 else
758 analyzer.Init(source, splitclosed, splitopen);
759 analyzer.Perform();
760 TopoDS_Compound closed, open;
761 BRep_Builder B;
762
763 Standard_Integer nb = analyzer.NbClosedFreeBounds();
764 di<<"\n";
586db386 765 di<<" \tArea mm2\tLength mm\tRatio L/W\tWidth mm\tNb noth\n";
7fd59977 766 B.MakeCompound(closed);
767 if (nb) {
586db386 768 di<<"Closed bounds properties\n";
7fd59977 769 for (Standard_Integer i=1; i <= nb; i++) {
770 Handle(ShapeAnalysis_FreeBoundData) fbd = analyzer.ClosedFreeBound(i);
771 PrintProps(i, fbd, di);
772 B.Add(closed,fbd->FreeBound());
773 }
774 }
775
776 nb = analyzer.NbOpenFreeBounds();
777 B.MakeCompound(open);
778 if (nb) {
586db386 779 di<<"Open bounds properties\n";
7fd59977 780 for (Standard_Integer i=1; i <= nb; i++) {
781 Handle(ShapeAnalysis_FreeBoundData) fbd = analyzer.OpenFreeBound(i);
782 PrintProps(i, fbd, di);
783 B.Add(open,fbd->FreeBound());
784 }
785 }
786
787 char name[100];
91322f44 788 Sprintf (name, "%s_c",a[1]);
7fd59977 789 di << name << " - closed wires, ";
790 DBRep::Set(name, closed);
91322f44 791 Sprintf (name, "%s_o",a[1]);
586db386 792 di << name << " - closed wires \n";
7fd59977 793 DBRep::Set(name, open);
794 return 0;
795}
796
797
798//=======================================================================
799//function : closefreebounds
800//purpose :
801//=======================================================================
802static Standard_Integer closefreebounds (Draw_Interpretor& di,
803 Standard_Integer n, const char** a)
804{
805 if ((n < 4) || (n > 6)) return 1;
806 TopoDS_Shape shape = DBRep::Get(a[1]);
807 if (shape.IsNull()) return 1;
91322f44 808 Standard_Real sewtoler = Draw::Atof (a[2]), closetoler = Draw::Atof (a[3]);
7fd59977 809 Standard_Boolean splitclosed = Standard_False, splitopen = Standard_False;
dde68833 810 if ( n > 4) splitclosed = Draw::Atoi (a[3]) != 0;
811 if ( n > 5) splitopen = Draw::Atoi (a[4]) != 0;
7fd59977 812
813 ShapeFix_FreeBounds F;
814 if (sewtoler <= 0)
815 F = ShapeFix_FreeBounds (shape, closetoler, splitclosed, splitopen);
816 else
817 F = ShapeFix_FreeBounds (shape, sewtoler, closetoler, splitclosed, splitopen);
818
819 char name[100];
820 TopoDS_Shape wires = F.GetClosedWires();
91322f44 821 Sprintf (name, "%s_c", a[1]);
7fd59977 822 DBRep::Set (name, wires);
586db386 823 di << name << " - closed wires\n";
7fd59977 824 wires = F.GetOpenWires();
91322f44 825 Sprintf (name, "%s_o", a[1]);
7fd59977 826 DBRep::Set (name, wires);
586db386 827 di << name << " - open wires\n";
7fd59977 828
829 return 0;
830}
831
832
7fd59977 833static Standard_Integer getareacontour (Draw_Interpretor& di,
834 Standard_Integer n, const char** a)
835{
836 if (n < 2) return 1;
837 TopoDS_Shape shape = DBRep::Get(a[1]);
838 if (shape.IsNull()) {
586db386 839 di<<"Shape is not defined\n";
7fd59977 840 return 1;
841 }
842 if(shape.ShapeType() != TopAbs_WIRE) {
586db386 843 di<<"invalid type of argument\n";
7fd59977 844 return 1;
845 }
846 //Handle(ShapeExtend_WireData) asewd = new ShapeExtend_WireData(TopoDS::Wire(shape));
847 Standard_Real anArea = ShapeAnalysis::ContourArea(TopoDS::Wire(shape));
848 di<<"Area = "<<anArea<<"\n";
849 return 0;
850}
851
ef57920a
V		 852
853static Standard_Integer checkselfintersection
854 (Draw_Interpretor& di, Standard_Integer argc, const char** argv)
855{
856 if (argc < 2)
857 {
586db386 858 di<<"Call please \"checkselfintersection wire [face]\"\n";
ef57920a
V		 859 return 1;
860 }
861
862 // Get wire.
863 const char* arg1 = argv[1];
864 TopoDS_Shape wire = DBRep::Get(arg1);
865 if (wire.IsNull() || wire.ShapeType() != TopAbs_WIRE)
866 {
586db386 867 di<<"A null shape or not a wire is used.\n";
ef57920a
V		 868 return 2;
869 }
870
871 // Get face if the user provided us with a face.
872 TopoDS_Shape face;
873 if (argc > 2)
874 {
875 const char* arg2 = argv[2];
876 face = DBRep::Get(arg2);
877 if (face.IsNull() || face.ShapeType() != TopAbs_FACE)
878 {
586db386 879 di<<"A null shape or not a face is used.\n";
ef57920a
V		 880 return 3;
881 }
882 }
883
884 // If the face is null, make a plane inside the wire.
885 if (face.IsNull())
886 {
887 BRepBuilderAPI_MakeFace mkFace(TopoDS::Wire(wire), true);
888 if (mkFace.IsDone())
889 face = mkFace.Face();
890 else
891 {
586db386 892 di<<"Can't make a face for the wire. Provide please a face for analysis.\n";
ef57920a
V		 893 return 4;
894 }
895 }
896
897 ShapeAnalysis_Wire analyser(TopoDS::Wire(wire), TopoDS::Face(face), Precision::Confusion());
898 Standard_Boolean result = analyser.CheckSelfIntersection();
899
900 if (result == Standard_True)
586db386 901 di<<"A self-intersecting wire.\n";
ef57920a 902 else
586db386 903 di<<"Not self-intersecting wire.\n";
ef57920a
V		 904 return 0;
905}
906
c31e0def 907static Standard_Integer checkedge(Draw_Interpretor& di, Standard_Integer argc, const char** argv)
908{
909 if (argc < 2)
910 {
586db386 911 di<<"Call please \"checkedge edge [face]\"\n";
c31e0def 912 return 1;
913 }
914
915 // Get edge.
916 const char* arg1 = argv[1];
917 TopoDS_Shape edge = DBRep::Get(arg1);
918 if (edge.IsNull() || edge.ShapeType() != TopAbs_EDGE)
919 {
586db386 920 di<<"A null shape or not an edge is used.\n";
c31e0def 921 return 2;
922 }
923
924 // Get face.
925 TopoDS_Shape face;
926 if (argc == 3)
927 {
928 const char* arg2 = argv[2];
929 face = DBRep::Get(arg2);
930 if (face.IsNull() || face.ShapeType() != TopAbs_FACE)
931 {
586db386 932 di<<"A null shape or not a face is used.\n";
c31e0def 933 return 3;
934 }
935 }
936
937 // Analysis of the edge.
938 ShapeAnalysis_Edge analyser;
939 Standard_Boolean isOk(Standard_True);
940
941 // Curve 3D.
942 if (analyser.HasCurve3d(TopoDS::Edge(edge)))
943 {
944 // Check vertices.
945 if (analyser.CheckVerticesWithCurve3d(TopoDS::Edge(edge)))
946 {
947 isOk = Standard_False;
948 di<<"Vertices of the edge don't coincide with start/end points of 3d-curve (using tolerance of the vertices).\n";
949 }
950 }
951 else
952 {
953 isOk = Standard_False;
954 di<<"Edge doesn't have a 3d-curve\n";
955 }
956
957 if (!face.IsNull())
958 {
959 // Curve 2D.
960 if (analyser.HasPCurve(TopoDS::Edge(edge), TopoDS::Face(face)))
961 {
962 // Check vertices.
963 if (analyser.CheckVerticesWithPCurve(TopoDS::Edge(edge), TopoDS::Face(face)))
964 {
965 isOk = Standard_False;
966 di<<"Vertices of the edge don't coincide with start/end points of 2d-curve (using tolerance of the vertices).\n";
967 }
968 }
969 else
970 {
971 isOk = Standard_False;
972 di<<"Edge doesn't have a 2d-curve on this face\n";
973 }
974 }
975
976 if (isOk)
977 di<<"Edge seems OK.\n";
978
979 return 0;
980}
981
7fd59977 982//=======================================================================
983//function : InitCommands
984//purpose :
985//=======================================================================
986
987 void SWDRAW_ShapeAnalysis::InitCommands(Draw_Interpretor& theCommands)
988{
989 static Standard_Integer initactor = 0;
c48e2889 990 if (initactor)
991 {
992 return;
993 }
994 initactor = 1;
995
7fd59977 996 Standard_CString g = SWDRAW::GroupName();
997
998 theCommands.Add ("tolerance","shape [tolmin tolmax:real]", __FILE__,tolerance,g);
999 theCommands.Add ("projface","nom_face X Y [Z]", __FILE__,projface,g);
1000 theCommands.Add ("projcurve","nom_edge | curve3d | curve3d first last + X Y Z",
1001 __FILE__,projcurve,g);
f6b08ecf 1002 theCommands.Add("projpcurve", "edge face tol x y z [start_param]",
1003 __FILE__, projpcurve, g);
1004 theCommands.Add("anaface", "nomface", __FILE__, anaface, g);
7fd59977 1005 theCommands.Add ("statshape","shape [particul] : stats/particularites",
1006 __FILE__,XSHAPE_statshape,g);
1007 theCommands.Add ("comptol","shape [nbpoints]",__FILE__,XSHAPE_comptoledge,g);
1008 theCommands.Add("freebounds",
1009 "shp toler [splitclosed [splitopen]] - free bounds; toler <= 0 for shells (no sewing call)",
1010 __FILE__, freebounds, g);
1011
1012 const char* groupold="DE: old";
1013 theCommands.Add("fbprops",
1014 "shp [toler [splitclosed [splitopen]]] - free bounds properties; toler <= 0 or not specified - for shells (no sewing call)",
1015 __FILE__, FreeBoundsProps, groupold);
1016 theCommands.Add("fbclose",
1017 "shp sewtoler closetoler [splitclosed [splitopen]] - closes free bounds; use sewtoler <= 0 for shells (no sewing call)",
1018 __FILE__, closefreebounds, groupold);
7fd59977 1019
1020 theCommands.Add("getareacontour","wire ",__FILE__, getareacontour, groupold);
ef57920a 1021 theCommands.Add ("checkselfintersection","wire [face]", __FILE__,checkselfintersection,g);
c31e0def 1022 theCommands.Add ("checkedge","edge [face]", __FILE__,checkedge,g);
7fd59977 1023}
Open CASCADE Technology repository