例子: 计算 梅西耶 M13 在 1998年8月10号在23点10分的 ALT 和AZ ,地点是 UT,英国伯明翰。M13的 RA 和 DEC 是:
RA = 16 h 41.7 min
DEC = 36 d 28 min
一、准备工作
英国伯明翰的经度和纬度作为;
LAT = 52 d 30 min North
LONG = 1 d 55 min West
转换成十进制
RA = 16 h 41.7 min = 16 + 41.7/60 = 16.695 hrs
DEC = 36 d 28 min = 36 + 28/60 = 36.466667 degs
Time = 2310 hrs = 23 + 10/60 = 23.166667 hrs
LAT = 52 d 30 min North = 52 + 30/60 = 52.5 degs
LONG = 1 d 55 min West = -(1 + 55/60) = -1.9166667 degs
一直保留小数点后面的数字,直至计算完成。 经纬度的 西方为负,东为正。 再把RA换算成度,
24H 转360度, 1h 等于15°
RA = 16.695 * 15 = 250.425 degrees
LST = 100.46 + 0.985647 * d + long + 15*UT
d J2000的天数
UT 时间 小时为单位
long is your longitude in decimal degrees, East positive.
Add or subtract multiples of 360 to bring LST in range 0 to 360
longitude 1 degree 55 minutes west).
举例 1998年8月10号 23点10分 Birmingham UK
I know that UT = 23.166667
d = -508.53472
long = -1.9166667 (西方为负数)
so
LST = 100.46 + 0.985647 * d + long + 15*UT
= 100.46 + 0.985647 * -508.53472 - 1.9166667 + 15 * 23.166667
= -55.192383 degrees
= 304.80762 degrees
LST 最后是值需要再0-360范围内,需要看看情况加一个360 或减
3、时角 HA
HA = LST - RA
例子
RA = 250.425 degs
LST = 304.80762
则
HA = LST - RA
= 304.80762 - 250.425
= 54.382617 degs
4、HA、DEC 计算出ALT、AZ
公式:
sin(ALT) = sin(DEC)*sin(LAT)+cos(DEC)*cos(LAT)*cos(HA)
ALT = asin(ALT)
sin(DEC) - sin(ALT)*sin(LAT)
cos(A) = ---------------------------------
cos(ALT)*cos(LAT)
A = acos(A)
If sin(HA) is negative, then AZ = A, otherwise
AZ = 360 - A
问题来了 S输出的 20位数据中 时间的含义 这是LST吗 ? 说明是说 从1970 年1月1号开始的 那是不是 LST= (2000年-1970年)- 2012年12月25号 TIME (8 bytes,integer): current time on the client computer in microseconds
since 1970.01.01 UT. Currently unused.
1、协议说明是这样的
RA (4 bytes,unsigned integer): right ascension of the telescope (J2000)
a value of 0x100000000 = 0x0 means 24h=0h,
a value of 0x80000000 means 12h
一秒等
详情回复
发表于 2012-12-27 17:54
1、协议说明是这样的
RA (4 bytes,unsigned integer): right ascension of the telescope (J2000)
a value of 0x100000000 = 0x0 means 24h=0h,
a value of 0x80000000 means 12h