本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。3 M# Y; m J# D/ t0 T% ]5 j6 o
WoS数据
& {' y2 A4 T/ T5 S, {5 n5 I根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,505
3 q" Q8 V7 C+ h# c# p. V* [! C8 INATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747
, Z! g# Y" _& xCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735/ ?" y8 c9 h4 v/ w) K& u
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.313
3 f" U5 {# a; S& M; c( b D3 \8 FCHINESE ACADEMY OF SCIENCES 629 7.396
: ~ B; h5 Q v1 Q: `8 H6 LUNIVERSITE DE TOULOUSE 517 6.079: {8 @9 O: A0 l9 e8 V
UNIVERSITE TOULOUSE III PAUL SABATIER 498 5.855
' t& ^ N$ f1 B, {NASA GODDARD SPACE FLIGHT CENTER 483 5.679: c' j: R3 q, `& }: S0 B! y" j( D- @
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162
' P3 \$ ]" B: l, nNASA JET PROPULSION LABORATORY 431 5.068
0 W' }% t7 Y5 M( `: Z...国家数据： USA 3118
g) X* t7 j0 C4 f" w3 sCHINA 1493: ^2 ]% n) q8 K( ^
FRANCE 1385
! ]( V Z7 b0 o* l$ i- I MENGLAND 804
# c4 w. } }8 K8 _$ X J$ UGERMANY 748
8 R6 e, ~4 H6 V! [. I' XAUSTRALIA 421" @0 Y2 k4 B1 V" j$ t
ITALY 421
, V1 ~, {+ I9 U. N1 j6 U/ LSPAIN 371
3 s) y1 F5 E6 lNETHERLANDS 340
& X7 R* |4 m2 A: M3 n# m! nCANADA 3236 f2 ^; A k. Z+ C! a. q6 y) [
1 u, R8 @7 @5 m1 X3 K! R从地名到坐标
" M: p r! G) _; r使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati
; }/ V& V4 W! e4 p% ?5 A* k! eNATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.883056! G9 l4 V3 q* q2 X3 L1 L) w0 [( h
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037+ Y9 O. x$ ]- _+ F
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.97 h% L) H) }/ ~
CHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327
2 e/ y# K q) M, U/ ~/ x0 `, gUNIVERSITE DE TOULOUSE 1.4 43.6
+ Z: P' p% ^6 ]7 w2 LUNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.56195868 b7 e! y) d& c. z3 s
NASA GODDARD SPACE FLIGHT CENTER -76.852344 38.994953 i/ ?/ v' |+ F: X" _" Z- h" q
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65
@) n9 H" y; g) H G8 ?" xNASA JET PROPULSION LABORATORY -118.1713944 34.2013081
) M- n; v+ L' Q1 [0 t! b# B+ D
$ }0 z( J. d$ E3 `- d0 C
8 W: Q* o% P8 S6 K绘图 - l" K6 X) g, R0 B
数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows V+ Q5 `8 j# R1 G& j ~
8 i7 g* N1 @6 Y8 \ [gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p* ~' t1 @) i1 M4 Y
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
) \8 k( C3 _1 A9 N/ N, S
7 m+ }/ N; x- s: Lset ps=3d.ps
( L. F/ V3 b# n, x& Vgmt 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%% R! G. O$ v% ] `2 m5 k
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%5 r( w$ X% n0 \2 e5 g4 X, B
* h! M/ p; P+ E, @) w2 ]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%) Z+ n# d# t7 g
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%5 _0 ~" K' f1 R" w
gmt psconvert %ps% -A -P -Tf; q7 V5 O* f' W! v: d1 O; R! A
6 M' ~8 U+ ~/ }5 _
结果
# I. E" W: n# b! M1 Y' K) B" d- ^2 ~图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。
1 J; e4 _5 L& D7 o
附赠平面图绘制代码 / C. z8 ?2 {! Y% ?# ]4 [' ?, ~+ J* h
- a& z+ C$ E7 k6 B& v* ~' a9 k
4 Z0 ?& m4 v1 p+ l' I/ }5 Z) y' ~
REM For windows
% N5 r5 m5 n3 M$ A; c$ s3 [7 J6 |" T6 n! s/ v6 h- [2 n2 u/ `
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
9 q% r4 H9 y3 n/ C0 \- H6 Agmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
- F. j2 W# w, U2 \% p0 ?* Y" ]7 E' Q w$ _/ G/ e
set ps=fig1.ps
; \* ~7 }/ K7 c( G. v V) v& w, ygmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%, I: G7 H+ a1 h9 Q m/ I7 d2 S
gawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%
4 d5 Z7 d& l" W( x, S" v/ ]gawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%
! N. l, ]. ?' B+ d
# l: m& P) W* F% @gmt psconvert %ps% -A -P -Tg
L. i: u8 L& D# e! a8 t/ f; i注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。 ' d8 U6 J3 W" b. N
样本数据下载 : X6 b( y4 ]5 p8 Q
https://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
9 |( H) V, N; u& e# z4 `, O: ? K0 y% O1 V& P6 S
|
