Hi Lester, Without further information, the simplest approach would be to write two nested loops. Now, regarding the computation of radially averaged PSDs, you may want to take a look in the web for published code. For example, look at the function raPsd2d, written by Ruzanski (2009) for Matlab. Regards, Rafael