Integrate ChambollePock with tomography package
Opened this issue · 1 comments
I maintain a tomography processing package (https://gitlab.com/aaherzing/tomotools) and I'm interested in integrating your C-P algorithm to my reconstruction module to add a TV minimization option. Does this sound reasonable?
I see that there are plans to integrate your algorithm with Astra, so this would simply be a stop-gap measure until that functionality is available.
Hello @AndrewHerzing
Wow, I did not even remember to have created this repository. Please go ahead if you find any of this code useful.
Regarding the Chambolle-Pock algorithm integration in Astra, I think it will not materialize. We had a discussion with Astra core developer, and the spirit is to bring "simple" building blocks for tomography. Total Variation stats to be a bit advanced.
In order to fully benefit from the performances of Astra when doing iterative reconstruction, the output of Projection/Backprojection have to "stay" on GPU. There is an option in parallel (back)projector to pass a Cuda memory pointer, but unfortunately this is currently available only for the 3D (back)projector. This feature should be added in 2D in the future, though.