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

气象家园

 找回密码
 立即注册

新浪微博登陆

只需一步, 快速开始

QQ登录

只需一步,快速开始

搜索
查看: 824|回复: 1

MeteoInfoLab脚本示例:区域填色图

[复制链接] |关注本帖

新浪微博达人勋

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

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

x
读取美国区域的温度格点数据和美国州的行政区域shape图层,给图层增加温度字段'temp',利用维度数组的maskout函数获取每个州的平均温度,并放入温度字段中,创建一个温度数据的图例,利用geoshow函数显示区域填色图。

  1. #Read temporature data from a surfer grid data
  2. f = addfile_surfer('D:/Temp/ascii/usgrid.dat')
  3. tdata = f['var'][:,:]

  4. #Read US states layer from shape file
  5. us = shaperead('D:/Temp/map/states.shp')

  6. #Add temp field
  7. us.addfield('temp', 'float')

  8. #Average temporature for each state and add to the temp field
  9. for i in range(us.shapenum()):
  10.     rpoly = us.shapes()
  11.     mdata = tdata.maskout(rpoly)
  12.     tave = mdata.ave()
  13.     us.setcellvalue('temp', i, tave)

  14. #Plot
  15. axesm()
  16. geoshow('country', facecolor='lightgray', edge=False)
  17. levs = arange(270, 300, 2)
  18. cols = makecolors(len(levs)+1)
  19. ls = makesymbolspec('polygon', field='temp', levels=levs, colors=cols,
  20.     edge=True, edgecolor='gray')
  21. geoshow(us, symbolspec=ls)
  22. xlim(-128, -65)
  23. ylim(24, 50)
  24. title('Average temporature of states')
  25. colorbar(us)


temp_ave_states.png
密码修改失败请联系qq:937062711

新浪微博达人勋

发表于 2018-3-10 08:15:54 来自手机 | 显示全部楼层 |取消关注该作者的回复
王老师,可以用站点资料吗?即每个县一个数据带表平均值
密码修改失败请联系qq:937062711
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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