读取/打开modis aqua .hdf文件并在gdal和matplotlib中显示/绘制输出

Question

我已经尝试并搜索了如何解决这个问题，但仍无法找到如何在给定的Modis Aqua .hdf文件中读取和绘制gdal和matplotlib中的方法。任何帮助深表感谢。顺便说一句，在Windows 7中使用Python 2.7.5。文件名是A2014037040000.L2_LAC.SeAHABS.hdf。在hdf文件的地球物理数据中，我只使用chlor_a。

更新

以下是示例文件的链接。

A2014037040500.L2_LAC.SeAHABS.hdf

Answer 1

HDF的诀窍是大多数时候你需要一个特定的子数据集。如果使用GDAL，则需要打开直接指向该子数据集的HDF：

import gdal
import matplotlib.pyplot as plt

ds = gdal.Open('HDF4_SDS:UNKNOWN:"MOD021KM.A2013048.0750.hdf":6')
data = ds.ReadAsArray()
ds = None

fig, ax = plt.subplots(figsize=(6,6))

ax.imshow(data[0,:,:], cmap=plt.cm.Greys, vmin=1000, vmax=6000)

你也可以打开'main'HDF文件并检查子数据集，然后从那里开始：

# open the main HDF
ds = gdal.Open('MOD021KM.A2013048.0750.hdf')

# get the path for a specific subdataset
subds = [sd for sd, descr in ds.GetSubDatasets() if descr.endswith('EV_250_Aggr1km_RefSB (16-bit unsigned integer)')][0]

# open and read it like normal
dssub = gdal.Open(subds)
data = dssub.ReadAsArray()
dssub = None

ds = None

Answer 2

您应该尝试为MODIS数据集设置数据类型。我猜它是16位无符号

ds = gdal.Open（hdfpath） data = ds.GetRasterBand（N）.ReadAsArray（）。astype（numpy.uint16）

N是您感兴趣的数据的带编号。您可以尝试使用QGIS或ENVI打开它以查看HDF文件的结构。

请记住，乐队从1开始而不是0.第一乐队是1。

希望有所帮助