X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=lib%2FcpExtensions%2FVisualization%2FImageSliceActors.cxx;h=024fbe09c9bc94051ac6652bc19092f7de7a2424;hb=e2fb8817731f6231d34941a208e46b36dad425b2;hp=ee4947e834d4394463a547de1e958d4752e66f58;hpb=ce7d718b88d08f559ca8328b335c4fe8fbfb9b9a;p=cpPlugins.git diff --git a/lib/cpExtensions/Visualization/ImageSliceActors.cxx b/lib/cpExtensions/Visualization/ImageSliceActors.cxx index ee4947e..024fbe0 100644 --- a/lib/cpExtensions/Visualization/ImageSliceActors.cxx +++ b/lib/cpExtensions/Visualization/ImageSliceActors.cxx @@ -1,11 +1,13 @@ #include +#include #include #include #include #include #include +#include #include #include #include @@ -14,6 +16,15 @@ #include #include +// ------------------------------------------------------------------------- +double cpExtensions::Visualization::ImageSliceActors:: +m_PlaneColors[ 3 ][ 3 ] = +{ + { 1, 0, 0 }, + { 0, 1, 0 }, + { 0, 0, 1 } +}; + // ------------------------------------------------------------------------- cpExtensions::Visualization::ImageSliceActors* cpExtensions::Visualization::ImageSliceActors:: @@ -22,166 +33,103 @@ New( ) return( new Self( ) ); } -// ------------------------------------------------------------------------- -cpExtensions::Visualization:: -ImageBlender* cpExtensions::Visualization::ImageSliceActors:: -GetBlender( ) -{ - return( this->m_ImageBlender ); -} - -// ------------------------------------------------------------------------- -const cpExtensions::Visualization:: -ImageBlender* cpExtensions::Visualization::ImageSliceActors:: -GetBlender( ) const -{ - return( this->m_ImageBlender ); -} - -// ------------------------------------------------------------------------- -void cpExtensions::Visualization::ImageSliceActors:: -SetBlender( ImageBlender* blender ) -{ - this->m_ImageBlender = blender; - - // Create mapper and actors - this->m_ImageMapper = vtkSmartPointer< vtkImageSliceMapper >::New( ); - this->m_ImageMapper->SetInputConnection( - this->m_ImageBlender->GetOutputPort( ) - ); - this->m_ImageMapper->SetOrientation( 0 ); - this->m_ImageMapper->Update( ); - - // Create actor - this->m_ImageActor = vtkSmartPointer< vtkImageActor >::New( ); - this->m_ImageActor->SetMapper( this->m_ImageMapper ); - this->m_ImageActor->SetInterpolate( this->m_Interpolate ); - this->m_ImageActor->Modified( ); - - if( this->m_Style.GetPointer( ) != NULL ) - this->m_Style->AssociateImageActor( this->m_ImageActor ); - this->AddItem( this->m_ImageActor ); - - this->SetSliceNumber( this->GetSliceNumberMinValue( ) ); - this->ResetCursor( ); - this->Modified( ); -} - // ------------------------------------------------------------------------- void cpExtensions::Visualization::ImageSliceActors:: SetAxis( int axis ) { this->m_ImageMapper->SetOrientation( axis ); + this->m_BlenderMapper->SetOrientation( axis ); this->m_ImageMapper->Update( ); + this->m_BlenderMapper->Update( ); this->SetSliceNumber( this->GetSliceNumberMinValue( ) ); this->m_ImageActor->Modified( ); + this->m_BlenderActor->Modified( ); this->Modified( ); - this->ResetCamera( ); } // ------------------------------------------------------------------------- void cpExtensions::Visualization::ImageSliceActors:: -AddInputConnection( vtkAlgorithmOutput* aout ) +SetInputConnection( vtkAlgorithmOutput* aout ) { // Get input vtkImageData if( aout == NULL ) return; - if( this->m_ImageBlender.GetPointer( ) == NULL ) - { - // Try to infere if input comes from a color mapping filter - vtkAlgorithm* producer = aout->GetProducer( ); - vtkSmartPointer< ImageBlender > blender = - dynamic_cast< ImageBlender* >( producer ); - if( blender.GetPointer( ) == NULL ) - { - // Ok, create a new blender with default window level - vtkImageData* data = dynamic_cast< vtkImageData* >( - producer->GetOutputDataObject( aout->GetIndex( ) ) - ); - double r[ 2 ]; - data->GetScalarRange( r ); - - blender = vtkSmartPointer< ImageBlender >::New( ); - blender->AddInputConnection( aout ); - blender->SetWindow( r[ 1 ] - r[ 0 ] ); - blender->SetLevel( ( r[ 1 ] + r[ 0 ] ) / double( 2 ) ); - - } // fi - this->SetBlender( blender ); - } - else - this->m_ImageBlender->AddInputConnection( aout ); + // Create mapper and actors + this->m_ImageMapper->SetInputConnection( aout ); + this->_ConfigureInputImage( ); } // ------------------------------------------------------------------------- void cpExtensions::Visualization::ImageSliceActors:: -AddInputData( vtkImageData* data ) +SetInputImage( vtkImageData* data ) { // Get input vtkImageData if( data == NULL ) return; - if( this->m_ImageBlender.GetPointer( ) == NULL ) - { - // Create a new blender with default window level - double r[ 2 ]; - data->GetScalarRange( r ); + // Create mapper and actors + this->m_ImageMapper->SetInputData( data ); + this->_ConfigureInputImage( ); +} - vtkSmartPointer< ImageBlender > blender = - vtkSmartPointer< ImageBlender >::New( ); - blender->AddInputData( data ); - blender->SetWindow( r[ 1 ] - r[ 0 ] ); - blender->SetLevel( ( r[ 1 ] + r[ 0 ] ) / double( 2 ) ); - this->SetBlender( blender ); - } - else - this->m_ImageBlender->AddInputData( data ); +// ------------------------------------------------------------------------- +int cpExtensions::Visualization::ImageSliceActors:: +AddBinaryConnection( + vtkAlgorithmOutput* aout, + const double& r, const double& g, const double& b + ) +{ + if( aout == NULL ) + return( -1 ); + this->m_Blender->AddInputConnection( aout ); + this->_ConfigureBinaryImage( r, g, b ); + return( this->m_Blender->GetNumberOfInputs( ) - 1 ); +} + +// ------------------------------------------------------------------------- +int cpExtensions::Visualization::ImageSliceActors:: +AddBinaryImage( + vtkImageData* data, + const double& r, const double& g, const double& b + ) +{ + if( data == NULL ) + return( -1 ); + this->m_Blender->AddInputData( data ); + this->_ConfigureBinaryImage( r, g, b ); + return( this->m_Blender->GetNumberOfInputs( ) - 1 ); } // ------------------------------------------------------------------------- void cpExtensions::Visualization::ImageSliceActors:: Clear( ) { - // Reset values - this->m_VisibleExtent[ 0 ] = - this->m_VisibleExtent[ 2 ] = - this->m_VisibleExtent[ 4 ] = -1; - this->m_VisibleExtent[ 1 ] = - this->m_VisibleExtent[ 3 ] = - this->m_VisibleExtent[ 5 ] = 0; - this->m_VisibleBounds[ 0 ] = - this->m_VisibleBounds[ 2 ] = - this->m_VisibleBounds[ 4 ] = double( 0 ); - this->m_VisibleBounds[ 1 ] = - this->m_VisibleBounds[ 3 ] = - this->m_VisibleBounds[ 5 ] = double( 0 ); - // Unbind from container this->RemoveAllItems( ); - // Delete all images - this->m_ImageActor = NULL; - this->m_ImageMapper = NULL; - this->m_ImageBlender = NULL; - - // Reconfigure unique objects - this->m_Cursor = vtkSmartPointer< vtkPolyData >::New( ); - this->m_CursorMapper = vtkSmartPointer< vtkPolyDataMapper >::New( ); - this->m_CursorActor = vtkSmartPointer< vtkActor >::New( ); - this->m_HorizontalLine = vtkSmartPointer< vtkPolyData >::New( ); - this->m_HorizontalLineMapper = vtkSmartPointer< vtkPolyDataMapper >::New( ); - this->m_HorizontalLineActor = vtkSmartPointer< vtkActor >::New( ); - this->m_VerticalLine = vtkSmartPointer< vtkPolyData >::New( ); - this->m_VerticalLineMapper = vtkSmartPointer< vtkPolyDataMapper >::New( ); - this->m_VerticalLineActor = vtkSmartPointer< vtkActor >::New( ); - this->m_Plane = vtkSmartPointer< vtkPolyData >::New( ); - this->m_PlaneMapper = vtkSmartPointer< vtkPolyDataMapper >::New( ); - this->m_TextActor = vtkSmartPointer< vtkTextActor >::New( ); - this->m_PlaneActor = vtkSmartPointer< vtkActor >::New( ); - this->m_TextBuffer[ 0 ] = '\0'; + // Filters and mappers + this->m_ImageMapper = vtkSmartPointer< vtkImageSliceMapper >::New( ); + this->m_ImageActor = vtkSmartPointer< vtkImageActor >::New( ); + this->m_Blender = vtkSmartPointer< TBlender >::New( ); + this->m_BlenderMapper = vtkSmartPointer< vtkImageSliceMapper >::New( ); + this->m_BlenderLUT = vtkSmartPointer< vtkLookupTable >::New( ); + this->m_BlenderActor = vtkSmartPointer< vtkImageActor >::New( ); + + this->m_ImageActor->SetMapper( this->m_ImageMapper ); + this->m_BlenderMapper-> + SetInputConnection( this->m_Blender->GetOutputPort( ) ); + + this->m_BlenderLUT->SetNumberOfTableValues( 1 ); + this->m_BlenderLUT->SetTableRange( 0, 1 ); + this->m_BlenderLUT->SetTableValue( 0, 0, 0, 0, 0 ); + this->m_BlenderLUT->Build( ); - // Unique objects configuration + this->m_BlenderActor->SetMapper( this->m_BlenderMapper ); + this->m_BlenderActor->GetProperty( )->SetLookupTable( this->m_BlenderLUT ); + this->m_BlenderActor->GetProperty( )->UseLookupTableScalarRangeOn( ); + + // White cursor vtkSmartPointer< vtkPoints > cursor_points = vtkSmartPointer< vtkPoints >::New( ); vtkSmartPointer< vtkCellArray > cursor_lines = @@ -190,26 +138,28 @@ Clear( ) cursor_points->InsertNextPoint( 0, 0, 0 ); cursor_points->InsertNextPoint( 0, 0, 0 ); cursor_points->InsertNextPoint( 0, 0, 0 ); - cursor_points->InsertNextPoint( 0, 0, 0 ); - cursor_points->InsertNextPoint( 0, 0, 0 ); - cursor_points->InsertNextPoint( 0, 0, 0 ); - cursor_points->InsertNextPoint( 0, 0, 0 ); cursor_lines->InsertNextCell( 2 ); cursor_lines->InsertCellPoint( 0 ); cursor_lines->InsertCellPoint( 1 ); cursor_lines->InsertNextCell( 2 ); cursor_lines->InsertCellPoint( 2 ); cursor_lines->InsertCellPoint( 3 ); - cursor_lines->InsertNextCell( 2 ); - cursor_lines->InsertCellPoint( 4 ); - cursor_lines->InsertCellPoint( 5 ); - cursor_lines->InsertNextCell( 2 ); - cursor_lines->InsertCellPoint( 6 ); - cursor_lines->InsertCellPoint( 7 ); + + this->m_Cursor = vtkSmartPointer< vtkPolyData >::New( ); + this->m_CursorMapper = vtkSmartPointer< vtkPolyDataMapper >::New( ); + this->m_CursorActor = vtkSmartPointer< vtkActor >::New( ); + this->m_Cursor->SetPoints( cursor_points ); this->m_Cursor->SetLines( cursor_lines ); this->m_CursorMapper->SetInputData( this->m_Cursor ); this->m_CursorActor->SetMapper( this->m_CursorMapper ); + this->m_CursorActor->GetProperty( )->SetColor( 1, 1, 0 ); + this->m_CursorActor->GetProperty( )->SetLineWidth( 2 ); + + // Plane intersections + this->m_Axis1 = vtkSmartPointer< vtkPolyData >::New( ); + this->m_Axis1Mapper = vtkSmartPointer< vtkPolyDataMapper >::New( ); + this->m_Axis1Actor = vtkSmartPointer< vtkActor >::New( ); vtkSmartPointer< vtkPoints > h_points = vtkSmartPointer< vtkPoints >::New( ); @@ -220,10 +170,14 @@ Clear( ) h_lines->InsertNextCell( 2 ); h_lines->InsertCellPoint( 0 ); h_lines->InsertCellPoint( 1 ); - this->m_HorizontalLine->SetPoints( h_points ); - this->m_HorizontalLine->SetLines( h_lines ); - this->m_HorizontalLineMapper->SetInputData( this->m_HorizontalLine ); - this->m_HorizontalLineActor->SetMapper( this->m_HorizontalLineMapper ); + this->m_Axis1->SetPoints( h_points ); + this->m_Axis1->SetLines( h_lines ); + this->m_Axis1Mapper->SetInputData( this->m_Axis1 ); + this->m_Axis1Actor->SetMapper( this->m_Axis1Mapper ); + + this->m_Axis2 = vtkSmartPointer< vtkPolyData >::New( ); + this->m_Axis2Mapper = vtkSmartPointer< vtkPolyDataMapper >::New( ); + this->m_Axis2Actor = vtkSmartPointer< vtkActor >::New( ); vtkSmartPointer< vtkPoints > v_points = vtkSmartPointer< vtkPoints >::New( ); @@ -234,16 +188,20 @@ Clear( ) v_lines->InsertNextCell( 2 ); v_lines->InsertCellPoint( 0 ); v_lines->InsertCellPoint( 1 ); - this->m_VerticalLine->SetPoints( v_points ); - this->m_VerticalLine->SetLines( v_lines ); - this->m_VerticalLineMapper->SetInputData( this->m_VerticalLine ); - this->m_VerticalLineActor->SetMapper( this->m_VerticalLineMapper ); + this->m_Axis2->SetPoints( v_points ); + this->m_Axis2->SetLines( v_lines ); + this->m_Axis2Mapper->SetInputData( this->m_Axis2 ); + this->m_Axis2Actor->SetMapper( this->m_Axis2Mapper ); + + // Plane + this->m_Plane = vtkSmartPointer< vtkPolyData >::New( ); + this->m_PlaneMapper = vtkSmartPointer< vtkPolyDataMapper >::New( ); + this->m_PlaneActor = vtkSmartPointer< vtkActor >::New( ); vtkSmartPointer< vtkPoints > plane_points = vtkSmartPointer< vtkPoints >::New( ); vtkSmartPointer< vtkCellArray > plane_lines = vtkSmartPointer< vtkCellArray >::New( ); - plane_points->InsertNextPoint( 0, 0, 0 ); plane_points->InsertNextPoint( 0, 1, 0 ); plane_points->InsertNextPoint( 1, 1, 0 ); @@ -260,9 +218,12 @@ Clear( ) this->m_PlaneMapper->SetInputData( this->m_Plane ); this->m_PlaneActor->SetMapper( this->m_PlaneMapper ); + // Text actor + this->m_TextActor = vtkSmartPointer< vtkTextActor >::New( ); + this->m_TextBuffer[ 0 ] = '\0'; this->m_TextActor->SetTextScaleModeToNone( ); vtkTextProperty* textprop = this->m_TextActor->GetTextProperty( ); - textprop->SetColor( 1, 1, 1 ); + textprop->SetColor( 1, 1, 0 ); textprop->SetFontFamilyToCourier( ); textprop->SetFontSize( 18 ); textprop->BoldOff( ); @@ -273,13 +234,6 @@ Clear( ) vtkCoordinate* coord = this->m_TextActor->GetPositionCoordinate( ); coord->SetCoordinateSystemToNormalizedViewport( ); coord->SetValue( 0.01, 0.01 ); - - // Update actor collection - this->AddItem( this->m_CursorActor ); - this->AddItem( this->m_HorizontalLineActor ); - this->AddItem( this->m_VerticalLineActor ); - this->AddItem( this->m_TextActor ); - this->AddItem( this->m_PlaneActor ); } // ------------------------------------------------------------------------- @@ -306,26 +260,20 @@ GetStyle( ) const // ------------------------------------------------------------------------- vtkImageData* cpExtensions::Visualization::ImageSliceActors:: -GetInputImage( unsigned int id ) +GetInputImage( ) { - if( this->m_ImageBlender.GetPointer( ) != NULL ) - return( - dynamic_cast< vtkImageData* >( this->m_ImageBlender->GetInput( id ) ) - ); + if( this->m_ImageMapper.GetPointer( ) != NULL ) + return( this->m_ImageMapper->GetInput( ) ); else return( NULL ); } // ------------------------------------------------------------------------- const vtkImageData* cpExtensions::Visualization::ImageSliceActors:: -GetInputImage( unsigned int id ) const +GetInputImage( ) const { - if( this->m_ImageBlender.GetPointer( ) != NULL ) - return( - dynamic_cast< const vtkImageData* >( - this->m_ImageBlender->GetInput( ) - ) - ); + if( this->m_ImageMapper.GetPointer( ) != NULL ) + return( this->m_ImageMapper->GetInput( ) ); else return( NULL ); } @@ -347,15 +295,15 @@ PushActorsInto( vtkRenderWindow* window, bool force_style ) rwi->SetInteractorStyle( this->m_Style ); // Update actors - vtkProp* prop; - this->InitTraversal( ); - while( prop = this->GetNextProp( ) ) - renderer->AddViewProp( prop ); - renderer->Modified( ); - if( !force_style ) + renderer->AddViewProp( this->m_ImageActor ); + renderer->AddViewProp( this->m_BlenderActor ); + if( force_style ) { - renderer->RemoveViewProp( this->m_CursorActor ); - renderer->RemoveViewProp( this->m_TextActor ); + renderer->AddViewProp( this->m_CursorActor ); + renderer->AddViewProp( this->m_PlaneActor ); + renderer->AddViewProp( this->m_TextActor ); + renderer->AddViewProp( this->m_Axis1Actor ); + renderer->AddViewProp( this->m_Axis2Actor ); } // fi @@ -379,14 +327,13 @@ PushActorsInto( vtkRenderWindow* window, bool force_style ) } else // if( axis == 2 ) { - camera->SetPosition( double( 0 ), double( 0 ), double( 1 ) ); - camera->SetViewUp ( double( 0 ), double( -1 ), double( 0 ) ); + camera->SetPosition( double( 0 ), double( 0 ), double( -1 ) ); + camera->SetViewUp ( double( 0 ), double( -1 ), double( 0 ) ); } // fi } // fi - this->ResetCamera( ); - rwi->Render( ); + renderer->ResetCamera( ); } // ------------------------------------------------------------------------- @@ -395,19 +342,24 @@ PopActorsFrom( vtkRenderWindow* window ) { vtkRenderWindowInteractor* rwi = window->GetInteractor( ); vtkRenderer* renderer = window->GetRenderers( )->GetFirstRenderer( ); - if( renderer != NULL ) { - // Update actors - vtkProp* prop; - this->InitTraversal( ); - while( prop = this->GetNextProp( ) ) - renderer->RemoveViewProp( prop ); - renderer->Modified( ); + renderer->RemoveViewProp( this->m_ImageActor ); + renderer->RemoveViewProp( this->m_BlenderActor ); + renderer->RemoveViewProp( this->m_CursorActor ); + renderer->RemoveViewProp( this->m_PlaneActor ); + renderer->RemoveViewProp( this->m_TextActor ); + renderer->RemoveViewProp( this->m_Axis1Actor ); + renderer->RemoveViewProp( this->m_Axis2Actor ); } // fi - if( rwi != NULL ) - rwi->Render( ); +} + +// ------------------------------------------------------------------------- +unsigned int cpExtensions::Visualization::ImageSliceActors:: +GetNumberOfImages( ) const +{ + return( this->m_Blender->GetNumberOfInputs( ) ); } // ------------------------------------------------------------------------- @@ -424,6 +376,20 @@ GetImageActor( ) const return( this->m_ImageActor ); } +// ------------------------------------------------------------------------- +vtkImageActor* cpExtensions::Visualization::ImageSliceActors:: +GetBinaryActor( ) +{ + return( this->m_BlenderActor ); +} + +// ------------------------------------------------------------------------- +const vtkImageActor* cpExtensions::Visualization::ImageSliceActors:: +GetBinaryActor( ) const +{ + return( this->m_BlenderActor ); +} + // ------------------------------------------------------------------------- vtkTextActor* cpExtensions::Visualization::ImageSliceActors:: GetTextActor( ) @@ -479,6 +445,7 @@ SetInterpolate( bool v ) if( this->m_Interpolate != v ) { this->m_ImageActor->SetInterpolate( v ); + this->m_BlenderActor->SetInterpolate( v ); this->m_Interpolate = v; this->Modified( ); @@ -515,7 +482,7 @@ GetDisplayBounds( double bounds[ 6 ] ) const { if( this->m_ImageActor.GetPointer( ) == NULL ) { - bounds[ 0 ] = bounds[ 2 ] = bounds[ 4 ] = double( -1 ); + bounds[ 0 ] = bounds[ 2 ] = bounds[ 4 ] = double( 0 ); bounds[ 1 ] = bounds[ 3 ] = bounds[ 5 ] = double( -1 ); } else @@ -528,8 +495,6 @@ ResetCursor( ) { if( this->m_ImageMapper.GetPointer( ) != NULL ) { - double bounds[ 6 ]; - this->m_ImageMapper->GetInput( )->GetBounds( bounds ); double pos[] = { this->m_VisibleBounds[ 0 ], @@ -545,10 +510,6 @@ ResetCursor( ) points->SetPoint( 1, 0, 0, 0 ); points->SetPoint( 2, 0, 0, 0 ); points->SetPoint( 3, 0, 0, 0 ); - points->SetPoint( 4, 0, 0, 0 ); - points->SetPoint( 5, 0, 0, 0 ); - points->SetPoint( 6, 0, 0, 0 ); - points->SetPoint( 7, 0, 0, 0 ); this->m_Cursor->Modified( ); this->m_CursorMapper->Modified( ); this->m_CursorActor->Modified( ); @@ -573,79 +534,137 @@ SetCursor( double pos[ 3 ] ) // Update cross double* bounds = this->m_VisibleBounds; - double - p0[ 3 ], p1[ 3 ], p2[ 3 ], p3[ 3 ], - p4[ 3 ], p5[ 3 ], p6[ 3 ], p7[ 3 ]; + double p0[ 3 ], p1[ 3 ], p2[ 3 ], p3[ 3 ]; - p0[ a2 ] = p1[ a2 ] = p4[ a2 ] = p5[ a2 ] = pos[ a2 ]; - p2[ a1 ] = p3[ a1 ] = p6[ a1 ] = p7[ a1 ] = pos[ a1 ]; + p0[ a2 ] = p1[ a2 ] = pos[ a2 ]; + p2[ a1 ] = p3[ a1 ] = pos[ a1 ]; p0[ a0 ] = p1[ a0 ] = p2[ a0 ] = p3[ a0 ] = bounds[ ma0 ]; - p4[ a0 ] = p5[ a0 ] = p6[ a0 ] = p7[ a0 ] = bounds[ ma0 + 1 ]; - p0[ a1 ] = p4[ a1 ] = bounds[ ma1 ]; - p1[ a1 ] = p5[ a1 ] = bounds[ ma1 + 1 ]; - p2[ a2 ] = p6[ a2 ] = bounds[ ma2 ]; - p3[ a2 ] = p7[ a2 ] = bounds[ ma2 + 1 ]; + p0[ a1 ] = bounds[ ma1 ]; + p1[ a1 ] = bounds[ ma1 + 1 ]; + p2[ a2 ] = bounds[ ma2 ]; + p3[ a2 ] = bounds[ ma2 + 1 ]; vtkPoints* points = this->m_Cursor->GetPoints( ); points->SetPoint( 0, p0 ); points->SetPoint( 1, p1 ); points->SetPoint( 2, p2 ); points->SetPoint( 3, p3 ); - points->SetPoint( 4, p4 ); - points->SetPoint( 5, p5 ); - points->SetPoint( 6, p6 ); - points->SetPoint( 7, p7 ); this->m_Cursor->Modified( ); this->m_CursorMapper->Modified( ); this->m_CursorActor->Modified( ); } +// ------------------------------------------------------------------------- +void cpExtensions::Visualization::ImageSliceActors:: +ResetAxesCursor( ) +{ + if( this->m_ImageMapper.GetPointer( ) != NULL ) + { + double pos[] = + { + this->m_VisibleBounds[ 0 ], + this->m_VisibleBounds[ 2 ], + this->m_VisibleBounds[ 4 ] + }; + this->SetAxesCursor( pos ); + } + else + { + vtkPoints* points = this->m_Axis1->GetPoints( ); + points->SetPoint( 0, 0, 0, 0 ); + points->SetPoint( 1, 0, 0, 0 ); + this->m_Axis1->Modified( ); + this->m_Axis1Mapper->Modified( ); + this->m_Axis1Actor->Modified( ); + + points = this->m_Axis2->GetPoints( ); + points->SetPoint( 0, 0, 0, 0 ); + points->SetPoint( 1, 0, 0, 0 ); + this->m_Axis2->Modified( ); + this->m_Axis2Mapper->Modified( ); + this->m_Axis2Actor->Modified( ); + + } // fi +} + +// ------------------------------------------------------------------------- +void cpExtensions::Visualization::ImageSliceActors:: +SetAxesCursor( double pos[ 3 ] ) +{ + if( this->m_ImageMapper.GetPointer( ) == NULL ) + return; + + // Get ordered axes + int a0 = this->GetAxis( ); + int a1 = ( a0 + 1 ) % 3; + int a2 = ( a0 + 2 ) % 3; + int ma0 = a0 << 1; + int ma1 = a1 << 1; + int ma2 = a2 << 1; + + // Update cross + double* bounds = this->m_VisibleBounds; + double p0[ 3 ], p1[ 3 ], p2[ 3 ], p3[ 3 ]; + + p0[ a2 ] = p1[ a2 ] = pos[ a2 ]; + p2[ a1 ] = p3[ a1 ] = pos[ a1 ]; + p0[ a0 ] = p1[ a0 ] = p2[ a0 ] = p3[ a0 ] = bounds[ ma0 ]; + p0[ a1 ] = bounds[ ma1 ]; + p1[ a1 ] = bounds[ ma1 + 1 ]; + p2[ a2 ] = bounds[ ma2 ]; + p3[ a2 ] = bounds[ ma2 + 1 ]; + + vtkPoints* points1 = this->m_Axis1->GetPoints( ); + points1->SetPoint( 0, p2 ); + points1->SetPoint( 1, p3 ); + + vtkPoints* points2 = this->m_Axis2->GetPoints( ); + points2->SetPoint( 0, p0 ); + points2->SetPoint( 1, p1 ); + + this->m_Axis1->Modified( ); + this->m_Axis1Mapper->Modified( ); + this->m_Axis1Actor->Modified( ); + + this->m_Axis2->Modified( ); + this->m_Axis2Mapper->Modified( ); + this->m_Axis2Actor->Modified( ); +} + // ------------------------------------------------------------------------- double cpExtensions::Visualization::ImageSliceActors:: GetMinWindow( ) const { - if( this->m_ImageBlender.GetPointer( ) != NULL ) - return( this->m_ImageBlender->GetMinWindow( ) ); - else - return( double( 0 ) ); + return( this->m_WLRange[ 0 ] ); } // ------------------------------------------------------------------------- double cpExtensions::Visualization::ImageSliceActors:: GetMaxWindow( ) const { - if( this->m_ImageBlender.GetPointer( ) != NULL ) - return( this->m_ImageBlender->GetMaxWindow( ) ); - else - return( double( 0 ) ); + return( this->m_WLRange[ 1 ] ); } // ------------------------------------------------------------------------- double cpExtensions::Visualization::ImageSliceActors:: GetMinLevel( ) const { - if( this->m_ImageBlender.GetPointer( ) != NULL ) - return( this->m_ImageBlender->GetMinLevel( ) ); - else - return( double( 0 ) ); + return( this->m_WLRange[ 2 ] ); } // ------------------------------------------------------------------------- double cpExtensions::Visualization::ImageSliceActors:: GetMaxLevel( ) const { - if( this->m_ImageBlender.GetPointer( ) != NULL ) - return( this->m_ImageBlender->GetMaxLevel( ) ); - else - return( double( 0 ) ); + return( this->m_WLRange[ 3 ] ); } // ------------------------------------------------------------------------- double cpExtensions::Visualization::ImageSliceActors:: GetWindow( ) const { - if( this->m_ImageBlender.GetPointer( ) != NULL ) - return( this->m_ImageBlender->GetWindow( ) ); + if( this->m_ImageActor.GetPointer( ) != NULL ) + return( this->m_ImageActor->GetProperty( )->GetColorWindow( ) ); else return( double( 0 ) ); } @@ -654,8 +673,8 @@ GetWindow( ) const double cpExtensions::Visualization::ImageSliceActors:: GetLevel( ) const { - if( this->m_ImageBlender.GetPointer( ) != NULL ) - return( this->m_ImageBlender->GetLevel( ) ); + if( this->m_ImageActor.GetPointer( ) != NULL ) + return( this->m_ImageActor->GetProperty( )->GetColorLevel( ) ); else return( double( 0 ) ); } @@ -664,45 +683,47 @@ GetLevel( ) const void cpExtensions::Visualization::ImageSliceActors:: SetWindow( double w ) { - if( this->m_ImageBlender.GetPointer( ) != NULL ) - { - this->m_ImageBlender->SetWindow( w ); - this->Modified( ); - - } // fi + if( this->m_ImageActor.GetPointer( ) == NULL ) + return; + double v = ( w < this->m_WLRange[ 0 ] )? this->m_WLRange[ 0 ]: w; + v = ( v > this->m_WLRange[ 1 ] )? this->m_WLRange[ 1 ]: v; + this->m_ImageActor->GetProperty( )->SetColorWindow( v ); } // ------------------------------------------------------------------------- void cpExtensions::Visualization::ImageSliceActors:: SetLevel( double l ) { - if( this->m_ImageBlender.GetPointer( ) != NULL ) - { - this->m_ImageBlender->SetLevel( l ); - this->Modified( ); - - } // fi + if( this->m_ImageActor.GetPointer( ) == NULL ) + return; + double v = ( l < this->m_WLRange[ 2 ] )? this->m_WLRange[ 2 ]: l; + v = ( v > this->m_WLRange[ 3 ] )? this->m_WLRange[ 3 ]: v; + this->m_ImageActor->GetProperty( )->SetColorLevel( v ); } // ------------------------------------------------------------------------- void cpExtensions::Visualization::ImageSliceActors:: SetWindowLevel( double w, double l ) { - if( this->m_ImageBlender.GetPointer( ) != NULL ) - { - this->m_ImageBlender->SetWindowLevel( w, l ); - this->UpdateText( w, l ); - this->Modified( ); - - } // fi + if( this->m_ImageActor.GetPointer( ) == NULL ) + return; + double a = ( w < this->m_WLRange[ 0 ] )? this->m_WLRange[ 0 ]: w; + a = ( a > this->m_WLRange[ 1 ] )? this->m_WLRange[ 1 ]: a; + double b = ( l < this->m_WLRange[ 2 ] )? this->m_WLRange[ 2 ]: l; + b = ( b > this->m_WLRange[ 3 ] )? this->m_WLRange[ 3 ]: b; + this->m_ImageActor->GetProperty( )->SetColorWindow( a ); + this->m_ImageActor->GetProperty( )->SetColorLevel( b ); + this->UpdateText( a, b ); } // ------------------------------------------------------------------------- void cpExtensions::Visualization::ImageSliceActors:: ResetWindowLevel( ) { - std::cerr << "ACA Resetear" << std::endl; - std::exit( 1 ); + this->SetWindowLevel( + this->m_WLRange[ 1 ] * double( 0.5 ), + ( this->m_WLRange[ 3 ] + this->m_WLRange[ 2 ] ) * double( 0.5 ) + ); } // ------------------------------------------------------------------------- @@ -757,9 +778,11 @@ SetSliceNumber( const int& slice ) // Update mappers and display bounds this->m_ImageMapper->SetSliceNumber( slice ); + this->m_BlenderMapper->SetSliceNumber( slice ); this->m_ImageMapper->Modified( ); + this->m_BlenderMapper->Modified( ); this->m_ImageActor->Modified( ); - this->m_ImageMapper->Update( ); + this->m_BlenderActor->Modified( ); // Update display extent (this isn't done automatically) this->m_ImageMapper->GetInput( )->GetExtent( this->m_VisibleExtent ); @@ -806,43 +829,36 @@ SetSliceNumber( const int& slice ) { plane_points->SetPoint( 1, bnds[ 1 ], bnds[ 2 ], bnds[ 5 ] ); plane_points->SetPoint( 3, bnds[ 0 ], bnds[ 3 ], bnds[ 4 ] ); - if( axis == 0 ) - this->m_PlaneActor->GetProperty( )->SetColor( 1, 0, 0 ); - else - this->m_PlaneActor->GetProperty( )->SetColor( 0, 0, 1 ); } else if( axis == 1 ) // ZX, y-normal { plane_points->SetPoint( 1, bnds[ 0 ], bnds[ 2 ], bnds[ 5 ] ); plane_points->SetPoint( 3, bnds[ 1 ], bnds[ 3 ], bnds[ 4 ] ); - this->m_PlaneActor->GetProperty( )->SetColor( 0, 1, 0 ); } // fi + + // Set plane colors + this->m_PlaneActor->GetProperty( )-> + SetColor( this->m_PlaneColors[ axis ] ); + this->m_Axis1Actor->GetProperty( )-> + SetColor( this->m_PlaneColors[ ( axis + 1 ) % 3 ] ); + this->m_Axis2Actor->GetProperty( )-> + SetColor( this->m_PlaneColors[ ( axis + 2 ) % 3 ] ); + this->m_Plane->Modified( ); this->m_PlaneMapper->Modified( ); this->m_PlaneActor->Modified( ); - // Update text - this->UpdateText( ); - - // Update lines from associated slices - /* TODO - auto sIt = this->m_AssociatedSlices.begin( ); - for( ; sIt != this->m_AssociatedSlices.end( ); ++sIt ) - { - Self* slice = *sIt; - for( unsigned int id = 0; id < slice->m_AssociatedSlices.size( ); ++id ) - { - std::cout << id << std::endl; - if( slice->m_AssociatedSlices[ id ] != this ) - continue; + this->m_Axis1->Modified( ); + this->m_Axis1Mapper->Modified( ); + this->m_Axis1Actor->Modified( ); - std::cout << "id : " << id << std::endl; + this->m_Axis2->Modified( ); + this->m_Axis2Mapper->Modified( ); + this->m_Axis2Actor->Modified( ); - } // rof - - } // rof - */ + // Update text + this->UpdateText( ); // Update camera position if( this->m_Window == NULL ) @@ -864,7 +880,7 @@ SetSliceNumber( const int& slice ) void cpExtensions::Visualization::ImageSliceActors:: SetSlice( double* pos ) { - vtkImageData* image = this->GetInputImage( 0 ); + vtkImageData* image = this->GetInputImage( ); if( image == NULL ) return; @@ -886,10 +902,17 @@ UpdateText( ) else if( axId == 1 ) axis = 'Y'; else if( axId == 2 ) axis = 'Z'; +#if defined(WIN32) + sprintf_s( + this->m_TextBuffer, 1024, "Axis: %c (%d)", + axis, this->m_ImageMapper->GetSliceNumber( ) + ); +#else // defined(WIN32) std::sprintf( this->m_TextBuffer, "Axis: %c (%d)", axis, this->m_ImageMapper->GetSliceNumber( ) ); +#endif // defined(WIN32) } else this->m_TextBuffer[ 0 ] = '\0'; @@ -911,7 +934,7 @@ UpdateText( double pos[ 3 ] ) else if( axId == 2 ) axis = 'Z'; int slice = this->GetSliceNumber( ); - vtkImageData* image = this->GetInputImage( 0 ); + vtkImageData* image = this->GetInputImage( ); int ijk[ 3 ]; double pcoords[ 3 ]; image->ComputeStructuredCoordinates( pos, ijk, pcoords ); @@ -940,11 +963,18 @@ UpdateText( double pos[ 3 ] ) ijk[ 0 ], ijk[ 1 ], ijk[ 2 ], n ); str << ")"; + +#if defined(WIN32) + sprintf_s( + this->m_TextBuffer, 1024, "Axis: %c (%d)\nPixel %s", + axis, slice, str.str( ).c_str( ) + ); +#else // defined(WIN32) std::sprintf( this->m_TextBuffer, "Axis: %c (%d)\nPixel %s", axis, slice, str.str( ).c_str( ) ); - +#endif // defined(WIN32) } // fi } else @@ -966,10 +996,17 @@ UpdateText( const double& w, const double& l ) else if( axId == 1 ) axis = 'Y'; else if( axId == 2 ) axis = 'Z'; +#if defined(WIN32) + sprintf_s( + this->m_TextBuffer, 1024, "Axis: %c (%d)\nW/L (%.2f/%.2f)", + axis, this->m_ImageMapper->GetSliceNumber( ), w, l + ); +#else // defined(WIN32) std::sprintf( this->m_TextBuffer, "Axis: %c (%d)\nW/L (%.2f/%.2f)", axis, this->m_ImageMapper->GetSliceNumber( ), w, l ); +#endif // defined(WIN32) } else this->m_TextBuffer[ 0 ] = '\0'; @@ -980,17 +1017,16 @@ UpdateText( const double& w, const double& l ) // ------------------------------------------------------------------------- void cpExtensions::Visualization::ImageSliceActors:: -Render( const double& t ) +Render( ) { - if( this->m_Window != NULL ) - { - vtkRenderer* renderer = - this->m_Window->GetRenderers( )->GetFirstRenderer( ); - if( renderer != NULL ) - renderer->SetAllocatedRenderTime( t ); - this->m_Window->Render( ); - - } // fi + if( this->m_Window == NULL ) + return; + vtkRenderer* renderer = + this->m_Window->GetRenderers( )->GetFirstRenderer( ); + if( renderer == NULL ) + return; + renderer->SetAllocatedRenderTime( 1e-3 ); + this->m_Window->Render( ); } // ------------------------------------------------------------------------- @@ -1001,8 +1037,9 @@ ResetCamera( ) return; vtkRenderer* renderer = this->m_Window->GetRenderers( )->GetFirstRenderer( ); - if( renderer != NULL ) - renderer->ResetCamera( this->m_VisibleBounds ); + if( renderer == NULL ) + return; + renderer->ResetCamera( ); } // ------------------------------------------------------------------------- @@ -1028,6 +1065,7 @@ _ConfigureStyle( ) { // Connect this view with a controller this->m_Style = vtkSmartPointer< TStyle >::New( ); + this->m_Style->SetAssociatedObject( this ); this->m_Style->AddMouseMoveCommand( Self::_MouseMoveCommand, this ); this->m_Style->AddMouseClickCommand( Self::_MouseClickCommand, this ); this->m_Style->AddMouseWheelCommand( Self::_MouseWheelCommand, this ); @@ -1036,6 +1074,72 @@ _ConfigureStyle( ) this->m_Style->AddLeaveCommand( Self::_LeaveCommand, this ); } +// ------------------------------------------------------------------------- +void cpExtensions::Visualization::ImageSliceActors:: +_ConfigureInputImage( ) +{ + this->m_ImageMapper->SetOrientation( 0 ); + this->m_ImageMapper->Update( ); + + // Create actor + this->m_ImageActor->SetInterpolate( this->m_Interpolate ); + this->m_ImageActor->Modified( ); + + if( this->m_Style.GetPointer( ) != NULL ) + this->m_Style->AssociateImageActor( this->m_ImageActor ); + this->AddItem( this->m_ImageActor ); + + this->SetSliceNumber( this->GetSliceNumberMinValue( ) ); + this->ResetCursor( ); + this->Modified( ); + + // Reset cursors + this->ResetCursor( ); + this->ResetAxesCursor( ); + + // Update window/level ranges + vtkImageData* data = this->GetInputImage( ); + if( data != NULL ) + { + double r[ 2 ]; + data->GetScalarRange( r ); + this->m_WLRange[ 0 ] = double( 0 ); + this->m_WLRange[ 1 ] = r[ 1 ] - r[ 0 ]; + this->m_WLRange[ 2 ] = r[ 0 ]; + this->m_WLRange[ 3 ] = r[ 1 ]; + this->ResetWindowLevel( ); + + // Configure blender + this->m_BlenderBase = vtkSmartPointer< vtkImageData >::New( ); + this->m_BlenderBase->ShallowCopy( data ); + this->m_BlenderBase->AllocateScalars( VTK_UNSIGNED_CHAR, 1 ); + std::memset( + this->m_BlenderBase->GetScalarPointer( ), 0, + this->m_BlenderBase->GetActualMemorySize( ) + ); + this->m_Blender->AddInputData( this->m_BlenderBase ); + + } // fi +} + +// ------------------------------------------------------------------------- +void cpExtensions::Visualization::ImageSliceActors:: +_ConfigureBinaryImage( const double& r, const double& g, const double& b ) +{ + unsigned int nValues = this->m_BlenderLUT->GetNumberOfTableValues( ); + this->m_BlenderLUT->SetNumberOfTableValues( nValues + 1 ); + this->m_BlenderLUT->SetTableRange( 0, nValues ); + this->m_BlenderLUT->SetTableValue( nValues, r, g, b, 0.5 ); + this->m_BlenderLUT->Build( ); + + this->m_BlenderLUT->Modified( ); + this->m_Blender->Modified( ); + this->m_Blender->Update( ); + + this->m_BlenderMapper->Modified( ); + this->m_BlenderActor->Modified( ); +} + // ------------------------------------------------------------------------- void cpExtensions::Visualization::ImageSliceActors:: _MouseMoveCommand( @@ -1052,17 +1156,20 @@ _MouseMoveCommand( // Just show the pixel information actors->SetCursor( pos ); actors->UpdateText( pos ); - actors->Render( 1e-3 ); + actors->Render( ); } else if( btn == TStyle::ButtonID_Left ) { if( !alt && ctr && !sft ) { + // Show axes in current renderer + actors->SetAxesCursor( pos ); + actors->Render( ); + // Interactively move slices auto i = actors->m_SlicesCommands.begin( ); for( ; i != actors->m_SlicesCommands.end( ); ++i ) i->first( pos, actors->GetAxis( ), i->second ); - actors->Render( 1e-3 ); } // fi } @@ -1087,13 +1194,13 @@ _MouseMoveCommand( dx += actors->m_StartWindowLevel[ 0 ]; dy += actors->m_StartWindowLevel[ 1 ]; actors->SetWindowLevel( dx, dy ); - actors->Render( 1e-3 ); + actors->Render( ); // Associate objects auto i = actors->m_WindowLevelCommands.begin( ); for( ; i != actors->m_WindowLevelCommands.end( ); ++i ) i->first( dx, dy, i->second ); - + } // fi } // fi @@ -1136,16 +1243,12 @@ _MouseWheelCommand( if( slice > actors->GetSliceNumberMaxValue( ) ) slice = actors->GetSliceNumberMaxValue( ); actors->SetSliceNumber( slice ); + actors->Render( ); auto a = actors->m_AssociatedSlices.begin( ); for( ; a != actors->m_AssociatedSlices.end( ); ++a ) - { ( *a )->SetSliceNumber( slice ); - ( *a )->Render( 1e-3 ); - } // rof - actors->Render( 1e-3 ); - // Associate objects auto i = actors->m_RenderCommands.begin( ); for( ; i != actors->m_RenderCommands.end( ); ++i ) @@ -1167,7 +1270,7 @@ _KeyCommand( void* data, const char& key ) case 'r': case 'R': { actors->ResetCamera( ); - actors->Render( 1e-3 ); + actors->Render( ); // Associate objects auto i = actors->m_RenderCommands.begin( ); @@ -1178,7 +1281,7 @@ _KeyCommand( void* data, const char& key ) case 'w': case 'W': { actors->ResetWindowLevel( ); - actors->Render( 1e-3 ); + actors->Render( ); // Associate objects auto i = actors->m_RenderCommands.begin( ); @@ -1201,7 +1304,7 @@ _EnterCommand( void* data ) actors->ResetCursor( ); actors->m_CursorActor->VisibilityOn( ); - actors->Render( 1e-3 ); + actors->Render( ); } // ------------------------------------------------------------------------- @@ -1214,7 +1317,7 @@ _LeaveCommand( void* data ) actors->ResetCursor( ); actors->m_CursorActor->VisibilityOff( ); - actors->Render( 1e-3 ); + actors->Render( ); } // eof - $RCSfile$