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 vote342 views1 answer
Django file field from base64
Python Django
1 vote347 views1 answer
Show jpg in jupyter notebook
Python
1 vote342 views1 answer
Scikit learn split data set
Python
0 vote321 views1 answer
Iniciar un projecto de python con pyenv
Python
4 votes282 views1 answer
Datetime has no attribute now
Python
0 vote341 views1 answer
Dictionary Get All Values
Python
0 vote356 views1 answer
Python inject into process
Python
54 votes396 views6 answers
Conda create environment python 3.6
Python
3 votes315 views3 answers
Pytest
Python
10 votes371 views1 answer
Https://docs.djangoproject.com/en/2.0/howto/static-files/
Python Django