我有 MODIS 大气产品。我使用下面的代码来读取数据。
%matplotlib inline
import numpy as np
from pyhdf import SD
import matplotlib.pyplot as plt
files = ['file1.hdf','file2.hdf','file3.hdf']
for n in files:
hdf=SD.SD(n)
lat = (hdf.select('Latitude'))[:]
lon = (hdf.select('Longitude'))[:]
sds=hdf.select('Deep_Blue_Aerosol_Optical_Depth_550_Land')
data=sds.get()
attributes = sds.attributes()
scale_factor = attributes['scale_factor']
data= data*scale_factor
plt.contourf(lon,lat,data)
问题是,在某些日子里,有 3 个数据集(在这种情况下,有些日子有四个数据集),所以我不能使用 hstack 或 vstack 来合并这些数据集。
我的意图是从三个不同的数据数组中获取单个数组。
我还附上了这个链接的数据文件:https ://drive.google.com/open?id=0B2rkXkOkG7ExYW9RNERaZU5lam8
您的帮助将不胜感激。