1 /*=========================================================================
2 Program: vv http://www.creatis.insa-lyon.fr/rio/vv
5 - University of LYON http://www.universite-lyon.fr/
6 - Léon Bérard cancer center http://oncora1.lyon.fnclcc.fr
7 - CREATIS CNRS laboratory http://www.creatis.insa-lyon.fr
9 This software is distributed WITHOUT ANY WARRANTY; without even
10 the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
11 PURPOSE. See the copyright notices for more information.
13 It is distributed under dual licence
15 - BSD See included LICENSE.txt file
16 - CeCILL-B http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
17 ======================================================================-====*/
18 #ifndef VVIMAGEFROMITK_H
19 #define VVIMAGEFROMITK_H
21 //#include <itkImage.h>
22 #include <itkExtractImageFilter.h>
23 #include <itkImageToVTKImageFilter.h>
25 /**Converts the itk image to vv, handling the 4D problem
26 * The time_sequence boolean specifies that the image is to be interpreted as a time sequence,
27 * even if its dim is < 4. */
28 template<unsigned int Dim, class PixelType> vvImage::Pointer vvImageFromITK(typename itk::Image<PixelType,Dim>::Pointer input, bool time_sequence=false)
30 assert(Dim < 5 && Dim > 0); // We don't handle anything higher than 4-dimensional (for the moment :-p)
31 vvImage::Pointer vv_image=vvImage::New();
32 vv_image->Init(); //Delete any existing images
33 typedef itk::Image< PixelType, Dim > InputImageType;
35 if (Dim == 4 || time_sequence) //The time sequence case: create a series of VTK images
37 typedef itk::Image< PixelType, Dim - 1 > ConnectorImageType;
38 typedef itk::ImageToVTKImageFilter <ConnectorImageType> ConnectorType;
39 typedef itk::ExtractImageFilter<InputImageType,ConnectorImageType> FilterType;
41 typename FilterType::Pointer filter = FilterType::New();
42 typename ConnectorType::Pointer connector = ConnectorType::New();
44 //extract the 3D slices and put them in a std::vector<vtkImageData*>
45 typename InputImageType::RegionType inputRegion = input->GetLargestPossibleRegion();
46 typename InputImageType::SizeType inputSize = inputRegion.GetSize();
48 typename InputImageType::SizeType extractedRegionSize = inputSize;
49 typename InputImageType::RegionType extractedRegion;
50 extractedRegionSize[Dim - 1] = 0;
51 extractedRegion.SetSize(extractedRegionSize);
53 filter->SetInput(input);
54 connector->SetInput(filter->GetOutput());
56 typename InputImageType::IndexType start = inputRegion.GetIndex();
58 for (unsigned int i = 0; i < inputSize[Dim - 1]; i++) {
60 extractedRegion.SetIndex(start);
61 filter->SetExtractionRegion(extractedRegion);
65 catch ( itk::ExceptionObject & err ) {
66 std::cerr << "Error while setting vvImage from ITK (Dim==4) [Extract phase]"
67 << " " << err << std::endl;
73 catch ( itk::ExceptionObject & err ) {
74 std::cerr << "Error while setting vvImage from ITK (Dim==4) [Connect phase]"
75 << " " << err << std::endl;
78 vtkImageData *image = vtkImageData::New();
79 image->DeepCopy(connector->GetOutput());
80 vv_image->AddImage(image);
82 vv_image->SetTimeSpacing(input->GetSpacing()[Dim-1]);
83 vv_image->SetTimeOrigin(input->GetOrigin()[Dim-1]);
85 else //Dim == 1,2,3 and not time_sequence
87 typedef itk::Image< PixelType, Dim > ConnectorImageType;
88 typedef itk::ImageToVTKImageFilter <ConnectorImageType> ConnectorType;
89 typename ConnectorType::Pointer connector = ConnectorType::New();
90 connector->SetInput(input);
95 catch ( itk::ExceptionObject & err ) {
96 std::cerr << "Error while setting vvImage from ITK (Dim==3)"
97 << " " << err << std::endl;
100 vtkImageData *image = vtkImageData::New();
101 image->DeepCopy(connector->GetOutput());
102 vv_image->AddImage(image);
107 #endif //vvImageFromITK