首页 > 解决方案 > Python:在纬度经度网格中绘制规则的数据网格

问题描述

我有一个数据点网格,每行/列的纬度/经度值。我想绘制这个来给出一个漂亮的地图。我该怎么做?我尝试将纬度/经度行/列“标题”转换为笛卡尔坐标,但无法使其正常工作。

from netCDF4 import Dataset
from mpl_toolkits.basemap import Basemap, cm
import matplotlib.pyplot as plt

general_file = 'data\ir2_20160815_114117_232_l3b_v10.nc'
file_data = Dataset(general_file)

lat = file_data.variables['latitude']
lon = file_data.variables['longitude']
rad = file_data.variables['radiance']

latitudes = lat[:]
longitudes = lon[:]
radiances = rad[0]

x= longitudes
y=latitudes
z = radiances

plt.contourf(x, y, z,25,cmap='Greys')
plt.title(general_file)
plt.colorbar()
plt.show()

数据文件从https://darts.isas.jaxa.jp/pub/pds3/extras/vco_ir2_l3_v1.0/vcoir2_7001/data/l3b/netcdf/r0024/下载

标签: pythonmatplotlibplotcontourf

解决方案


你可以使用我的包 ncplot ( https://pypi.org/project/ncplot/ ) 来做到这一点。

import ncplot
ncplot.view('data\ir2_20160815_114117_232_l3b_v10.nc', "radiance")

推荐阅读