收藏本站 劰载中...网站公告 | 吾爱海洋论坛交流QQ群:835383472

水深之处—地形数据格式Arc ASCII与GMT netCDF之间的转换 - 海洋测绘数据转化

[复制链接]
9 A" z% D8 y3 S

AcrGIS在地理信息系统空间分析有着得天独厚的优势,因此成为了我研究目前课题的主要工具。但是在交流过程中,由于大家处理数据和画图的主力工具不同,往往会遇到数据格式转换的问题。本文将介绍如何使用GMT 4的命令实现地形数据格式Arc ASCII与GMT可识别的netCDF之间的转换。

将Arc ASCII转换为GMT netCDF格式

我的专栏文章海洋地球科学开放数据库列举了一些实用的数据库。如下图,我们首先从GMRT选取下载一段高精度的Arc ASCII格式的多波束文件MAR.asc。

& p! @: G; z5 ?+ m
) v: V: p8 U6 {

然后使用xyz2grd 命令将Arc ASCII格式的数据直接转换为GMT可识别的netCDF格式:

* S3 \! \' i4 Q9 S& u" S; U
xyz2grd MAR.asc -GMAR.grd -E -V
6 W4 ^% w1 H0 l% ]/ d) e

如果是从JAMSTEC等数据库下载的ASCII文件,个人习惯使用nearneighbor插值:

5 t o# ^9 K2 d6 T' {2 N
blockmean MAR.dat -R-45.6/-44.2/25.3/27 -I0.001 > MAR.xyz: g! b- c8 z' v0 |7 L, Z3 W) } nearneighbor MAR.xyz -V -R -I -S0.002 -N2 -GMAR.grd
. n: |4 \* m# Z4 \

输出的MAR.grd文件即为GMT netCDF格式。

将GMT netCDF转换为Arc ASCII格式

将GMT netCDF格式的数据转换为Arc ASCII格式时,首先需要对netCDF文件重采样,使经、纬度方向的采样率完全相同,然后再使用grd2xyz命令将其转换为Arc ASCII格式:

7 A! x- m$ k; Y* W
grdsample MAR.grd -I0.001= -GMAR.nc6 G$ q1 P2 W0 \* }8 K {7 p grd2xyz MAR.nc -E > MAR_G.asc
# p! \, N; J, X

这里需要注意的是,重采样命令中的-I选项必须加上"=",使输出的MAR.nc文件在经、纬度方向的采样率完全相同,从而可以用于grd2xyz命令转换数据格式。

5 ?9 L& u. z( g& i* k

最终输出的MAR_G.asc文件即为Arc ASCII格式。使用ArcMAP计算一下hillshade,就可以看到下图的高精度海底地貌啦:

4 ~8 n: k: Z! q6 U
总结容易踩的坑

a. grd2xyz和xyz2grd命令处理涉及Arc ASCII格式数据时,必须加上"-E"选项!

, l5 g) t" d4 _# }

b. GMT netCDF转Arc ASCII时,第一步的重采样命令中的"-I"选项必须加上"="!

/ y7 a6 ~5 a+ H& [/ q& M$ \

最后声明:本文中涉及的GMT代码均为GMT 4版本,不适用于GMT更高版本哦!

P6 g N1 Z! \ 2 ?3 A' Z. u; d B8 N6 }2 g7 e; v5 z2 k, f & D2 X: E# V2 A k( U; u 3 a2 N. P# }& X0 l- N" Q+ d
回复

举报 使用道具

全部回帖
暂无回帖,快来参与回复吧
懒得打字?点击右侧快捷回复 【吾爱海洋论坛发文有奖】
您需要登录后才可以回帖 登录 | 立即注册
湘湘外婆
活跃在昨天 20:32
快速回复 返回顶部 返回列表