2

マウナケアから特定の時間に特定の RA/Dec で (有名でない) 星の alt/az を決定しようとしています。pyephem を使用してこれらのパラメーターを計算しようとしていますが、結果の alt/az は他のソースと一致しません。Keck による HAT-P-32 の計算は次のとおりです。

import ephem
telescope = ephem.Observer()
telescope.lat =  '19.8210'
telescope.long = '-155.4683'
telescope.elevation = 4154
telescope.date = '2013/1/18 10:04:14'
star = ephem.FixedBody()
star._ra = ephem.degrees('02:04:10.278')
star._dec = ephem.degrees('+46:41:16.21')
star.compute(telescope)
print star.alt, star.az

これは を返しますが-28:43:54.0 73:22:55.3、Stellarium によると、適切な alt/az は次のようになります62:26:03 349:15:13。私は何を間違っていますか？

EDIT : 以前は逆だった緯度と経度を修正しました。

4

2 に答える 2

1

まず、後方に長い緯度があります。次に、文字列を 16 進形式で指定する必要があります。第三に、RA を度ではなく時間で指定する必要があります。

import ephem
telescope = ephem.Observer()
# Reversed longitude and latitude for Mauna Kea
telescope.lat =    '19:49:28' # from Wikipedia
telescope.long = '-155:28:24'
telescope.elevation = 4154.
telescope.date = '2013/1/18 00:04:14'
star = ephem.FixedBody()
star._ra  = ephem.hours('02:04:10.278') # in hours for RA
star._dec = ephem.degrees('+46:41:16.21')
star.compute(telescope)

このようにして、次のようになります。

>>> print star.alt, star.az
29:11:57.2 46:43:19.6
于 2013-02-22T15:37:10.780 に答える
0

PyEphem は常に時刻に UTC を使用するため、プログラムがどこで実行されても同じように動作し、同じ出力が得られます。ローカル タイム ゾーンを使用する代わりに、使用している日付を UTC に変換するだけで、結果は Stellarium とほぼ一致します。使用する：

telescope.date = '2013/1/18 05:04:14'

結果はこの alt/az です:

62:27:19.0 349:26:19.4

残りのわずかな違いがどこから来るのかを知るには、2 つのプログラムが計算の各ステップをどのように処理するかを調べる必要があります。しかし、これはあなたを十分に近づけますか？

于 2013-02-23T15:56:39.313 に答える