+
+ // force index to zero
+ typename ImageType::RegionType region;
+ region = output->GetLargestPossibleRegion();
+ typename ImageType::IndexType index;
+ index.Fill(itk::NumericTraits<double>::Zero);
+ region.SetIndex(index);
+ output->SetRegions(region);