private void RotateImage(RotationOperation operation)
{
vtkImageReslice reslice = vtkImageReslice.New();
vtkTransform transform = vtkTransform.New();
transform.PostMultiply();
double[] center = { 75, 100, 0 };
transform.Translate(-center[0], -center[1], -center[2]);
if (operation == RotationOperation.Forward)
{
transform.RotateZ(90);
}
else if (operation == RotationOperation.Back)
{
transform.RotateZ(-90);
}
transform.Translate(+center[0], +center[1], +center[2]);
transform.Update();
reslice.SetInput(_viewer.GetInput());
reslice.SetResliceTransform(transform);
reslice.Update();
_viewer.SetInput(reslice.GetOutput());
UpdateViewer();
}
/// <summary>
/// Update the 2D visualization when the plane moved.
/// </summary>
/// <param name="plane">PlaneWidget which changed the coordinates.</param>
public void PlaneMoved(vtkImagePlaneWidget plane)
{
vtkImageReslice reslice = vtkImageReslice.New();
vtkTransform transform = vtkTransform.New();
transform.PostMultiply();
//TODO wyznaczenie centrum okna
double[] center = { 75, 100, 0 };
transform.Translate(-center[0], -center[1], -center[2]);
transform.RotateZ(_orientation);
transform.Translate(+center[0], +center[1], +center[2]);
transform.Update();
reslice.SetInput(plane.GetResliceOutput());
reslice.SetResliceTransform(transform);
reslice.Update();
//
vtkImageData viewerInputData;
viewerInputData = reslice.GetOutput();
_lastData = viewerInputData;
_viewer.SetInput(viewerInputData);
//
//_viewer.SetInput(reslice.GetOutput());
UpdateViewer();
}
/// <summary>
/// Rotates volume along given axis.
/// </summary>
/// <param name="volume"></param>
/// <param name="angles"></param>
public static void RotateData(ref vtkImageData volume, double[] angles)
{
// Rotation along image center
double[] center = volume.GetExtent().Divide(2);
// Dimensions of rotated image
int[] outExtent = volume.GetExtent().Multiply(1.1).Round().ToInt32();
// Rotation parameters
var rotate = new vtkTransform();
rotate.Translate(center[1], center[3], center[5]);
rotate.RotateX(angles[0]);
rotate.RotateY(angles[1]);
rotate.RotateZ(angles[2]); // z angle should be 0
rotate.Translate(-center[1], -center[3], -center[5]);
// Perform rotation
var slice = new vtkImageReslice();
slice.SetInput(volume);
slice.SetResliceTransform(rotate);
slice.SetInterpolationModeToCubic();
slice.SetOutputSpacing(volume.GetSpacing()[0], volume.GetSpacing()[1], volume.GetSpacing()[2]);
slice.SetOutputOrigin(volume.GetOrigin()[0], volume.GetOrigin()[1], volume.GetOrigin()[2]);
slice.SetOutputExtent(outExtent[0], outExtent[1], outExtent[2], outExtent[3], outExtent[4], outExtent[5]);
}
//TODO
public void DrawingModeRepaint()
{
if (_lastData == null)
{
return;
}
vtkImageReslice reslice = vtkImageReslice.New();
vtkTransform transform = vtkTransform.New();
transform.PostMultiply();
//TODO wyznaczenie centrum okna
double[] center = { 75, 100, 0 };
transform.Translate(-center[0], -center[1], -center[2]);
transform.RotateZ(_orientation);
transform.Translate(+center[0], +center[1], +center[2]);
transform.Update();
reslice.SetInput(_lastData);
reslice.SetResliceTransform(transform);
reslice.Update();
vtkImageData viewerInputData;
if (!DrawingEnabled)
{
viewerInputData = reslice.GetOutput();
}
else
{
var resliceOutput = reslice.GetOutput();
var sizeArray = resliceOutput.GetDimensions();
vtkImageCanvasSource2D imageCanvas = vtkImageCanvasSource2D.New();
imageCanvas.SetScalarTypeToUnsignedChar();
imageCanvas.SetExtent(0, sizeArray[0], 0, sizeArray[1], 0, 0);
imageCanvas.DrawImage(0, 0, resliceOutput);
DrawingUtils.Draw(imageCanvas, DrawingInfo);
imageCanvas.Update();
viewerInputData = imageCanvas.GetOutput();
}
_viewer.SetInput(viewerInputData);
UpdateViewer();
}
public void SetRotation()
{
vtkTransform.RotateX(rotateX);
vtkTransform.RotateY(rotateY);
vtkTransform.RotateZ(rotateZ);
}