ITK  5.0.0
Insight Segmentation and Registration Toolkit
WikiExamples/ImageProcessing/ProcessingNthImageElement.cxx
using VectorImageType = itk::Image<itk::CovariantVector< float, 3>, 2>;
static void CreateImage(VectorImageType::Pointer image);
int main(int, char *[])
{
VectorImageType::Pointer image = VectorImageType::New();
CreateImage(image);
using ImageAdaptorType = itk::NthElementImageAdaptor<VectorImageType,
float>;
ImageAdaptorType::Pointer adaptor = ImageAdaptorType::New();
adaptor->SelectNthElement(0);
adaptor->SetImage(image);
BlurFilterType::Pointer blurFilter = BlurFilterType::New();
blurFilter->SetInput(adaptor);
blurFilter->Update();
return EXIT_SUCCESS;
}
void CreateImage(VectorImageType::Pointer image)
{
start.Fill(0);
size.Fill(2);
region.SetSize(size);
region.SetIndex(start);
image->SetRegions(region);
image->Allocate();
itk::ImageRegionIterator<VectorImageType> imageIterator(image,image->GetLargestPossibleRegion());
vec[0] = 1;
vec[1] = 2;
vec[2] = 3;
while(!imageIterator.IsAtEnd())
{
imageIterator.Set(vec);
++imageIterator;
}
}