Python: 如何在带有散点图的地图上绘制数组？

Question

让我介绍一下背景。我有一个二进制文件，我必须在其中提取数据。这是一个巨大的矩阵，每个像素的深度为longitude/latitude。代码的第一部分是一小段代码，用于读取和存储每个数据（在预定义的 lon/lat 间隔中）。但是，我想在世界地图上绘制这些数据，但我被阻止了。我可以在带有网格的空白图形上绘制数据，但我不知道如何插入世界地图（使用底图）。

# -*- coding: utf-8 -*-
import numpy as np
from matplotlib import pyplot as plt
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap
from matplotlib.colors import LinearSegmentedColormap

# Definition of the grid
GLDB_res = 1./120.
GLDB_lon = np.linspace(-180.+0.5*GLDB_res, 180.-0.5*GLDB_res, 43200)
GLDB_lat = np.linspace(  90.-0.5*GLDB_res, -90.+0.5*GLDB_res, 21600)

j_Asia = pl.find(pl.logical_and(GLDB_lon>=0.,GLDB_lon<=40.))
i_Asia= pl.find(pl.logical_and(GLDB_lat>=50.,GLDB_lat<=90.))

# Read depth

GLDB_depth_file = 'fileURL/GlobalLakeDepth.dat'
GLDB_depth = np.fromfile(GLDB_depth_file,'<i2')
GLDB_depth_r = GLDB_depth.reshape((21600,43200))

GLDB_depth_Asia = GLDB_depth_r[i_Asia[0]:i_Asia[-1]+1,j_Asia[0]:j_Asia[-1]+1].astype(float)
GLDB_depth_Asia[np.where(GLDB_depth_Asia==0.)] = np.nan

# Création de la Basemap

map=Basemap(projection='cyl', llcrnrlat=50,urcrnrlat=90,\
        llcrnrlon=0,urcrnrlon=180, resolution='l')
map.drawcoastlines()
map.fillcontinents(color='#cc9955')

#Colormap
cmap_depth  = LinearSegmentedColormap.from_list("my_colormap", ((200./255,200/255.,1.), (20./255,20./255,165./255)), N=255, gamma=1.0)

#Graphique

plt.close()
fig=plt.figure(1,facecolor='w',figsize=(11,7))
ax=plt.subplot(1,2,1)
x,y=map(j_Asia,i_Asia) #associe un x et un y au lon et lat 
plt.scatter(x,y, c='my_colormap') #trace les données lacs avec Basemap

c=plt.colorbar(orientation='horizontal')
for i in np.arange(0,len(i_Asia),120/2): #res:1/120th deg
map=plt.plot(np.array([0,len(i_Asia)]),np.array([i,i]),'-',c=(150./255,150./255,150./255))
map=plt.plot(np.array([i,i]),np.array([0,len(i_Asia)]),'-',c=(150./255,150./255,150./255))

plt.show()

总而言之，我想插入我的数组 Global_Depth_Asia，它表示区域地图上的所有深度。问题是我有一个无效的 rgb 参数的问题！

谢谢

Answer 1

我显然不能运行测试代码，但是里面有些奇怪的东西。

plt.close() 关闭底图。你为什么要那样做？
plt.figure 创建一个没有底图的新图。这可能是不希望的。如果要使用底图绘图，也请删除 plt.subplot。
使用map.scatter()将数据绘制到地图上。
为 scatter 的 c 参数提供选择颜色所依据的数据

为 cmap 参数提供颜色图。

m.scatter(x,y, c=data, cmap=cmap_depth)

Python: 如何在带有散点图的地图上绘制数组？

Python: How can I draw an array over a map with scatter?

python

arrays

matplotlib

ipython

matplotlib-basemap