本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。' X! S8 t4 [" }5 v& ^* b
WoS数据 ( Q' H6 |3 n: \5 @) R' g
根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,505
/ o* J5 D2 {( t& B1 G' h+ NNATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747
' y8 \0 k/ ^6 D* _9 G, ]CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735
$ o. H6 g/ w8 D$ K2 s0 HINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.313
0 Q4 V. O6 k* F$ @: t/ aCHINESE ACADEMY OF SCIENCES 629 7.396) ~ I6 h& X2 d% x
UNIVERSITE DE TOULOUSE 517 6.079
6 ?( q9 s g% G( IUNIVERSITE TOULOUSE III PAUL SABATIER 498 5.855
# j# N# v! |, ^4 ]NASA GODDARD SPACE FLIGHT CENTER 483 5.679
$ M7 ^2 }5 T5 b0 E' `& ?LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162, U9 w+ a) R9 P, _6 ^$ X8 e3 h
NASA JET PROPULSION LABORATORY 431 5.068' }3 \3 S4 L1 ?! P/ { r; ]* ~
...国家数据： USA 3118
8 Y" T) F) z* ~2 I( v7 \3 FCHINA 14932 R" F. V9 K$ z& x# Y2 V, N$ z) e
FRANCE 13858 S3 J: _9 P& t0 o; S2 T) ] l
ENGLAND 804
) E; y- h5 F/ b# b- RGERMANY 748
+ f( a: _! V6 a! V9 e% XAUSTRALIA 4213 ]6 h3 F8 B( o6 K3 {6 V
ITALY 421
: I* Y' g5 D/ d' S# n- k9 ESPAIN 371
' w: y. [ n% b3 [$ vNETHERLANDS 340! y4 y3 g# d- t5 X/ y6 Q- m
CANADA 3234 }- z& m8 P! ?
9 E. D" z/ q Y5 X7 J8 ?
从地名到坐标
: g0 t1 j3 E) R使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati" d" i1 X% Y- g3 p. b
NATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.8830562 ^( S( ^7 F/ S$ z. \+ n) e
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037
% c8 P& V( a, \0 Z, q$ HINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.9% }: l$ c4 l% w' G
CHINESE ACADEMY OF SCIENCES 116.3356321 39.91063274 W. ]1 w: ~* W/ g. f$ w: S
UNIVERSITE DE TOULOUSE 1.4 43.6* j: H3 S# g& J1 H& E7 N) m; }( {
UNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.56195861 q$ x& |6 H* C6 _
NASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495
) V/ x& \0 }' ~* P" H3 s5 tLABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65
0 D0 g; G5 N2 \ N4 ?1 {NASA JET PROPULSION LABORATORY -118.1713944 34.2013081
8 o) O. J& W- }. f& Y# B: h8 W! ]( u- ]: G5 W a
Z+ {7 Q( n( I2 Q+ ^. k
绘图
8 G) y, c# v5 V( W数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows
6 h9 ^3 Y9 Z X1 q+ ^$ U& O& l2 B6 |) o- J* E
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
% E* @' h( s' r) l' T% g" wgmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black + q, L0 l% }( w, Q
9 V- m1 W& w6 S- S( kset ps=3d.ps J# ?6 B* z9 M w
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%
/ O2 |! _8 \3 i9 G, b4 H" ~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%% O* ~) D, Q; U0 |
, @( o& B5 T+ P6 fgmt 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 G. R/ C" g9 r
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%; |3 k; X( Q4 D! A( u, r; R
gmt psconvert %ps% -A -P -Tf$ y- l1 U$ T& }) A0 }, X- R
9 @% ]( h8 ?! O! x( e+ }. y0 W结果
* r! v4 ^- q- D4 H图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。 ! j/ ?6 S3 }1 D
附赠平面图绘制代码
( T7 z5 l/ i. B# O- k& O( o, U: S' g7 s0 ^
+ t" M- S7 Q) m" q( }5 C' J
* p% I$ _# m* Y" q. O
REM For windows' q8 {: F C% E" S6 H q n( z
8 ~! g U3 f2 J n
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
# o3 Y) }; t" e# o- }gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
8 L$ W; }: ^% b/ f( Q; q, I, m
* Q5 T. t' T# U1 Q2 l& Jset ps=fig1.ps- h; e2 M) J1 U4 Q# q
gmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%' b9 z( w: h7 m) @7 y: o4 n% n
gawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%
$ y, G5 F3 z! b' C$ ]gawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%
8 P9 ?' t; _5 h' b- V- _" ~* M8 z; t! \( I& N
gmt psconvert %ps% -A -P -Tg) M& J0 n& j7 S5 J/ v& S
注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。 2 o# }: b+ U2 [( J
样本数据下载
8 t* h; a8 f4 d- N: z4 G2 {https://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
% D2 y& E2 P! J6 ?: \
$ b" V) P M+ M3 d8 y
|
