1 // Copyright (C) 2007-2016 CEA/DEN, EDF R&D, OPEN CASCADE
3 // Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
6 // This library is free software; you can redistribute it and/or
7 // modify it under the terms of the GNU Lesser General Public
8 // License as published by the Free Software Foundation; either
9 // version 2.1 of the License, or (at your option) any later version.
11 // This library is distributed in the hope that it will be useful,
12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 // Lesser General Public License for more details.
16 // You should have received a copy of the GNU Lesser General Public
17 // License along with this library; if not, write to the Free Software
18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
23 #include "VTKViewer_Utilities.h"
24 #include "VTKViewer_Actor.h"
25 #include "VTKViewer_Algorithm.h"
31 #include <vtkCamera.h>
32 #include <vtkRenderer.h>
33 #include <vtkRenderWindow.h>
35 /*!@see vtkRenderer::ResetCamera(double bounds[6]) method*/
37 ResetCamera(vtkRenderer* theRenderer,
38 int theUsingZeroFocalPoint)
43 vtkCamera* aCamera = theRenderer->GetActiveCamera();
48 int aCount = ComputeVisiblePropBounds(theRenderer,aBounds);
50 if(theUsingZeroFocalPoint || aCount){
51 static double MIN_DISTANCE = 1.0 / VTK_FLOAT_MAX;
53 double aLength = aBounds[1]-aBounds[0];
54 aLength = std::max((aBounds[3]-aBounds[2]),aLength);
55 aLength = std::max((aBounds[5]-aBounds[4]),aLength);
57 if(aLength < MIN_DISTANCE)
61 sqrt((aBounds[1]-aBounds[0])*(aBounds[1]-aBounds[0]) +
62 (aBounds[3]-aBounds[2])*(aBounds[3]-aBounds[2]) +
63 (aBounds[5]-aBounds[4])*(aBounds[5]-aBounds[4]));
65 if(aWidth < MIN_DISTANCE)
68 double aViewPlaneNormal[3];
69 aCamera->GetViewPlaneNormal(aViewPlaneNormal);
71 double aCenter[3] = {0.0, 0.0, 0.0};
72 if(!theUsingZeroFocalPoint){
73 aCenter[0] = (aBounds[0] + aBounds[1])/2.0;
74 aCenter[1] = (aBounds[2] + aBounds[3])/2.0;
75 aCenter[2] = (aBounds[4] + aBounds[5])/2.0;
77 aCamera->SetFocalPoint(aCenter[0],aCenter[1],aCenter[2]);
79 double aViewAngle = aCamera->GetViewAngle();
80 double aDistance = 2.0*aWidth/tan(aViewAngle*vtkMath::Pi()/360.0);
82 // check view-up vector against view plane normal
84 aCamera->GetViewUp(aViewUp);
85 if(fabs(vtkMath::Dot(aViewUp,aViewPlaneNormal)) > 0.999)
86 aCamera->SetViewUp(-aViewUp[2], aViewUp[0], aViewUp[1]);
89 aCamera->SetPosition(aCenter[0]+aDistance*aViewPlaneNormal[0],
90 aCenter[1]+aDistance*aViewPlaneNormal[1],
91 aCenter[2]+aDistance*aViewPlaneNormal[2]);
93 // find size of the window
94 int* aWinSize = theRenderer->GetSize();
95 if(aWinSize[0] < aWinSize[1])
96 aWidth *= double(aWinSize[1])/double(aWinSize[0]);
98 if(theUsingZeroFocalPoint)
101 aCamera->SetParallelScale(aWidth/2.0);
104 ResetCameraClippingRange(theRenderer);
107 /*! Compute the bounds of the visible props*/
109 ComputeVisiblePropBounds(vtkRenderer* theRenderer,
112 VTK::ActorCollectionCopy aCopy(theRenderer->GetActors());
113 vtkActorCollection* aCollection = aCopy.GetActors();
114 return ComputeBounds( aCollection, theBounds );
117 /*! Compute the bounds of actors*/
119 ComputeBounds(vtkActorCollection* theCollection, double theBounds[6])
123 theBounds[0] = theBounds[2] = theBounds[4] = VTK_FLOAT_MAX;
124 theBounds[1] = theBounds[3] = theBounds[5] = -VTK_FLOAT_MAX;
126 // loop through all props
127 theCollection->InitTraversal();
128 while (vtkActor* aProp = theCollection->GetNextActor()) {
129 // if it's invisible, or has no geometry, we can skip the rest
130 if(aProp->GetVisibility() && aProp->GetMapper() && vtkMath::AreBoundsInitialized(aProp->GetBounds())){
131 if(VTKViewer_Actor* anActor = VTKViewer_Actor::SafeDownCast(aProp))
132 if(anActor->IsInfinitive())
134 double *aBounds = aProp->GetBounds();
135 static double MIN_DISTANCE = 1./VTK_FLOAT_MAX;
136 static double MAX_DISTANCE = 0.9*VTK_FLOAT_MAX;
138 if(fabs(aBounds[1] - aBounds[0]) < MIN_DISTANCE) {
143 if(fabs(aBounds[3] - aBounds[2]) < MIN_DISTANCE) {
148 if(fabs(aBounds[5] - aBounds[4]) < MIN_DISTANCE) {
153 // make sure we haven't got bogus bounds
154 if ( aBounds != NULL &&
155 aBounds[0] > -MAX_DISTANCE && aBounds[1] < MAX_DISTANCE &&
156 aBounds[2] > -MAX_DISTANCE && aBounds[3] < MAX_DISTANCE &&
157 aBounds[4] > -MAX_DISTANCE && aBounds[5] < MAX_DISTANCE)
161 theBounds[0] = std::min(aBounds[0],theBounds[0]);
162 theBounds[2] = std::min(aBounds[2],theBounds[2]);
163 theBounds[4] = std::min(aBounds[4],theBounds[4]);
165 theBounds[1] = std::max(aBounds[1],theBounds[1]);
166 theBounds[3] = std::max(aBounds[3],theBounds[3]);
167 theBounds[5] = std::max(aBounds[5],theBounds[5]);
175 /*!@see vtkRenderer::ResetCameraClippingRange(double bounds[6]) method*/
177 ResetCameraClippingRange(vtkRenderer* theRenderer)
179 if(!theRenderer || !theRenderer->VisibleActorCount()) return;
181 vtkCamera* anActiveCamera = theRenderer->GetActiveCamera();
182 if( anActiveCamera == NULL ){
186 // Find the plane equation for the camera view plane
188 anActiveCamera->GetViewPlaneNormal(vn);
190 anActiveCamera->GetPosition(position);
193 //theRenderer->ComputeVisiblePropBounds(bounds);
194 ComputeVisiblePropBounds(theRenderer, bounds);
197 center[0] = (bounds[0] + bounds[1])/2.0;
198 center[1] = (bounds[2] + bounds[3])/2.0;
199 center[2] = (bounds[4] + bounds[5])/2.0;
201 double width = sqrt((bounds[1]-bounds[0])*(bounds[1]-bounds[0]) +
202 (bounds[3]-bounds[2])*(bounds[3]-bounds[2]) +
203 (bounds[5]-bounds[4])*(bounds[5]-bounds[4]));
205 double distance = sqrt((position[0]-center[0])*(position[0]-center[0]) +
206 (position[1]-center[1])*(position[1]-center[1]) +
207 (position[2]-center[2])*(position[2]-center[2]));
209 double range[2] = {distance - width/2.0, distance + width/2.0};
211 // Do not let the range behind the camera throw off the calculation.
212 if (range[0] < 0.0) range[0] = 0.0;
214 anActiveCamera->SetClippingRange( range );
217 /*!Compute trihedron size.*/
219 ComputeTrihedronSize( vtkRenderer* theRenderer,
221 const double theSize,
222 const double theSizeInPercents )
224 // calculating diagonal of visible props of the renderer
226 if ( ComputeVisiblePropBounds( theRenderer, bnd ) == 0 )
228 bnd[ 1 ] = bnd[ 3 ] = bnd[ 5 ] = 100;
229 bnd[ 0 ] = bnd[ 2 ] = bnd[ 4 ] = 0;
233 aLength = bnd[ 1 ]-bnd[ 0 ];
234 aLength = std::max( ( bnd[ 3 ] - bnd[ 2 ] ),aLength );
235 aLength = std::max( ( bnd[ 5 ] - bnd[ 4 ] ),aLength );
237 static double EPS_SIZE = 5.0E-3;
238 theNewSize = aLength * theSizeInPercents / 100.0;
240 // if the new trihedron size have sufficient difference, then apply the value
241 return fabs( theNewSize - theSize) > theSize * EPS_SIZE ||
242 fabs( theNewSize-theSize ) > theNewSize * EPS_SIZE;
245 bool IsBBEmpty(vtkRenderer* theRenderer)
250 double aNewBndBox[6];
251 aNewBndBox[ 0 ] = aNewBndBox[ 2 ] = aNewBndBox[ 4 ] = VTK_FLOAT_MAX;
252 aNewBndBox[ 1 ] = aNewBndBox[ 3 ] = aNewBndBox[ 5 ] = -VTK_FLOAT_MAX;
254 // iterate through displayed objects and set size if necessary
255 VTK::ActorCollectionCopy aCopy(theRenderer->GetActors());
256 vtkActorCollection* anActors = aCopy.GetActors();
257 anActors->InitTraversal();
259 while(vtkActor* anAct = anActors->GetNextActor())
260 //if(SALOME_Actor* anActor = dynamic_cast<SALOME_Actor*>(anAct))
261 if(VTKViewer_Actor* anActor = VTKViewer_Actor::SafeDownCast(anAct))
262 if(anActor->GetVisibility() && !anActor->IsInfinitive())
264 double *aBounds = anActor->GetBounds();
265 if(aBounds[0] > -VTK_FLOAT_MAX && aBounds[1] < VTK_FLOAT_MAX &&
266 aBounds[2] > -VTK_FLOAT_MAX && aBounds[3] < VTK_FLOAT_MAX &&
267 aBounds[4] > -VTK_FLOAT_MAX && aBounds[5] < VTK_FLOAT_MAX)
275 Check that the given bounding box is valid, i.e each min bound < each max bound
278 bool isBoundValid(double* theBounds) {
279 if(theBounds[0] > theBounds[1] ||
280 theBounds[2] > theBounds[3] ||
281 theBounds[4] > theBounds[5])
287 bool ComputeBBCenter(vtkRenderer* theRenderer, double theCenter[3])
289 theCenter[0] = theCenter[1] = theCenter[2] = 0.0;
294 double aNewBndBox[6];
295 aNewBndBox[ 0 ] = aNewBndBox[ 2 ] = aNewBndBox[ 4 ] = VTK_FLOAT_MAX;
296 aNewBndBox[ 1 ] = aNewBndBox[ 3 ] = aNewBndBox[ 5 ] = -VTK_FLOAT_MAX;
298 // iterate through displayed objects and set size if necessary
299 VTK::ActorCollectionCopy aCopy(theRenderer->GetActors());
300 vtkActorCollection* anActors = aCopy.GetActors();
301 anActors->InitTraversal();
303 while(vtkActor* anAct = anActors->GetNextActor())
305 //if(SALOME_Actor* anActor = dynamic_cast<SALOME_Actor*>(anAct))
306 if(VTKViewer_Actor* anActor = VTKViewer_Actor::SafeDownCast(anAct))
308 if(anActor->GetVisibility() && !anActor->IsInfinitive())
310 double *aBounds = anActor->GetBounds();
312 //Ignore invalid bounds
313 if(!isBoundValid(aBounds)) continue;
315 if(aBounds[0] > -VTK_FLOAT_MAX && aBounds[1] < VTK_FLOAT_MAX &&
316 aBounds[2] > -VTK_FLOAT_MAX && aBounds[3] < VTK_FLOAT_MAX &&
317 aBounds[4] > -VTK_FLOAT_MAX && aBounds[5] < VTK_FLOAT_MAX)
319 for(int i = 0; i < 5; i = i + 2){
320 if(aBounds[i] < aNewBndBox[i])
321 aNewBndBox[i] = aBounds[i];
322 if(aBounds[i+1] > aNewBndBox[i+1])
323 aNewBndBox[i+1] = aBounds[i+1];
333 // null bounding box => the center is (0,0,0)
337 if(aNewBndBox[0] > -VTK_FLOAT_MAX && aNewBndBox[1] < VTK_FLOAT_MAX &&
338 aNewBndBox[2] > -VTK_FLOAT_MAX && aNewBndBox[3] < VTK_FLOAT_MAX &&
339 aNewBndBox[4] > -VTK_FLOAT_MAX && aNewBndBox[5] < VTK_FLOAT_MAX)
341 static double MIN_DISTANCE = 1.0 / VTK_FLOAT_MAX;
343 double aLength = aNewBndBox[1]-aNewBndBox[0];
344 aLength = std::max((aNewBndBox[3]-aNewBndBox[2]),aLength);
345 aLength = std::max((aNewBndBox[5]-aNewBndBox[4]),aLength);
347 if(aLength < MIN_DISTANCE)
351 sqrt((aNewBndBox[1]-aNewBndBox[0])*(aNewBndBox[1]-aNewBndBox[0]) +
352 (aNewBndBox[3]-aNewBndBox[2])*(aNewBndBox[3]-aNewBndBox[2]) +
353 (aNewBndBox[5]-aNewBndBox[4])*(aNewBndBox[5]-aNewBndBox[4]));
355 if(aWidth < MIN_DISTANCE)
358 theCenter[0] = (aNewBndBox[0] + aNewBndBox[1])/2.0;
359 theCenter[1] = (aNewBndBox[2] + aNewBndBox[3])/2.0;
360 theCenter[2] = (aNewBndBox[4] + aNewBndBox[5])/2.0;
368 int aCount = ComputeVisiblePropBounds(theRenderer,aBounds);
369 printf("aNewBndBox[0] = %f, aNewBndBox[1] = %f,\naNewBndBox[2] = %f, aNewBndBox[3] = %f,\naNewBndBox[4] = %f, aNewBndBox[5] = %f\n",
370 aBounds[0],aBounds[1],aBounds[2],aBounds[3],aBounds[4],aBounds[5]);
371 printf("aCount = %d\n",aCount);
374 static double MIN_DISTANCE = 1.0 / VTK_FLOAT_MAX;
376 double aLength = aBounds[1]-aBounds[0];
377 aLength = max((aBounds[3]-aBounds[2]),aLength);
378 aLength = max((aBounds[5]-aBounds[4]),aLength);
380 if(aLength < MIN_DISTANCE)
384 sqrt((aBounds[1]-aBounds[0])*(aBounds[1]-aBounds[0]) +
385 (aBounds[3]-aBounds[2])*(aBounds[3]-aBounds[2]) +
386 (aBounds[5]-aBounds[4])*(aBounds[5]-aBounds[4]));
388 if(aWidth < MIN_DISTANCE)
391 theCenter[0] = (aBounds[0] + aBounds[1])/2.0;
392 theCenter[1] = (aBounds[2] + aBounds[3])/2.0;
393 theCenter[2] = (aBounds[4] + aBounds[5])/2.0;