Question

让我介绍一下背景。我有一个二进制文件，我必须提取数据。这是一个巨大的矩阵，每个经度/纬度的像素深度。代码的第一部分是一个小代码，用于读取和存储每个数据（在预定义的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的{{1}}参数
将色彩映射提供给c参数。
```
cmap
```

Python：如何使用散点在地图上绘制数组？

1 个答案: