TRMM 3B43是衛(wèi)星觀測月平均降水量產(chǎn)品，是HDF的格點(diǎn)數(shù)據(jù)。需要注意的是數(shù)據(jù)中降水變量維的順序里經(jīng)度維在前緯度維在后，這與通常的設(shè)置（緯度維在前經(jīng)度維在后）相反，需要對獲取的二維數(shù)組進(jìn)行轉(zhuǎn)置，使用transpose函數(shù)。數(shù)據(jù)的經(jīng)緯度數(shù)組也需要定義（文件中沒有）。數(shù)據(jù)中的0值為缺測值（fill_value），并不是一個(gè)好的缺測值設(shè)置，因此改為-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()

轉(zhuǎn)載于:https://www.cnblogs.com/yaqiang/p/4629276.html

總結(jié)

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

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。