私の軌道上の各伝搬ステップのデカルト座標と速度ベクトルを取得しようとしています。私はOrbitalPy http://pythonhosted.org/OrbitalPy/を使って、古典的なケプラーの要素で軌道を生成しています。私はclass orbital.utilities.StateVector
から状態ベクトル(位置と速度の両方)を得ることができる必要がありますが、私はタイプエラー取得のドキュメントによるとPython OrbitalPyトレースバックエラーデカルト状態ベクトルケプラー要素からの位置と速度
:新しい()は正確に3つの引数(2与えられた)を取り
ここにコードがあります:
from scipy.constants import kilo
import orbital
from orbital import earth, KeplerianElements, Maneuver, plot, utilities
from orbital.utilities import Position, Velocity
import matplotlib.pyplot as plt
import numpy as np
#Orbit Setup
orbitPineapple = KeplerianElements.with_period(96 * 60, body=earth, i=(np.deg2rad(51.6)))
plot(orbitPineapple)
plt.show()
orbitPineapple
Out[23]: KeplerianElements(a=6945033.343911132,
e=0,
i=0.90058989402907408,
raan=0,
arg_pe=0,
M0=0.0,
body=orbital.bodies.earth,
ref_epoch=<Time object: scale='utc' format='jyear_str' value=J2000.000>)
prop1 = orbital.maneuver.PropagateAnomalyTo(M=1.00)
orbitX = orbitPineapple.apply_maneuver(prop1)
plot(orbitPineapple, title='Go Pineapple!')
plt.show()
orbital.utilities.StateVector(orbitPineapple)
TypeError Traceback (most recent call last)
<ipython-input-53-91fb5303082b> in <module>()
4 #print(orbital.utilities.StateVector.velocity(orbitPineapple))
5
----> 6 orbital.utilities.StateVector(orbitPineapple)
7 #orbital.utilities.StateVector.position(orbitPineapple())
8
TypeError: __new__() takes exactly 3 arguments (2 given)
位置ベクトルは、座標を他のシステムに変換しようとする人のために、AstroPyのGCRSであるGeocentric Equatorial Inertial(GEI)参照フレーム内にあります。 'gcrs。transform_to() 'はAstroPyで素晴らしい – Rose