X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=lib%2FmaracasVisuLib%2Fsrc%2Finterface%2FwxWindows%2Fwidgets%2FwxVtkMPR3DView.cxx;h=108b5913765e923e3c7e8d1d0eea9b9868dbc690;hb=b80b4dcc301c1c2c63c0f76cd06b0a434350f7a7;hp=28c4ca778c3eb9885e99a99a0571a78c5984dd11;hpb=5c61ebf06cedf398ed1b3412899366d16a8b079c;p=creaMaracasVisu.git

diff --git a/lib/maracasVisuLib/src/interface/wxWindows/widgets/wxVtkMPR3DView.cxx b/lib/maracasVisuLib/src/interface/wxWindows/widgets/wxVtkMPR3DView.cxx
index 28c4ca7..108b591 100644
--- a/lib/maracasVisuLib/src/interface/wxWindows/widgets/wxVtkMPR3DView.cxx
+++ b/lib/maracasVisuLib/src/interface/wxWindows/widgets/wxVtkMPR3DView.cxx
@@ -1,3 +1,28 @@
+/*# ---------------------------------------------------------------------
+#
+# 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.
+# ------------------------------------------------------------------------ */
+
 #include "wxVtkMPR3DView.h"
 
 #include "vtkDataSetMapper.h"
@@ -5,7 +30,6 @@
 #include "vtkProbeFilter.h"
 #include "vtkStructuredPoints.h"
 #include "vtkMetaImageWriter.h"
-#include "vtkDataSetMapper.h"
 #include "vtkImageChangeInformation.h"
 #include "vtkPointData.h"
 #include "vtkWindowLevelLookupTable.h"
@@ -13,89 +37,174 @@
 #include "vtkCutter.h"
 #include "vtkTransform.h"
 #include "vtkmyPWCallback_3DPointWidget.h"
+#include "vtkProperty.h"
+#include "vtkCellPicker.h"
+#include "vtkInteractorStyleBaseView.h"
+#include <vtkImageReslice.h>
+
+
+
+//----------------------------------------------------------------------------
+//----------------------------------------------------------------------------
+//----------------------------------------------------------------------------
+class vtkWindowLevelImagePlaneWidgetCallback : public vtkCommand
+{
+	 public:
+	   static vtkWindowLevelImagePlaneWidgetCallback *New()
+	   { 
+		   return new vtkWindowLevelImagePlaneWidgetCallback;
+	   }
+	
+	   void Execute( vtkObject *caller, unsigned long ev, void *callData )
+		{
+			bool needRefresh=false;
+			double slicePosition;
+			vtkImagePlaneWidget* ipw =dynamic_cast< vtkImagePlaneWidget* >( caller );
+		    if ( ev == vtkCommand::WindowLevelEvent)
+			{
+				double info[2];
+				ipw->GetWindowLevel(info);
+				_wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->SetColorWindow(info[0]);
+				_wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->SetColorLevel(info[1]);
+				needRefresh=true;
+			} // if ev
+		   
+		   if ( ev == vtkCommand::InteractionEvent)
+		   {
+			   needRefresh=true;
+			   if (ipw->GetCursorDataStatus())
+			   {
+				   double xyzv[4];
+				   ipw->GetCursorData(xyzv);
+				   _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->SetX(xyzv[0]);
+				   _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->SetY(xyzv[1]);
+				   _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->SetZ(xyzv[2]);
+			   } else {
+				   
+				   
+				   vtkImageData *image = _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData();
+				   double spc[3];
+				   if(image)
+				   {
+					   image->GetSpacing(spc);
+				   }
+				   
+				   double normal[3];
+				   slicePosition = ipw->GetSlicePosition();
+				   ipw->GetNormal(normal);
+				   
+				   // FreePlaneX
+				   if ((normal[0]==1)&&(normal[1]==0)&&(normal[2]==0))
+				   {
+					   _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->SetX( slicePosition/spc[0] );
+				   }
+				   // FreePlaneY
+				   if ((normal[0]==0)&&(normal[1]==1)&&(normal[2]==0))
+				   {
+					   _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->SetY( slicePosition/spc[1] );
+				   }
+				   // FreePlaneZ
+				   if ((normal[0]==0)&&(normal[1]==0)&&(normal[2]==1))
+				   {
+					   _wxvtkmpr3Dview->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->SetZ( slicePosition/spc[2] );
+				   }
+			   }// ipw->GetCursorDataStatus
+		   }// ev
+			
+		   if (needRefresh)	
+		   {
+			   vtkInteractorStyleBaseView *isbv = (vtkInteractorStyleBaseView*)_wxvtkmpr3Dview->GetWxvtk3Dbaseview()->GetInteractorStyleBaseView();
+			   isbv->SetParent_refresh_waiting();
+			   isbv->EvaluateToRefresh();
+		   }
+			
+	   } // Execute
+	
+		vtkWindowLevelImagePlaneWidgetCallback(){}
+		wxVtkMPR3DView *_wxvtkmpr3Dview;
+};
+
+
+
+//-------------------------------------------------------------------
+//-------------------------------------------------------------------
+//-------------------------------------------------------------------
+
+
 
 wxVtkMPR3DView::wxVtkMPR3DView( wxVtk3DBaseView *wxvtk3Dbaseview )
 {
 	_wxvtk3Dbaseview			=	wxvtk3Dbaseview;
 	_vtkmpr3Ddataviewer			=	NULL;
-	_wxvtkmpr3DviewCntrlPanel	=	NULL;	
-	//_myCallback=NULL;
-	_pointWidget=NULL;
-	_planeWidget=NULL;
-	_vtkplane=NULL;
-	_probe=NULL;
-	_contourMapper=NULL;
+	_wxvtkmpr3DviewCntrlPanel	=	NULL;
+	//_myCallback				=	NULL;
+	_pointWidget				=	NULL;
+	_planeWidget				=	NULL;
+	_vtkplane					=	NULL;
+	_probe						=	NULL;
+	_contourMapper				=	NULL;
+
+	_planeWidgetX				=	NULL;
+	_planeWidgetY				=	NULL;
+	_planeWidgetZ				=	NULL;
+	//EED Fev 3 2010
+	//Free Planes
 }
-void wxVtkMPR3DView::RemoveActor(vtkActor* actor){	
 
-	_wxvtk3Dbaseview->GetRenderer()->RemoveActor(actor);	
 
+//-------------------------------------------------------------------
+wxVtkMPR3DView::~wxVtkMPR3DView()
+{
+	VisiblePointWidget( false );
+	VisiblePlaneWidget( false );
 
-}
-void wxVtkMPR3DView::ResetCamera(int* ext, double* origin,double* spc){
-	if(ext == NULL){
+	if(_pointWidget!=NULL){
+		_pointWidget->Delete();
+	}
+	if(_planeWidget!=NULL){
+		_planeWidget->Delete();
+	}
 
-		_wxvtk3Dbaseview->GetRenderer()->ResetCamera ();	
-		_wxvtk3Dbaseview->GetCamera()->Dolly(1.5);	
-	}else{
-		/*double x = (spc[0])*(origin[0]+(((double)ext[1]-(double)ext[0])/2.0));		
-		double y = (spc[1])*(origin[1]+(double)ext[3]);
-		double z = (spc[2])*(origin[2]+(((double)ext[5]-(double)ext[4])/2.0));*/
-		/*double x0=(spc[0])*((double)ext[0]+origin[0]);
-		double x1=(spc[0])*((double)ext[1]+origin[0]);
-		double y0=(spc[1])*((double)ext[2]+origin[1]);
-		double y1=(spc[1])*((double)ext[3]+origin[1]);
-		double z0=(spc[2])*((double)ext[4]+origin[2]);
-		double z1=(spc[2])*((double)ext[5]+origin[2]);*/
-		double x0=(spc[0])*((double)ext[0]);
-		double x1=(spc[0])*((double)ext[1]);
-		double y0=(spc[1])*((double)ext[2]);
-		double y1=(spc[1])*((double)ext[3]);
-		double z0=(spc[2])*((double)ext[4]);
-		double z1=(spc[2])*((double)ext[5]);
-		
-		_wxvtk3Dbaseview->GetRenderer()->ResetCamera(x0,x1,y0,y1,z0,z1);
-		//_wxvtk3Dbaseview->GetCamera()->SetPosition(x,y,z);
-		_wxvtk3Dbaseview->GetCamera()->Dolly(1.5);	
+	if(_vtkplane!=NULL){
+		_vtkplane->Delete();
+		_probe->Delete();
+		_contourMapper->Delete();
+		_contourPlaneActor->Delete();
 	}
+
+}
+
+//-------------------------------------------------------------------
+void wxVtkMPR3DView::RemoveActor(vtkActor* actor)
+{
+	_wxvtk3Dbaseview->GetRenderer()->RemoveActor(actor);
 }
-void wxVtkMPR3DView::Configure(){
 
 
+
+//-------------------------------------------------------------------
+void wxVtkMPR3DView::Configure()
+{
+	vtkImageData *imageData = GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData();
+
 	_wxvtk3Dbaseview->Configure();
 	_wxvtkmpr3DviewCntrlPanel->UpdateControlPanel();
-	
-//EED 27 Mai 2009 
-//	_wxvtk3Dbaseview->GetRenderer()->Clear();    
-	
-	// Actors are added to the renderer. 
+
+//EED 27 Mai 2009
+//	_wxvtk3Dbaseview->GetRenderer()->Clear();
+
+	// Actors are added to the renderer.
 	vtkActor* _outlineActor = _vtkmpr3Ddataviewer->GetOutlineActor();
-	_wxvtk3Dbaseview->GetRenderer()->AddActor( _outlineActor );			
+
+
+	_wxvtk3Dbaseview->GetRenderer()->AddActor( _outlineActor );
+
 //	_wxvtk3Dbaseview->GetRenderer()->AddActor( _vtkmpr3Ddataviewer->GetImageActor(0)  );	// _saggital
 //	_wxvtk3Dbaseview->GetRenderer()->AddActor( _vtkmpr3Ddataviewer->GetImageActor(1)  );	// _axial
 //	_wxvtk3Dbaseview->GetRenderer()->AddActor( _vtkmpr3Ddataviewer->GetImageActor(2)  );	// _coronal
 
 
-	// An initial camera view is created.  The Dolly() method moves 
-	// the camera towards the FocalPoint, thereby enlarging the image.	
-	_wxvtk3Dbaseview->GetRenderer()->SetActiveCamera(_wxvtk3Dbaseview->GetCamera());
-	this->ResetCamera();
 	
-
-	// Set a background color for the renderer and set the size of the
-	// render window (expressed in pixels).
-	_wxvtk3Dbaseview->GetRenderer()->SetBackground( 0.36 , 0.36 , 0.36 );
-	_wxvtk3Dbaseview->GetRenWin()->SetSize(400, 400);
-
-	// Note that when camera movement occurs (as it does in the Dolly()
-	// method), the clipping planes often need adjusting. Clipping planes
-	// consist of two planes: near and far along the view direction. The 
-	// near plane clips out objects in front of the plane; the far plane
-	// clips out objects behind the plane. This way only what is drawn
-	// between the planes is actually rendered.
-	_wxvtk3Dbaseview->GetRenderer()->ResetCameraClippingRange();
-
-
 	// vtkPointWidget
 	//if(_myCallback!=NULL){
 		//_myCallback->Delete();
@@ -103,162 +212,376 @@ void wxVtkMPR3DView::Configure(){
 	//}
 	vtkmyPWCallback_3DPointWidget *_myCallback = vtkmyPWCallback_3DPointWidget::New();
 	_myCallback->SetWxVtkMPR3DView(this);
- 
-	if(_pointWidget!=NULL){
-		_pointWidget->Delete();
+
+
+	if(imageData){
+		if(_pointWidget==NULL){
+			_pointWidget = vtkPointWidget::New();
+		}
+		_pointWidget->SetInput( imageData );
+		_myCallback->SetVtkPointWidget(_pointWidget);
+		_pointWidget->SetInteractor( GetWxvtk3Dbaseview()->GetWxVTKRenderWindowInteractor() );
+		_pointWidget->AllOff();
+		_pointWidget->PlaceWidget();
+
+		_pointWidget->AddObserver(vtkCommand::InteractionEvent,_myCallback);
+
+		if(_planeWidget==NULL)
+		{
+			_planeWidget = vtkPlaneWidget::New();
+		}
+		_planeWidget->SetInput( imageData );
+
+		_myCallback->SetVtkPlaneWidget(_planeWidget);
+
+		_planeWidget->NormalToXAxisOn();
+		_planeWidget->SetResolution(50);
+		_planeWidget->SetRepresentationToOutline();
+
+		int dim[3];
+
+		imageData->GetDimensions(dim);
+		int px=(dim[0]/2);
+		int py=(dim[1]/2);
+		int pz=(dim[2]/2);
+		int dd=20;
+		_planeWidget->PlaceWidget( px-dd , px+dd , py-dd , py+dd , pz-dd , pz+dd );
+
+		if(_vtkplane==NULL){
+			_vtkplane = vtkPolyData::New();
+
+			_probe = vtkProbeFilter::New();
+			_probe->SetInput(_vtkplane);
+
+			_contourMapper = vtkPolyDataMapper::New();
+
+			_contourMapper->SetInput( _probe->GetPolyDataOutput() );
+
+			_contourPlaneActor = vtkActor::New();
+			_contourPlaneActor->SetMapper(_contourMapper);
+
+			_contourPlaneActor->VisibilityOff();
+
+			_planeWidget->SetInteractor( GetWxvtk3Dbaseview()->GetWxVTKRenderWindowInteractor() );
+			_planeWidget->AddObserver(vtkCommand::InteractionEvent,_myCallback);
+
+			_planeWidget->Off();
+			_wxvtk3Dbaseview->GetRenderer()->AddActor( _contourPlaneActor );
+
+		}
+
+		_planeWidget->GetPolyData(_vtkplane);
+		_probe->SetSource( imageData );
+		_contourMapper->SetScalarRange( imageData->GetScalarRange() );
+
+	ConfigureFreePlanes();
 	}
-	_pointWidget = vtkPointWidget::New();
-	_myCallback->SetVtkPointWidget(_pointWidget);
-	_pointWidget->SetInteractor( GetWxvtk3Dbaseview()->GetWxVTKRenderWindowInteractor() ); 
-	_pointWidget->SetInput( GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData() );
-	
-	_pointWidget->AllOff();
-	
-	_pointWidget->PlaceWidget();
-//	_pointWidget->SetPosition( x,y,z ); 
-	
+
+//	_pointWidget->SetPosition( x,y,z );
+
 //EED01 29Mars2009
-// MACOS probleme vtk-window out of wx-window 	
+// MACOS probleme vtk-window out of wx-window
 //	_pointWidget->On();
-	
-	_pointWidget->AddObserver(vtkCommand::InteractionEvent,_myCallback);
-
 
 // PlaneWidget
-	if(_planeWidget!=NULL){
-		_planeWidget->Delete();
-	}
-	_planeWidget = vtkPlaneWidget::New();
-	_myCallback->SetVtkPlaneWidget(_planeWidget);
-	_planeWidget->SetInput( GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData() );
-	_planeWidget->NormalToXAxisOn();
-	_planeWidget->SetResolution(50);
-	_planeWidget->SetRepresentationToOutline();
-	int dim[3];
-	GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData()->GetDimensions(dim);
-	int px=(dim[0]/2);
-	int py=(dim[1]/2);
-	int pz=(dim[2]/2);
-	int dd=20;
-	_planeWidget->PlaceWidget( px-dd , px+dd , py-dd , py+dd , pz-dd , pz+dd );
-
-
-	if(_vtkplane==NULL){
-		_vtkplane = vtkPolyData::New();
-		
-		_probe = vtkProbeFilter::New();
-		_probe->SetInput(_vtkplane);
 
-		_contourMapper = vtkPolyDataMapper::New();
-		
-		_contourMapper->SetInput( _probe->GetPolyDataOutput() );
 
-		_contourPlaneActor = vtkActor::New();
-		_contourPlaneActor->SetMapper(_contourMapper);
+// An initial camera view is created.  The Dolly() method moves
+	// the camera towards the FocalPoint, thereby enlarging the image.
+//	_wxvtk3Dbaseview->GetRenderer()->SetActiveCamera(_wxvtk3Dbaseview->GetCamera());
 
-		_contourPlaneActor->VisibilityOff();
+	// Note that when camera movement occurs (as it does in the Dolly()
+	// method), the clipping planes often need adjusting. Clipping planes
+	// consist of two planes: near and far along the view direction. The
+	// near plane clips out objects in front of the plane; the far plane
+	// clips out objects behind the plane. This way only what is drawn
+	// between the planes is actually rendered.
 
-		_planeWidget->SetInteractor( GetWxvtk3Dbaseview()->GetWxVTKRenderWindowInteractor() );
-		_planeWidget->AddObserver(vtkCommand::InteractionEvent,_myCallback);
+//EED 27/05/2013
+//	_wxvtk3Dbaseview->GetRenderer()->ResetCameraClippingRange();
 
-		_planeWidget->Off();
-		_wxvtk3Dbaseview->GetRenderer()->AddActor( _contourPlaneActor );	
+//EED 27/05/2013
+//      this->ResetCamera();
+	_wxvtk3Dbaseview->ResetCamera();
 
-	}else{
-		
-	}
-	_planeWidget->GetPolyData(_vtkplane);
+	// Set a background color for the renderer and set the size of the
+	// render window (expressed in pixels).
+	_wxvtk3Dbaseview->GetRenderer()->SetBackground( 0.36 , 0.36 , 0.36 );
+//EED 23oct2010	_wxvtk3Dbaseview->GetRenWin()->SetSize(400, 400);
 
-	
-	vtkImageData *imageData = GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData();	
-	_probe->SetSource( imageData );
+}
 
+//-------------------------------------------------------------------
+void wxVtkMPR3DView::ConfigureFreePlanes()
+{
+
+	// The shared picker enables us to use 3 planes at one time
+	// and gets the picking order right
+	vtkCellPicker* picker = vtkCellPicker::New();
+	picker->SetTolerance(0.005);
+
+	// The 3 image plane widgets
+	_planeWidgetX = GetPlaneWidget('x', 1, 0, 0, picker);
+	_planeWidgetY = GetPlaneWidget('y', 1, 1, 0, picker);
+	_planeWidgetZ = GetPlaneWidget('z', 0, 0, 1, picker);
+
+	_planeWidgetY->SetLookupTable(_planeWidgetX->GetLookupTable());
+	_planeWidgetZ->SetLookupTable(_planeWidgetX->GetLookupTable());
+
+	picker->UnRegister(NULL);
+
+	// ColorWindow ColorLevel Callback
+	vtkWindowLevelImagePlaneWidgetCallback *wlipwc	= vtkWindowLevelImagePlaneWidgetCallback::New();
+	wlipwc->_wxvtkmpr3Dview							= this;
+	_planeWidgetX->AddObserver(vtkCommand::WindowLevelEvent, wlipwc );
+	_planeWidgetY->AddObserver(vtkCommand::WindowLevelEvent, wlipwc );
+	_planeWidgetZ->AddObserver(vtkCommand::WindowLevelEvent, wlipwc );
 	
-	_contourMapper->SetScalarRange( imageData->GetScalarRange() );
-    
-	
+	// Cursor Callback
+	_planeWidgetX->AddObserver(vtkCommand::InteractionEvent, wlipwc );
+	_planeWidgetY->AddObserver(vtkCommand::InteractionEvent, wlipwc );
+	_planeWidgetZ->AddObserver(vtkCommand::InteractionEvent, wlipwc );	
+}
 
-			
+
+//-------------------------------------------------------------------
+void wxVtkMPR3DView::SetImage()
+{
+	vtkImageData *imageData					= GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData();
+	// Orthogonal planes B&W
+	_planeWidgetX->SetInput( imageData );
+	_planeWidgetY->SetInput( imageData );
+	_planeWidgetZ->SetInput( imageData );
+	// -- Plane widget
+	_probe->SetSource( imageData );
+	_vtkmpr3Ddataviewer->SetImage();
 }
+
 //-------------------------------------------------------------------
-wxVtkMPR3DView::~wxVtkMPR3DView()
+vtkImagePlaneWidget* wxVtkMPR3DView::GetPlaneWidget(unsigned char activationkey, double r, double g, double b, vtkCellPicker* picker)
 {
-	printf("EED wxVtkMPR3DView::~wxVtkMPR3DView\n");
-	VisiblePointWidget( false );
-	VisiblePlaneWidget( false );
-	
-	if(_pointWidget!=NULL){
-		_pointWidget->Delete();
+	vtkProperty* prop1					= 0;
+	vtkImagePlaneWidget* planeWidget	= 0;
+	planeWidget							= vtkImagePlaneWidget::New();
+	planeWidget->SetInteractor( GetWxvtk3Dbaseview()->GetWxVTKRenderWindowInteractor() );
+	planeWidget->EnabledOff();
+	vtkImageData *image					= GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData();
+	int xMin = 0, xMax = 0, yMin = 0, yMax = 0, zMin = 0, zMax = 0;
+	double xSpacing = 0, ySpacing = 0, zSpacing = 0;
+	if(image)
+	{
+		planeWidget->SetInput( image );
+		image->GetExtent(xMin, xMax, yMin, yMax, zMin, zMax);
+		image->GetSpacing(xSpacing, ySpacing, zSpacing);
 	}
-	if(_planeWidget!=NULL){
-		_planeWidget->Delete();
+
+	if (activationkey=='x')
+	{
+		planeWidget->SetPlaneOrientationToXAxes();
+		planeWidget->SetSlicePosition((xMax+xMin)/2.*xSpacing);
 	}
 
-	if(_vtkplane!=NULL){
-		_vtkplane->Delete();		
-		_probe->Delete();
-		_contourMapper->Delete();			
-		_contourPlaneActor->Delete();		
+	if (activationkey=='y')
+	{
+		planeWidget->SetPlaneOrientationToYAxes();
+		planeWidget->SetSlicePosition((yMax+yMin)/2.*ySpacing);
+	}
+
+	if (activationkey=='z')
+	{
+		planeWidget->SetPlaneOrientationToZAxes();
+		planeWidget->SetSlicePosition((zMax+zMin)/2.*zSpacing);
 	}
 
+	planeWidget->DisplayTextOn();
+	planeWidget->SetPicker(picker);
+	planeWidget->SetKeyPressActivationValue(activationkey);
+	prop1 = planeWidget->GetPlaneProperty();
+	prop1->SetColor(r, g, b);
+	return planeWidget;
+}
+
+//-------------------------------------------------------------------
+void wxVtkMPR3DView::SetFreePlanesOrtho()
+{
+	_planeWidgetX->SetPlaneOrientationToXAxes();
+	_planeWidgetY->SetPlaneOrientationToYAxes();
+	_planeWidgetZ->SetPlaneOrientationToZAxes();
+	vtkInteractorStyleBaseView *isbv = (vtkInteractorStyleBaseView*)GetWxvtk3Dbaseview()->GetInteractorStyleBaseView();
+	isbv->SetParent_refresh_waiting();
+	isbv->EvaluateToRefresh();
+}
+
+//-------------------------------------------------------------------
+void wxVtkMPR3DView::FreePlaneVisible(int plane, bool ok)
+{
+	if(_planeWidgetX)
+	{
+              	if (ok==true)
+              	{
+//EED 21/07/2013
+//			1. Make backup of actors
+// 			2. Putting out all the actors
+//			3. Add vtkImagageActor in first place
+//			4. Putting again the originals actors
+//			5. Refresh
+
+			// step 1, step 2
+        		vtkActorCollection *tmpCollection    = vtkActorCollection::New();
+	 		vtkActorCollection *collectionActors =  _wxvtk3Dbaseview->GetRenderer()->GetActors();
+        		int i,size=collectionActors->GetNumberOfItems();
+			for (i=0; i<size; i++)
+			{
+				tmpCollection->AddItem ( (vtkActor*)collectionActors->GetItemAsObject(0)  );           		
+			        _wxvtk3Dbaseview->GetRenderer()->RemoveActor( (vtkProp*)tmpCollection->GetItemAsObject(i) );
+			} // for
+
+			// step 3
+ 	      		if ((plane==1) && (ok==true)) _planeWidgetX->EnabledOn();
+	      		if ((plane==2) && (ok==true)) _planeWidgetY->EnabledOn();
+	      		if ((plane==3) && (ok==true)) _planeWidgetZ->EnabledOn();
+
+			// step 4
+			 for (i=0; i<size; i++)
+			{
+			        _wxvtk3Dbaseview->GetRenderer()->AddActor( (vtkProp*)tmpCollection->GetItemAsObject(i) );
+			} // for
+			tmpCollection->Delete();
+
+			// step 5 Refresh
+			_wxvtk3Dbaseview->GetRenderer()->GetRenderWindow()->Render();
+              	} // if ok
+		if ((plane==1) && (ok==false)) _planeWidgetX->EnabledOff();
+		if ((plane==2) && (ok==false)) _planeWidgetY->EnabledOff();
+		if ((plane==3) && (ok==false)) _planeWidgetZ->EnabledOff();
+	} // if  _planeWidgetX
+}
+
+//-------------------------------------------------------------------
+void wxVtkMPR3DView::FreePlaneInteraction(bool ok)
+{
+	if(_planeWidgetX)
+	{
+		if (ok==true)
+		{
+			_planeWidgetX->InteractionOn();
+			_planeWidgetY->InteractionOn();
+			_planeWidgetZ->InteractionOn();
+		} else {
+			_planeWidgetX->InteractionOff();
+			_planeWidgetY->InteractionOff();
+			_planeWidgetZ->InteractionOff();
+		}
+	} // if _planeWidgetX
+}
+
+
+
+//-------------------------------------------------------------------
+void wxVtkMPR3DView::setColorTransferFunction(vtkColorTransferFunction* colortable)
+{
+	GetVtkMPR3DDataViewer ()->setColorTransferFunction(colortable);
 }
 
 //-------------------------------------------------------------------
 void wxVtkMPR3DView::Refresh()   // virtual
 {
-	_vtkmpr3Ddataviewer	->	Refresh();
+	_vtkmpr3Ddataviewer->Refresh();
 	if (_wxvtkmpr3DviewCntrlPanel!=NULL)
 	{
 		_wxvtkmpr3DviewCntrlPanel->Refresh();
 	}
-
 }
+
 //-------------------------------------------------------------------
 void wxVtkMPR3DView::RefreshView()   // virtual
 {
 	double spc[3];
-	this->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData()->GetSpacing(spc);
+	vtkImageData *image = GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData();
+
+	if(image)
+	{
+		image->GetSpacing(spc);
+	}else{
+		spc[0] = 0;
+		spc[1] = 0;
+		spc[2] = 0;
+	}
+	//this->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData()->GetSpacing(spc);
 
 	//double nx=1,ny=0,nz=0; // JPRx
 	Refresh();
-	double x=GetVtkMPR3DDataViewer ()->GetVtkMPRBaseData()->GetX();
-	double y=GetVtkMPR3DDataViewer ()->GetVtkMPRBaseData()->GetY();
-	double z=GetVtkMPR3DDataViewer ()->GetVtkMPRBaseData()->GetZ();
+	double x=GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetX();
+	double y=GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetY();
+	double z=GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetZ();
 
-	x=x*spc[0];
-	y=y*spc[1];
-	z=z*spc[2];
-
-	_pointWidget->SetPosition( x,y,z ); 
+	SetColorWindow( _vtkmpr3Ddataviewer->GetVtkMPRBaseData()->GetColorWindow() );
+	SetColorLevel( _vtkmpr3Ddataviewer->GetVtkMPRBaseData()->GetColorLevel() );
 
 
-	double in[3]; 
 	double normal[3];
-	vtkTransform *transform = GetVtkMPR3DDataViewer ()->GetVtkMPRBaseData()->GetTransformOrientation();
-	in[0]=1;		in[1]=0;			in[2]=0;
-	transform->TransformPoint(in,normal);
-	_planeWidget->SetNormal( normal[0],normal[1],normal[2] );
-
+	// FreePlaneX
+	_planeWidgetX->GetNormal(normal);
+	if ((normal[0]==1)&&(normal[1]==0)&&(normal[2]==0))
+	{
+		_planeWidgetX->SetSlicePosition( x*spc[0] );
+//EED 2016/02/19
+//		_planeWidgetX->GetReslice()->SetInterpolate( _vtkmpr3Ddataviewer->GetVtkMPRBaseData()->GetInterpolate() );
+		_planeWidgetX->GetTexture()->SetInterpolate( _vtkmpr3Ddataviewer->GetVtkMPRBaseData()->GetInterpolate() );
+	}
+	// FreePlaneY
+	_planeWidgetY->GetNormal(normal);
+	if ((normal[0]==0)&&(normal[1]==1)&&(normal[2]==0))
+	{
+		_planeWidgetY->SetSlicePosition( y*spc[1] );
+//EED 2016/02/19
+		_planeWidgetY->GetReslice()->SetInterpolate( _vtkmpr3Ddataviewer->GetVtkMPRBaseData()->GetInterpolate() );
+		_planeWidgetY->GetTexture()->SetInterpolate( _vtkmpr3Ddataviewer->GetVtkMPRBaseData()->GetInterpolate() );
+	}
+	// FreePlaneZ
+	_planeWidgetZ->GetNormal(normal);
+	if ((normal[0]==0)&&(normal[1]==0)&&(normal[2]==1))
+	{
+		_planeWidgetZ->SetSlicePosition( z*spc[2] );
+//EED 2016/02/19
+		_planeWidgetZ->GetReslice()->SetInterpolate( _vtkmpr3Ddataviewer->GetVtkMPRBaseData()->GetInterpolate() );
+		_planeWidgetZ->GetTexture()->SetInterpolate( _vtkmpr3Ddataviewer->GetVtkMPRBaseData()->GetInterpolate() );
+	}
+		
+	x=x*spc[0];
+	y=y*spc[1];
+	z=z*spc[2];
 
-	_planeWidget->SetCenter( x,y,z );
-	_planeWidget->UpdatePlacement();
-	_planeWidget->GetPolyData(_vtkplane);
+	if(_pointWidget){
+		_pointWidget->SetPosition( x,y,z );
+		double in[3];
+		double normal[3];
+		vtkTransform *transform = GetVtkMPR3DDataViewer ()->GetVtkMPRBaseData()->GetTransformOrientation();
+		in[0]=1;		in[1]=0;			in[2]=0;
+		transform->TransformPoint(in,normal);
+		_planeWidget->SetNormal( normal[0],normal[1],normal[2] );
+		_planeWidget->SetCenter( x,y,z );
+		_planeWidget->UpdatePlacement();
+		_planeWidget->GetPolyData(_vtkplane);
+	}
 }
+
 //-------------------------------------------------------------------
-wxPanel* wxVtkMPR3DView::CreateControlPanel(wxWindow *parent)
+wxPanel* wxVtkMPR3DView::CreateControlPanel(wxWindow *parent, bool align)
 {
-	_wxvtkmpr3DviewCntrlPanel = new wxVtkMPR3DViewCntrlPanel(parent,this);
+	_wxvtkmpr3DviewCntrlPanel = new wxVtkMPR3DViewCntrlPanel(parent,this, align);
 	return _wxvtkmpr3DviewCntrlPanel;
 }
 //-------------------------------------------------------------------
-vtkMPR3DDataViewer* wxVtkMPR3DView::GetVtkMPR3DDataViewer() 
+vtkMPR3DDataViewer* wxVtkMPR3DView::GetVtkMPR3DDataViewer()
 {
-	return _vtkmpr3Ddataviewer; 
+	return _vtkmpr3Ddataviewer;
 }
 //-------------------------------------------------------------------
 wxVtk3DBaseView* wxVtkMPR3DView::GetWxvtk3Dbaseview() throw (char*)
 {
-	if(_wxvtk3Dbaseview==NULL){
+	if(_wxvtk3Dbaseview==NULL)
+	{
 		throw "wxVtk3DBaseView* wxVtkMPR3DView::GetWxvtk3Dbaseview() _wxvtk3Dbaseview=NULL";
 	}
 	return _wxvtk3Dbaseview;
@@ -267,11 +590,54 @@ wxVtk3DBaseView* wxVtkMPR3DView::GetWxvtk3Dbaseview() throw (char*)
 
 //-------------------------------------------------------------------
 void wxVtkMPR3DView::VisibleImageActor(int idPosition, bool visible){
+	/*if(visible == false)
+		printf("CPR: wxVtkMPR3DView::VisibleImageActor-> visibleActor == false \n");
+	else
+	{
+		printf("CPR: wxVtkMPR3DView::VisibleImageActor-> visibleActor == true\n");
+		printf("CPR: wxVtkMPR3DView::VisibleImageActor-> idPosition == %d\n", idPosition);
+	}
+	*/
 	if (visible!=_vtkmpr3Ddataviewer->GetVisiblePosition(idPosition)){
 		if (visible==false){
-			_wxvtk3Dbaseview->GetRenderer()->RemoveActor( _vtkmpr3Ddataviewer->GetImageActor(idPosition)  );	
+			_wxvtk3Dbaseview->GetRenderer()->RemoveActor( _vtkmpr3Ddataviewer->GetImageActor(idPosition)  );
 		} else {
-			_wxvtk3Dbaseview->GetRenderer()->AddActor( _vtkmpr3Ddataviewer->GetImageActor(idPosition)  );	
+//EED 21/07/2013
+//			1. Make backup of actors
+// 			2. Putting out all the actors
+//			3. Add vtkImagageActor in first place
+//			4. Putting again the originals actors
+//			5. Refresh
+
+//EED 4 nov 2015  (For Estelle)
+// This works in ViewerNV with nTymeView=6
+// Because the others modules (5) generate also more actors with an OutLineFilter -> Actor
+//showOutlineActor(false);
+
+			// step 1, step 2
+			vtkActorCollection *tmpCollection    = vtkActorCollection::New();
+	 		vtkActorCollection *collectionActors =  _wxvtk3Dbaseview->GetRenderer()->GetActors();
+        		int i,size=collectionActors->GetNumberOfItems();
+			for (i=0; i<size; i++)
+			{
+				printf(" wxVtkMPR3DView::VisibleImageActor %d\n" , i);
+				tmpCollection->AddItem ( (vtkActor*)collectionActors->GetItemAsObject(0)  );           		
+			        _wxvtk3Dbaseview->GetRenderer()->RemoveActor( (vtkProp*)tmpCollection->GetItemAsObject(i) );
+			} // for
+
+			// step 3
+			_wxvtk3Dbaseview->GetRenderer()->AddActor( _vtkmpr3Ddataviewer->GetImageActor(idPosition)  );
+
+			// step 4
+			 for (i=0; i<size; i++)
+			{
+			        _wxvtk3Dbaseview->GetRenderer()->AddActor( (vtkProp*)tmpCollection->GetItemAsObject(i) );
+			} // for
+			tmpCollection->Delete();
+
+			// step 5 Refresh
+			_wxvtk3Dbaseview->GetRenderer()->GetRenderWindow()->Render();
+
 		}
 		_vtkmpr3Ddataviewer->SetVisiblePosition(idPosition,visible);
 	}
@@ -281,28 +647,48 @@ void wxVtkMPR3DView::VisibleImageActor(int idPosition, bool visible){
 //-------------------------------------------------------------------
 void wxVtkMPR3DView::VisiblePointWidget( bool visible )
 {
-	if (visible==true)
-	{
-		_pointWidget->On();
-	} else {
-		_pointWidget->Off();
+	if(_pointWidget){
+		if (visible==true)
+		{
+			_pointWidget->On();
+		} else {
+			_pointWidget->Off();
+		}
 	}
 }
 
 //-------------------------------------------------------------------
 void wxVtkMPR3DView::VisiblePlaneWidget( bool visible )
 {
-	if (visible==true)
-	{
-		_planeWidget->On();
-		_contourPlaneActor->VisibilityOn();
+	if(_planeWidget){
+		if (visible==true)
+		{
+			_planeWidget->On();
+			_contourPlaneActor->VisibilityOn();
 
-	} else {
-		_planeWidget->Off();
-		_contourPlaneActor->VisibilityOff();
+		} else {
+			_planeWidget->Off();
+			_contourPlaneActor->VisibilityOff();
+		}
 	}
 }
 
+//CPR: Method added 30 Nov 2009
+void wxVtkMPR3DView::showOutlineActor(bool value)
+{
+
+	vtkActor* _outlineActor = _vtkmpr3Ddataviewer->GetOutlineActor();
+	if(value == true)
+	{
+		_wxvtk3Dbaseview->GetRenderer()->AddActor( _outlineActor );
+	}
+	else
+	{
+		_wxvtk3Dbaseview->GetRenderer()->RemoveActor( _outlineActor );
+	}
+
+}
+
 //-------------------------------------------------------------------
 void wxVtkMPR3DView::SetVtkMPR3DDataViewer(vtkMPR3DDataViewer *vtkmpr3Ddataviewer)
 {
@@ -312,11 +698,22 @@ void wxVtkMPR3DView::SetVtkMPR3DDataViewer(vtkMPR3DDataViewer *vtkmpr3Ddataviewe
 //-------------------------------------------------------------------
 void wxVtkMPR3DView::InitOrientationPointWidget()
 {
-	this->_planeWidget->SetNormal(1,0,0);   	
+	this->_planeWidget->SetNormal(1,0,0);
 }
 
+//-------------------------------------------------------------------
+void wxVtkMPR3DView::SetColorWindow(double colorWindow)
+{
+	double colorLevel = _planeWidgetX->GetLevel();
+	_planeWidgetX->SetWindowLevel(colorWindow,colorLevel);
+}
 
-
+//-------------------------------------------------------------------
+void wxVtkMPR3DView::SetColorLevel(double colorLevel)
+{
+	double colorWindow = _planeWidgetX->GetWindow();
+	_planeWidgetX->SetWindowLevel(colorWindow,colorLevel);
+}
 
 //-------------------------------------------------------------------
 // EED 25 Janvier 2007 testLoic
@@ -364,7 +761,7 @@ void wxVtkMPR3DView::TestLoic1()
 	stPoints -> Update();
 
 	vtkImageChangeInformation  *change = vtkImageChangeInformation ::New();
-	change	  -> SetInput( stPoints );  
+	change	  -> SetInput( stPoints );
 	change	  -> Update();    //important
 
 	double _range[2];
@@ -377,7 +774,7 @@ void wxVtkMPR3DView::TestLoic1()
    _bwlookup->SetValueRange( 0 , 1 );
    _bwlookup->SetAlphaRange( 1 , 1 );
    _bwlookup->Build( );
- 
+
    vtkLookupTable * _collookup = vtkLookupTable::New( );
    _collookup->SetNumberOfColors( 256 );
    _collookup->SetTableRange( 0 , 255 );
@@ -420,7 +817,7 @@ void wxVtkMPR3DView::TestLoic1()
       renderer->SetActiveCamera(camera);
       renderer->ResetCamera();
       renderer->SetBackground(1,1,1);
-  
+
   renWin->SetSize(300,300);
 
   // interact with data
@@ -447,7 +844,6 @@ void wxVtkMPR3DView::TestLoic2()
         vtkPlane *slicePlane = vtkPlane::New();
 				this->_planeWidget->GetPlane( slicePlane );
 
-                
         vtkCutter* sliceCutter = vtkCutter::New();
 		vtkImageData *vtkimagedata =  this->GetVtkMPR3DDataViewer()->GetVtkMPRBaseData()->GetImageData();
                 sliceCutter->SetInput( vtkimagedata );
@@ -456,7 +852,6 @@ void wxVtkMPR3DView::TestLoic2()
 //EED
 //        vtkLookupTable *lut = BuildHueWeightBaseMap();
 
-
         vtkPolyDataMapper *slice = vtkPolyDataMapper::New();
                 slice->SetInput( sliceCutter->GetOutput() );
         double range[2];
@@ -476,12 +871,9 @@ void wxVtkMPR3DView::TestLoic2()
                 contourMapper->SetScalarRange( range );
   //              contourMapper->SetLookupTable( lut );
 
-                
         vtkActor *contourActor = vtkActor::New();
                 contourActor->SetMapper( contourMapper );
 
-
-
  // The usual rendering stuff.
   vtkCamera *camera = vtkCamera::New();
       camera->SetPosition(1,1,1);
@@ -499,13 +891,33 @@ void wxVtkMPR3DView::TestLoic2()
       renderer->SetActiveCamera(camera);
       renderer->ResetCamera();
       renderer->SetBackground(1,1,1);
-  
+
   renWin->SetSize(300,300);
 
   // interact with data
   renWin->Render();
   iren->Start();
+}
 
-
+//-------------------------------------------------------------------
+void wxVtkMPR3DView::SetBackGroundType(int type)
+{
+	if (_wxvtk3Dbaseview!=NULL)
+	{
+		if (type==0)
+		{
+			_wxvtk3Dbaseview->GetRenderer()->GradientBackgroundOff();
+			_wxvtk3Dbaseview->GetRenderer()->SetBackground( 0 , 0 , 0 );
+		} // if type
+		if (type==1)
+		{
+			_wxvtk3Dbaseview->GetRenderer()->GradientBackgroundOn();
+			_wxvtk3Dbaseview->GetRenderer()->SetBackground( 0.33 , 0.33 , 0.33 );
+			_wxvtk3Dbaseview->GetRenderer()->SetBackground2( 0.66 , 0.66 , 0.66 );
+		} // if type
+		_wxvtk3Dbaseview->Refresh();
+	}// if _wxvtk3Dbaseview
+	
 }
 
+