本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。
' j9 Q2 c2 o+ @, V7 ]; dWoS数据
, D8 y: O! t0 d, ^8 L) M3 l, S! M根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,505
# N- m0 j& K5 S8 G: FNATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747
@1 z" C% X2 t( p1 b% LCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735
' r. e( Z X8 l5 Z" eINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.313* W! E' T- A$ j i5 C6 H
CHINESE ACADEMY OF SCIENCES 629 7.396& u3 D( A$ `$ k( d
UNIVERSITE DE TOULOUSE 517 6.079( G3 Q4 c% u$ E9 b( K) P
UNIVERSITE TOULOUSE III PAUL SABATIER 498 5.855
) W* A1 G F5 Q7 ^) J" QNASA GODDARD SPACE FLIGHT CENTER 483 5.679
5 A U3 b6 p }( k% w/ hLABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162; p2 ?7 l0 a, Z- x. E" q4 v
NASA JET PROPULSION LABORATORY 431 5.068
' c1 I4 w2 C2 ~8 W2 g8 K7 d...国家数据： USA 3118" C+ U5 ~# N7 B$ J0 Q
CHINA 1493 W F$ V3 O0 q7 f7 c) [" u! `
FRANCE 1385
) U0 }$ i: m) n, B$ ]1 K1 w! ^ENGLAND 804; f- [3 R, K! p H+ l- g2 v
GERMANY 748
7 W8 F. t/ v1 Z& {6 qAUSTRALIA 421
* ~4 j. h s5 O% bITALY 421
7 e+ y* m: Y; ]$ MSPAIN 371
/ h1 d3 C3 Z2 WNETHERLANDS 3402 F* _$ R, ^$ \* v7 |3 @* I: S
CANADA 323: h3 B; Q8 r, R* C2 Y) W
9 O6 L. F3 ]! s; N0 v8 ]从地名到坐标
8 I( i( }2 _6 d9 `7 M& X& O; Y& g使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati& t# j# a$ N" R
NATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.883056
. {7 h4 i# s6 Z7 C$ R% GCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037
# t3 M2 m- _/ i$ v9 z+ f; K+ B3 \" iINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.9+ J! ]0 [& c6 R8 x+ Y" d
CHINESE ACADEMY OF SCIENCES 116.3356321 39.91063277 {" O4 c: e7 S5 P2 N2 Y* s; K
UNIVERSITE DE TOULOUSE 1.4 43.6
5 c2 s7 b/ v, K7 M* YUNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586
) {" E+ z5 E- x6 V/ vNASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495
0 d( A1 f0 Q M; H' VLABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65( ^& a% L3 f& i, p# v. x: a
NASA JET PROPULSION LABORATORY -118.1713944 34.2013081
; c/ y0 |( }1 n% B: {, ?7 Q E
# s7 Q& U; ?1 ]/ H& o# x0 B! p4 @ [8 D& l, R2 ^; y7 p7 F. `4 X
绘图
' }+ `* Z- @1 k5 n9 q; L数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows
- ~- Q3 ]' q9 H
; r" ?# e) E5 K M8 Tgmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
7 s1 B+ _- V2 L) B5 Mgmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black 7 a! K2 m3 K R) h) O4 i
! X; J! F* v, k0 G3 F+ c2 X2 @3 rset ps=3d.ps
; Q, j; [! f9 ?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%
P% m8 b) V" v" G4 @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%
2 _+ [, G5 [1 {) F0 J P
$ D f* U' S9 a: o3 a1 S: u! q$ f. w) @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 b/ d4 @) A& p' a8 M- Q* a, A
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%
y, Y% E1 O1 M: |+ }5 c1 l0 }# Dgmt psconvert %ps% -A -P -Tf7 g7 s: A% x) j, ?' G @7 ]/ V
/ l j& {; [+ C2 f5 E$ ?. u9 A
结果
7 S& b' j3 t1 l( {; m1 `图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。
1 H& I# \' o( n0 y# j4 `2 E1 r
附赠平面图绘制代码
2 T/ u. O& m# u7 H$ A3 ]$ m) [: j; y3 I$ E( `6 n7 E
r/ b' N" y- V' y' I6 ?2 ?0 m( `6 i. K: ]2 t! y. S- v; W; a9 j
REM For windows/ q% M) {7 c* B% {
- A3 o, j% r) B& s: F' Q& w
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
6 G- I8 ^) G# Rgmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
% g5 t k1 n3 c6 L; w0 a8 V% r* r# c# h7 a% `- u
set ps=fig1.ps
) h" r/ ]9 Q, }4 V$ q! ^) h7 d' Bgmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%
& O( \+ m8 ~+ s1 ~gawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%
! A5 B. t# r6 O8 i8 Hgawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%! Y; |$ U, d" D: Q; k& h
: D2 e& M) A# L9 P2 Z
gmt psconvert %ps% -A -P -Tg
/ F, c* y" N+ x" ~ Y: T1 A* ^, v# _/ h注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。
/ z5 d% h) j$ L' B4 l& ?4 z/ \
样本数据下载
, \" w R) z0 a3 ?7 g* a/ Rhttps://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
+ ]- `- H3 N; f1 s
- ?! i, x! Q K" h |
