本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。, h2 Y: a2 b9 K. Q
WoS数据
, ^, ~' o" x# \( |( `4 @根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,505! m2 V' u- s$ u; L& y) [& @( H1 j
NATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747
/ }1 v2 T& m7 Q* H; [5 {" VCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735
8 s. h2 T5 @+ S, kINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.313% M* Q8 ^5 J9 h+ q
CHINESE ACADEMY OF SCIENCES 629 7.396/ ?) ~5 C* B1 P# r
UNIVERSITE DE TOULOUSE 517 6.0791 G' W Z9 |3 J4 f
UNIVERSITE TOULOUSE III PAUL SABATIER 498 5.855' @+ _) {) F4 h( z
NASA GODDARD SPACE FLIGHT CENTER 483 5.679" j; w |9 S* V- S& [
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162! L+ y4 W" w1 `3 ]/ h
NASA JET PROPULSION LABORATORY 431 5.0685 B. C. `) _% `/ J3 c1 b7 C
...国家数据： USA 3118
5 c2 w6 x" Y+ ?2 D7 F* ?4 ECHINA 1493
) }& a1 |+ s) K7 T4 u. L1 HFRANCE 1385& f% W! j$ K& e+ T4 Z3 k3 E1 U
ENGLAND 804- s, a: C' z# f4 Q4 T) ~* N% ^
GERMANY 748
0 x( F$ E& y0 k3 S& cAUSTRALIA 421: T/ n% ^% S3 @ r
ITALY 421
3 @( T/ P7 T; d n1 C/ YSPAIN 371
6 r( C+ A( m2 H7 R( A; }NETHERLANDS 340
1 _: K6 z; g7 i+ X+ UCANADA 3236 y/ Q, E' Y6 }0 O; u
' {! Y9 ]: v! U$ U8 V; x3 ^. g5 T从地名到坐标 # N$ |6 t5 Z- S
使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati
' Y/ ~: Q. s1 c' iNATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.883056; J! z6 }7 ~* }/ }- V% E. R; D
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037
- w: [' t( U; zINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.9, f) d# f! ~, ]' A, R, N5 k! y
CHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327; v6 f5 L+ ^7 i- C* V
UNIVERSITE DE TOULOUSE 1.4 43.6
+ f0 Z* b4 c! f( q' ~& _: h5 eUNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586
8 k5 e% N6 w2 _% ]" N6 ] t' d ~2 kNASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495
6 R u6 d6 F s7 ]3 E8 rLABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65
4 C# v( H8 ]8 n0 c h3 [NASA JET PROPULSION LABORATORY -118.1713944 34.2013081
1 M- W0 B5 `6 u" v- ~% h% h3 T1 H% @7 i1 w0 R
$ a! ]2 u7 C7 z$ ^$ ~" d绘图 + b2 \$ ^9 K) J2 S
数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows
& o* n+ R4 ] J0 [2 q6 d9 s% d, {& b9 t9 E/ P/ Z8 H
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
; r! T' x; {1 w, ^gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
) x% o& G. r' V- v( N* C* I( z" G% L% R& R4 v& E; Y+ A
set ps=3d.ps
1 s! z. D+ a+ F$ wgmt 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%
/ @- Q; v! {- r' s# G0 a' x* B- ?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%4 |9 L( s% F( Y5 N, _5 d8 Z9 W4 \
' X( B4 H9 d! s1 |) ^7 ^
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%6 e( r! S- h: A- ^3 J v& r# h# y
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 f" z& m; b5 @% @gmt psconvert %ps% -A -P -Tf
+ K% i3 Q, P" [( s# b$ M0 g8 e0 W0 ^
结果 / b* [9 o+ ]) d
图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。 7 Z; }2 g q2 j H' W
附赠平面图绘制代码
0 D! K+ a' b4 U& X) G& Y' i& B1 S4 P2 ~) P5 z. V
. O7 e ]; a: n4 X
& m0 Q, x8 c0 p5 D' n' Q
REM For windows
) F; U0 G$ x; B( x; W2 F v0 s; X: M6 \, ^+ B' n5 `) V3 Y
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
: a4 f' |* @$ L5 d8 ^/ Sgmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black & B+ {8 m' f c+ ]
: B ~' S1 J0 d1 U. {0 _8 lset ps=fig1.ps
5 q r3 K- [8 u' j$ K. Mgmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%3 h# p' U8 P& y1 b% g( Q$ `! ?# N
gawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%
5 J# `/ ^5 N' x1 `gawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%* V3 `5 L3 \ P8 O/ i9 ^
: c5 n$ U6 f" f) N$ |) l. [gmt psconvert %ps% -A -P -Tg
/ ^4 V, L, u- m; }. A8 C注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。 4 |0 R g; d/ t& C Y& _0 z6 o
样本数据下载
" ?8 V, o# w4 ~3 e8 lhttps://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
; h) Y r+ {$ g |: C2 T7 ^( _, @
9 d( U' y( Z' Z5 z+ E, H% v; x |
