MeteoInfo介绍

MeteoInfo · 发表于 2012-10-3 14:15:22

haobang008 发表于 2012-10-3 09:44

登录/注册后可看大图

这次真心麻烦你了，呵呵

写了一个脚本将你提供的coastline数据生成一个线图层，基本可以肯定coastline数据不适合做mask，线条很零碎（见下图中的彩色线，每条零碎的线用不同颜色显示），想要形成封闭的区域几乎不可能。不知道你为什么这么执着要用coastline数据？我之前说过，MeteoInfo带的110m-land.shp很适合做海陆mask，你如果愿意的话可以试试。

具体代码：

#--------------------------------------------------------
# Author: Yaqiang Wang
# Date: 2012-10-3
# Purpose: Read and plot coastline data
# Note: Sample
#-----------------------------------------------------------
#---- Import clr and classes
import clr
clr.AddReferenceByPartialName("System")
clr.AddReferenceByPartialName("System.Drawing")
from System import *
from System.Drawing import *
from System.Collections.Generic import *
clr.AddReference("MeteoInfoC.dll")
from MeteoInfoC import *
from MeteoInfoC.Data.MapData import *
from MeteoInfoC.Geoprocess import *
from MeteoInfoC.Layer import *
from MeteoInfoC.Shape import *
from MeteoInfoC.Legend import *
from MeteoInfoC.Drawing import *
import os.path
#---- Create a line layer
print 'Create polyline layer...'
lLayer = VectorLayer(ShapeTypes.Polyline)
lLayer.LayerName = "Coast_Line"
lLayer.Visible = True
lLayer.EditAddField("ID", String)
#---- Read data file
print 'Read data file...'
fn = "D:\\Temp\\ascii\\data_1_bad.dat"
tf = open(fn)
plist = List[PointD]()
tf.readline()
i = 0
for aline in tf:
if aline.strip()[0:3] == 'nan':
aPolyline = PolylineShape()
aPolyline.Points = plist
shapeNum = lLayer.ShapeNum
if lLayer.EditInsertShape(aPolyline, shapeNum):
lLayer.EditCellValue("ID", shapeNum, shapeNum)
plist = List[PointD]()
print "Line %i" % (shapeNum)
else:
datalist = aline.split()
lat = float(datalist[1])
lon = float(datalist[0])
aPoint = PointD(lon, lat)
plist.Add(aPoint)
lLayer.UpdateLegendScheme(LegendType.SingleSymbol, "ID")
mipy.MapDocument.ActiveMapFrame.AddLayer(lLayer)
print 'Finished!'

复制代码

haobang008 · 发表于 2012-10-5 10:53:38

MeteoInfo 发表于 2012-10-3 14:15

登录/注册后可看大图

写了一个脚本将你提供的coastline数据生成一个线图层，基本可以肯定coastline数据不适合做mask，线条很零 ...

十分感谢版主，之所以执着着用coastline 的数据，是因为我手上的数据只有它的精度是最高的，可以精确到 1''左右，我的数据也是比较精确，不知道你说的110m-land.shp的精度有多大？

haobang008 · 发表于 2012-10-5 11:00:23

MeteoInfo 发表于 2012-10-3 14:15

登录/注册后可看大图

写了一个脚本将你提供的coastline数据生成一个线图层，基本可以肯定coastline数据不适合做mask，线条很零 ...

Yaqiang Wang 大神，十分感谢你的算法，但是我的Meteoinfo四最新版的，但是为什么找不到你说的自带的110m-land 文件呢？是不是你已经下载好了的？

MeteoInfo · 发表于 2012-10-5 14:55:28

haobang008 发表于 2012-10-5 11:00

登录/注册后可看大图

Yaqiang Wang 大神，十分感谢你的算法，但是我的Meteoinfo四最新版的，但是为什么找不到你说的自带的110 ...

在MeteoInfo安装路径的Map目录下。

haobang008 · 发表于 2012-10-5 15:03:29

MeteoInfo 发表于 2012-10-5 14:55

登录/注册后可看大图

在MeteoInfo安装路径的Map目录下。

嗯，找到陆地说，我先尝试一下呢

haobang008 · 发表于 2012-10-5 15:27:40

按照你的编码已经画出来了，但是我想把陆地和海洋分别设置为 0 和1 矩阵，不知道如何去做啊？

haobang008 · 发表于 2012-10-5 15:28:34

哦，哦对了，这是我做的图，想标记都不行呢

MeteoInfo · 发表于 2012-10-5 16:37:35

haobang008 发表于 2012-10-5 15:28

登录/注册后可看大图

哦，哦对了，这是我做的图，想标记都不行呢

可以参考此帖：MeteoInfo脚本示例：创建点图层并判断点是否在某区域内
http://bbs.06climate.com/forum.p ... 100&fromuid=106

haobang008 · 发表于 2012-10-5 19:26:56

MeteoInfo 发表于 2012-10-5 16:37

登录/注册后可看大图

可以参考此帖：MeteoInfo脚本示例：创建点图层并判断点是否在某区域内
http://bbs.06climate.com/forum. ...

正如前辈所说，是否在闭合图形之内，这个算法是有的，但是之前绘制的图形你想必也看到了，狠毒都是不闭合的线段，那我应该如何处理？还有就国界和湖泊的问题，那该如何处理呢？请赐教

cold_wq · 发表于 2012-10-24 15:43:43

楼主，wCcontour,TrajStat,EMDIS 这些下载了放在哪里？就可以直接用？（希望这个问题不太白痴）

		自动登录	找回密码
密码			立即注册