// Set input image iterator
typedef itk::ImageRegionIterator<ImageType> IteratorType;
IteratorType it(input1, input1->GetLargestPossibleRegion());
// typedef input2
// Set input image iterator
typedef itk::ImageRegionIterator<ImageType> IteratorType;
IteratorType it(input1, input1->GetLargestPossibleRegion());
// typedef input2
output->SetRegions(input1->GetLargestPossibleRegion());
output->SetOrigin(input1->GetOrigin());
output->SetSpacing(input1->GetSpacing());
output->SetRegions(input1->GetLargestPossibleRegion());
output->SetOrigin(input1->GetOrigin());
output->SetSpacing(input1->GetSpacing());