Acrgis在GIS空间分析方面具有独特的优势,因此它已成为我目前研究的主要工具。然而,在通信过程中,由于处理数据和绘图的主要工具不同,我们经常会遇到数据格式转换的问题。本文将介绍如何使用GMT4命令实现地形数据格式arc ASCII与GMT可识别NetCDF之间的转换。 我的专栏《海洋地球科学开放数据库》列出了一些实用数据库。如下图所示,我们首先从gmrt选择并下载一个高精度arc ASCII多波束文件mar.asc。 然后使用xyz2grd 命令将Arc ASCII格式的数据直接转换为GMT可识别的netCDF格式: 如果是从JAMSTEC等数据库下载的ASCII文件,个人习惯使用nearneighbor插值: 输出的MAR.grd文件即为GMT netCDF格式。 将GMT NetCDF格式的数据转换为arc ASCII格式时,首先对NetCDF文件重新采样,使经度和纬度的采样率完全相同,然后使用grd2xyz命令将其转换为arc ASCII格式: 这里需要注意的是,重采样命令中的-I选项必须加上"=",使输出的MAR.nc文件在经、纬度方向的采样率完全相同,从而可以用于grd2xyz命令转换数据格式。 最终输出的MAR_G.asc文件即为Arc ASCII格式。使用ArcMAP计算一下hillshade,就可以看到下图的高精度海底地貌啦: 总结容易踩的坑 , r6 r6 ?3 `. l, B: s) Q$ X+ a$ Y
最后声明:本文涉及的所有GMT代码均为GMT 4版本,不适用于GMT更高版本! |