有趣的问题。我试着基本上遵循@Skeletor概述的想法，并Map图像，以便它可以用plot_surface显示：

import PIL
import matplotlib.pyplot as plt
import numpy as np
from mpl_toolkits.mplot3d import Axes3D

# load bluemarble with PIL
bm = PIL.Image.open('bluemarble.jpg')
# it's big, so I'll rescale it, convert to array, and divide by 256 to get RGB values that matplotlib accept 
bm = np.array(bm.resize([d/5 for d in bm.size]))/256.

# coordinates of the image - don't know if this is entirely accurate, but probably close
lons = np.linspace(-180, 180, bm.shape[1]) * np.pi/180 
lats = np.linspace(-90, 90, bm.shape[0])[::-1] * np.pi/180 

# repeat code from one of the examples linked to in the question, except for specifying facecolors:
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')

x = np.outer(np.cos(lons), np.cos(lats)).T
y = np.outer(np.sin(lons), np.cos(lats)).T
z = np.outer(np.ones(np.size(lons)), np.sin(lats)).T
ax.plot_surface(x, y, z, rstride=4, cstride=4, facecolors = bm)

plt.show()

结果：x1c 0d1x

这是我几个小时前做的：
首先，我们导入所需的库：

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
import imageio

其次，我们制作这些图，并将它们作为png存储在我们的目录中：注意我写了range(0,330,20)

for i in range(0,330,20):
    my_map = Basemap(projection='ortho', lat_0=0, lon_0=i, resolution='l', area_thresh=1000.0)
    my_map.bluemarble()
    my_map.etopo()
    name=str(i)
    path='/path/to/your/directory/'+name
    plt.savefig(path+'.png')
    plt.show()
    plt.clf()
    plt.cla()
    plt.close()

最后，我们可以将所有图像加入一个动画GIF：

images = []
for f in range(0,330,20):
    images.append(imageio.imread("/path/to/your/directory/"+str(f)+".png"))
    imageio.mimsave('movie.gif', images, duration=0.5)

然后享受结果：

我可以想象以下解决方案：使用numpy.roll，你可以在每次调用时将数组移动一列（或更多）。因此，您可以将地球表面的图像加载到numpy数组中作为模板，并将旋转后的图像导出为jpg。您将按照扭曲示例中所示进行绘制。