Current code:
return this->GetIndexToWorldTransform()->TransformPoint(m_BoundingBox->GetCenter());
Proposed:
if (m_ImageGeometry)
{
Point3D c = m_BoundingBox->GetCenter();
c[0] -= 0.5;
c[1] -= 0.5;
c[2] -= 0.5;
IndexToWorld(c, c);
return c;
}