/src/contrib/geom-5.1.2.7/src/BlockFix/BlockFix.cpp

  1//  Copyright (C) 2007-2008  CEA/DEN, EDF R&D, OPEN CASCADE
  2//
  3//  Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
  4//  CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
  5//
  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.
 10//
 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.
 15//
 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
 19//
 20//  See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
 21//
 22// File:	BlockFix.cxx
 23// Created:	Tue Dec  7 11:59:05 2004
 24// Author:	Pavel DURANDIN
 25//
 26#include <BlockFix.hxx>
 27#include <TopoDS_Shape.hxx>
 28#include <TopTools_DataMapOfShapeShape.hxx>
 29#include <ShapeCustom.hxx>
 30#include <BRepTools.hxx>
 31#include <ShapeBuild_ReShape.hxx>
 32#include <TopoDS_Face.hxx>
 33#include <TopExp_Explorer.hxx>
 34#include <TopoDS.hxx>
 35#include <TopLoc_Location.hxx>
 36#include <Geom_Surface.hxx>
 37#include <Geom_CylindricalSurface.hxx>
 38#include <Geom_ConicalSurface.hxx>
 39#include <ShapeFix_Wire.hxx>
 40#include <TopoDS_Wire.hxx>
 41#include <BRepTools_Modifier.hxx>
 42#include <Geom_SphericalSurface.hxx>
 43#include <Geom_ToroidalSurface.hxx>
 44#include <BRep_Tool.hxx>
 45#include <TopoDS_Edge.hxx>
 46#include <Geom2d_Curve.hxx>
 47#include <BRep_Builder.hxx>
 48#include <ShapeAnalysis_Edge.hxx>
 49#include <ShapeFix_Edge.hxx>
 50#include <ShapeFix.hxx>
 51#include <ShapeFix_Face.hxx>
 52#include <ShapeAnalysis.hxx>
 53
 54#include <TColgp_SequenceOfPnt2d.hxx>
 55#include <ShapeAnalysis_Curve.hxx>
 56#include <TopoDS_Vertex.hxx>
 57#include <ShapeBuild_Edge.hxx>
 58
 59#include <BlockFix_SphereSpaceModifier.hxx>
 60#include <TopTools_DataMapIteratorOfDataMapOfShapeShape.hxx>
 61#include <TopTools_MapOfShape.hxx>
 62#include <BlockFix_PeriodicSurfaceModifier.hxx>
 63
 64#include <TopoDS_Solid.hxx>
 65
 66
 67//=======================================================================
 68//function : FixResult
 69//purpose  : auxilary
 70//=======================================================================
 71static void FixResult(const TopoDS_Shape& result,
 72                      Handle(ShapeBuild_ReShape)& Context,
 73                      const Standard_Real Tol)
 74{
 75  for (TopExp_Explorer ex_f(result,TopAbs_FACE); ex_f.More(); ex_f.Next()) {
 76    TopoDS_Shape aShape = Context->Apply(ex_f.Current().Oriented(TopAbs_FORWARD));
 77    // face coud not be dropped or splitted on this step
 78    TopoDS_Face aFace = TopoDS::Face(aShape);
 79    TopLoc_Location L;
 80    Handle(Geom_Surface) Surf = BRep_Tool::Surface(aFace,L);
 81    
 82    if( Surf->IsKind(STANDARD_TYPE(Geom_SphericalSurface)) ||
 83        Surf->IsKind(STANDARD_TYPE(Geom_CylindricalSurface)) ) {
 84    
 85      Standard_Integer nbWires = 0;
 86      for (TopExp_Explorer ex_w(aFace,TopAbs_WIRE); ex_w.More(); ex_w.Next()) {
 87        nbWires++;
 88        Handle(ShapeFix_Wire) sfw = new ShapeFix_Wire(TopoDS::Wire(ex_w.Current()), 
 89                                                      aFace, 
 90                                                      Precision::Confusion());
 91        sfw->FixReorder();
 92        if(sfw->StatusReorder ( ShapeExtend_FAIL ))
 93          continue;
 94        
 95        sfw->SetPrecision(2.*Tol);
 96        sfw->FixShifted();
 97        
 98        Standard_Boolean isDone = sfw->LastFixStatus ( ShapeExtend_DONE );
 99        isDone |= sfw->FixDegenerated();
100        
101        // remove degenerated edges from not degenerated points
102        ShapeAnalysis_Edge sae;
103        Handle(ShapeExtend_WireData) sewd = sfw->WireData();
104        Standard_Integer i;
105        for( i = 1; i<=sewd->NbEdges();i++) {
106          TopoDS_Edge E = sewd->Edge(i);
107          if(BRep_Tool::Degenerated(E)&&!sae.HasPCurve(E,aFace)) {
108            sewd->Remove(i);
109            isDone = Standard_True;
110            i--;
111          }
112        }
113        
114        //isDone |= sfw->FixLacking(); // commented by skl 22.03.2005 (PAL8395)
115        
116        // remove neighbour seam edges 
117        if(isDone) {
118          for( i = 1; i<sewd->NbEdges();i++) {
119            if(sewd->IsSeam(i) && sewd->IsSeam(i+1)) {
120              isDone = Standard_True;
121              sewd->Remove(i);
122              sewd->Remove(i);
123              i--;
124            }
125          }
126          if(sewd->IsSeam(1) && sewd->IsSeam(sewd->NbEdges())) {
127            sewd->Remove(1);
128            sewd->Remove(sewd->NbEdges());
129          }
130        }
131        
132               
133        if(isDone) {
134          TopoDS_Wire ResWire = sfw->Wire();
135          Context->Replace(ex_w.Current(), ResWire);
136        };
137      }
138      // Implement fix orientation in case of several wires
139      if(nbWires > 1) {
140        TopoDS_Face aFixedFace = TopoDS::Face(Context->Apply(aFace));
141        Handle(ShapeFix_Face) sff = new ShapeFix_Face(aFixedFace);
142        if(sff->FixOrientation())
143          Context->Replace(aFixedFace,sff->Face());
144      }
145      
146    }
147  }
148}
149
150
151
152
153
154//=======================================================================
155//function : ConvertToAnalytical
156//purpose  : 
157//=======================================================================
158
159TopoDS_Shape BlockFix::RotateSphereSpace (const TopoDS_Shape& S,
160                                          const Standard_Real Tol) 
161{
162
163  // Create a modification description
164  Handle(BlockFix_SphereSpaceModifier) SR = new BlockFix_SphereSpaceModifier;
165  SR->SetTolerance(Tol);
166
167  TopTools_DataMapOfShapeShape context;
168  BRepTools_Modifier MD;
169  TopoDS_Shape result = ShapeCustom::ApplyModifier ( S, SR, context,MD );
170  
171  Handle(ShapeBuild_ReShape) RS = new ShapeBuild_ReShape;
172  FixResult(result,RS,Tol);
173  result = RS->Apply(result);
174  
175  ShapeFix_Edge sfe;
176  for(TopExp_Explorer exp(result,TopAbs_EDGE); exp.More(); exp.Next()) {
177    TopoDS_Edge E = TopoDS::Edge(exp.Current());
178    sfe.FixVertexTolerance (E);
179  }
180  
181  ShapeFix::SameParameter(result,Standard_False);
182  return result;
183}
184
185
186//=======================================================================
187//function : FixRanges
188//purpose  : 
189//=======================================================================
190
191TopoDS_Shape BlockFix::FixRanges (const TopoDS_Shape& S,
192                                  const Standard_Real Tol) 
193{
194  // Create a modification description
195  Handle(BlockFix_PeriodicSurfaceModifier) SR = new BlockFix_PeriodicSurfaceModifier;
196  SR->SetTolerance(Tol);
197
198  TopTools_DataMapOfShapeShape context;
199  BRepTools_Modifier MD;
200  TopoDS_Shape result = ShapeCustom::ApplyModifier ( S, SR, context,MD );
201  
202  Handle(ShapeBuild_ReShape) RS = new ShapeBuild_ReShape;
203  FixResult(result,RS,Tol);
204  result = RS->Apply(result);
205  
206  ShapeFix_Edge sfe;
207  for(TopExp_Explorer exp(result,TopAbs_EDGE); exp.More(); exp.Next()) {
208    TopoDS_Edge E = TopoDS::Edge(exp.Current());
209    sfe.FixVertexTolerance (E);
210  }
211  
212  ShapeFix::SameParameter(result,Standard_False);
213
214  return result;
215}