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

2 votes134 views2 answers
Inspect rows in dictionary pthon
Python
0 vote133 views1 answer
Filter nulla values only pandas
Python
1 vote132 views1 answer
Python type hint for a string
Python
14 votes160 views6 answers
Binary search python
Python
0 vote147 views1 answer
How to handle response from tkinter messagebox.askquestion() function in Python
Python
21 votes136 views2 answers
How to stop running code in python
Python
0 vote152 views1 answer
Field 'id' expected a number but got <django.db.models.fields.related_descriptors.create_forward_many_to_many_manager.<locals>.ManyRelatedManager object at 0x7fb02044de80>
Python
0 vote146 views1 answer
Get list of all document in django-elasticsearch-dsl
Python
0 vote141 views1 answer
Python read from stdin pipe
Python
0 vote183 views1 answer
Nested dictionary to dataframe
Python