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 vote177 views2 answers
Pandas sub dataframe
Python
0 vote182 views1 answer
Python data structures 9.4
Python
0 vote193 views1 answer
Print elements without print function in python
Python
2 votes192 views1 answer
Private variables python
Python
2 votes168 views1 answer
Python only numbers in string
Python
0 vote185 views1 answer
Python googledriver download
Python
3 votes204 views2 answers
ModuleNotFoundError: No module named 'skvideo'
Python
27 votes185 views3 answers
Python date and time to datetime
Python
0 vote196 views2 answers
Instabot source code python library
Python
80 votes183 views4 answers
Save pandas into csv
Python