本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。
- q! J: K# C; `. ^- kWoS数据
$ l! v( G7 ?3 f根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,5053 ]3 ]7 D" N$ t% x* l& r9 f B
NATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747
D9 v6 Q/ k8 _+ _CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735
' s \: P8 J% q0 g3 bINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.3137 C8 o w o- Q* u
CHINESE ACADEMY OF SCIENCES 629 7.396
]% ~% m$ D" [9 nUNIVERSITE DE TOULOUSE 517 6.079
3 Y. I' Z+ b0 A) h- y& p8 k: NUNIVERSITE TOULOUSE III PAUL SABATIER 498 5.8557 v. V" \& H% B3 S, M
NASA GODDARD SPACE FLIGHT CENTER 483 5.6791 j+ g8 n7 ~7 t: i5 E- ?
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.1624 Z7 ~1 k5 d' Q# L2 w1 {9 n4 I5 a
NASA JET PROPULSION LABORATORY 431 5.0689 N5 f4 F2 z2 q
...国家数据： USA 3118
. b/ x( J5 V! D$ _CHINA 1493' i3 C0 V& F& l' B1 k3 q
FRANCE 1385: E% e2 o! X2 I% e9 u& t, G$ O' _; U
ENGLAND 8046 H) ?2 o f- z" m0 y: D M, C- O* i
GERMANY 7482 Z3 m y0 l/ G; W" b
AUSTRALIA 421" a& p3 ?; R, B
ITALY 421
/ U R' l. S. I) g4 C2 GSPAIN 371
( \* ?9 M9 j' _& O3 Q# L5 [NETHERLANDS 340
, m/ \( M3 s1 s. r) YCANADA 3237 l! ^# @4 `9 u5 P# {- ?6 r
. |2 u1 `) j/ h; L. o
从地名到坐标
' d9 a/ P! R8 l+ S6 ]: Q- G使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati
: D; b* s. G$ O( [' v f. R2 DNATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.883056
0 P0 ?, b% `/ U; {. VCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037' `1 t9 W [: P
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.93 H" J7 P# x1 Q+ V$ V( ^
CHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327
/ |, C- O0 Y( p5 o! }9 ~UNIVERSITE DE TOULOUSE 1.4 43.6' @0 |: a# ^+ V2 H
UNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586
+ K! H/ x" Y) E- c0 v1 r0 [NASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495
' Q/ e. F& @. _' V* Z/ [2 hLABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65
7 x$ D( X" Y7 n5 [) h* D" \2 FNASA JET PROPULSION LABORATORY -118.1713944 34.2013081
7 Q8 ?0 c: y/ H
% Q4 ?& n- r! ?" p8 l, v7 n- y" A% d/ N5 K W$ ~
绘图
' ^6 }) E. L3 w, c8 r1 j: t数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows
' j7 G6 x) ~5 C! @- B, A* E7 }& c* {4 {5 [+ \4 k8 g0 `( M" S/ t1 j
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p; u2 k* N1 N4 t. e$ a9 H' W, `/ x
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
M2 R4 A- I6 l; `5 t8 H z0 m, M! S1 k+ ~3 M# y% o# \8 E8 f$ j, y
set ps=3d.ps: K* i3 h) \: }4 C2 }
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%6 s! ~+ R( |% s# R4 _
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%
2 Z' z- W# W! g- \7 s7 ]
6 [6 k+ Q, C; r* n9 R7 dgmt 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%
: x& ^9 {" m2 a4 S* Hgawk "{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%
' ^: C$ Q/ @& r& g" @$ Z) I9 [gmt psconvert %ps% -A -P -Tf- S8 U* w T. Y, {+ X
7 j4 x3 T& w; E# G4 j# D( M7 h
结果
3 f* s" q! I7 K6 s& ~: [图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。
) D2 c- j- p, G3 O
附赠平面图绘制代码
" x( _2 X7 s* @( Y. B; L
% Z% i1 ~- e) h a8 U0 o
0 `4 Q5 }5 N2 q* u: B& ?' `( D8 g% m g* Y4 H
REM For windows7 @; F. d4 v% o% z' z0 p+ M! A
+ \8 z' x$ g9 q+ n$ o/ G0 W& z
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p v3 u" B( U, @4 z8 t X/ R A6 ]
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
* I& Z) T i1 f% R; {7 ~$ s7 k+ H4 l, h O# J
set ps=fig1.ps
# E+ Y% ]& [; A- |gmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%
+ v& X* f' |1 N/ H; l/ egawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%: j: z0 {, X0 W5 e# V9 d$ X
gawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%
3 `, n" G3 f! I: n$ J! Q! A. W# [2 U8 g% W: A! A8 F
gmt psconvert %ps% -A -P -Tg
8 S: b: w" o' K6 F5 ?1 t/ Q注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。
: b. y/ H* [/ f0 H9 C( m4 c
样本数据下载
; k/ R/ ?2 b% G/ t; ?9 s6 U& ihttps://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
l" h1 v: f% I
/ _- m, D6 p" k& Q2 z8 z8 G% ^, U |
