物探论坛

 找回密码
 立即注册

QQ登录

只需一步,快速开始

查看: 1307|回复: 1

[GMT] GMT格网化方法

[复制链接]
发表于 2013-8-8 22:52:09 | 显示全部楼层 |阅读模式
GMT的很多命令中都会涉及到格网文件,在上一节中已经较为详细的介绍了格网文件的相关内容。本章将详细介绍如何生成格网文件。
在我们所获取的数据里,通常都是(x,y,z)形式,而这些数据有的已经是按规则的格网排列好的,而大多数都是不规则分布的数据形式。对于按规则格网排列好的(x,y,z)数据,只需用xyz2grd命令即可生成格网文件;对于不规则分布或者随意分布的数据,首先需要将这些数据进行格网化,生成按规则格网排列的(x,y,z)数据,然后再转化为格网数据。这里主要讨论随意分布的数据如何格网化并生成格网数据。
数据格网化,常用有两种算法:nearest neighbor griddinggridding with splines in tenson
1.       nearest neighbor gridding
该方法对应的GMT命令为nearneighbor,该命令的详细用法这里不做介绍,只讲述这种格网化算法。
nearneighbor命令将会指定数据范围和格网间隔,从而即可确定格网结点。此外,该命令还需要给出搜索半径,-Sradius选项。有了结点,有了搜索半径,即可搜索出该结点附近分布的点值了,根据这些点值,利用距离做权重,计算出结点处的值,即为格网点的值了。该算法要求在每个结点的搜索范围内只要有一个点,如果没有点,则该结点值被设为NaN值。
例:nearneighbor –R245/255/20/30 –I5m –S40k –Gship.nc –V ship.xyz
注:当数据较密时,适合用该方法。


2.       gridding with spline in tension
利用全局的数据进行格网化。具体说来就是,将所有的数据投影到一个surface上,surface上对应的格网点值就是格网值。因此,相比这两种方法,nearneighbor是利用局部范围的数据进行格网化,而该方法是用全局的数据进行格网化。
在利用全局的数据进行格网化之前,为了消除aliasing效应,需要用blockmeanblockmedianblockmode命令对数据进行预处理。blockmean命令适用于较为平滑的数据;blockmedian适用于起伏较大的数据。预处理后的数据,即可用surface命令进行格网化并生成格网数据了。
例:blockmean –R245/255/20/30 –I5m –V ship.xyz > ship_5m.xyz
surface ship_5m.xyz –R245/255/20/30 –I5m –Gship.nc –V
注:相同的数据,surfacenearneighbor命令格网化的结果差别较大。
此外,GMT还提供了两个格网数据提取和信息查询的命令,分别是grdrastergrdinfo,这里就不在详细介绍了。
回复

使用道具 举报

发表于 2014-1-27 15:17:52 | 显示全部楼层
请问gmt是开源的吗?
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|Archiver|手机版|小黑屋|物探论坛 ( 鄂ICP备12002012号 微信号:iwutan )

GMT+8, 2022-7-6 08:40 , Processed in 0.066756 second(s), 15 queries .

Powered by Discuz! X3.4

© 2001-2017 Comsenz Inc.

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