def plot_map():
lon=grd.lon; lat=grd.lat
lat_ts=((lat[0,0]+lat[-1,0])*0.5)
cbox=[lon[0, 0],lat[0, 0],lon[-1, -1],lat[-1, -1]]
map = Basemap(projection='merc', resolution='l',
lat_ts=lat_ts,
llcrnrlon=cbox[0], llcrnrlat=cbox[1],
urcrnrlon=cbox[2], urcrnrlat=cbox[3])
map.drawmapboundary(fill_color='white')
map.drawcoastlines()
map.fillcontinents(color='#cc9966', lake_color='#99ffff')
lon_tcks = np.arange(-100.,30.,20.)
lat_tcks = np.arange(0.,70.,10.)
map.drawmeridians(lon_tcks, labels=[True, False, False, True])
map.drawparallels(lat_tcks, labels=[False, True, True, False])
# get projected coordinates
x,y=map(lon,lat)
return map, x, y
评论列表
文章目录