I'm trying to implement an optimization algorithm in Python for solving a computerized tomography problem with TV regularization. I know I can use the function "radon" from scikit-image, but the point it that I also need the transpose (or adjoint operator) of the Radon transform as well.
Does anyone know if there is a way to compute the transpose matrix of the Radon Transform in scikit-image?
By reading the documentation, I don't think that Scikit Radon transform has this feature. I would therefore not use it on my own, since the adjoint is the most needed operation in scientific computing. Browsing through the web, I saw at least one library (pylops) offering this feature: https://pylops.readthedocs.io/en/latest/gallery/plot_radon.html
The other option is to reprogram it as a sparse matrix and use standard matrix-vector products...