pythonにて、計算結果をグラフに表す方法について

Question

前提・実現したいこと

人工衛星の軌道計算を行っており、時刻に対する高度の推移のグラフを作りたいです。

該当のソースコード

from skyfield.api import Topos, load, EarthSatellite
from datetime import datetime, timedelta
from pytz import timezone
ts = load.timescale()
tokyo = Topos('35.65870 N', '139.74539 E')
line1 = '1 25544U 98067A   20226.06311231  .00000634  00000-0  19556-4 0  9992'
line2 = '2 25544  51.6462  66.9823 0001637  29.7739 108.2756 15.49160058240839'
iss = EarthSatellite(line1, line2, 'ISS (ZARYA)', ts)
t0 = ts.utc(2020, 11, 1, 0, 0, 0)
with open("resalt.txt", mode='w') as f:
    for s in range(600):
        t1 = ts.utc(t0.utc_datetime() + timedelta(seconds=s))
        alt, az, distance = (iss - tokyo).at(t1).altaz()
        print('時刻:', t1.utc_strftime("%Y %b %d %H:%M:%S"), '高度:{0:.1f} 度'.format(alt.degrees), file=f)

試したこと

計算結果をテキストファイルとして保存し、それを読み込みmatplotlibでグラフ化しようとしています。

計算結果のテキストは下記のようにテキストファイルに書き込まれています。
時刻: 2020 Nov 01 00:32:51 高度:0.0 度
時刻: 2020 Nov 01 00:32:52 高度:0.1 度
時刻: 2020 Nov 01 00:32:53 高度:0.1 度
時刻: 2020 Nov 01 00:32:54 高度:0.2 度

カンマ区切りの方がやりやすいのでしょうか

Answer 1

skyfieldを初めて知り面白そうだったので書いてみました。
ISSは約９０分で地球を一周するので、１分おきに９０分間の方位角と高度をプロットしました。
satellites（ISS [ZARYA]） の取得を変更しました。

◎以下追加
変数gtypeに以下のように設定すると、それぞれのグラフを出力します。
az-alt : x軸=方位角 y軸=高度
time-alt : x軸=時刻 y軸=高度（先頭３０個のみ）

from skyfield.api import Topos, load
from datetime import timedelta
import matplotlib.pyplot as plt
import pandas as pd
import japanize_matplotlib

ts = load.timescale()
tokyo = Topos('35.65870 N', '139.74539 E')
satellites = load.tle('http://celestrak.com/NORAD/elements/stations.txt')
iss = satellites['ISS (ZARYA)']
t0 = ts.utc(2020, 11, 1, 0, 0, 0)
with open("result.csv", mode='w') as f:
    print('time,alt,az', file=f)
    for s in range(0, 5400, 60):
        t1 = ts.utc(t0.utc_datetime() + timedelta(seconds=s))
        alt, az, distance = (iss - tokyo).at(t1).altaz()
        t = t1.utc_strftime("%Y/%m/%d %H:%M:%S")
        alt_degree = alt.degrees  # 高度（度）
        az_degree = az.degrees    # 方位角（度）
        print('{},{},{}'.format(t, alt_degree, az_degree), file=f)

df = pd.read_csv("result.csv")
gtype = 'time-alt'  # az-alt = x:方位角-y:高度 、 time-alt = x:time-y:高度
if gtype == 'az-alt':
    plt.axis([0, 360, -90, 90])
    x = df['az']
    y = df['alt']
    plt.scatter(x, y)
    plt.plot(x, y)
    plt.xlabel('方位角')
    plt.ylabel('高度')
    plt.title('ISS (ZARYA)の軌跡：2020/11/01 00:00:00(UTC) から90分間(60秒毎)')
elif gtype == 'time-alt':
    x = df['time'][0:30]
    y = df['alt'][0:30]
    plt.subplots_adjust(bottom=0.5)
    plt.xticks(rotation=90)
    plt.scatter(x, y)
    plt.plot(x, y)
    plt.xlabel('時刻')
    plt.ylabel('高度')
    plt.title('ISS (ZARYA)の時刻・高度：2020/11/01 00:00:00(UTC) から30個(60秒毎)')
plt.show()

結果は以下のようになりました。
gtypeがaz-altの時

gtypeがtime-altの時

Comments

[sa-to]

ありがとうございます！
時間と高度のグラフなどににしたい場合は、
x = df['az']
y = df['alt']
この部分を変更すれば良いのでしょうか？