1 // Created on: 1999-05-13
2 // Created by: data exchange team
3 // Copyright (c) 1999-1999 Matra Datavision
4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
6 // This file is part of Open CASCADE Technology software library.
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
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.
14 // Alternatively, this file may be used under the terms of Open CASCADE
15 // commercial license or contractual agreement.
18 #include <Geom_BezierCurve.hxx>
19 #include <Geom_BSplineCurve.hxx>
20 #include <Geom_Circle.hxx>
21 #include <Geom_Conic.hxx>
22 #include <Geom_Curve.hxx>
23 #include <Geom_Line.hxx>
24 #include <Geom_TrimmedCurve.hxx>
25 #include <GeomConvert.hxx>
26 #include <GeomConvert_ApproxCurve.hxx>
27 #include <GeomConvert_BSplineCurveToBezierCurve.hxx>
28 #include <Precision.hxx>
29 #include <ShapeExtend.hxx>
30 #include <ShapeUpgrade_ConvertCurve3dToBezier.hxx>
31 #include <Standard_Type.hxx>
32 #include <TColGeom_HArray1OfCurve.hxx>
33 #include <TColgp_Array1OfPnt.hxx>
34 #include <TColStd_Array1OfReal.hxx>
36 IMPLEMENT_STANDARD_RTTIEXT(ShapeUpgrade_ConvertCurve3dToBezier,ShapeUpgrade_SplitCurve3d)
38 ShapeUpgrade_ConvertCurve3dToBezier::ShapeUpgrade_ConvertCurve3dToBezier()
40 mySegments = new TColGeom_HSequenceOfCurve;
41 mySplitParams = new TColStd_HSequenceOfReal;
42 myLineMode = Standard_True;
43 myCircleMode = Standard_True;
44 myConicMode = Standard_True;
47 //=======================================================================
50 //=======================================================================
52 void ShapeUpgrade_ConvertCurve3dToBezier::Compute()
55 mySplitParams->Clear();
56 Standard_Real precision = Precision::PConfusion();
57 Standard_Real First = mySplitValues->Value(1);
58 Standard_Real Last = mySplitValues->Value(mySplitValues->Length());
59 if(myCurve->IsKind(STANDARD_TYPE(Geom_TrimmedCurve))) {
60 Handle(Geom_TrimmedCurve) tmp = Handle(Geom_TrimmedCurve)::DownCast (myCurve);
61 Handle(Geom_Curve) BasCurve = tmp->BasisCurve();
62 ShapeUpgrade_ConvertCurve3dToBezier converter;
63 converter.Init(BasCurve,First,Last);
64 converter.SetSplitValues(mySplitValues);
66 mySplitValues->ChangeSequence() = converter.SplitValues()->Sequence();
67 myNbCurves = mySplitValues->Length()-1;
68 myStatus |= converter.myStatus;
69 mySegments->ChangeSequence() = converter.Segments()->Sequence();
70 mySplitParams->ChangeSequence() = converter.SplitParams()->Sequence();
73 else if(myCurve->IsKind(STANDARD_TYPE(Geom_BezierCurve))) {
74 Handle(Geom_BezierCurve) bezier = Handle(Geom_BezierCurve)::DownCast (myCurve);
75 myNbCurves = mySplitValues->Length()-1;
76 mySplitParams->Append(First);
77 mySplitParams->Append(Last);
78 if(First < precision && Last > 1 - precision) {
79 mySegments->Append(bezier);
80 myStatus = ShapeExtend::EncodeStatus (ShapeExtend_OK);
82 Handle(Geom_BezierCurve) besNew = Handle(Geom_BezierCurve)::DownCast(bezier->Copy());
83 besNew->Segment(First,Last);
84 mySegments->Append(besNew);
85 myStatus = ShapeExtend::EncodeStatus (ShapeExtend_DONE2);
89 else if(myCurve->IsKind(STANDARD_TYPE(Geom_Line))) {
90 Handle(Geom_Line) aLine = Handle(Geom_Line)::DownCast(myCurve);
91 myNbCurves = mySplitValues->Length()-1;
92 mySplitParams->Append(First);
93 mySplitParams->Append(Last);
95 mySegments->Append(aLine);
96 myStatus = ShapeExtend::EncodeStatus (ShapeExtend_OK);
99 TColgp_Array1OfPnt poles(1,2);
100 poles(1) = aLine->Value(First);
101 poles(2) = aLine->Value(Last);
102 Handle(Geom_BezierCurve) bezier = new Geom_BezierCurve(poles);
103 mySegments->Append(bezier);
104 myStatus = ShapeExtend::EncodeStatus (ShapeExtend_DONE1);
107 else if((myCurve->IsKind(STANDARD_TYPE(Geom_Conic))&&!myConicMode)||
108 (myCurve->IsKind(STANDARD_TYPE(Geom_Circle))&&!myCircleMode)) {
109 myNbCurves = mySplitValues->Length()-1;
110 mySplitParams->Append(First);
111 mySplitParams->Append(Last);
112 mySegments->Append(myCurve);
113 myStatus = ShapeExtend::EncodeStatus (ShapeExtend_OK);
117 Handle(Geom_BSplineCurve) aBSpline;
118 Standard_Real Shift = 0.;
119 if(myCurve->IsKind(STANDARD_TYPE(Geom_Conic))) {
120 Handle(Geom_Curve) tcurve = new Geom_TrimmedCurve(myCurve,First,Last); //protection agains parabols ets
121 GeomConvert_ApproxCurve approx (tcurve, Precision::Approximation(),
122 GeomAbs_C1, 100, 6 );
123 if ( approx.HasResult() )
124 aBSpline = approx.Curve();
126 Handle(Geom_TrimmedCurve) t3d = new Geom_TrimmedCurve(myCurve,First,Last);
127 aBSpline = GeomConvert::CurveToBSplineCurve(t3d,Convert_QuasiAngular);
129 Shift = First - aBSpline->FirstParameter();
130 First = aBSpline->FirstParameter();
131 Last = aBSpline->LastParameter();
133 else if(!myCurve->IsKind(STANDARD_TYPE(Geom_BSplineCurve))) {
134 aBSpline = GeomConvert::CurveToBSplineCurve(myCurve,Convert_QuasiAngular);
137 aBSpline = Handle(Geom_BSplineCurve)::DownCast(myCurve);
139 Standard_Real bf = aBSpline->FirstParameter();
140 Standard_Real bl = aBSpline->LastParameter();
141 if(Abs(First-bf) < precision)
143 if(Abs(Last-bl) < precision)
147 cout <<"Warning: The range of the edge exceeds the curve domain" <<endl;
150 mySplitValues->SetValue(1,First);
154 cout <<"Warning: The range of the edge exceeds the curve domain" <<endl;
157 mySplitValues->SetValue(mySplitValues->Length(),Last);
160 GeomConvert_BSplineCurveToBezierCurve tool(aBSpline,First,Last,precision);
161 Standard_Integer nbArcs = tool.NbArcs();
162 TColStd_Array1OfReal knots(1,nbArcs+1);
164 mySplitParams->Append(First+Shift);
165 Standard_Integer j; // svv Jan 10 2000 : porting on DEC
166 for(j = 1; j <=nbArcs; j++) {
167 Standard_Real nextKnot = knots(j+1)+Shift;
168 if(nextKnot - mySplitParams->Value(mySplitParams->Length()) > precision) {
169 mySegments->Append(tool.Arc(j));
170 mySplitParams->Append(knots(j+1)+Shift);
174 First = mySplitValues->Value(1);
175 for(j = 2; j <= mySplitValues->Length(); j++) {
176 Last = mySplitValues->Value(j);
177 for(Standard_Integer i = 2; i <= nbArcs+1; i++) {
178 Standard_Real valknot = knots(i)+Shift;
179 if(valknot <= First + precision) continue;
180 if(valknot >= Last - precision) break;
181 mySplitValues->InsertBefore(j++,valknot);
185 myNbCurves = mySplitValues->Length()-1;
187 myStatus = ShapeExtend::EncodeStatus (ShapeExtend_DONE1);
190 //=======================================================================
193 //=======================================================================
195 void ShapeUpgrade_ConvertCurve3dToBezier::Build(const Standard_Boolean /*Segment*/)
197 Standard_Real prec = Precision::PConfusion();
198 Standard_Integer nb = mySplitValues->Length();
199 myResultingCurves = new TColGeom_HArray1OfCurve (1,nb-1);
200 Standard_Real prevPar = 0.;
201 Standard_Integer j=2;
202 for(Standard_Integer i = 2; i <= nb; i++) {
203 Standard_Real par = mySplitValues->Value(i);
204 for(; j<= mySplitParams->Length(); j++)
205 if(mySplitParams->Value(j)+prec > par)
210 Handle(Geom_Curve) crv = Handle(Geom_Curve)::DownCast(mySegments->Value(j-1)->Copy());
211 if(crv->IsKind(STANDARD_TYPE(Geom_BezierCurve))) {
212 Handle(Geom_BezierCurve) bes = Handle(Geom_BezierCurve)::DownCast(crv);
213 Standard_Real uFact = mySplitParams->Value(j) - mySplitParams->Value(j-1);
214 Standard_Real pp = mySplitValues->Value(i-1);
215 Standard_Real length = (par - pp)/uFact;
216 bes->Segment(prevPar, prevPar+length);
218 myResultingCurves->SetValue(i-1,bes);
221 myResultingCurves->SetValue(i-1,crv);
225 //=======================================================================
226 //function : Segments
228 //=======================================================================
230 Handle(TColGeom_HSequenceOfCurve) ShapeUpgrade_ConvertCurve3dToBezier::Segments() const
235 //=======================================================================
236 //function : SplitParams
238 //=======================================================================
240 Handle(TColStd_HSequenceOfReal) ShapeUpgrade_ConvertCurve3dToBezier::SplitParams() const
242 return mySplitParams;