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 vote129 views1 answer
Pandas write to csv without first line
Python
72 votes133 views5 answers
How to run for loop in python
Python
5 votes130 views1 answer
Mongodb group by having
Python
1 vote122 views1 answer
How to make a use of list in python to make your own length function
Python
0 vote121 views1 answer
Flask gunicorn get ip
Python
2 votes124 views2 answers
Python how to sum two lists
Python
0 vote134 views1 answer
Ros python service server
Python
1 vote187 views2 answers
Divide a value by all values in a list
Python
0 vote129 views1 answer
Matplotlib overlapping labels
Python
0 vote127 views1 answer
Nested ternary operator python
Python