散点图不上绘制底图

Question

任何点我有一个数组anomalies_ind这是以这种方式创建：

data_path = r"C:\Users\matth\Downloads\TRMM_3B42RT\3B42RT_Daily.201001.7.nc4" 
f = Dataset(data_path) 

latbounds = [ -45 , -10 ] 
lonbounds = [ 105, 160 ] 
lats = f.variables['lat'][:] 
lons = f.variables['lon'][:] 

# latitude lower and upper index 
latli = np.argmin(np.abs(lats - latbounds[0])) 
latui = np.argmin(np.abs(lats - latbounds[1])) 

# longitude lower and upper index 
lonli = np.argmin(np.abs(lons - lonbounds[0])) 
lonui = np.argmin(np.abs(lons - lonbounds[1])) 

precip_subset = f.variables['precipitation'][ : , lonli:lonui , latli:latui ] 

data_low_indices1 = np.where((precip_subset > 0) & (precip_subset < 1)) 
data_low_indices2 = np.array(np.where((precip_subset > 0) & (precip_subset < 1))).T 
anomalies_ind = [] 
for ind in data_low_indices2: 
    anomalies_ind.append(ind) 
    print(np.asarray(anomalies_ind)) 

这样做的输出是这样的：

[[1, 23, 45] 
[3, 45, 56] 
... 
[31, 45, 89]] 

第一元素表示天在1月份，第二和第三元素分别代表经度和纬度。我想在地图上给出像这样的经度和纬度积点：

foo = np.asarray(anomalies_ind) 
longs = foo[:,1] 
lat = foo[:,2] 
m = Basemap(llcrnrlon=105.,llcrnrlat=-45,urcrnrlon=160,urcrnrlat=-10) 
m.drawcoastlines() 
m.fillcontinents(color = 'lightgray', zorder = 0) 
m.scatter(longs, lat, marker = 'o', color = 'k', zorder=10) 
plt.show() 

然而，没有点在地图上。有谁知道什么是错的？

编辑：这里是真正的foo阵列”

[[ 0 0 0] 
[ 0 0 16] 
[ 0 0 17] 
..., 
[ 30 219 113] 
[ 30 219 114] 
[ 30 219 116]] 

Answer 1

我碰到这个问题，我自己

底图绘图功能有一个latlon关键字，并改变这对我有用。默认值是latlon=False，以便将x和y值解释为投影坐标。添加latlon=True告诉底图将x和y值解释为地图坐标。

看到这里上分散的底图文档： http://matplotlib.org/basemap/api/basemap_api.html#mpl_toolkits.basemap.Basemap.scatter

在原始的绘图语法，你要改线：

m.scatter(longs, lat, marker = 'o', color = 'k', zorder=10) 

到：

m.scatter(longs, lat, marker = 'o', color = 'k', zorder=10, latlon=True) 

Answer 2

你有你的纬度转换ANS经度调用散射之前地图投影的一些值：

x,y=m(longs,lat) 
m.scatter(x, y, marker = 'o', color = 'k', zorder=10)