请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 32735|回复: 61

MeteoInfoLab脚本:TXT站点数据绘制等值线填色图

  [复制链接]

新浪微博达人勋

发表于 2017-5-16 08:43:21 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 静言_GRMC 于 2017-5-16 09:21 编辑

参考了版主发的脚本示例(http://bbs.06climate.com/forum.p ... 6422&extra=page%3D1),版主的是用micaps资料,我也依葫芦画瓢做了一个,数据存在txt(经、纬度、PM2.5浓度),分享给大家,有代码多余或者不规范的请多多指教~

代码入下:
  1. #Get file names
  2. fn = 'f:/meteoinfoJS/AQI.txt'

  3. ncol = numasciicol(fn)
  4. nrow = numasciirow(fn)
  5. a = asciiread(fn,shape=(nrow,ncol))
  6. lon = a[:,0]
  7. lat = a[:,1]
  8. pm= a[:,2]

  9. #To grid data
  10. x = arange(112.9, 114.3, 0.2)
  11. y = arange(22.5, 24.2, 0.2)
  12. gtemp,gx,gy = griddata((lon, lat), pm, xi=(x, y), method='idw', radius=0.8)
  13. #Plot
  14. axesm()
  15. bou1_layer = shaperead('f:/MeteoInfo/map/guangzhou_county_new.shp')
  16. mlayer = shaperead('f:/MeteoInfo/map/guangzhou_county_new.shp')
  17. geoshow(bou1_layer, edgecolor='lightgray')
  18. geoshow(mlayer, visible=False)

  19. levs = [0, 25, 35, 50, 75, 95, 115, 130, 150,200]
  20. cols = [(255,255,255),(0,255,0),(127,255,0),(255,255,0),(255,215,0),(255,128,0),(255,97,0), \
  21.     (255,0,0),(176,23,31),(135,38,87),(255,0,255)]
  22. #layer = contourfm(x, y, gtemp,20)  

  23. layer = contourfm(x, y, gtemp,levs,colors=cols)
  24. #slayer = scatterm(lon, lat,pm,colors=['k'], size=10)
  25. slayer = scatterm(lon, lat,pm,levs,colors=cols, size=8)

  26. masklayer(mlayer, [layer])
  27. xlim(112.9, 114.1)
  28. ylim(22.5, 24)
  29. title(u'广州市空气质量站点PM2.5浓度',fontname=u'黑体',fontsize=20,bold=False,color='blue')   
  30. text(113.3, 23.9, u'2017年5月6日13时', fontname=u'黑体', fontsize=16)
  31. colorbar(layer)
复制代码
PM2.5.png

评分

参与人数 1威望 +2 金钱 +10 贡献 +8 体力 +20 收起 理由
MeteoInfo + 2 + 10 + 8 + 20 赞一个!

查看全部评分

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-5-16 09:53:03 | 显示全部楼层
哈哈,不错不错,加油加油。
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2017-5-16 09:30:11 | 显示全部楼层
欢迎这样的经验分享贴!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-5-16 10:00:15 | 显示全部楼层
MeteoInfo 发表于 2017-5-16 09:30
欢迎这样的经验分享贴!

谢谢王老师点赞~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-5-16 10:07:02 | 显示全部楼层
很好很强大
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-5-26 14:49:07 | 显示全部楼层
请问广州市或者广东省的地图去哪里下载呢?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-5-26 14:52:04 | 显示全部楼层
跟大牛要到啦
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-5-27 10:36:46 | 显示全部楼层

那就好啊   ~~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-5-27 10:36:51 | 显示全部楼层

那就好啊   ~~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-5-27 10:40:31 | 显示全部楼层
AQI.txt'能不能把AQI.txt'给我们看下
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表