radialavg.zip

Radially averaqe 2D square matrix z into m bins and demonstration routine

[Zr, R] = RADIALAVG(z,m,xo,yo) computes the average along the radius of a unit circle inscribed in the square matrix z. The average is computed in M bins. The radial average is computed 1/2 bin beyond the unit circle, but all values beyond that (i.e. corners of the matrix z) are excluded from the calculation. The radial average is returned in Zr and the mid-points of the M bins are returned in vector R. RADIALAVG will fail if Not a Number (NaN) values are in the data z; they are not excluded from the calculation. If offset values xo,yo are used, the origin (0,0) of the unit circle about which the RADIALAVG is computed is offset by xo and yo relative to the origin of the unit square of the input z matrix.