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 vote156 views1 answer
Ordered dict
Python
2 votes136 views1 answer
Django models distinct
Python
0 vote139 views1 answer
\r\n python
Python
1 vote135 views1 answer
Module "django.contrib.auth.hashers" does not define a "BcryptPasswordHasher" attribute/class
Python
0 vote140 views1 answer
Python get github file content
Python
0 vote143 views1 answer
Python get ids from array of objects
Python
0 vote141 views1 answer
Function continuity python
Python
0 vote191 views1 answer
Plot bar chart python with resulting numbers
Python
1 vote177 views1 answer
Django proper capitalization case jinja
Python Django
1 vote155 views2 answers
How to find 1 st digit in python
Python