使用库pyephem,我想找到一种方法来绘制表示卫星地面轨迹的漂亮绘图(经度、纬度)曲线。我已经计算了 (longitude,latitude) 但是当经度超过 +180° 时,计算的下一个值是例如 -178°,这样就绘制了一个线段:这就是糟糕的数字。
在我目前拥有的代码和获得的数字下方。
currentDate = date.datetime(2018,12,1,0,0,0);
for i in range(nPoints):
iss.compute(currentDate)
# compute latitude
posLat[i] = iss.sublat*(180/math.pi)
# compute longitude
posLong[i] = iss.sublong*(180/math.pi)
currentDate = currentDate + date.timedelta(seconds = (time3orbits/nPoints*3600))
print currentDate, posLong[i], posLat[i]
plt.plot(posLong,posLat)
plt.show()
这是我得到的图(绕地球 3 圈):
如何防止绘制的线段?问候
iss.sublat
和iss.sublong
分别是 iss 位置的当前纬度和经度(从一个currentDate
固定的)。
牛魔王的故事
相关分类