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

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

• 将Arc ASCII转换为GMT netCDF格式

  
  ! Z9 ~% S7 w& S8 \

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

( o/ E; P9 G+ A# f' F1 ^2 }                               
登录/注册后可看大图

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

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

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

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

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

• 将GMT netCDF转换为Arc ASCII格式

  
  

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

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

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

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

                               
登录/注册后可看大图

• 总结容易踩的坑

  
  , N1 |6 L0 x# n

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

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