WRF / WPS > Jupyter Notebook > OI SST(NetCDFファイル)を読込みMatplotlibで画像表示

動作環境

Xeon E5-2620 v4 (8コア) x 2 
32GB RAM 
CentOS Linux release 7.7.1908 (Core)
Python 3.6.8 (default, Aug  7 2019, 17:28:10) 
IPython 7.9.0 -- An enhanced Interactive Python. Type '?' for help.

関連

• 2020-01-15 WRF > NCEP SST Analysisデータ終了後のSSTデータの選択 > OI SSTデータの使用

実装

OI SSTファイル(NetCDF形式)を読込んで画像表示する。

ReadNetCDF_200116.ipynb

import netCDF4
import numpy as np
import matplotlib.pyplot as plt
import matplotlib.cm as cm

INFILE = '20200111120000-NCEI-L4_GHRSST-SSTblend-AVHRR_OI-GLOB-v02.0-fv02.0.nc'
nc = netCDF4.Dataset(INFILE, 'r')
nlat = len(nc.dimensions['lat'])
nlon = len(nc.dimensions['lon'])
lons = nc.variables['lon'][:]
lats = nc.variables['lat'][:]
ssts = nc.variables['analysed_sst'][:]
print(nlat, nlon)
# print(lons)
# print(lats)
# print(ssts)
nc.close()

grid = ssts.reshape((nlat, nlon))
plt.imshow(grid, interpolation='nearest', cmap=cm.gist_rainbow)
plt.colorbar()

結果

上下反対になったが、実際にはlat, lon情報をもとに後続の処理で使う予定。
lats, lonsは0.25刻みのデータであった。

参考

• PythonでNetCDFの読み書き

情報感謝です。