Difference between revisions of "ITK/Examples/Broken/Images/GradientImageFilter"

From KitwarePublic
< ITK‎ | Examples
Jump to: navigation, search
(Deprecated content that is moved to sphinx)
 
Line 1: Line 1:
==GradientImageFilter.cxx==
+
{{warning|1=The media wiki content on this page is no longer maintainedThe examples presented on the https://itk.org/Wiki/*  pages likely require ITK version 4.13 or earlier releasesIn many cases, the examples on this page no longer conform to the best practices for modern ITK versions.}}
<source lang="cpp">
 
#include "itkImage.h"
 
#include "itkCovariantVector.h"
 
#include "itkImageFileReader.h"
 
#include "itkImageFileWriter.h"
 
#include "itkGradientImageFilter.h"
 
 
 
#include <itkImageToVTKImageFilter.h>
 
 
 
#include "vtkPointData.h"
 
#include "vtkFloatArray.h"
 
#include "vtkImageData.h"
 
#include "vtkRenderWindowInteractor.h"
 
#include "vtkRenderWindow.h"
 
#include "vtkSmartPointer.h"
 
#include "vtkImageActor.h"
 
#include "vtkActor.h"
 
#include "vtkInteractorStyleImage.h"
 
#include "vtkRenderer.h"
 
#include "vtkGlyph3DMapper.h"
 
#include "vtkArrowSource.h"
 
 
 
static void VectorImageToVTKImage(itk::Image<itk::CovariantVector<float, 2>, 2>::Pointer vectorImage, vtkImageData* VTKImage);
 
 
 
int main(int argc, char * argv[])
 
{
 
  // Verify command line arguments
 
  if( argc < 2 )
 
    {
 
    std::cerr << "Usage: " << std::endl;
 
    std::cerr << argv[0] << "inputImageFile" << std::endl;
 
    return EXIT_FAILURE;
 
    }
 
 
 
  // Parse command line arguments
 
  std::string inputFilename = argv[1];
 
 
 
  // Setup types
 
  typedef itk::Image< float,  2 >  FloatImageType;
 
  typedef itk::Image< unsigned char, 2 >  UnsignedCharImageType;
 
 
 
  // Create and setup a reader
 
  typedef itk::ImageFileReader< UnsignedCharImageType >  ReaderType;
 
  ReaderType::Pointer reader = ReaderType::New();
 
  reader->SetFileName( inputFilename.c_str() );
 
 
 
  // Create and setup a gradient filter
 
  typedef itk::GradientImageFilter<
 
      UnsignedCharImageType, float> GradientFilterType;
 
  GradientFilterType::Pointer gradientFilter = GradientFilterType::New();
 
  gradientFilter->SetInput( reader->GetOutput() );
 
  gradientFilter->Update();
 
 
 
  // Visualize original image
 
  typedef itk::ImageToVTKImageFilter<UnsignedCharImageType> ConnectorType;
 
  ConnectorType::Pointer originalConnector = ConnectorType::New();
 
  originalConnector->SetInput(reader->GetOutput());
 
 
 
  vtkSmartPointer<vtkImageActor> originalActor =
 
    vtkSmartPointer<vtkImageActor>::New();
 
  originalActor->SetInput(originalConnector->GetOutput());
 
 
 
  // Visualize gradient
 
  vtkSmartPointer<vtkImageData> gradientImage =
 
    vtkSmartPointer<vtkImageData>::New();
 
  VectorImageToVTKImage(gradientFilter->GetOutput(), gradientImage);
 
 
 
  vtkSmartPointer<vtkArrowSource> arrowSource =
 
    vtkSmartPointer<vtkArrowSource>::New();
 
 
 
  vtkSmartPointer<vtkGlyph3DMapper> gradientMapper =
 
    vtkSmartPointer<vtkGlyph3DMapper>::New();
 
  gradientMapper->ScalingOn();
 
  gradientMapper->SetScaleFactor(.05);
 
  gradientMapper->SetSourceConnection(arrowSource->GetOutputPort());
 
  gradientMapper->SetInputConnection(gradientImage->GetProducerPort());
 
  gradientMapper->Update();
 
 
 
  vtkSmartPointer<vtkActor> gradientActor =
 
    vtkSmartPointer<vtkActor>::New();
 
  gradientActor->SetMapper(gradientMapper);
 
 
 
  // Visualize
 
  // 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<vtkRenderWindow> renderWindow =
 
    vtkSmartPointer<vtkRenderWindow>::New();
 
  renderWindow->SetSize(600,300);
 
 
 
  vtkSmartPointer<vtkRenderer> leftRenderer =
 
    vtkSmartPointer<vtkRenderer>::New();
 
  renderWindow->AddRenderer(leftRenderer);
 
  leftRenderer->SetViewport(leftViewport);
 
 
 
  vtkSmartPointer<vtkRenderer> rightRenderer =
 
    vtkSmartPointer<vtkRenderer>::New();
 
  renderWindow->AddRenderer(rightRenderer);
 
  rightRenderer->SetViewport(rightViewport);
 
  rightRenderer->SetBackground(1,0,0);
 
 
 
  leftRenderer->AddActor(originalActor);
 
  rightRenderer->AddActor(gradientActor);
 
 
 
  vtkSmartPointer<vtkRenderWindowInteractor> renderWindowInteractor =
 
    vtkSmartPointer<vtkRenderWindowInteractor>::New();
 
 
 
  vtkSmartPointer<vtkInteractorStyleImage> style =
 
    vtkSmartPointer<vtkInteractorStyleImage>::New();
 
  renderWindowInteractor->SetInteractorStyle(style);
 
 
 
  renderWindowInteractor->SetRenderWindow(renderWindow);
 
  renderWindowInteractor->Initialize();
 
 
 
  renderWindowInteractor->Start();
 
 
 
  return EXIT_SUCCESS;
 
}
 
 
 
void VectorImageToVTKImage(itk::Image<itk::CovariantVector<float, 2>, 2>::Pointer vectorImage, vtkImageData* VTKImage)
 
{
 
  itk::Image<itk::CovariantVector<float, 2>, 2>::RegionType region = vectorImage->GetLargestPossibleRegion();
 
  itk::Image<itk::CovariantVector<float, 2>, 2>::SizeType imageSize = region.GetSize();
 
  VTKImage->SetExtent(0, imageSize[0] -1, 0, imageSize[1] - 1, 0, 0);
 
 
 
  vtkSmartPointer<vtkFloatArray> vectors =
 
    vtkSmartPointer<vtkFloatArray>::New();
 
  vectors->SetNumberOfComponents(3);
 
  vectors->SetNumberOfTuples(imageSize[0] * imageSize[1]);
 
  vectors->SetName("GradientVectors");
 
 
 
  int counter = 0;
 
  for(unsigned int j = 0; j < imageSize[1]; j++)
 
    {
 
    for(unsigned int i = 0; i < imageSize[0]; i++)
 
      {
 
      itk::Image<itk::CovariantVector<float, 2>, 2>::IndexType index;
 
      index[0] = i;
 
      index[1] = j;
 
 
 
      itk::Image<itk::CovariantVector<float, 2>, 2>::PixelType pixel = vectorImage->GetPixel(index);
 
   
 
      float v[2];
 
      v[0] = pixel[0];
 
      v[1] = pixel[1];
 
      v[2] = 0;
 
      vectors->InsertTupleValue(counter, v);
 
      counter++;
 
      }
 
    }
 
  //std::cout << region << std::endl;
 
 
 
  VTKImage->GetPointData()->SetVectors(vectors);
 
}
 
</source>
 
 
 
{{ITKVTKCMakeLists|{{SUBPAGENAME}}}}
 

Latest revision as of 17:41, 7 June 2019

Warning: The media wiki content on this page is no longer maintained. The examples presented on the https://itk.org/Wiki/* pages likely require ITK version 4.13 or earlier releases. In many cases, the examples on this page no longer conform to the best practices for modern ITK versions.