（知识点）天体的时角计算公式，你知道几个？

t = LST - α (1)

上式中，LST为当地恒星时，α为天体赤经。

#

当地恒星时可以通过经度与格林威治恒星时联系起来。 公式如下：

#

LST= 商品及服务税 + \ (2) #

格林威治平均恒星时等于地球自转角（ERA）与二分差（EE）之和。 GMST的公式如下（直到三次项）： #

GMST=\theta(UT1)+0".+4613".+1".^2 - 0".^3 (3)

二部差分仅保留交角章动第一项，可写为：

EE=\Delta \psi\cos \ (4) #

式(3)中第一项为地球自转角，计算公式为： #

\theta(UT1)=2\pi(0. + 1.d) (5) #

上式中赤纬角计算公式，d为从标准历元J2000.0开始的世界时（UT1时标）儒略日数。 #

有了这么一堆公式，就可以算出任意天体的时角了。 #

---------------------------- 太长的分割线 ------------ --- ------------------

如果你的机器上有它并下载了软件包，你可以简单地使用以下程序来计算： #

from astropy.coordinates import EarthLocation, get_body
from astropy.time import Time
lon = '125d19m22.0s'  # 经度
lat = '43d52m48.0s'   # 纬度
alt = '206m'          # 海拔
time = Time ('2007-01-01 00:00:00')   # UTC时间
site = EarthLocation ( lon, lat, alt )
venus = get_body('venus', time = time, location = site)
lmst = time.sidereal_time ('mean', longitude = lon)
ha = lmst - venus.ra
print (ha) 
#
 

只需自行修改经度、纬度、海拔和时间即可。 #

最终输出的小时角可能小于0赤纬角计算公式第一考试网，所以手动添加24小时即可。 #

或者如果你想计算其他主要行星的小时角度，只需将程序中所有金星搜索替换为其他主要行星的英文名称即可。 #