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 단위의 데이터였다.

참고

파이썬에서 NetCDF 읽기 및 쓰기

정보 감사입니다.