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

10 votes207 views1 answer
Python 3 text file leng
Python
8 votes182 views2 answers
How to move a column in pandas dataframe
Python
0 vote182 views2 answers
Scikit learn introduction
Python
1 vote188 views1 answer
Title tikinter
Python
0 vote228 views1 answer
How to close windows in selenium python without quitting the browser
Python
0 vote234 views1 answer
What is a good django orm cookbook
Python
0 vote229 views1 answer
How to set text in qdateedit pyqt5
Python
1 vote226 views1 answer
Join pyspark stackoverflow
Python
0 vote199 views1 answer
'djdt' is not a registered namespace
Python
28 votes190 views3 answers
UnicodeDecodeError: 'utf-8' codec can't decode byte 0xa5 in position 10: invalid start byte
Python