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

57 votes233 views6 answers
Print current date and time in python
Python
9 votes220 views6 answers
Python sort list
Python
12 votes261 views4 answers
Get index of all element in list python
Python
3 votes229 views2 answers
Multiply all values in column pandas
Python
1 vote244 views1 answer
Django ModelChoiceField value not id
Python Django
0 vote247 views1 answer
Ensure string length 2 python
Python
0 vote229 views1 answer
Requests.packages.urllib3.util.retry could not be resolved from source
Python
20 votes244 views2 answers
How to create string in python
Python
5 votes233 views3 answers
How to make addition in python
Python
1 vote226 views1 answer
Python write text file on the next line
Python