如何绘制在基础地图投影的散点图？

Question

下面是我的代码试图在一些其他数据的基础地图投影的顶部绘制散点。我在我的代码中的实际路径数据。

该散点不会显示在基础地图投影的顶部。

from mpl_toolkits.basemap import Basemap, shiftgrid, cm
import numpy as np
import matplotlib.pyplot as plt
from netCDF4 import Dataset
import matplotlib as mpl
samosdata = Dataset()
d = Dataset()


sst = d.variables["analysed_sst"][:][0, ::-1, :]  
lon = d.variables["lon"][:]
lat = d.variables["lat"][::-1]
lats = samosdata.variables['lat']
lons = samosdata.variables['lon']
time = samosdata.variables['time']
ts = samosdata.variables['TS']
ts = np.array(ts)

lons = np.array(lons)
lats = np.array(lats)
fig = plt.figure()
ax = fig.add_axes([0.1, 0.1, 0.8, 0.8])

m = Basemap(projection='mill', llcrnrlat=-80, urcrnrlat=80,
        llcrnrlon=-180, urcrnrlon=180, lat_ts=20, resolution='c')

nx = int((m.xmax - m.xmin)/11113.2); ny = int((m.ymax m.ymin)/11113.2)
sst = m.transform_scalar(sst, lon, lat, nx, ny)
im = m.imshow(sst, interpolation = "none")


x, y =m(list(lons), list(lats)) 
pts = plt.scatter(x,y, marker='o', s=5, c=ts, cmap='plasma')
plt.colorbar(pts)
m.drawcoastlines()
parallels = np.arange(-90, 90, 30)
meridians = np.arange(-180, 180, 60)
m.drawparallels(parallels, labels = [1, 0, 0, 1])
m.drawmeridians(meridians, labels = [1, 0, 0, 1])

cb = m.colorbar(im, "right", size = "5%", pad = "2%")
ax.set_title("SST 2010 01 01")
plt.show()

当我输入此为界限，将绘制散点，但不会accuratly显示背景数据。

m = Basemap(projection='mill', llcrnrlat=min(lats)-30, 
urcrnrlat=max(lats)+30,
        llcrnrlon=min(lons)-30, urcrnrlon=max(lons)+30, lat_ts=20, 
resolution='c')

Answer 1

更换功能plt.scatter()与m.scatter()。和使用选项latlon=True应该解决您的问题

...
pts = m.scatter(x,y, marker='o', s=5, c=ts, cmap='plasma', latlon=True)
...

希望这是有用的。