最近下载水深数据时候，发现NGDC不提供自定义区域下载，但可以下载全球的数据，这样以来，用起来还是不方便。同事给我发一个博客，是一位知乎作者，叫任昱，德国亥姆霍兹基尔海洋研究中心 海洋地球科学博士在读，以下内容转载自他的文章，https://www.zhihu.com/column/c_184222208
本文的第二部分，netcdf和arc ASCII码的转换的方法，非常的实用。
对于大多数科研工作者来说，利用开放数据库是非常重要的一个技能。本文将会列举一些我在科研工作中使用的优秀的海洋地球科学开放数据库（open access geo-databases）。
- 海底地形起伏数据
1.1 全球高精度水深数据
JAMSTEC - 记得勾选“Bathymetry (MBES)”
1.2 区域高精度水深数据
IHO Data Centre for Digital Bathymetry
Baltic Sea Bathymetry Database
AusSeabed Bathymetry and Backscatter Data, Australia
1.3 15弧秒全球地形起伏数据
GEBCO 2019 Gridded Bathymetry Data
GEBCO - Undersea Feature Names Gazetteer
2. 地震数据
2.1 地震目录
Data Centers Supporting FDSN Web Services
2.2 震源机制解
ISC Bulletin - Focal mechanism
NIED F-net Earthquake Mechanism Information
Istituto Nazionale Di Geofisica E Vulcanologia - Time Domain Moment Tensor Catalogue, Italy
2.3 地球模型
GLiM - Global Lithological Map
CRUST 1.0: Global Crustal Model at 1x1 Degrees
IRIS Earth Model Collaboration
3. 重力
4. 地磁
The Global Seafloor Fabric and Magnetic Lineation Data Base Project
Geomagnetic and Electric Field Models
5. 地热
6. 火山
Volcano Location Database Search
水深之处—地形数据格式Arc ASCII与GMT netCDF之间的转换，https://zhuanlan.zhihu.com/p/138883556
AcrGIS在地理信息系统空间分析有着得天独厚的优势，因此成为了我研究目前课题的主要工具。但是在交流过程中，由于大家处理数据和画图的主力工具不同，往往会遇到数据格式转换的问题。本文将介绍如何使用GMT 4的命令实现地形数据格式Arc ASCII与GMT可识别的netCDF之间的转换。
- 将Arc ASCII转换为GMT netCDF格式
我的专栏文章海洋地球科学开放数据库列举了一些实用的数据库。如下图，我们首先从GMRT选取下载一段高精度的Arc ASCII格式的多波束文件MAR.asc。
然后使用xyz2grd 命令将Arc ASCII格式的数据直接转换为GMT可识别的netCDF格式：
xyz2grd MAR.asc -GMAR.grd -E -V如果是从JAMSTEC等数据库下载的ASCII文件，个人习惯使用nearneighbor插值：
blockmean MAR.dat -R-45.6/-44.2/25.3/27 -I0.001 > MAR.xyz
nearneighbor MAR.xyz -V -R -I -S0.002 -N2 -GMAR.grd输出的MAR.grd文件即为GMT netCDF格式。
- 将GMT netCDF转换为Arc ASCII格式
将GMT netCDF格式的数据转换为Arc ASCII格式时，首先需要对netCDF文件重采样，使经、纬度方向的采样率完全相同，然后再使用grd2xyz命令将其转换为Arc ASCII格式：
grdsample MAR.grd -I0.001= -GMAR.nc
grd2xyz MAR.nc -E > MAR_G.asc这里需要注意的是，重采样命令中的-I选项必须加上"="，使输出的MAR.nc文件在经、纬度方向的采样率完全相同，从而可以用于grd2xyz命令转换数据格式。
最终输出的MAR_G.asc文件即为Arc ASCII格式。使用ArcMAP计算一下hillshade，就可以看到下图的高精度海底地貌啦：
- 总结容易踩的坑
a. grd2xyz和xyz2grd命令处理涉及Arc ASCII格式数据时，必须加上"-E"选项！
b. GMT netCDF转Arc ASCII时，第一步的重采样命令中的"-I"选项必须加上"="!
最后声明：本文中涉及的GMT代码均为GMT 4版本，不适用于GMT更高版本哦！