0

私は 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 つのデータ セットがある日がある (この場合、4 つのデータセットがある日がある) ため、hstack または vstack を使用してこれらのデータセットをマージできないことです。

私の意図は、3 つの異なるデータ配列から 1 つの配列を取得することです。

このリンクとともにデータファイルも添付しました: https://drive.google.com/open?id=0B2rkXkOkG7ExYW9RNERaZU5lam8

あなたの助けは非常に高く評価されます。

4

0 に答える 0