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

8 votes152 views6 answers
How to convert categorical data to numerical data in python
Python
1 vote232 views2 answers
Run flask app windows
Python Flask
12 votes155 views1 answer
Make a zero list python
Python
0 vote152 views1 answer
How to make a button open a new window in python
Python
33 votes154 views5 answers
Python check if string contains symbols
Python
1 vote154 views1 answer
Keep the user logged in even though user changes password django
Python
2 votes166 views3 answers
Combine list of dicts
Python
0 vote138 views1 answer
Linux pyspark select java version
Python
0 vote151 views1 answer
Open tar file pandas
Python
2 votes154 views1 answer
Get median of column pandas
Python