MeteoInfoLab脚本示例：SMAP L4 HDF数据

YF2015

读取数据过程中，名字出错了！

MeteoInfo

YF2015 发表于 2016-5-27 14:09
run script...
Read data...
Traceback (most recent call last):

你把代码贴出来看看

YF2015

#Add data file
print 'Read data...'
fn = 'C:/SMAP/SMAP_L4_SM_gph_20151015T013000_Vb1010_001.h5'
f = addfile(fn)
lon = f['cell_lon'][:,:]
lat = f['cell_lat'][::-1,:]
vname = 'sm_surface'
v = f[vname]
data = v[::-1,:]
#Interpolate data to grid
print 'Interpolate data to grid...'
lon1 = linspace(lon.min(), lon.max(), lon.dimlen(1))
lat1 = linspace(lat.min(), lat.max(), lat.dimlen(0))
data1 = griddata((lon, lat), data, xi=(lon1, lat1), method='surface')[0]
#Plot
print 'Plot...'
axesm()
mlayer = shaperead('C:/MeteoInfo/map/country1.shp')
geoshow(mlayer)
layer = imshowm(lon1, lat1, data1, 20)
colorbar(layer)
title('{0}\n {1}'.format(fn, vname))
print 'Finish!'
出现下面的提示
run script...
Read data...
Interpolate data to grid...
Traceback (most recent call last):
  File "<iostream>", line 14, in <module>
TypeError: 'NoneType' object is unsubscriptable

MeteoInfo

YF2015 发表于 2016-5-27 15:05
#Add data file
print 'Read data...'
fn = 'C:/SMAP/SMAP_L4_SM_gph_20151015T013000_Vb1010_001.h5'

你的脚步程序在我这里运行没有问题，你用SMAP_L4_SM_gph_20150414T013000_Vb1010_001.h5数据文件试试，如果这个文件能行可能是不同时期数据文件差异导致的。

YF2015

没问题了，刚才是因为我.shp的路径出问题了，所以才出错，很不错做出来的图

YF2015

MeteoInfo 发表于 2016-5-27 15:38
你的脚步程序在我这里运行没有问题，你用SMAP_L4_SM_gph_20150414T013000_Vb1010_001.h5数据文件试试，如 ...

王老师，还有一个要问一下，就是把出图的时间加上去，不要文件路径那个名字要怎么设置呢！

MeteoInfo

YF2015 发表于 2016-5-27 16:25
王老师，还有一个要问一下，就是把出图的时间加上去，不要文件路径那个名字要怎么设置呢！

用一个字符串就行了，参考这里：http://www.meteothinker.com/exam ... llite/fy3a_aod.html

YF2015

王老师，问题都解决，请求老师帮助一下，怎么MASK出广西地区出来，看了那个MASK中国的，做不出广西图像！！！

MeteoInfo

YF2015 发表于 2016-5-27 18:57
王老师，问题都解决，请求老师帮助一下，怎么MASK出广西地区出来，看了那个MASK中国的，做不出广西图像！！ ...

在MeteoInfo中利用bou2_4p.shp图层，选中广西，输出为一个只包含广西的shape文件，然后用此文件来做maskout

YF2015

MeteoInfo 发表于 2016-5-29 12:20
在MeteoInfo中利用bou2_4p.shp图层，选中广西，输出为一个只包含广西的shape文件，然后用此文件来做masko ...

别的软件的shp它还不认