网上找的:
! B/ P% B' Z+ l% l/ o* } V* d2 J
已知两天体赤经赤纬,如何计算角距离? - u+ f) C8 u# H1 |' N" i/ q
/ Q" H" d) C( ^0 {" @
如果你了解球坐标和向量代数, 化为球极坐标问题来解决思路比较简明清晰. (下面, 赤经---R.A. 赤纬---Dec), v. S9 P2 u3 G/ E
8 j& h/ W( l d每颗恒星在单位球上的投影的球坐标为(r, theta, phi), 这里r=1, theta=-(Dec-90), phi=R.A.
) }1 _( A9 x+ y, v" w _4 B8 a+ A8 p8 W2 j
变成投影点的直角坐标:
& F4 S9 r9 T6 \1 ^ s8 X5 `% l" ^1 Jx=sin(theta)cos(phi) / X) q, R+ I L$ b/ X8 t/ ]1 O0 r
y=sin(theta)sin(phi) 1 r; h+ u" P' _, i: U
z=cos(theta) ( v( k) z5 w1 p- C* W' ^/ a) H
. J) s, h+ `8 Z! g
这样可以把两颗星投影的角距离变为求两点间对球心的张角的简单问题. (点x0, y0, z0与点x1, y1, z1). 如果你用的角距离不大于180度(=Pi), 那么这个张角的值A可以直接求出: A=arc cos (x0*x1+y0*y1+z0*z1)
; d# f; E+ s( D0 y7 |# j$ b' N Y" l$ N# n7 Q" q
==============================================================
. R' q6 c0 A% D1 Y! o; b: @" q令:
+ u# K% q; i) y5 s0 M赤纬 = a1, a21 }! V8 k( [$ Y. g% v% V/ {9 f
赤经 = b1, b2) B7 Y% v( @* ^( T/ M, b5 Z
9 I8 s7 m5 h2 F2 ^' z0 B3 e3 m则:
2 _. L$ X6 }8 r# e. q5 z
' D1 N* p- E& m: q: H3 J0 k. J角距离 = arccos(sin(a1) * sin(a2) + cos(a1) * cos(a2) * cos(b2 - b1)) |