爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4884|回复: 9

[vbscript] 只网格化有数据区域

[复制链接]

新浪微博达人勋

发表于 2014-6-9 09:05:50 | 显示全部楼层 |阅读模式

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

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

x
几年前,因为遇到大量有非规则边界的数据,每次要做白化都头痛,因为做个bln文件太烦了。譬如 20140606135925.jpg


于是跑到goldensoftware问,有没有快速、自动的办法按数据生成白化边界。有大神说白化就是个屁,边界文件算什么东西,你穿越到未来用 surfer10 就有“Blank grid outside contvex hull of data”选项了。


20140609082809.jpg

当时是surfer8,所以享受不到这个以凸包检测来自动白化无数据区的先进功能。

然后高手们群策群力,说白化边界文件是没必要的,用一个三角剖分可以解决。先正常网格化(就是你希望用什么方法、什么参数),再用三角剖分做一次网格化,接着用网格数学,把正常网格化的网格选为A,把三角剖分网格化的网格选为B,公式呢就写 “A+B-B”,结果很漂亮。

20140606140007.jpg

但是我当时的数据,都是有凹多边形边界啊,上面的方法只能搞定凸多边形。

好吧,我最终选用的遮罩法横空出世了。大神们说了,管你什么形状,还是一个 Mask 解决啊。跟上面一样,先还是正常网格化,得到网格A。然后用 Data Metrics 做第二次网格化,这次要选择点位计数,得到网格B。
20140609084952.jpg


网格 B 还不能直接用,用网格数学处理一下,公式“if(a<1,1.70141e+038,0)”,因为网格 B 的节点值就是该节点周围的数据个数,所以前面的公式就是把没数据的节点白化(1.70141e38是surfer默认的白化值)。

再来一次网格数学,这次选择网格A和处理过的网格B,用公式 A+B 解决。
20140606151726.jpg


根据这些讨论,我写了脚本自动处理(因为当时有大量需要白化的数据)。结果的好坏要看选用的搜索半径,当然网格A的质量是至关重要的。

GridInside.zip (7.29 KB, 下载次数: 16)

评分

参与人数 3金钱 +58 贡献 +7 收起 理由
陌小默的 + 18 + 2
chengxf + 20 很给力!
言深深 + 20 + 5 赞一个!

查看全部评分

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

新浪微博达人勋

发表于 2014-6-12 09:32:24 | 显示全部楼层
我们自己开发的国产地震解释软件GEOEAST轻松可以自己生成数据边界并自动处理边界外网格,算法是我开发的。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-6-12 10:31:23 | 显示全部楼层
好人一僧平安啊
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-6-12 13:51:53 | 显示全部楼层
本帖最后由 郭小侠V 于 2014-6-12 21:02 编辑

求分享啊,2楼!!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-6-14 08:02:17 | 显示全部楼层
这个不错,赞一个
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-6-14 09:43:30 | 显示全部楼层
求分享一下!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-1-27 09:12:18 | 显示全部楼层
看上去不错,试试
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-1-27 09:54:24 | 显示全部楼层
棒啊~攒够钱了过来下载!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-5-8 08:34:24 | 显示全部楼层
学习了,赞一下
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-11-26 14:47:41 | 显示全部楼层
好厉害,真是高手啊,学习了
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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