ITK/Examples/ImageProcessing/RegionOfInterestImageFilter: Difference between revisions

From KitwarePublic
< ITK‎ | Examples
Jump to navigationJump to search
(Use QuickView)
Line 4: Line 4:
#include "itkImageFileReader.h"
#include "itkImageFileReader.h"
#include "itkRegionOfInterestImageFilter.h"
#include "itkRegionOfInterestImageFilter.h"
#include "itkRGBPixel.h"


#include "itkImageToVTKImageFilter.h"
#include "QuickView.h"
 
#include "vtkImageViewer.h"
#include "vtkRenderWindowInteractor.h"
#include "vtkSmartPointer.h"
#include "vtkImageActor.h"
#include "vtkInteractorStyleImage.h"
#include "vtkRenderer.h"


#include <iostream>
#include <iostream>
Line 19: Line 13:
int main(int argc, char *argv[])
int main(int argc, char *argv[])
{
{
if(argc < 2)
  if(argc < 2)
{
    {
std::cerr << "Required: filename" << std::endl;
    std::cerr << "Usage: ";
return EXIT_FAILURE;
    std::cerr << argv[0] << " inputImageFile [start] [size]" << std::endl;
}
    return EXIT_FAILURE;
 
    }
std::string filename = argv[1];
 
typedef itk::Image<unsigned char, 2> ImageType;
typedef itk::ImageFileReader<ImageType> ReaderType;
typedef itk::ImageToVTKImageFilter<ImageType> ConnectorType;
 
ReaderType::Pointer reader = ReaderType::New();
ConnectorType::Pointer originalConnector = ConnectorType::New();
ConnectorType::Pointer extractedConnector = ConnectorType::New();
 
reader->SetFileName(filename.c_str());
reader->Update();
 
ImageType::SizeType inSize = reader->GetOutput()->GetLargestPossibleRegion().GetSize();
 
originalConnector->SetInput(reader->GetOutput());
 
vtkSmartPointer<vtkImageActor> originalActor =
vtkSmartPointer<vtkImageActor>::New();
originalActor->SetInput(originalConnector->GetOutput());
 
typedef itk::RegionOfInterestImageFilter< ImageType, ImageType > FilterType;
FilterType::Pointer filter = FilterType::New();
 
ImageType::IndexType start;
start[0] = inSize[0]/4;
start[1] = inSize[1]/4;
 
ImageType::SizeType size;
size[0] = inSize[0]/4;
size[1] = inSize[1]/4;
 
ImageType::RegionType desiredRegion;
desiredRegion.SetSize(size);
desiredRegion.SetIndex(start);


filter->SetRegionOfInterest(desiredRegion);
  std::string filename = argv[1];
filter->SetInput(reader->GetOutput());


extractedConnector->SetInput(filter->GetOutput());
  typedef itk::Image<itk::RGBPixel<unsigned char>, 2> ImageType;
  typedef itk::ImageFileReader<ImageType>            ReaderType;


vtkSmartPointer<vtkImageActor> extractedActor =
  ReaderType::Pointer reader = ReaderType::New();
vtkSmartPointer<vtkImageActor>::New();
  reader->SetFileName(filename.c_str());
extractedActor->SetInput(extractedConnector->GetOutput());
  reader->Update();


// There will be one render window
  ImageType::SizeType inSize = reader->GetOutput()->GetLargestPossibleRegion().GetSize();
vtkSmartPointer<vtkRenderWindow> renderWindow =
vtkSmartPointer<vtkRenderWindow>::New();
renderWindow->SetSize(600, 300);


vtkSmartPointer<vtkRenderWindowInteractor> interactor =
  typedef itk::RegionOfInterestImageFilter< ImageType, ImageType > FilterType;
vtkSmartPointer<vtkRenderWindowInteractor>::New();
  FilterType::Pointer filter = FilterType::New();
interactor->SetRenderWindow(renderWindow);


// Define viewport ranges
  ImageType::IndexType start;
// (xmin, ymin, xmax, ymax)
  start[0] = inSize[0]/4;
double leftViewport[4] = {0.0, 0.0, 0.5, 1.0};
  start[1] = inSize[1]/4;
double rightViewport[4] = {0.5, 0.0, 1.0, 1.0};


// Setup both renderers
  ImageType::SizeType size;
vtkSmartPointer<vtkRenderer> leftRenderer =
  size[0] = inSize[0]/4;
vtkSmartPointer<vtkRenderer>::New();
  size[1] = inSize[1]/4;
renderWindow->AddRenderer(leftRenderer);
leftRenderer->SetViewport(leftViewport);
leftRenderer->SetBackground(.6, .5, .4);


vtkSmartPointer<vtkRenderer> rightRenderer =
  if (argc > 2)
vtkSmartPointer<vtkRenderer>::New();
    {
renderWindow->AddRenderer(rightRenderer);
    start.Fill(atoi(argv[2]));
rightRenderer->SetViewport(rightViewport);
    }
rightRenderer->SetBackground(.4, .5, .6);
  if (argc > 3)
    {
    size.Fill(atoi(argv[3]));
    }


// Add the sphere to the left and the cube to the right
leftRenderer->AddActor(originalActor);
rightRenderer->AddActor(extractedActor);


leftRenderer->ResetCamera();
  ImageType::RegionType desiredRegion;
rightRenderer->ResetCamera();
  desiredRegion.SetSize(size);
  desiredRegion.SetIndex(start);


renderWindow->Render();
  filter->SetRegionOfInterest(desiredRegion);
  filter->SetInput(reader->GetOutput());


vtkSmartPointer<vtkInteractorStyleImage> style =
  QuickView viewer;
vtkSmartPointer<vtkInteractorStyleImage>::New();
  viewer.AddRGBImage(
    reader->GetOutput(),true,
    itksys::SystemTools::GetFilenameName(argv[1]));


interactor->SetInteractorStyle(style);
  std::stringstream desc;
  desc << "RegionOfInterestImageFilter\nstart = " << start
      << " size = " << size;
  viewer.AddRGBImage(
    filter->GetOutput(),
    true,
    desc.str());


interactor->Start();
  viewer.Visualize();


return EXIT_SUCCESS;
  return EXIT_SUCCESS;
}
}
</source>
</source>

Revision as of 05:25, 23 December 2010

RegionOfInterestImageFilter.cxx

<source lang="cpp">

  1. include "itkImage.h"
  2. include "itkImageFileReader.h"
  3. include "itkRegionOfInterestImageFilter.h"
  4. include "itkRGBPixel.h"
  1. include "QuickView.h"
  1. include <iostream>
  2. include <string>

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

 if(argc < 2)
   {
   std::cerr << "Usage: ";
   std::cerr << argv[0] << " inputImageFile [start] [size]" << std::endl;
   return EXIT_FAILURE;
   }

 std::string filename = argv[1];

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

 ReaderType::Pointer reader = ReaderType::New();
 reader->SetFileName(filename.c_str());
 reader->Update();

 ImageType::SizeType inSize = reader->GetOutput()->GetLargestPossibleRegion().GetSize();

 typedef itk::RegionOfInterestImageFilter< ImageType, ImageType > FilterType;
 FilterType::Pointer filter = FilterType::New();

 ImageType::IndexType start;
 start[0] = inSize[0]/4;
 start[1] = inSize[1]/4;

 ImageType::SizeType size;
 size[0] = inSize[0]/4;
 size[1] = inSize[1]/4;

 if (argc > 2)
   {
   start.Fill(atoi(argv[2]));
   }
 if (argc > 3)
   {
   size.Fill(atoi(argv[3]));
   }



 ImageType::RegionType desiredRegion;
 desiredRegion.SetSize(size);
 desiredRegion.SetIndex(start);

 filter->SetRegionOfInterest(desiredRegion);
 filter->SetInput(reader->GetOutput());

 QuickView viewer;
 viewer.AddRGBImage(
   reader->GetOutput(),true,
   itksys::SystemTools::GetFilenameName(argv[1]));  

 std::stringstream desc;
 desc << "RegionOfInterestImageFilter\nstart = " << start
      << " size = " << size;
 viewer.AddRGBImage(
   filter->GetOutput(),
   true,
   desc.str());  

 viewer.Visualize();

 return EXIT_SUCCESS;

} </source>

CMakeLists.txt

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

PROJECT(RegionOfInterestImageFilter)

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

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

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

ADD_EXECUTABLE(RegionOfInterestImageFilter RegionOfInterestImageFilter.cpp) TARGET_LINK_LIBRARIES(RegionOfInterestImageFilter vtkHybrid ITKBasicFilters ITKCommon ITKIO)

</source>

Retrieved from "https://public.kitware.com/Wiki/index.php?title=ITK/Examples/ImageProcessing/RegionOfInterestImageFilter&oldid=35923"