IteratorType it(input1, input1->GetLargestPossibleRegion());
// typedef input2
- typename ImageType::Pointer input2 = NULL;
+ typename ImageType::Pointer input2 = ITK_NULLPTR;
IteratorType it2;
// Special case for normalisation
output->SetRegions(input1->GetLargestPossibleRegion());
output->SetOrigin(input1->GetOrigin());
output->SetSpacing(input1->GetSpacing());
+ output->SetDirection(input1->GetDirection());
output->Allocate();
// Set output iterator
typedef itk::ImageRegionIterator<OutputImageType> IteratorOutputType;