1 // Copyright (c) 1998-1999 Matra Datavision
2 // Copyright (c) 1999-2014 OPEN CASCADE SAS
4 // This file is part of Open CASCADE Technology software library.
6 // This library is free software; you can redistribute it and/or modify it under
7 // the terms of the GNU Lesser General Public License version 2.1 as published
8 // by the Free Software Foundation, with special exception defined in the file
9 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
10 // distribution for complete text of the license and disclaimer of any warranty.
12 // Alternatively, this file may be used under the terms of Open CASCADE
13 // commercial license or contractual agreement.
16 #include <Bnd_Box.hxx>
17 #include <BRep_Builder.hxx>
18 #include <BRep_Tool.hxx>
19 #include <BRepAlgo_Section.hxx>
20 #include <BRepBndLib.hxx>
21 #include <BRepFill_Generator.hxx>
22 #include <BRepLib_MakeEdge.hxx>
23 #include <BRepLib_MakeVertex.hxx>
24 #include <BRepLib_MakeWire.hxx>
25 #include <BRepProj_Projection.hxx>
26 #include <BRepSweep_Prism.hxx>
27 #include <BRepTools_Modifier.hxx>
28 #include <BRepTools_TrsfModification.hxx>
31 #include <gp_Trsf.hxx>
33 #include <Precision.hxx>
34 #include <ShapeAnalysis_FreeBounds.hxx>
35 #include <Standard_ConstructionError.hxx>
36 #include <Standard_NoSuchObject.hxx>
37 #include <Standard_NullObject.hxx>
39 #include <TopExp_Explorer.hxx>
40 #include <TopLoc_Location.hxx>
42 #include <TopoDS_Iterator.hxx>
43 #include <TopoDS_Shape.hxx>
44 #include <TopTools_ListIteratorOfListOfShape.hxx>
45 #include <TopTools_ListOfShape.hxx>
47 //=======================================================================
48 //function : DistanceOut
49 //purpose : Compute the minimum distance between input shapes
50 // (using Bounding Boxes of each Shape)
51 //=======================================================================
52 static Standard_Real DistanceOut (const TopoDS_Shape& S1, const TopoDS_Shape& S2)
55 BRepBndLib::Add(S1,BBox1);
56 BRepBndLib::Add(S2,BBox2);
57 return BBox1.Distance(BBox2);
60 //=======================================================================
61 //function : DistanceIn
62 //purpose : Compute the maximum distance between input Shapes
63 // we compute the maximum dimension of each Bounding Box and then
64 // add each other with the minimum distance of shapes.
65 //=======================================================================
67 static Standard_Real DistanceIn (const TopoDS_Shape& S1, const TopoDS_Shape& S2)
70 BRepBndLib::Add(S1,SBBox);
71 BRepBndLib::Add(S2,LBBox);
73 Standard_Real LXmin, LYmin, LZmin, LXmax, LYmax, LZmax,
74 SXmin, SYmin, SZmin, SXmax, SYmax, SZmax;
75 SBBox.Get(SXmin, SYmin, SZmin,
77 LBBox.Get(LXmin, LYmin, LZmin,
80 //Compute the max distance between input shapes------------//
81 gp_XYZ Lmin(LXmin, LYmin, LZmin),
82 Lmax(LXmax, LYmax, LZmax);
83 gp_XYZ Smin(SXmin, SYmin, SZmin),
84 Smax(SXmax, SYmax, SZmax);
87 return Lmax.Modulus() + Smax.Modulus() + DistanceOut(S1, S2);
90 //=======================================================================
91 //function : BuildSection
92 //purpose : Cuts theShape by theTool using BRepAlgoAPI_Section and
93 // stores result as set of connected wires and compound
94 //=======================================================================
96 void BRepProj_Projection::BuildSection (const TopoDS_Shape& theShape,
97 const TopoDS_Shape& theTool)
99 myIsDone = Standard_False;
104 // if theShape is compound, extract only faces -- section algorithm
105 // may refuse to work if e.g. vertex is present
107 if (theShape.ShapeType() == TopAbs_FACE ||
108 theShape.ShapeType() == TopAbs_SHELL ||
109 theShape.ShapeType() == TopAbs_SOLID ||
110 theShape.ShapeType() == TopAbs_COMPSOLID)
112 else if (theShape.ShapeType() == TopAbs_COMPOUND)
116 TopExp_Explorer exp (theShape, TopAbs_FACE);
117 for (; exp.More(); exp.Next())
121 B.Add (C, exp.Current());
125 if ( aShape.IsNull() )
126 throw Standard_ConstructionError(__FILE__": target shape has no faces");
128 // build section computing pcurves on the shape
129 // BRepAlgoAPI_Section aSectionTool (aShape, theTool, Standard_False);
130 BRepAlgo_Section aSectionTool (aShape, theTool, Standard_False);
131 aSectionTool.Approximation (Standard_True);
132 aSectionTool.ComputePCurveOn1 (Standard_True);
133 aSectionTool.Build();
135 // check for successful work of the section tool
136 if (! aSectionTool.IsDone())
139 // get edges of the result
140 Handle(TopTools_HSequenceOfShape) anEdges = new TopTools_HSequenceOfShape;
141 TopExp_Explorer exp(aSectionTool.Shape(), TopAbs_EDGE);
142 for (; exp.More(); exp.Next())
143 anEdges->Append (exp.Current());
145 // if no edges are found, this means that this section yields no result
146 if (anEdges->Length() <= 0)
149 // connect edges to wires using ShapeAnalysis functionality
150 ShapeAnalysis_FreeBounds::ConnectEdgesToWires (anEdges, Precision::Confusion(),
151 Standard_True, mySection);
152 myIsDone = (! mySection.IsNull() && mySection->Length() > 0);
154 // collect all resulting wires to compound
158 B.MakeCompound (myShape);
159 for (Standard_Integer i=1; i <= mySection->Length(); i++)
160 B.Add (myShape, mySection->Value(i));
162 // initialize iteration (for compatibility with previous versions)
167 //=======================================================================
168 //function : BRepProj_Projection
169 //purpose : Cylindrical Projection
170 //=======================================================================
172 BRepProj_Projection::BRepProj_Projection(const TopoDS_Shape& Wire,
173 const TopoDS_Shape& Shape,
175 : myIsDone(Standard_False), myItr(0)
178 Standard_NullObject_Raise_if((Wire.IsNull() || Shape.IsNull()),__FILE__": null input shape");
179 if (Wire.ShapeType() != TopAbs_EDGE &&
180 Wire.ShapeType() != TopAbs_WIRE )
181 throw Standard_ConstructionError(__FILE__": projected shape is neither wire nor edge");
183 // compute the "length" of the cylindrical surface to build
184 Standard_Real mdis = DistanceIn(Wire, Shape);
185 gp_Vec Vsup (D.XYZ() * 2 * mdis);
186 gp_Vec Vinf (D.XYZ() * -mdis);
188 // move the base of the cylindrical surface by translating it by -mdis
190 T.SetTranslation(Vinf);
191 // Note: it is necessary to create copy of wire to avoid adding new pcurves into it
192 Handle(BRepTools_TrsfModification) Trsf = new BRepTools_TrsfModification(T);
193 BRepTools_Modifier Modif (Wire, Trsf);
194 TopoDS_Shape WireBase = Modif.ModifiedShape(Wire);
196 // Creation of a cylindrical surface
197 BRepSweep_Prism CylSurf (WireBase, Vsup, Standard_False);
200 BuildSection (Shape, CylSurf.Shape());
203 //=======================================================================
204 //function : BRepProj_Projection
205 //purpose : Conical projection
206 //=======================================================================
208 BRepProj_Projection::BRepProj_Projection (const TopoDS_Shape& Wire,
209 const TopoDS_Shape& Shape,
211 : myIsDone(Standard_False), myItr(0)
214 Standard_NullObject_Raise_if((Wire.IsNull() || Shape.IsNull()),__FILE__": null input shape");
215 if (Wire.ShapeType() != TopAbs_EDGE &&
216 Wire.ShapeType() != TopAbs_WIRE )
217 throw Standard_ConstructionError(__FILE__": projected shape is neither wire nor edge");
219 // if Wire is only an edge, transform it into a Wire
221 if (Wire.ShapeType() == TopAbs_EDGE)
228 aWire = TopoDS::Wire(Wire);
230 // compute the "length" of the conical surface to build
231 Standard_Real mdis = DistanceIn(Wire, Shape);
233 // Initialize iterator to get first sub-shape of Wire
234 TopExp_Explorer ExpWire;
235 ExpWire.Init (aWire, TopAbs_VERTEX);
237 // get the first Point of the first sub-shape os the Wire
238 gp_Pnt PC = BRep_Tool::Pnt(TopoDS::Vertex(ExpWire.Current()));
240 // compute the ratio of the scale transformation
241 Standard_Real Scale = PC.Distance(P);
242 if ( Abs (Scale) < Precision::Confusion() )
243 throw Standard_ConstructionError("Projection");
244 Scale = 1. + mdis / Scale;
246 // move the base of the conical surface by scaling it with ratio Scale
248 T.SetScale(P, Scale);
249 Handle(BRepTools_TrsfModification) Tsca = new BRepTools_TrsfModification(T);
250 BRepTools_Modifier ModifScale(aWire,Tsca);
251 TopoDS_Shape ShapeGen1 = ModifScale.ModifiedShape(aWire);
253 TopoDS_Vertex aVertex = BRepLib_MakeVertex(P);
256 BB.MakeEdge( DegEdge );
257 BB.Add( DegEdge, aVertex.Oriented(TopAbs_FORWARD) );
258 BB.Add( DegEdge, aVertex.Oriented(TopAbs_REVERSED) );
259 BB.Degenerated( DegEdge, Standard_True );
262 BB.MakeWire( DegWire );
263 BB.Add( DegWire, DegEdge );
264 DegWire.Closed( Standard_True );
266 // Build the Ruled surface based shape
267 BRepFill_Generator RuledSurf;
268 RuledSurf.AddWire(DegWire);
269 RuledSurf.AddWire(TopoDS::Wire(ShapeGen1));
271 TopoDS_Shell SurfShell = RuledSurf.Shell();
274 BuildSection (Shape, SurfShell);