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¶
Danger
The generation of tar files has not been ported to the new VTKExamples. Some tarballs may be missing or out-of-date.
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 PLEASE NOTE: Be sure to add the VTK bin directory to your path. This will resolve the VTK dll's at run time.