爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 8602|回复: 11

MeteoInfoLab脚本示例:Maskout图形

[复制链接]

新浪微博达人勋

发表于 2015-6-18 15:14:35 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 MeteoInfo 于 2020-5-30 21:59 编辑

Maskout通常有两种类型:Maskout图形和Maskout数据。这里是Maskout图形的示例。需要用shaperead读取地图数据形成图层作为Maskout图层(这里是中国的行政区域china.shp),然后用masklayer函数对数据图层进行屏蔽,第一个参数是用来Maskout的图层,第二次参数是一个图层列表(可以同时屏蔽多个图层)。
Image00843.png

脚本程序:
  1. f = addfile('D:/Temp/GrADS/model.ctl')
  2. ps = f['PS'][0,'10:60','60:140']
  3. axesm()
  4. m_china = geoshow('china', visible=False)
  5. geoshow('country', linecolor=(0,0,255))
  6. layer = contourfm(ps, 12)
  7. masklayer(m_china, [layer])
  8. title('Pressure')
  9. ylabel('Latitude')
  10. xlabel('Longitude')
  11. colorbar(layer)

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

新浪微博达人勋

 楼主| 发表于 2016-5-18 08:40:55 | 显示全部楼层
YF2015 发表于 2016-5-17 22:04
还想问一下怎么在软件中加点的名字比如下面这个图,加文字怎么改不来!只是text

点击工具栏中的“选择”工具,用鼠标选中文字,双击,在弹出的对话框中可以修改。
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2016-5-16 13:06:42 | 显示全部楼层
这个程序可以读写哪些文件,除了ctl{:eb513:}
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-5-16 13:34:43 | 显示全部楼层
YF2015 发表于 2016-5-16 13:06
这个程序可以读写哪些文件,除了ctl

参考这里:http://www.meteothinker.com/docs ... /data_tutorial.html
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-5-16 16:43:51 | 显示全部楼层
f = addfile('C:/AE_Land3.002/2008/AMSR_E_L3_DailyLand_V06_20080202.grd')
ps = f['PS'][0,(10,60),(60,140)]
axesm()
m_china = shaperead('C:/MeteoInfo/map/china.shp')
m_world = shaperead('C:/MeteoInfo/map/country1.shp')
geoshow(m_china, visible=False)
geoshow(m_world, linecolor=(0,0,255))
layer = contourfm(ps, 12)
masklayer(m_china, [layer])
title('Pressure')
ylabel('Latitude')
xlabel('Longitude')
colorbar(layer)
show()
我利用生成的GRD文件,不能运行
Traceback (most recent call last):
  File "<iostream>", line 1, in <module>
  File "C:\MeteoInfo\pylib\mipylib\minum.py", line 317, in addfile
    meteodata.openData(fname)
at java.lang.NumberFormatException.forInputString(Unknown Source)
at java.lang.Integer.parseInt(Unknown Source)
at java.lang.Short.parseShort(Unknown Source)
at java.lang.Short.parseShort(Unknown Source)
at org.meteoinfo.data.meteodata.arl.ARLDataInfo.readDataLabel(ARLDataInfo.java:527)
at org.meteoinfo.data.meteodata.arl.ARLDataInfo.canOpen(ARLDataInfo.java:195)
at org.meteoinfo.data.meteodata.MeteoDataInfo.openData(MeteoDataInfo.java:437)
at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)
at sun.reflect.NativeMethodAccessorImpl.invoke(Unknown Source)
at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)
at java.lang.reflect.Method.invoke(Unknown Source)
java.lang.NumberFormatException: java.lang.NumberFormatException: For input string: "DS"
出现这几行代码,问一下王老师,这是什么原因
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-5-16 17:30:23 | 显示全部楼层
meteinfo 可以打开shp文件吧,在哪里打开
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-5-16 18:12:56 | 显示全部楼层
YF2015 发表于 2016-5-16 16:43
f = addfile('C:/AE_Land3.002/2008/AMSR_E_L3_DailyLand_V06_20080202.grd')
ps = f['PS'][0,(10,60),(60 ...

GRD文件是什么格式?你怎么生成的?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-5-16 18:13:30 | 显示全部楼层
YF2015 发表于 2016-5-16 17:30
meteinfo 可以打开shp文件吧,在哪里打开

参考MeteoInfo的帮助文档:http://www.meteothinker.com/docs/meteoinfo/desktop/index.html
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-5-17 22:04:18 | 显示全部楼层
MeteoInfo 发表于 2016-5-16 18:13
参考MeteoInfo的帮助文档:http://www.meteothinker.com/docs/meteoinfo/desktop/index.html

还想问一下怎么在软件中加点的名字比如下面这个图,加文字怎么改不来!只是text
NMTXJDBH`]RHTB6ILTWJUQB.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-6-11 20:16:14 | 显示全部楼层
加文字的问题,我也很碰到以上麻烦,请王老师帮忙指点下好吗,谢市长
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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