- 积分
- 5574
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-8-12
- 最后登录
- 1970-1-1
|

楼主 |
发表于 2015-5-15 11:45:51
|
显示全部楼层
试了下改了路径和变量名,但是有报错,说是没有getName这个属性?
#Add data file
fn = 'F:/paper/fy2f-ctt/FY2F_CTT_MLT_NOM_20140717_1730.hdf'
f = addfile(fn)
#Get data variable
v = f['ctt']
#Set x/y
x = linspace(-5750000.0,5750000.0, 2288)
y = linspace(-5750000.0,5750000.0, 2288)
#Add x/y dimensions
v.adddim('Y', y)
v.adddim('X', x)
#Get data array
ctt = v[:,:]
#Plotax = worldmap()
proj = axesm(proj='geos', origin=(0,104.5,0), h=35785864, gridlabel=False)
mlayer = shaperead('C:/Users/Administrator/AppData/Local/Temp/Rar$EXa0.954/MeteoInfo/map/country1.shp')
geoshow(ax, mlayer)
layer = imshowm(ax, olr, 20, proj=proj)
colorbar(layer)
show() |
|