TRMM 3B43是衛星觀測月平均降水量產品，是HDF的格點數據。需要注意的是數據中降水變量維的順序里經度維在前緯度維在后，這與通常的設置（緯度維在前經度維在后）相反，需要對獲取的二維數組進行轉置，使用transpose函數。數據的經緯度數組也需要定義（文件中沒有）。數據中的0值為缺測值（fill_value），并不是一個好的缺測值設置，因此改為-9999.0。

腳本程序：

#Add data file folder = 'D:/Temp/hdf/' fns = '3B43.100301.6A.HDF' fn = folder + fns f = addfile(fn) vname = 'precipitation' t = f[vname] rain = t[0,:,:] rain = transpose(rain) rain[rain==0] = -9999.0 rain.fill_value = -9999.0 lat = arange(-49.875, 49.875, 0.249375) lon = arange(-179.875, 179.876, 0.25) #Plot axesm() mlayer = shaperead('D:/Temp/map/country1.shp') geoshow(mlayer) layer = imshowm(lon, lat, rain, 20) colorbar(layer, orientation='horizontal') title(vname + 'at scan=0 (mm/hr)') axism()

轉載于:https://www.cnblogs.com/yaqiang/p/4629276.html

總結

以上是生活随笔為你收集整理的MeteoInfoLab脚本示例：TRMM 3B43 HDF数据的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。