+ typedef itk::Image< PixelType, Dim > InputImageType;
+ typedef itk::Image< PixelType, Dim - 1 > 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[Dim - 1] = 0;
+ extractedRegion.SetSize(extractedRegionSize);
+
+ for (unsigned int i = 0; i < inputSize[Dim - 1]; i++) {
+ start[Dim - 1] = i;
+ extractedRegion.SetIndex(start);
+
+ typename FilterType::Pointer filter = FilterType::New();
+ filter->SetDirectionCollapseToSubmatrix();
+ filter->SetExtractionRegion(extractedRegion);
+ filter->SetInput(input);
+ filter->ReleaseDataFlagOn();
+ vv_image->AddItkImage<ItkImageType>(filter->GetOutput());
+ vv_image->ComputeScalarRangeBase<PixelType, Dim-1>(filter->GetOutput());
+ }
+ vv_image->SetTimeSpacing(input->GetSpacing()[Dim-1]);
+ vv_image->SetTimeOrigin(input->GetOrigin()[Dim-1]);
+}
+
+template<unsigned int Dim, class PixelType>
+struct vvImageFromITK_Impl
+{
+ static vvImage::Pointer Do (typename itk::Image<PixelType,Dim>::Pointer input, bool time_sequence=false)
+ {