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();
107 _imageChangeInformation = vtkImageChangeInformation::New();
110 //----------------------------------------------------------------------------
111 LayerImageBase::~LayerImageBase()
115 //----------------------------------------------------------------------------
116 void LayerImageBase::SetX(int x)
121 //----------------------------------------------------------------------------
122 void LayerImageBase::SetY(int y)
127 //----------------------------------------------------------------------------
128 void LayerImageBase::SetZ(int z)
133 //----------------------------------------------------------------------------
134 int LayerImageBase::GetX() // virtual
139 //----------------------------------------------------------------------------
140 int LayerImageBase::GetY() // virtual
144 //----------------------------------------------------------------------------
145 int LayerImageBase::GetZ() // virtual
151 //----------------------------------------------------------------------------
152 vtkImageData* LayerImageBase::GetImage()
155 return _imageChangeInformation->GetOutput();
158 //----------------------------------------------------------------------------
159 bool LayerImageBase::GetActorPresent()
161 return _actorPresent;
165 //----------------------------------------------------------------------------
166 void LayerImageBase::SetImage(vtkImageData* image)
169 int dimensionOriginalLayer[3];
171 _image->GetWholeExtent(ext);
172 dimensionOriginalLayer[0] = ext[1]-ext[0]+1;
173 dimensionOriginalLayer[1] = ext[3]-ext[2]+1;
174 dimensionOriginalLayer[2] = ext[5]-ext[4]+1;
175 SetDimensionOriginalLayer(dimensionOriginalLayer);
176 double spcOriginalLayer[3];
177 _image->GetSpacing(spcOriginalLayer);
178 SetSpcOriginalLayer(spcOriginalLayer);
179 if (GetwxVtkBaseView()!=NULL)
181 vtkImageData *imagebase = GetwxVtkBaseView()->GetVtkBaseData()->GetImageData();
185 imagebase->GetWholeExtent(ext);
186 dimBase[0]=ext[1]-ext[0]+1;
187 dimBase[1]=ext[3]-ext[2]+1;
188 dimBase[2]=ext[5]-ext[4]+1;
189 SetDimensionBase(dimBase);
191 imagebase->GetSpacing(spcBase);
193 SetX( dimBase[0]*spcBase[0]/2 );
194 SetY( dimBase[1]*spcBase[1]/2 );
195 SetZ( dimBase[2]*spcBase[2]/2 );
197 } // if wxvtkbaseview
199 _image->GetScalarRange( _range );
200 _thresholdTable = vtkLookupTable::New();
202 _imageChangeInformation->SetInput(_image);
205 //------------------------------------------------------------------------------
206 void LayerImageBase::SetSpcOriginalLayer(double spc[3])
208 _spcOriginalLayer[0] = spc[0];
209 _spcOriginalLayer[1] = spc[1];
210 _spcOriginalLayer[2] = spc[2];
213 //------------------------------------------------------------------------------
214 void LayerImageBase::SetNewSpacingLayer(double spc[3])
216 _newSpcLayer[0]=spc[0];
217 _newSpcLayer[1]=spc[1];
218 _newSpcLayer[2]=spc[2];
219 _imageChangeInformation->SetOutputSpacing( _newSpcLayer );
224 //------------------------------------------------------------------------------
225 void LayerImageBase::SetSpcBase(double spc[3])
227 _spcBase[0] = spc[0];
228 _spcBase[1] = spc[1];
229 _spcBase[2] = spc[2];
232 //------------------------------------------------------------------------------
233 void LayerImageBase::SetDimensionOriginalLayer(int dim[3])
235 _dimOriginalLayer[0] = dim[0];
236 _dimOriginalLayer[1] = dim[1];
237 _dimOriginalLayer[2] = dim[2];
240 //------------------------------------------------------------------------------
241 void LayerImageBase::SetDimensionBase(int dim[3])
243 _dimBase[0] = dim[0];
244 _dimBase[1] = dim[1];
245 _dimBase[2] = dim[2];
250 //------------------------------------------------------------------------------
251 void LayerImageBase::GetSpcOriginalLayer(double *spc)
253 spc[0] = _spcOriginalLayer[0];
254 spc[1] = _spcOriginalLayer[1];
255 spc[2] = _spcOriginalLayer[2];
258 //------------------------------------------------------------------------------
259 void LayerImageBase::GetSpcBase(double *spc)
261 spc[0] = _spcBase[0];
262 spc[1] = _spcBase[1];
263 spc[2] = _spcBase[2];
266 //------------------------------------------------------------------------------
267 void LayerImageBase::GetDimensionOriginalLayer(int *dim)
269 dim[0] = _dimOriginalLayer[0];
270 dim[1] = _dimOriginalLayer[1];
271 dim[2] = _dimOriginalLayer[2];
274 //------------------------------------------------------------------------------
275 void LayerImageBase::GetDimensionBase(int *dim)
277 dim[0] = _dimBase[0];
278 dim[1] = _dimBase[1];
279 dim[2] = _dimBase[2];
283 //------------------------------------------------------------------------------
284 void LayerImageBase::SetRangeForColorBar(std::vector<double> &range)
286 _rangeForColorBar=range;
290 //------------------------------------------------------------------------------
291 void LayerImageBase::GetRangeForColorBar(double &minRange, double &maxRange)
293 if (_rangeForColorBar.size()==2)
295 minRange = _rangeForColorBar[0];
296 maxRange = _rangeForColorBar[1];
298 minRange = _range[0];
299 maxRange = _range[1];
305 //----------------------------------------------------------------------------
306 void LayerImageBase::SetwxVtkBaseView(wxVtkBaseView *baseview)
310 _baseView = baseview;
311 ColorLayerObserver *colorlayerobserver = ColorLayerObserver::New();
312 _baseView->GetVtkBaseData()->AddObserver( 1, colorlayerobserver );
313 colorlayerobserver->SetLayerImageBase(this);
318 //----------------------------------------------------------------------------
319 wxVtkBaseView *LayerImageBase::GetwxVtkBaseView()
324 //----------------------------------------------------------------------------
325 void LayerImageBase::Refresh()
329 _baseView->Refresh();
334 //----------------------------------------------------------------------------
335 vtkLookupTable* LayerImageBase::GetThresholdTable()
337 return _thresholdTable;
340 //----------------------------------------------------------------------------
341 void LayerImageBase::CleanXYZ(double &x, double &y, double &z)
361 maxSize=_dimOriginalLayer[0]*_newSpcLayer[0];
367 maxSize=_dimOriginalLayer[1]*_newSpcLayer[1];
373 maxSize=_dimOriginalLayer[2]*_newSpcLayer[2];
382 if ( x > (ext[1]-ext[0]) )
387 if ( y > (ext[3]-ext[2]) )
392 if ( z > (ext[5]-ext[4]) )
400 //----------------------------------------------------------------------------
401 void LayerImageBase::onThreshold()
403 if ((_image!=NULL) && (_baseView!=NULL))
406 GetImage()-> GetSpacing(spc);
407 printf("EED LayerImageBase::onThreshold A img-spc %f %f %f\n", spc[0],spc[1],spc[2]);
425 vtkCamera *camera = _baseView->GetRenderer()->GetActiveCamera();
426 int directionViewer=0;
427 if (camera->GetParallelProjection()==1)
429 wxVtk2DBaseView *wxvtk2dbasevie = (wxVtk2DBaseView *)_baseView;
430 directionViewer = wxvtk2dbasevie->GetDirection();
431 } // ParallelProjection
434 if (!GetActorPresent())
436 if (_thresholdTable==NULL)
439 _thresholdTable = vtkLookupTable::New();
442 if (_thresholdMapper==NULL)
444 _thresholdMapper = vtkImageMapToColors::New( );
447 if (_thresholdActor==NULL)
449 _thresholdActor = vtkImageActor::New( );
450 _scalarbarActor = vtkScalarBarActor::New();
451 _thresholdActor->SetOpacity( 0.6 );
452 _thresholdActor->InterpolateOn( );
453 if (directionViewer==0)
455 _thresholdActor->SetPosition( 900-1,0,0 );
456 _thresholdActor->RotateY( -90 );
458 if (directionViewer==1)
460 _thresholdActor->SetPosition( 0,-900-1, 0 );
461 _thresholdActor->RotateX( 90 );
463 if (directionViewer==2)
465 _thresholdActor->SetPosition( 0,0, -900-1 );
468 _baseView->GetRenderer()->AddActor( _thresholdActor );
469 _baseView->GetRenderer()->AddActor2D( _scalarbarActor );
470 _actorPresent = true;
471 } // !GetActorPresent()
473 ConfigLookupTable(); // virtual method
474 _imageReslicer->SetInput( GetImage() );
475 _imageReslicer->SetInformationInput( GetImage() );
477 if (directionViewer==0)
479 _imageReslicer->SetResliceAxesDirectionCosines( 0,0,1, 0,1,0, -1,0,0 );
480 _imageReslicer->SetResliceAxesOrigin(x,0,0);
483 if (directionViewer==1)
485 _imageReslicer->SetResliceAxesDirectionCosines( 1,0,0 ,0,0,1, 0,-1,0 );
486 _imageReslicer->SetResliceAxesOrigin(0,y,0);
489 if (directionViewer==2)
491 _imageReslicer->SetResliceAxesDirectionCosines( 1,0,0, 0,1,0 ,0,0,1 );
492 _imageReslicer->SetResliceAxesOrigin(0,0,z);
495 _imageReslicer->SetOutputDimensionality(2);
496 // _imageReslicer->SetInterpolationModeToLinear();
497 _imageReslicer->SetInterpolationModeToNearestNeighbor();
498 _imageReslicer->Modified();
501 vtkImageData *img = _imageReslicer->GetOutput();
503 // img->UpdateInformation();
505 // _thresholdTable->Update();
506 _thresholdMapper->SetInput( img );
507 _thresholdMapper->SetLookupTable( _thresholdTable );
508 _thresholdActor->SetInput( _thresholdMapper->GetOutput() );
510 // _scalarbarActor->SetLookupTable( _thresholdTable );
511 _scalarbarActor->SetLookupTable( _thresholdMapper->GetLookupTable() );
512 _scalarbarActor->SetTitle("Value");
513 _scalarbarActor->SetNumberOfLabels(4);
514 int fontsize = _scalarbarActor->GetLabelTextProperty()->GetFontSize();
515 _scalarbarActor->GetLabelTextProperty()->SetFontSize(fontsize/2);
517 // _scalarbarActor->SetTextPad(4); ??
519 GetImage()-> GetSpacing(spc);
520 printf("EED LayerImageBase::onThreshold B img-spc %f %f %f\n", spc[0],spc[1],spc[2]);
528 //----------------------------------------------------------------------------
529 void LayerImageBase::onThresholdChange()
537 //----------------------------------------------------------------------------
538 void LayerImageBase::onThresholdInterpolation(bool interpolate)
540 if (_thresholdActor!=NULL)
544 _thresholdActor->InterpolateOn( );
548 _thresholdActor->InterpolateOff( );
553 //----------------------------------------------------------------------------
554 void LayerImageBase::onThresholdChangeOpacity (int opacity)
558 _thresholdActor->SetOpacity(opacity*0.1);
562 //----------------------------------------------------------------------------
563 void LayerImageBase::onThresholdRemove()
567 wxVtkBaseView * baseView = _baseView;
568 baseView->GetRenderer()->RemoveActor( _thresholdActor );
569 baseView->GetRenderer()->RemoveActor( _scalarbarActor );
571 _actorPresent = false;
575 //----------------------------------------------------------------------------
576 vtkLookupTable *LayerImageBase::GetvtkLookupTable()
578 return _thresholdTable;
581 //----------------------------------------------------------------------------
582 vtkImageReslice *LayerImageBase::GetvtkImageReslice()
584 return _imageReslicer;
587 //----------------------------------------------------------------------------
588 void LayerImageBase::GetImageScalarRange()
591 GetRangeForColorBar(min,max);
592 // _image->GetScalarRange(_range);