夜空をのぞく仮想望遠鏡のようなものを書きたいと思っています。
今まで私は星のカタログを持っていて、それらを平面に投影してモック写真を撮りたいと思っています。
射影は であると推測します。gnomonic projectionこれは、こことここにあります。
2番目のリンクでは、星のピクセル位置の計算に関するalg。
Forward:
Define
scale: number of pixels per degree in the map
alpha, delta: Equatorial coordinates of a given position
alpha0, delta0: Equatorial coordinates of the map center
A = cos(delta) x cos(alpha - alpha0)
F = scale x (180/pi)/[sin(delta0) x sin(delta) + A x cos(delta0)]
then the pixel coordinates in the image are
LINE = -F x [cos(delta0) x sin(delta) - A x sin(delta0)]
SAMPLE = -F x cos(delta) x sin(alpha - alpha0)
Reverse:
Define
X = SAMPLE/(scale x 180/pi)
Y = LINE/(scale x 180/pi)
D = arctan[(X^2 + Y^2)^0.5]
B = arctan(-X/Y)
XX = sin(delta0) x sin(D) x cos(B) + cos(delta0) x cos(D)
YY = sin(D) x sin(B)
then the right ascension and declination are
alpha = alpha0 + arctan(YY/XX)
delta = arcsin[sin(delta0) x cos(D) - cos(delta0) x sin(D) x cos(B)]
NOTE: The arctangent functions for B and alpha must be four-quadrant arctangents.
ただし、角度を deg にするか rad にするか、SAMPLE と LINE の意味がわかりません。
また、ノーモニック プロジェクションの使用についても確信が持てません。
ヘルプやディスカッションは大歓迎です。