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 vote187 views1 answer
Float value in regression expression python
Python
20 votes204 views3 answers
How to convert all items in a list to integer python
Python
1 vote179 views1 answer
How to append leading zeros in python
Python
0 vote173 views1 answer
How to join models from another app
Python
41 votes187 views6 answers
How to fill nan values in pandas
Python
1 vote207 views1 answer
Keras import optimizer adam
Python
1 vote193 views1 answer
The operands of the logical operators should be boolean expressions, but python is not very strict. any nonzero number is interpreted as true.
Python
-2 votes215 views1 answer
Quotes generator python
Python
0 vote207 views1 answer
Back of list
Python
2 votes281 views1 answer
How to check how much time elapsed Python
Python