for (unsigned int pixelDim=0; pixelDim<Dimension; ++pixelDim)
{
size[pixelDim]=adaptor->GetLargestPossibleRegion().GetSize(pixelDim);
spacing[pixelDim]=input->GetSpacing()[pixelDim];
origin[pixelDim]=input->GetOrigin()[pixelDim];
for (unsigned int pixelDim=0; pixelDim<Dimension; ++pixelDim)
{
size[pixelDim]=adaptor->GetLargestPossibleRegion().GetSize(pixelDim);
spacing[pixelDim]=input->GetSpacing()[pixelDim];
origin[pixelDim]=input->GetOrigin()[pixelDim];
}
typename OutputImageType::RegionType region;
region.SetSize(size);
region.SetIndex(index);
output->SetRegions(region);
output->SetOrigin(origin);
}
typename OutputImageType::RegionType region;
region.SetSize(size);
region.SetIndex(index);
output->SetRegions(region);
output->SetOrigin(origin);