typename ResampleFilterType::Pointer resampler = ResampleFilterType::New();
resampler->SetInput(input);
resampler->SetOutputOrigin(input->GetOrigin());
typename ResampleFilterType::Pointer resampler = ResampleFilterType::New();
resampler->SetInput(input);
resampler->SetOutputOrigin(input->GetOrigin());