计算从地球上某个位置看到的拉格朗日点的RA DEC/AZ EL

如果可以将任意(地心或日心)位置向量转换为RA/DEC,则解决方案很简单,因为可以相对容易地计算拉格朗日点的位置向量。

如果你有太阳的地心向量, Rsun ,可以执行以下操作:

因为L1,L2和L3在日地线上,它们只是不同比例的版本 Rsun公司 . 这里有一些很好的近似方法,因为太阳的质量比地球的质量大得多(精确的公式要复杂得多):

L1 = Rsun * (m / (3*M))**(1/3)

L2 = -Rsun * (m / (3*M))**(1/3)

L3 = Rsun * (2 + (5*m / (12*M))

质量单位无关紧要,所以我们可以用地球质量作为单位 m == 1 M == 333000

L4和L5点是与太阳和地球等边三角形的角。所以你可以通过旋转得到它们 Rsun公司

# Get the three components of the original vector
x,y,z = Rsun 

# Compute the rotated vector
L4 = [
  x * cos(60°) - y * sin(60°),
  x * sin(60°) + y * cos(60°),
  z 
]