- typedef itk::Image< PixelType, Dim > ConnectorImageType;
- typedef itk::ImageToVTKImageFilter <ConnectorImageType> ConnectorType;
- typename ConnectorType::Pointer connector = ConnectorType::New();
- connector->SetInput(input);
-
- try {
- connector->Update();
- }
- catch ( itk::ExceptionObject & err ) {
- std::cerr << "Error while setting vvImage from ITK (Dim==3)"
- << " " << err << std::endl;
- return vv_image;
- }
- vtkImageData *image = vtkImageData::New();
- image->DeepCopy(connector->GetOutput());
- vv_image->AddImage(image);
+ vv_image->AddItkImage<InputImageType>(input);
+ }
+ return vv_image;
+ }
+};
+
+template<class PixelType>
+struct vvImageFromITK_Impl<4u, PixelType>
+{
+ static vvImage::Pointer Do (typename itk::Image<PixelType,4u>::Pointer input, bool time_sequence=false)
+ {
+ vvImage::Pointer vv_image=vvImage::New();
+ typedef itk::Image< PixelType, 4 > InputImageType;
+
+ typedef itk::Image< PixelType, 3 > ItkImageType;
+ typedef itk::ExtractImageFilter<InputImageType, ItkImageType> FilterType;
+
+ //extract the 3D slices and put them in a std::vector<vtkImageData*>
+ input->UpdateOutputInformation();
+ typename InputImageType::RegionType inputRegion = input->GetLargestPossibleRegion();
+ typename InputImageType::SizeType inputSize = inputRegion.GetSize();
+ typename InputImageType::IndexType start = inputRegion.GetIndex();
+ typename InputImageType::SizeType extractedRegionSize = inputSize;
+ typename InputImageType::RegionType extractedRegion;
+ extractedRegionSize[3] = 0;
+ extractedRegion.SetSize(extractedRegionSize);
+
+ for (unsigned int i = 0; i < inputSize[3]; i++) {
+ start[3] = i;
+ extractedRegion.SetIndex(start);
+
+ typename FilterType::Pointer filter = FilterType::New();
+#if ITK_VERSION_MAJOR == 4
+ filter->SetDirectionCollapseToSubmatrix();
+#endif
+ filter->SetExtractionRegion(extractedRegion);
+ filter->SetInput(input);
+ filter->ReleaseDataFlagOn();
+ vv_image->AddItkImage<ItkImageType>(filter->GetOutput());