b311480e |
1 | // Created on: 1994-09-21 |
2 | // Created by: Bruno DUMORTIER |
3 | // Copyright (c) 1994-1999 Matra Datavision |
4 | // Copyright (c) 1999-2012 OPEN CASCADE SAS |
5 | // |
6 | // The content of this file is subject to the Open CASCADE Technology Public |
7 | // License Version 6.5 (the "License"). You may not use the content of this file |
8 | // except in compliance with the License. Please obtain a copy of the License |
9 | // at http://www.opencascade.org and read it completely before using this file. |
10 | // |
11 | // The Initial Developer of the Original Code is Open CASCADE S.A.S., having its |
12 | // main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France. |
13 | // |
14 | // The Original Code and all software distributed under the License is |
15 | // distributed on an "AS IS" basis, without warranty of any kind, and the |
16 | // Initial Developer hereby disclaims all such warranties, including without |
17 | // limitation, any warranties of merchantability, fitness for a particular |
18 | // purpose or non-infringement. Please see the License for the specific terms |
19 | // and conditions governing the rights and limitations under the License. |
20 | |
7fd59977 |
21 | |
22 | // Modified by skv - Wed Aug 11 17:26:03 2004 OCC6272 |
23 | |
24 | #include <stdio.h> |
25 | |
26 | #include <GeomProjLib.ixx> |
27 | |
28 | #include <ProjLib_ProjectedCurve.hxx> |
29 | #include <ProjLib_CompProjectedCurve.hxx> |
30 | #include <ProjLib_ProjectOnPlane.hxx> |
31 | #include <ProjLib_ProjectOnSurface.hxx> |
32 | |
33 | #include <Approx_CurveOnSurface.hxx> |
34 | |
35 | #include <GeomAdaptor_Curve.hxx> |
36 | #include <GeomAdaptor_Surface.hxx> |
37 | #include <GeomAdaptor_HCurve.hxx> |
38 | #include <GeomAdaptor_HSurface.hxx> |
39 | |
40 | #include <Geom_Line.hxx> |
41 | #include <Geom_Circle.hxx> |
42 | #include <Geom_Ellipse.hxx> |
43 | #include <Geom_Parabola.hxx> |
44 | #include <Geom_Hyperbola.hxx> |
45 | #include <Geom_BezierCurve.hxx> |
46 | #include <Geom_BSplineCurve.hxx> |
47 | #include <Geom_TrimmedCurve.hxx> |
48 | #include <Geom2d_Line.hxx> |
49 | #include <Geom2d_Circle.hxx> |
50 | #include <Geom2d_Ellipse.hxx> |
51 | #include <Geom2d_Parabola.hxx> |
52 | #include <Geom2d_Hyperbola.hxx> |
53 | #include <Geom2d_BezierCurve.hxx> |
54 | #include <Geom2d_BSplineCurve.hxx> |
55 | #include <Geom2d_TrimmedCurve.hxx> |
56 | |
57 | #include <gp_Pln.hxx> |
58 | |
59 | #include <Precision.hxx> |
60 | #include <TColStd_Array1OfReal.hxx> |
61 | #include <TColStd_Array1OfInteger.hxx> |
62 | #include <TColgp_Array1OfPnt.hxx> |
63 | #include <TColgp_Array1OfPnt2d.hxx> |
64 | |
65 | #ifdef DRAW |
66 | #include <DrawTrSurf.hxx> |
67 | #endif |
68 | #ifdef DEB |
69 | static Standard_Boolean Affich = Standard_False; |
70 | static Standard_Integer NBPROJ = 1; |
71 | #endif |
72 | |
73 | |
74 | //======================================================================= |
75 | //function : Curve2d |
76 | //purpose : |
77 | //======================================================================= |
78 | |
79 | Handle(Geom2d_Curve) GeomProjLib::Curve2d(const Handle(Geom_Curve)& C, |
80 | const Standard_Real First, |
81 | const Standard_Real Last, |
82 | const Handle(Geom_Surface)& S, |
83 | const Standard_Real UDeb, |
84 | const Standard_Real UFin, |
85 | const Standard_Real VDeb, |
86 | const Standard_Real VFin, |
87 | Standard_Real& Tolerance) |
88 | { |
89 | #ifdef DRAW |
90 | if ( Affich) { |
91 | //POP pour NT |
92 | // char name[100]; |
93 | char* name = new char[100]; |
94 | sprintf(name,"PROJCURV_%d",NBPROJ); |
95 | DrawTrSurf::Set(name,C); |
96 | sprintf(name,"PROJSURF_%d",NBPROJ); |
97 | DrawTrSurf::Set(name,S); |
98 | NBPROJ++; |
99 | } |
100 | #endif |
101 | |
102 | Tolerance = Max(Precision::PConfusion(),Tolerance); |
103 | |
104 | GeomAdaptor_Curve AC(C,First,Last); |
105 | GeomAdaptor_Surface AS(S, |
106 | UDeb, |
107 | UFin, |
108 | VDeb, |
109 | VFin); |
110 | |
111 | Handle(GeomAdaptor_HSurface) HS = new GeomAdaptor_HSurface(AS); |
112 | Handle(GeomAdaptor_HCurve) HC = new GeomAdaptor_HCurve(AC); |
113 | |
114 | ProjLib_ProjectedCurve Proj(HS,HC,Tolerance); |
115 | |
116 | Handle(Geom2d_Curve) G2dC; |
117 | |
118 | switch ( Proj.GetType()) { |
119 | |
120 | case GeomAbs_Line: |
121 | G2dC = new Geom2d_Line(Proj.Line()); |
122 | break; |
123 | |
124 | case GeomAbs_Circle: |
125 | G2dC = new Geom2d_Circle(Proj.Circle()); |
126 | break; |
127 | |
128 | case GeomAbs_Ellipse: |
129 | G2dC = new Geom2d_Ellipse(Proj.Ellipse()); |
130 | break; |
131 | |
132 | case GeomAbs_Parabola: |
133 | G2dC = new Geom2d_Parabola(Proj.Parabola()); |
134 | break; |
135 | |
136 | case GeomAbs_Hyperbola: |
137 | G2dC = new Geom2d_Hyperbola(Proj.Hyperbola()); |
138 | break; |
139 | |
140 | case GeomAbs_BezierCurve: |
141 | G2dC = Proj.Bezier(); |
142 | break; |
143 | |
144 | case GeomAbs_BSplineCurve: |
145 | G2dC = Proj.BSpline(); |
146 | break; |
147 | |
148 | default: |
149 | return G2dC; |
150 | |
151 | } |
152 | |
153 | if(G2dC.IsNull()) { |
154 | Tolerance = Proj.GetTolerance(); |
155 | return G2dC; |
156 | } |
157 | |
158 | if ( C->IsKind(STANDARD_TYPE(Geom_TrimmedCurve)) ) { |
159 | Handle(Geom_TrimmedCurve) CTrim = Handle(Geom_TrimmedCurve)::DownCast(C); |
160 | Standard_Real U1 = CTrim->FirstParameter(); |
161 | Standard_Real U2 = CTrim->LastParameter(); |
162 | G2dC = new Geom2d_TrimmedCurve( G2dC, U1, U2); |
163 | } |
164 | |
165 | #ifdef DRAW |
166 | if ( Affich) { |
167 | static Standard_CString aprojcurv = "projcurv" ; |
168 | DrawTrSurf::Set(aprojcurv,G2dC); |
169 | } |
170 | #endif |
171 | Tolerance = Proj.GetTolerance(); |
172 | return G2dC; |
173 | } |
174 | |
175 | //======================================================================= |
176 | //function : Curve2d |
177 | //purpose : |
178 | //======================================================================= |
179 | |
180 | Handle(Geom2d_Curve) GeomProjLib::Curve2d(const Handle(Geom_Curve)& C, |
181 | const Standard_Real First, |
182 | const Standard_Real Last, |
183 | const Handle(Geom_Surface)& S, |
184 | Standard_Real& Tolerance) |
185 | { |
186 | Standard_Real UFirst, |
187 | ULast, |
188 | VFirst, |
189 | VLast ; |
190 | |
191 | S->Bounds(UFirst, |
192 | ULast, |
193 | VFirst, |
194 | VLast) ; |
195 | return Curve2d(C, |
196 | First, |
197 | Last, |
198 | S, |
199 | UFirst, |
200 | ULast, |
201 | VFirst, |
202 | VLast, |
203 | Tolerance) ; |
204 | } |
205 | // Modified by skv - Wed Aug 11 17:26:03 2004 OCC6272 Begin |
206 | // Add not implemented method. |
207 | //======================================================================= |
208 | //function : Curve2d |
209 | //purpose : |
210 | //======================================================================= |
211 | |
212 | Handle(Geom2d_Curve) GeomProjLib::Curve2d( const Handle(Geom_Curve)& C, |
213 | const Handle(Geom_Surface)& S, |
214 | const Standard_Real UDeb, |
215 | const Standard_Real UFin, |
216 | const Standard_Real VDeb, |
217 | const Standard_Real VFin) |
218 | { |
219 | Standard_Real First = C->FirstParameter(); |
220 | Standard_Real Last = C->LastParameter(); |
221 | Standard_Real Tol = Precision::PConfusion(); |
222 | return GeomProjLib::Curve2d(C,First,Last,S,UDeb,UFin,VDeb,VFin,Tol); |
223 | } |
224 | // Modified by skv - Wed Aug 11 17:26:03 2004 OCC6272 End |
225 | |
226 | //======================================================================= |
227 | //function : Curve2d |
228 | //purpose : |
229 | //======================================================================= |
230 | |
231 | Handle(Geom2d_Curve) GeomProjLib::Curve2d( const Handle(Geom_Curve)& C, |
232 | const Handle(Geom_Surface)& S, |
233 | const Standard_Real UDeb, |
234 | const Standard_Real UFin, |
235 | const Standard_Real VDeb, |
236 | const Standard_Real VFin, |
237 | Standard_Real& Tolerance) |
238 | { |
239 | Standard_Real First = C->FirstParameter(); |
240 | Standard_Real Last = C->LastParameter(); |
241 | return GeomProjLib::Curve2d(C,First,Last,S,UDeb,UFin,VDeb,VFin,Tolerance); |
242 | } |
243 | |
244 | |
245 | //======================================================================= |
246 | //function : Curve2d |
247 | //purpose : |
248 | //======================================================================= |
249 | |
250 | Handle(Geom2d_Curve) GeomProjLib::Curve2d( const Handle(Geom_Curve)& C, |
251 | const Handle(Geom_Surface)& S) |
252 | { |
253 | Standard_Real First = C->FirstParameter(); |
254 | Standard_Real Last = C->LastParameter(); |
255 | Standard_Real Tol = Precision::PConfusion(); |
256 | return GeomProjLib::Curve2d(C,First,Last,S,Tol); |
257 | } |
258 | |
259 | |
260 | //======================================================================= |
261 | //function : Curve2d |
262 | //purpose : |
263 | //======================================================================= |
264 | |
265 | Handle(Geom2d_Curve) GeomProjLib::Curve2d( const Handle(Geom_Curve)& C, |
266 | const Standard_Real First, |
267 | const Standard_Real Last, |
268 | const Handle(Geom_Surface)& S) |
269 | { |
270 | Standard_Real Tol = Precision::PConfusion(); |
271 | return GeomProjLib::Curve2d(C,First,Last,S,Tol); |
272 | } |
273 | |
274 | |
275 | //======================================================================= |
276 | //function : Project |
277 | //purpose : |
278 | //======================================================================= |
279 | |
280 | Handle(Geom_Curve) GeomProjLib::Project( const Handle(Geom_Curve)& C, |
281 | const Handle(Geom_Surface)& S) |
282 | { |
283 | GeomAdaptor_Curve AC(C); |
284 | GeomAdaptor_Surface AS(S); |
285 | |
286 | Handle(Geom_Curve) GC; |
287 | |
288 | if ( AS.GetType() == GeomAbs_Plane) { |
289 | ProjLib_ProjectOnPlane Proj( AS.Plane().Position()); |
290 | Handle(GeomAdaptor_HCurve) HC = new GeomAdaptor_HCurve(AC); |
291 | Proj.Load(HC,Precision::PApproximation()); |
292 | |
293 | switch ( Proj.GetType()) { |
294 | case GeomAbs_Line: |
295 | GC = new Geom_Line(Proj.Line()); |
296 | break; |
297 | |
298 | case GeomAbs_Circle: |
299 | GC = new Geom_Circle(Proj.Circle()); |
300 | break; |
301 | |
302 | case GeomAbs_Ellipse: |
303 | GC = new Geom_Ellipse(Proj.Ellipse()); |
304 | break; |
305 | |
306 | case GeomAbs_Parabola: |
307 | GC = new Geom_Parabola(Proj.Parabola()); |
308 | break; |
309 | |
310 | case GeomAbs_Hyperbola: |
311 | GC = new Geom_Hyperbola(Proj.Hyperbola()); |
312 | break; |
313 | |
314 | case GeomAbs_BezierCurve: |
315 | GC = Proj.Bezier(); |
316 | break; |
317 | |
318 | case GeomAbs_BSplineCurve: |
319 | GC = Proj.BSpline(); |
320 | break; |
321 | |
322 | default: |
323 | return GC; |
324 | |
325 | } |
326 | |
327 | if ( C->IsKind(STANDARD_TYPE(Geom_TrimmedCurve))) { |
328 | Handle(Geom_TrimmedCurve) CTrim = Handle(Geom_TrimmedCurve)::DownCast(C); |
329 | Standard_Real U1 = CTrim->FirstParameter(); |
330 | Standard_Real U2 = CTrim->LastParameter(); |
331 | GC = new Geom_TrimmedCurve( GC, U1, U2); |
332 | } |
333 | |
334 | } |
335 | else { |
336 | Handle(GeomAdaptor_HSurface) HS = new GeomAdaptor_HSurface(AS); |
337 | Handle(GeomAdaptor_HCurve) HC = new GeomAdaptor_HCurve(AC); |
338 | // Standard_Real Tol = Precision::Approximation(); |
339 | // Standard_Real TolU = Precision::PApproximation(); |
340 | // Standard_Real TolV = Precision::PApproximation(); |
341 | Standard_Real Tol = 0.0001; |
342 | Standard_Real TolU = Pow(Tol, 2./3); |
343 | Standard_Real TolV = Pow(Tol, 2./3); |
344 | ProjLib_CompProjectedCurve Proj(HS,HC,TolU,TolV,-1.); |
345 | |
7fd59977 |
346 | Standard_Real f,l; |
347 | Proj.Bounds(1,f,l); |
348 | Handle(Adaptor2d_HCurve2d) HC2d = Proj.Trim(f,l,TolU); |
349 | Approx_CurveOnSurface Approx(HC2d, HS, f, l, Tol, |
350 | GeomAbs_C2,14,16,Standard_True); |
351 | |
352 | // ici, on a toujours un type BSpline. |
353 | if (Approx.IsDone() && Approx.HasResult()) |
354 | GC = Approx.Curve3d(); |
355 | } |
356 | |
357 | return GC; |
358 | } |
359 | |
360 | //======================================================================= |
361 | //function : ProjectOnPlane |
362 | //purpose : |
363 | //======================================================================= |
364 | |
365 | Handle(Geom_Curve) GeomProjLib::ProjectOnPlane |
366 | (const Handle(Geom_Curve)& Curve, |
367 | const Handle(Geom_Plane)& Plane, |
368 | const gp_Dir& Dir, |
369 | const Standard_Boolean KeepParametrization) |
370 | { |
371 | GeomAdaptor_Curve AC(Curve); |
372 | Handle(GeomAdaptor_HCurve) HC = new GeomAdaptor_HCurve(AC); |
373 | |
374 | ProjLib_ProjectOnPlane Proj(Plane->Position(), Dir); |
375 | Proj.Load(HC,Precision::Approximation(), KeepParametrization); |
376 | |
377 | Handle(Geom_Curve) GC; |
378 | |
379 | switch ( Proj.GetType()) { |
380 | case GeomAbs_Line: |
381 | GC = new Geom_Line(Proj.Line()); |
382 | break; |
383 | |
384 | case GeomAbs_Circle: |
385 | GC = new Geom_Circle(Proj.Circle()); |
386 | break; |
387 | |
388 | case GeomAbs_Ellipse: |
389 | GC = new Geom_Ellipse(Proj.Ellipse()); |
390 | break; |
391 | |
392 | case GeomAbs_Parabola: |
393 | GC = new Geom_Parabola(Proj.Parabola()); |
394 | break; |
395 | |
396 | case GeomAbs_Hyperbola: |
397 | GC = new Geom_Hyperbola(Proj.Hyperbola()); |
398 | break; |
399 | |
400 | case GeomAbs_BezierCurve: |
401 | GC = Proj.Bezier(); |
402 | break; |
403 | |
404 | case GeomAbs_BSplineCurve: |
405 | GC = Proj.BSpline(); |
406 | break; |
407 | default: |
408 | return GC; |
409 | |
410 | } |
411 | |
412 | if ( Curve->IsKind(STANDARD_TYPE(Geom_TrimmedCurve)) ) { |
413 | Handle(Geom_TrimmedCurve) CTrim |
414 | = Handle(Geom_TrimmedCurve)::DownCast(Curve); |
415 | GC = new Geom_TrimmedCurve( GC, Proj.FirstParameter(), |
416 | Proj.LastParameter()); |
417 | } |
418 | |
419 | return GC; |
420 | |
421 | } |
422 | |
423 | |