Dear FSL experts,
Is there any easy way of calculating the full Jacobian matrices from the deformation fields?
I'm trying to add a function to FnirtFileReader similar to ::Jacobian, whose loop does this:
jac(i,j,k,0) = (a11+(1.0/dx.Vxs_x())*dxdx->element(indx));
jac(i,j,k,1) = (a12+(1.0/dx.Vxs_y())*dxdy->element(indx));
jac(i,j,k,2) = (a13+(1.0/dx.Vxs_z())*dxdz->element(indx));
jac(i,j,k,3) = (a21+(1.0/dy.Vxs_x())*dydx->element(indx));
jac(i,j,k,4) = (a22+(1.0/dy.Vxs_y())*dydy->element(indx));
jac(i,j,k,5) = (a23+(1.0/dy.Vxs_z())*dydz->element(indx));
jac(i,j,k,6) = (a31+(1.0/dz.Vxs_x())*dzdx->element(indx));
jac(i,j,k,7) = (a32+(1.0/dz.Vxs_y())*dzdy->element(indx));
jac(i,j,k,8) = (a33+(1.0/dz.Vxs_z())*dzdz->element(indx));
(jac is a volume4D)
Would this be correct? Or am I killing a fly with a sledgehammer?
Thanks!
My kind regards,
Alexandre Manhães Savio
Grupo de Inteligencia Computacional (GIC)
Universidad del País Vasco
|