Matplotlib radial averaging

import numpy as np
import matplotlib.pyplot as plt
from skimage import data

# get some image
image = data.coins()
image = image[:,0:303]

# create array of radii
x,y = np.meshgrid(np.arange(image.shape[1]),np.arange(image.shape[0]))
R = np.sqrt(x**2+y**2)

# calculate the mean
f = lambda r : image[(R >= r-.5) & (R < r+.5)].mean()
r  = np.linspace(1,302,num=302)
mean = np.vectorize(f)(r)

# plot it
fig,ax=plt.subplots()
ax.plot(r,mean)
plt.show()

Source: stackoverflow.com

Similar codes

1 vote270 views1 answer
Flask app route in different files
Python
0 vote334 views1 answer
Change order of barh matplotlib
Python
3 votes298 views1 answer
Visual studio code import library python
Python
0 vote328 views1 answer
How to kill a script if error is hit python
Python
0 vote335 views1 answer
Stop a subprocess python
Python
0 vote300 views1 answer
Installing intel-numpy
Python
1 vote346 views5 answers
Embed python discord
Python
1 vote311 views1 answer
How to change the values of a column in numpy array
Python
0 vote325 views1 answer
Cos inverse in python numpy
Python
0 vote328 views1 answer
R - scale data from 0 to 1
Python