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

40 votes209 views3 answers
Python push to dataframe pandas
Python
-2 votes217 views2 answers
Create virtual environment python stack overflow
Python
1 vote220 views1 answer
Require http method django view
Python
23 votes204 views2 answers
How to replace zero value in python dataframe
Python
1 vote237 views1 answer
Error popup in django not visible
Python Django
63 votes254 views4 answers
Python find in string
Python
-1 vote216 views2 answers
Pandas dataframe apply function with multiple arguments
Python
7 votes235 views5 answers
Convert categorical variable to numeric python
Python
0 vote211 views1 answer
Python tableau a 2 dimensions
Python
1 vote218 views1 answer
Bezier curve python
Python