m-pilia / disptools

Generate displacement fields with known volume changes
https://martinopilia.com/disptools
MIT License
23 stars 3 forks source link

Possibly implement "geometric" based jacobian calculation #6

Open gdevenyi opened 5 years ago

gdevenyi commented 5 years ago

ANTs has a "geometric" jacobian implementation, which computes the jacobian based on the change in volume of a tetrahedra centered at each voxel.

https://github.com/ANTsX/ANTs/blob/master/Utilities/itkGeometricJacobianDeterminantImageFilter.hxx

https://github.com/ANTsX/ANTs/blob/master/Examples/CreateJacobianDeterminantImage.cxx

gdevenyi commented 5 years ago

Related paper: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2823935/

m-pilia commented 5 years ago

Thanks for the suggestion! I had not seen that paper before, and it seems an interesting approach. In this context it would take some additional effort to implement, however, because the tricky part is to get a good cost function for gradient descent. It is an interesting possibility, however I have little time budget for disptools nowadays, so I am not sure it is something I can tackle in the short term.