请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 38088|回复: 48

[资源分享] Surfer软件图形的体积面积计算

  [复制链接]

新浪微博达人勋

发表于 2015-6-6 22:12:58 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 chengxf 于 2015-6-6 22:17 编辑

Surfer软件图形的体积面积计算
等值线图是Surfer软件功能的最主要应用。对于一个已知的X、Y、Z数据集,利用丰富的插值算法,可以得到该数据集的等值线分布。例如下图,就是四川省和重庆市(当时尚未分开)的一次降水分布(降水数据是伪造的)。现在我们的题目是计算这次降水过程中雨量大于等于70mm的区域面积有多大。图上我们直观地看到被绿色部分包围的就是我们要求的区域。
                               2015-6-6 10-30-29.jpg
首先我们来了解一下Surfer软件中关于面积的定义,面积分平面面积和曲面面积,每一种面积里面又分正(裁切)面积和负(填充)面积。Surfer 求面积时需要两个面,称为上表面和下表面,这两个表面其中一个一定是曲面,另一个可以是曲面也可以是平面。例如你有了降水的分布等值线图,这就是一个曲面,然后你确定一个平面的值,就是 Z 值,假定你需要求大于70mm降水的面积,这个平面的 Z 值就是70mm那根等值线的值(70)。然后求这个平面和曲面之间的正(裁切)面积。
曲面面积很好理解,它就是曲面整个的面积,平面面积是被定义为曲面投影到一个平面上的面积,平面面积总是小于曲面面积。我们这个题目中,只要求求平面面积。左下图是上述降水分布的3D曲面图,注意不要将它看成是地形图,这里越是高的地方表示雨量越大罢了。
2015-6-6 11-04-07.jpg     2015-6-6 13-26-26.jpg
上右图是加入Z=70的平面图以后的3D曲面图,这个平面切割出来的部分,如果将它们投影到这个平面上,那么就称之为正平面面积,也叫裁切平面面积。除此以外的这个平面以下的低于70mm的部分,它的(向上的)投影,就称之为负平面面积,也叫填充平面面积。我们将等值线分布以及边界图叠加上来,可以看到下面的图形,就更加清楚了。实际上,我们将第一张图稍加修正,就可以看到绿色部分是正平面面积部分,其他白色部分是负平面面积。平面面积里面还有一个白化面积的概念,就是说当上表面和下表面重叠的时候,这一部分就是白化面积。在我们这个实例中,四川省和重庆市边界以外的地区,就是白化面积,一般情况下,这个面积是没有什么用途的。
   2015-6-6 13-47-57.jpg 2015-6-6 13-54-48.jpg    
说到这里也许会有人问,到底该如何计算面积呢?我要告诉大家的是计算工作是很烦琐的,应该交由计算机去完成,我们需要的是搞懂这些术语的定义以便将来可以看懂Surfer软件提供的报告。
搞懂了面积的定义,基本上也就搞懂了体积的定义,体积就没有平面曲面之分了,对应上面的正面积,就是正体积,也称裁切体积;对应负面积,就是负体积,也称填充体积。关于体积的论述,我们不能以上述降水来举例,那是不恰当的。当地震山体塌方形成堰塞湖的时候,我们根据上述原理,立刻可以求出负体积、填充体积来,体积出来了,蓄水量多少吨也就知道了。当我们需要削去某个山头,开凿一条高速公路时,立刻可以计算出这个山头的正体积、裁切体积,工程的土方量也就出来了。这里关于体积的叙述就不多说了,原理和概念是一样的。
在Surfer中,运行“网格(Grid)|体积(Volume)”,在“打开文件”对话框中,选中我们上面作为范例的网格文件,这时弹出一个“网格体积”对话框,可以看到,刚刚选中的文件已经作为“上表面”列在对话框中,而“下表面”我们选中“常数”,在Z值的编辑框中输入“70”,就是说将Z=70这个平面作为下表面。如下图。如果你希望了解一下上表面的网格状况,可以点击上表面网格文件名右边的 2015-6-6 19-26-40.jpg   按钮,将显示网格的信息如下。
2015-6-6 19-27-48.jpg 2015-6-6 19-28-38.jpg
其他取默认值,点击“确认”按钮,一份网格面积体积计算报告会呈现在你面前。
这份报告就是我们要的结果,下面将这个报告分成两部分贴在下面,由于我这里已经翻译成汉字,所以比较好理解。我用的是Surfer11.0.642的32位版本,已经将这一部分汉化。抱歉,其他版本还没有汉化,尤其是Surfer12的版本,原先辛苦做的汉化工作全部因为它的bug打了水漂。和我使用相同版本的Surfer的网友,可以下载这个压缩文件,解压后替换掉原来的两个文件即可。用其他版本的朋友请不要做这事,搞乱了程序我可不管。
Surfer11更新.rar (1.95 MB, 下载次数: 7)

评分

参与人数 2威望 +4 金钱 +20 贡献 +12 体力 +50 收起 理由
kongfeng0824 + 10 + 8 + 50 很给力!
陌小默的 + 4 + 10 + 4 赞一个!

查看全部评分

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

新浪微博达人勋

发表于 2015-6-7 00:13:39 | 显示全部楼层
老师这篇文章很有启发性
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-6-7 08:23:12 | 显示全部楼层
程老师给力啊
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-6-7 08:56:57 | 显示全部楼层
估测面降雨量有极大地现实意义,程老师V5!,为程老师点赞!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-6-7 09:49:02 | 显示全部楼层
基本上看明白了,程老师辛苦,费那么大劲就是为了让大家学会,我由衷的感动。
关于知识,求精度永无止境,(100square-85)×100/10000,是个例子,但也说明了通过鼠标点击求算的困难性,这之中应该有个批量处理的过程,我可到简单操作,因为surfer每一步都可以用程序实现,但是目前只能一般性的调用,还没有接触到您那深入。同时插值本身也是个大概值...,这里的这个度确实不好掌握。但是至少要做到在插值的基础上以毫米为单位的等值线一个一个的计算才有说服力。
最后,再次感谢程老师,祝福好人一生幸福,平安!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-6-7 12:27:13 | 显示全部楼层
谢谢程老师     
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-6-7 13:32:12 | 显示全部楼层
程老师又出教学贴了,太无私了,感谢!!可惜我的surfer是12版的
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-6-8 08:53:00 | 显示全部楼层
下载pdf文档,慢慢品味。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-6-9 22:39:10 | 显示全部楼层
程老师真是好人,之前的教程从头到尾认真学习了,现在正好需要进行面积计算就看到此贴,真是太激动了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-6-10 08:39:23 | 显示全部楼层
做得好,画的挺好看的
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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