1

で画像処理をしていscipy.ndimageます。リング状のオブジェクトが与えられた場合、その円周に「プロファイル」を生成したいと思います。プロファイルは、リングの周りのさまざまなポイントでの厚さの測定値、またはリングの「厚さ」に沿った平均信号のようなものです。

ndimage.mean最初に良いラベル画像を取得できれば使用できるように思えます。

私の指輪がこんな感じなら、

 A = array([[0, 0, 0, 0, 1, 0, 0, 0, 0, 0],
            [0, 0, 1, 1, 1, 1, 1, 1, 0, 0],
            [0, 0, 1, 1, 1, 1, 1, 1, 0, 0],
            [0, 0, 1, 1, 0, 0, 1, 1, 0, 0],
            [0, 0, 1, 1, 0, 0, 1, 1, 0, 0],
            [0, 0, 1, 1, 0, 0, 1, 1, 0, 0],
            [0, 0, 1, 1, 0, 0, 1, 1, 1, 0],
            [0, 0, 1, 1, 0, 0, 1, 1, 0, 0],
            [0, 0, 1, 1, 1, 1, 1, 1, 0, 0],
            [0, 0, 0, 0, 1, 0, 0, 0, 0, 0]])

numpy.mean( A, labels )where labelsisで「手段のプロファイル」を取得できました。

array([[0,  0,  0,  0,  2,  0,  0,  0,  0,  0],
       [0,  0,  1,  1,  2,  3,  4,  4,  0,  0],
       [0,  0,  1,  1,  2,  3,  4,  5,  0,  0],
       [0,  0, 16, 16,  0,  0,  5,  5,  0,  0],
       [0,  0, 15, 15,  0,  0,  6,  6,  0,  0],
       [0,  0, 14, 14,  0,  0,  7,  7,  0,  0],
       [0,  0, 13, 13,  0,  0,  8,  8,  8,  0],
       [0,  0, 13, 12,  0,  0,  9,  9,  0,  0],
       [0,  0, 12, 12, 11, 10,  9,  9,  0,  0],
       [0,  0,  0,  0, 11,  0,  0,  0,  0,  0]])

これで見落とされる補間要素がいくつかあるに違いありませんが、それは私が自分で思いつくことができるすべてです.

提案されたlabels画像を生成する方法はありますか? プロファイルを生成するためのより良い方法はありますか?

4

1 に答える 1

1

円 (または球) 上の標準的な確率分布は、フォン ミーゼス フィッシャー分布です。

Scipy は次のディストリビューションをサポートしています: http://docs.scipy.org/doc/scipy/reference/generated/scipy.stats.vonmises.html

したがって、fit 関数を使用して、データの最尤パラメーターを見つけることができるはずです。

于 2012-09-21T17:27:32.977 に答える