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')