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 votes286 views1 answer
Numpy where
Python
0 vote318 views1 answer
Calculate Euclidean Distance in Python using norm()
Python
0 vote299 views2 answers
For schleife python
Python
0 vote295 views1 answer
Python file back to beginning
Python
0 vote271 views1 answer
Python codes for counting the occurrence of a letters in dictionary excluding digits
Python
0 vote317 views1 answer
How to convert csv columns to text python
Python
5 votes300 views2 answers
Pyaudio not installing ubuntu
Python
1 vote287 views1 answer
Matplotlib change bar color under threshold
Python
2 votes296 views4 answers
Python dict get random key
Python
82 votes286 views6 answers
How to create a for loop in python
Python