python - 在 Python 中绘制 3D 球体上的圆柱地图数据

Question

假设我有一张带有行星圆柱图的图像，比如说其中一个：

http://www.johnstonsarchive.net/spaceart/cylmaps.html

我想将它绘制在一个 3D 球体上以恢复地球的原始图形。

有没有办法使用像 matplotlib、mayavi、底图或类似的 Python 包来做到这一点？

Answer 1

更新：这是使用 Cartopy 的新版本，因为底图是 EOL。以下是原答案。

import matplotlib.pyplot as plt
import cartopy.crs as ccrs

img = plt.imread("/tmp/venuscyl4.tif")

plt.figure(figsize=(3, 3))

ax = plt.axes(projection=ccrs.Orthographic(-10, 45))
ax.gridlines(color='black', linestyle='dotted')
ax.imshow(img, origin="upper", extent=(-180, 180, -90, 90),
          transform=ccrs.PlateCarree())  # Important

plt.show()

---

感谢 Raphael Roth 的回答，我终于找到了我想要的东西：底图方法warpimage。

这是一个非常简单的例子。使用这张金星的圆柱图，并基于食谱的简单示例：

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
import numpy as np
# set up orthographic map projection with
# perspective of satellite looking down at 50N, 100W.
# use low resolution coastlines.
# don't plot features that are smaller than 1000 square km.
bmap = Basemap(projection='ortho', lat_0 = 50, lon_0 = -100,
              resolution = 'l', area_thresh = 1000.)
# plot surface
bmap.warpimage(image='venuscyl4.jpg')
# draw the edge of the map projection region (the projection limb)
bmap.drawmapboundary()
# draw lat/lon grid lines every 30 degrees.
bmap.drawmeridians(np.arange(0, 360, 30))
bmap.drawparallels(np.arange(-90, 90, 30))
plt.show()

产生以下输出：

Answer 2

底图工具包非常适合这项任务。

问题似乎是你没有真正的数据，即每个像素的纬度、经度、值。问题还在于，对于给定的图像，您通常不知道用于创建图像的投影，因此您无法进行反向转换来获取原始数据。

如果您只想绘制地球图像，请使用底图工具包中的 bluemarble() 函数：

http://wiki.scipy.org/Cookbook/Matplotlib/Maps