[Insight-users] Jacobian determinant image of the deformation field

[Mathieu De Craene]

Hi ITK folks,


I found this email on the mailing list

        http://public.kitware.com/pipermail/insight-users/2006-August/019246.html

and I think these two users are right.

I was wondering if we shouldn't add the following line at line 250 of
itkDeformationFieldJacobianDeterminantFilter.h to correct the problem

        if (i == j) J[i][j] += 1.;


so that we effectively compute


        det[ dT/dx ] = det[ I + du/dx ]
        
        
        
Thanks for any comment,



Mathieu