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 vote131 views1 answer
Maximum count of replacements in python
Python
0 vote133 views1 answer
ImportError: cannot import name 'StringIO'
Python
7 votes153 views1 answer
From string to time python dataframe
Python
1 vote129 views1 answer
Isistance exmaple
Python
1 vote136 views1 answer
Django loop index
Python Django
0 vote165 views1 answer
Parser.add_argument array python
Python
6 votes147 views3 answers
Python how to save a Seaborn plot into a file
Python
0 vote120 views1 answer
Setting flask app environment variable
Python
0 vote140 views1 answer
What data types store value by index in python
Python
42 votes129 views4 answers
Replace column values pandas
Python