MeteoInfoLab脚本示例：MICAPS第一类数据

MeteoInfo

这里演示从MICAPS第一类数据中读取站点经纬度及某一个变量的值（比如温度），然后绘制散点图以及插值为格点数据绘制等值线图。打开MICAPS数据的命令是addfile_micaps()，然后就可以从文件变量里读取各个变量的数据。

1. f = addfile_micaps('D:/Temp/micaps/10101414.000')
   
2. data = f['Temperature'][:]
   
3. lon = f['Longitude'][:]
   
4. lat = f['Latitude'][:]
   
5. t = f.gettime(0)
   
6. #Plot
   
7. axesm(bgcolor=(204,255,255))
   
8. geoshow('country', edgecolor='k', facecolor=(255,251,195))
   
9. geoshow('cn_province', edgecolor='k')
   
10. levs = arange(0, 35, 2)
    
11. layer = scatterm(lon, lat, data, levs)
    
12. colorbar(layer)
    
13. yticks([20,30,40,50])
    
14. title('Temperature (' + t.strftime('%Y-%m-%d %H:00') + ')')
    
15. xlim(72, 136)
    
16. ylim(16, 55)

f = addfile_micaps('D:/Temp/micaps/10101414.000')<br /> data = f['Temperature'][:]<br /> lon = f['Longitude'][:]<br /> lat = f['Latitude'][:]<br /> t = f.gettime(0)<br /> #Plot<br /> axesm(bgcolor=(204,255,255))<br /> geoshow('country', edgecolor='k', facecolor=(255,251,195))<br /> geoshow('cn_province', edgecolor='k')<br /> levs = arange(0, 35, 2)<br /> layer = scatterm(lon, lat, data, levs)<br /> colorbar(layer)<br /> yticks([20,30,40,50])<br /> title('Temperature (' + t.strftime('%Y-%m-%d %H:00') + ')')<br /> xlim(72, 136)<br /> ylim(16, 55)

用站点数据绘制等值线图需要先将站点数据插值为格点数据。

1. f = addfile_micaps('D:/Temp/micaps/10101414.000')
   
2. data = f['Temperature'][:]
   
3. lon = f['Longitude'][:]
   
4. lat = f['Latitude'][:]
   
5. t = f.gettime(0)
   
6. #To grid data
   
7. x = arange(70, 140, 0.5)
   
8. y = arange(15, 58, 0.5)
   
9. gdata,gx,gy = griddata((lon, lat), data, xi=(x, y), method='idw')
   
10. #Plot
    
11. axesm()
    
12. geoshow('country')
    
13. geoshow('cn_province')
    
14. levs = arange(0, 35, 2)
    
15. layer = contourfm(x, y, gdata, levs)
    
16. colorbar(layer)
    
17. yticks([20,30,40,50])
    
18. title('Temperature (' + t.strftime('%Y-%m-%d %H:00') + ')')
    
19. xlim(72, 136)
    
20. ylim(16, 55)

f = addfile_micaps('D:/Temp/micaps/10101414.000')<br /> data = f['Temperature'][:]<br /> lon = f['Longitude'][:]<br /> lat = f['Latitude'][:]<br /> t = f.gettime(0)<br /> #To grid data<br /> x = arange(70, 140, 0.5)<br /> y = arange(15, 58, 0.5)<br /> gdata,gx,gy = griddata((lon, lat), data, xi=(x, y), method='idw')<br /> #Plot<br /> axesm()<br /> geoshow('country')<br /> geoshow('cn_province')<br /> levs = arange(0, 35, 2)<br /> layer = contourfm(x, y, gdata, levs)<br /> colorbar(layer)<br /> yticks([20,30,40,50])<br /> title('Temperature (' + t.strftime('%Y-%m-%d %H:00') + ')')<br /> xlim(72, 136)<br /> ylim(16, 55)

运行了f = addfile_micaps('D:/Temp/micaps/10101414.000')语句后，可以在Console里输入f获取数据文件信息，包括所有的变量名。


>>> f
File Name: D:/Temp/micaps/10101414.000
Description: diamond 1 10年10月14日14时地面填图
Time: 2010-10-14 14:00
Station Number: 2180
Number of Variables = 22
        Stid
        Longitude
        Latitude
        Altitude
        Grade
        CloudCover
        WindDirection
        WindSpeed
        Pressure
        PressVar3h
        WeatherPast1
        WeatherPast2
        Precipitation6h
        LowCloudShape
        LowCloudAmount
        LowCloudHeight
        DewPoint
        Visibility
        WeatherNow
        Temperature
        MiddleCloudShape
        HighCloudShape

virl

谢谢，太有用了

qxjy009

楼主我想用第一类数据格式画风矢量，是要用mipy.MeteoDataset.SetDrawType("barb")?吗，按照您这上面的格式读取了WindDirection和WindSpeed之后，要用什么命令画出barb呢

MeteoInfo

qxjy009 发表于 2016-12-14 16:33
楼主我想用第一类数据格式画风矢量，是要用mipy.MeteoDataset.SetDrawType("barb")?吗，按照您这上面的格式 ...

参考此贴：http://bbs.06climate.com/forum.php?mod=viewthread&tid=36171

风向杆用barbsm函数

媛女

第二类数据格式 高空的数据怎么读呢？

黑色青春

你好 micaps 中国一般有多少站点呢

Linda.W

王老师，请问Micaps第五类探空资料如何用MI读取？谢谢老师

MeteoInfo

Linda.W 发表于 2018-10-12 22:09
王老师，请问Micaps第五类探空资料如何用MI读取？谢谢老师

参考这个帖子：http://bbs.06climate.com/forum.php?mod=viewthread&tid=46470

Linda.W

谢谢王老师，，学习了。

chongzika

王老师我附件里面南京站点6小时降水有156之多，为什么画完图只有125不到（忽略标题）。
58238  118.80   32.00   12    2    8   90    3   26  -10    8    8
  156.0 9999    8  300   21.7    1.9   81   22.0 9999 9999    1    2 9999 9999