ITK/Examples/SimpleOperations/TranslationTransform: Difference between revisions

From KitwarePublic
< ITK‎ | Examples
Jump to navigationJump to search
No edit summary
Line 1: Line 1:
This example currently freezes. The idea is to display an image along with a sphere at (0,0). Then display the translated image with the sphere still at (0,0). This should be enough to see the translation that occured.
==TranslationTransform.cxx==
 
==TranslateAnImage.cxx==
<source lang="cpp">
<source lang="cpp">
#include "itkImage.h"
#include "itkImage.h"
Line 8: Line 6:
#include "itkNormalizeImageFilter.h"
#include "itkNormalizeImageFilter.h"
#include "itkResampleImageFilter.h"
#include "itkResampleImageFilter.h"
#include "itkImageFileWriter.h"


#include <itkImageToVTKImageFilter.h>
#include "QuickView.h"


#include "vtkImageViewer.h"
typedef itk::Image<unsigned char, 2>  ImageType;
#include "vtkRenderWindowInteractor.h"
 
#include "vtkSmartPointer.h"
void CreateImage(ImageType::Pointer image);
#include "vtkImageActor.h"
#include "vtkInteractorStyleImage.h"
#include "vtkRenderer.h"
#include "vtkPolyDataMapper.h"
#include "vtkActor.h"
#include "vtkProperty.h"
#include "vtkSphereSource.h"


int main(int argc, char *argv[])
int main(int argc, char *argv[])
{
{
   if(argc < 2)
   ImageType::Pointer image = ImageType::New();
    {
   CreateImage(image);
    std::cerr << "Required: filename" << std::endl;
    return EXIT_FAILURE;
    }
 
  typedef itk::Image<unsigned char, 2>  ImageType;
 
  typedef itk::ImageFileReader<ImageType> ReaderType;
  ReaderType::Pointer reader = ReaderType::New();
  reader->SetFileName(argv[1]);
 
  typedef itk::ImageToVTKImageFilter<ImageType> ConnectorType;
  ConnectorType::Pointer originalConnector = ConnectorType::New();
   originalConnector->SetInput(reader->GetOutput());
  vtkSmartPointer<vtkImageActor> originalActor =
    vtkSmartPointer<vtkImageActor>::New();
  originalActor->SetInput(originalConnector->GetOutput());


   itk::TranslationTransform<double,2>::Pointer transform =
   typedef itk::TranslationTransform<double,2> TranslationTransformType;
     itk::TranslationTransform<double,2>::New();
  TranslationTransformType::Pointer transform =
   itk::TranslationTransform<double,2>::OutputVectorType translation;
     TranslationTransformType::New();
   TranslationTransformType::OutputVectorType translation;
   translation[0] = 10;
   translation[0] = 10;
   translation[0] = 20;
   translation[1] = 20;
   transform->Translate(translation);
   transform->Translate(translation);


   vtkSmartPointer<vtkSphereSource> sphereSource =
   typedef itk::ResampleImageFilter<ImageType, ImageType> ResampleImageFilterType;
    vtkSmartPointer<vtkSphereSource>::New();
  ResampleImageFilterType::Pointer resampleFilter = ResampleImageFilterType::New();
   sphereSource->SetRadius(10);
   resampleFilter->SetTransform(transform.GetPointer());
  sphereSource->SetCenter(0, 0, 0);
   resampleFilter->SetInput(image);
   sphereSource->Update();


   vtkSmartPointer<vtkPolyDataMapper> mapper =
   /*
    vtkSmartPointer<vtkPolyDataMapper>::New();
  // These are the defaults
   mapper->SetInputConnection(sphereSource->GetOutputPort());
   double spacing[ 2 ];
   vtkSmartPointer<vtkActor> sphereActor =
   spacing[0] = 1.0;
    vtkSmartPointer<vtkActor>::New();
   spacing[1] = 1.0;
  sphereActor->SetMapper(mapper);
   resampleFilter->SetOutputSpacing( spacing );
  sphereActor->GetProperty()->SetColor(1,0,0);
    
  typedef itk::ResampleImageFilter< ImageType, ImageType > ResampleFilterType;
  ResampleFilterType::Pointer resampleFilter = ResampleFilterType::New();
  resampleFilter->SetInput(reader->GetOutput());
   resampleFilter->SetTransform(transform);


   ConnectorType::Pointer resampledConnector = ConnectorType::New();
   double origin[ 2 ];
   resampledConnector->SetInput(resampleFilter->GetOutput());
  origin[0] = 0.0;
  origin[1] = 0.0;
   resampleFilter->SetOutputOrigin( origin );
  */


   vtkSmartPointer<vtkImageActor> resampledActor =
   // Without this, the program crashes
    vtkSmartPointer<vtkImageActor>::New();
  ImageType::SizeType  size = image->GetLargestPossibleRegion().GetSize();
   resampledActor->SetInput(resampledConnector->GetOutput());
   resampleFilter->SetSize( size );
 
  resampleFilter->Update();


   // There will be one render window
   QuickView viewer;
  vtkSmartPointer<vtkRenderWindow> renderWindow =
   viewer.AddImage(image.GetPointer());
    vtkSmartPointer<vtkRenderWindow>::New();
   viewer.AddImage(resampleFilter->GetOutput());
   renderWindow->SetSize(600, 300);
   viewer.Visualize();
 
    
  vtkSmartPointer<vtkRenderWindowInteractor> interactor =
    vtkSmartPointer<vtkRenderWindowInteractor>::New();
  interactor->SetRenderWindow(renderWindow);
 
  // Define viewport ranges
  // (xmin, ymin, xmax, ymax)
  double leftViewport[4] = {0.0, 0.0, 0.5, 1.0};
  double rightViewport[4] = {0.5, 0.0, 1.0, 1.0};
 
  // Setup both renderers
  vtkSmartPointer<vtkRenderer> leftRenderer =
    vtkSmartPointer<vtkRenderer>::New();
  renderWindow->AddRenderer(leftRenderer);
  leftRenderer->SetViewport(leftViewport);
   leftRenderer->SetBackground(.6, .5, .4);
 
  vtkSmartPointer<vtkRenderer> rightRenderer =
    vtkSmartPointer<vtkRenderer>::New();
  renderWindow->AddRenderer(rightRenderer);
  rightRenderer->SetViewport(rightViewport);
   rightRenderer->SetBackground(.4, .5, .6);
 
  // Add the sphere to the left and the cube to the right
  leftRenderer->AddActor(originalActor);
   leftRenderer->AddActor(sphereActor);
 
  rightRenderer->AddActor(resampledActor);
  rightRenderer->AddActor(sphereActor);
    
    
   leftRenderer->ResetCamera();
   return EXIT_SUCCESS;
  rightRenderer->ResetCamera();
}


   renderWindow->Render();
void CreateImage(ImageType::Pointer image)
{
  ImageType::IndexType start;
   start.Fill(0);


   vtkSmartPointer<vtkInteractorStyleImage> style =
   ImageType::SizeType size;
    vtkSmartPointer<vtkInteractorStyleImage>::New();
  size.Fill(100);


   interactor->SetInteractorStyle(style);
   ImageType::RegionType region(start, size);
  image->SetRegions(region);
  image->Allocate();
  image->FillBuffer(0);


   interactor->Start();
   // Make a square
  for(unsigned int r = 40; r < 60; r++)
    {
    for(unsigned int c = 40; c < 60; c++)
      {
      ImageType::IndexType pixelIndex;
      pixelIndex[0] = r;
      pixelIndex[1] = c;


  return EXIT_SUCCESS;
      image->SetPixel(pixelIndex, 255);
      }
    }
}
}
</source>
</source>
Line 130: Line 92:
cmake_minimum_required(VERSION 2.6)
cmake_minimum_required(VERSION 2.6)


PROJECT(ResampleImageFilter)
PROJECT(TranslationTransform)


include_directories(/home/doriad/src/ITK/Wrapping/WrapITK/ExternalProjects/ItkVtkGlue/src/)
include_directories(/home/doriad/ITKWikiExamples/ItkVtkGlue)


FIND_PACKAGE(VTK REQUIRED)
FIND_PACKAGE(VTK REQUIRED)
Line 140: Line 102:
INCLUDE(${ITK_USE_FILE})
INCLUDE(${ITK_USE_FILE})


ADD_EXECUTABLE(ResampleImageFilter ResampleImageFilter.cxx)
ADD_EXECUTABLE(TranslationTransform TranslationTransform.cxx /home/doriad/ITKWikiExamples/ItkVtkGlue/QuickView.cxx)
TARGET_LINK_LIBRARIES(ResampleImageFilter vtkHybrid
TARGET_LINK_LIBRARIES(TranslationTransform
ITKIO ITKNumerics ITKBasicFilters ITKCommon
vtkHybrid
)
ITKBasicFilters ITKCommon ITKIO)
 
</source>
</source>

Revision as of 17:30, 27 January 2011

TranslationTransform.cxx

<source lang="cpp">

  1. include "itkImage.h"
  2. include "itkTranslationTransform.h"
  3. include "itkImageFileReader.h"
  4. include "itkNormalizeImageFilter.h"
  5. include "itkResampleImageFilter.h"
  6. include "itkImageFileWriter.h"
  1. include "QuickView.h"

typedef itk::Image<unsigned char, 2> ImageType;

void CreateImage(ImageType::Pointer image);

int main(int argc, char *argv[]) { 

 ImageType::Pointer image = ImageType::New();
 CreateImage(image);

 typedef itk::TranslationTransform<double,2> TranslationTransformType;
 TranslationTransformType::Pointer transform =
   TranslationTransformType::New();
 TranslationTransformType::OutputVectorType translation;
 translation[0] = 10;
 translation[1] = 20;
 transform->Translate(translation);

 typedef itk::ResampleImageFilter<ImageType, ImageType> ResampleImageFilterType;
 ResampleImageFilterType::Pointer resampleFilter = ResampleImageFilterType::New();
 resampleFilter->SetTransform(transform.GetPointer());
 resampleFilter->SetInput(image);

 /*
 // These are the defaults
 double spacing[ 2 ];
 spacing[0] = 1.0;
 spacing[1] = 1.0;
 resampleFilter->SetOutputSpacing( spacing );

 double origin[ 2 ];
 origin[0] = 0.0;
 origin[1] = 0.0;
 resampleFilter->SetOutputOrigin( origin );
 */

 // Without this, the program crashes
 ImageType::SizeType   size = image->GetLargestPossibleRegion().GetSize();
 resampleFilter->SetSize( size );
 
 resampleFilter->Update();

 QuickView viewer;
 viewer.AddImage(image.GetPointer());
 viewer.AddImage(resampleFilter->GetOutput());
 viewer.Visualize();
 
 
 return EXIT_SUCCESS;

}

void CreateImage(ImageType::Pointer image) { 

 ImageType::IndexType start;
 start.Fill(0);

 ImageType::SizeType size;
 size.Fill(100);

 ImageType::RegionType region(start, size);
 image->SetRegions(region);
 image->Allocate();
 image->FillBuffer(0);

 // Make a square
 for(unsigned int r = 40; r < 60; r++)
   {
   for(unsigned int c = 40; c < 60; c++)
     {
     ImageType::IndexType pixelIndex;
     pixelIndex[0] = r;
     pixelIndex[1] = c;

     image->SetPixel(pixelIndex, 255);
     }
   }

} </source>

CMakeLists.txt

<source lang="cmake"> cmake_minimum_required(VERSION 2.6)

PROJECT(TranslationTransform)

include_directories(/home/doriad/ITKWikiExamples/ItkVtkGlue)

FIND_PACKAGE(VTK REQUIRED) INCLUDE(${VTK_USE_FILE})

FIND_PACKAGE(ITK REQUIRED) INCLUDE(${ITK_USE_FILE})

ADD_EXECUTABLE(TranslationTransform TranslationTransform.cxx /home/doriad/ITKWikiExamples/ItkVtkGlue/QuickView.cxx) TARGET_LINK_LIBRARIES(TranslationTransform vtkHybrid ITKBasicFilters ITKCommon ITKIO)

</source>

Retrieved from "https://public.kitware.com/Wiki/index.php?title=ITK/Examples/SimpleOperations/TranslationTransform&oldid=36840"