I’m trying to create a c program that implements the radon transform algorithm. I know that for an image f(x,y) a collection of g(phi,s), at all phi is the randon transform of the image, where g(phi,s) is defined as:
Now, I don’t have a lot of experience with C and I always used an external library with Java and C# to perform complex math operations. I’m having a hard time finding one for c, I’m also having trouble creating a function to do it, I’ve been looking into numerical integration but that is for definite integrals. Any help would be appreciated
Disclaimer: there may be specialized implementations of Radon transform known to the image processing community. I don’t know about those, I just approach the problem as an integration problem. I see eg. http://takinginitiative.net/2008/04/02/radon-transform-c-implementation-update/ that there may be some specific method for computing Radon transforms of images.
What you are really doing here is a 1D integral along a line. See the 4th formula in the Wikipedia article http://en.wikipedia.org/wiki/Radon_transform .
If you are doing this over an image, the method you use for computing the integral will be quite dependent on the interpolation you choose. If you go for bilinear interpolation, then a plain (adaptive) trapezoidal rule will likely give you good results.
For C libraries to assist you, you can have a look at the GNU Scientific Library.