Oleksandr Huziy

Here is a small workaround for you: ``` python from mpl_toolkits import basemap from matplotlib.lines import Line2D import matplotlib.pyplot as plt m1 = basemap.Basemap(projection='cea', lon_0=0) m1.drawcoastlines() lon = [-135, -45,...

Actually this also works fine: ``` python from mpl_toolkits import basemap m1 = basemap.Basemap(projection='cea', lon_0=0) m1.drawcoastlines() lon = [-135, -45, 45, 135] lat = [45, 45, 45, 45] m1.plot(*m1(lon[:2], lat[:2]))...

I think a good fix would be to check the latlon=True and the length of the arrays and throw an exception if the length is less than 3, with an...

Actually it is also calling shiftdata and the shiftdata is kind of limited, here is the doc to shiftdata: ``` def shiftdata(self,lonsin,datain=None,lon_0=None): """ Shift longitudes (and optionally data) so that...

Duplicating the code here, I could not make it highlighted above, probably because I sent it using email... ``` python def shiftdata(self,lonsin,datain=None,lon_0=None): """ Shift longitudes (and optionally data) so that...

I am pretty sure that shiftdata is not needed for scatter, so I've created a PR for that, but I am not sure about plot though...

Hi @neillewis96: you have to convert (lon, lat) to the projection coordinates ```python xx, yy = map(lon, lat) map.pcolormesh(xx,yy,convrain_mean) ``` Cheers

Hi @neillewis96: I've tried to reproduce your issue with random data, but it seems to be OK for me... Which version of basemap are you using. Here is the notebook...

Hi @abc123634: I think I might have a solution for you, because readshapefile actually saves the info into fields of the basemap object, so you could reuse it without calling...

Hi Daryl, It works like this, pcolormesh fills space between the point defined by X,Y, so this way you get the number of intervals less than the number of points......