MeteoInfoLab脚本示例：数据投影-FLEXPART

MeteoInfo

FLEXPART是一个类似HYSPLIT的扩散模式，它输出的netcdf文件参照了WRF，可惜全局属性没有写全，比如只有一个投影名称（例如Lambert），没有相关的投影参数：中央经度，标准纬度等等。必须查阅WRF的头文件才能重建投影（为什么不照猫画虎把属性写全呢？）。数据的经纬度坐标是有的，但在Lambert投影下的坐标没有，可以通过projectxy函数获得投影下的x, y坐标，其中的lon, lat是数据左下角的经纬度。


脚本程序：

1. f = addfile('D:/Temp/nc/header_d01.nc')
   
2. topov = f['TOPOGRAPHY']
   
3. topo = topov[:,:]
   
4. #Set projection
   
5. fromproj = projinfo(proj='longlat')
   
6. toproj = projinfo(proj='lcc', lon_0=117.0092, lat_0=40.11041, lat_1=30.0, lat_2=60.0)
   
7. #Set x/y
   
8. lon = 112.14053344726562
   
9. lat = 35.62180709838867
   
10. xn = 216
    
11. yn = 240
    
12. dx = 4000
    
13. dy = 4000
    
14. x, y = projectxy(lon, lat, xn, yn, dx, dy, toproj, fromproj)
    
15. #Plot
    
16. axesm(projinfo=toproj, griddx=2, griddy=2, gridline=True)
    
17. mlayer = shaperead('D:/Temp/map/bou2_4p.shp')
    
18. geoshow(mlayer)
    
19. layer = imshowm(x, y, topo, 50, proj=toproj)
    
20. title('FLEXPART Topography')
    
21. colorbar(layer)
    
22. show()

f = addfile('D:/Temp/nc/header_d01.nc')<br /> topov = f['TOPOGRAPHY']<br /> topo = topov[:,:]<br /> #Set projection<br /> fromproj = projinfo(proj='longlat')<br /> toproj = projinfo(proj='lcc', lon_0=117.0092, lat_0=40.11041, lat_1=30.0, lat_2=60.0)<br /> #Set x/y<br /> lon = 112.14053344726562<br /> lat = 35.62180709838867<br /> xn = 216<br /> yn = 240<br /> dx = 4000<br /> dy = 4000<br /> x, y = projectxy(lon, lat, xn, yn, dx, dy, toproj, fromproj)<br /> #Plot<br /> axesm(projinfo=toproj, griddx=2, griddy=2, gridline=True)<br /> mlayer = shaperead('D:/Temp/map/bou2_4p.shp')<br /> geoshow(mlayer)<br /> layer = imshowm(x, y, topo, 50, proj=toproj)<br /> title('FLEXPART Topography')<br /> colorbar(layer)<br /> show()

MeteoInfo

kakame08 发表于 2016-1-8 09:12
才看到 这个真的是好东西啊
请问这个是对FLEXPART输出结果的显示吗？

一楼不是说得很清楚了吗

MeteoInfo

kakame08 发表于 2016-1-8 11:36
抱歉 没仔细看图
是不是只要是.nc格式的 都可以用MeteoInfo你这个脚本画图呢？
这个脚本是针对 heade ...

MeteoInfo原则上支持所有的netCDF格式数据。

lqhlqhlqh

haodongxi1!!!

刘勇的下铺

楼主你好，我输出的这个header文件好像不是nc文件，想问一下这样应该怎样处理，还是在程序的哪个地方能控制输出的格式呢？谢谢。

猪妹妹要去北京

刘勇的下铺 发表于 2017-4-28 12:00
楼主你好，我输出的这个header文件好像不是nc文件，想问一下这样应该怎样处理，还是在程序的哪个地方能控制 ...

好像FLEXPART——WRF耦合版本的输出格式是nc的，FLEXPART直接输出的是二进制文件，官网新出了一个补丁小程序，加进去可以控制输出文件类型为nc，但是其他的一些程序貌似也要改动，你可以试试

对方的凤飞飞

猪妹妹要去北京 发表于 2018-3-30 08:54
好像FLEXPART——WRF耦合版本的输出格式是nc的，FLEXPART直接输出的是二进制文件，官网新出了一个补丁小 ...

请问是什么补丁，能详细介绍一下，我怎么没有看到呢？