水深之处—地形数据格式Arc ASCII与GMT netCDF之间的转换

[复制链接]

Acrgis在GIS空间分析方面具有独特的优势,因此它已成为我目前研究的主要工具。然而,在通信过程中,由于处理数据和绘图的主要工具不同,我们经常会遇到数据格式转换的问题。本文将介绍如何使用GMT4命令实现地形数据格式arc ASCII与GMT可识别NetCDF之间的转换。

  • 将Arc ASCII转换为GMT netCDF格式


    6 s2 D# ?- g4 a; \: Y* [5 s

我的专栏《海洋地球科学开放数据库》列出了一些实用数据库。如下图所示,我们首先从gmrt选择并下载一个高精度arc ASCII多波束文件mar.asc。


, x3 i; |2 I7 m' c* h                               
登录/注册后可看大图

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

游客,如果您要查看本帖隐藏内容请回复

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

游客,如果您要查看本帖隐藏内容请回复

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

  • 将GMT netCDF转换为Arc ASCII格式


    9 c( I4 }: F. f# t8 F3 a

将GMT NetCDF格式的数据转换为arc ASCII格式时,首先对NetCDF文件重新采样,使经度和纬度的采样率完全相同,然后使用grd2xyz命令将其转换为arc ASCII格式:

游客,如果您要查看本帖隐藏内容请回复

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

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


) M# e6 K# Z+ O- w& `                               
登录/注册后可看大图

  • 总结容易踩的坑

    , r6 r6 ?3 `. l, B: s) Q$ X+ a$ Y

游客,如果您要查看本帖隐藏内容请回复

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

回复

举报 使用道具

相关帖子

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