本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。* N1 k# h0 @5 f) c8 A0 \
WoS数据
& y7 A$ I4 _3 B0 ^: e' i9 I9 D根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,505 t6 ?0 b$ W* u1 F6 g: e) ~/ z
NATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747
: O0 r3 M/ K/ PCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735
4 o' V" f/ b' I L" ]5 x3 n; RINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.313
: p# d, e" y) `" ]3 q4 [CHINESE ACADEMY OF SCIENCES 629 7.396
/ t9 h. L. m: a, v; l! m' L; MUNIVERSITE DE TOULOUSE 517 6.079
: p# D7 _& l4 X( `; c( ~( `UNIVERSITE TOULOUSE III PAUL SABATIER 498 5.8552 D2 ^+ H& \5 z6 G) o8 C3 k
NASA GODDARD SPACE FLIGHT CENTER 483 5.679
# \6 Y* J- g& ~2 P+ g6 NLABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.1621 \8 b0 m* `! \$ X1 o
NASA JET PROPULSION LABORATORY 431 5.068
* }/ ?. \" m$ M( A0 N( D9 E+ @...国家数据： USA 3118
2 @$ M: w6 S0 }! { ~CHINA 14934 {5 d9 y7 ~8 C5 q
FRANCE 1385' m& D3 A7 R& w- R
ENGLAND 8046 f* P& i& m: t' ^. B5 p) ^% @, @
GERMANY 748
: D6 ?, O Q9 E8 I2 @% WAUSTRALIA 421
$ X" y" F$ m8 _( b) B# a. O- QITALY 4213 ]6 t2 B9 K# N
SPAIN 371, o. u( K* L; d. f+ S% K8 r: t
NETHERLANDS 340
K/ q' \; z: x6 E2 {$ \. DCANADA 323
9 l3 r; n5 I' _2 Y8 G8 i$ P s2 z" i N3 Z' i% _2 g
从地名到坐标 4 \% K7 a4 v3 r% p& A! E( H' H
使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati
3 H- X, {$ s, t2 f3 U# v N* {$ e4 INATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.8830569 A4 P4 X: `- q6 v
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037
' a+ q9 _* D5 U3 `5 G! [# [( \INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.9
0 A% |/ E) w2 Q% L3 G& CCHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327
0 o$ U- Y0 c; p O( b' m+ \UNIVERSITE DE TOULOUSE 1.4 43.6
$ r- X! k0 O: K7 [- U3 cUNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586- o5 F) Z+ ]7 g3 I: P" \# i1 k
NASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495
$ e/ x1 J4 L( ?, @9 SLABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65( w! O3 ? l! p7 T0 i% ^
NASA JET PROPULSION LABORATORY -118.1713944 34.2013081# Z* i( ]' s) Z o
6 W7 k# X t# X8 g5 ]
! B( B- h m4 C; Y0 l( w% D n: |绘图 , U C$ h! j/ n7 z
数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows) d# v4 ?* {8 z6 [
+ g$ v9 }5 E3 I$ X* M+ P
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
2 E8 X" \* g$ c4 k7 S$ A4 u2 Cgmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
( x2 ^# w2 w; B- X- H
5 ^' W, u# M/ }1 ^set ps=3d.ps! r, p+ F" K) r: A
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%
' K6 i+ K9 U- J7 L8 \$ t$ O2 Igawk "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%( ]4 Y* ]8 l1 ^& u& f
. D0 i- K& _* r5 r8 H. Wgmt 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 N1 o, }- v) b) \
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%
/ L6 g- ]. k8 U2 N5 u' N! Q1 qgmt psconvert %ps% -A -P -Tf
) Y, y# [6 G# q9 T* o: d4 P" h
/ C$ v' O# o4 o; D' U. [结果 4 T q: `6 f9 h0 f3 l+ v
图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。
5 t4 |) r1 w, C7 D: J
附赠平面图绘制代码 \0 P& s; ^/ M% X3 }. |
: v' j% \ A: x N
7 J) h5 ^) f7 T) ]- b* K" d+ @/ a& g1 \: \
REM For windows
! ?. |, [' h" O
7 q5 Z' w! C$ r, i3 a; r9 ?gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
6 L& ?# {4 K+ e! Xgmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black 2 Z7 b/ o$ ]) ?9 m8 Z
[" C; U, v% ?2 u3 `- Rset ps=fig1.ps* Y: k& f2 k5 g. h8 L
gmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%
* r! \7 G- w0 g8 e7 [" agawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%
@, I; M/ a, H' a" j5 [gawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%" H0 a/ @- u+ _& ^) c) e
, F1 ^6 V7 v; N) r% F
gmt psconvert %ps% -A -P -Tg8 }8 t1 @6 @; U0 Y7 w
注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。
/ g" m8 T! `5 y' ?
样本数据下载
7 Z' |9 X$ m: ^( b! I; lhttps://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
, {4 K# a. g# \* K' @: q. }& Q' d; Z# S, ^
6 r7 j c" p7 Y! g `1 o) [ |
