9dfed5f9e0f5b7f2702966844357b22329ad380e
[occt.git] / src / PrsDim / PrsDim_OffsetDimension.cxx
1 // Created on: 1997-03-04
2 // Created by: Jean-Pierre COMBE
3 // Copyright (c) 1997-1999 Matra Datavision
4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
5 //
6 // This file is part of Open CASCADE Technology software library.
7 //
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.
13 //
14 // Alternatively, this file may be used under the terms of Open CASCADE
15 // commercial license or contractual agreement.
16
17 #include <PrsDim_OffsetDimension.hxx>
18
19 #include <PrsDim.hxx>
20 #include <BRep_Tool.hxx>
21 #include <BRepAdaptor_Surface.hxx>
22 #include <BRepBuilderAPI_Transform.hxx>
23 #include <DsgPrs_OffsetPresentation.hxx>
24 #include <ElCLib.hxx>
25 #include <ElSLib.hxx>
26 #include <gce_MakeLin.hxx>
27 #include <GeomAbs_SurfaceType.hxx>
28 #include <gp_Ax1.hxx>
29 #include <gp_Ax2.hxx>
30 #include <gp_Ax3.hxx>
31 #include <gp_Lin.hxx>
32 #include <gp_Pln.hxx>
33 #include <gp_Pnt.hxx>
34 #include <gp_Trsf.hxx>
35 #include <gp_Vec.hxx>
36 #include <Graphic3d_Structure.hxx>
37 #include <Precision.hxx>
38 #include <Prs3d_ArrowAspect.hxx>
39 #include <Prs3d_DimensionAspect.hxx>
40 #include <Prs3d_Drawer.hxx>
41 #include <Prs3d_LineAspect.hxx>
42 #include <Prs3d_Presentation.hxx>
43 #include <Select3D_SensitiveBox.hxx>
44 #include <Select3D_SensitiveSegment.hxx>
45 #include <SelectMgr_EntityOwner.hxx>
46 #include <SelectMgr_Selection.hxx>
47 #include <Standard_NotImplemented.hxx>
48 #include <StdPrs_WFShape.hxx>
49 #include <TCollection_ExtendedString.hxx>
50 #include <TColStd_Array2OfReal.hxx>
51 #include <TopExp_Explorer.hxx>
52 #include <TopoDS.hxx>
53 #include <TopoDS_Shape.hxx>
54 #include <TopoDS_Vertex.hxx>
55
56 IMPLEMENT_STANDARD_RTTIEXT(PrsDim_OffsetDimension, PrsDim_Relation)
57
58 //=======================================================================
59 //function : PrsDim_OffsetDimension
60 //purpose  : 
61 //=======================================================================
62 PrsDim_OffsetDimension::PrsDim_OffsetDimension(const TopoDS_Shape& FistShape, 
63                                          const TopoDS_Shape& SecondShape,
64                                          const Standard_Real aVal,
65                                          const TCollection_ExtendedString& aText)
66 :PrsDim_Relation(),
67 myFAttach(0.,0.,0.),
68 mySAttach(0.,0.,0.)
69 {
70   myFShape = FistShape;
71   mySShape = SecondShape;
72   mySymbolPrs = DsgPrs_AS_BOTHAR;
73   myVal = aVal;
74   myText = aText;
75   //myArrowSize = fabs (myVal/5.);
76   myArrowSize = fabs (myVal/10.0);
77   if (myArrowSize > 30.) myArrowSize = 30.;
78   if (myArrowSize < 15.) myArrowSize = 15.;
79   //std::cout<<"PrsDim_OffsetDimension::PrsDim_OffsetDimension " <<  myArrowSize << " myArrowSize"<<std::endl;
80 }
81
82 //=======================================================================
83 //function : Compute
84 //purpose  : 
85 //=======================================================================
86 void PrsDim_OffsetDimension::Compute(const Handle(PrsMgr_PresentationManager3d)&,
87                                   const Handle(Prs3d_Presentation)& aprs,
88                                   const Standard_Integer)
89 {
90   gp_Trsf aInvertTrsf = myRelativePos;
91   //myArrowSize = fabs (myVal/5.);
92   myArrowSize = fabs (myVal/10.0);
93   if (myArrowSize > 30.) myArrowSize = 30.;
94   if (myArrowSize < 15.) myArrowSize = 15.;
95   //std::cout<<"PrsDim_OffsetDimension::PrsDim_OffsetDimension " <<  myArrowSize << " myArrowSize"<<std::endl;
96   
97   BRepAdaptor_Surface surf1(TopoDS::Face(myFShape));
98   BRepAdaptor_Surface surf2(TopoDS::Face(mySShape));
99   
100   if (surf1.GetType() == GeomAbs_Cylinder || 
101       surf1.GetType() == GeomAbs_Cone     || 
102       surf1.GetType() == GeomAbs_Torus) {
103     if (surf2.GetType() == GeomAbs_Cylinder ||
104         surf2.GetType() == GeomAbs_Cone     ||
105         surf2.GetType() == GeomAbs_Torus) {
106       ComputeTwoAxesOffset(aprs, aInvertTrsf);
107     } else {
108       ComputeAxeFaceOffset(aprs, aInvertTrsf);
109     }
110   }
111   else {
112     //myDirAttach : oriente de myFShape vers mySShape
113     gp_Pln aPln = surf1.Plane();
114     gp_Pnt aPnt = aPln.Location();
115     
116     gp_Pln bPln = surf2.Plane();
117     
118     Standard_Real uPnt, vPnt;
119     ElSLib::Parameters (bPln , aPnt , uPnt, vPnt);
120     gp_Pnt bPnt = ElSLib::Value (uPnt, vPnt, bPln);
121     if (aPnt.IsEqual(bPnt,Precision::Confusion())) {
122       gp_Ax1 aAx1 = aPln.Axis();
123       myDirAttach = aAx1.Direction();
124     } else {
125       gp_Vec aVec (aPnt,bPnt);
126       myDirAttach.SetCoord(aVec.X(),aVec.Y(),aVec.Z());
127     }
128     ComputeTwoFacesOffset(aprs, aInvertTrsf);
129   }  
130 }
131
132 //=======================================================================
133 //function : ComputeSelection
134 //purpose  : 
135 //=======================================================================
136 void PrsDim_OffsetDimension::ComputeSelection(const Handle(SelectMgr_Selection)& aSel, 
137                                            const Standard_Integer)
138 {
139   //myArrowSize = fabs (myVal/5.);
140   myArrowSize = fabs (myVal/10.0);
141   if (myArrowSize > 30.) myArrowSize = 30.;
142   if (myArrowSize < 15.) myArrowSize = 15.;
143   //std::cout<<"PrsDim_OffsetDimension::PrsDim_OffsetDimension " <<  myArrowSize << " myArrowSize"<<std::endl;
144   gp_Pnt myTFAttach = myFAttach.Transformed (myRelativePos);
145   gp_Pnt myTSAttach = mySAttach.Transformed (myRelativePos);
146   gp_Dir myTDirAttach = myDirAttach.Transformed (myRelativePos);
147   gp_Dir myTDirAttach2 = myDirAttach2.Transformed (myRelativePos);
148   gp_Pnt Tcurpos = myPosition.Transformed (myRelativePos);
149
150   gp_Lin L1 (myTFAttach,myTDirAttach);
151   gp_Lin L2 (myTSAttach,myTDirAttach2);
152   gp_Pnt Proj1 = ElCLib::Value(ElCLib::Parameter(L1,Tcurpos),L1);
153   gp_Pnt Proj2 = ElCLib::Value(ElCLib::Parameter(L2,Tcurpos),L2);
154   gp_Lin L3;
155
156   Handle(SelectMgr_EntityOwner) own = new SelectMgr_EntityOwner(this,7);
157
158   if (!Proj1.IsEqual(Proj2,Precision::Confusion())) {
159     L3 = gce_MakeLin(Proj1,Proj2);
160   }
161   else {    // cas ou la dimension est nulle
162     if (!Proj1.IsEqual(Tcurpos,Precision::Confusion())) {
163       gp_Vec v3 (Proj1,Tcurpos);
164       gp_Dir d3 (v3);
165       L3 = gce_MakeLin(Proj1,d3);
166     } else {
167       L3 = gce_MakeLin(Proj1,myTDirAttach);
168     }
169
170   // Text
171     Standard_Real size(Min(myVal/100.+1.e-6,myArrowSize+1.e-6));
172     Handle( Select3D_SensitiveBox ) box = new Select3D_SensitiveBox( own,
173                                                                      Tcurpos.X(),
174                                                                      Tcurpos.Y(),
175                                                                      Tcurpos.Z(),
176                                                                      Tcurpos.X() + size,
177                                                                      Tcurpos.Y() + size,
178                                                                      Tcurpos.Z() + size);
179     aSel->Add(box);
180   }
181
182   Standard_Real parmin,parmax,parcur;
183   parmin = ElCLib::Parameter(L3,Proj1);
184   parmax = parmin;
185
186   parcur = ElCLib::Parameter(L3,Proj2);
187   parmin = Min(parmin,parcur);
188   parmax = Max(parmax,parcur);
189
190   parcur = ElCLib::Parameter(L3,Tcurpos);
191   parmin = Min(parmin,parcur);
192   parmax = Max(parmax,parcur);
193
194   gp_Pnt PointMin = ElCLib::Value(parmin,L3);
195   gp_Pnt PointMax = ElCLib::Value(parmax,L3);
196
197   Handle(Select3D_SensitiveSegment) seg;
198   if (!PointMin.IsEqual(PointMax,Precision::Confusion())) {
199     seg = new Select3D_SensitiveSegment(own,
200                                         PointMin,
201                                         PointMax);
202     aSel->Add(seg);
203   }
204   
205   if (!myTFAttach.IsEqual(Proj1,Precision::Confusion())) {
206     seg = new Select3D_SensitiveSegment(own,
207                                         myTFAttach,
208                                         Proj1);
209     aSel->Add(seg);
210   }
211   if (!myTSAttach.IsEqual(Proj2,Precision::Confusion())) {
212     seg = new Select3D_SensitiveSegment(own,
213                                         myTSAttach,
214                                         Proj2);
215     aSel->Add(seg);
216   }
217 }
218
219 //=======================================================================
220 //function : ComputeTwoAxesOffset
221 //purpose  : 
222 //=======================================================================
223 void PrsDim_OffsetDimension::ComputeTwoAxesOffset(const Handle(Prs3d_Presentation)& aprs, 
224                                                 const gp_Trsf& aTrsf)
225 {
226   BRepAdaptor_Surface surf1(TopoDS::Face(myFShape));
227   BRepAdaptor_Surface surf2(TopoDS::Face(mySShape));
228
229   gp_Ax1 Ax1Surf1, Ax1Surf2;
230
231   if (surf1.GetType() == GeomAbs_Cylinder) {
232     gp_Cylinder aCyl= surf1.Cylinder();
233     Ax1Surf1 = aCyl.Axis();    
234   } else if (surf1.GetType() == GeomAbs_Cone) {
235     gp_Cone aCone= surf1.Cone();
236     Ax1Surf1 = aCone.Axis();    
237   } else if (surf1.GetType() == GeomAbs_Torus) {
238     gp_Torus aTore= surf1.Torus();
239     Ax1Surf1 = aTore.Axis();    
240   }
241   Standard_Real FirstUParam = surf1.FirstUParameter();
242   Standard_Real FirstVParam = surf1.FirstVParameter();
243   Standard_Real LastVParam  = surf1.LastVParameter();
244   gp_Pnt P1First = surf1.Value(FirstUParam,FirstVParam);
245   gp_Pnt P1Last  = surf1.Value(FirstUParam,LastVParam);
246
247
248   if (surf2.GetType() == GeomAbs_Cylinder) {
249     gp_Cylinder aCyl= surf2.Cylinder();
250     Ax1Surf2 = aCyl.Axis();    
251   } else if (surf2.GetType() == GeomAbs_Cone) {
252     gp_Cone aCone= surf2.Cone();
253     Ax1Surf2 = aCone.Axis();    
254   } else if (surf2.GetType() == GeomAbs_Torus) {
255     gp_Torus aTore= surf2.Torus();
256     Ax1Surf2 = aTore.Axis();    
257   }
258   FirstUParam = surf2.FirstUParameter();
259   FirstVParam = surf2.FirstVParameter();
260   LastVParam  = surf2.LastVParameter();
261   gp_Pnt P2First = surf2.Value(FirstUParam,FirstVParam);
262   gp_Pnt P2Last  = surf2.Value(FirstUParam,LastVParam);
263
264   
265   
266   myFAttach = Ax1Surf1.Location();
267   mySAttach = Ax1Surf2.Location();
268   myDirAttach = Ax1Surf1.Direction();
269   myDirAttach2 = myDirAttach;
270   gp_Pnt curpos;
271   gp_Lin aProjLine  = gce_MakeLin(myFAttach,myDirAttach);
272
273   if (myAutomaticPosition) {
274     curpos.SetX ( (myFAttach.X() +  mySAttach.X()) /2. + 0.01);
275     curpos.SetY ( (myFAttach.Y() +  mySAttach.Y()) /2. + 0.01);
276     curpos.SetZ ( (myFAttach.Z() +  mySAttach.Z()) /2. + 0.01);
277     // + 0.01 pour eviter un raise de ComputeSelection...
278
279     myPosition = curpos;
280   } 
281   else {
282     curpos = myPosition;
283   }
284   
285   curpos = ElCLib::Value(ElCLib::Parameter(aProjLine,curpos),aProjLine);
286   // on projette pour la presentation
287
288   gp_Pnt P1FirstProj  = ElCLib::Value(ElCLib::Parameter(aProjLine,P1First),aProjLine);
289   gp_Pnt P1LastProj   = ElCLib::Value(ElCLib::Parameter(aProjLine,P1Last),aProjLine);
290   if (P1FirstProj.Distance(curpos) > P1LastProj.Distance(curpos))
291     myFAttach = P1FirstProj;
292   else
293     myFAttach = P1LastProj;
294   
295   gp_Pnt P2FirstProj  = ElCLib::Value(ElCLib::Parameter(aProjLine,P2First),aProjLine);
296   gp_Pnt P2LastProj   = ElCLib::Value(ElCLib::Parameter(aProjLine,P2Last),aProjLine);
297   if (P2FirstProj.Distance(curpos) > P2LastProj.Distance(curpos))
298     mySAttach = P2FirstProj;
299   else
300     mySAttach = P2LastProj;
301
302
303   Handle(Prs3d_DimensionAspect) la = myDrawer->DimensionAspect();
304   Handle(Prs3d_ArrowAspect) arr = la->ArrowAspect();
305   //std::cout<<"PrsDim_OffsetDimension::PrsDim_OffsetDimension " <<  myArrowSize << " myArrowSize"<<std::endl;
306   arr->SetLength(myArrowSize);
307   arr = la->ArrowAspect();
308   arr->SetLength(myArrowSize);
309
310   gp_Pnt myTFAttach = myFAttach.Transformed (aTrsf);
311   gp_Pnt myTSAttach = mySAttach.Transformed (aTrsf);
312   gp_Dir myTDirAttach = myDirAttach.Transformed (aTrsf);
313   gp_Dir myTDirAttach2 = myTDirAttach;
314   gp_Pnt Tcurpos = curpos.Transformed (aTrsf);
315  
316   if (myIsSetBndBox)
317     Tcurpos = PrsDim::TranslatePointToBound( Tcurpos, myDirAttach, myBndBox );
318
319   DsgPrs_OffsetPresentation::AddAxes(aprs,
320                                      myDrawer,
321                                      myText,
322                                      myTFAttach,
323                                      myTSAttach,
324                                      myTDirAttach,
325                                      myTDirAttach2,
326                                      Tcurpos);
327   
328   BRepBuilderAPI_Transform transform1 (myFShape, aTrsf, Standard_True);
329   TopoDS_Shape myTFShape = transform1.Shape();
330   BRepBuilderAPI_Transform transform2 (mySShape, aTrsf, Standard_True);
331   TopoDS_Shape myTSShape = transform2.Shape();
332
333   StdPrs_WFShape::Add (aprs, myTFShape, myDrawer);
334   StdPrs_WFShape::Add (aprs, myTSShape, myDrawer);
335 }
336
337 //=======================================================================
338 //function : ComputeTwoFacesOffset
339 //purpose  : 
340 //=======================================================================
341 void PrsDim_OffsetDimension::ComputeTwoFacesOffset(const Handle(Prs3d_Presentation)& aprs, 
342                                                 const gp_Trsf& aTrsf)
343 {
344   gp_Dir norm1 = myDirAttach;
345   gp_Pnt curpos;
346   gp_Ax2 myax2;
347   if (myAutomaticPosition && ! myIsSetBndBox) {
348     TopExp_Explorer explo(myFShape,TopAbs_VERTEX);
349     if (explo.More()) {
350       TopoDS_Vertex vertref = TopoDS::Vertex(explo.Current());
351       myFAttach = BRep_Tool::Pnt(vertref);
352       gp_Vec trans = norm1.XYZ()*fabs(myVal/2);
353       gp_Ax2 ax2(myFAttach,norm1);
354       myDirAttach = ax2.XDirection();
355       curpos = myFAttach.Translated(trans);
356       if (myVal <= Precision::Confusion()) {
357         gp_Vec vecnorm1 = norm1.XYZ()*.001;
358         curpos.Translate(vecnorm1);
359       }
360       myPosition = curpos;
361       myax2 = ax2;
362     }
363   }
364   else {
365     if (myAutomaticPosition && myIsSetBndBox)
366       {
367         Standard_Real aXmin, aYmin, aZmin, aXmax, aYmax, aZmax;
368         myBndBox.Get( aXmin, aYmin, aZmin, aXmax, aYmax, aZmax );
369         myPosition.SetCoord( aXmax, aYmax, aZmax );
370       }
371      
372     curpos = myPosition;
373     myFAttach = PrsDim::Nearest(myFShape,curpos);
374     if (myFAttach.Distance(curpos) <= Precision::Confusion()) {
375       gp_Ax2 ax2(myFAttach,norm1);
376       myDirAttach = ax2.XDirection();
377       myax2 = ax2;
378     }
379     else {
380       gp_Dir orient(myFAttach.XYZ()-curpos.XYZ());
381       gp_Ax2 ax2(myFAttach,norm1);
382       if (orient.Angle(norm1) <= Precision::Angular()) {
383         myDirAttach = ax2.XDirection();
384       }
385       else {
386         gp_Dir adir = norm1 ^ orient;
387         myDirAttach = adir ^ norm1;
388       }
389       myax2 = ax2;
390     }
391   }
392   // en attendant mieux
393   mySAttach = PrsDim::Nearest(mySShape,curpos);
394   gp_Ax3 anax3 (myax2);
395   gp_Pln apln (anax3);
396   
397   //gp_Pnt proj2;
398   Standard_Real u2,v2, uatt, vatt;
399   ElSLib::Parameters (apln , mySAttach, uatt,vatt);
400   ElSLib::Parameters (apln , curpos   , u2,v2);
401   
402   if (uatt== u2 && vatt == v2) {
403     myDirAttach2 = myDirAttach;
404   } else {
405     gp_Vec avec (ElSLib::Value (uatt,vatt, apln) , ElSLib::Value (u2,v2, apln));
406     myDirAttach2.SetCoord (avec.X(),avec.Y(),avec.Z());
407   }
408   
409   Handle(Prs3d_DimensionAspect) la = myDrawer->DimensionAspect();
410   Handle(Prs3d_ArrowAspect) arr = la->ArrowAspect();
411   //std::cout<<"PrsDim_OffsetDimension::PrsDim_OffsetDimension " <<  myArrowSize << " myArrowSize"<<std::endl;
412   arr->SetLength(myArrowSize);
413   arr = la->ArrowAspect();
414   arr->SetLength(myArrowSize);
415
416   gp_Pnt myTFAttach = myFAttach.Transformed (aTrsf);
417   gp_Pnt myTSAttach = mySAttach.Transformed (aTrsf);
418   gp_Dir myTDirAttach = myDirAttach.Transformed (aTrsf);
419   gp_Dir myTDirAttach2 = myDirAttach2.Transformed (aTrsf);
420   gp_Pnt Tcurpos = curpos.Transformed (aTrsf);
421
422 /*
423   if (myIsSetBndBox)
424     {
425       BRepAdaptor_Surface surf1(TopoDS::Face(myFShape));
426       Tcurpos = PrsDim::TranslatePointToBound( Tcurpos, surf1.Plane().XAxis().Direction(), myBndBox );
427     }
428 */
429   DsgPrs_OffsetPresentation::Add(aprs,
430                                  myDrawer,
431                                  myText,
432                                  myTFAttach,
433                                  myTSAttach,
434                                  myTDirAttach,
435                                  myTDirAttach2,
436                                  Tcurpos);
437   
438
439   BRepBuilderAPI_Transform transform1 (myFShape, aTrsf, Standard_True);
440   TopoDS_Shape myTFShape = transform1.Shape();
441   BRepBuilderAPI_Transform transform2 (mySShape, aTrsf, Standard_True);
442   TopoDS_Shape myTSShape = transform2.Shape();
443
444   StdPrs_WFShape::Add (aprs, myTFShape, myDrawer);
445   StdPrs_WFShape::Add (aprs, myTSShape, myDrawer);
446 }
447
448 //=======================================================================
449 //function : ComputeAxeFaceOffset
450 //purpose  : 
451 //=======================================================================
452 void PrsDim_OffsetDimension::ComputeAxeFaceOffset(const Handle(Prs3d_Presentation)& aprs, 
453                                                const gp_Trsf& aTrsf)
454 {
455   BRepBuilderAPI_Transform transform1 (myFShape, aTrsf, Standard_True);
456   TopoDS_Shape myTFShape = transform1.Shape();
457   BRepBuilderAPI_Transform transform2 (mySShape, aTrsf, Standard_True);
458   TopoDS_Shape myTSShape = transform2.Shape();
459
460   StdPrs_WFShape::Add (aprs, myTFShape, myDrawer);
461   StdPrs_WFShape::Add (aprs, myTSShape, myDrawer);
462 }
463