GMT绘制海底地形和剖面图

---

前言

---

使用GMT绘制海底地形和剖面图，主要参考GMT中文社区里面的案例。

一、GMT是什么？

GMT（Generic Mapping Tools）是一个用于绘制地图、地震图、海洋图等的开源软件套件，非常适合用于制作海底地形图（也称为水深图或海底地形图）。

二、使用步骤

注：文心一言给出的步骤。需要进一步学习可以到GMT中文社区去查看手册。

**1，获取数据：**首先，你需要获取海底地形数据。这些数据通常以网格（grid）的形式提供，其中包含了海底不同点的深度信息。一个常用的数据源是ETOPO5，它提供了全球范围内的海底地形数据。

**2，安装GMT：**在你的计算机上安装GMT软件。GMT支持多种操作系统，包括Windows、Linux和Mac OS X。你可以从GMT的官方网站下载并安装它。

**3，准备数据：**将你的海底地形数据转换为GMT可以识别的格式。GMT支持多种网格文件格式，但你可能需要将你的数据转换为GMT的默认格式之一（如.grd）。

**4，绘制地图：**使用GMT的命令行工具或脚本语言（如bash、Python等）来绘制海底地形图。你可以指定地图的投影方式、颜色方案、图例等。例如，你可以使用grdimage命令来绘制网格数据，并使用-C选项来指定颜色方案。

**5，添加图层：**除了海底地形图之外，你还可以添加其他图层来丰富你的地图。例如，你可以添加海岸线、国家边界、海洋特征（如洋流、海底山脉等）等。GMT提供了许多命令和选项来添加这些图层。

**6，导出地图：**完成地图绘制后，你可以将地图导出为多种格式，包括PNG、JPEG、PDF等。你可以使用psconvert命令来导出地图，并指定所需的输出格式和选项。

三、绘图

在准备好之后就可以按照GMT中文社区的例子来修改自己想要的区域和地形了

gmt begin sat_plot png
gmt set MAP_FRAME_TYPE plain
gmt set FONT 10p,Times-Roman,black
gmt grdcut @earth_relief_15s -R128/132/15/18 -Gspac_33.nc  
gmt makecpt -Cbroc -T-8000/-3000  
gmt grdimage spac_33.nc -I+a15+ne0.75 -JM15c -B --FORMAT_GEO_MAP=dddF  
gmt colorbar -Bxa500f100 -DjcR+w8c/0.3c+o1.3c/0c+m+ebf  -By+l(m) -F+gwhite+p0.1p+c0.1c
gmt basemap -TdjLT+w3c+l+o0.5c
gmt basemap -Lg131.5/15.2+c2+w50k+f+u -F+p0.1p+gazure+r 
echo 129.7166667 16.93333333 c | gmt plot  -Sc0.3c -W0.5p,white -Gred 
echo 129.4 17 WPS-1/2 | gmt text -F+f11p,1,red
gmt inset begin -DjBL+w1.5i+o0i/0i -F+gwhite+p0.1p+c0c
        gmt basemap -JG140/20/? -Rg -Bg 
        gmt makecpt -Cearth -T-8000/5000
        gmt grdimage @earth_relief_01m -I+d
        gmt coast -JG140/20/? -Rg -Bg -W0.2p,black -Ggrey
        echo 120 20 140 0 | gmt plot -Sr+s -W1p,red
        echo 150 20 Pacific |gmt text -F+f7p,0.8,white
    gmt inset end  
echo 129.7 17 > ridge.txt
echo 129.9 16.8 >> ridge.txt
gmt plot -Rspac_33.nc -W1.5p,blue ridge.txt
gmt plot -Sc0.2c -Gblue ridge.txt
gmt grdtrack ridge.txt -Gspac_33.nc -C50k/4k/5k+v -Sm+sstack.txt > table.txt
gmt plot -W0.8p table.txt
gmt plot -R-30/30/-8000/-4000 -JX15c/7.5c -W3p stack.txt -i0,1,5,6 -L+b -Glightgray -Y-h+2.8c
echo -4 -7700 | gmt plot  -Sc0.3c -W0.5p,white -Gred 
echo 0 -7750 WPS-1/2 | gmt text -F+f11p,1,red
gmt basemap  -Bxafg500+l"Distance from ridge (km)" -Byaf+l"Depth (m)" -BWSne
echo "MEDIAN STACKED PROFILE" | gmt text -Gwhite -F+cTC+f14p -Dj8p

rm -f ridge.txt table.txt stack.txt spac_33.nc
gmt end show

绘图结果