本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。7 l7 X* B. M8 X( _' I8 |% }
WoS数据 {4 B* M+ ^+ Y
根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,5059 s8 S( C# s& O! X9 T
NATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747
1 _! S5 o8 R' [; T! r) XCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735' t) p, o, s, ~/ E+ m
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.313" _5 w+ l- c1 T8 t
CHINESE ACADEMY OF SCIENCES 629 7.3967 ~ {' Q& s6 M7 a/ t
UNIVERSITE DE TOULOUSE 517 6.079 A- j" Q! N' O* v1 u) R2 s
UNIVERSITE TOULOUSE III PAUL SABATIER 498 5.855
; K( w7 {) O5 e9 P3 `: Y! JNASA GODDARD SPACE FLIGHT CENTER 483 5.679. f; u* l9 v4 f
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162' m l" s( |1 g# E7 E8 v
NASA JET PROPULSION LABORATORY 431 5.0686 ~9 y' F( r3 M$ R+ M5 A
...国家数据： USA 3118% {+ l3 t9 N$ N* J4 J* ]
CHINA 1493
: ?& E% C. p9 U/ {- VFRANCE 13857 y6 e7 V6 f( |& V5 A9 Q, [
ENGLAND 804: d8 X D# y' o0 ?' T
GERMANY 748
, c6 u8 x. p. D+ Q# ^AUSTRALIA 4216 W' t# k; K5 R
ITALY 421
! Q* o' W$ r: Z0 z; ^2 Q, a( _SPAIN 371
2 ~4 _, n& Z; f+ Z4 c0 R6 nNETHERLANDS 340 K, X) c5 i! G
CANADA 323
" }0 M) i8 ]! n* Z9 p, y, Q d0 x3 b0 E5 M9 g( `+ H! o
从地名到坐标 8 T, z Q. x9 e( L( g. Z3 e
使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati
3 q( @8 c. H0 a/ V- l3 yNATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.883056/ ^% O' r0 H2 j9 `
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037
. `% P4 h1 y0 `( w# e, AINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.9
+ Z0 n7 I% U; O' X5 X' a& a% }CHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327
. R$ n% p) w( ` PUNIVERSITE DE TOULOUSE 1.4 43.6
8 F( x# s: p' X! N% N, b6 B0 KUNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586
* L2 Q6 g- }8 {) E' z9 |NASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495
! C" r5 e1 u; w# bLABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.651 j; j1 T# }, t/ O- {, }9 }: o( c
NASA JET PROPULSION LABORATORY -118.1713944 34.2013081* ] j: n. ` L" o, h' T* ~" k
) W- P3 s5 D6 T1 k; R8 s2 \* T3 v" s! ^9 e
绘图
3 c5 n! h- T2 r+ V数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows
) R- w4 @# B& w& F) v8 r5 K) g0 X" _: j5 b8 t
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
( s( r# R8 c. B0 B, Mgmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black 4 t% g$ ]* W' Z5 d$ Q% A+ r
/ S3 d9 e* L) M* E% {set ps=3d.ps4 r) F& b5 l+ j U7 U5 T5 e
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%) A" L; n% q3 z) W: @7 ^
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%9 m9 l7 d5 E: y n4 r
7 {; c: n- i! M1 S
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%) D1 k1 J; c" ^: {* @& f; d8 |# Q
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%1 Z; f( s# C( F9 a# \
gmt psconvert %ps% -A -P -Tf8 j5 i8 y/ F/ J1 y
1 t8 I, K. i* {. O2 i' a结果 $ G( s2 g) U! K; x5 @
图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。
. G% g6 }& \+ j: c6 Z+ ?* V
附赠平面图绘制代码 2 |4 q3 I0 R; g7 \. Y; u
* u9 s# E7 \: u2 V; v$ g
+ L( f/ o r2 w. T( R& a( k1 ]% K% Q$ R" S/ A7 w
REM For windows7 V# s; a; ^ H/ N3 i" @# w
$ e; g8 p) q5 \. b+ Q& t2 t6 Bgmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
; @" ?9 U# g1 Ogmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
C8 O+ d2 L# p& W2 L5 r/ e
3 h9 g1 L+ S) l2 s6 S$ m* z, Eset ps=fig1.ps( F5 f* r) S3 q2 _; _4 H
gmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%
1 ?4 ]8 a( V% N" \/ `# S& F% ggawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%: n* ]4 y2 ]; N
gawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%3 o7 m6 h3 n' D$ J" s0 s X
9 \7 b% |, ]* S9 U1 ~gmt psconvert %ps% -A -P -Tg6 g% ]* x/ w0 p7 J2 Y
注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。 5 f* ^$ y. ~; P* ], z
样本数据下载
+ T2 y/ R' g+ Z1 R0 x, Thttps://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
0 R7 d/ j$ @$ M% p
) n9 c0 O+ m: y0 w3 ` |
