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 vote126 views1 answer
Open google chrome from python
Python
0 vote180 views1 answer
Return a sorted copy of the list. Does not modify original list.
Python
0 vote173 views1 answer
Compile python folder
Python
10 votes169 views2 answers
Pygame get mouse position
Python
6 votes213 views4 answers
Flask get with parameters
Python Flask
0 vote223 views1 answer
Max value indices
Python
0 vote202 views1 answer
Area of a right angle triangle in python
Python
0 vote159 views1 answer
Python ssl module is not available
Python
0 vote176 views1 answer
Is root node an internal node
Python
0 vote176 views1 answer
Comment voir les requêtes SQL brutes exécutées par Django?
Python