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 vote228 views2 answers
Flask return html
Python
0 vote239 views1 answer
Pyhton apend to list
Python
3 votes249 views1 answer
N unique random numbers in python
Python
26 votes245 views1 answer
Dict from two lists
Python
0 vote240 views1 answer
Linkedin dynamic scrolling using selenium python
Python
10 votes237 views3 answers
Python infinite loop
Python
4 votes192 views2 answers
How to remove all elements from dictionary python
Python
2 votes258 views1 answer
Downgrade python 3.8 to 3.7 ubuntu
Python
1 vote253 views1 answer
Python selenium itemprop
Python
-1 vote250 views1 answer
Python convert integer to signed base 2 complement
Python