本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。
: v* s" z! M% s0 z: U& GWoS数据
# w; l4 i5 Z# k! u根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,505
5 y: h/ s {7 Z, M: ~- kNATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.7479 L& v3 X# F; t. O( d
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735% b' k" ~0 R8 o4 T: |
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.3132 W, h) O: Z7 c
CHINESE ACADEMY OF SCIENCES 629 7.396% }' `, w4 B+ E; K) b
UNIVERSITE DE TOULOUSE 517 6.079
2 I$ f' M* r) ~8 B4 mUNIVERSITE TOULOUSE III PAUL SABATIER 498 5.8553 |, x3 t; Y# M" @9 S6 T0 ^
NASA GODDARD SPACE FLIGHT CENTER 483 5.679
' e& T' h9 w4 ^. ELABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162
2 B# [3 ?7 i" F9 Y- ]" FNASA JET PROPULSION LABORATORY 431 5.0681 ]( ]' e1 r0 m) Y$ Y
...国家数据： USA 3118
/ b* r' t7 M7 Y' X( ICHINA 1493
% S( V: s6 v9 q* _: v1 y/ V x9 sFRANCE 1385, C5 `! O! t9 S! O; h
ENGLAND 8047 n6 G; x. a9 |+ E6 X
GERMANY 748
/ @( Q$ k! a, j+ b/ fAUSTRALIA 421, o4 a. F9 j b4 n# @
ITALY 4218 M4 N( z) ?5 J6 @5 n9 R
SPAIN 371+ Q% @8 G4 j/ g
NETHERLANDS 3406 r" e0 ? s3 X/ f8 G
CANADA 323
0 M& x L9 q. z1 `: Z; L, I" Y4 C
从地名到坐标 1 I3 T* V, ]: x( t% h
使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati) q( x0 }" ^/ b% @
NATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.883056
. n: y4 p* ^% I1 l' iCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037
: u+ P! Q3 m {0 s7 a3 `INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.9
; V3 i8 k8 ^- `5 `CHINESE ACADEMY OF SCIENCES 116.3356321 39.91063272 u+ O# q! I- R3 Y$ }2 H- v
UNIVERSITE DE TOULOUSE 1.4 43.69 ^6 M) Y0 @ `/ t& _5 V
UNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586* ?1 d2 F# `7 \% x/ i9 H, k9 N; z' A
NASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495
1 G1 |& z4 r; S$ pLABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65$ W9 A/ F$ O& A( T0 T
NASA JET PROPULSION LABORATORY -118.1713944 34.2013081
4 i, E6 y" S1 | G, E3 C0 v6 S6 W6 m5 i5 v8 N; |# S; @% ?% u0 ~
]. q9 _9 ^. A$ l! b
绘图 & ^" {+ s9 u" }. P" K% |& S. H
数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows' o; X1 w6 x* k2 v5 X" a9 w. g9 \
: ~* Z1 T' a# E" k ^
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
" x5 N n% ?. a& Ugmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
; v1 K1 S$ T6 a. Q- \% F+ B/ @0 f! _; B
set ps=3d.ps
" [6 _& g! U3 S- B {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%- y; M) |, W! ` c, x6 a
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%* {9 S* f Z6 w* w3 y3 r
1 ^ W, X" l! L! t0 b
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%
5 ?9 n* \6 V! q# ngawk "{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% h% A5 P4 p+ e# g
gmt psconvert %ps% -A -P -Tf# @' U# M' l1 H8 Z! x# E- _! ]- Y& J
. i! s* @/ F: L7 L" M结果
& t$ P+ ]" D5 a5 p* i- I图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。
1 S% E% j W9 T( ^; |% Q# C
附赠平面图绘制代码
. v' M3 _5 J5 v1 q2 t0 ?/ [: n. l5 l; O( W! ?
) w! w z9 J3 X( t" i- u5 B7 \$ n) _1 p! U9 I' f/ M' N
REM For windows
0 j! e* i7 P' O+ b0 \6 ~1 N/ o
5 ^( A8 h/ i* L" D5 b: T7 Jgmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p0 c, T' D2 _7 o! C
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black ' ?8 p: s3 {: p7 f' { s
% G ~: O2 U2 _0 Zset ps=fig1.ps
. |- J3 a1 }* M S$ kgmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%
5 C, ~5 H1 Z1 k2 Q: lgawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%* S) U# U6 H$ G, H% ~
gawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%: f5 K( ]! J6 a4 C2 i2 b" G$ p
9 U. M& k) x8 E. M; V- Z" _1 R
gmt psconvert %ps% -A -P -Tg% q# a5 Y$ `" f% G
注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。
% @; B" A5 x! B' r. J+ `5 {' E
样本数据下载 ) A" J1 {" A1 K, Y! O
https://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
( O [ d* o+ W) @3 H2 Q; w
5 Q: b N8 {2 f. I+ K% p, J' ? |
