| 90 | |
| 91 | template <typename TPixel, unsigned int VDimensions> |
| 92 | static void SplitVectorImage(const itk::VectorImage<TPixel, VDimensions>* image, |
| 93 | std::vector<mitk::Image::Pointer>& result) |
| 94 | { |
| 95 | typedef itk::VectorImage<TPixel, VDimensions> VectorImageType; |
| 96 | typedef itk::Image<TPixel, VDimensions> ImageType; |
| 97 | typedef itk::VectorIndexSelectionCastImageFilter<VectorImageType, ImageType> VectorIndexSelectorType; |
| 98 | |
| 99 | auto numberOfLayers = image->GetVectorLength(); |
| 100 | for (decltype(numberOfLayers) layer = 0; layer < numberOfLayers; ++layer) |
| 101 | { |
| 102 | auto layerSelector = VectorIndexSelectorType::New(); |
| 103 | layerSelector->SetInput(image); |
| 104 | layerSelector->SetIndex(layer); |
| 105 | layerSelector->Update(); |
| 106 | |
| 107 | mitk::Image::Pointer layerImage = mitk::GrabItkImageMemoryChannel(layerSelector->GetOutput(), nullptr, nullptr, false); |
| 108 | result.push_back(layerImage); |
| 109 | } |
| 110 | } |
| 111 | |
| 112 | std::vector<mitk::Image::Pointer> mitk::SplitVectorImage(const Image* vecImage) |
| 113 | { |