我在繪制多邊形時遇到了困難，例如，在使用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]

創作挑戰賽新人創作獎勵來咯，堅持創作打卡瓜分現金大獎

總結

以上是生活随笔為你收集整理的python plot map_使用matplotlibbasemap在边界打印的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。