我试图按照维汶提反问题的方法实现维基上这里描述的内容。
问题在于lambda根本无法收敛。如果我尝试迭代公式序列,它的值就会保持不变,我真的不确定为什么。也许我只是对一个明显的问题盯得太久而看不到。
需要指出的是,我是Python的新手,还在学习语言,因此我不确定是否使用了错误的语言可能导致问题,或者我在执行某些计算时有错误。我只是似乎找不到公式中的任何错误。
基本上,我尽可能地按照维基文章的格式编写了代码,结果如下:
import math
# Length of radius at equator of the ellipsoid
a = 6378137.0
# Flattening of the ellipsoid
f = 1/298.257223563
# Length of radius at the poles of the ellipsoid
b = (1 - f) * a
# Latitude points
la1, la2 = 10, 60
# Longitude points
lo1, lo2 = 5, 150
# For the inverse problem, we calculate U1, U2 and L.
# We set the initial value of lamb = L
u1 = math.atan( (1 - f) * math.tan(la1) )
u2 = math.atan( (1 - f) * math.tan(la2) )
L = (lo2 - lo1) * 0.0174532925
lamb = L
while True:
sinArc = math.sqrt( math.pow(math.cos(u2) * math.sin(lamb),2) + math.pow(math.cos(u1) * math.sin(u2) - math.sin(u1) * math.cos(u2) * math.cos(lamb),2) )
cosArc = math.sin(u1) * math.sin(u2) + math.cos(u1) * math.cos(u2) * math.cos(lamb)
arc = math.atan2(sinArc, cosArc)
sinAzimuth = ( math.cos(u1) * math.cos(u2) * math.sin(lamb) ) // ( sinArc )
cosAzimuthSqr = 1 - math.pow(sinAzimuth, 2)
cosProduct = cosArc - ((2 * math.sin(u1) * math.sin(u2) ) // (cosAzimuthSqr))
C = (f//16) * cosAzimuthSqr * (4 + f * (4 - 3 * cosAzimuthSqr))
lamb = L + (1 - C) * f * sinAzimuth * ( arc + C * sinArc * ( cosProduct + C * cosArc * (-1 + 2 * math.pow(cosProduct, 2))))
print(lamb)
如上所述,问题在于变量“lamb”(lambda)的值不会变小。我甚至尝试将我的代码与其他实现进行比较,但它们看起来几乎一样。
我在这里做错了什么? :-)
谢谢大家!
lamb
应该收敛到哪个值? - Jivanfrom __future__ import division
。 - Jonathan Daviessin α = blabla
时,它并不意味着sinAlpha = blabla
,而是意味着α = arcsin blablabla
。 - JivanL = (lo2 - lo1) * 0.0174532925
中,为什么要乘以* 0.017...
?维基百科的文章中没有提到这一点。 - Jivan