使用GDAL模块打开HDF文件，并显示其图像

# -*- coding: utf-8 -*-
from osgeo import gdal
import pylab
from get_filename import file_name
 
ds=gdal.Open("h:/Mark/Mark_1/data/VIPPHEN_EVI2.A1981.004.2016180133348.hdf")
#subdatasets = ds.GetSubDatasets()
#print(subdatasets)


path = 'HDF4_EOS:EOS_GRID:"h:/Mark/Mark_1/data/'
path += 'VIPPHEN_EVI2.A1981.004.2016180133348.hdf'
path += '":VIP_CMG_GRID:"Start of Season 1"'

#start_of_season_1 = 

ds = gdal.Open(path)
band = ds.GetRasterBand(1)
data = band.ReadAsArray()
nodata_value = band.GetNoDataValue() 

#pylab.imshow()
pylab.imshow(data)
pylab.show()