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

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

WoS数据

7 X6 a. p2 D7 S5 p$ [- `2 o5 _

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

机构数据：

Affiliations Record Count % of 8,505
7 P* X& n( y5 E9 d. O7 ENATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747
# F. J2 R& Q% X- C5 JCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE  828 9.735
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT  707 8.313
CHINESE ACADEMY OF SCIENCES 629 7.396
3 H$ ~; F3 N6 rUNIVERSITE DE TOULOUSE 517 6.079
) T. D; L3 w$ k9 X' ?0 i& yUNIVERSITE TOULOUSE III PAUL SABATIER 498 5.855
- i8 M. c6 d+ i% \' l! {NASA GODDARD SPACE FLIGHT CENTER 483 5.679
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162
NASA JET PROPULSION LABORATORY  431 5.068
...

国家数据：

USA 3118
CHINA 1493
. B; F0 [% Q. g% c" E2 d- wFRANCE 1385
, Q+ e7 W+ R; hENGLAND 804
" o/ ?1 T  E  N* G# c2 @GERMANY 748
1 \- s, _% E" ^AUSTRALIA 421
1 }+ P# `4 Q+ f7 n4 b* qITALY 421
SPAIN 371
NETHERLANDS 340
! B) B& a+ k# d, \CANADA 323

从地名到坐标

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

参考：www.52ocean.cn

得到：

Affiliations long lati
NATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.883056
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
& c0 S/ P! j3 W5 IUNIVERSITE DE TOULOUSE 1.4 43.6
UNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586
. }- R0 J9 l) ?0 z' z5 MNASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495
* G2 H6 Z# e, `4 v+ i8 XLABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65
NASA JET PROPULSION LABORATORY  -118.1713944 34.2013081
/ S" [4 W. _! D. p( T, `' W
. z: ?% x3 W# U( h' s6 Z

绘图

9 j. v& r) C6 S, \: b+ t8 t; W& Z

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

REM For windows
; n: O8 S! q( ]
- O6 _% w& _; J/ I+ v5 B4 ]gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
: Q0 S# x2 Z* X& |3 Sgmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black

5 K: a4 Q$ _. N; R- h# Z$ hset ps=3d.ps
, m- z6 ^2 N: L1 D- o0 ~/ ggmt 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%
gawk "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%
  B, M0 i' p% k8 D( A
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%
7 A  G# j6 }3 g0 Egawk "{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%
4 P7 q) V+ I' K2 J5 mgmt psconvert %ps%  -A -P -Tf

! `" T: V* I3 ^: A0 H0 {5 Y

结果

  b( `5 Q; E3 p8 V5 a6 o

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

5 E& [" V9 w- T; N$ I

附赠平面图绘制代码

  K. y2 X0 x! }; F$ e% {, k; q, ~3 O

REM For windows

. c( d9 m  [7 {8 }  [/ P, Ygmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
6 ?: C8 t$ J4 O, v
: `9 ?4 M4 P  P! Y$ _7 _2 q3 iset ps=fig1.ps
gmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K  -W0.1p --MAP_ANNOT_OBLIQUE=45  > %ps%
gawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J  -Sc -Gred -K -O -t55 >>%ps%
( G* R$ N/ V: l' ~1 _& ^! Qgawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J  -Sc -Gblue -O -t75 >>%ps%
# }5 J" s" ^+ S, {& `0 u5 h: A9 ^
. o% Y+ j1 S7 h# j, ?( jgmt psconvert %ps%  -A -P -Tg
+ _0 i* @2 U& H  m1 O. [- F# p1 K) P

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

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

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

样本数据下载

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