我需要使用从 txt 文件中读取的数据来绘制地图(纬度、经度)。
数据写成:
lat lon var
32 12 55.65
(..)
每个向量(lat
、lon
和 var)
)的值个数为 16384。
所以我认为我需要将 var
(16384) 转换为 var_matrix
(128,128)。
为此,我尝试了以下代码:
em_m=[]
n=0
for i in range (ncol):
for j in range (nrow):
var_m[i,j]=var[n]
n=n+1
但 var_m
始终等于 0。
print(var[0])
1.12
var_m[0,0]=var[0]
print(var_m[0,0])
0
var_m
的最佳方式是什么?mp = Basemap(projection='merc', ,
llcrnrlat=lat_south, #15, #lower left corner lat
llcrnrlon=lon_west, #-130, #lower left conrner lon
urcrnrlat=lat_north, #75, #uper righ corner lat
urcrnrlon=lon_east, #50, #uper right corner lon
resolution='c')
lons,lats = np.meshgrid(lon,lat)
x,y = mp(lons,lats)
mp.drawcoastlines(linewidth=1.)
c_scheme = mp.imshow (np.squeeze(var_m[:,:]),cmap=gist_nacr, interpolation = 'lanczos')