打开包含站点填图的站点数据文件(比如micaps 1)之后,用文件对象的smodeldata函数获取StationModel数据对象,然后用stationmodel函数绘制站点填图图层。

脚本程序:

#Set data folders
basedir = 'D:/MyProgram/Distribution/java/MeteoInfo/MeteoInfo'
datadir = os.path.join(basedir, 'sample/MICAPS')
mapdir = os.path.join(basedir, 'map')
#Read shape files
lworld = shaperead(os.path.join(mapdir, 'country1.shp'))
lchina = shaperead(os.path.join(mapdir, 'bou2_4p.shp'))
#Read station data
f = addfile_micaps(os.path.join(datadir, '10101414.000'))
data = f.smodeldata()
#Plot
axesm(bgcolor=(204,255,255))
geoshow(lworld, facecolor=(255,251,195))
geoshow(lchina, edgecolor='k')
layer = stationmodel(data)
#colorbar(layer)
xlim(72, 136)
ylim(16, 55)

转载于:https://www.cnblogs.com/yaqiang/p/4869324.html

MeteoInfoLab脚本示例:站点填图相关推荐

  1. MeteoInfoLab脚本示例:多Y轴图

    数据范围相差比较大的数据序列进行对比的时候多Y轴图就很重要了.MeteoInfoLab中提供了一个twinx函数来根据已有的坐标系(Axes)生成一个新的Axes,这个命令会使得已有的Axes不绘制右 ...

  2. MeteoInfoLab脚本示例:CloudSAT Swath HDF数据

    读取CloudSAT HDF Swath数据,绘图分上下两部分,上面是时间和高度维的Radar Reflectivity Factor二维图,下面是卫星轨迹图. 示例程序: # Add file f ...

  3. MeteoInfoLab脚本示例:Streamline流线图

    绘制Stramline流线图的函数是streamline,需要两个变量(U/V分量或者风向/风速). 脚本程序: f = addfile('D:/Temp/GrADS/model.ctl') u = ...

  4. MeteoInfoLab脚本示例:TRMM 3B43 HDF数据

    TRMM 3B43是卫星观测月平均降水量产品,是HDF的格点数据.需要注意的是数据中降水变量维的顺序里经度维在前纬度维在后,这与通常的设置(纬度维在前经度维在后)相反,需要对获取的二维数组进行转置,使 ...

  5. MeteoInfoLab脚本示例:OMI Grid HDF数据

    OMI卫星格点数据的例子,全球臭氧柱总量分布. 脚本程序: #Add data file folder = 'D:/Temp/hdf/' fns = 'OMI-Aura_L3-OMTO3e_2005m ...

  6. MeteoInfoLab脚本示例:MODIS Sinusoidal投影HDF数据

    MODIS卫星很多陆面数据都是Sinusoidal投影,数据被分为一个个10*10度(赤道地区)的瓦片(http://modis-land.gsfc.nasa.gov/MODLAND_grid.htm ...

  7. MeteoInfoLab脚本示例:获取气团轨迹每个节点的气象数据

    读取HYSPLIT输出的轨迹数据文件和相应时间的气象数据文件,生成轨迹图层,循环每条轨迹的节点,读出该节点的经度.纬度.气压.时间,通过对气象数据插值获得该节点的气象数据. 脚本程序: #------ ...

  8. MeteoInfoLab脚本示例:OMI Swath HDF数据

    这个例子读取OMI卫星Swath数据中的CloudFaction变量并绘图. 脚本程序: #Add data file folder = 'D:/Temp/hdf/' fns = 'OMI-Aura_ ...

  9. MeteoInfoLab脚本示例:AIRS Swath HDF数据

    例子中的AIRS Swath HDF数据在Polar Stereographic(南极)投影中接近矩形,需要先从数据中读出经纬度及相关数据数组,利用surfacem函数绘制Swath数据(散点),在s ...

最新文章

  1. IDEA IntelliJ 如何设置web项目的热部署
  2. 12.QT线程的两种启动方式
  3. 如何让你用 Python 年薪 40 万?答案早就写到 JD 上了
  4. 函数 单片机glint_GL 库函数
  5. 1.2鼠标移入移出改变背景色和其他大小样式
  6. CRC32算法详细推导(2)
  7. Android\OPhone自定义视图(View)
  8. Thinkphp3.2+ 微信小程序图片上传
  9. 为什么.class文件查看不了_恕我直言,这可能是你见过最详细的class文件结构分析
  10. 第九届蓝桥杯c语言b组试题+部分答案
  11. UE4.26官方文档网页浏览录屏打包下载版
  12. C#,入门教程——一个优秀程序员的修炼之道
  13. linux系统支持网银吗,Linux系统能使用网银吗?
  14. Ubuntu下安装Genymotion安卓模拟器
  15. contents()与children()的用法和区别
  16. 如何做好一个产品经理
  17. 百度地图画扇形区域覆盖(大小方向颜色透明图可调)
  18. 电子行业测试软件,CAM软件技术在电子行业的应用测试验证
  19. 使用FreeMarker自动生成Word文档
  20. 展现AI与自动化测试技术之间的神奇化学反应

热门文章

  1. 假期第一天,第一次见丈母娘~
  2. Eth-Trunk负载分担不均怎么办,如何通过Hash算法实现负载分担?
  3. vue.runtime.esm.js?2b0e:619 [Vue warn]: Duplicate keys detected: ‘xxx‘. This may cause an update err
  4. Python PIL Image的使用
  5. 量子信息-学习记录13
  6. WordCloud库的安装教程
  7. 莫纳什大学计算机专业研究生在哪个校区,2020年莫纳什大学信息技术学院具体地址及在哪个校区...
  8. 88steamCSGO即开即取回的开箱网站,CSGO皮肤交易平台
  9. NO.ONE进程、进程、线程、线程——阿古兽高级超级终极究极进化暴龙兽喷火暴龙兽机械暴龙兽战斗暴龙兽
  10. jeeplus框架代码生成器功能的使用