# ImageRotate

VTKExamples/Cxx/Images/ImageRotate

### Code

ImageRotate.cxx

    #include <vtkSmartPointer.h>
    #include <vtkImageReader2Factory.h>
    #include <vtkImageReader2.h>
    #include <vtkImageData.h>
    #include <vtkRenderWindow.h>
    #include <vtkRenderWindowInteractor.h>
    #include <vtkInteractorStyleImage.h>
    #include <vtkRenderer.h>
    #include <vtkImageActor.h>
    #include <vtkImageReslice.h>
    #include <vtkImageViewer2.h>
    #include <vtkTransform.h>

    int main(int argc, char *argv[])
    {
      // Verify command line arguments
      if(argc < 2)
      {
        std::cout << "Usage: " << argv[0]
                  << " InputFilename" << std::endl;
        return EXIT_FAILURE;
      }

      double angle = 45;
      if (argc > 2)
      {
        angle = atof(argv[2]);
      }

      // Read file
      vtkSmartPointer<vtkImageReader2Factory> readerFactory =
        vtkSmartPointer<vtkImageReader2Factory>::New();

      vtkImageReader2 *reader =
        readerFactory->CreateImageReader2(argv[1]);
      reader->SetFileName(argv[1]);
      reader->Update();
      double bounds[6];
      reader->GetOutput()->GetBounds(bounds);

      // Rotate about the center of the image
      vtkSmartPointer<vtkTransform> transform =
        vtkSmartPointer<vtkTransform>::New();

      // Compute the center of the image
      double center[3];
      center[0] = (bounds[1] + bounds[0]) / 2.0;
      center[1] = (bounds[3] + bounds[2]) / 2.0;
      center[2] = (bounds[5] + bounds[4]) / 2.0;

      // Rotate about the center
      transform->Translate(center[0], center[1], center[2]);
      transform->RotateWXYZ(angle, 0, 0, 1);
      transform->Translate(-center[0], -center[1], -center[2]);

      // Reslice does all of the work
      vtkSmartPointer<vtkImageReslice> reslice =
        vtkSmartPointer<vtkImageReslice>::New();
      reslice->SetInputConnection(reader->GetOutputPort());
      reslice->SetResliceTransform(transform);
      reslice->SetInterpolationModeToCubic();
      reslice->SetOutputSpacing(
        reader->GetOutput()->GetSpacing()[0],
        reader->GetOutput()->GetSpacing()[1],
        reader->GetOutput()->GetSpacing()[2]);
      reslice->SetOutputOrigin(
        reader->GetOutput()->GetOrigin()[0],
        reader->GetOutput()->GetOrigin()[1],
        reader->GetOutput()->GetOrigin()[2]);
      reslice->SetOutputExtent(reader->GetOutput()->GetExtent()); // Use a larger extent than the original image's to prevent clipping

      // Visualize
      vtkSmartPointer<vtkImageViewer2> imageViewer =
        vtkSmartPointer<vtkImageViewer2>::New();
      imageViewer->SetInputConnection(reslice->GetOutputPort());
      vtkSmartPointer<vtkRenderWindowInteractor> renderWindowInteractor =
        vtkSmartPointer<vtkRenderWindowInteractor>::New();
      imageViewer->SetupInteractor(renderWindowInteractor);
      imageViewer->Render();
      imageViewer->GetRenderer()->ResetCamera();
      imageViewer->Render();

      renderWindowInteractor->Start();

      reader->Delete();

      return EXIT_SUCCESS;
    }

### CMakeLists.txt

    cmake_minimum_required(VERSION 2.8)

    PROJECT(ImageRotate)

    find_package(VTK REQUIRED)
    include(${VTK_USE_FILE})

    add_executable(ImageRotate MACOSX_BUNDLE ImageRotate.cxx )

    target_link_libraries(ImageRotate ${VTK_LIBRARIES})

### Download and Build ImageRotate

Click here to download ImageRotate and its CMakeLists.txt file. Once the tarball ImageRotate.tar has been downloaded and extracted,

    cd ImageRotate/build

If VTK is installed:

    cmake ..

If VTK is not installed but compiled on your system, you will need to specify the path to your VTK build:

    cmake -DVTK_DIR:PATH=/home/me/vtk_build ..

Build the project:

    make

and run it:

    ./ImageRotate

WINDOWS USERS

Be sure to add the VTK bin directory to your path. This will resolve the VTK dll's at run time.