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

12 votes260 views4 answers
'chromedriver' executable needs to be in PATH
Python
0 vote270 views1 answer
Python random email generator
Python
0 vote254 views1 answer
Python + credit-german.csv + class
Python
3 votes290 views1 answer
Python split lines
Python
0 vote271 views1 answer
Python break list into n lists
Python
3 votes254 views3 answers
Python check if variable has value
Python
0 vote266 views1 answer
Python select columns names from dataframe
Python
0 vote270 views1 answer
Figure in matplotlib
Python
0 vote274 views2 answers
How to print the freq of each char by using dict in python
Python
0 vote275 views1 answer
Read ms word with python
Python