MeteoInfoLab脚本示例：SMAP L4 HDF数据

MeteoInfo

SMAP (Soil Moisture Active Passive) Level 4 HDF数据的处理示例。数据并非等经纬度投影，从数据中读出经纬度二维数组（cell_lon, cell_lat），注意纬度数组要进行反序（从由北向南改为由南向北），相应读取某要素数据（这里示例sruface_temp）时也要进行纬度方向反序处理。然后将数据再次格点化为等经纬度数据，用griddata()函数，注意插值方法要选为method=surface。

1. #Add data file
   
2. print 'Read data...'
   
3. fn = 'D:/Temp/hdf/SMAP_L4_SM_gph_20150414T013000_Vb1010_001.h5'
   
4. f = addfile(fn)
   
5. lon = f['cell_lon'][:,:]
   
6. lat = f['cell_lat'][::-1,:]
   
7. vname = 'surface_temp'
   
8. v = f[vname]
   
9. data = v[::-1,:]
   
10. #Interpolate data to grid
    
11. print 'Interpolate data to grid...'
    
12. lon1 = linspace(lon.min(), lon.max(), lon.dimlen(1))
    
13. lat1 = linspace(lat.min(), lat.max(), lat.dimlen(0))
    
14. data1 = griddata((lon, lat), data, xi=(lon1, lat1), method='surface')[0]
    
15. #Plot
    
16. print 'Plot...'
    
17. axesm()
    
18. mlayer = shaperead('D:/Temp/map/country1.shp')
    
19. geoshow(mlayer)
    
20. layer = imshowm(lon1, lat1, data1, 20)
    
21. colorbar(layer)
    
22. title('{0}\n {1}'.format(fn, vname))
    
23. print 'Finish!'

#Add data file<br /> print 'Read data...'<br /> fn = 'D:/Temp/hdf/SMAP_L4_SM_gph_20150414T013000_Vb1010_001.h5'<br /> f = addfile(fn)<br /> lon = f['cell_lon'][:,:]<br /> lat = f['cell_lat'][::-1,:]<br /> vname = 'surface_temp'<br /> v = f[vname]<br /> data = v[::-1,:]<br /> #Interpolate data to grid<br /> print 'Interpolate data to grid...'<br /> lon1 = linspace(lon.min(), lon.max(), lon.dimlen(1))<br /> lat1 = linspace(lat.min(), lat.max(), lat.dimlen(0))<br /> data1 = griddata((lon, lat), data, xi=(lon1, lat1), method='surface')[0]<br /> #Plot<br /> print 'Plot...'<br /> axesm()<br /> mlayer = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(mlayer)<br /> layer = imshowm(lon1, lat1, data1, 20)<br /> colorbar(layer)<br /> title('{0}\n {1}'.format(fn, vname))<br /> print 'Finish!'

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数据文件试试，如果这个文件能行可能是不同时期数据文件差异导致的。

MeteoInfo

YF2015 发表于 2016-5-27 12:56
土壤湿度数据要怎么改呢，我改了一下还是不行，提示行列数出错，请王老师帮忙了

vname = 'surface_temp'
改为：
vname = 'sm_surface_wetness'

MeteoInfo

YF2015 发表于 2016-5-26 19:42
刚才试了一下读取温度数据，很卡，然后就出来这样的提示：是什么问题
run script...
Read data...

重新下载MeteoInfo Java 1.3.2R13版试试

YF2015

王老师 太感谢

YF2015

如果做广西地区和中国的，用哪个.shp文件，我的.shp不能再软件里用

MeteoInfo

YF2015 发表于 2016-5-26 16:33
如果做广西地区和中国的，用哪个.shp文件，我的.shp不能再软件里用

MeteoInfo软件里有一些shape文件（在map目录中），你自己打开看看。

YF2015

王老师，要是读取里面的土壤湿度数据要改动哪个地方比如是这个“土壤湿度”，

YF2015

刚才试了一下读取温度数据，很卡，然后就出来这样的提示：是什么问题
run script...
Read data...
Interpolate data to grid...
Traceback (most recent call last):
  File "<iostream>", line 14, in <module>
  File "C:\MeteoInfo\pylib\mipylib\minum.py", line 1447, in griddata
    r = ArrayUtil.interpolation_Surface(x_s.asarray(), y_s.asarray(), values, x_g.asarray(), y_g.asarray())
at java.util.Arrays.copyOf(Unknown Source)
at java.util.ArrayList.grow(Unknown Source)
at java.util.ArrayList.ensureExplicitCapacity(Unknown Source)
at java.util.ArrayList.ensureCapacityInternal(Unknown Source)
at java.util.ArrayList.add(Unknown Source)
at org.meteoinfo.shape.PolygonShape.updatePolygons(PolygonShape.java:278)
at org.meteoinfo.shape.PolygonShape.setPoints(PolygonShape.java:188)
at org.meteoinfo.data.ArrayUtil.interpolation_Surface(ArrayUtil.java:1521)
at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)
at sun.reflect.NativeMethodAccessorImpl.invoke(Unknown Source)
at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)
at java.lang.reflect.Method.invoke(Unknown Source)
java.lang.OutOfMemoryError: java.lang.OutOfMemoryError: Java heap space

YF2015

MeteoInfo 发表于 2016-5-26 21:43
重新下载MeteoInfo Java 1.3.2R13版试试

土壤湿度数据要怎么改呢，我改了一下还是不行，提示行列数出错，请王老师帮忙了

YF2015

MeteoInfo 发表于 2016-5-27 13:34
vname = 'surface_temp'
改为：
vname = 'sm_surface_wetness'

run script...
Read data...
Traceback (most recent call last):
  File "<iostream>", line 8, in <module>
TypeError: str indices must be integers
出现这样的提示