二十八、在 Basemap 中绘制坐标

欢迎阅读另一个 Matplotlib Basemap 教程。 在本教程中，我们将介绍如何绘制单个坐标，以及如何在地理区域中连接这些坐标。

首先，我们将从一些基本的起始数据开始：

  1. from mpl_toolkits.basemap import Basemap
  2. import matplotlib.pyplot as plt
  4. m = Basemap(projection='mill',
  5.             llcrnrlat = 25,
  6.             llcrnrlon = -130,
  7.             urcrnrlat = 50,
  8.             urcrnrlon = -60,
  9.             resolution='l')
  11. m.drawcoastlines()
  12. m.drawcountries(linewidth=2)
  13. m.drawstates(color='b')

接下来，我们可以绘制坐标，从获得它们的实际坐标开始。 记住，南纬和西经坐标需要转换为负值。 例如，纽约市是北纬40.7127西经74.0059。 我们可以在我们的程序中定义这些坐标，如：

  1. NYClat, NYClon = 40.7127, -74.0059

之后我们将这些转换为要绘制的x和y坐标。

  1. xpt, ypt = m(NYClon, NYClat)

注意这里，我们现在已经将坐标顺序翻转为lon, lat（纬度，经度）。 坐标通常以lat, lon顺序给出。 然而，在图形中，lat, long转换为y, x，我们显然不需要。 在某些时候，你必须翻转它们。 不要忘记这部分！

最后，我们可以绘制如下的坐标：

  1. m.plot(xpt, ypt, 'c*', markersize=15)

这个图表上有一个青色的星，大小为 15。更多标记类型请参阅：Matplotlib 标记文档。

接下来，让我们再画一个位置，洛杉矶，加利福尼亚：

  1. LAlat, LAlon = 34.05, -118.25
  2. xpt, ypt = m(LAlon, LAlat)
  3. m.plot(xpt, ypt, 'g^', markersize=15)

这次我们画出一个绿色三角，执行代码会生成：

二十八、在 Basemap 中绘制坐标 - 图1

如果我们想连接这些图块怎么办？原来，我们可以像其它 Matplotlib 图表那样实现它。

首先，我们将那些xpt和ypt坐标保存到列表，类似这样的东西：

  1. xs = []
  2. ys = []
  4. NYClat, NYClon = 40.7127, -74.0059
  5. xpt, ypt = m(NYClon, NYClat)
  6. xs.append(xpt)
  7. ys.append(ypt)
  8. m.plot(xpt, ypt, 'c*', markersize=15)
  10. LAlat, LAlon = 34.05, -118.25
  11. xpt, ypt = m(LAlon, LAlat)
  12. xs.append(xpt)
  13. ys.append(ypt)
  14. m.plot(xpt, ypt, 'g^', markersize=15)
  16. m.plot(xs, ys, color='r', linewidth=3, label='Flight 98')

会生成：

二十八、在 Basemap 中绘制坐标 - 图2

太棒了。有时我们需要以圆弧连接图上的两个坐标。如何实现呢？

  1. m.drawgreatcircle(NYClon, NYClat, LAlon, LAlat, color='c', linewidth=3, label='Arc')

我们的完整代码为：

  1. from mpl_toolkits.basemap import Basemap
  2. import matplotlib.pyplot as plt
  4. m = Basemap(projection='mill',
  5.             llcrnrlat = 25,
  6.             llcrnrlon = -130,
  7.             urcrnrlat = 50,
  8.             urcrnrlon = -60,
  9.             resolution='l')
  11. m.drawcoastlines()
  12. m.drawcountries(linewidth=2)
  13. m.drawstates(color='b')
  14. #m.drawcounties(color='darkred')
  15. #m.fillcontinents()
  16. #m.etopo()
  17. #m.bluemarble()
  19. xs = []
  20. ys = []
  22. NYClat, NYClon = 40.7127, -74.0059
  23. xpt, ypt = m(NYClon, NYClat)
  24. xs.append(xpt)
  25. ys.append(ypt)
  26. m.plot(xpt, ypt, 'c*', markersize=15)
  28. LAlat, LAlon = 34.05, -118.25
  29. xpt, ypt = m(LAlon, LAlat)
  30. xs.append(xpt)
  31. ys.append(ypt)
  32. m.plot(xpt, ypt, 'g^', markersize=15)
  34. m.plot(xs, ys, color='r', linewidth=3, label='Flight 98')
  35. m.drawgreatcircle(NYClon, NYClat, LAlon, LAlat, color='c', linewidth=3, label='Arc')
  38. plt.legend(loc=4)
  39. plt.title('Basemap Tutorial')
  40. plt.show()

结果为：

二十八、在 Basemap 中绘制坐标 - 图3

这就是Basemap的全部了，下一章关于 Matplotlib 的 3D 绘图。