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 vote307 views1 answer
Pandas series example
Python
0 vote331 views1 answer
Python - How To Check if a String Is a Palindrome
Python
0 vote361 views1 answer
Django create object from dict
Python Django
3 votes329 views2 answers
Dict_values to list python
Python
4 votes300 views4 answers
Python save output to file
Python
34 votes343 views6 answers
Round python
Python
2 votes300 views1 answer
How to clear an array python
Python
0 vote308 views1 answer
16 bit floating point numpy
Python
-1 vote325 views2 answers
Get tuple value python
Python
0 vote379 views1 answer
Using claudinary in django
Python Django