Hi guys,<br><br>I'm planing to use vtkImageReslice to get a slice from an 3D vtkImageData in two different ways:<br><br>- the first one based on rotation matrices<br>i just used the rotation angles as the directions cosines and it works as expected. In this implementation i'm using<br>
<br>SetResliceAxesDirectionCosines(...); //rotation matrix<br>
SetResliceAxesOrigin(...); // rotation point<br><br>
- the second one is based on the normal to a plane<br>
a have two points (in world coordinates) and i want to get the a plane that is orthogonal to that vector. I've made some calculations and I got the plane, the normal and two other points inside that plane. As far a i know, i must use the same methods to get the slice<br>
<br>slicer->SetResliceAxesDirectionCosines(...);<br>slicer->SetResliceAxesOrigin(...);<br><br>plus the following:<br><br>slicer->SetOutputDimensionality(2);<br>slicer->SetOutputSpacing(...);<br>slicer->SetOutputExtent(...);<br>
<br>The result slice seems to have the correct pixels, however the data is skewed 45 degrees!! <a href="http://tinypic.com/r/qpj6kz/5">http://tinypic.com/r/qpj6kz/5</a>.<br><br>So my questions are:<br>- what are the actual parameters of 'SetResliceAxesDirectionCosines'?<br>
<br>- this post: <a href="http://www.vtk.org/pipermail/vtkusers/2004-May/074177.html">http://www.vtk.org/pipermail/vtkusers/2004-May/074177.html</a>, depicts that the best example to look at is 'vtkImagePlaneWidget', however after making those changes i'm getting a completely black image. So i think the problem must be the parameters to SetResliceAxesDirectionCosines that I'm invoquing in this way:<br>
<br>SetResliceAxesDirectionCosines (<br> v1[0], v1[1], v1[2], //based on the first point in the plane<br> v2[0], v2[1], v2[2], //based on the second point in the plane<br> normal[0], normal[1], normal[2]<br>
);<br><a href="http://www.vtk.org/pipermail/vtkusers/2004-May/074177.html" target="_blank"></a><a href="http://web.archiveorange.com/archive/v/50zHjsWfslvTDBylAOk7" target="_blank"></a><br>Apart from normalization, what else must be done to the vectors that represent the plane? do they must the something in particular, like the corners? or can be just anywhere in the plane?<br>
<br>I'am also aware that for images with different scales a workaround could be using a ImageChangeInformation as a wrapper:<br><br>vtkSmartPointer<itkImageChangeInformation> change = vtkImageChangeInformation::New();<br>
change->SetInput(reader->GetOutput());<br><div class="gmail_extra">change->SetOutputSpacing(1,1,1);<br>change->Update();<br><br>but i still get the same results. Thanks for your help.<br></div><div class="gmail_extra">
<br>Best regards,<br><br>Fabian</div><br>