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

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

[复制链接]
0 R% x" z M6 H: a

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

将Arc ASCII转换为GMT netCDF格式

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

) m, i4 X2 B, k$ m4 L7 H
; T6 {. o' E2 Y- F# P

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

% _9 T7 [5 G" `3 f0 v% z9 o `
xyz2grd MAR.asc -GMAR.grd -E -V
6 a5 J+ Q& \! C+ d# C8 {; g1 X

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

* ?& n4 d: Y) @6 A# F" ?& y. V; i
blockmean MAR.dat -R-45.6/-44.2/25.3/27 -I0.001 > MAR.xyz* m$ W6 }' B- P; c6 F nearneighbor MAR.xyz -V -R -I -S0.002 -N2 -GMAR.grd
: D0 W) @5 v& `

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

将GMT netCDF转换为Arc ASCII格式

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

2 o8 `* ], ~$ O+ W1 D$ b/ y
grdsample MAR.grd -I0.001= -GMAR.nc & E$ s7 p8 ]. J1 A* C9 t grd2xyz MAR.nc -E > MAR_G.asc
3 O+ w4 [" `- P; w- Q

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

, v) Z& N: Y6 F' l" R$ p W2 W

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

! V9 S; e+ u. Y, K& F" K! a
总结容易踩的坑

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

# h+ A7 O4 n1 l. n& S) U* g

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

9 k$ p0 g' ]+ ~1 p) w

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

8 u4 n8 a! E) z" K" N( z) U K5 l5 M2 M# S, g, f9 k/ s- H3 @ ; O0 W+ I6 l3 v. E7 Q6 K; B) H$ `- H $ o5 [7 o8 l* w; y y' j" [
回复

举报 使用道具

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