本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。
* p7 Y: J) [+ q! d$ P8 ]WoS数据 6 p; `5 k! r$ b& a9 f) ?
根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,505
4 k U- N3 y# b! Q( w3 P- SNATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747
! i7 M5 q1 P) H, zCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735$ B3 G9 C3 \* Z0 K4 A
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.313
8 G; t4 n) M: G% L. ICHINESE ACADEMY OF SCIENCES 629 7.396$ h. }3 u8 p e$ l( [* y" ?
UNIVERSITE DE TOULOUSE 517 6.079
& {6 V* i" t8 u: \7 R4 MUNIVERSITE TOULOUSE III PAUL SABATIER 498 5.855
( `4 u8 `$ E' f7 j. |NASA GODDARD SPACE FLIGHT CENTER 483 5.679
0 R# A* N; Q9 z& a! l6 y" r' tLABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162) R0 R4 j6 \4 r5 ^
NASA JET PROPULSION LABORATORY 431 5.068
' T$ i% [0 m! o0 K8 w6 I9 h! \...国家数据： USA 31188 u; g! m9 Z* d! t' b. o7 c( @
CHINA 14938 y0 r5 k" t4 r8 d8 m, f0 `# |
FRANCE 13855 R( [* L( _! ~$ z# W7 F1 A) w3 D
ENGLAND 804
2 Z6 V$ e! B6 O. eGERMANY 748
, \8 ?& v/ Q& X N v" y( {AUSTRALIA 4215 O( k, L& Q; o2 N; Q8 Z7 `
ITALY 421
. D& u- V0 e. m; S" T4 SSPAIN 371
" o/ ]2 \2 d6 m1 |" ZNETHERLANDS 340 J. `2 u0 L6 _) L+ ]& v0 g
CANADA 323
. o* x9 K, V9 I" T) ] z$ o6 u: `2 h" J/ ^* Y
从地名到坐标
! ~0 q6 G9 Y* j+ S/ C使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati
8 Q9 @6 y$ y: Z VNATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.883056
6 e2 c- }+ p2 M- [7 V( }' ?$ YCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037
. V6 W' W- D; X. |INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.9
+ f7 A0 |$ F0 p6 eCHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327
4 _, l* B- M3 [6 o* rUNIVERSITE DE TOULOUSE 1.4 43.6
9 E' d& g/ q- R9 \9 _& }2 |UNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586% i6 z$ K2 B6 a" y# h& q" }
NASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495
6 Y) q$ L* l4 K. C: V0 D7 |LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.651 n9 K; t1 @4 i, ?7 R. t9 ?3 c8 ^ P% o4 S
NASA JET PROPULSION LABORATORY -118.1713944 34.2013081
9 u# r8 \) f# g, ?$ d( q! \
) u0 V4 L& _+ N2 g }
2 A/ z! |0 p, V绘图 # D7 h* {( j& b; h# E
数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows! {) ]( Q y% `5 n8 a
% Y! V" l- r& \. R3 S4 ~
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p+ X8 G; b, h$ R# Y
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black # d5 E& b* Q6 a2 @$ b
7 C0 ?: d ~! h
set ps=3d.ps
' X- ~' v8 f2 s6 z- U9 g6 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%
; Q3 o: b& b, J# E M$ ~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 q/ ^ j/ Y* }, x
$ {+ V3 X1 i' F. 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%0 P0 q( U! U2 l' @& I
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%
% C) B' Q3 @' e6 r% ~gmt psconvert %ps% -A -P -Tf# K: T5 U8 s7 ]" x# n/ g$ j
/ r8 j" U, D- ^/ G) {" P
结果 ; K( X9 N: d' i! X& L l2 W
图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。
0 {! T# r5 c3 f" X9 ]- s
附赠平面图绘制代码
( {" x. W. h, ~: ^' n7 u, V" `! Q2 c1 C' c, E6 A
* Q, U" t3 S1 B
4 L% A0 U6 S9 ?/ x# _) F
REM For windows% y& T9 m9 J5 Z; z
# p6 e+ Q: d0 L8 ]. W/ u! K$ V$ Bgmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p3 J4 I/ L4 b& F; s; D5 T7 o
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black # ?+ Q$ J6 D, f( N+ k: V
, m. o1 M9 ? `
set ps=fig1.ps, R4 ~" n5 x4 i& t9 e
gmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%
! T! z, w3 y$ X5 y* |4 `/ sgawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%( b- ?+ u/ g5 _+ v9 d# ?4 D, D
gawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%
/ \- j6 k" N; }, X/ D" @
% ~) O9 [& M; B C! U% ~gmt psconvert %ps% -A -P -Tg
# k& D+ k. h- B7 J注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。 ! V; U# G2 N& M9 M4 C2 S$ @/ o
样本数据下载 3 _$ i3 `" p1 I8 ~% n: K! b& f# H
https://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
9 \$ e) m+ e( K# Q5 ]' p7 G, K. D3 H
|
