#include <cpExtensions/Visualization/MPRActors.h>
+#include <cmath>
+
#include <vtkAlgorithmOutput.h>
#include <vtkImageData.h>
#include <vtkInformation.h>
+#include <vtkLookupTable.h>
+#include <vtkOutlineSource.h>
+#include <vtkProperty.h>
#include <vtkRenderer.h>
#include <vtkWindowLevelLookupTable.h>
this->ImageMaps[ 0 ]->SetInputConnection( aout );
this->SetLookupTableToWindowLevel( 0 );
this->_Update( 0 );
+ this->ResetWindowLevel( 0 );
return( 0 );
}
else
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 );
}
this->ImageMaps[ 0 ]->SetInputData( image );
this->SetLookupTableToWindowLevel( 0 );
this->_Update( 0 );
+ this->ResetWindowLevel( 0 );
return( 0 );
}
else
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 );
}
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 )
} // 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 )
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->SetAboveRangeColor( _0, _0, _0, _0 );
+ lut->SetBelowRangeColor( _0, _0, _0, _0 );
+ */
+ 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 );
+}
+
// -------------------------------------------------------------------------
int cpExtensions::Visualization::MPRActors::
GetSliceNumberMinValue( const int& axis ) const
void cpExtensions::Visualization::MPRActors::
ResetSlices( )
{
- // TODO
+ for( unsigned int i = 0; i < 2; ++i )
+ for( unsigned int j = 0; j < 3; ++j )
+ this->Slices[ i ][ j ]->SetSliceNumber(
+ this->Slices[ i ][ j ]->GetSliceNumberMinValue( )
+ );
}
// -------------------------------------------------------------------------
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
void cpExtensions::Visualization::MPRActors::
_Update( unsigned int i )
{
- /*
// Check if the input has been configured
- vtkImageData* image = this->_InputImage( );
+ vtkImageData* image = this->_Image( i );
if( image == NULL )
return;
- this->ImageMapToColors->Update( );
+ this->ImageMaps[ i ]->Update( );
for( int j = 0; j < 2; ++j )
{
- for( int i = 0; i < 3; ++i )
+ for( int k = 0; k < 3; ++k )
{
- this->Slices[ j ][ i ]->AddInputConnection(
- this->ImageMapToColors->GetOutputPort( ), i
+ this->Slices[ j ][ k ]->AddInputConnection(
+ this->ImageMaps[ i ]->GetOutputPort( ), k
);
- this->Slices[ j ][ i ]->UpdateText( );
+ 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
- // Create 3D outline
- double bb[ 6 ];
- image->GetBounds( bb );
+ if( i == 0 )
+ {
+ // Create 3D outline
+ double bb[ 6 ];
+ image->GetBounds( bb );
- vtkSmartPointer< vtkOutlineSource > img_ol =
- vtkSmartPointer< vtkOutlineSource >::New( );
- img_ol->SetBounds( bb );
+ vtkSmartPointer< vtkOutlineSource > img_ol =
+ vtkSmartPointer< vtkOutlineSource >::New( );
+ img_ol->SetBounds( bb );
- vtkSmartPointer< vtkPolyDataMapper > img_ol_mapper =
- vtkSmartPointer< vtkPolyDataMapper >::New( );
- img_ol_mapper->SetInputConnection( img_ol->GetOutputPort( ) );
- this->ImageOutlineActor->SetMapper( img_ol_mapper );
- this->ImageOutlineActor->GetProperty( )->SetColor( 1, 1, 1 );
- this->ImageOutlineActor->GetProperty( )->SetLineWidth( 1 );
+ vtkSmartPointer< vtkPolyDataMapper > img_ol_mapper =
+ vtkSmartPointer< vtkPolyDataMapper >::New( );
+ img_ol_mapper->SetInputConnection( img_ol->GetOutputPort( ) );
+ this->ImageOutlineActor->SetMapper( img_ol_mapper );
+ this->ImageOutlineActor->GetProperty( )->SetColor( 1, 1, 1 );
+ this->ImageOutlineActor->GetProperty( )->SetLineWidth( 1 );
- this->AddItem( this->ImageOutlineActor );
- */
+ this->AddItem( this->ImageOutlineActor );
+
+ } // fi
}
// eof - $RCSfile$