X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=lib%2FcpExtensions%2FVisualization%2FImageSliceActors.cxx;h=d6542efbd25152f0870fa8b2a2645827d67407c4;hb=8b6eaedb25048704cb0a1a3fa1dd751c4fb25db4;hp=52be1f1ddeb099e82990b032f4e002038356c40b;hpb=1d735deae0ae5e8bff643a59370a03c961ef780c;p=cpPlugins.git diff --git a/lib/cpExtensions/Visualization/ImageSliceActors.cxx b/lib/cpExtensions/Visualization/ImageSliceActors.cxx index 52be1f1..d6542ef 100644 --- a/lib/cpExtensions/Visualization/ImageSliceActors.cxx +++ b/lib/cpExtensions/Visualization/ImageSliceActors.cxx @@ -1,128 +1,44 @@ #include -#include -#include +#include +#include #include -#include +#include #include -#include #include -#include +#include +#include #include // ------------------------------------------------------------------------- -cpExtensions::Visualization::ImageSliceActors* cpExtensions::Visualization::ImageSliceActors:: +Self* cpExtensions::Visualization::ImageSliceActors:: New( ) { return( new Self( ) ); } -// ------------------------------------------------------------------------- -void cpExtensions::Visualization::ImageSliceActors:: -AddInputConnection( vtkAlgorithmOutput* aout, int axis ) -{ - unsigned int nImages = this->SliceMappers.size( ); - - vtkSmartPointer< vtkImageSliceMapper > mapper = - vtkSmartPointer< vtkImageSliceMapper >::New( ); - this->SliceMappers.push_back( mapper ); - mapper->SetInputConnection( aout ); - mapper->SetOrientation( - ( nImages == 0 )? axis: this->SliceMappers[ 0 ]->GetOrientation( ) - ); - mapper->Update( ); - - vtkSmartPointer< vtkImageActor > actor = - vtkSmartPointer< vtkImageActor >::New( ); - this->ImageActors.push_back( actor ); - actor->SetMapper( mapper ); - actor->Modified( ); - - if( nImages == 0 ) - { - this->AddItem( this->TextActor ); - this->AddItem( this->PlaneActor ); - - } // fi - this->AddItem( actor ); - this->Modified( ); -} - -// ------------------------------------------------------------------------- -void cpExtensions::Visualization::ImageSliceActors:: -AddInputData( vtkImageData* data, int axis ) -{ - unsigned int nImages = this->SliceMappers.size( ); - - vtkSmartPointer< vtkImageSliceMapper > mapper = - vtkSmartPointer< vtkImageSliceMapper >::New( ); - this->SliceMappers.push_back( mapper ); - mapper->SetInputData( data ); - mapper->SetOrientation( - ( nImages == 0 )? axis: this->SliceMappers[ 0 ]->GetOrientation( ) - ); - mapper->Update( ); - - vtkSmartPointer< vtkImageActor > actor = - vtkSmartPointer< vtkImageActor >::New( ); - this->ImageActors.push_back( actor ); - actor->SetMapper( mapper ); - actor->Modified( ); - - if( nImages == 0 ) - { - this->AddItem( this->TextActor ); - this->AddItem( this->PlaneActor ); - - } // fi - this->AddItem( actor ); - this->Modified( ); -} - // ------------------------------------------------------------------------- void cpExtensions::Visualization::ImageSliceActors:: Clear( ) { - // Unbind from container this->RemoveAllItems( ); - // Delete all images - this->SliceMappers.clear( ); - this->ImageActors.clear( ); - - // Reconfigure unique objects - this->PlaneSource = vtkSmartPointer< vtkPolyData >::New( ); - this->PlaneMapper = vtkSmartPointer< vtkPolyDataMapper >::New( ); - this->TextActor = vtkSmartPointer< vtkTextActor >::New( ); - this->PlaneActor = vtkSmartPointer< vtkActor >::New( ); - this->TextBuffer[ 0 ] = '\0'; - - // Unique objects configuration - 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 ); - plane_points->InsertNextPoint( 1, 0, 0 ); - plane_lines->InsertNextCell( 5 ); - plane_lines->InsertCellPoint( 0 ); - plane_lines->InsertCellPoint( 1 ); - plane_lines->InsertCellPoint( 2 ); - plane_lines->InsertCellPoint( 3 ); - plane_lines->InsertCellPoint( 0 ); - this->PlaneSource->SetPoints( plane_points ); - this->PlaneSource->SetLines( plane_lines ); - - this->PlaneMapper->SetInputData( this->PlaneSource ); - this->PlaneActor->SetMapper( this->PlaneMapper ); - - this->TextActor->SetTextScaleModeToNone( ); - vtkTextProperty* textprop = this->TextActor->GetTextProperty( ); - textprop->SetColor( 1, 1, 1 ); + // Main image + this->m_Mapper = vtkSmartPointer< vtkImageSliceMapper >::New( ); + this->m_Actor = vtkSmartPointer< vtkImageActor >::New( ); + this->m_Actor->SetMapper( this->m_Mapper ); + + // Cursor + this->m_Cursor.Create( ); + this->m_Cursor.Actor->GetProperty( )->SetColor( 1, 1, 0 ); + + // Text + this->m_TextBuffer[ 0 ] = '\0'; + this->m_TextActor = vtkSmartPointer< vtkTextActor >::New( ); + this->m_TextActor->SetTextScaleModeToNone( ); + vtkTextProperty* textprop = this->m_TextActor->GetTextProperty( ); + textprop->SetColor( 1, 1, 0 ); textprop->SetFontFamilyToCourier( ); textprop->SetFontSize( 18 ); textprop->BoldOff( ); @@ -130,261 +46,582 @@ Clear( ) textprop->ShadowOff( ); textprop->SetJustificationToLeft( ); textprop->SetVerticalJustificationToBottom( ); - vtkCoordinate* coord = this->TextActor->GetPositionCoordinate( ); + vtkCoordinate* coord = this->m_TextActor->GetPositionCoordinate( ); coord->SetCoordinateSystemToNormalizedViewport( ); - coord->SetValue( 0.01, 0.01 ); + coord->SetValue( 0.01, 0.05 ); + + // Configure style + vtkSmartPointer< TStyle > st = vtkSmartPointer< TStyle >::New( ); + st->SetAssociatedObject( this ); + st->AddMouseMoveCommand( Self::_MouseMoveCommand, this ); + st->AddMouseClickCommand( Self::_MouseClickCommand, this ); + st->AddMouseWheelCommand( Self::_MouseWheelCommand, this ); + st->AddKeyCommand( Self::_KeyCommand, this ); + st->AddEnterCommand( Self::_EnterCommand, this ); + st->AddLeaveCommand( Self::_LeaveCommand, this ); + this->m_Style = st; } // ------------------------------------------------------------------------- -unsigned int cpExtensions::Visualization::ImageSliceActors:: -GetNumberOfImageActors( ) const +void cpExtensions::Visualization::ImageSliceActors:: +SetInputConnection( vtkAlgorithmOutput* aout, int orientation ) { - return( this->ImageActors.size( ) ); + if( aout == NULL ) + return; + this->m_Mapper->SetInputConnection( aout ); + this->_ConfigureInput( orientation ); } // ------------------------------------------------------------------------- -vtkImageActor* cpExtensions::Visualization::ImageSliceActors:: -GetImageActor( unsigned int id ) +void cpExtensions::Visualization::ImageSliceActors:: +SetInputData( vtkImageData* data, int orientation ) { - if( id < this->ImageActors.size( ) ) - return( this->ImageActors[ id ] ); - else - return( NULL ); + if( data == NULL ) + return; + this->m_Mapper->SetInputData( data ); + this->_ConfigureInput( orientation ); } // ------------------------------------------------------------------------- -const vtkImageActor* cpExtensions::Visualization::ImageSliceActors:: -GetImageActor( unsigned int id ) const +vtkImageData* cpExtensions::Visualization::ImageSliceActors:: +GetInputData( ) { - if( id < this->ImageActors.size( ) ) - return( this->ImageActors[ id ] ); + if( this->m_Mapper.GetPointer( ) != NULL ) + return( this->m_Mapper->GetInput( ) ); else return( NULL ); } // ------------------------------------------------------------------------- -vtkTextActor* cpExtensions::Visualization::ImageSliceActors:: -GetTextActor( ) +void cpExtensions::Visualization::ImageSliceActors:: +AddMesh( vtkPolyData* mesh ) { - return( this->TextActor ); + SourceActor< vtkCutter > a; + a.Create( ); + a.Source->SetInputData( mesh ); + a.Source->SetCutFunction( this->m_Mapper->GetSlicePlane( ) ); + a.Source->SetValue( 0, 0 ); + a.Source->GenerateTrianglesOff( ); + a.Source->Update( ); + a.Modified( ); + this->m_Meshes[ mesh ] = a; + this->AddItem( a.Actor ); + + auto ren = this->m_Style->GetCurrentRenderer( ); + if( ren != NULL ) + if( ren->HasViewProp( this->m_Actor ) ) + ren->AddViewProp( a.Actor ); } // ------------------------------------------------------------------------- -const vtkTextActor* cpExtensions::Visualization::ImageSliceActors:: -GetTextActor( ) const +void cpExtensions::Visualization::ImageSliceActors:: +AssociateSlice( Self* slice ) { - return( this->TextActor ); + this->m_AssociatedSlices.push_back( slice ); + this->Modified( ); } // ------------------------------------------------------------------------- -vtkActor* cpExtensions::Visualization::ImageSliceActors:: -GetPlaneActor( ) +void cpExtensions::Visualization::ImageSliceActors:: +CleanAssociatedSlices( ) { - return( this->PlaneActor ); + this->m_AssociatedSlices.clear( ); + this->Modified( ); } // ------------------------------------------------------------------------- -const vtkActor* cpExtensions::Visualization::ImageSliceActors:: -GetPlaneActor( ) const +cpExtensions::Interaction::ImageInteractorStyle* +cpExtensions::Visualization::ImageSliceActors:: +GetStyle( ) { - return( this->PlaneActor ); + return( dynamic_cast< TStyle* >( this->m_Style.GetPointer( ) ) ); } // ------------------------------------------------------------------------- -double* cpExtensions::Visualization::ImageSliceActors:: -GetDisplayBounds( ) const +void cpExtensions::Visualization::ImageSliceActors:: +SetStyle( vtkInteractorStyle* st ) { - if( this->ImageActors.size( ) > 0 ) - return( this->ImageActors[ 0 ]->GetDisplayBounds( ) ); - else - return( NULL ); + this->m_Style = st; + this->Modified( ); } // ------------------------------------------------------------------------- void cpExtensions::Visualization::ImageSliceActors:: -GetDisplayBounds( double bounds[ 6 ] ) const +PushInto( vtkRenderer* ren ) { - if( this->ImageActors.size( ) > 0 ) - this->ImageActors[ 0 ]->GetDisplayBounds( bounds ); + this->InitTraversal( ); + vtkProp* prop; + while( ( prop = this->GetNextProp( ) ) != NULL ) + ren->AddViewProp( prop ); } // ------------------------------------------------------------------------- -int cpExtensions::Visualization::ImageSliceActors:: -GetAxis( ) const +void cpExtensions::Visualization::ImageSliceActors:: +PopFrom( vtkRenderer* ren ) { - if( this->SliceMappers.size( ) > 0 ) - return( this->SliceMappers[ 0 ]->GetOrientation( ) ); - else - return( -1 ); + this->InitTraversal( ); + vtkProp* prop; + while( ( prop = this->GetNextProp( ) ) != NULL ) + ren->RemoveViewProp( prop ); } // ------------------------------------------------------------------------- -int cpExtensions::Visualization::ImageSliceActors:: +long cpExtensions::Visualization::ImageSliceActors:: GetSliceNumber( ) const { - if( this->SliceMappers.size( ) > 0 ) - return( this->SliceMappers[ 0 ]->GetSliceNumber( ) ); + if( this->m_Mapper.GetPointer( ) != NULL ) + return( this->m_Mapper->GetSliceNumber( ) ); else return( -1 ); } // ------------------------------------------------------------------------- -int cpExtensions::Visualization::ImageSliceActors:: -GetSliceNumberMinValue( ) const +void cpExtensions::Visualization::ImageSliceActors:: +SetSliceNumber( long slice ) { - if( this->SliceMappers.size( ) > 0 ) - return( this->SliceMappers[ 0 ]->GetSliceNumberMinValue( ) ); - else - return( -1 ); + if( this->m_Mapper.GetPointer( ) == NULL ) + return; + + // Get orientation + int a = this->m_Mapper->GetOrientation( ); + + // Check extent + int ext[ 6 ]; + this->m_Mapper->GetInput( )->GetExtent( ext ); + long rs = slice; + rs = ( rs > ext[ a << 1 ] )? rs: ext[ a << 1 ]; + rs = ( rs < ext[ ( a << 1 ) + 1 ] )? rs: ext[ ( a << 1 ) + 1 ]; + + // Update pipeline + this->m_Mapper->SetSliceNumber( rs ); + + // Update display extent (this isn't done automatically) + ext[ a << 1 ] = ext[ ( a << 1 ) + 1 ] = rs; + this->m_Actor->SetDisplayExtent( ext ); + + this->m_Mapper->Modified( ); + this->m_Actor->Modified( ); + this->_ConfigureCursor( ); + this->Modified( ); } // ------------------------------------------------------------------------- -int cpExtensions::Visualization::ImageSliceActors:: -GetSliceNumberMaxValue( ) const +void cpExtensions::Visualization::ImageSliceActors:: +ShowPixelText( double* pos ) { - if( this->SliceMappers.size( ) > 0 ) - return( this->SliceMappers[ 0 ]->GetSliceNumberMaxValue( ) ); + if( this->m_Mapper.GetPointer( ) != NULL ) + { + char axis; + int axId = this->m_Mapper->GetOrientation( ); + if ( axId == 0 ) axis = 'X'; + else if( axId == 1 ) axis = 'Y'; + else if( axId == 2 ) axis = 'Z'; + long slice = this->GetSliceNumber( ); + + vtkImageData* image = this->m_Actor->GetInput( ); + int ijk[ 3 ]; + double pcoords[ 3 ]; + image->ComputeStructuredCoordinates( pos, ijk, pcoords ); + ijk[ axId ] = slice; + + int ext[ 6 ]; + image->GetExtent( ext ); + if( + ext[ 0 ] <= ijk[ 0 ] && ijk[ 0 ] <= ext[ 1 ] && + ext[ 2 ] <= ijk[ 1 ] && ijk[ 1 ] <= ext[ 3 ] && + ext[ 4 ] <= ijk[ 2 ] && ijk[ 2 ] <= ext[ 5 ] + ) + { + int nScl = image->GetNumberOfScalarComponents( ); + std::stringstream str; + str + << "[" << ijk[ 0 ] + << "," << ijk[ 1 ] + << "," << ijk[ 2 ] << "]=("; + str << + image->GetScalarComponentAsFloat( ijk[ 0 ], ijk[ 1 ], ijk[ 2 ], 0 ); + for( int n = 1; n < nScl; ++n ) + str + << " " + << image->GetScalarComponentAsFloat( + ijk[ 0 ], ijk[ 1 ], ijk[ 2 ], n + ); + str << ")"; + +#if defined(WIN32) + sprintf_s( + this->m_TextBuffer, MAX_TEXT_BUFFER, "Axis: %c (%d)\nPixel %s", + axis, slice, str.str( ).c_str( ) + ); +#else // defined(WIN32) + std::sprintf( + this->m_TextBuffer, "Axis: %c (%ld)\nPixel %s", + axis, slice, str.str( ).c_str( ) + ); +#endif // defined(WIN32) + } // fi + } else - return( -1 ); + this->m_TextBuffer[ 0 ] = '\0'; + this->m_TextActor->SetInput( this->m_TextBuffer ); + this->m_TextActor->Modified( ); + this->Modified( ); } // ------------------------------------------------------------------------- void cpExtensions::Visualization::ImageSliceActors:: -SetSliceNumber( const int& slice ) +GetScalarRange( double r[ 2 ] ) const { - unsigned int nImages = this->SliceMappers.size( ); - if( nImages == 0 ) - return; + r[ 0 ] = this->m_ScalarRange[ 0 ]; + r[ 1 ] = this->m_ScalarRange[ 1 ]; +} + +// ------------------------------------------------------------------------- +void cpExtensions::Visualization::ImageSliceActors:: +SetScalarRange( const double& a, const double& b ) +{ + this->m_ManualScalarRange = true; + this->m_ScalarRange[ 0 ] = a; + this->m_ScalarRange[ 1 ] = b; + this->SetWindowLevel( 1, 0.5 ); +} - // Change visualization extent - for( unsigned int i = 0; i < nImages; ++i ) +// ------------------------------------------------------------------------- +void cpExtensions::Visualization::ImageSliceActors:: +UnsetScalarRange( ) +{ + auto image = this->m_Actor->GetInput( ); + if( image != NULL ) { - this->SliceMappers[ i ]->SetSliceNumber( slice ); - this->SliceMappers[ i ]->Modified( ); - this->ImageActors[ i ]->Modified( ); - this->SliceMappers[ i ]->Update( ); + double r[ 2 ]; + image->GetScalarRange( r ); + this->SetScalarRange( r[ 0 ], r[ 1 ] ); - } // rof + } // fi + this->m_ManualScalarRange = false; +} - // Compute plane - vtkAlgorithm* algo = this->SliceMappers[ 0 ]->GetInputAlgorithm( ); - vtkInformation* info = algo->GetOutputInformation( 0 ); - int ext[ 6 ]; - double ori[ 3 ], spac[ 3 ], pos[ 3 ]; - info->Get( vtkStreamingDemandDrivenPipeline::WHOLE_EXTENT( ), ext ); - info->Get( vtkDataObject::ORIGIN( ), ori ); - info->Get( vtkDataObject::SPACING( ), spac ); - this->SliceMappers[ 0 ]->GetSlicePlane( )->GetOrigin( pos ); - - // Prevent obscuring voxels by offsetting the plane geometry - double xbnds[ ] = - { - ori[ 0 ] + ( spac[ 0 ] * double( ext[ 0 ] ) ), - ori[ 0 ] + ( spac[ 0 ] * double( ext[ 1 ] ) ) - }; - double ybnds[ ] = - { - ori[ 1 ] + ( spac[ 1 ] * double( ext[ 2 ] ) ), - ori[ 1 ] + ( spac[ 1 ] * double( ext[ 3 ] ) ) - }; - double zbnds[ ] = - { - ori[ 2 ] + ( spac[ 2 ] * double( ext[ 4 ] ) ), - ori[ 2 ] + ( spac[ 2 ] * double( ext[ 5 ] ) ) - }; +// ------------------------------------------------------------------------- +void cpExtensions::Visualization::ImageSliceActors:: +SetWindowLevel( const double& w, const double& l ) +{ + static const double _0 = double( 0 ); + static const double _1 = double( 1 ); - if( spac[ 0 ] < double( 0 ) ) - { - double t = xbnds[ 0 ]; - xbnds[ 0 ] = xbnds[ 1 ]; - xbnds[ 1 ] = t; + double rw = ( w < _0 )? _0: ( ( w > _1 )? _1: w ); + double rl = ( l < _0 )? _0: ( ( l > _1 )? _1: l ); - } // fi - if( spac[ 1 ] < double( 0 ) ) - { - double t = ybnds[ 0 ]; - ybnds[ 0 ] = ybnds[ 1 ]; - ybnds[ 1 ] = t; + double d = this->m_ScalarRange[ 1 ] - this->m_ScalarRange[ 0 ]; + rw *= d; + rl *= d; + rl += this->m_ScalarRange[ 0 ]; - } // fi - if( spac[ 2 ] < double( 0 ) ) + if( this->m_Actor.GetPointer( ) != NULL ) { - double t = zbnds[ 0 ]; - zbnds[ 0 ] = zbnds[ 1 ]; - zbnds[ 1 ] = t; + this->m_Actor->GetProperty( )->SetColorWindow( rw ); + this->m_Actor->GetProperty( )->SetColorLevel( rl ); + this->Modified( ); } // fi +} - int axis = this->SliceMappers[ 0 ]->GetOrientation( ); - this->PlaneActor->GetProperty( )->SetRepresentationToWireframe( ); - this->PlaneActor->GetProperty( )->SetLineWidth( 2 ); - vtkPoints* plane_points = this->PlaneSource->GetPoints( ); - if( axis == 0 ) // YZ, x-normal +// ------------------------------------------------------------------------- +double cpExtensions::Visualization::ImageSliceActors:: +GetWindow( ) const +{ + if( this->m_Actor.GetPointer( ) != NULL ) + return( this->m_Actor->GetProperty( )->GetColorWindow( ) ); + else + return( double( 0 ) ); +} + +// ------------------------------------------------------------------------- +double cpExtensions::Visualization::ImageSliceActors:: +GetLevel( ) const +{ + if( this->m_Actor.GetPointer( ) != NULL ) + return( this->m_Actor->GetProperty( )->GetColorLevel( ) ); + else + return( double( 0 ) ); +} + +// ------------------------------------------------------------------------- +void cpExtensions::Visualization::ImageSliceActors:: +Render( ) +{ + auto ren = this->m_Style->GetCurrentRenderer( ); + if( ren == NULL ) + return; + auto win = ren->GetRenderWindow( ); + if( win == NULL ) + return; + win->Render( ); +} + +// ------------------------------------------------------------------------- +cpExtensions::Visualization::ImageSliceActors:: +ImageSliceActors( ) + : Superclass( ), + m_ManualScalarRange( false ) +{ + this->Clear( ); +} + +// ------------------------------------------------------------------------- +cpExtensions::Visualization::ImageSliceActors:: +~ImageSliceActors( ) +{ +} + +// ------------------------------------------------------------------------- +void cpExtensions::Visualization::ImageSliceActors:: +_ConfigureInput( int orientation ) +{ + // Prepare main image + this->m_Mapper->SetOrientation( orientation ); + this->m_Mapper->Update( ); + this->m_Actor->GetProperty( )->SetOpacity( 0.8 ); + this->m_Actor->Modified( ); + auto st = this->GetStyle( ); + if( st != NULL ) + st->AssociateImageActor( this->m_Actor ); + this->AddItem( this->m_Actor ); + + // Prepare main cursor + this->_ConfigureCursor( ); + this->AddItem( this->m_Cursor.Actor ); + + // Prepare window/level + if( !( this->m_ManualScalarRange ) ) + this->UnsetScalarRange( ); + + // Add all other actors + this->AddItem( this->m_TextActor ); + + // Put it in the first slice + auto image = this->m_Mapper->GetInput( ); + int ext[ 6 ]; + image->GetExtent( ext ); + this->SetSliceNumber( ext[ orientation << 1 ] ); + + /* + if( this->m_Style.GetPointer( ) != NULL ) + this->AddItem( this->m_Actor ); + + 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:: +_ConfigureCursor( ) +{ + int a = this->m_Mapper->GetOrientation( ); + double bounds[ 6 ]; + this->m_Actor->GetBounds( bounds ); + bounds[ a << 1 ] += 1e-5; + bounds[ ( a << 1 ) + 1 ] += 1e-5; + + this->m_Cursor.Source->SetModelBounds( bounds ); + this->m_Cursor.Source->AllOn( ); + this->m_Cursor.Actor->GetProperty( )->SetOpacity( 1 ); + this->m_Cursor.Actor->GetProperty( )->SetLineWidth( 3 ); + this->m_Cursor.Modified( ); +} + +// ------------------------------------------------------------------------- +void cpExtensions::Visualization::ImageSliceActors:: +_MouseMoveCommand( + void* data, + const TStyle::ButtonID& btn, + int* idx, double* pos, + bool alt, bool ctr, bool sft + ) +{ + ImageSliceActors* self = reinterpret_cast< ImageSliceActors* >( data ); + if( self == NULL ) + return; + + if( btn == TStyle::ButtonID_None ) { - plane_points->SetPoint( 0, pos[ 0 ], ybnds[ 0 ], zbnds[ 0 ] ); - plane_points->SetPoint( 1, pos[ 0 ], ybnds[ 1 ], zbnds[ 0 ] ); - plane_points->SetPoint( 2, pos[ 0 ], ybnds[ 1 ], zbnds[ 1 ] ); - plane_points->SetPoint( 3, pos[ 0 ], ybnds[ 0 ], zbnds[ 1 ] ); - this->PlaneActor->GetProperty( )->SetColor( 1, 0, 0 ); + self->ShowPixelText( pos ); + self->m_Cursor.Source->SetFocalPoint( pos ); + self->m_Cursor.Source->Modified( ); + self->m_Cursor.Modified( ); + self->Render( ); } - else if( axis == 1 ) // ZX, y-normal + else if( btn == TStyle::ButtonID_Right ) { - plane_points->SetPoint( 0, xbnds[ 0 ], pos[ 1 ], zbnds[ 0 ] ); - plane_points->SetPoint( 1, xbnds[ 0 ], pos[ 1 ], zbnds[ 1 ] ); - plane_points->SetPoint( 2, xbnds[ 1 ], pos[ 1 ], zbnds[ 1 ] ); - plane_points->SetPoint( 3, xbnds[ 1 ], pos[ 1 ], zbnds[ 0 ] ); - this->PlaneActor->GetProperty( )->SetColor( 0, 1, 0 ); - } - else // XY, z-normal + if( !alt && !ctr && sft ) + { + int a0 = self->m_Mapper->GetOrientation( ); + int a1 = ( a0 + 1 ) % 3; + int a2 = ( a0 + 2 ) % 3; + double dx = pos[ a1 ] - self->m_StartMouseEvent[ a1 ]; + double dy = pos[ a2 ] - self->m_StartMouseEvent[ a2 ]; + double bounds[ 6 ]; + self->m_Actor->GetBounds( bounds ); + dx /= bounds[ ( a1 << 1 ) + 1 ] - bounds[ a1 << 1 ]; + dy /= bounds[ ( a2 << 1 ) + 1 ] - bounds[ a2 << 1 ]; + double w = self->m_StartWindow + dx; + double l = self->m_StartLevel + dy; + + self->SetWindowLevel( w, l ); + self->Render( ); + for( + auto aIt = self->m_AssociatedSlices.begin( ); + aIt != self->m_AssociatedSlices.end( ); + ++aIt + ) + { + ( *aIt )->SetWindowLevel( w, l ); + ( *aIt )->Render( ); + + } // rof + + } // fi + + } // fi +} + +// ------------------------------------------------------------------------- +void cpExtensions::Visualization::ImageSliceActors:: +_MouseClickCommand( + void* data, + const TStyle::ButtonID& btn, + int* idx, double* pos, + bool alt, bool ctr, bool sft + ) +{ + ImageSliceActors* self = reinterpret_cast< ImageSliceActors* >( data ); + if( self == NULL ) + return; + + self->m_StartMouseEvent[ 0 ] = pos[ 0 ]; + self->m_StartMouseEvent[ 1 ] = pos[ 1 ]; + self->m_StartMouseEvent[ 2 ] = pos[ 2 ]; + + if( btn == TStyle::ButtonID_Right ) { - plane_points->SetPoint( 0, xbnds[ 0 ], ybnds[ 0 ], pos[ 2 ] ); - plane_points->SetPoint( 1, xbnds[ 1 ], ybnds[ 0 ], pos[ 2 ] ); - plane_points->SetPoint( 2, xbnds[ 1 ], ybnds[ 1 ], pos[ 2 ] ); - plane_points->SetPoint( 3, xbnds[ 0 ], ybnds[ 1 ], pos[ 2 ] ); - this->PlaneActor->GetProperty( )->SetColor( 0, 0, 1 ); + if( !alt && !ctr && sft ) + { + double r[ 2 ]; + self->GetScalarRange( r ); + double d = r[ 1 ] - r[ 0 ]; + self->m_StartWindow = self->GetWindow( ) / d; + self->m_StartLevel = ( self->GetLevel( ) - r[ 0 ] ) / d; + + } // fi } // fi - this->PlaneSource->Modified( ); - this->PlaneMapper->Modified( ); - this->PlaneActor->Modified( ); - this->Modified( ); } // ------------------------------------------------------------------------- void cpExtensions::Visualization::ImageSliceActors:: -UpdateText( ) +_MouseWheelCommand( + void* data, + const int& dir, bool alt, bool ctr, bool sft + ) { - if( this->SliceMappers.size( ) > 0 ) + ImageSliceActors* self = reinterpret_cast< ImageSliceActors* >( data ); + if( self == NULL ) + return; + auto ren = self->m_Style->GetCurrentRenderer( ); + if( ren == NULL ) + return; + auto cam = ren->GetActiveCamera( ); + if( cam == NULL ) + return; + auto win = ren->GetRenderWindow( ); + if( win == NULL ) + return; + + // Get previous values + int a = self->m_Mapper->GetOrientation( ); + double bounds[ 6 ]; + self->m_Actor->GetBounds( bounds ); + double prev_pos = bounds[ a << 1 ]; + + // Move slice + long slice = self->GetSliceNumber( ) + dir * ( ( sft )? 10: 1 ); + self->SetSliceNumber( slice ); + self->m_Actor->GetBounds( bounds ); + + // Reset camera + double cam_pos[ 3 ]; + cam->GetPosition( cam_pos ); + cam_pos[ a ] += bounds[ a << 1 ] - prev_pos; + cam->SetPosition( cam_pos ); + + // Update text + self->ShowPixelText( self->m_Cursor.Source->GetFocalPoint( ) ); + + // Update all + win->Render( ); + for( + auto aIt = self->m_AssociatedSlices.begin( ); + aIt != self->m_AssociatedSlices.end( ); + ++aIt + ) { - char axis; - int axId = this->SliceMappers[ 0 ]->GetOrientation( ); - if ( axId == 0 ) axis = 'X'; - else if( axId == 1 ) axis = 'Y'; - else if( axId == 2 ) axis = 'Z'; + if( ( *aIt )->m_Mapper->GetOrientation( ) == a ) + { + ( *aIt )->SetSliceNumber( slice ); + ( *aIt )->Render( ); - std::sprintf( - this->TextBuffer, "Axis: %c (%d)", - axis, this->SliceMappers[ 0 ]->GetSliceNumber( ) - ); - } - else - this->TextBuffer[ 0 ] = '\0'; - this->TextActor->SetInput( this->TextBuffer ); - this->TextActor->Modified( ); - this->Modified( ); + } // fi + + } // fi } // ------------------------------------------------------------------------- -cpExtensions::Visualization::ImageSliceActors:: -ImageSliceActors( ) - : Superclass( ) +void cpExtensions::Visualization::ImageSliceActors:: +_KeyCommand( + void* data, + const char& key + ) { - this->Clear( ); } // ------------------------------------------------------------------------- -cpExtensions::Visualization::ImageSliceActors:: -~ImageSliceActors( ) +void cpExtensions::Visualization::ImageSliceActors:: +_EnterCommand( void* data ) +{ +} + +// ------------------------------------------------------------------------- +void cpExtensions::Visualization::ImageSliceActors:: +_LeaveCommand( void* data ) { }