1 // Copyright (C) 2013-2016 CEA/DEN, EDF R&D, OPEN CASCADE
3 // This library is free software; you can redistribute it and/or
4 // modify it under the terms of the GNU Lesser General Public
5 // License as published by the Free Software Foundation; either
6 // version 2.1 of the License, or (at your option) any later version.
8 // This library is distributed in the hope that it will be useful,
9 // but WITHOUT ANY WARRANTY; without even the implied warranty of
10 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
11 // Lesser General Public License for more details.
13 // You should have received a copy of the GNU Lesser General Public
14 // License along with this library; if not, write to the Free Software
15 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
17 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
20 #include <Standard_Stream.hxx>
22 #include "AdvancedEngine_SmoothingSurfaceDriver.hxx"
23 #include "AdvancedEngine_ISmoothingSurface.hxx"
24 #include "AdvancedEngine_Types.hxx"
26 #include "GEOM_Function.hxx"
28 #include <StdFail_NotDone.hxx>
30 #include <TopoDS_Vertex.hxx>
31 #include <TopoDS_Face.hxx>
32 #include <TopoDS_Compound.hxx>
34 #include <TopTools_MapOfShape.hxx>
35 #include <TopTools_MapIteratorOfMapOfShape.hxx>
36 #include <TopExp_Explorer.hxx>
38 #include <TColgp_SequenceOfPnt.hxx>
39 #include <TColgp_Array2OfPnt.hxx>
40 #include <TColgp_SequenceOfXY.hxx>
41 #include <TColgp_SequenceOfXYZ.hxx>
42 #include <TColStd_Array1OfInteger.hxx>
44 #include <BRepAdaptor_HSurface.hxx>
46 #include <BRep_Builder.hxx>
47 #include <BRepGProp.hxx>
48 #include <BRep_Tool.hxx>
49 #include <BRepPrimAPI_MakeSphere.hxx>
50 #include <BRepPrimAPI_MakeBox.hxx>
51 #include <BRepBuilderAPI_MakeFace.hxx>
52 #include <BRepAlgoAPI_Cut.hxx>
54 #include <GeomPlate_Surface.hxx>
55 #include <GeomPlate_BuildPlateSurface.hxx>
56 #include <GeomPlate_PointConstraint.hxx>
57 #include <GeomPlate_MakeApprox.hxx>
58 #include <GeomPlate_PlateG0Criterion.hxx>
59 #include <GeomPlate_BuildAveragePlane.hxx>
61 #include <Geom_BSplineSurface.hxx>
62 #include <Geom_Surface.hxx>
63 #include <Geom_Plane.hxx>
65 #include <GProp_GProps.hxx>
66 #include <Bnd_Box.hxx>
67 #include <BRepBndLib.hxx>
75 #include <GC_MakePlane.hxx>
77 //=======================================================================
80 //=======================================================================
81 const Standard_GUID& AdvancedEngine_SmoothingSurfaceDriver::GetID()
83 static Standard_GUID aGUID("1C3A0F30-729D-4E83-8232-78E74FC5637C");
87 //=======================================================================
88 //function : AdvancedEngine_SmoothingSurfaceDriver
90 //=======================================================================
91 AdvancedEngine_SmoothingSurfaceDriver::AdvancedEngine_SmoothingSurfaceDriver()
95 //=======================================================================
96 //function : MakeSmoothingSurfaceUnClosed
98 //=======================================================================
99 TopoDS_Shape AdvancedEngine_SmoothingSurfaceDriver::MakeSmoothingSurfaceUnClosed
100 (const Handle(TColgp_HArray1OfPnt) &theListOfPoints,
101 const Standard_Integer theNbMax,
102 const Standard_Integer theDegMax,
103 const Standard_Real theDMax) const
105 TopoDS_Face aInitShape;
107 // Create an average Plane
108 //Handle(TColgp_HArray1OfPnt) HAP = new TColgp_HArray1OfPnt(1,theListOfPoints.Length())
109 GeomPlate_BuildAveragePlane gpbap(theListOfPoints,theListOfPoints->Length(),Precision::Confusion(),1,1);
110 Handle(Geom_Plane) plane(gpbap.Plane());
111 Standard_Real Umin, Umax, Vmin, Vmax;
112 gpbap.MinMaxBox(Umin,Umax,Vmin,Vmax);
113 // cout << "Vals : " << Umin << ", " << Umax << ", " << Vmin << ", " << Vmax << endl;
114 BRepBuilderAPI_MakeFace mf(plane,Umin,Umax,Vmin,Vmax,Precision::Confusion());
115 aInitShape = mf.Face();
118 GeomPlate_BuildPlateSurface aBuilder(3,10);
119 // ** Initialization of surface
120 Handle(BRepAdaptor_HSurface) HSI = new BRepAdaptor_HSurface();
121 HSI->ChangeSurface().Initialize(aInitShape);
122 aBuilder.LoadInitSurface( BRep_Tool::Surface(HSI->ChangeSurface().Face()));
124 Standard_Integer j, j1, j2;
125 // cout << "Init surface" << endl;
126 j1 = theListOfPoints->Lower();
127 j2 = theListOfPoints->Upper();
128 for (j=j1; j<=j2 ; j++)
130 gp_Pnt aPnt = theListOfPoints->Value(j);
131 Handle(GeomPlate_PointConstraint) PCont = new GeomPlate_PointConstraint(aPnt,0);
134 // cout << "avant Perform surface" << endl;
136 // cout << "Perform surface" << endl;
138 // A ce niveau : surface algo
139 Handle(GeomPlate_Surface) gpPlate = aBuilder.Surface();
141 Standard_Integer nbcarreau = (theNbMax > 0 ? theNbMax : 2);
142 Standard_Integer degmax = (theDegMax > 0 ? theDegMax : 8);
143 Standard_Real seuil = theDMax;
146 seuil = Max(0.0001,10*aBuilder.G0Error());
149 GeomPlate_MakeApprox Mapp(gpPlate,0.0001,nbcarreau,degmax,seuil);
150 // cout << "Approx surface" << endl;
152 Handle (Geom_Surface) Surf (Mapp.Surface());
154 aBuilder.Surface()->Bounds( Umin, Umax, Vmin, Vmax);
156 BRepBuilderAPI_MakeFace MF(Surf,Umin, Umax, Vmin, Vmax, Precision::Confusion());
157 TopoDS_Shape aShape = MF.Shape();
163 //=======================================================================
166 //=======================================================================
167 Standard_Integer AdvancedEngine_SmoothingSurfaceDriver::Execute(Handle(TFunction_Logbook)& log) const
169 if (Label().IsNull()) return 0;
170 Handle(GEOM_Function) aFunction = GEOM_Function::GetFunction(Label());
171 if (aFunction.IsNull()) return 0;
173 Standard_Integer aType = aFunction->GetType();
175 if (aType != SMOOTHINGSURFACE_LPOINTS) return 0;
176 // cout << "Youhou : " << aType << endl;
178 AdvancedEngine_ISmoothingSurface aData (aFunction);
180 // Fill the map of vertices.
181 Standard_Integer aNbShapes = aData.GetLength();
182 TopTools_MapOfShape aMapPoints;
185 for (i = 1; i <= aNbShapes; i++) {
186 Handle(GEOM_Function) aFShape = aData.GetPntOrComp(i);
187 TopoDS_Shape aShape = aFShape->GetValue();
189 if (aShape.ShapeType() == TopAbs_VERTEX) {
190 aMapPoints.Add(aShape);
192 TopExp_Explorer anExp(aShape, TopAbs_VERTEX);
194 for (; anExp.More(); anExp.Next()) {
195 aMapPoints.Add(anExp.Current());
200 // Add points to the array of points.
201 const Standard_Integer aNbPoints = aMapPoints.Extent();
202 Handle(TColgp_HArray1OfPnt) anArrayofPnt =
203 new TColgp_HArray1OfPnt(1, aNbPoints);
204 TopTools_MapIteratorOfMapOfShape anIter(aMapPoints);
206 for (i = 1; anIter.More(); anIter.Next(), i++) {
207 TopoDS_Vertex aPoint = TopoDS::Vertex(anIter.Key());
208 gp_Pnt aPnt = BRep_Tool::Pnt(aPoint);
210 anArrayofPnt->SetValue(i, aPnt);
213 const Standard_Integer aNbMax = aData.GetNbMax();
214 const Standard_Integer aDegMax = aData.GetDegMax();
215 const Standard_Real aDMax = aData.GetDMax();
217 // Make smoothing surface.
218 TopoDS_Shape aShape = AdvancedEngine_SmoothingSurfaceDriver::
219 MakeSmoothingSurfaceUnClosed(anArrayofPnt, aNbMax, aDegMax, aDMax);
221 if (aShape.IsNull()) return 0;
223 aFunction->SetValue(aShape);
225 log->SetTouched(Label());
230 //================================================================================
232 * \brief Returns a name of creation operation and names and values of creation parameters
234 //================================================================================
236 bool AdvancedEngine_SmoothingSurfaceDriver::
237 GetCreationInformation(std::string& theOperationName,
238 std::vector<GEOM_Param>& theParams)
240 if (Label().IsNull()) return 0;
241 Handle(GEOM_Function) function = GEOM_Function::GetFunction(Label());
243 AdvancedEngine_ISmoothingSurface aCI( function );
244 Standard_Integer aType = function->GetType();
246 theOperationName = "SMOOTHINGSURFACE";
249 case SMOOTHINGSURFACE_LPOINTS:
251 AddParam( theParams, "Points" );
252 if ( aCI.GetLength() > 1 )
253 theParams[0] << aCI.GetLength() << " points: ";
254 for ( int i = 1, nb = aCI.GetLength(); i <= nb; ++i )
255 theParams[0] << aCI.GetPntOrComp( i ) << " ";
257 const Standard_Integer aNbMax = aCI.GetNbMax();
258 const Standard_Integer aDegMax = aCI.GetDegMax();
259 const Standard_Real aDMax = aCI.GetDMax();
261 AddParam(theParams, "Max nbr of Bezier pieces", aCI.GetNbMax());
262 AddParam(theParams, "Max BSpline surface degree", aCI.GetDegMax());
263 AddParam(theParams, "3D tolerance of initial approximation", aCI.GetDMax());
273 IMPLEMENT_STANDARD_RTTIEXT (AdvancedEngine_SmoothingSurfaceDriver,GEOM_BaseDriver);