Pythonのサークルからデータを取る

Question

角度によってリングの強さがどのように変化するかを調べています。私がやりたい何

はそのドーナツの中心部内から値の円を取り、角度対してプロットしたものである。ここでは画像の例があります。私が現在行っているのは、scipy.ndimage.interpolation.rotateを使ってリングを半径方向にスライスし、2つのピークの最大値を抽出し、それらの角度をプロットすることです。

crop = np.ones((width,width)) #this is my image 
    slices = np.arange(0,width,1) 
    stack = np.zeros((2*width,len(slices))) 
    angles = np.linspace(0,2*np.pi,len(crop2)) 

    for j in range(len(slices2)): # take slices 
      stack[:,j] = rotate(crop,slices[j],reshape=False)[:,width] 

しかし、これは私が実際に探していることではないと思います。私は、私が望むデータをどのように抽出するのかとほとんど苦労しています。私はまたこのようなマスクを適用しようとしました。

画像へ

、その後私は、正しい順序で、マスク内の値を取得する方法がわからない（すなわち角度0増加のために - 。2PI）を

その他アイデアは大きな助けになるでしょう！

Answer 1

私は正しさの検証を支援するために、異なる入力画像を作っ：画像内の円形の経路からのサンプル値に

import numpy as np 
import scipy as sp 
import scipy.interpolate 
import matplotlib.pyplot as plt 

# Mock up an image. 
W = 100 
x = np.arange(W) 
y = np.arange(W) 
xx,yy = np.meshgrid(x,y) 

image = xx//5*5 + yy//5*5 
image = image/np.max(image) # scale into [0,1] 

plt.imshow(image, interpolation='nearest', cmap='gray') 
plt.show() 

を、我々は最初のインターポレータを構築し、我々は、任意の場所にアクセスしたいので。私たちはまたそれをより高速にベクトル化します。
次に円のパラメトリック定義x(t) = sin(t), y(t) = cos(t)を使用して円の円周上に座標Nの座標を生成します。
Nは、円周の少なくとも2倍（ナイキストシャノンサンプリング定理）でなければなりません。

interp = sp.interpolate.interp2d(x, y, image) 
vinterp = np.vectorize(interp) 

for r in (15, 30, 45): # radii for circles around image's center 
    xcenter = len(x)/2 
    ycenter = len(y)/2 
    arclen = 2*np.pi*r 
    angle = np.linspace(0, 2*np.pi, arclen*2, endpoint=False) 
    value = vinterp(xcenter + r*np.sin(angle), 
        ycenter + r*np.cos(angle)) 
    plt.plot(angle, value, label='r={}'.format(r)) 

plt.legend() 
plt.show()