b311480e |
1 | // Created on: 1995-09-22 |
2 | // Created by: Bruno DUMORTIER |
3 | // Copyright (c) 1995-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. |
7fd59977 |
16 | |
7fd59977 |
17 | |
7fd59977 |
18 | #include <AppDef_Compute.hxx> |
19 | #include <AppDef_MultiLine.hxx> |
20 | #include <AppDef_MultiPointConstraint.hxx> |
42cf5bc1 |
21 | #include <AppParCurves_MultiCurve.hxx> |
22 | #include <BRepFill_ApproxSeewing.hxx> |
7fd59977 |
23 | #include <BSplCLib.hxx> |
7fd59977 |
24 | #include <Geom2d_BSplineCurve.hxx> |
42cf5bc1 |
25 | #include <Geom2d_Curve.hxx> |
26 | #include <Geom_BSplineCurve.hxx> |
27 | #include <Geom_Curve.hxx> |
28 | #include <PLib.hxx> |
29 | #include <StdFail_NotDone.hxx> |
7fd59977 |
30 | #include <TColgp_Array1OfPnt.hxx> |
31 | #include <TColgp_Array1OfPnt2d.hxx> |
7fd59977 |
32 | #include <TColStd_Array1OfInteger.hxx> |
42cf5bc1 |
33 | #include <TColStd_Array1OfReal.hxx> |
7fd59977 |
34 | |
35 | //======================================================================= |
36 | //function : BRepFill_ApproxSeewing |
37 | //purpose : |
38 | //======================================================================= |
7fd59977 |
39 | BRepFill_ApproxSeewing::BRepFill_ApproxSeewing() |
40 | :myIsDone(Standard_False) |
41 | { |
42 | } |
43 | |
44 | |
45 | //======================================================================= |
46 | //function : BRepFill_ApproxSeewing |
47 | //purpose : |
48 | //======================================================================= |
49 | |
50 | BRepFill_ApproxSeewing::BRepFill_ApproxSeewing(const BRepFill_MultiLine& ML) |
51 | :myIsDone(Standard_False) |
52 | { |
53 | Perform(ML); |
54 | } |
55 | |
56 | |
57 | //======================================================================= |
58 | //function : Perform |
59 | //purpose : |
60 | //======================================================================= |
61 | |
62 | void BRepFill_ApproxSeewing::Perform(const BRepFill_MultiLine& ML) |
63 | { |
64 | myML = ML; |
65 | |
66 | // evaluate the approximative length of the 3dCurve |
67 | Standard_Integer i; |
68 | Standard_Real Length = 0.; |
69 | Standard_Real U1 = myML.FirstParameter(); |
70 | Standard_Real U2 = myML.LastParameter(); |
71 | Standard_Integer NbPoints = 50; |
72 | Standard_Real Dist, dU = (U2 - U1) / ( 2*NbPoints - 1); |
73 | |
74 | TColgp_Array1OfPnt2d LP(1,2*NbPoints); // tableau Longueur <-> Param |
51740958 |
75 | gp_Pnt aPnt1, aPnt2; |
76 | aPnt1 = myML.Value(U1); |
7fd59977 |
77 | |
78 | for ( i = 0; i < 2*NbPoints ; i++) { |
51740958 |
79 | aPnt2 = myML.Value(U1 + i*dU); |
80 | Dist = aPnt1.Distance(aPnt2); |
7fd59977 |
81 | Length += Dist; |
82 | LP(i+1) = gp_Pnt2d( Length, U1 + (i*dU)); |
51740958 |
83 | aPnt1 = aPnt2; |
7fd59977 |
84 | } |
85 | |
86 | // On cherche a mettre NbPoints dans la curve. |
87 | // on met les points environ a Length/NbPoints. |
88 | |
89 | AppDef_MultiLine MLS ( NbPoints); |
90 | AppDef_MultiPointConstraint MP ( 1, 2); |
91 | gp_Pnt P3d; |
92 | gp_Pnt2d PF1,PF2; |
93 | |
94 | ML.Value3dOnF1OnF2(U1,P3d,PF1,PF2); |
95 | MP.SetPoint (1, P3d); |
96 | MP.SetPoint2d(2, PF1); |
97 | MP.SetPoint2d(3, PF2); |
98 | MLS.SetValue (1, MP); |
99 | |
100 | Standard_Real DCorde = Length / ( NbPoints - 1); |
101 | Standard_Real Corde = DCorde; |
102 | Standard_Integer Index = 1; |
103 | Standard_Real U, Alpha; |
104 | for ( i = 2; i < NbPoints; i++) { |
105 | while ( LP(Index).X() < Corde) Index ++; |
106 | Alpha = (Corde - LP(Index-1).X()) / (LP(Index).X() - LP(Index-1).X()); |
107 | U = LP(Index-1).Y() + Alpha * ( LP(Index).Y() - LP(Index-1).Y()); |
108 | AppDef_MultiPointConstraint MPC( 1, 2); |
109 | ML.Value3dOnF1OnF2(U,P3d,PF1,PF2); |
110 | MPC.SetPoint (1, P3d); |
111 | MPC.SetPoint2d(2, PF1); |
112 | MPC.SetPoint2d(3, PF2); |
113 | MLS.SetValue (i, MPC); |
114 | Corde = i*DCorde; |
115 | } |
116 | AppDef_MultiPointConstraint MPE( 1, 2); |
117 | ML.Value3dOnF1OnF2(U2,P3d,PF1,PF2); |
118 | MPE.SetPoint (1, P3d); |
119 | MPE.SetPoint2d(2, PF1); |
120 | MPE.SetPoint2d(3, PF2); |
121 | MLS.SetValue (NbPoints, MPE); |
122 | |
123 | AppDef_Compute Fit(MLS); |
124 | |
125 | Standard_Integer NbCurves = Fit.NbMultiCurves(); |
126 | // Standard_Integer MaxDeg = 0; |
127 | |
128 | if ( NbCurves == 0) { |
0797d9d3 |
129 | #ifdef OCCT_DEBUG |
7fd59977 |
130 | cout << " TrimSurfaceTool : Approx echoue, on met les polygones" << endl; |
131 | #endif |
132 | |
133 | TColStd_Array1OfReal Knots(1,NbPoints); |
134 | TColStd_Array1OfInteger Mults(1,NbPoints); |
135 | Mults.Init(1); |
136 | Mults(1) = Mults(NbPoints) = 2; |
137 | TColgp_Array1OfPnt P (1,NbPoints); |
138 | TColgp_Array1OfPnt2d P1(1,NbPoints); |
139 | TColgp_Array1OfPnt2d P2(1,NbPoints); |
140 | |
141 | Standard_Real Uf = ML.FirstParameter(); |
142 | Standard_Real Ul = ML.LastParameter(); |
51740958 |
143 | Standard_Real dUlf = (Ul-Uf)/(NbPoints-1); |
7fd59977 |
144 | AppDef_MultiPointConstraint MPC; |
145 | for ( i = 1; i<= NbPoints-1; i++) { |
146 | MPC = MLS.Value(i); |
51740958 |
147 | U = Uf + (i-1) * dUlf; |
7fd59977 |
148 | P (i) = MPC.Point(1); |
149 | P1(i) = MPC.Point2d(2); |
150 | P2(i) = MPC.Point2d(3); |
151 | Knots(i) = U; |
152 | } |
153 | // eval the last point on Ul |
154 | MPC = MLS.Value(NbPoints); |
155 | P (NbPoints) = MPC.Point(1); |
156 | P1(NbPoints) = MPC.Point2d(2); |
157 | P2(NbPoints) = MPC.Point2d(3); |
158 | Knots(NbPoints) = Ul; |
159 | |
160 | myCurve = new Geom_BSplineCurve ( P , Knots, Mults, 1); |
161 | myPCurve1 = new Geom2d_BSplineCurve( P1, Knots, Mults, 1); |
162 | myPCurve2 = new Geom2d_BSplineCurve( P2, Knots, Mults, 1); |
163 | |
164 | myIsDone = Standard_True; |
165 | |
166 | return; |
167 | } |
168 | |
169 | // Les approx sont a priori OK. |
170 | |
171 | const AppParCurves_MultiBSpCurve& MBSp = |
172 | Fit.SplineValue(); |
173 | Standard_Integer NbPoles = MBSp.NbPoles(); |
174 | TColgp_Array1OfPnt Poles (1 , NbPoles); |
175 | TColgp_Array1OfPnt2d Poles2d1(1 , NbPoles); |
176 | TColgp_Array1OfPnt2d Poles2d2(1 , NbPoles); |
177 | |
178 | MBSp.Curve(1, Poles); |
179 | MBSp.Curve(2, Poles2d1); |
180 | MBSp.Curve(3, Poles2d2); |
181 | |
182 | const TColStd_Array1OfReal& Knots = MBSp.Knots(); |
183 | const TColStd_Array1OfInteger& Mults = MBSp.Multiplicities(); |
184 | Standard_Integer Degree = MBSp.Degree(); |
185 | |
186 | myCurve = new Geom_BSplineCurve (Poles, Knots,Mults,Degree); |
187 | myPCurve1 = new Geom2d_BSplineCurve(Poles2d1,Knots,Mults,Degree); |
188 | myPCurve2 = new Geom2d_BSplineCurve(Poles2d2,Knots,Mults,Degree); |
189 | |
190 | myIsDone = Standard_True; |
191 | } |
192 | |
193 | |
194 | //======================================================================= |
195 | //function : IsDone |
196 | //purpose : |
197 | //======================================================================= |
198 | |
199 | Standard_Boolean BRepFill_ApproxSeewing::IsDone() const |
200 | { |
201 | return myIsDone; |
202 | } |
203 | |
204 | |
205 | //======================================================================= |
857ffd5e |
206 | //function : Handle(Geom_Curve)& |
7fd59977 |
207 | //purpose : |
208 | //======================================================================= |
209 | |
210 | const Handle(Geom_Curve)& BRepFill_ApproxSeewing::Curve() const |
211 | { |
212 | StdFail_NotDone_Raise_if( !myIsDone, |
213 | "BRepFill_ApproxSeewing::Curve"); |
214 | return myCurve; |
215 | } |
216 | |
217 | |
218 | //======================================================================= |
857ffd5e |
219 | //function : Handle(Geom2d_Curve)& |
7fd59977 |
220 | //purpose : |
221 | //======================================================================= |
222 | |
223 | const Handle(Geom2d_Curve)& BRepFill_ApproxSeewing::CurveOnF1() const |
224 | { |
225 | StdFail_NotDone_Raise_if( !myIsDone, |
226 | "BRepFill_ApproxSeewing::CurveOnF1"); |
227 | return myPCurve1; |
228 | } |
229 | |
230 | |
231 | //======================================================================= |
857ffd5e |
232 | //function : Handle(Geom2d_Curve)& |
7fd59977 |
233 | //purpose : |
234 | //======================================================================= |
235 | |
236 | const Handle(Geom2d_Curve)& BRepFill_ApproxSeewing::CurveOnF2() const |
237 | { |
238 | StdFail_NotDone_Raise_if( !myIsDone, |
239 | "BRepFill_ApproxSeewing::CurveOnF2"); |
240 | return myPCurve2; |
241 | } |
242 | |
243 | |