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

3 votes282 views1 answer
What is cpython
Python
0 vote282 views1 answer
Numpy generlized ufunc
Python
0 vote282 views1 answer
Get return value of contract call in brownie
Python
11 votes283 views3 answers
Pandas count values by column
Python
15 votes305 views5 answers
Python http server
Python
2 votes282 views2 answers
Np float to int
Python
0 vote286 views1 answer
How to copy list python
Python
19 votes278 views4 answers
Convert column series to datetime in pandas dataframe
Python
0 vote263 views1 answer
Use an async check function for discord.py wait_for?
Python
2 votes279 views3 answers
Python logger get level
Python