Блог

У меня есть сетка точек данных со значением широты / долготы для каждой строки / столбца. Я хочу построить это, чтобы получить красивую карту. Как мне это сделать? Я попытался преобразовать заголовки строк / столбцов широты / долготы в декартовы, но не могу заставить это работать.

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

general_file = 'datair2_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 /

Вы могли бы сделать это с помощью моего пакета ncplot (https://pypi.org/project/ncplot /).

 import ncplot
ncplot.view('datair2_20160815_114117_232_l3b_v10.nc', "radiance")