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 vote171 views1 answer
Convert numpy array to tfrecord and back
Python
0 vote190 views1 answer
How to run mac terminal from python script
Python
0 vote172 views1 answer
How to write a table from 1 to 10 with for loop in fython in 3 lines
Python
0 vote283 views1 answer
Django 3 create async rest api
Python Django
0 vote172 views3 answers
Duck typing in python
Python
0 vote173 views1 answer
How to pre populate field flask wtforms
Python
1 vote183 views2 answers
Operators in python
Python
3 votes167 views2 answers
How to make python code run faster
Python
6 votes150 views3 answers
Pi in python
Python
25 votes177 views2 answers
Check if a list of strings is in a string python
Python