回答編集履歴
1
質問者の要望
test
CHANGED
|
@@ -3,6 +3,18 @@
|
|
|
3
3
|
ISSは約９０分で地球を一周するので、１分おきに９０分間の方位角と高度をプロットしました。
|
|
4
4
|
|
|
5
5
|
satellites（ISS [ZARYA]） の取得を変更しました。
|
|
6
|
+
|
|
7
|
+
|
|
8
|
+
|
|
9
|
+
◎以下追加
|
|
10
|
+
|
|
11
|
+
変数gtypeに以下のように設定すると、それぞれのグラフを出力します。
|
|
12
|
+
|
|
13
|
+
az-alt : x軸=方位角 y軸=高度
|
|
14
|
+
|
|
15
|
+
time-alt : x軸=時刻 y軸=高度（先頭３０個のみ）
|
|
16
|
+
|
|
17
|
+
|
|
6
18
|
|
|
7
19
|
```python
|
|
8
20
|
|
|
@@ -50,25 +62,47 @@
|
|
|
50
62
|
|
|
51
63
|
df = pd.read_csv("result.csv")
|
|
52
64
|
|
|
53
|
-
plt
|
|
65
|
+
gtype = 'time-alt' # az-alt = x:方位角-y:高度 、 time-alt = x:time-y:高度
|
|
54
66
|
|
|
55
|
-
|
|
67
|
+
if gtype == 'az-alt':
|
|
56
68
|
|
|
57
|
-
|
|
69
|
+
plt.axis([0, 360, -90, 90])
|
|
58
70
|
|
|
59
|
-
|
|
71
|
+
x = df['az']
|
|
60
72
|
|
|
61
|
-
|
|
73
|
+
y = df['alt']
|
|
62
74
|
|
|
63
|
-
plt.
|
|
75
|
+
plt.scatter(x, y)
|
|
64
76
|
|
|
65
|
-
plt.
|
|
77
|
+
plt.plot(x, y)
|
|
66
78
|
|
|
79
|
+
plt.xlabel('方位角')
|
|
80
|
+
|
|
81
|
+
plt.ylabel('高度')
|
|
82
|
+
|
|
67
|
-
plt.title('ISS (ZARYA)の軌跡：2020/11/01 00:00:00(UTC) から90分間(60秒毎)')
|
|
83
|
+
plt.title('ISS (ZARYA)の軌跡：2020/11/01 00:00:00(UTC) から90分間(60秒毎)')
|
|
84
|
+
|
|
85
|
+
elif gtype == 'time-alt':
|
|
86
|
+
|
|
87
|
+
x = df['time'][0:30]
|
|
88
|
+
|
|
89
|
+
y = df['alt'][0:30]
|
|
90
|
+
|
|
91
|
+
plt.subplots_adjust(bottom=0.5)
|
|
92
|
+
|
|
93
|
+
plt.xticks(rotation=90)
|
|
94
|
+
|
|
95
|
+
plt.scatter(x, y)
|
|
96
|
+
|
|
97
|
+
plt.plot(x, y)
|
|
98
|
+
|
|
99
|
+
plt.xlabel('時刻')
|
|
100
|
+
|
|
101
|
+
plt.ylabel('高度')
|
|
102
|
+
|
|
103
|
+
plt.title('ISS (ZARYA)の時刻・高度：2020/11/01 00:00:00(UTC) から30個(60秒毎)')
|
|
68
104
|
|
|
69
105
|
plt.show()
|
|
70
|
-
|
|
71
|
-
|
|
72
106
|
|
|
73
107
|
```
|
|
74
108
|
|
|
@@ -76,4 +110,12 @@
|
|
|
76
110
|
|
|
77
111
|
結果は以下のようになりました。
|
|
78
112
|
|
|
113
|
+
gtypeがaz-altの時
|
|
114
|
+
|
|
79
|
-
|
|
115
|
+
|
|
116
|
+
|
|
117
|
+
|
|
118
|
+
|
|
119
|
+
gtypeがtime-altの時
|
|
120
|
+
|
|
121
|
+
|