- 积分
- 57092
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-6-21
- 最后登录
- 1970-1-1
|
发表于 2019-8-8 21:12:43
|
显示全部楼层
根据你提供的正确投影信息可以在脚本程序中重新计算投影和数据坐标:
 - fn = 'D:/Temp/grib/Z_NAFP_C_BCSH_20190724000000_P_surface-warms-f36.grib'
- f = addfile(fn)
- vproj = f['LambertConformal_Projection']
- proj = projinfo(proj='lcc', lat_0=35, lon_0=105, lat_1=30, lat_2=60)
- llproj = projinfo()
- lon00 = 78.30664
- lat00 = 7.29969
- x0, y0 = project(lon00, lat00, llproj, proj)
- nx = 759
- ny = 599
- delta = 9000 # 9km
- x = arange1(x0, nx, delta)
- y = arange1(y0, ny, delta)
- data = f['Temperature_surface'][0,:,:]
- #Plot
- axesm(projinfo=proj)
- geoshow('country', edgecolor='k')
- layer = imshowm(x, y, data, 20, proj=proj)
- colorbar(layer)
|
|