本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。
' J- v# J$ _ P2 HWoS数据
6 T( ]+ v3 A+ @! W; a根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,505
! @, l: C6 L8 E6 `9 LNATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747% n. C, u/ I% x' i! d
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735
7 }$ ~" I6 M+ J. o \INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.313- ~2 P# S" J5 W' I5 j, W
CHINESE ACADEMY OF SCIENCES 629 7.396, o H! J' g0 h0 Q4 d9 a
UNIVERSITE DE TOULOUSE 517 6.079
+ m3 t) J4 ?# PUNIVERSITE TOULOUSE III PAUL SABATIER 498 5.855# I% _9 i6 `& e6 i( B8 K& X
NASA GODDARD SPACE FLIGHT CENTER 483 5.679
6 A1 s0 @( h/ oLABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162
3 P4 k2 W: o, o5 _* D) ]NASA JET PROPULSION LABORATORY 431 5.068* l* L+ t* Z4 f' k
...国家数据： USA 3118
8 g c+ V& t9 V& y; Z- n# M% nCHINA 1493
! P1 l- U5 ?0 Q, j' v8 ]FRANCE 1385
# o1 K# o" Z1 A( wENGLAND 804$ g! g0 G' q& u7 `! [- N1 ?
GERMANY 748
% T+ V; L' \0 ^8 O- }4 W! ?2 R/ o- |AUSTRALIA 421
! c9 ~5 A$ T* n! C1 [5 B7 \ITALY 4218 x. H7 j Z0 [* c; S: Q
SPAIN 371
. k6 `% R* z6 sNETHERLANDS 340
3 ]0 @2 l$ `9 ?$ P0 QCANADA 3231 l T! j" n- }+ u2 ^6 H
) t4 c$ R/ f% N! {; K- Q
从地名到坐标 T/ ^+ ?4 R# X$ F
使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati
: v! X' T I! }" T9 C. N6 P: hNATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.883056& Q8 t P3 B2 a
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037
* D4 ]' w2 G1 Y4 ]4 {8 ~INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.9, F/ z/ \8 [1 s; k2 Z* \
CHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327+ s3 Y+ E3 _. M/ U# j+ k
UNIVERSITE DE TOULOUSE 1.4 43.6
. }) ^( Z& ^- R0 Y6 M! R6 FUNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586
: f, A" y6 ~8 y0 O2 q' g. wNASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495. _0 t+ I4 F3 i* }+ V7 X
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65; y; b9 p" T5 h' L& B A" L
NASA JET PROPULSION LABORATORY -118.1713944 34.2013081: \+ E# u% d% _, l! P+ ?$ [
( Z' s& o' `: O3 N* }/ V/ z
+ z/ b5 b/ J# W3 J% d绘图 : e9 J; k3 y) V9 Q7 l
数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows* `5 J. a+ U; q
* ~2 [+ f5 M: a$ y! X7 y) i- _
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
+ V% X0 \& u( [" u. M8 i2 _gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black 9 Y+ M$ \. ]0 s t3 v6 @7 Q
& P6 S q$ Z' F. N6 z
set ps=3d.ps* u) ~ `4 K2 N5 o* O8 R& m
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%
1 i1 z9 u% l7 j6 J4 Pgawk "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 }" i: n' b4 G5 W
1 v5 Q6 O0 i; \' A
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%
; g2 ]6 E/ C" \8 a: j6 Tgawk "{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%$ m4 X3 d2 h7 S! P, F; T" A
gmt psconvert %ps% -A -P -Tf+ Y! L" P4 L/ \1 ?/ l1 z
) @/ g( ]) j, Y1 d% ~0 j6 z$ \
结果 ) N5 |# Y, i/ V) ?3 }7 v
图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。
$ a: @1 K4 Q+ W
附赠平面图绘制代码 0 {2 b) t5 Q" y& s; _1 v: h+ H1 C0 T5 D
$ c( R3 a- G3 F' g v1 V
9 S; t% y; k. }* k6 Y2 j) ~
# F& L- o. K& n9 x7 h" f! Y5 e CREM For windows0 v& W7 s0 h3 {# X, G
3 S# t& y/ ~% q! S% d" Q$ D
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
0 a/ j2 C5 _, n& f! z7 W& `; [gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black 9 q0 P0 E+ k4 @8 x: j& Y) h
( D2 D3 w- @$ c& v0 N: oset ps=fig1.ps4 D4 A4 W. }. X8 x
gmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%! Y& U$ s) q, L h% q5 y+ V
gawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%
4 t% V3 ]& f$ B; Bgawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%6 k5 Q7 O2 D4 K" F' j. }
. \& @ A4 j7 s/ f
gmt psconvert %ps% -A -P -Tg/ L( X0 ]: O. p
注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。 $ H7 q' y( c* t |
样本数据下载 9 f( C5 ?. k% Q6 i" O& g
https://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
. S8 z" P# j( A3 f k; x1 F
( ^6 P) Q8 T! i9 I, S |
