2

matplotlibを使用して、一部のデータをグリッドに適合させ、極座標投影としてプロットしています。以下の例のようなもの。ただし、プロットのエッジが0/360度で交わる場所を滑らかにしたいと思います。誰もが私がこれを行う方法を知っていますか？

from pylab import *
import random

x = linspace(0, 360, 361).astype(int)
x = x*pi/180
y = linspace(0.05, 0.5, 800)
xgrid, ygrid = meshgrid(x, y)


baz = []
for c in range(2000): baz.append(random.randint(0,360))


freq = rand(len(baz))
pwr = rand(len(baz))

zgrid = griddata(baz,freq,pwr, xgrid, ygrid)
subplot(111,  polar=True)
pcolormesh(xgrid, ygrid, zgrid)
show()

また、私が使用しているデータには、griddataによって作成されたマスクのためにギャップがあります（上記のようにgriddataを使用しますが、ループ内の多くのグリッドを合計します）。不足しているセグメントを埋めたいのですが（添付の図を参照）、これを行う方法を知っている人はいますか？

ここに画像の説明を入力してください ありがとうデイブ

4

1 に答える 1

1

0/360度の位置でどのグリッドが結合するかがわかっている場合は、それらを連結してスプライン補間（scipy補間）を実行できます。2番目の問題についてはわかりませんが、極座標でグリッドを作成するのはどうですか？これで問題は解決しますか？

敬具

于 2012-11-23T15:16:01.043 に答える