使用matplotlib-basemap绘制边界

Question

我在绘图方面遇到困难，例如跨使用matplotlib-basemap生成的地图边界的多边形。在下面的示例中，地图边界由日期线指定。我尝试通过指定三角形顶点的坐标在日期线上绘制一个三角形。当所有坐标都在地图内时，这可以正常工作，但如果它们穿过地图边界，则底图执行奇怪的外推，因为它似乎不知道如何以正确的方式绘制矩形。

在我看来，正确的方式意味着绘制三角形直到地图边界，然后继续在地图的另一边。

以下是最小代码示例和说明一般问题的图。 如何以一般方式解决这个问题的任何想法都非常受欢迎。

from mpl_toolkits.basemap import Basemap
import matplotlib.pylab as plt
import numpy as np
import matplotlib.path as mpath
import matplotlib.patches as mpatches
import matplotlib as mpl
from matplotlib.collections import PatchCollection

![plt.close('all')
Path = mpath.Path
fig=plt.figure(); ax=fig.add_subplot(121); ax1=fig.add_subplot(122)

def do_plot(ax,lons,lats,title):

    patches=\[\]
    m = Basemap(projection='robin',  resolution='c',lon_0=0.,ax=ax) #todo: how to make it properly work for other projections ???
    m.drawmapboundary(fill_color='grey')
    m.drawcoastlines()

    #--- generate first sample with no problem
    x,y=m(lons,lats)
    verts = np.asarray(\[x,y\]).T
    codes = \[Path.MOVETO,Path.LINETO,Path.LINETO\]
    patches.append(mpatches.PathPatch(mpath.Path(verts, codes,closed=True)))

    #--- generate collection
    cmap = plt.cm.get_cmap('jet', 50); norm = mpl.colors.Normalize(vmin=None, vmax=None) #colorbar mapping
    collection = PatchCollection(patches, cmap=cmap,norm=norm, alpha=1.,match_original=False)  #construct library of all objects
    colors = np.asarray(np.random.random(len(patches)))
    collection.set_array(np.array(colors)) #assign data values here

    #--- do actual plotting
    im=m.ax.add_collection(collection)

    ax.set_title(title)

do_plot(ax,\[-10.,0.,20.\],\[30.,50.,20.\],'This works')
do_plot(ax1,\[170,180,-175\],\[30.,50.,20.\],'... and here is the boundary problem')

plt.show()][1]

Answer 1

您无法以简单的方式解决Basemap的这个问题。在您的行x,y=m(lons,lats)中，您已将点转换为地图坐标，并且绘制多边形只是在这些投影点之间绘制。

您可以尝试使用Cartopy，这可以执行此操作。 This example可能有帮助。