我已经尝试并搜索了如何解决此问题,但仍然找不到如何从给定的Modis Aqua .hdf文件在gdal和matplotlib中读取和绘制此图的方法。任何帮助深表感谢。顺便说一句,我在Windows 7中使用Python 2.7.5。文件名是A2014037040000.L2_LAC.SeAHABS.hdf
。在hdf文件的地球物理数据中,我将仅使用chlor_a。
更新:
这是示例文件的链接。
A2014037040500.L2_LAC.SeAHABS.hdf
最佳答案
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)
您还可以打开“主” 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
关于python - 读取/打开modis aqua .hdf文件,并在gdal和matplotlib中显示/绘制输出,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/22703498/