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

22 votes188 views2 answers
Python count repeated elements in a list
Python
4 votes288 views5 answers
Python booleans
Python
4 votes178 views2 answers
Check installed django version
Python Django
0 vote264 views1 answer
Django-filter field name greater than
Python
6 votes192 views5 answers
Python weekday
Python
1 vote198 views2 answers
Split a string with 2 char each in python
Python
0 vote201 views1 answer
Flask - store object directly in a session [duplicate]
Python
0 vote186 views1 answer
Merge all mp4 video files into one file python
Python
4 votes199 views1 answer
How to make a list string in python
Python
3 votes232 views4 answers
Csv to python dictionary
Python