Pythonで衛星データの.nc画像を読み込んで表示する

from osgeo import gdal
import numpy as np
from netCDF4 import Dataset
import matplotlib.pyplot as plt

data = Dataset('image/AT_316_S2_10m_256.nc')
print (data)
#B2バンドを抽出
B2 = data['B2']
#Timeスライス(0~5の6枚から選ぶ)
im = B2[3,:,:]
#与えられた.ncが上下反転していたので戻す
im = np.flipud(im)
plt.figure(figsize=(6,6))
plt.imshow(im)
plt.axis('off')