地球上の特定のオブザーバー (経度、緯度、時間) と特定の銀河座標 (GLON、GLAT) について、PyEphem を使用して空の対応する (alt、az) ポイントを計算するにはどうすればよいですか?
質問する
2220 次
1 に答える
3
PyEphem の現在の動作方法を考えると、質問に答えるには 2 つのステップがあります。まず、代わりに、銀河座標のペアを赤道 RA/dec 座標のペアに変換する必要があります。
import ephem
# Convert a Galactic coordinate to RA and dec
galactic_center = ephem.Galactic(0, 0)
eq = ephem.Equatorial(galactic_center)
print 'RA:', eq.ra, 'dec:', eq.dec
→ RA: 17:45:37.20 dec: -28:56:10.2
これらの座標は、ウィキペディアが銀河の中心について示しているものにかなり近いです。
通常の赤経/赤緯座標が得られたので、現在の空の位置を把握する必要があります。PyEphem は、星や惑星などの天体の「天体」についてのみ知っているライブラリの上に構築されているため、その場所に偽の「星」を作成し、その方位角と高度を尋ねるだけです。
# So where is that RA and dec above Boston?
# Pretend that a star or other fixed body is there.
body = ephem.FixedBody()
body._ra = eq.ra
body._dec = eq.dec
body._epoch = eq.epoch
obs = ephem.city('Boston')
obs.date = '2012/6/24 02:00' # 10pm EDT
body.compute(obs)
print 'Az:', body.az, 'Alt:', body.alt
→ Az: 149:07:25.6 Alt: 11:48:43.0
そして、その夜遅くにボストンの天空図を見ることで、この答えが妥当であることを確認できます: 射手座 — 銀河の中心の位置 — は空の南東の縁の上にちょうど昇っています。 149°、および 11° のようなまだ低い高度の場合。
于 2012-06-24T22:12:25.473 に答える