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

1 vote76 views1 answer
Python add word to list
Python
0 vote106 views1 answer
Keep calm and carry on memes
Python
5 votes112 views1 answer
Python undefine variable
Python
0 vote174 views1 answer
Subprocess.popen no output
Python
0 vote114 views1 answer
Python check file format
Python
38 votes110 views2 answers
Python check if path does not exist
Python
23 votes110 views6 answers
Python even or odd
Python
0 vote120 views1 answer
Unpack a sequence into variables python
Python
1 vote119 views1 answer
Conda specify multiple channels
Python
6 votes112 views1 answer
Send data through tcp sockets python
Python