1 // Copyright (C) 2007-2019 CEA/DEN, EDF R&D, OPEN CASCADE
3 // Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
6 // This library is free software; you can redistribute it and/or
7 // modify it under the terms of the GNU Lesser General Public
8 // License as published by the Free Software Foundation; either
9 // version 2.1 of the License, or (at your option) any later version.
11 // This library is distributed in the hope that it will be useful,
12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 // Lesser General Public License for more details.
16 // You should have received a copy of the GNU Lesser General Public
17 // License along with this library; if not, write to the Free Software
18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
24 // Created: Tue Dec 7 11:59:05 2004
25 // Author: Pavel DURANDIN
27 #include <BlockFix.hxx>
29 #include <BlockFix_SphereSpaceModifier.hxx>
30 #include <BlockFix_PeriodicSurfaceModifier.hxx>
33 #include <TopExp_Explorer.hxx>
35 #include <TopLoc_Location.hxx>
38 #include <TopoDS_Edge.hxx>
39 #include <TopoDS_Face.hxx>
40 #include <TopoDS_Wire.hxx>
41 #include <TopoDS_Shape.hxx>
42 #include <TopoDS_Solid.hxx>
43 #include <TopoDS_Vertex.hxx>
45 #include <TopTools_ListOfShape.hxx>
46 #include <TopTools_ListIteratorOfListOfShape.hxx>
47 #include <TopTools_MapOfShape.hxx>
48 #include <TopTools_DataMapOfShapeShape.hxx>
49 #include <TopTools_DataMapIteratorOfDataMapOfShapeShape.hxx>
51 #include <BRep_Tool.hxx>
52 #include <BRep_Builder.hxx>
54 #include <BRepAdaptor_Surface.hxx>
56 #include <BRepTools.hxx>
57 #include <BRepTools_Modifier.hxx>
58 #include <BRepTools_Substitution.hxx>
60 #include <BRepOffsetAPI_MakeFilling.hxx>
62 #include <ShapeFix.hxx>
63 #include <ShapeFix_Edge.hxx>
64 #include <ShapeFix_Face.hxx>
66 #include <ShapeAnalysis.hxx>
67 #include <ShapeAnalysis_Edge.hxx>
68 #include <ShapeAnalysis_Curve.hxx>
69 #include <ShapeAnalysis_Surface.hxx>
71 #include <ShapeCustom.hxx>
73 #include <ShapeBuild_Edge.hxx>
74 #include <ShapeBuild_ReShape.hxx>
76 #include <ShapeFix_Wire.hxx>
78 #include <Geom_Surface.hxx>
79 #include <Geom_CylindricalSurface.hxx>
80 #include <Geom_ConicalSurface.hxx>
81 #include <Geom_SphericalSurface.hxx>
82 #include <Geom_ToroidalSurface.hxx>
84 #include <Geom2d_Curve.hxx>
86 #include <TColgp_SequenceOfPnt2d.hxx>
88 //=======================================================================
89 //function : FixResult
91 //=======================================================================
92 static void FixResult(const TopoDS_Shape& result,
93 Handle(ShapeBuild_ReShape)& Context,
94 const Standard_Real Tol)
96 for (TopExp_Explorer ex_f(result,TopAbs_FACE); ex_f.More(); ex_f.Next()) {
97 TopoDS_Shape aShape = Context->Apply(ex_f.Current().Oriented(TopAbs_FORWARD));
98 // face could not be dropped or split on this step
99 TopoDS_Face aFace = TopoDS::Face(aShape);
101 Handle(Geom_Surface) Surf = BRep_Tool::Surface(aFace,L);
103 if( Surf->IsKind(STANDARD_TYPE(Geom_SphericalSurface)) ||
104 Surf->IsKind(STANDARD_TYPE(Geom_CylindricalSurface)) ) {
106 Standard_Integer nbWires = 0;
107 for (TopExp_Explorer ex_w(aFace,TopAbs_WIRE); ex_w.More(); ex_w.Next()) {
109 Handle(ShapeFix_Wire) sfw = new ShapeFix_Wire(TopoDS::Wire(ex_w.Current()),
111 Precision::Confusion());
113 if(sfw->StatusReorder ( ShapeExtend_FAIL ))
116 sfw->SetPrecision(2.*Tol);
119 Standard_Boolean isDone = sfw->LastFixStatus ( ShapeExtend_DONE );
120 isDone |= sfw->FixDegenerated();
122 // remove degenerated edges from not degenerated points
123 ShapeAnalysis_Edge sae;
124 Handle(ShapeExtend_WireData) sewd = sfw->WireData();
126 for( i = 1; i<=sewd->NbEdges();i++) {
127 TopoDS_Edge E = sewd->Edge(i);
128 if(BRep_Tool::Degenerated(E)&&!sae.HasPCurve(E,aFace)) {
130 isDone = Standard_True;
135 //isDone |= sfw->FixLacking(); // commented by skl 22.03.2005 (PAL8395)
137 // remove neighbour seam edges
139 for( i = 1; i<sewd->NbEdges();i++) {
140 if(sewd->IsSeam(i) && sewd->IsSeam(i+1)) {
141 isDone = Standard_True;
147 if(sewd->IsSeam(1) && sewd->IsSeam(sewd->NbEdges())) {
149 sewd->Remove(sewd->NbEdges());
154 TopoDS_Wire ResWire = sfw->Wire();
155 Context->Replace(ex_w.Current(), ResWire);
158 // Implement fix orientation in case of several wires
160 TopoDS_Face aFixedFace = TopoDS::Face(Context->Apply(aFace));
161 Handle(ShapeFix_Face) sff = new ShapeFix_Face(aFixedFace);
162 if(sff->FixOrientation())
163 Context->Replace(aFixedFace,sff->Face());
169 //=======================================================================
170 //function : RotateSphereSpace
172 //=======================================================================
173 TopoDS_Shape BlockFix::RotateSphereSpace (const TopoDS_Shape& S,
174 const Standard_Real Tol)
176 // Create a modification description
177 Handle(BlockFix_SphereSpaceModifier) SR = new BlockFix_SphereSpaceModifier;
178 SR->SetTolerance(Tol);
180 TopTools_DataMapOfShapeShape context;
181 BRepTools_Modifier MD;
182 TopoDS_Shape result = ShapeCustom::ApplyModifier ( S, SR, context,MD );
184 Handle(ShapeBuild_ReShape) RS = new ShapeBuild_ReShape;
185 FixResult(result,RS,Tol);
186 result = RS->Apply(result);
189 for(TopExp_Explorer exp(result,TopAbs_EDGE); exp.More(); exp.Next()) {
190 TopoDS_Edge E = TopoDS::Edge(exp.Current());
191 sfe.FixVertexTolerance (E);
194 ShapeFix::SameParameter(result,Standard_False);
198 //=======================================================================
199 //function : RefillProblemFaces
201 //=======================================================================
202 TopoDS_Shape BlockFix::RefillProblemFaces (const TopoDS_Shape& aShape)
204 Standard_Integer NbSamples = 10;
206 TopTools_ListOfShape theFaces;
208 TopExp_Explorer Explo(aShape, TopAbs_FACE);
209 for (; Explo.More(); Explo.Next())
211 TopoDS_Face aFace = TopoDS::Face(Explo.Current());
212 BRepAdaptor_Surface BAsurf(aFace);
213 GeomAbs_SurfaceType SurfType = BAsurf.GetType();
214 if (SurfType >= GeomAbs_BezierSurface)
216 TopExp_Explorer fexp(aFace, TopAbs_EDGE);
217 for (; fexp.More(); fexp.Next())
219 const TopoDS_Edge& anEdge = TopoDS::Edge(fexp.Current());
220 if (BRep_Tool::Degenerated(anEdge))
222 TopoDS_Vertex V1, V2;
223 TopExp::Vertices(anEdge, V1, V2);
226 gp_Pnt aPnt = BRep_Tool::Pnt(V1);
227 Standard_Real TolV = BRep_Tool::Tolerance(V1);
228 Handle(Geom_Surface) aSurf = BRep_Tool::Surface(aFace);
229 Handle(ShapeAnalysis_Surface) Analyser = new ShapeAnalysis_Surface(aSurf);
230 if (Analyser->IsDegenerated(aPnt, TolV))
232 theFaces.Append(aFace);
241 //Now all problem faces are collected in the list "theFaces"
242 BRepTools_Substitution aSubst;
243 TopTools_ListIteratorOfListOfShape itl(theFaces);
244 for (; itl.More(); itl.Next())
246 const TopoDS_Face& aFace = TopoDS::Face(itl.Value());
247 BRepOffsetAPI_MakeFilling Filler;
248 for (Explo.Init(aFace, TopAbs_EDGE); Explo.More(); Explo.Next())
250 const TopoDS_Edge& anEdge = TopoDS::Edge(Explo.Current());
251 if (!BRep_Tool::Degenerated(anEdge))
252 Filler.Add(anEdge, GeomAbs_C0);
254 Standard_Real Umin, Umax, Vmin, Vmax;
255 BRepTools::UVBounds(aFace, Umin, Umax, Vmin, Vmax);
256 //Handle(Geom_Surface) aSurf = BRep_Tool::Surface(aFace);
257 Standard_Integer i, j;
258 for (i = 1; i < NbSamples; i++)
259 for (j = 1; j < NbSamples; j++) {
261 gp_Pnt aPoint = aSurf->Value(Umin + i*(Umax-Umin)/NbSamples,
262 Vmin + j*(Vmax-Vmin)/NbSamples);
265 Filler.Add(Umin + i*(Umax-Umin)/NbSamples,
266 Vmin + j*(Vmax-Vmin)/NbSamples,
273 for (Explo.Init(aFace, TopAbs_EDGE); Explo.More(); Explo.Next())
275 const TopoDS_Edge& anEdge = TopoDS::Edge(Explo.Current());
276 TopTools_ListOfShape Ledge;
277 if (!BRep_Tool::Degenerated(anEdge))
279 const TopTools_ListOfShape& Ledges = Filler.Generated(anEdge);
280 if (!Ledges.IsEmpty()) {
281 TopoDS_Shape NewEdge = Ledges.First();
282 Ledge.Append(NewEdge.Oriented(TopAbs_FORWARD));
285 aSubst.Substitute(anEdge, Ledge);
287 TopTools_ListOfShape Lface;
288 TopoDS_Face NewFace = TopoDS::Face(Filler.Shape());
289 NewFace.Orientation(TopAbs_FORWARD);
290 BRepAdaptor_Surface NewBAsurf(NewFace);
292 gp_Vec D1U, D1V, Normal, NewNormal;
293 Handle(Geom_Surface) aSurf = BRep_Tool::Surface(aFace);
294 aSurf->D1((Umin+Umax)*0.5, (Vmin+Vmax)*0.5, MidPnt, D1U, D1V);
296 NewBAsurf.D1((NewBAsurf.FirstUParameter() + NewBAsurf.LastUParameter())*0.5,
297 (NewBAsurf.FirstVParameter() + NewBAsurf.LastVParameter())*0.5,
299 NewNormal = D1U ^ D1V;
300 if (Normal * NewNormal < 0.)
302 Lface.Append(NewFace);
303 aSubst.Substitute(aFace, Lface);
306 aSubst.Build(aShape);
308 TopoDS_Shape Result = aShape;
309 if (aSubst.IsCopied(aShape))
310 Result = aSubst.Copy(aShape).First();
312 BRepTools::RemoveUnusedPCurves(Result);
317 //=======================================================================
318 //function : FixRanges
320 //=======================================================================
321 TopoDS_Shape BlockFix::FixRanges (const TopoDS_Shape& S,
322 const Standard_Real Tol)
324 // Create a modification description
325 Handle(BlockFix_PeriodicSurfaceModifier) SR = new BlockFix_PeriodicSurfaceModifier;
326 SR->SetTolerance(Tol);
328 TopTools_DataMapOfShapeShape context;
329 BRepTools_Modifier MD;
330 TopoDS_Shape result = ShapeCustom::ApplyModifier ( S, SR, context,MD );
332 Handle(ShapeBuild_ReShape) RS = new ShapeBuild_ReShape;
333 FixResult(result,RS,Tol);
334 result = RS->Apply(result);
337 for(TopExp_Explorer exp(result,TopAbs_EDGE); exp.More(); exp.Next()) {
338 TopoDS_Edge E = TopoDS::Edge(exp.Current());
339 sfe.FixVertexTolerance (E);
342 ShapeFix::SameParameter(result,Standard_False);