GMT加载shp文件

当我们需要在GMT中绘制来自外部的数据时，比如，自带的河流数据不能够满足要求，我们通常使用shp文件。

而我一般的做法是在Arcgis中加载shp数据，然后采用【要素折点转点】工具，可能还要使用【添加xy坐标】工具，最后导出shp文件，利用matlab工具读取点数据的经纬度数据。接着在gmt读取保存的txt文件即可，但是这样存在一个问题：因为我们是采用散点的方式，如果要绘制一个区域边界，那么如果采样点过于系数，那么显示的效果就不是很好，像下图所示：

那么有没有解决办法呢？最近我明白可以将shp文件先转成gmt文件，这样就可以直接绘制连续的线性元素。我这里借助了GDAL的ogr2ogr转换工具，具体下载安装见：

https://www.cnblogs.com/gisalameda/p/10409683.html

安装好后，我们可以在C盘下看见这样的一个文件夹：

进入文件夹，运行OSGeo4W.bat。需要注意的是，你需要提前把需要转换的shp文件放于此目录下，见下图中的changjiang.shp。

可输入ogr2ogr，弹出以下信息，说明配置成功。接着输入以下的命令即可。

然后运行得到转换后的gmt文件，直接在gmt中读取，得到完整的长江线状要素。