MeteoInfo介绍

yaohw123456

有中文版的吗，表示英文学的不好，看起来很吃力

w443079344

谢谢分享.....

designwuying

初学者，受教了

letme

王老师，请问我用MeteoInfo -Default为什么打不开TRMM3B43数据呢？谢谢

MeteoInfo

letme 发表于 2016-8-3 20:50
王老师，请问我用MeteoInfo -Default为什么打不开TRMM3B43数据呢？谢谢

参考此贴：MeteoInfoLab脚本汇总贴
http://bbs.06climate.com/forum.p ... 151&fromuid=106
(出处: 气象家园)

letme

MeteoInfo 发表于 2016-8-3 22:13
参考此贴：MeteoInfoLab脚本汇总贴
http://bbs.06climate.com/forum.php?mod=viewthread&tid=36151&from ...

使用您介绍的MeteoInfoLab脚本，改了之后显示
run script...
Traceback (most recent call last):
  File "<iostream>", line 4, in <module>
NameError: name 'folder' is not defined

脚本内容：
#Add data file
#folder = 'E:/TRMM'
fns = '3B43.20010801.7A.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)
#Add data file
folder ='E:/TRMM'
fns = '3B43.20010801.7A.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('E:/TRMM/LYFW.shp')
geoshow(mlayer)
layer = imshowm(lon, lat, rain, 20)
colorbar(layer, orientation='horizontal')
title(vname + 'at scan=0')
axism()
lon = arange(-179.875, 179.876, 0.25)
#Plot
axesm()
mlayer = shaperead('E:/TRMM/LYFW.shp')
geoshow(mlayer)
layer = imshowm(lon, lat, rain, 20)
colorbar(layer, orientation='horizontal')
title(vname + 'at scan=0')
axism()


请问怎么解决呢？谢谢您

MeteoInfo

letme 发表于 2016-8-3 22:28
使用您介绍的MeteoInfoLab脚本，改了之后显示
run script...
Traceback (most recent call last):

#folder = 'E:/TRMM'
改为：
folder = 'E:/TRMM'

语句前面有#表示被注释掉了，不会执行。

letme

MeteoInfo 发表于 2016-8-3 22:30
#folder = 'E:/TRMM'
改为：
folder = 'E:/TRMM'

谢谢您。万忙之中回复。真的非常感谢

letme

letme 发表于 2016-8-3 22:36
谢谢您。万忙之中回复。真的非常感谢

王老师，还要请教下

脚本如下：

#Add data file
folder = 'E:/TRMM/'
fns = '3B43.20010801.7A.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('E:/TRMM/LYFW.shp')
geoshow(mlayer)
layer = imshowm(lon, lat, rain, 20)
colorbar(layer, orientation='horizontal')
title(vname + 'at scan=0')
axism()

显示结果是：
run script...
indices must be 2 dimensions!
Traceback (most recent call last):
  File "<iostream>", line 10, in <module>
  File "D:\MeteoInfo_Java_1.3.4R2_Files\MeteoInfo\pylib\mipylib\minum.py", line 1264, in transpose
    r = ArrayMath.transpose(a.asarray(), dim1, dim2)
AttributeError: 'NoneType' object has no attribute 'asarray'

请问如何处理。谢谢

MeteoInfo

letme 发表于 2016-8-3 22:46
王老师，还要请教下

脚本如下：

t可能是2维变量，
rain = t[0,:,:]
改为：
rain = t[:,:]
试试