This method attempts a direct integration of the Abel transform integral. It makes no assumptions about the data (apart from cylindrical symmetry), but it typically requires fine sampling to converge. Such methods are typically inefficient, but thanks to this Cython implementation (by Roman Yurchuk), this ‘direct’ method is competitive with the other methods.
How it works¶
Information about the algorithm and the numerical optimizations is contained in PR #52
When to use it¶
When a robust forward transform is required, this method works quite well. It is not typically recommended for the inverse transform, but it can work well for smooth functions that are finely sampled.
How to use it¶
To complete the forward or inverse transform of a full image with the direct method, simply use the
abel.Transform(myImage, method='direct', direction='forward').transform abel.Transform(myImage, method='direct', direction='inverse').transform
If you would like to access the Direct algorithm directly (to transform a right-side half-image), you can use