]> Creatis software - cpPlugins.git/blobdiff - lib/cpExtensions/Visualization/MPRActors.cxx
Merge branch 'master' of ssh://git.creatis.insa-lyon.fr/cpPlugins
[cpPlugins.git] / lib / cpExtensions / Visualization / MPRActors.cxx
index 90fdb7d83c7f68842d4d356191b61bdeedbe03c0..72bf719729b46078d34948f7aacf5477e859a71c 100644 (file)
@@ -1,12 +1,9 @@
 #include <cpExtensions/Visualization/MPRActors.h>
 
-#include <vtkAlgorithmOutput.h>
 #include <vtkImageData.h>
-#include <vtkInformation.h>
-#include <vtkOutlineSource.h>
-#include <vtkProperty.h>
 #include <vtkRenderer.h>
-#include <vtkWindowLevelLookupTable.h>
+#include <vtkRendererCollection.h>
+#include <vtkRenderWindow.h>
 
 // -------------------------------------------------------------------------
 cpExtensions::Visualization::MPRActors*
@@ -29,21 +26,26 @@ GetSliceActors( const int& i ) const
 
 // -------------------------------------------------------------------------
 int cpExtensions::Visualization::MPRActors::
-AddInputConnection( vtkAlgorithmOutput* aout )
+AddInputConnection(
+  vtkAlgorithmOutput* aout, ImageSliceActors::LUTType lut
+  )
 {
-  int N = this->ImageMaps.size( );
+  int N = this->Slices[ 0 ][ 0 ]->GetNumberOfImageActors( );
   if( N == 0 )
   {
-    this->ImageMaps.push_back(
-      vtkSmartPointer< vtkImageMapToColors >::New( )
-      );
-    this->ImageMaps[ 0 ]->SetInputConnection( aout );
-    this->SetLookupTableToWindowLevel( 0 );
-    this->_Update( 0 );
-    return( 0 );
+    this->Slices[ 0 ][ 0 ]->AddInputConnection( aout, 0, lut );
+    vtkImageMapToColors* imap = this->Slices[ 0 ][ 0 ]->GetImageMap( 0 );
+    for( unsigned int i = 0; i < 2; ++i )
+      for( unsigned int j = 0; j < 3; ++j )
+        if( i != 0 || j != 0 )
+          this->Slices[ i ][ j ]->AddInputConnection(
+            ( ( imap != NULL )? imap->GetOutputPort( ): aout ), j,
+            ImageSliceActors::LUTType_None
+            );
   }
   else
   {
+    /*
     // Check if the image share the same space
     vtkImageData* ref_image = this->_Image( 0 );
     vtkImageData* new_image =
@@ -64,20 +66,50 @@ AddInputConnection( vtkAlgorithmOutput* aout )
         vtkSmartPointer< vtkImageMapToColors >::New( )
         );
       this->ImageMaps[ N ]->SetInputConnection( aout );
-      this->SetLookupTableToWindowLevel( N );
+      this->SetLookupTableToColor(
+        N, double( 1 ), double( 0 ), double( 0 )
+        );
       this->_Update( N );
       return( N );
     }
     else
-      return( -1 );
+    N = -1;
+  */
 
   } // fi
+  return( N );
 }
 
 // -------------------------------------------------------------------------
 int cpExtensions::Visualization::MPRActors::
-AddInputData( vtkImageData* image )
+AddInputData( vtkImageData* image, ImageSliceActors::LUTType lut )
 {
+  int N = this->Slices[ 0 ][ 0 ]->GetNumberOfImageActors( );
+  if( N == 0 )
+  {
+    this->Slices[ 0 ][ 0 ]->AddInputData( image, 0, lut );
+    vtkImageMapToColors* imap = this->Slices[ 0 ][ 0 ]->GetImageMap( 0 );
+    for( unsigned int i = 0; i < 2; ++i )
+      for( unsigned int j = 0; j < 3; ++j )
+        if( i != 0 || j != 0 )
+        {
+          if( imap != NULL )
+            this->Slices[ i ][ j ]->AddInputConnection(
+              imap->GetOutputPort( ), j,
+              ImageSliceActors::LUTType_None
+              );
+          else
+            this->Slices[ i ][ j ]->AddInputData(
+              image, j,
+              ImageSliceActors::LUTType_None
+              );
+        } // fi
+  }
+  else
+  {
+  } // fi
+  return( N );
+  /*
   int N = this->ImageMaps.size( );
   if( N == 0 )
   {
@@ -87,6 +119,7 @@ AddInputData( vtkImageData* image )
     this->ImageMaps[ 0 ]->SetInputData( image );
     this->SetLookupTableToWindowLevel( 0 );
     this->_Update( 0 );
+    this->ResetWindowLevel( 0 );
     return( 0 );
   }
   else
@@ -107,7 +140,9 @@ AddInputData( vtkImageData* image )
         vtkSmartPointer< vtkImageMapToColors >::New( )
         );
       this->ImageMaps[ N ]->SetInputData( image );
-      this->SetLookupTableToWindowLevel( N );
+      this->SetLookupTableToColor(
+        N, double( 1 ), double( 0 ), double( 0 )
+        );
       this->_Update( N );
       return( N );
     }
@@ -115,43 +150,49 @@ AddInputData( vtkImageData* image )
       return( -1 );
 
   } // fi
+  */
 }
 
 // -------------------------------------------------------------------------
 void cpExtensions::Visualization::MPRActors::
-PushDataInto( vtkRenderer* x, vtkRenderer* y, vtkRenderer* z, vtkRenderer* w )
+PushActorsInto(
+  vtkRenderWindow* x, vtkRenderWindow* y, vtkRenderWindow* z,
+  vtkRenderWindow* w
+  )
 {
-  vtkRenderer* rends[] = { x, y, z };
+  this->Slices[ 0 ][ 0 ]->PushActorsInto( x, true );
+  this->Slices[ 0 ][ 1 ]->PushActorsInto( y, true );
+  this->Slices[ 0 ][ 2 ]->PushActorsInto( z, true );
+
+  vtkRenderer* wren =
+    ( w != NULL )? w->GetRenderers( )->GetFirstRenderer( ): NULL;
+  vtkRenderer* rends[ ] =
+    {
+      ( x != NULL )? x->GetRenderers( )->GetFirstRenderer( ): NULL,
+      ( y != NULL )? y->GetRenderers( )->GetFirstRenderer( ): NULL,
+      ( z != NULL )? z->GetRenderers( )->GetFirstRenderer( ): NULL
+    };
   for( int i = 0; i < 3; ++i )
   {
     if( rends[ i ] != NULL )
-    {
-      for(
-        unsigned int k = 0;
-        k < this->Slices[ 0 ][ i ]->GetNumberOfImageActors( );
-        ++k
-        )
-        rends[ i ]->AddActor( this->Slices[ 0 ][ i ]->GetImageActor( k ) );
-      rends[ i ]->AddActor( this->Slices[ 0 ][ i ]->GetTextActor( ) );
       for( int j = 0; j < 3; ++j )
-        rends[ i ]->AddActor( this->Slices[ 0 ][ j ]->GetPlaneActor( ) );
-
-    } // fi
-    if( w != NULL )
+        if( i != j )
+          rends[ i ]->AddActor( this->Slices[ 0 ][ j ]->GetPlaneActor( ) );
+    if( wren != NULL )
     {
       for(
         unsigned int k = 0;
         k < this->Slices[ 1 ][ i ]->GetNumberOfImageActors( );
         ++k
         )
-        w->AddActor( this->Slices[ 1 ][ i ]->GetImageActor( k ) );
-      w->AddActor( this->Slices[ 1 ][ i ]->GetPlaneActor( ) );
+        wren->AddActor( this->Slices[ 1 ][ i ]->GetImageActor( k ) );
+      wren->AddActor( this->Slices[ 1 ][ i ]->GetPlaneActor( ) );
 
     } // fi
 
   } // rof
-  if( w != NULL )
-    w->AddActor( this->ImageOutlineActor );
+  if( wren != NULL )
+    wren->AddActor( this->ImageOutlineActor );
 }
 
 // -------------------------------------------------------------------------
@@ -163,6 +204,7 @@ PopDataFrom(
   vtkRenderer* w
   )
 {
+  /*
   vtkRenderer* rends[] = { x, y, z };
   for( int i = 0; i < 3; ++i )
   {
@@ -194,12 +236,25 @@ PopDataFrom(
   } // rof
   if( w != NULL )
     w->RemoveActor( this->ImageOutlineActor );
+  */
 }
 
+// -------------------------------------------------------------------------
+  /*
+void cpExtensions::Visualization::MPRActors::
+LinkInteractors( )
+{
+  this->_Update( 0 );
+  this->_Update( 1 );
+  this->_Update( 2 );
+}
+  */
+
 // -------------------------------------------------------------------------
 void cpExtensions::Visualization::MPRActors::
 SetLookupTable( unsigned int i, vtkScalarsToColors* lut )
 {
+  /*
   if( i < this->ImageMaps.size( ) )
   {
     this->ImageMaps[ i ]->SetLookupTable( lut );
@@ -207,22 +262,26 @@ SetLookupTable( unsigned int i, vtkScalarsToColors* lut )
     this->Modified( );
 
   } // fi
+  */
 }
 
 // -------------------------------------------------------------------------
 vtkScalarsToColors* cpExtensions::Visualization::MPRActors::
 GetLookupTable( unsigned int i ) const
 {
+  /*
   if( i < this->ImageMaps.size( ) )
     return( this->ImageMaps[ i ]->GetLookupTable( ) );
   else
     return( NULL );
+    */
 }
 
 // -------------------------------------------------------------------------
 void cpExtensions::Visualization::MPRActors::
 SetLookupTableToWindowLevel( unsigned int i )
 {
+  /*
   // Check if the input has been configured
   vtkImageData* image = this->_Image( i );
   if( image == NULL )
@@ -238,6 +297,7 @@ SetLookupTableToWindowLevel( unsigned int i )
   lut->Build( );
 
   this->SetLookupTable( i, lut );
+  */
 }
 
 // -------------------------------------------------------------------------
@@ -251,6 +311,7 @@ GetMinWindow( unsigned int i ) const
 double cpExtensions::Visualization::MPRActors::
 GetMaxWindow( unsigned int i ) const
 {
+  /*
   // Check if the input has been configured
   vtkImageData* image = this->_Image( i );
   if( image == NULL )
@@ -259,12 +320,14 @@ GetMaxWindow( unsigned int i ) const
   double r[ 2 ];
   image->GetScalarRange( r );
   return( r[ 1 ] - r[ 0 ] );
+  */
 }
 
 // -------------------------------------------------------------------------
 double cpExtensions::Visualization::MPRActors::
 GetMinLevel( unsigned int i ) const
 {
+  /*
   // Check if the input has been configured
   vtkImageData* image = this->_Image( i );
   if( image == NULL )
@@ -273,12 +336,14 @@ GetMinLevel( unsigned int i ) const
   double r[ 2 ];
   image->GetScalarRange( r );
   return( r[ 0 ] );
+  */
 }
 
 // -------------------------------------------------------------------------
 double cpExtensions::Visualization::MPRActors::
 GetMaxLevel( unsigned int i ) const
 {
+  /*
   // Check if the input has been configured
   vtkImageData* image = this->_Image( i );
   if( image == NULL )
@@ -287,36 +352,42 @@ GetMaxLevel( unsigned int i ) const
   double r[ 2 ];
   image->GetScalarRange( r );
   return( r[ 1 ] );
+  */
 }
 
 // -------------------------------------------------------------------------
 double cpExtensions::Visualization::MPRActors::
 GetWindow( unsigned int i ) const
 {
+  /*
   vtkWindowLevelLookupTable* lut =
     dynamic_cast< vtkWindowLevelLookupTable* >( this->GetLookupTable( i ) );
   if( lut != NULL )
     return( lut->GetWindow( ) );
   else
     return( double( 0 ) );
+  */
 }
 
 // -------------------------------------------------------------------------
 double cpExtensions::Visualization::MPRActors::
 GetLevel( unsigned int i ) const
 {
+  /*
   vtkWindowLevelLookupTable* lut =
     dynamic_cast< vtkWindowLevelLookupTable* >( this->GetLookupTable( i ) );
   if( lut != NULL )
     return( lut->GetLevel( ) );
   else
     return( double( 0 ) );
+  */
 }
 
 // -------------------------------------------------------------------------
 void cpExtensions::Visualization::MPRActors::
 SetWindow( unsigned int i, const double& w )
 {
+  /*
   vtkWindowLevelLookupTable* lut =
     dynamic_cast< vtkWindowLevelLookupTable* >( this->GetLookupTable( i ) );
   if( lut != NULL )
@@ -327,28 +398,55 @@ SetWindow( unsigned int i, const double& w )
     this->Modified( );
 
   } // fi
+  */
 }
 
 // -------------------------------------------------------------------------
 void cpExtensions::Visualization::MPRActors::
 SetLevel( unsigned int i, const double& l )
 {
+  /*
+  vtkWindowLevelLookupTable* lut =
+    dynamic_cast< vtkWindowLevelLookupTable* >( this->GetLookupTable( i ) );
+  if( lut != NULL )
+  {
+    lut->SetLevel( l );
+    lut->Build( );
+    this->ImageMaps[ i ]->Modified( );
+    this->Modified( );
+
+  } // fi
+  */
+}
+
+// -------------------------------------------------------------------------
+void cpExtensions::Visualization::MPRActors::
+SetWindowLevel( unsigned int i, const double& w, const double& l )
+{
+  /*
   vtkWindowLevelLookupTable* lut =
     dynamic_cast< vtkWindowLevelLookupTable* >( this->GetLookupTable( i ) );
   if( lut != NULL )
   {
+    lut->SetWindow( w );
     lut->SetLevel( l );
     lut->Build( );
     this->ImageMaps[ i ]->Modified( );
+
+    for( unsigned int j = 0; j < 3; ++j )
+      this->Slices[ 0 ][ j ]->UpdateText( w, l );
+
     this->Modified( );
 
   } // fi
+  */
 }
 
 // -------------------------------------------------------------------------
 void cpExtensions::Visualization::MPRActors::
 ResetWindowLevel( unsigned int i )
 {
+  /*
   vtkImageData* image = this->_Image( i );
   vtkWindowLevelLookupTable* lut =
     dynamic_cast< vtkWindowLevelLookupTable* >( this->GetLookupTable( i ) );
@@ -357,11 +455,72 @@ ResetWindowLevel( unsigned int i )
     double r[ 2 ];
     image->GetScalarRange( r );
     lut->SetTableRange( r );
+    lut->SetWindow( r[ 1 ] - r[ 0 ] );
+    lut->SetLevel( ( r[ 1 ] + r[ 0 ] ) / double( 2 ) );
     lut->Build( );
     this->ImageMaps[ i ]->Modified( );
     this->Modified( );
 
   } // fi
+  */
+}
+
+// -------------------------------------------------------------------------
+void cpExtensions::Visualization::MPRActors::
+SetLookupTableToColor(
+  unsigned int i, const double& r, const double& g, const double& b
+  )
+{
+  /*
+  static const double _0 = double( 0 );
+  static const double _1 = double( 1 );
+  static const double _2 = double( 2 );
+  static const double _4 = double( 4 );
+  static const double _6 = double( 6 );
+  static const double _OPACITY = double( 0.6 );
+
+  // Check if the input has been configured
+  vtkImageData* image = this->_Image( i );
+  if( image == NULL )
+    return;
+
+  double range[ 2 ];
+  image->GetScalarRange( range );
+
+  // Get HSV from display color
+  double cmax = ( r > g )? r: g; cmax = ( b > cmax )? b: cmax;
+  double cmin = ( r < g )? r: g; cmin = ( b < cmin )? b: cmin;
+  double d = cmax - cmin;
+
+  double saturation = ( std::fabs( cmax ) > _0 )? d / cmax: _0;
+  double value = cmax;
+  double hue = _0;
+  if( d > _0 )
+  {
+    if( r == cmax )
+      hue = std::fmod( ( g - b ) / d, _6 );
+    else if( g == cmax )
+      hue = ( ( b - r ) / d ) + _2;
+    else if( b == cmax )
+      hue = ( ( r - g ) / d ) + _4;
+    hue *= _1 / _6;
+
+  } // fi
+
+  // Define new lookup table
+  vtkSmartPointer< vtkLookupTable > lut =
+    vtkSmartPointer< vtkLookupTable >::New( );
+  lut->SetScaleToLinear( );
+  lut->SetNanColor( _0, _0, _0, _0 );
+  lut->SetTableRange( range[ 0 ], range[ 1 ] );
+  lut->SetAlphaRange( _0, _OPACITY );
+  lut->SetHueRange( _0, hue );
+  lut->SetSaturationRange( _0, saturation );
+  lut->SetValueRange( _0, value );
+  lut->Build( );
+
+  this->SetLookupTable( i, lut );
+  */
 }
 
 // -------------------------------------------------------------------------
@@ -459,22 +618,27 @@ cpExtensions::Visualization::MPRActors::
 vtkImageData* cpExtensions::Visualization::MPRActors::
 _Image( unsigned int i ) const
 {
+  /*
   if( i < this->ImageMaps.size( ) )
   {
     vtkAlgorithm* algo = this->ImageMaps[ i ]->GetInputAlgorithm( );
     vtkInformation* info = algo->GetOutputInformation( 0 );
     return(
-      vtkImageData::SafeDownCast( info->Get( vtkDataObject::DATA_OBJECT( ) ) )
+      vtkImageData::SafeDownCast(
+        info->Get( vtkDataObject::DATA_OBJECT( ) )
+        )
       );
   }
   else
-    return( NULL );
+  */
+  return( NULL );
 }
 
 // -------------------------------------------------------------------------
 void cpExtensions::Visualization::MPRActors::
 _Update( unsigned int i )
 {
+  /*
   // Check if the input has been configured
   vtkImageData* image = this->_Image( i );
   if( image == NULL )
@@ -490,6 +654,16 @@ _Update( unsigned int i )
         );
       this->Slices[ j ][ k ]->UpdateText( );
 
+      // Add all of slice's props
+      this->Slices[ j ][ k ]->InitTraversal( );
+      vtkProp* prop = this->Slices[ j ][ k ]->GetNextProp( );
+      while( prop != NULL )
+      {
+        this->AddItem( prop );
+        prop = this->Slices[ j ][ k ]->GetNextProp( );
+
+      } // elihw
+
     } // rof
 
   } // rof
@@ -514,6 +688,7 @@ _Update( unsigned int i )
     this->AddItem( this->ImageOutlineActor );
 
   } // fi
+  */
 }
 
 // eof - $RCSfile$