Page MenuHomePhabricator

BackTransform in Geometry3D still uses the unstable ITK Matirx inversion
Closed, ResolvedPublic

Description

BackTransform for a mitk::Vector is not affected by the world coordinate voxel center shift.
Still, it uses the ITK method Transform->GetInverseMatrix(), which has shown to be unstable on matrix inversion for 3D matrices.

It should be replaced as shown in T1995. (By Transform->GetInverse->GetMatrix) which seems to be stable.

Daniel M., again, please comment on this.

Event Timeline

I'd try to look at this together tomorrow afternoon.

[SVN revision 20058]
FIX (#2884): Use stable matrix inversion; relocated multiplication to separate function.

Fixed by using ->GetInverse.
Relocated Multiplication with inverse matrix to seperate method
MultiplyWithWorldToIndex.