1 /*# ---------------------------------------------------------------------
3 # Copyright (c) CREATIS (Centre de Recherche en Acquisition et Traitement de l'Image
5 # Authors : Eduardo Davila, Frederic Cervenansky, Claire Mouton
6 # Previous Authors : Laurent Guigues, Jean-Pierre Roux
7 # CreaTools website : www.creatis.insa-lyon.fr/site/fr/creatools_accueil
9 # This software is governed by the CeCILL-B license under French law and
10 # abiding by the rules of distribution of free software. You can use,
11 # modify and/ or redistribute the software under the terms of the CeCILL-B
12 # license as circulated by CEA, CNRS and INRIA at the following URL
13 # http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
14 # or in the file LICENSE.txt.
16 # As a counterpart to the access to the source code and rights to copy,
17 # modify and redistribute granted by the license, users are provided only
18 # with a limited warranty and the software's author, the holder of the
19 # economic rights, and the successive licensors have only limited
22 # The fact that you are presently reading this means that you have had
23 # knowledge of the CeCILL-B license and that you accept its terms.
24 # ------------------------------------------------------------------------ */
28 * \brief Class bbtk::ThresholdImageView .
31 #include "LayerImageBase.h"
33 #include "wxMPRBaseData.h"
34 #include "wxVtk2DBaseView.h"
35 #include "wxVTKRenderWindowInteractor.h"
36 #include "vtkTextProperty.h"
39 //---------------------------------------------------------------------------------------------
40 //---------------------------------------------------------------------------------------------
41 //---------------------------------------------------------------------------------------------
43 // Callback for the interaction
44 class ColorLayerObserver : public vtkCommand{
46 virtual char const *GetClassName() const { return "ColorLayerObserver";}
48 static ColorLayerObserver *New() { return new ColorLayerObserver; }
52 _layerimagebase = NULL;
55 virtual void Execute(vtkObject *wdg, unsigned long event, void* calldata)
57 if (_layerimagebase!= NULL)
60 wxVtkBaseView *wxvtkbaseview = _layerimagebase->GetwxVtkBaseView();
61 if (wxvtkbaseview!=NULL)
64 vtkMPRBaseData *vtkmprbasedata = (vtkMPRBaseData*)wxvtkbaseview->GetVtkBaseData();
65 if (vtkmprbasedata!=NULL)
67 _layerimagebase->SetX( vtkmprbasedata->GetX() );
68 _layerimagebase->SetY( vtkmprbasedata->GetY() );
69 _layerimagebase->SetZ( vtkmprbasedata->GetZ() );
70 _layerimagebase->onThresholdChange();
77 void SetLayerImageBase(LayerImageBase* layerimagebase)
79 _layerimagebase = layerimagebase;
83 LayerImageBase* _layerimagebase;
88 //---------------------------------------------------------------------------------------------
89 //---------------------------------------------------------------------------------------------
90 //---------------------------------------------------------------------------------------------
94 LayerImageBase::LayerImageBase()
96 _actorPresent = false;
100 _thresholdTable = NULL;
101 _thresholdMapper = NULL;
102 _scalarbarActor = NULL;
103 _thresholdActor = NULL;
106 _imageReslicer = vtkImageReslice::New();
109 //----------------------------------------------------------------------------
110 LayerImageBase::~LayerImageBase()
114 //----------------------------------------------------------------------------
115 void LayerImageBase::SetX(int x)
120 //----------------------------------------------------------------------------
121 void LayerImageBase::SetY(int y)
126 //----------------------------------------------------------------------------
127 void LayerImageBase::SetZ(int z)
132 //----------------------------------------------------------------------------
133 int LayerImageBase::GetX() // virtual
138 //----------------------------------------------------------------------------
139 int LayerImageBase::GetY() // virtual
143 //----------------------------------------------------------------------------
144 int LayerImageBase::GetZ() // virtual
150 //----------------------------------------------------------------------------
151 vtkImageData* LayerImageBase::GetImage()
156 //----------------------------------------------------------------------------
157 bool LayerImageBase::GetActorPresent()
159 return _actorPresent;
164 //----------------------------------------------------------------------------
165 void LayerImageBase::SetImage(vtkImageData* image)
170 _image->GetSpacing(spc);
171 _image->GetWholeExtent(ext);
172 int sx=ext[1]-ext[0]+1;
173 int sy=ext[3]-ext[2]+1;
174 int sz=ext[5]-ext[4]+1;
175 int ssx=(sx/2)*spc[0];
176 int ssy=(sy/2)*spc[1];
177 int ssz=(sz/2)*spc[2];
182 _image->GetScalarRange( _range );
183 _thresholdTable = vtkLookupTable::New();
186 void LayerImageBase::SetRangeForColorBar(std::vector<double> &range)
188 _rangeForColorBar=range;
191 void LayerImageBase::GetRangeForColorBar(double &minRange, double &maxRange)
193 if (_rangeForColorBar.size()==2)
195 minRange = _rangeForColorBar[0];
196 maxRange = _rangeForColorBar[1];
198 minRange = _range[0];
199 maxRange = _range[1];
205 //----------------------------------------------------------------------------
206 void LayerImageBase::SetwxVtkBaseView(wxVtkBaseView *baseview)
210 _baseView = baseview;
211 ColorLayerObserver *colorlayerobserver = ColorLayerObserver::New();
212 _baseView->GetVtkBaseData()->AddObserver( 1, colorlayerobserver );
213 colorlayerobserver->SetLayerImageBase(this);
218 //----------------------------------------------------------------------------
219 wxVtkBaseView *LayerImageBase::GetwxVtkBaseView()
224 //----------------------------------------------------------------------------
225 void LayerImageBase::Refresh()
229 _baseView->Refresh();
234 //----------------------------------------------------------------------------
235 vtkLookupTable* LayerImageBase::GetThresholdTable()
237 return _thresholdTable;
240 //----------------------------------------------------------------------------
241 void LayerImageBase::CleanXYZ(double &x, double &y, double &z)
244 _image->GetWholeExtent(ext);
261 if ( x > (ext[1]-ext[0]) )
266 if ( y > (ext[3]-ext[2]) )
271 if ( z > (ext[5]-ext[4]) )
279 //----------------------------------------------------------------------------
280 void LayerImageBase::onThreshold()
282 if ((_image!=NULL) && (_baseView!=NULL))
285 GetImage()-> GetSpacing(spc);
295 vtkCamera *camera = _baseView->GetRenderer()->GetActiveCamera();
296 int directionViewer=0;
297 if (camera->GetParallelProjection()==1)
299 wxVtk2DBaseView *wxvtk2dbasevie = (wxVtk2DBaseView *)_baseView;
300 directionViewer = wxvtk2dbasevie->GetDirection();
301 } // ParallelProjection
304 if (!GetActorPresent())
306 if (_thresholdTable==NULL)
309 _thresholdTable = vtkLookupTable::New();
312 if (_thresholdMapper==NULL)
314 _thresholdMapper = vtkImageMapToColors::New( );
317 if (_thresholdActor==NULL)
319 _thresholdActor = vtkImageActor::New( );
320 _scalarbarActor = vtkScalarBarActor::New();
321 _thresholdActor->SetOpacity( 0.6 );
322 _thresholdActor->InterpolateOn( );
323 if (directionViewer==0)
325 _thresholdActor->SetPosition( 900-1,0,0 );
326 _thresholdActor->RotateY( -90 );
328 if (directionViewer==1)
330 _thresholdActor->SetPosition( 0,-900-1, 0 );
331 _thresholdActor->RotateX( 90 );
333 if (directionViewer==2)
335 _thresholdActor->SetPosition( 0,0, -900-1 );
338 _baseView->GetRenderer()->AddActor( _thresholdActor );
339 _baseView->GetRenderer()->AddActor2D( _scalarbarActor );
340 _actorPresent = true;
341 } // !GetActorPresent()
343 ConfigLookupTable(); // virtual method
344 _imageReslicer->SetInput( GetImage() );
345 _imageReslicer->SetInformationInput( GetImage() );
347 if (directionViewer==0)
349 _imageReslicer->SetResliceAxesDirectionCosines( 0,0,1, 0,1,0, -1,0,0 );
350 _imageReslicer->SetResliceAxesOrigin(x,0,0);
353 if (directionViewer==1)
355 _imageReslicer->SetResliceAxesDirectionCosines( 1,0,0 ,0,0,1, 0,-1,0 );
356 _imageReslicer->SetResliceAxesOrigin(0,y,0);
359 if (directionViewer==2)
361 _imageReslicer->SetResliceAxesDirectionCosines( 1,0,0, 0,1,0 ,0,0,1 );
362 _imageReslicer->SetResliceAxesOrigin(0,0,z);
366 _imageReslicer->SetOutputDimensionality(2);
367 // _imageReslicer->SetInterpolationModeToLinear();
368 _imageReslicer->SetInterpolationModeToNearestNeighbor();
370 vtkImageData *img = _imageReslicer->GetOutput();
372 // img->UpdateInformation();
374 // _thresholdTable->Update();
375 _thresholdMapper->SetInput( img );
376 _thresholdMapper->SetLookupTable( _thresholdTable );
377 _thresholdActor->SetInput( _thresholdMapper->GetOutput() );
379 // _scalarbarActor->SetLookupTable( _thresholdTable );
380 _scalarbarActor->SetLookupTable( _thresholdMapper->GetLookupTable() );
381 _scalarbarActor->SetTitle("Value");
382 _scalarbarActor->SetNumberOfLabels(4);
383 int fontsize = _scalarbarActor->GetLabelTextProperty()->GetFontSize();
384 _scalarbarActor->GetLabelTextProperty()->SetFontSize(fontsize/2);
386 // _scalarbarActor->SetTextPad(4); ??
393 //----------------------------------------------------------------------------
394 void LayerImageBase::onThresholdChange()
402 //----------------------------------------------------------------------------
403 void LayerImageBase::onThresholdInterpolation(bool interpolate)
405 if (_thresholdActor!=NULL)
409 _thresholdActor->InterpolateOn( );
413 _thresholdActor->InterpolateOff( );
418 //----------------------------------------------------------------------------
419 void LayerImageBase::onThresholdChangeOpacity (int opacity)
423 _thresholdActor->SetOpacity(opacity*0.1);
427 //----------------------------------------------------------------------------
428 void LayerImageBase::onThresholdRemove()
432 wxVtkBaseView * baseView = _baseView;
433 baseView->GetRenderer()->RemoveActor( _thresholdActor );
434 baseView->GetRenderer()->RemoveActor( _scalarbarActor );
436 _actorPresent = false;
440 //----------------------------------------------------------------------------
441 vtkLookupTable *LayerImageBase::GetvtkLookupTable()
443 return _thresholdTable;
446 //----------------------------------------------------------------------------
447 vtkImageReslice *LayerImageBase::GetvtkImageReslice()
449 return _imageReslicer;
452 //----------------------------------------------------------------------------
453 void LayerImageBase::GetImageScalarRange()
456 GetRangeForColorBar(min,max);
457 // _image->GetScalarRange(_range);