- 积分
- 56367
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-6-21
- 最后登录
- 1970-1-1

|
发表于 2015-6-30 14:55:13
|
显示全部楼层
本帖最后由 MeteoInfo 于 2015-7-3 16:28 编辑
也可以把经纬线和标注加上:
 - #Add data file
- fn = 'D:/Temp/hdf/FY2C_CLC_MLT_NOM_20070730_1800.hdf'
- f = addfile(fn)
- #Get data variable
- v = f['FY2C_Hourly_Cloud_Classification']
- #Set x/y
- x = linspace(-5750000.0,5750000.0, 2288)
- y = linspace(-5750000.0,5750000.0, 2288)
- #Set x/y dimensions
- v.setdim('Y', y, 0)
- v.setdim('X', x, 1)
- #Get data array
- cloud = v[:,:]
- #Plot
- ax, proj = axesm(proj='geos', lon_0=104.5, h=35785864, gridlabel=True, gridline=True, frameon=False)
- mlayer = shaperead('D:/Temp/map/country1.shp')
- geoshow(mlayer, edgecolor=(0,0,255))
- layer = imshowm(x, y, cloud, proj=proj, cmap='hotcold_18lev')
- colorbar(layer)
可以用鼠标在图形区域拉框放大,在图形区域点击鼠标右键,在右键菜单中点击‘Undo zoom’返回原图形范围。
|
|