MeteoInfoLab脚本示例：AMSR-E卫星数据投影

MeteoInfo

AMSR-E（http://nsidc.org/data/amsre/index.html）数据中的Land3数据是HDF-EOS4格式，投影是Cylindrical_Equal_Area。这里示例读取数据并投影至等经纬度投影。

脚本程序：

1. #Add data file
   
2. fn = 'D:/Temp/hdf/AMSR_E_L3_DailyLand_V06_20091231.hdf'
   
3. f = addfile(fn)
   
4. vname = 'D_Soil_Moisture'
   
5. data = f[vname][::-1,:]
   
6. data[data==0] = -9999
   
7. data.fill_value = -9999
   
8. yn = data.dimlen(0)
   
9. xn = data.dimlen(1)
   
10. #Project data
    
11. toproj = projinfo()    #longlat projection
    
12. x = linspace(0, 360, xn)
    
13. y = linspace(-90, 90, yn)
    
14. lon, lat = meshgrid(x, y)
    
15. data = data.project(lon, lat, toproj=toproj)
    
16. #Plot
    
17. axesm()
    
18. mlayer = shaperead('D:/Temp/map/country1.shp')
    
19. geoshow(mlayer, edgecolor='k')
    
20. levs = arange(0, 201, 20)
    
21. layer = imshowm(x, y, data, levs, fill_value=0)
    
22. colorbar(layer)
    
23. xticks(arange(0, 361, 60))
    
24. yticks(arange(-90, 91, 60))
    
25. title(vname)

#Add data file<br /> fn = 'D:/Temp/hdf/AMSR_E_L3_DailyLand_V06_20091231.hdf'<br /> f = addfile(fn)<br /> vname = 'D_Soil_Moisture'<br /> data = f[vname][::-1,:]<br /> data[data==0] = -9999<br /> data.fill_value = -9999<br /> yn = data.dimlen(0)<br /> xn = data.dimlen(1)<br /> #Project data<br /> toproj = projinfo()    #longlat projection<br /> x = linspace(0, 360, xn)<br /> y = linspace(-90, 90, yn)<br /> lon, lat = meshgrid(x, y)<br /> data = data.project(lon, lat, toproj=toproj)<br /> #Plot<br /> axesm()<br /> mlayer = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(mlayer, edgecolor='k')<br /> levs = arange(0, 201, 20)<br /> layer = imshowm(x, y, data, levs, fill_value=0)<br /> colorbar(layer)<br /> xticks(arange(0, 361, 60))<br /> yticks(arange(-90, 91, 60))<br /> title(vname)

MeteoInfo

YF2015 发表于 2016-6-1 18:59
dat,nitf,dted,这几个格式也不可以吗

可以写成基本的二进制格式或者netCDF

MeteoInfo

YF2015 发表于 2016-5-24 19:00
SMAP的四级产品，希望老师能给个解读的程序，很感谢！链接：http://pan.baidu.com/s/1pK9LjIj 密码：1k8j

参考此贴：MeteoInfoLab脚本示例：SMAP L4 HDF数据
http://bbs.06climate.com/forum.p ... 292&fromuid=106
(出处: 气象家园)

需要下载最新的MeteoInfo软件。

MeteoInfo

YF2015 发表于 2016-5-16 11:53
王老师，怎么从出来的图像裁剪出自己想要的区域，比如说中国或者一个省的区域呢！

参考此贴：MeteoInfoLab脚本示例：Maskout图形
http://bbs.06climate.com/forum.p ... 169&fromuid=106
(出处: 气象家园)

ziqiangbuxi

王老师，您好，再请教您个问题，就是转换好的经纬的影响，怎么能在meteoinfo中保存为tiff格式导出来，或者是img的格式？能用于GIS直接处理的

MeteoInfo

ziqiangbuxi 发表于 2015-12-29 22:01
王老师，您好，再请教您个问题，就是转换好的经纬的影响，怎么能在meteoinfo中保存为tiff格式导出来，或者 ...

可以加一条语句输出数组为Surfer ASCII Grid格式数据：
data.savegrid(x, y, 'D:/Temp/test.grd')

ziqiangbuxi

MeteoInfo 发表于 2015-12-29 22:48
可以加一条语句输出数组为Surfer ASCII Grid格式数据：
data.savegrid(x, y, 'D:/Temp/test.grd')

真的很感谢王老师，这几天查了很多资料，看到hdf-eos上说明，经过您的指导和摸索，在GIS里面也实现了，但是操作较为繁琐。您教的方法是最为方便的，后面认真学习meteoinfo。再次感谢您的耐心指导。

ziqiangbuxi

王老师，再请教您一个问题，在MeteoInfoLab脚本中，我查询了资料，看到可以插值，但不知道能否进行重采样，就是说，原始数据的分辨率是25km，那能不能把分辨率重采样为5km？谢谢您。

MeteoInfo

ziqiangbuxi 发表于 2016-1-1 23:21
王老师，再请教您一个问题，在MeteoInfoLab脚本中，我查询了资料，看到可以插值，但不知道能否进行重采样， ...

下面的语句其实就是格点数据重采样，用的双线性插值。
data = data.project(lon, lat, toproj=toproj)

当然这里同时包含了投影的计算，如果将toproj参数去掉就是典型的重采样了（不进行投影转换）。

aas314

学习了

ziqiangbuxi

MeteoInfo 发表于 2016-1-3 10:27
下面的语句其实就是格点数据重采样，用的双线性插值。
data = data.project(lon, lat, toproj=toproj)
...

好的，谢谢王老师，正在学习中！

YF2015

王老师我的软件打不开Land3数据

MeteoInfo

YF2015 发表于 2016-5-15 20:25
王老师我的软件打不开Land3数据

你用的什么软件？