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 vote194 views1 answer
How to instal django cities
Python
-1 vote185 views1 answer
How to split a string into a dictionary in Python
Python
1 vote201 views1 answer
Python write txt utf8
Python
0 vote193 views1 answer
Make a function that accepts any nuber of arguments python
Python
12 votes225 views6 answers
Pandas drop duplicates
Python
0 vote200 views1 answer
Read a csv and plot in python
Python
4 votes201 views2 answers
Tty escape
Python
0 vote219 views3 answers
Removing duplicates from django models data
Python
7 votes185 views1 answer
Program to find even numbers in python
Python
0 vote202 views1 answer
Unique mark boolean django model field
Python