當(dāng)我們需要在GMT中繪制來自外部的數(shù)據(jù)時，比如，自帶的河流數(shù)據(jù)不能夠滿足要求，我們通常使用shp文件。

而我一般的做法是在Arcgis中加載shp數(shù)據(jù)，然后采用【要素折點轉(zhuǎn)點】工具，可能還要使用【添加xy坐標(biāo)】工具，最后導(dǎo)出shp文件，利用matlab工具讀取點數(shù)據(jù)的經(jīng)緯度數(shù)據(jù)。接著在gmt讀取保存的txt文件即可，但是這樣存在一個問題：因為我們是采用散點的方式，如果要繪制一個區(qū)域邊界，那么如果采樣點過于系數(shù)，那么顯示的效果就不是很好，像下圖所示：

那么有沒有解決辦法呢？最近我明白可以將shp文件先轉(zhuǎn)成gmt文件，這樣就可以直接繪制連續(xù)的線性元素。我這里借助了GDAL的ogr2ogr轉(zhuǎn)換工具，具體下載安裝見：

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

安裝好后，我們可以在C盤下看見這樣的一個文件夾：

進(jìn)入文件夾，運行OSGeo4W.bat。需要注意的是，你需要提前把需要轉(zhuǎn)換的shp文件放于此目錄下，見下圖中的changjiang.shp。

可輸入ogr2ogr，彈出以下信息，說明配置成功。接著輸入以下的命令即可。

然后運行得到轉(zhuǎn)換后的gmt文件，直接在gmt中讀取，得到完整的長江線狀要素。