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