Very Fast Radial Profile

Radial profile of an image without a "for loop". At size 2048x1792 it works 50 times faster than prevalent "for loop" variants
Updated 5 Nov 2021

View License

Over the years there were several questions posted here at MathWorks about how to compute a radial profile (like average spectral density). The same answer was repeated every time: go in a for loop over all radii and average the data which lays in a circle. When dealing with an image of just 1024x1024 it starts to be faster to deconvolve an image than to compute a radial profile!
MatLab has a very fast own function accumarray which can do the trick with zero for loops. I have written a basic radial profiler which calls accumarray twice: first time to create radial ticks, and second time to create the average profile. In my application with 2048x1792 data I have ahieved 50x accelearation compared with the "for loop" approach.
If you need to speed-up this function even more, consider using bwdist instead of meshgrid.
This function can be further improved further in terms of accuracy of averaging by taking into account subpixel subtelties.

Cite As

Sergey Loginov (2024). Very Fast Radial Profile (, MATLAB Central File Exchange. Retrieved .

MATLAB Release Compatibility
Created with R2021b
Compatible with any release
Platform Compatibility
Windows macOS Linux

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!
Version Published Release Notes

typos in the summary