我在绘制多边形时遇到了困难，例如，在使用matplotlib basemap生成的地图边界上绘制多边形。在下面的示例中，地图边界由日期线指定。我试图通过指定三角形顶点的坐标来绘制一个跨越日期线的三角形。当所有的坐标都在地图内时，这种方法很有效，但是如果它们穿过地图边界，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]