GMT 中文社区

GMT脚本执行利器——VS Code

Wed, 19 Feb 2020 00:00:00 +0000

Visual Studio Code，简称vscode，是微软推出的一款免费、轻量且功能强大的编辑器。vscode除一键执行、语法高亮、分屏显示外，还可实现bat/bash脚本的顺滑切换运行、jpg/gif/pdf等成图结果实时预览、脚本报错信息显示，非常适合全平台下执行GMT脚本。

注：vscode是跨平台软件，支持linux/Mac/Windows平台。本文以windows为例介绍其使用方法。

vscode安装

vscode下载地址为：官网

安装中建议勾选下图中两项，目的是添加右键打开的快捷方式。

安装完成之后有三种打开vscode的方式：

开始菜单中，直接找到Visual Studio Code应用程序打开。
直接在对应脚本或工作文件夹点击右键，选择“通过Code打开（Open with Code）”。
Ctrl+R打开运行窗口，输入code (前提是重启后，vscode已加入PATH路径)。

插件安装

vscode通过不同功能的插件实现代码的开发运行、图片文件的预览等等一系列功能。执行GMT脚本命令，推荐安装的插件包括：

Code Runner 支持代码一键运行
vscode-pdf pdf预览插件

插件安装方式：打开Extensions扩展插件菜单，搜索对应插件，点击安装。安装成功之后，Install图标将变为齿轮，可点击该齿轮，修改插件属性。已安装启用的插件都可在Extension中查看设置。

在 Extensions 菜单中，点击 Code Runner 插件的齿轮图标，弹出菜单中选择“Extension Settings”，勾选以下三个选项。不勾选这三个选项可能会导致无法正确执行 GMT 脚本。

vscode执行GMT脚本过程

GMT脚本建议在独立工作空间（文件夹）中执行，有以下两种方式：

在vscode中的Explore资源管理器菜单，选择“打开文件夹“：

或者，在文件夹上右键，选择“通过Code打开”：

含有中文的脚本需要设置编码为GB2312，否则会出现乱码：点击右下角的编码按钮，再选择“重新打开编码（Reopen with Encoding）”，选择中文简体GB2312编码。

执行脚本可选择脚本页面右上方三角符号（Alt+F1），脚本运行信息显示在下方Terminal终端窗体。若想同时预览生成结果，可增加编辑窗口（Split Editor），两个窗口分别显示脚本和运行结果。

若本机安装了Bash终端，vscode也可通过手工切换终端执行 bash 脚本。

如果你运行的是 bash 脚本，则必须在 gmt begin 之前添加一行 export GMT_SESSION_NAME=$$ ，否则脚本运行将出现错误。

详情请参阅GMT中文手册的begin模块说明页面

为提高编码效率，推荐使用快捷键进行代码注释、代码运行、窗格切换等操作。键盘单击F1，调出vscode的命令面板，输入”shortcut“检索出快捷键编辑命令：

我们可以查询或自定义众多快捷键，包括但不限于行注释（Add Line Comment/Remove Line Comment）、窗口分割（View: Split Editor Down/Right/Up/Left）、窗口关闭（View：Close Editor）等。

2021版1:100万全国标准基础地理数据！你值得拥有

Sun, 15 Apr 2018 00:00:00 +0000

本文介绍国家标准地理数据的获取以及处理为 GMT 支持格式的方法。

具体包括：

国家标准1:100万（非涉密）基础地理数据的申请、下载
数据集内容简介
数据集要素合并、提取及格式转换

本文仅供个人科研作图参考。若使用者编制的地图需向社会公开的，请务必遵守《地图管理条例》有关规定依法履行地图审核程序。

数据概况

国家基础地理信息中心于2017年11月份开始免费向公众提供全国1:100万非涉密的全国基础地理数据库。目前数据最新为2021版。该数据库：

共77幅DLG数据集，包括全国陆地范围以及台湾岛、海南岛、钓鱼岛、南海诸岛在内的主要岛屿及其临近海域。
包括水系、居民地及设施、公路铁路、行政边界、地名等5个数据集12个数据层。

数据更多介绍可访问全国地理信息资源目录服务系统。

数据获取

进入全国地理信息资源目录服务系统，点右上角“注册”。已注册用户直接登陆。
在“数据下载”栏目中，点击进入“1：100万全国基础地理信息数据库”。
该数据集下有77个图幅产品。可选择一两个图幅，或全部选中，加入成果车。
结算，填写订单中的个人信息，验证手机号，下载数据。 **注：**若未提示验证手机，换个firefox浏览器试试。

数据集内容简介

该数据集包括5类12个要素集，如下图所示：

全国基础地理数据库数据集

每个要素集内又包括不同等级的自然地理要素，以BOUL要素集（图层）为例：

BOUL 中国及邻区各级行政界限（线）图层

该数据集包括各级行政区界线，以GB码进行区分：

行政区界类型	GB码
海岸线	250200
水系交汇处	260100
外国地区国界线	610200
国界线(已定)	620201
国界线(未定)	620202
省级界线	630200
特别行政区界	630400
地、市、州级行政区界	640200
县级行政区界	650200
特殊地区界限	670202

其他各图层内不同要素的分类编码可参考基础地理信息要素分类与代码（GB/T 13923-2006）。

此外，常用的还有 AGNP 各级行政地名和城乡居民地坐标（点）图层

行政地名点以政府驻地来定位，无GB码，以CLASS类别区分：

行政地名点类型	CLASS类别
国名	AA
省（直辖市、自治区、特别行政区）行政地名	AB
自治州、盟、地区行政地名	AC
地级市行政地名	AD
县级市行政地名	AE
县（自治县、旗、自治旗、地级市市辖区）级市行政地名	AF
县辖区及县级行政区域的派出机构地名	AG
街道办事处地名	AH
镇行政地名	AI
乡行政地名	AJ
建制村地名	AK
城镇区片、小区名	BA
自然村、屯、片村、村民小组名	BB
牧点、渔点、棚房名	BC
其它	BD
党政机关、党派团体名	CA
企事业单位名	CB
农、林、牧、渔场	CC

数据处理

从网站上直接下载的数据是分图幅的 geodatabase(gdb) 数据库文件，如果需要用 GMT画“全国一张图”，需要：

合并77个gdb文件中的同类要素
提取某类要素中的某个等级数据
转为gmt识别的文本格式

gdb数据库文件合并

采用本社区一贯推荐的GDAL/ORG库将分幅的gdb数据进行批量合并。首先解压所有下载的压缩包，得到以.gdb结尾的目录。目录名类似H51.gdb、I49.gdb等。然后在这些目录的上级目录执行以下脚本：

Linux系统：

#!/bin/bash
file="./final"
for i in $(ls -d *.gdb)
do
 if [ ! -d "$file" ]
      then
           echo "creating final filefolder" 
           ogr2ogr -f 'ESRI Shapefile' $file $i -lco encoding=UTF-8
      else
           echo "merging……"
      ogr2ogr -f 'ESRI Shapefile' $file $i -update -append
fi
done

Windows系统。注意必须保存为 bat 文件执行，直接复制粘贴到 cmd 中会出现报错：

REM 将A49分幅的gdb数据转为一组shp格式的文件，文件夹命名为Merged
REM 编码为UTF-8以正确显示属性字段中的中文
ogr2ogr -f "ESRI Shapefile" -lco encoding=UTF-8 Merged A49.gdb

REM 将当前目录下剩余所有的gdb数据追加Merged，同类shp文件合并。
for /f %%i in ('dir /b "*.gdb"') do (
ogr2ogr -f "ESRI Shapefile" -update -append Merged %%i
)

pause

若安装有python环境，还可选择ogrmerge.py工具进行批量的格式转换和合并，详见ogrmerge使用手册。

该节将分幅的交通、水系、边界等12类要素分别合并为全国一张图，最后转为12个shp文件。

查询GB代码与CLASS类别

使用ogrinfo可以查询shp文件内包含哪些要素。例如以下命令查询BOUL.shp中包含的GB代码：

ogrinfo -sql "SELECT distinct GB FROM BOUL" BOUL.shp

而以下命令则查询AGNP.shp中包含的CLASS类别：

ogrinfo -sql "SELECT distinct CLASS FROM AGNP" AGNP.shp

提取某类要素

由于每个数据集（一个shp文件）中存在不同等级的要素，而实际可能只使用其中一种，则仍可利用ogr2ogr工具提供的sql查询选项进行提取。

例如：

行政区划BOUL图层中GB代码620201、620202、250200分别为已定国界、未定国界和海岸线，那么提取数据并转化为GMT格式的命令为：

ogr2ogr -f GMT -where "GB=620201 or GB=620202 or GB=250200" BOUL_sub.gmt BOUL.shp

最后GMT绘图如下：

gmt begin ChinaMap
    gmt plot BOUL_sub.gmt -R70/145/3/55 -JM105/35/6.5i 
gmt end show

绘图结果：

如果想要提取行政地名点AGNP图层中的建制村信息（对应CLASS为AK），则使用：

ogr2ogr -f GMT -where "CLASS='AK'" AGNP.gmt AGNP.shp

注意因为CLASS类别为字符串，所以此处AK需要加上引号。GB代码为整型数因而不需要加引号。

最后转换成的GMT格式数据如下：

...
# @NCLASS|NAME|PINYIN|GNID|XZNAME
# @Tstring|string|string|string|string
# FEATURE_DATA
# @DAK|岙里王村|Aoliwangcun|330226114037|胡陈乡
121.68517089 29.388652295
# @DAK|坝头村|Batoucun|331022000000|亭旁镇
121.310746425 29.0109724900001
# @DAK|白岩村|Baiyancun|331082111028|尤溪镇
120.971417165 28.7226002200001
# @DAK|柏树下村|Baishuxiacun|331082116042|邵家渡街道
121.257111875 28.8522451100001
...

以 Windows 系统为例，首先将以上数据的编码转换为GB2312，再使用以下代码绘制：

gmt begin AGNP jpg
  gmt set PS_CHAR_ENCODING Standard+

  gmt basemap -R120/121/28/29 -JM10c -Baf
  gmt grdimage @earth_relief_15s -Cgeo
  gmt plot -Sc0.2c -Gred -W0p,black AGNP.gmt
  chcp 936
  gmt convert AGNP.gmt -a2="NAME" | gmt text -F+f10p,40+jTL -D0.15c/0.3c -Gwhite
gmt end

绘图结果：

GMT5在Windows下提示无法创建用户目录的解决办法

Fri, 27 Oct 2017 00:00:00 +0000

GMT5的某些版本（<5.4.4）在某些特定的Windows下会出现如下警告信息：

gmt: Unable to create GMT User directory : /Users/用户名/.gmt
gmt: Auto-downloading of earth_relief_##m|s.grd files has been disabled.

这是一个Windows特有，但又不是所有人都会遇到的一个问题。出现该警告信息的本质原因在于，当前系统没有定义环境变量 HOME，进而导致GMT无法知道当前计算机的主目录的具体位置。

虽然大部分情况下并不影响程序的正常运行，但总是出现这些警告信息真的有点难受。解决办法是：

1、打开 “我的电脑”->“属性”->“高级系统设置”->“环境变量”
2、新建 “系统变量”
     变量名： HOME
     变量值： C:\Users\用户名

重新运行GMT检查问题是否得到解决，可能需要重启计算机。

经测试，HOME 变量的值可以取任意盘符中任何已存在的文件夹，正常情况下，执行脚本后会生成 %HOME%/.gmt/cache 空文件夹。

GMT 中文社区第三次线下聚餐

Tue, 17 Oct 2017 00:00:00 +0000

金秋十月，秋高气爽，雁过留声。GMT 中文社区在北京成功举办了第三次线下聚餐活动。

2017 年 CGU 年会届时举行，全国同行汇集首都北京国际会议中心，我们借此机会齐聚一堂。个别同学本来表示由于时间原因或距离问题无法参会，最终还是觉得机会难得克服各种困难先后赶到。最终，13 位 GMTers 进行了一次成功的“网友见面会”。

在聚会上，本来活跃在 QQ 群里的大家面对面相互认识和了解了彼此的学习经历和科研方向，除了探讨了 GMT 学习、社区发展，还聊到了地震预报预警、防震减灾和地震学家如何对公众进行地震科普等问题。不仅如此，本次聚会大家还进行了一些简单的跨学科交流，桌上言谈之间对于自然科学研究及探索的热情和勇气溢于言表。

希望通过本次聚会我们能更加熟悉和了解该领域各个研究方向的同学，彼此保持更密切的联络，在科研道路上相互交流、积极互助。

科研道路上，知己天涯若比邻，期待下次再聚~

从左至右依次为：章浩东（中科院地质与地球物理所）、王志伟（地震局地质研究所）、王少卿（地震局地球物理所）、张琰（地震局地球物理所）、张周（广州地化所）、徐荟（地震局地球物理所）、王亮（成都理工）、杨晓瑜（地震局地球物理所）、李佳威（地震局地球物理所+北大）、神秘人、田冬冬（中科大）、任昱（南海所）、方洪健（中科大）、胡少乾（中科大）

利用 Google Maps 获取超高精度高程数据

Thu, 14 Sep 2017 00:00:00 +0000

在绘制超小区域的地形图时，常需要超高精度的高程数据。

目前已知的公开高程数据中，ASTER GDEM的精度最高，为1弧秒，即约30米一个数据点。可以通过 Google Maps Elevation API 获取米级精度的数据。

数据简介

Google Maps Elevation API 提供地球表面所有位置的海拔高度数据，包括海床上的深度位置（此类位置返回负值）。

官方文档：

获取密钥 API KEY

想要从Google Maps中获取地形数据，首先需要一个密钥。可以从 get API key 页面开始申请秘钥（需要提供信用卡等信息）。

限额说明：

每个用户最多可能只能申请10个免费密钥
每个密钥每天2500次免费请求
每次请求最多512个位置
每秒最多50次请求

使用说明

获取单点高程

通过简单的网络请求即可获取单点的高程：

https://maps.googleapis.com/maps/api/elevation/json?locations=39.7391536,-104.9847034&key=YOUR_API_KEY

其中：

locations=39.7391536,-104.9847034 指定了要获取高程的经纬度，纬度和经度之间用逗号分隔；
key=YOU_API_KEY 指定了用户密钥，用户需要将自己申请的密钥替换 YOU_API_KEY

将上面的链接复制到浏览器中，即可得到返回值为::

{
   "results" : [
        {
            "elevation" : 1608.637939453125,
            "location" : {
                "lat" : 39.73915360,
                "lng" : -104.98470340
            },
            "resolution" : 4.771975994110107
        }
   ],
   "status" : "OK"
}

从中可以看到这一点的海拔高度是 1608.637939453125 米（小数点后几位可能都没有意义）。 resolution给出了该点高程的分辨率值，即内插海拔高度的数据点之间的最大距离，以米为单位。

获取多点高程

也可以在一次网络请求中获取多个点的高程，多个位置之间用 | 分隔即可。例如：

https://maps.googleapis.com/maps/api/elevation/json?locations=39.7391536,-104.9847034|36.455556,-116.866667&key=YOUR_API_KEY

但一次请求中的位置不能超过512个。

批量获取高程

Google Maps 提供了多种客户端以实现批量获取高程。此处以Python客户端为例。

安装 Google Map API 的 Python模块:
```
 $ pip install -U googlemaps
```
获取密钥
参考如下脚本批量获取高程数据

#!/usr/bin/env python
# -*- coding: utf-8 -*-

import googlemaps
import numpy as np


def request_elevation(xmin, xmax, xinc, ymin, ymax, yinc, maxnum=512):
    """Request elevation and return them with a list

    Parameter
    =========
    xmin : float
        minimum longitude, unit is degree
    xmax : float
        maximum longitude, unit is degree
    xinc : float
        increment of longitude

    ymin : float
        minimum latitude, unit is degree
    ymax : float
        maximum latitude, unit is degree
    yinc : float
        increment of latitude
    maxnum : int
        point number of one request, maximum value is 512
    """

    xlist = np.arange(xmin, xmax, xinc)
    ylist = np.arange(ymin, ymax, yinc)
    points = [(y, x) for y in ylist for x in xlist]

    requestnum = int(np.floor(len(points) / maxnum) + 1)
    npoints = np.array_split(points, requestnum)
    if requestnum > 2500:
        print("Error: Request number exceed! Change grid or region.")

    return [gmaps.elevation(locations=loc.tolist()) for loc in npoints]


# 设置要下载高程数据的经度范围、纬度范围以及网格间隔
xmin, xmax, xinc = 120.000, 129.005, 0.00005
ymin, ymax, yinc = 30.000, 30.005, 0.00005

gmaps = googlemaps.Client(key='YOU_API_KEY')

# Obtain all elevations
# 请设置 maxnum 确保每一次请求的数据点数不超过 512 个以控制请求字符串长度
elevations = request_elevation(xmin, xmax, xinc, ymin, ymax, yinc, maxnum=512)

# output data with format `longitude   latitude elevation resolution`
with open("Elevations.dat", "w") as f:
    for loclist in elevations:
        for subloc in loclist:
            outputstr = "{:.5f}  {:.5f}  {:.5f}  {:.5f}\n".format(
                subloc['location']['lng'],
                subloc['location']['lat'],
                subloc['elevation'],
                subloc['resolution'])
            f.writelines(outputstr)

GMT 中文社区 QQ 群入群须知

Wed, 23 Aug 2017 00:00:00 +0000

欢迎加入“地学GMT中文社区群”：

一群：218905582
二群：791856541

GMT 中文社区群是一个相对活跃且学术氛围浓厚的群。本文是 GMT 中文社区群的入群须知，也是 GMT 中文社区的提问规范。希望所有群友可以阅读并遵守以下规范来营造更好的群环境。

群内可以做的事

GMT 相关问题
其他地学软件相关问题
地学科研相关问题
求助文献
发科研相关的培训、会议信息
介绍并推广自己的最新研究成果

群内禁止做的事：

发广告
签到、灌水等无意义行为
人身攻击
讨论政治问题

安装类问题

如果你在安装 GMT 的过程中遇到问题，请参照下面的模板提问。

我在安装 GMT 的时候遇到了问题，请问该如何解决？谢谢。

操作系统是: Windows XP | Windows 7 | Ubuntu | CentOS | macOS …
要安装的GMT版本是: GMT 4.5.16 | GMT 5.4.2 …
安装GMT参考的文档是: 此处可贴文档的网址
安装过程中显示的错误信息为: 此处可贴图

思路类问题

想要实现某一类图，但不知GMT能否实现或没有实现的思路。

如果别人文章中有自己想要实现的效果图就直接贴图，并指出自己在实现图中哪一部分时没有思路。如果找不到已有的效果图，可以考虑手绘示意图，或精准描述自己想要的效果，以免引起误会，造成不必要的时间浪费。

此处先贴一张图

请问图中的XXXXX能否用GMT实现？该如何实现？谢谢。

请问我想用GMT实现XXX效果。此处详细介绍自己想要的效果 。不知GMT中该如何实现？谢谢。

如果遇到自己画出来的图与自己预期的图不一致，请参考如下模板。

当前效果：截图
预期效果：预期效果截图，或对预期效果的精确描述
我的代码是：给出自己使用的代码，可以截图可以贴代码。最好是给出最小示例而不是一堆代码。
GMT版本是: GMT 5.4.2

脚本报错类

遇到绘图脚本或数据处理脚本出错时，请参考如下模板。

我想要实现的是 ：…
我的命令/脚本是: 可截图或贴脚本
终端报错是: 终端报错截图

文献求助

需要他人帮忙下载文献时，请先使用 sci-hub 工具自行尝试下载（sci-hub教程）。

对于中文文献或 sci-hub 无法下载的文献，方可在群中寻找帮助。求助文献时，请先明确告知“使用 sci-hub无法下载该文章”并提供文章的官方链接。若找不到官方链接，则提供文章的标题、作者、期刊、年份等信息。

求助文献示例如下：

文献求助！我已经尝试用 sci-hub 下载但依然下载不到。文章链接为: XXXXXXX。谢谢。

提问的智慧

时间对任何人都是一笔宝贵的财富。没有人有义务回答你的问题。在期望别人回答你的问题之前先要问一问自己针对这个问题做了哪些努力。

提问前需要做的事情

大多数的问题，文档都可以为你提供标准的解答。因此在你提问之前，你最好先去阅读一下相关的文档，至少确保你阅读过一份完整的入门文档。如果你能够在你提问的同时表明自己已经阅读过文档，但是依旧留有困惑，潜在的回答者会更加愿意为你解答。

推荐的一些文档：

提问技巧

精确描述自己做了什么，遇到了什么问题，要提供足够的信息给潜在的回答者，尽量避免回答者多次要求提问者补充额外的信息的情况。
提供当前所使用的代码。贴完图之后如果自行修正了代码中的某些错误，一定要更新图。
提供截图以及文本。有些时候需要准确看到代码的样子，此时截图更合适；有些时候需要复制代码在自己电脑上跑一遍，此时提供文本更合适。
提问时说一说自己为了解决这个问题做了哪些事情，没人喜欢“伸手党”。
时间对任何人都是一笔宝贵的财富，请心怀感恩之心。
当你解决了某个问题（不管是别人帮助还是自己茅塞顿开），请不要吝啬告诉大家你的解决方案
当别人正在讨论某个问题时，请不要中间插入自己的问题，待别人的问题告一段落之后再提问。
除非必要，不要在首次发问的时候@某个人，这样给所有人的感觉是你在向这个人提问，从而失去了更多的潜在回答者。
有问题就直接问，不必管有没有人用过。很多时候，你在A软件中遇到的问题可能跟A软件完全没有关系，完全没接触过A软件的人可能也可以帮你解决该问题。

投稿须知

Thu, 17 Aug 2017 00:00:00 +0000

GMT 中文社区欢迎读者投稿，稿件主要包括社区博文和绘图实例两大类。

稿件类型

社区博文

社区博文是与 GMT 官方及社区活动有关的简短文章，包括但不限于：

GMT 官方新闻（比如新版本发布）
GMT 社区新闻（比如聚餐与培训）
GMT 使用技巧等

绘图实例

绘图实例是指使用 GMT 绘制的具有典型性的图件。目前社区维护的绘图实例位于 GMT 中文手册：绘图实例中。

请参考GMT 中文手册：贡献指南，直接向中文手册提交绘图实例。

投稿方式

可以通过如下几种方式投稿：

直接在社区主页或 GMT中文手册的 GitHub 仓库中提交 Issue 或 Pull Request。（不熟悉 Pull Request 的读者可以参考 Pull Request 流程
将稿件发送至 admin@gmt-china.org

作者权益

请确保你有将稿件投稿到社区的权利，特别是在图件修改自他人图件的情况下。你所拥有的权益是稿件公开后，你有署名权（如果维护者对你的稿件有大改，则维护者可能会在你的名字后面署名）。接受的稿件在公开后，任何人对你的稿件都可以以任何媒介以任何形式进行复制、发行、修改、转换或以其为基础进行再创作用于任何可能的用途，包括商业目的。别人在使用你的作品之后，应当提及你和 GMT 中文社区，但社区无法保证第三方会尊重你的权益。

使用sci-hub下载文献

Sat, 12 Aug 2017 00:00:00 +0000

科研人员经常要看一些全文，但由于单位不可能花钱购买所有数据库的下载权限，因而会经常需要需要付费下载文章的情况。本文将介绍文献下载的好工具 sci-hub，可以下载到绝大多数外文期刊的文章甚至书籍。

sci-hub 网站

sci-hub 是一个在线网站，可以下载众多本该付费下载的文章，因而本质上是非法的，用户请自行判断。原始域名已经被封禁，用户可以访问如下网站获取可用的网址：

https://whereisscihub.now.sh/

使用方法

下面将用多种方法演示如何下载地球科学引用率非常高的一篇文章 Preliminary reference Earth model。

要下载一篇文章，需要知道如下两个信息中的至少一个：

文献的官方网址，比如 https://www.sciencedirect.com/science/article/pii/0031920181900467, 可以通过 Google Scholar 等网站搜索得到
文献的DOI，比如 10.1016/0031-9201(81)90046-7，可以从文章的官方网站上找到，现在大多数文章的参考文献也要求给出DOI信息

DOI 是 Digital Object Identifier 的缩写，也叫唯一标识符，是国际通用、全球唯一、终身不变的数字资源标识符。正如我们每人都有自己独一无二的身份证号一样， DOI代码具有唯一性，一旦分配，终身不变。已知DOI，将 DOI 加到 https://dx.doi.org/ 后（例如 https://dx.doi.org/10.1016/0031-9201(81)90046-7 ）即可跳转到该文章的真实网址。

到 sci-hub 网站下载

访问 sci-hub 网址，比如 https://sci.hubg.org/ ，如下图所示：

将文献网址或DOI粘贴到搜索框中，按下确认键或点击搜索框右边的open按钮即可。

直接利用网址下载全文

打开文章对应的网址，在域名的后面加上 sci-hub 的域名，比如 sci-hub.org 即可。

例如，这篇文章的网址是 https://www.sciencedirect.com/science/article/pii/0031920181900467 ，在域名的后面加上 sci-hub 的域名，变成： https://www.sciencedirect.com.sci-hub.org/science/article/pii/0031920181900467 即可直接下载文章全文。

直接利用DOI下载全文

已知文章DOI的前提下，直接在 sci-hub 域名后面加上 DOI 即可，比如直接访问 https://sci-hub.org/10.1016/0031-9201(81)90046-7 即可下载文章全文。

注意事项

当文献下载过频繁会要求用户输入验证码，如下图：
DOI是如 10.1016/0031-9201(81)90046-7 这样的字符串，前面不要带 doi: 等其他字符。
sci-hub 不支持下载国内期刊的文章

GMT中文社区第二次线下聚餐

Tue, 27 Sep 2016 00:00:00 +0000

2016年9月21日，借GMT中文社区首次线下培训契机，在社区创始人田冬冬博士和GMT群管理员刘珠妹老师的号召下，GMT中文社区在湖北武汉组织了第二次线下聚餐（GMT User Dinner，简称GUD），四位群管理员和来自河北、福建、安徽、甘肃、湖北等地的共16位社区成员参加了此次聚会。

聚会伊始，杨萍带来了兰州似火的热情，为大家奉上美味糖果牦牛肉，大红色的包装纸被单斌老师认成喜糖。红色代表着喜庆，象征着今天大家有缘聚在一起，也是一种皆大欢喜的团圆。随后大家自我介绍了自己的群内ID和真实姓名，社区成员对彼此的学习情况、工作经历、研究领域等有了初步认识。稍后，社区成员商讨了GMT社区网站的美观设计、腾讯网盘容量已扩容性问题、相同问题反复提问和解答等。针对上述内容，田博士希望将来可以重新设计一个漂亮的社区logo，建议把GMT社区内容重新移植回QQ群文件中，并希冀整理群问题的工作可以尽快开展，因此再次号召更多志愿者加入到GMT社区的管理和维护工作。稍后，大家谈到了亘古不变的话题“虎哥究竟是男是女”，众人仁者见仁智者见智，各持己见，瞬间聚会气氛达到了高潮，可惜的是虎哥本人并不在现场！

随后田博士在师弟和美女的陪同下，参观了位于东湖之滨、珞珈山下、拥有百年历史的全国重点高校—国立武汉大学。走出武大校园后，随后乘船游览了国家5A级旅游景区—东湖。此次聚会活动虎哥没有参加，但虎哥热情依在，虎哥看完聚会照片，和大家调侃了“桂花山药” 美味的一盘菜，并押宝5毛，最后在虎哥强烈邀约下，筹备了第三次线下聚餐活动。

GMT社区线下聚会在轻松愉快的氛围中结束，此次聚会谈论话题广泛，但又不失主题性，为GMT社区成员之间深入交流提供了一个很好的平台，促进了GMT社区的稳健成长。遥祝GMT中文社区像东湖一样，可以永葆原创魅力，又像武汉大学一样，亦可承载一路的历史、厚重、辉煌！ “相见时难别亦难”，美好的聚会终有一别，这次聚会虽然短暂，但短暂的聚会让我们会加倍的去珍惜，希望下一次聚会再团聚！

参会人员从左往右依次为：章婷（湖北省地震局）、刘珠妹（中国地震局地震研究所）、王丽霞（中国地震局兰州地震研究所）、杨萍（中国地震局兰州地震研究所）、郑海刚（安徽省地震局）、谢小晋（中国地质大学-武汉）、贾媛（山东省地震局）、杨锦玲（福建省地震局）、田冬冬（中国科学技术大学）、肖晓（武汉大学）、赵志远（河北省地震局）、刘旭耀（中国地质大学-武汉）、罗娜（河北省地震局）、邓明莉（中国科学院测量与地球物理研究所）、刘成利（中国地质大学-武汉）、单斌（中国地质大学-武汉）

GMT中文社区首次线下培训

Sun, 25 Sep 2016 00:00:00 +0000

2016年9月21至22日，GMT中文社区协同中国地震局地震研究所和中国地震局地壳形变学科技术管理组举行了首次GMT线下培训活动。

此次培训邀请了田冬冬（中国科学技术大学在读博士，社区 ID：知名不具）、单斌（中国地质大学特任教授，社区 ID：Bins）和刘珠妹（中国地震局地震研究所，社区 ID：忆尤）三位老师在湖北省地震局21楼报告厅做了《GMT入门、基础及进阶》、《GMT应用与科技论文写作》、《地学常见GMT图件模板讲解》专题培训讲座，来自全国地震系统的学员和慕名人士共计120余人参加了此次培训。

此次培训考虑到培训人员GMT操作能力差异巨大，为了尽可能地照顾到所有培训人员的需求，经过一个月精心调研准备，从内容、思路、绘图技巧等方面综合考虑，确定了此次GMT专题讲座的最终主题。

21日上午田冬冬博士首先介绍了GMT的概念、开发团队、GMT的版本等，随后对GMT基础知识、配置文件、常用命令进行了详细讲解，最后总结了GMT学习方法及学习资源，对GMT中文社区创建初衷、社区功能、未来发展等做了阐述，并呼吁各位仁人志士加入到GMT社区的设计和日常维护中。

单斌教授对GMT做了简单介绍，随后深入浅出的介绍了GMT常用选项、常用命令、shell脚本编写等，并对GMT在科技论文中的应用做了具体实例讲解，提供了科技论文写作技巧。最后具神指出了GMT5和GMT4的语法命令区别，对学习兴趣、专业优势等和大家做了深入分享交流。

22日下午刘珠妹老师对地学常见GMT图件模板进行了讲解，刘老师首先讲授了怎样读懂脚本的使用平台、脚本版本、脚本结构等内容，然后讲授了编写绘图脚本的6个基本步骤，介绍了自由空气异常图、布格异常图等几种常用绘图模板，讨论了中文支持、脚本批量运行、二次开发等GMT绘图技巧，刘老师强调在学习过程中，应对照样例修改参数，具体问题具体分析，带着问题去找答案。

此次GMT专题培训授课内容全面、专题性比较强、课程安排比较合理，通过此次培训，推广了GMT的使用范围，提高了GMT软件在日常工作中实际操作应用能力，促进了GMT群和GMT中文社区的未来更快更好发展，加强了GMT群成员的交流，此次培训取得了预期效果，期待下一次培训的召开。

GMT中文社区首次线下聚餐

Sat, 24 Sep 2016 00:00:00 +0000

2016年6月28日，GMT中文社区在北京举行了首次GMT用户线下聚餐（GMT User Dinner，简称GUD）。

6月27日至29日，很多GMT用户参加了在中国科学院地质与地球物理研究所举办的 “首届中国区域地震学参考模型交流会”。社区创始人田冬冬提议借此机会，举行社区的第一次线下聚会。

6月28日晚，参与聚会的GMT社区用户相约前往北京北辰西路的芝麻花餐厅聚会。在聚会上，社区成员相互认识，对彼此的研究领域、GMT 学习和社区的发展进行了交流和讨论。

左至右依次为：李佳威（中国地震局地球物理研究所）、王志伟（中国地震局地质研究所）、田冬冬（中国科学技术大学）、王亮（桂林理工大学）、李君（中国地震局第二监测中心

GMT国内镜像正式建成

Fri, 08 Jul 2016 00:00:00 +0000

经过本社区积极联系，GMT中国区镜像站点已于2016年7月6日正式建立，并已加入GMT官网镜像列表中。该镜像站点依托中国科学技术大学LUG服务器，是中国唯一的一个 GMT 镜像站。 GMT国内镜像的建立，可以大大加快下载GMT源码及相关数据的速度。

GMT 源码包镜像

镜像地址：http://mirrors.ustc.edu.cn/gmt/

GMT 源码包镜像中包含了 GMT 源码软件包以及 GSHHG 和 DCW 数据压缩包。进入 GMT 镜像即可直接下载 GMT 源码包，Windows/Mac 安装包位于镜像的 bin 目录下。

GMT 远程数据镜像

镜像地址：https://mirrors.ustc.edu.cn/gmtdata/

GMT 远程数据镜像中提供了全球地形起伏等地学数据。该镜像的使用方法可以参考 https://docs.gmt-china.org/latest/dataset/usage/。

GMT中文社区正式成立

Wed, 06 Apr 2016 00:00:00 +0000

GMT的中文用户们，大家好，我是 seisman 。

GMT中文社区于今日正式成立了！

GMT中文社区是针对GMT中文用户而成立的社区。整个社区目前包含了三个主要部分：社区主站、GMT中文手册和问答社区。

在社区的第一篇博文中，我想谈谈建立这个社区的原因及一些想法。

GMT中文手册

GMT中文手册的项目始于2014年，到2016年初已经初步成形。

在2014年，网络上能够找到的几本GMT中文文档，大多基于GMT4甚至GMT3。这几本文档存在一些明显的缺陷：

排版质量一般，大多是在Word里写好然后转换成PDF的
缺乏维护，几乎都是写好发布之后就没有再更新
不完整，很多有价值的内容都没有被包含在文档中
不严谨，很多文档声称是基于GMT4，但实际上某些命令却在使用GMT3的语法
无法适配最新版的GMT，某些时候对于新用户而言是一种障碍
重复劳动：不同的作者，花费了大量的时间和精力，翻译整理出了几份相似的文档，这其中很多工作都是在重复劳动

因而，我于2014年开始写针对GMT5的中文手册。新的GMT中文手册具有如下特点：

开源：源代码用rST语言写成，并在GitHub上开源
高质量排版：文档有网页版和PDF版，保证高质量的排版（目前尚存在不足，见第5条）
完整性：尽可能覆盖GMT的方方面面，成为一个完整的参考手册（没有最完整，只有更完整，见第5条）
协同合作：依托强大的GitHub和Git，使得多人合作共同编辑同一份文档变得很简单
持续维护：不仅适配GMT的新版本，而且不断修正文档中存在的错误
严谨性：尽可能保证所有命令均通过GMT最新版本的测试

经过两年努力，文档的整体结构已经大体完成。我希望文档今后能够持续更新，不断完善，不管是接下来的GMT 5.2.2、5.3.0还是未来的GMT 6.0甚至7.0，都能随着GMT版本的升级而不断更新，但这些已非我一人之力所能完成。众人拾柴火焰高。我希望有更多的GMT中文用户加入到文档维护的队伍中来，这便是我建立GMT中文社区的目的。

问答社区

问答社区试运行了一段时间，但最终还是因为用户参与度较低而决定关闭问答社区。

GMT问答社区的前身是地学GMT学习群（群号：218905582）。QQ群属于即时通讯，通常在群中提的问题都可以快速得到答案，但与此同时，QQ群存在如下缺点：

即时性导致的随意性：由于QQ本质上是一种聊天工具，导致很多人在提问时太过随意。遇到问题时，不加思索在QQ群中提问，却不去思考如何提问，如何在提问时提供完整有效的信息，如何在提问时解释清楚自己想要实现什么和遇到了什么问题。 QQ群中的很多提问都是如此，需要回答者多次向提问者发问才能搞清楚提问者的问题。这不仅浪费了回答者的时间，也不利于提问者能力的提高。
多话题引起的错乱：经常遇到几个人同时提问，几个话题互相交叉，使得搞不清到底谁在回复谁
问答记录无法长期保存：只有用户接收到的消息会被保存到本地磁盘中，而对于新加入的成员而言，则无法看到历史聊天记录
问答记录无法搜索：仅当聊天记录保存在本地磁盘中才能被搜索，无法通过搜索引擎检索
问答记录无法被非群成员搜索：有大量的GMT用户未加入群中，这些非成员则无法获取这些优质的问答资源
Linux用户提问不便：Linux下几乎没有可用的QQ客户端，导致众多Linux下的GMT用户只能通过手机打字在群中提问，且只能通过手机拍照的方式上传屏幕截图

鉴于QQ群的以上几个缺点，萌生了建立问答社区的想法。希望这样一个问答社区可以解决QQ 群存在的问题。问答社区将努力保持高信噪比，专注于GMT的问与答。

GMT中文社区

最初是计划把中文手册以及问答社区都放在我自己的博客域名之下的，但这样做一方面不够正式，另一方面也会影响到我对于自己博客的整体定位。最终经过一番考量之后，决定单独购买域名，建立GMT中文社区，然后将文档以及问答社区作为GMT中文社区的子项目。

现在，整个GMT中文社区的基本框架已经搭建完毕。我希望它能顺利运行下去，也欢迎更多的人提意见和建议，加入到整个社区的维护中来。

最后，我要感谢王亮在整个社区的建立过程中和我交流与讨论，并给出了很多有意义的建议。我还要感谢各种开源服务，我的工作是建立在它们之上的。我也希望华语地学界有更多的人了解开源、认同开源、热爱开源并最终参与到开源中来。