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 vote315 views1 answer
Python dash bootstrap buttons with icons
Python
6 votes276 views2 answers
How to limit a command to a permission in discord.py
Python
4 votes269 views1 answer
Np.random.uniform
Python
134 votes277 views6 answers
Python install pip
Python
0 vote276 views1 answer
Modbusfc03 python
Python
-1 vote259 views1 answer
What will be the output of the following python code? print("abc. def".capitalize())
Python
0 vote269 views1 answer
Django data from many to many field in template
Python
0 vote282 views1 answer
Waitress serve
Python
0 vote252 views1 answer
Multiply every nth element
Python
1 vote302 views1 answer
How to open dialog box to select folder in python
Python