/*# --------------------------------------------------------------------- # # Copyright (c) CREATIS (Centre de Recherche en Acquisition et Traitement de l'Image # pour la Sant�) # Authors : Eduardo Davila, Frederic Cervenansky, Claire Mouton # Previous Authors : Laurent Guigues, Jean-Pierre Roux # CreaTools website : www.creatis.insa-lyon.fr/site/fr/creatools_accueil # # This software is governed by the CeCILL-B license under French law and # abiding by the rules of distribution of free software. You can use, # modify and/ or redistribute the software under the terms of the CeCILL-B # license as circulated by CEA, CNRS and INRIA at the following URL # http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html # or in the file LICENSE.txt. # # As a counterpart to the access to the source code and rights to copy, # modify and redistribute granted by the license, users are provided only # with a limited warranty and the software's author, the holder of the # economic rights, and the successive licensors have only limited # liability. # # The fact that you are presently reading this means that you have had # knowledge of the CeCILL-B license and that you accept its terms. # ------------------------------------------------------------------------ */ /** * \file * \brief Class bbtk::ThresholdImageView . */ #include "LayerImageBase.h" #include "wxMPRBaseData.h" #include "wxVtk2DBaseView.h" #include "wxVTKRenderWindowInteractor.h" //--------------------------------------------------------------------------------------------- //--------------------------------------------------------------------------------------------- //--------------------------------------------------------------------------------------------- // Callback for the interaction class ColorLayerObserver : public vtkCommand{ public: virtual char const *GetClassName() const { return "ColorLayerObserver";} static ColorLayerObserver *New() { return new ColorLayerObserver; } ColorLayerObserver() { _layerimagebase = NULL; } virtual void Execute(vtkObject *wdg, unsigned long event, void* calldata) { if (_layerimagebase!= NULL) { wxVtkBaseView *wxvtkbaseview = _layerimagebase->GetwxVtkBaseView(); if (wxvtkbaseview!=NULL) { vtkMPRBaseData *vtkmprbasedata = (vtkMPRBaseData*)wxvtkbaseview->GetVtkBaseData(); if (vtkmprbasedata!=NULL) { _layerimagebase->SetX( vtkmprbasedata->GetX() ); _layerimagebase->SetY( vtkmprbasedata->GetY() ); _layerimagebase->SetZ( vtkmprbasedata->GetZ() ); _layerimagebase->onThresholdChange(); } // if basedata } // if wxvtkbaseview } //if layerimagebase } void SetLayerImageBase(LayerImageBase* layerimagebase) { _layerimagebase = layerimagebase; } private: LayerImageBase* _layerimagebase; }; //--------------------------------------------------------------------------------------------- //--------------------------------------------------------------------------------------------- //--------------------------------------------------------------------------------------------- ColorLayerImageBasevtkInteractor::ColorLayerImageBasevtkInteractor(LayerImageBase* layerImageBase) { _state = false; _stateKey = false; _layerImageBase = layerImageBase; } //--------------------------------------------------------------------------------------------- ColorLayerImageBasevtkInteractor::~ColorLayerImageBasevtkInteractor() { } //--------------------------------------------------------------------------------------------- bool ColorLayerImageBasevtkInteractor::OnChar() { return true; } //--------------------------------------------------------------------------------------------- bool ColorLayerImageBasevtkInteractor::OnMouseMove() // vitual { printf("EED ColorLayerImageBasevtkInteractor::OnMouseMove start \n"); vtkRenderWindowInteractor *interactor = _vtkInteractorStyleBaseView->GetInteractor(); // For example in ManualPaint Action // This is not working, because the order of insertion events mecanism (_vtkInteractorStyleBaseView) // Layer first Paint after. // the only way to make it works is to take the axes an paint. this force a first SetParent_refresh_waiting() // if ( this->_vtkInteractorStyleBaseView->GetParent_refresh_waiting() == true ) if ( interactor->GetControlKey()==1 ) // EED This if for ManaulPaint Interaction compatibility { printf("EED ColorLayerImageBasevtkInteractor::OnMouseMove need Refresh \n" ); _layerImageBase->onThreshold(); _layerImageBase->GetvtkImageReslice()->Modified(); } else { printf("EED ColorLayerImageBasevtkInteractor::OnMouseMove NOT need Refresh \n" ); } if ( (_state==true) && (_layerImageBase->GetActorPresent()) ) { wxVtk2DBaseView *wxvtk2Dbaseview = (wxVtk2DBaseView*) _vtkInteractorStyleBaseView->GetWxVtk2DBaseView(); wxVTKRenderWindowInteractor *wxVTKiren = wxvtk2Dbaseview->GetWxVTKRenderWindowInteractor(); int px, py; wxVTKiren->GetEventPosition(px, py); // int typeView = 1; double X = (double) px; double Y = (double) py; double Z = wxvtk2Dbaseview->GetActualSlice(); // wxvtk2Dbaseview->TransFromCoordScreenToWorld(X, Y, Z,true, typeView); wxvtk2Dbaseview->TransFromCoordScreenToWorld(X, Y, Z, true); // _manualPaintControler->PaintImage((int) X, (int) Y, (int) Z); _layerImageBase->SetX( (int)X ); _layerImageBase->SetY( (int)Y ); _layerImageBase->SetZ( (int)Z ); _layerImageBase->onThreshold(); // _layerImageBase->Refresh(); this->_vtkInteractorStyleBaseView->SetParent_refresh_waiting(); if (_vtkInteractorStyleBaseView->GetRefresh_waiting()==true) { _layerImageBase->GetvtkImageReslice()->Modified(); } printf("EED ColorLayerImageBasevtkInteractor::YESSSSSS \n"); } return true; } //--------------------------------------------------------------------------------------------- bool ColorLayerImageBasevtkInteractor::OnRightButtonDown() // virtual { printf("EED ColorLayerImageBasevtkInteractor::OnRightButtonDown \n"); _state=true; return true; } //--------------------------------------------------------------------------------------------- bool ColorLayerImageBasevtkInteractor::OnRightButtonUp() // virtual { printf("EED ColorLayerImageBasevtkInteractor::OnRightButtonUp \n"); _state=false; return true; } //--------------------------------------------------------------------------------------------- //--------------------------------------------------------------------------------------------- //--------------------------------------------------------------------------------------------- LayerImageBase::LayerImageBase() { _actorPresent = false; _X = 0; _Y = 0; _Z = 0; _thresholdTable = NULL; _thresholdMapper = NULL; _scalarbarActor = NULL; _thresholdActor = NULL; _image = NULL; _baseView = NULL; _imageReslicer = vtkImageReslice::New(); } //---------------------------------------------------------------------------- LayerImageBase::~LayerImageBase() { } //---------------------------------------------------------------------------- void LayerImageBase::SetX(int x) { _X = x; } //---------------------------------------------------------------------------- void LayerImageBase::SetY(int y) { _Y = y; } //---------------------------------------------------------------------------- void LayerImageBase::SetZ(int z) { _Z = z; } //---------------------------------------------------------------------------- int LayerImageBase::GetX() // virtual { return _X; } //---------------------------------------------------------------------------- int LayerImageBase::GetY() // virtual { return _Y; } //---------------------------------------------------------------------------- int LayerImageBase::GetZ() // virtual { return _Z; } //---------------------------------------------------------------------------- vtkImageData* LayerImageBase::GetImage() { return _image; } //---------------------------------------------------------------------------- bool LayerImageBase::GetActorPresent() { return _actorPresent; } //---------------------------------------------------------------------------- void LayerImageBase::SetImage(vtkImageData* image) { _image = image; } //---------------------------------------------------------------------------- void LayerImageBase::SetwxVtkBaseView(wxVtkBaseView *baseview) { if (_baseView==NULL) { _baseView = baseview; // Old EED 2014 11 13 // vtkInteractorStyleBaseView *isbv = (vtkInteractorStyleBaseView*)(_baseView->GetInteractorStyleBaseView()); // isbv->AddInteractorStyleMaracas( new ColorLayerImageBasevtkInteractor(this) ); ColorLayerObserver *colorlayerobserver = ColorLayerObserver::New(); _baseView->GetVtkBaseData()->AddObserver( 1, colorlayerobserver ); colorlayerobserver->SetLayerImageBase(this); } // if _baseView } //---------------------------------------------------------------------------- wxVtkBaseView *LayerImageBase::GetwxVtkBaseView() { return _baseView; } //---------------------------------------------------------------------------- void LayerImageBase::Refresh() { if (_baseView!=NULL) { _baseView->Refresh(); } } //---------------------------------------------------------------------------- vtkLookupTable* LayerImageBase::GetThresholdTable() { return _thresholdTable; } //---------------------------------------------------------------------------- void LayerImageBase::CleanXYZ(double &x, double &y, double &z) { int ext[6]; _image->GetWholeExtent(ext); if (x<0) { x=0; } if (y<0) { y=0; } if (z<0) { z=0; } if ( x > (ext[1]-ext[0]) ) { z=ext[1]-ext[0]; } if ( y > (ext[3]-ext[2]) ) { z=ext[3]-ext[2]; } if ( z > (ext[5]-ext[4]) ) { z=ext[5]-ext[4]; } } //---------------------------------------------------------------------------- void LayerImageBase::onThreshold() { if ((_image!=NULL) && (_baseView!=NULL)) { double spc[3]; GetImage()-> GetSpacing(spc); double x=GetX(); double y=GetY(); double z=GetZ(); CleanXYZ(x,y,z); x = x*spc[0]; y = y*spc[1]; z = z*spc[2]; vtkCamera *camera = _baseView->GetRenderer()->GetActiveCamera(); int directionViewer=0; if (camera->GetParallelProjection()==1) { wxVtk2DBaseView *wxvtk2dbasevie = (wxVtk2DBaseView *)_baseView; directionViewer = wxvtk2dbasevie->GetDirection(); } // ParallelProjection if (!GetActorPresent()) { if (_thresholdTable==NULL) { //Lookup Table _thresholdTable = vtkLookupTable::New(); } // _thresholdTable if (_thresholdMapper==NULL) { _thresholdMapper = vtkImageMapToColors::New( ); } if (_thresholdActor==NULL) { _thresholdActor = vtkImageActor::New( ); _scalarbarActor = vtkScalarBarActor::New(); _thresholdActor->SetOpacity( 0.6 ); _thresholdActor->InterpolateOn( ); if (directionViewer==0) { _thresholdActor->SetPosition( 900-1,0,0 ); _thresholdActor->RotateY( -90 ); } if (directionViewer==1) { _thresholdActor->SetPosition( 0,-900-1, 0 ); _thresholdActor->RotateX( 90 ); } if (directionViewer==2) { _thresholdActor->SetPosition( 0,0, -900-1 ); } } // _thresholdActor _baseView->GetRenderer()->AddActor( _thresholdActor ); _baseView->GetRenderer()->AddActor2D( _scalarbarActor ); _actorPresent = true; } // !GetActorPresent() ConfigLookupTable(); // virtual method _imageReslicer->SetInput( GetImage() ); _imageReslicer->SetInformationInput( GetImage() ); if (directionViewer==0) { _imageReslicer->SetResliceAxesDirectionCosines( 0,0,1, 0,1,0, -1,0,0 ); _imageReslicer->SetResliceAxesOrigin(x,0,0); } if (directionViewer==1) { _imageReslicer->SetResliceAxesDirectionCosines( 1,0,0 ,0,0,1, 0,-1,0 ); _imageReslicer->SetResliceAxesOrigin(0,y,0); } if (directionViewer==2) { _imageReslicer->SetResliceAxesDirectionCosines( 1,0,0, 0,1,0 ,0,0,1 ); _imageReslicer->SetResliceAxesOrigin(0,0,z); } _imageReslicer->SetOutputDimensionality(2); // _imageReslicer->SetInterpolationModeToLinear(); _imageReslicer->SetInterpolationModeToNearestNeighbor(); vtkImageData *img = _imageReslicer->GetOutput(); // img->Update(); // img->UpdateInformation(); // _thresholdTable->Update(); _thresholdMapper->SetInput( img ); _thresholdMapper->SetLookupTable( _thresholdTable ); _thresholdActor->SetInput( _thresholdMapper->GetOutput() ); // _scalarbarActor->SetLookupTable( _thresholdTable ); _scalarbarActor->SetLookupTable( _thresholdMapper->GetLookupTable() ); _scalarbarActor->SetTitle("Value"); _scalarbarActor->SetNumberOfLabels(5); // _scalarbarActor->SetTextPad(4); ?? } // _image } //---------------------------------------------------------------------------- void LayerImageBase::onThresholdChange() { if (_actorPresent) { onThreshold(); } } //---------------------------------------------------------------------------- void LayerImageBase::onThresholdInterpolation(bool interpolate) { if (_thresholdActor!=NULL) { if (interpolate) { _thresholdActor->InterpolateOn( ); } else { _thresholdActor->InterpolateOff( ); } } } //---------------------------------------------------------------------------- void LayerImageBase::onThresholdChangeOpacity (int opacity) { if (_actorPresent) { _thresholdActor->SetOpacity(opacity*0.1); } } //---------------------------------------------------------------------------- void LayerImageBase::onThresholdRemove() { if (_actorPresent) { wxVtkBaseView * baseView = _baseView; baseView->GetRenderer()->RemoveActor( _thresholdActor ); baseView->GetRenderer()->RemoveActor( _scalarbarActor ); _actorPresent = false; } } //---------------------------------------------------------------------------- vtkLookupTable *LayerImageBase::GetvtkLookupTable() { return _thresholdTable; } //---------------------------------------------------------------------------- vtkImageReslice *LayerImageBase::GetvtkImageReslice() { return _imageReslicer; } // EOF