Hi, I found this to be faster, hope it helps, S.
%% code for fast radial averaging:
[~,~,iir]=unique(ri); %find unique values
Seems really useful set of functions. The compiling of the mex files can be a bit tricky though.
I would like to ask: why is only single precision supported? This really creates problem for me because I have an application that sensitive on the accuracy.