请问一下华东的模式数据不能正常打开是什么情况呢？

onekisy

数据文件用的兰伯特投影
取到的参数
int LambertConformal_Projection);
                LambertConformal_Projection: grid_mapping_name = "lambert_conformal_conic"
                LambertConformal_Projection: latitude_of_projection_origin = 0.029999999329447746
                LambertConformal_Projection: longitude_of_central_meridian = 0.10499999672174454
                LambertConformal_Projection: standard_parallel = 0.029999999329447746, 0.05999999865889549
                LambertConformal_Projection: earth_radius = 6371229.0

和华东模式的说明参数不一致呢
兰伯特投影(LambertConformal)
CEN_LAT=35.0；CEN_LON=105.0；TRUELAT1=30.0;TRUELAT2=60.0;STAND_LON=105.0;
XLAT(0,0)=7.29969,XLONG(0,0)=78.30664

MeteoInfo

这是数据本身的问题，可能是写文件的程序有问题

MeteoInfo

根据你提供的正确投影信息可以在脚本程序中重新计算投影和数据坐标：

1. fn = 'D:/Temp/grib/Z_NAFP_C_BCSH_20190724000000_P_surface-warms-f36.grib'
   
2. f = addfile(fn)
   
3. vproj = f['LambertConformal_Projection']
   
4. proj = projinfo(proj='lcc', lat_0=35, lon_0=105, lat_1=30, lat_2=60)
   
5. llproj = projinfo()
   
6. lon00 = 78.30664
   
7. lat00 = 7.29969
   
8. x0, y0 = project(lon00, lat00, llproj, proj)
   
9. nx = 759
   
10. ny = 599
    
11. delta = 9000    # 9km
    
12. x = arange1(x0, nx, delta)
    
13. y = arange1(y0, ny, delta)
    
14. data = f['Temperature_surface'][0,:,:]
    
15. 
    
16. #Plot
    
17. axesm(projinfo=proj)
    
18. geoshow('country', edgecolor='k')
    
19. layer = imshowm(x, y, data, 20, proj=proj)
    
20. colorbar(layer)

fn = 'D:/Temp/grib/Z_NAFP_C_BCSH_20190724000000_P_surface-warms-f36.grib'<br /> f = addfile(fn)<br /> vproj = f['LambertConformal_Projection']<br /> proj = projinfo(proj='lcc', lat_0=35, lon_0=105, lat_1=30, lat_2=60)<br /> llproj = projinfo()<br /> lon00 = 78.30664<br /> lat00 = 7.29969<br /> x0, y0 = project(lon00, lat00, llproj, proj)<br /> nx = 759<br /> ny = 599<br /> delta = 9000    # 9km<br /> x = arange1(x0, nx, delta)<br /> y = arange1(y0, ny, delta)<br /> data = f['Temperature_surface'][0,:,:]<br /> <br /> #Plot<br /> axesm(projinfo=proj)<br /> geoshow('country', edgecolor='k')<br /> layer = imshowm(x, y, data, 20, proj=proj)<br /> colorbar(layer)

onekisy

好的 谢谢。

onekisy

MeteoInfo 发表于 2019-8-8 21:12
根据你提供的正确投影信息可以在脚本程序中重新计算投影和数据坐标：

请问一下这种格点数据meteoinfo有没有方法可以直接转换成经纬度格点