如何映射netcdf数据ob基本映射

时间:2019-04-17 12:36:44

标签: geospatial google-colaboratory matplotlib-basemap netcdf4

文件包含回声值(经度/纬度)到经纬度,我必须在底图上绘制完整范围的回声。

from netCDF4 import Dataset
import numpy as np
import pandas as pd
from google.colab import files
upload = files.upload()
my_example_nc_file = 'a.nc'
fh = Dataset(my_example_nc_file, mode='r')
lons = fh.variables['longitude'][:]
lats = fh.variables['latitude'][:]
ech= fh.variables['echos'][:]

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
%matplotlib inline
m = Basemap(width=5000000,height=3500000,
        resolution='l',projection='stere',\
        lat_ts=40,lat_0=lat_0,lon_0=lon_0)

xi, yi = m(lons, lats)
#simple plot
#m.plot(xi, yi, 'co')
m.scatter(rge,yi, marker = 'o', color='r', zorder=5)

当前代码执行结果如下。

enter image description here 我要绘制总回波,其颜色由变化表示,如以下屏幕简短显示

enter image description here

0 个答案:

没有答案