ITK  5.0.0
Insight Segmentation and Registration Toolkit
WikiExamples/VectorImages/JoinImageFilter.cxx
#include "itkImage.h"
using ImageType = itk::Image<unsigned char, 2>;
static void CreateImage(ImageType::Pointer image, unsigned char value);
int main(int, char *[])
{
ImageType::Pointer image1 = ImageType::New();
CreateImage(image1, 0);
ImageType::Pointer image2 = ImageType::New();
CreateImage(image2, 10);
using JoinImageFilterType = itk::JoinImageFilter<ImageType, ImageType>;
JoinImageFilterType::Pointer joinFilter = JoinImageFilterType::New();
joinFilter->SetInput1(image1);
joinFilter->SetInput2(image2);
joinFilter->Update();
index[0] = 0;
index[1] = 0;
std::cout << static_cast<int>(joinFilter->GetOutput()->GetPixel(index)[0]) << std::endl;
std::cout << static_cast<int>(joinFilter->GetOutput()->GetPixel(index)[1]) << std::endl;
return EXIT_SUCCESS;
}
void CreateImage(ImageType::Pointer image, unsigned char value)
{
// Create an image
start.Fill(0);
size.Fill(100);
ImageType::RegionType region(start,size);
image->SetRegions(region);
image->Allocate();
image->FillBuffer(value);
}