- if (Dim == 4 || time_sequence) //The time sequence case: create a series of VTK images
- {
- 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();
+ for (unsigned int i = 0; i < inputSize[Dim - 1]; i++) {
+ start[Dim - 1] = i;
+ extractedRegion.SetIndex(start);
+
+ typename FilterType::Pointer filter = FilterType::New();