标签：python matplotlib plot

使用图书馆pyephem,我想找到一种方法来做一个很好的情节(经度,纬度)曲线,代表卫星的地面轨迹.我已经计算了(经度,纬度),但是当经度超过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分别是当前经常位置的经度和纬度(来自currentDate固定).

解决方法:

您可以将符号更改时的数据数组从正数分割为负数,然后绘制单个子数组：

orbits = np.argwhere(
    (data[:-1, 0] * data[1:, 0] < 0)
    & (data[1:, 0] < 0)
).ravel() + 1

for orbit in np.split(data, orbits, axis=0):
    plt.plot(orbit[:, 0], orbit[:, 1])

标签：python,matplotlib,plot
来源： https://codeday.me/bug/20190710/1424463.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。