radialAverage¶
-
image_tools.radialAverage(image, center=None, stddev=False, returnAz=False, return_naz=False, binsize=1.0, weights=None, steps=False, interpnan=False, left=None, right=None, mask=None, symmetric=None)[source] [edit on github]¶ Calculate the radially averaged azimuthal profile. (this code has not been optimized; it could be speed boosted by ~20x)
image - The 2D image center - The [x,y] pixel coordinates used as the center. The default is
None, which then uses the center of the image (including fractional pixels).stddev - if specified, return the radial standard deviation instead of the average returnAz - if specified, return (azimuthArray,azimuthal_profile) return_naz - if specified, return number of pixels per azimuth and azimuth binsize - size of the averaging bin. Can lead to strange results if
non-binsize factors are used to specify the center and the binsize is too large- weights - can do a weighted average instead of a simple average if this keyword parameter
- is set. weights.shape must = image.shape. weighted stddev is undefined, so don’t set weights and stddev.
- steps - if specified, will return a double-length bin array and azimuthal
- profile so you can plot a step-form azimuthal profile (which more accurately represents what’s going on)
- interpnan - Interpolate over NAN values, i.e. bins where there is no data?
- left,right - passed to interpnan; they set the extrapolated values
- mask - can supply a mask (boolean array same size as image with True for OK and False for not)
- to average over only select data.
If a bin contains NO DATA, it will have a NAN value because of the divide-by-sum-of-weights component. I think this is a useful way to denote lack of data, but users let me know if an alternative is prefered...