本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。2 K) j' c5 x( ?: L) W; i
WoS数据 % L# E% ~3 A; v2 w" ^, S' V
根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,5058 Y: Y. _1 R" O5 e/ ]1 P% _. D; s- B
NATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.7478 g% B& w& I3 @1 J" Q1 b9 V% L
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735
; Q2 h9 o3 S0 O1 CINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.313
; x+ R) m" U) H4 ZCHINESE ACADEMY OF SCIENCES 629 7.396
' P- B# [, b* l6 GUNIVERSITE DE TOULOUSE 517 6.079) v, p! |6 G7 ^0 M, c3 ~* ?6 m
UNIVERSITE TOULOUSE III PAUL SABATIER 498 5.855
1 L, r M$ U G# q, aNASA GODDARD SPACE FLIGHT CENTER 483 5.679" j4 h- a; @! q9 v$ J) a8 V
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162
- }) l8 N/ M- u8 f! y N0 uNASA JET PROPULSION LABORATORY 431 5.068) q% n5 D5 l# `+ G6 T
...国家数据： USA 3118
. }% |3 K" U6 F7 v ^: ~3 p5 p- dCHINA 1493
4 J' j- j% u* l/ |7 B1 rFRANCE 13854 @2 R4 c# @, N! W3 z, v8 x
ENGLAND 804
; L5 ]7 q/ V& E, P, `7 @GERMANY 748
8 T7 P1 J1 X7 O$ k( ?( X; YAUSTRALIA 421/ r. w) z. f5 T! F
ITALY 421
8 A8 V4 @2 M2 q9 ]; E( USPAIN 3710 i$ \: I" R+ Y4 x8 D- M5 R G! I& |" w% \4 p
NETHERLANDS 340( o# y1 V; ]# c j b0 K4 j" [7 [
CANADA 323
4 H( N7 V+ S: i! y- M, x B
1 p, P1 _7 s; A2 A7 ~+ M3 A+ y从地名到坐标
7 T" X$ x Q! D; ?7 {7 }' f使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati/ ?5 Q. C8 S% V6 }7 _& W
NATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.883056
l( b$ H8 N+ y( y& b! FCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037
: E9 S; D; N" Q: bINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.9
: \) x& n8 g. X/ nCHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327: e4 w% g2 u7 |8 l+ z# Q
UNIVERSITE DE TOULOUSE 1.4 43.6
* u9 p* X. k! j& G5 g L" {UNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586" N6 m g; v# M( q
NASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495 S( N/ ^; A0 e7 X! {# a
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65
) M8 r. N4 M+ Y5 e- f/ ^1 VNASA JET PROPULSION LABORATORY -118.1713944 34.2013081' @/ f w7 I- h* ?
2 i, c7 Q0 C* B% A
5 O+ o" |5 e5 k0 D5 b8 d1 X- \4 J
绘图
1 `& F' y6 S# a! ]. G3 ~数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows0 {, q( k5 k3 a/ m$ _1 F; t( u
0 g% W% F u' R
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
% G$ C. y( E u- s0 b6 sgmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black # ^: i1 [% q- \( [# D. ]" J: r8 b
2 k4 W# n* c8 ?9 r8 z+ Mset ps=3d.ps
4 W" l9 W; h. L3 Q- P4 Agmt 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%
' D ?9 Y! R' N7 L! c' I- wgawk "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%
2 c' t: c0 a' a, t- J
: E j G( {0 o& e- dgmt 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%1 q5 C9 M8 c' R$ t' l! A4 M, A$ 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%$ k. v4 k+ J" B. V
gmt psconvert %ps% -A -P -Tf" [. \: c" a% h4 l- {$ q3 p. w
3 v; a6 F2 y, r. q, f
结果 $ ?% b* M3 g8 n' x! a0 S
图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。 ! Y. }% a: e7 P+ s) ?& m
附赠平面图绘制代码 ' ?$ Y6 A) L8 h; X
& O# \1 L( O* P6 h9 V: c9 ~1 E) G0 I1 q+ { A& W' Q- S5 Y
5 c9 `# J* o; @/ m
REM For windows/ Q# F6 w6 }1 T9 i6 v9 Y) }2 h) e
: k* g8 n4 C: o- ~gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p1 Y m& t2 X: Z$ W6 x+ u
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
7 k. K2 W$ V- M; p: K. \, x3 N6 Q3 C- F) g
set ps=fig1.ps: @" k% W1 N: {+ N9 l5 O; X" L
gmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%/ t+ s& k; j$ a/ U
gawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%8 j/ R: n" ]0 ~1 ?- C% D- l
gawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%
' I" }3 H. e0 l1 O: r: H" w: Q$ h5 E$ q+ h6 Z- I1 _. O j6 |
gmt psconvert %ps% -A -P -Tg- l3 V3 j: t! ?" E6 ], L* f6 @
注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。
+ p y5 g3 j4 i8 a5 ^- v# L4 ]8 }/ r- p8 {' P
样本数据下载
3 T# Z' h; d( S; u! R7 h5 chttps://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
: `3 \4 n0 a2 L8 a+ m
/ g+ S% C: o: A: Z3 ?+ W& K/ a; C8 W$ h/ P
|
