MeteoInfo介绍

MeteoInfo · 发表于 2012-10-1 22:32:30

haobang008 发表于 2012-10-1 21:07

登录/注册后可看大图

请问你有什么好的办法吗？

格式是可以转换的，参考此帖：全球高精度海岸线数据的处理
http://bbs.06climate.com/forum.p ... 115&fromuid=106

重要的是数据必须是封闭的。不知道你要做多大范围的mask，如果范围大的话（比如全球的），没必要用精度这么高的数据，计算起来会很慢。MeteoInfo软件带的110m-land.shp就是很好的做海陆mask的数据。

haobang008 · 发表于 2012-10-2 16:20:15

MeteoInfo 发表于 2012-10-1 22:32

登录/注册后可看大图

格式是可以转换的，参考此帖：全球高精度海岸线数据的处理
http://bbs.06climate.com/forum.php?mod=vie ...

你好，我下载了数据以后发现一个问题，每组数据都是用nan 分隔开，有些是封闭的，有些则不是的，如果都是闭合的话，我自己完全可以用 Fortran来做，但一半是闭合的数据，而另一半则不是闭合的数据，而且其中还掺杂了国界线呢，所以请问楼主的具体思路是如何来处的呢？

haobang008 · 发表于 2012-10-2 16:21:34

至于 meteoinfo 我的兴趣倒不是很大，只是想知道你的具体思路是如何处理coastline extractor 的数据的呢？拜托了，多谢

MeteoInfo · 发表于 2012-10-2 16:43:32

haobang008 发表于 2012-10-2 16:21

登录/注册后可看大图

至于 meteoinfo 我的兴趣倒不是很大，只是想知道你的具体思路是如何处理coastline extractor 的数据的呢？拜 ...

相关的帖子已经说得很清楚了，如果不愿意尝试，也无可奈何。

haobang008 · 发表于 2012-10-2 16:55:29

MeteoInfo 发表于 2012-10-2 16:43

登录/注册后可看大图

相关的帖子已经说得很清楚了，如果不愿意尝试，也无可奈何。

我只是想知道你的具体思路，而不是说MeteoInfo 的问题，就是你的思路是如何做的，如何将不是闭合的图形设置的呢？

haobang008 · 发表于 2012-10-2 16:56:12

你的算法是如何的啊，其实无论 Meteoinfo还是Fortran，本质上不都是一样的吗？

MeteoInfo · 发表于 2012-10-3 09:15:39

haobang008 发表于 2012-10-2 16:56

登录/注册后可看大图

你的算法是如何的啊，其实无论 Meteoinfo还是Fortran，本质上不都是一样的吗？

算法其实就是判断一个点是否在多边形内，可以用射线法判断，具体你可以google一下，应该有很多相关的代码。

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

嗯，这才说道点子上，我现在纠结的是我要处理的点并不是闭合的图形，我把具体的数据粘到下面，你可以看下：文件中每一组数据用nan 隔开，但是至于为什么隔开，可能是为了画图方便吧，每组nan 隔开的书不会组成一个闭合的图形，那我应该怎么处理了？（数据参见附加呢）

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

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

登录/注册后可看大图

算法其实就是判断一个点是否在多边形内，可以用射线法判断，具体你可以google一下，应该有很多相关的代码 ...

这次真心麻烦你了，呵呵

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!'

复制代码

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