MeteoInfoLab脚本示例：MODIS Sinusoidal投影HDF数据

MeteoInfo

MODIS卫星很多陆面数据都是Sinusoidal投影，数据被分为一个个10*10度（赤道地区）的瓦片（http://modis-land.gsfc.nasa.gov/MODLAND_grid.html），数据格式通常是HDF-EOS格式，作为EOS格式MeteoInfo可以自动读出数据的坐标信息（包括投影等）。

脚本程序：

1. f = addfile('D:/Temp/Hdf/MOD10A1.A2010365.h27v05.005.2011002103013.hdf')
   
2. vname = 'snow_Albedo_Daily_Tile'
   
3. data = f[vname][:,:]
   
4. axesm(projinfo=f.proj, griddx=5, griddy=5, gridline=True)
   
5. mlayer = shaperead('D:/Temp/map/country1.shp')
   
6. geoshow(mlayer)
   
7. levs = arange(0, 100, 5)
   
8. layer = imshowm(data, 20, proj=f.proj)
   
9. title(vname)
   
10. colorbar(layer)
    
11. axism()

f = addfile('D:/Temp/Hdf/MOD10A1.A2010365.h27v05.005.2011002103013.hdf') <br /> vname = 'snow_Albedo_Daily_Tile' <br /> data = f[vname][:,:] <br /> axesm(projinfo=f.proj, griddx=5, griddy=5, gridline=True) <br /> mlayer = shaperead('D:/Temp/map/country1.shp') <br /> geoshow(mlayer) <br /> levs = arange(0, 100, 5) <br /> layer = imshowm(data, 20, proj=f.proj) <br /> title(vname) <br /> colorbar(layer) <br /> axism()

放大到数据范围：


在创建axesm坐标系是不给出投影参数，缺省是经纬度投影，数据图形会自动投影。

1. axesm(griddx=30, griddy=30, gridline=True)

axesm(griddx=30, griddy=30, gridline=True)

放大到数据范围：

spi

{:5_236:}

837078493

楼主你好 我想问一下  这个投影出来的数据 就是正常的经纬度对应的500m*500m的图像了吗？这个 Sinusoidal投影 那么和WRF所用的兰波托投影 岂不是对应不上？因为我想 和WRF结合起来  这个是1200km*1200km的  我的WRF是44Km*44Km的  谢谢楼主了

MeteoInfo

837078493 发表于 2015-10-24 09:55
楼主你好 我想问一下  这个投影出来的数据 就是正常的经纬度对应的500m*500m的图像了吗？这个 Sinusoidal投 ...

这两天很忙（正在单位i加班），等有空了可以帮你处理一下。

837078493

MeteoInfo 发表于 2015-10-24 12:38
这两天很忙（正在单位i加班），等有空了可以帮你处理一下。

哈哈  谢谢你啊   

rceclx

MeteoInfo 发表于 2015-10-24 12:38
这两天很忙（正在单位i加班），等有空了可以帮你处理一下。

王老师，这么辛苦啊！

MeteoInfo

837078493 发表于 2015-10-24 09:55
楼主你好 我想问一下  这个投影出来的数据 就是正常的经纬度对应的500m*500m的图像了吗？这个 Sinusoidal投 ...

你的WRF domain的投影参数，中心经纬度（或者左下角），X/Y方向网格数等是什么？最好提供一下Sinuoidal投影的数据文件（或者下载地址），我抽空试试。其实就是投影转换的问题。

837078493

MeteoInfo 发表于 2015-10-25 09:23
你的WRF domain的投影参数，中心经纬度（或者左下角），X/Y方向网格数等是什么？最好提供一下Sinuoidal投 ...

您好  我的WRF namelist 设置：parent_id         =   1,   1,   2
parent_grid_ratio =   1,   5,   5
i_parent_start    =   1,  12,  48
j_parent_start    =   1,  12,  48
e_we              =  45, 111,  81
e_sn              =  45, 111,  81
geog_data_res     = '10m','2m','30s'
dx = 25000,
dy = 25000,
map_proj = 'lambert',
ref_lat   =  35.946,
ref_lon   =  104.137,
truelat1  =  30.0,
truelat2  =  60.0,
stand_lon =  104.137,
MODIS 数据太大 穿不上来  链接是：http://e4ftl01.cr.usgs.gov/MOLT/MOD09A1.006/2008.08.04/
这个链接里面的MOD09A1.A2008217.h26v05.006.2015177090630.hdf  数据  我想在MODIS数据里面提取WRF的最内层区域的数据  谢谢你啊

MeteoInfo

837078493 发表于 2015-10-25 09:51
您好  我的WRF namelist 设置：parent_id         =   1,   1,   2
parent_grid_ratio =   1,   5,   5 ...

下载起来太慢了，建议你放在百度云之类的地方方便下载。

837078493

MeteoInfo 发表于 2015-10-26 17:13
下载起来太慢了，建议你放在百度云之类的地方方便下载。

http://pan.baidu.com/s/1o6irQJ4