GMT绘制3D柱状地理分布图

本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。

WoS数据

根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。

机构数据：

Affiliations Record Count % of 8,505
% |6 u" p- y9 ^, t# `* JNATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747
' m  v$ W- u7 m% ECENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE  828 9.735
' }* e; R* `1 l1 _9 }7 n5 P. sINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT  707 8.313
CHINESE ACADEMY OF SCIENCES 629 7.396
UNIVERSITE DE TOULOUSE 517 6.079
UNIVERSITE TOULOUSE III PAUL SABATIER 498 5.855
NASA GODDARD SPACE FLIGHT CENTER 483 5.679
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162
) K% x) \- p; ^' q8 E3 ~1 PNASA JET PROPULSION LABORATORY  431 5.068
0 Y/ W7 L6 @7 Q3 V3 @- C...

国家数据：

USA 3118
' K' F: ^- Q$ z- g; m# yCHINA 1493
FRANCE 1385
5 c& y; l: x7 @0 X. Y; fENGLAND 804
% J" X! c6 `5 p* `* l( c, R$ r. ~GERMANY 748
. w$ ^' E+ @) k) L1 R7 |# xAUSTRALIA 421
8 y. P. c/ @2 V  [7 U6 `9 k6 \3 TITALY 421
3 h  W9 G7 ?) A+ T$ P; zSPAIN 371
  T" p2 g& w8 HNETHERLANDS 340
9 G  J4 P- r0 DCANADA 323

从地名到坐标

( W0 S: B$ v2 R& E! [/ w: d

使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。

参考：www.52ocean.cn

得到：

Affiliations long lati
/ n/ {: J' c0 @# V% Q0 k1 ENATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.883056
0 E1 f6 [2 [4 p8 X! ]9 e3 J+ a. |CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE  2.2639934 48.8476037
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT  2.3 48.9
CHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327
UNIVERSITE DE TOULOUSE 1.4 43.6
UNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586
7 }5 \& g& a. `5 [. N* ]- ANASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65
NASA JET PROPULSION LABORATORY  -118.1713944 34.2013081
. O! I  R* J, u9 g$ v

绘图

1 K2 V- r$ b. u: F2 |4 N6 ]* m$ _/ G1 p

数据整理完毕后，使用GMT绘制3D柱状图。

REM For windows
3 I- B" l# G' p2 ~2 d5 \
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
# K0 S; F: ?0 Q5 agmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black

set ps=3d.ps
gmt pscoast -JQ0/37.5/10c -R-180/180/-90/90 -Dc -A10000/0/1 -Sazure2 -Gwheat -Bxa60 -Bya30 --FONT_TITLE=10p -K  --MAP_ANNOT_OBLIQUE=45 -p200/30 > %ps%
9 G# D3 S8 l3 {) z% jgawk "NF>1 {print $1, $2, $3}" loc | gmt psxyz  -R-180/180/-90/90/0/1000 -J -JZ1.5i -So0.1c -Gred -O -K -p -Bza200+l"Publication numbers by affiliation" -BWSneZ -t20 >>%ps%
7 ~8 N6 D$ O2 Z! H- \
gmt pscoast -JQ0/37.5/10c -R-180/180/-90/90 -Dc -A10000/0/1 -Sazure2 -Gwheat -Bxa60 -Bya30 --FONT_TITLE=10p -K -O --MAP_ANNOT_OBLIQUE=45 -p200/30 -Y6c >> %ps%
gawk "{print $2, $3, $4}" loc_c | gmt psxyz -R-180/180/-90/90/0/3200 -J -JZ1.5i -So0.2c -Gblue -Bza500+l"Publication numbers by country/region" -BWSneZ -O  -p -hi1 >>%ps%
gmt psconvert %ps%  -A -P -Tf
) q/ E* t' E+ z" o7 s

结果

2 l4 S& K$ D3 l' {' o* t9 `5 s

图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。

9 @8 ~% E3 d1 {5 x  k) u( A

附赠平面图绘制代码

& [0 w/ B" R: d: X
! j- u& B: O3 P1 _6 e+ c0 C9 |
  o. B5 [3 y; F* q4 zREM For windows
3 ^; |% ~6 q% D& V, F
( e7 S5 S# r2 t8 q  Hgmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
( R4 d! ~* A# ?% \gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black

: |# L' O1 Y2 H1 s7 qset ps=fig1.ps
; ~! s8 ^# a* V, i1 Qgmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K  -W0.1p --MAP_ANNOT_OBLIQUE=45  > %ps%
7 z8 _, r  j8 a6 e4 l* Ygawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J  -Sc -Gred -K -O -t55 >>%ps%
gawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J  -Sc -Gblue -O -t75 >>%ps%
" l$ Q+ o1 Y4 f; R1 [1 k2 D
3 l2 V* g: C  {- V; L( X1 ~4 Ygmt psconvert %ps%  -A -P -Tg

注意国家发文数量做了处理（log）。

图件缺少比例尺，有兴趣的同学可以修改代码加上。

至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。

& a9 u3 h0 I% \, u) H% r9 s' r

样本数据下载

4 `8 B) h$ {! k  [* o& c3 d- o# \

https://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA