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

0 vote349 views1 answer
Cm to m in python
Python
0 vote369 views1 answer
App is not a registered namespace django
Python Django
0 vote332 views1 answer
Add comma to csv in pandas
Python
0 vote385 views1 answer
Read dict from text
Python
1 vote759 views1 answer
How to check which submit button is clicked in flask wtf
Python
0 vote361 views1 answer
Pivot tables pandas explicación
Python
17 votes313 views4 answers
Flask render_template
Python
1 vote348 views2 answers
Multiple values in python loop for x,y
Python
1 vote337 views1 answer
What happens when you use the built-in function any() on a list?
Python
0 vote385 views1 answer
Plural name django
Python Django