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 vote217 views1 answer
Extract value from tensor pytorch
Python
4 votes222 views2 answers
The list of prime number in a given range python
Python
9 votes212 views3 answers
Create qr code python
Python
1 vote222 views1 answer
Display max number of columns pandas
Python
53 votes226 views3 answers
How to find length of list python
Python
1 vote224 views1 answer
Stop procedure python
Python
2 votes203 views1 answer
Vowel and consonant list python
Python
0 vote253 views1 answer
Activate inherit function django
Python
5 votes265 views1 answer
Python reverse linked list
Python
0 vote255 views1 answer
Count proportion pandas
Python