本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。; X U: u1 R2 h$ W# i7 C5 M2 e0 n& s
WoS数据 , n3 D. z( {7 J e# H( L+ |
根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,5050 q) e0 t' v# J" K
NATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.7472 i( }+ X; l" V# v" y$ T" M$ g
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735
w7 m; h- w1 N D _2 xINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.3133 v$ K( N# u |6 d
CHINESE ACADEMY OF SCIENCES 629 7.3961 g4 v6 ^. Q) G# a5 [1 B
UNIVERSITE DE TOULOUSE 517 6.079' |' _7 E1 f9 y8 N+ @6 A; Q# r* X
UNIVERSITE TOULOUSE III PAUL SABATIER 498 5.855* P1 @4 y( N/ B7 K I) o
NASA GODDARD SPACE FLIGHT CENTER 483 5.679 _' S9 R8 c( w- ? p
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162
7 {! R& A |" i, F7 n' }NASA JET PROPULSION LABORATORY 431 5.068
7 u2 d3 P. T4 p, O. z5 F( k1 f...国家数据： USA 3118
1 K4 s2 K5 Y" ^5 h6 a' UCHINA 1493 C% T, L/ L3 c4 x
FRANCE 1385
- T# F% _+ ]& {( ZENGLAND 804
, ^) B+ b5 X1 F+ j0 h; G0 g: ]GERMANY 748. c+ v* S2 c; J* Q ]) U T6 t; N
AUSTRALIA 421
) M# S$ h: ~6 PITALY 421
% m& m' e; {1 d# S; }SPAIN 371
7 l& ^' w5 F& r, G# j2 `8 G0 o* ]NETHERLANDS 340+ v0 x, e; s% d% L. V+ U( }$ ?- }
CANADA 323" W8 n; {5 y' u2 V
- e+ ^- T% Q7 Y$ v' y从地名到坐标 ' z. k0 f; u, s8 U/ U: Q
使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati! c) j/ p: R3 }2 W6 d
NATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.8830563 h( D' T% g* D
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.84760373 j( \: P2 P L {$ v
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.9
5 ]) |; H, ] v0 b$ cCHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327
7 T! A) |- W' d) OUNIVERSITE DE TOULOUSE 1.4 43.6
. }# Y& S' S4 j% d6 HUNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586
) l6 d2 M# P" }, t, vNASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495$ N1 n* @9 p5 B3 ^7 Y9 S4 W
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.657 A. p6 Q' d. q- T5 r. A; q0 J) B
NASA JET PROPULSION LABORATORY -118.1713944 34.2013081
) @8 u, a8 O. x
' [3 J8 p l( [; w% J5 H
2 ]6 Q3 j5 R7 N9 O绘图 % L! r; ^2 o' k, N
数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows9 R- [8 x& C/ G. V% c" T x8 n
7 G% l/ O6 [( G8 g, d, m
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
0 S& T% z/ m* j. D! J% g6 agmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
9 D; m4 y3 u) V1 H: E5 H0 g/ J& @0 D, o$ |* c' o
set ps=3d.ps
9 i; {" @) T5 B: u* K/ hgmt 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%
) C: [4 o1 u; S$ w2 S3 b, j! Q/ ygawk "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%, w7 {& n! C; d# T! b: M' K! v
" z) h4 T: r) {% b4 rgmt 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%' x9 |' i+ ~. a ]4 j6 u
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%
+ Y) C+ }* A/ _* n" p" Pgmt psconvert %ps% -A -P -Tf
4 I& t0 C: t- @) Z
* ]7 b& n" e7 ]% w& B结果
3 f3 F4 ?3 q" p' B+ F4 Q. H0 Q图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。 ! i( l3 i0 k. p7 p! u
附赠平面图绘制代码
5 S; s* L/ n, S* U5 D: S1 F. |' I" G0 Z' E% I6 z6 F/ E i' N9 l+ X
5 n* w. f5 V0 H- U; a% ] H) \
, b$ B- O5 R# D! l. u" d
REM For windows# }' p' F1 O0 N) Q% ?2 o
$ M. l3 b% {0 u4 f, M, c7 v
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p4 M# K6 w# ~: M: z b
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
. u! E# ^2 I Z. d' n
# U5 J. M; `6 h6 Oset ps=fig1.ps
" Q, r$ t# }' Y7 Ygmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%
* h R8 g {5 Cgawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%% N6 d8 ~ V _- B) Z" S
gawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%4 S% `) z4 s8 i& P' ^
' @8 C! u4 z& J9 v* {5 I6 J$ h
gmt psconvert %ps% -A -P -Tg/ d- i+ R& ^0 [/ x" }
注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。
& a# F+ w5 l6 v* ?. x
样本数据下载
+ ~/ Z" M R( k) U0 ^* k7 _https://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
) W8 b+ Q" W" X& h0 R
8 ]& F5 P0 u" U; d- E |
