2020-02-21

• 反射率数据：MOD09
• 温度数据：MOD11_L2

def open_gdal(fileName):
  hdf = gdal.Open(fileName)
  # subDataset
  subdatasets = hdf.GetSubDatasets()
  for subDataset in subdatasets:
      print(subDataset)
  # size
  print(hdf.RasterXSize)
  print(hdf.RasterYSize)
  # SDSs
  sdsdict = hdf.GetMetadata('SUBDATASETS')
  sdslist = [sdsdict[k] for k in sdsdict.keys() if '_NAME' in k]
  sds = []
  for n in sdslist:
      sds.append(gdal.Open(n))
  print(len(sds))

  # Dataset
  data = sds[0].ReadAsArray()