[vtkusers] How can i calculate my new extent in Z after vtkImageReslice

[marcelo oliveira]

Hi user !!!!

I´m very in doubt !!!!

I have a stack of axial CT image

x = 512
y = 512
Z = 81

space x = 1
space y = 1
space z = 13.65

How can i caculate my new stack in Z when i Rotate Y
by 64 degrees

i try 

double ZExtent = 
80.00/sin(1.11);                

                  //90-64 degrees 
spaceZ = 13.65 / sin(1.54-1.11);
spaceZ = sqrt( pow( (13.65 * 80), 2 ) + 
pow( spaceZ, 2 ) ) / ZExtent ;

then

reslice->SetOutputOrigin( 0.0, 0.000, 0.00 );
reslice->SetOutputExtent( -2000.00, 2000.00, -2000.00,
2000.00, 000.00, ZExtent);
reslice->SetOutputSpacing( 1.00, 1.00, spaceZ );

reslice->SetResliceAxesDirectionCosines(  
cos( -rotacaoImagem ), 0, sin(-rotacaoImagem),
0, 1, 0,
-sin(-rotacaoImagem), 0, cos(-rotacaoImagem)
);

Is right ????


Tanks a lot for help me !!!!!!1

 


_______________________________________________________________________
Desafio AntiZona: participe do jogo de perguntas e respostas que vai
dar um Renault Clio, computadores, câmeras digitais, videogames e muito
mais! www.cade.com.br/antizona