}
// Compute origin based on image corner
- typename FilterType::OriginPointType origin = input->GetOrigin();
for(unsigned int i=0; i<OutputImageType::ImageDimension; i++) {
- origin[i] -= 0.5 * input->GetSpacing()[i];
- origin[i] += 0.5 * m_OutputSpacing[i];
+ m_OutputOrigin[i] -= 0.5 * input->GetSpacing()[i];
+ m_OutputOrigin[i] += 0.5 * m_OutputSpacing[i];
}
// Instance of the transform object to be passed to the resample
filter->SetTransform(m_Transform);
filter->SetSize(m_OutputSize);
filter->SetOutputSpacing(m_OutputSpacing);
- filter->SetOutputOrigin(origin);
+ filter->SetOutputOrigin(m_OutputOrigin);
filter->SetDefaultPixelValue(m_DefaultPixelValue);
filter->SetNumberOfThreads(this->GetNumberOfThreads());
filter->SetOutputDirection(m_OutputDirection); // <-- NEEDED if we want to keep orientation (in case of PermutAxes for example)