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 vote155 views1 answer
How to make a button open a new window in python
Python
2 votes139 views2 answers
Get rid of unnamed column pandas
Python
0 vote158 views1 answer
Text generate gpt 2 huggingface
Python
0 vote150 views1 answer
Scatter matrix plot
Python
8 votes155 views6 answers
How to convert categorical data to numerical data in python
Python
9 votes168 views3 answers
Python 2.7 datetime to timestamp
Python
0 vote166 views3 answers
Pandas assign multiple columns at once
Python
0 vote157 views1 answer
Video timestamp opencv python
Python
0 vote140 views1 answer
Dataframe look at every second column
Python
0 vote157 views1 answer
Hewwo world
Python