如何结合numpy ndarray?

时间:2016-08-01 14:37:27

标签: python numpy multidimensional-array pyhdf

我有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)

enter image description here

问题是,在某些日子里,有3个数据集(在这种情况下,有些日子有4个数据集)所以我不能使用hstack或vstack来合并这些数据集。

我的目的是从三个不同的数据阵列中获取单个数组。

我还附上了数据文件以及此链接:https://drive.google.com/open?id=0B2rkXkOkG7ExYW9RNERaZU5lam8

您的帮助将受到高度赞赏。

0 个答案:

没有答案