Python读取micaps mdfs数据

eeeee

一直没有看到有人读取mdfs的站点数据，于是今天就按照数据说明自己写了一个小程序来读MDFS数据。
传入文件路径实例化即可。

1. mdfs = MDFS_Station(文件路径)
   

复制代码

MDFS_Station类含有的属性：
data_dsc：描述信息
level：数据层面
level_dsc：层面描述信息
utc_time：数据时间（世界时），为datetime.datetime类型
data：主体数据，为dict类型，含有的键值有'ID'（站号），'Lon'（经度），'Lat'（纬度）以及数据里对应的物理量数据的编号（见dtype_link.xml）
演示：

1. from read_mdfs import MDFS_Station
   
2. import matplotlib.pyplot as plt
   
3. import matplotlib.colors as ccolor
   
4. import matplotlib.cm as cmx
   
5. import cartopy.crs as ccrs
   
6. import cartopy.feature as cfeature
   
7. 
   
8. x = MDFS_Station(r'D:\20181028160000.000')
   
9. lon = x.data['Lon']
   
10. lat = x.data['Lat']
    
11. var = x.data[603]
    
12. 
    
13. cm = plt.get_cmap('jet')
    
14. cNorm = ccolor.Normalize(vmin=min(var), vmax=max(var))
    
15. scalarMap = cmx.ScalarMappable(norm=cNorm, cmap=cm)
    
16. c = scalarMap.to_rgba(var)
    
17. ax = plt.axes(projection=ccrs.PlateCarree())
    
18. ax.scatter(lon, lat, s=5, c=c)
    
19. ax.coastlines(resolution='10m')
    
20. ax.add_feature(cfeature.BORDERS)
    
21. plt.show()

复制代码

————————————更新（2018.11.4）——————————————
加入格点数据的读取

1. from read_mdfs import MDFS_Grid
   
2. import numpy as np
   
3. import matplotlib.pyplot as plt
   
4. import cartopy.crs as ccrs
   
5. import cartopy.feature as cfeature
   
6. import metpy.calc as mpcalc
   
7. from metpy.units import units
   
8. filepath = r'D:\18110408.000'
   
9. f = MDFS_Grid(filepath)
   
10. lon = f.data['Lon']
    
11. lat = f.data['Lat']
    
12. norm = f.data['Norm']
    
13. dir = f.data['Direction']
    
14. u, v = mpcalc.wind_components(norm * units('m/s'), dir * units.degree)
    
15. ax = plt.axes(projection=ccrs.PlateCarree())
    
16. ax.streamplot(lon, lat, u, v, density=3)
    
17. ax.coastlines(resolution='10m')
    
18. ax.add_feature(cfeature.BORDERS)
    
19. plt.show()

复制代码

注：data里面向量数据的键值是Norm和Direction，网格数据是Grid


————————————更新（2018.11.9）——————————————
修复了读取站点数据时返回要素长度不同的bug。

————————————更新（2018.12.7）——————————————
添加存储为NetCDF文件的功能

1. from read_mdfs import MDFS_Station, NetCDFWriter
   
2. x = MDFS_Station(r'D:\\20181028160000.000')
   
3. a = NetCDFWriter('D:\\1.nc')
   
4. a.load_data(x)
   
5. a.close()

复制代码

游客，如果您要查看本帖隐藏内容请回复

苏小眠

请问读取mdfs的那个库在哪里下载啊

fujianbeige

学习学习，谢谢分享

eeeee

xiahou90 发表于 2019-2-26 11:45
怎么才能批量下载导出mdfs的数据到本地呢？

你可以在micaps里面操作，或者自己写程序下载

飞雪残剑

请问楼主，怎么才能下载mdfs的数据呢？

xiaohu1643

支持楼主！！非常实用{:5_213:}

淡水河

我看看哈

tulalang

楼主这个写的很简便啊，看看

随缘

kankankankankan

月漫星空

多谢楼主分享！！

ysc930802

很给力的楼主，最近也在学这个，有现成的很好。

tulalang

楼主实例中的20181028160000.000这个数据是怎么获得的呢？GDS么

道狼52203

O(∩_∩)O哈哈~，好东西，学习学习！！

zhouling

学习一个