1

Pythonとmatplotlibを使用して、可視衛星を表すいくつかの値を極座標グラフにプロットしたいと思います。matplotlibの例に従ってコードを記述しました。極座標チャートが表示されていますが、指定されたポイントはプロットされていません。

import matplotlib
from matplotlib.pyplot import figure, show, rc, grid
from math import pi

# radar green, solid grid lines
rc('grid', color='#316931', linewidth=1, linestyle='-')
rc('xtick', labelsize=15)
rc('ytick', labelsize=15)

# force square figure and square axes looks better for polar, IMO
width, height = matplotlib.rcParams['figure.figsize']
size = min(width, height)
# make a square figure
fig = figure(figsize=(size, size))
ax = fig.add_axes([0.1, 0.1, 0.8, 0.8], polar=True, axisbg='#d5de9c')

# Satellite info [PRN, E, Az, Ss, Used (1 yes, 0 no)]
sat = [ [1, 62, 255, 46, 1],
        [14, 62, 26, 46, 1],
        [31, 42, 158, 36, 1],
        [22, 40, 76, 50, 1],
        [11, 29, 308, 0, 0],
        [19, 26, 243, 36, 1],
        [3, 13, 217, 0, 0],
        [18, 10, 93, 0, 0],
        [20, 6, 291, 0, 0],
        [5, 1, 72, 0, 0],
        [122, 43, 216, 0, 0],
        [135, 47, 203, 43, 0] ]

for index in (0, len(sat)-1):
    if(sat[index][4]>0):
        ax.plot(sat[index][2], sat[index][1], color='green', marker='s', markersize=12)
    else:
        ax.plot(sat[index][2], sat[index][1], color='gray', marker='s', markersize=12)

ax.set_rmax(2.0)
grid(True)

ax.set_title("Visible satellites", fontsize=20)
show()

私は何が間違っているのですか?

4

2 に答える 2

1

次のように、satのリストを反復処理する方が簡単です。

for s in sat:
    if(s[4]>0):
        ax.plot(s[2], s[1],color='green', marker='s', markersize=5)
    else:
        ax.plot(s[2], s[1],color='gray', marker='s', markersize=5)

それに加えて、set_rmaxを使用して方位角を2未満に制限しています。これにより、1つの衛星のみが表示され、コメントを外してすべてを表示します。

編集:リストを直接解凍すると、読みやすさがさらに向上する可能性があります。

for (PRN, E, Az, Ss, Used) in sat:
    if(Used>0):
        ax.plot(Ss, Az,color='green', marker='s', markersize=5)
    else:
        ax.plot(Ss, Az,color='gray', marker='s', markersize=5)
于 2012-10-18T11:24:40.790 に答える
1

この行で:

for index in (0, len(sat)-1):

インデックスは、値のペアに対してのみ実行されます。あなたは範囲を意味しましたfor index in range(0, len(sat))か?

于 2012-10-18T11:20:55.177 に答える