本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。6 b6 f* m6 K1 Q" ] P8 I( T; E
WoS数据
" {' |9 ~8 }/ H7 |: g' z! f' Y& Z根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,5056 z1 X& u# s2 g. o- ~$ t& r
NATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747
' {6 x' e4 i4 U |7 G6 cCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735' v6 m" y% [% h5 ]
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.313
( C1 E6 B% p- I2 lCHINESE ACADEMY OF SCIENCES 629 7.396: T" `. h I/ @! p1 j
UNIVERSITE DE TOULOUSE 517 6.079' |! a2 m; s( C: W4 W0 K! l
UNIVERSITE TOULOUSE III PAUL SABATIER 498 5.855
9 e0 ?3 }7 |' G ~ A9 SNASA GODDARD SPACE FLIGHT CENTER 483 5.679& y% l1 x/ i9 Y
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162
0 R- |- H+ Y! w- g" b( B( PNASA JET PROPULSION LABORATORY 431 5.0688 a0 c r( U- E- j2 @
...国家数据： USA 31180 V# X8 @9 @, t2 k2 N+ e
CHINA 1493
0 F! g2 _( i) X; ^, ~FRANCE 13854 j, ~" w% G2 u$ p. V/ Q
ENGLAND 8045 ~& I" x+ }7 x7 v1 L1 E
GERMANY 748
4 Q1 d: G6 N8 y2 ?( e" m4 rAUSTRALIA 421# P/ T2 p. Q& A
ITALY 421) J! V5 N& y$ I4 f
SPAIN 371! `- W9 W& E; z) l2 y
NETHERLANDS 340' }- b; s0 v, h
CANADA 323
, T' x- n( k# r# q" Z# J1 U* ]6 n: W6 K4 j
从地名到坐标
* Y0 U2 p( d3 p# \! ^* ]; N: s使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati s1 L( O/ p" P& z
NATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.8830564 ?1 k: C5 u* a
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037, c4 l5 i9 L- |5 H$ q/ O
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.9
5 N" s/ b0 r" ?, hCHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327
% ^. r h% k i9 ^. n% i" [UNIVERSITE DE TOULOUSE 1.4 43.6
$ n# S! t& R0 C- A3 V2 rUNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586& m6 W" e( R t/ Y3 B8 d
NASA GODDARD SPACE FLIGHT CENTER -76.852344 38.994950 |: M; n" W' \) c# B
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65
( k, ~0 Y$ k% tNASA JET PROPULSION LABORATORY -118.1713944 34.2013081# F% b$ I- O* A$ r4 t2 h
$ Y5 b, L( a' N! S0 t1 l, d1 g+ e
( k$ M+ x+ C( ?" ]
绘图 4 d8 g- h3 K3 Z- t! p* c
数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows7 |: Q0 f" W$ @, a
( i+ u( ]# {' o
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
5 E8 Q' m! c- t3 }gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black 5 \0 u- D' L% D3 V
Y( ~' N& ~, [) i6 pset ps=3d.ps
0 z, u! l" ^* m; ?8 `0 rgmt 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%
4 M& G( E* u4 z x3 r; _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%. C% X8 f- u9 Q; m2 w# Y( R& I
) }" P' r5 l: _4 I: N. x4 Vgmt 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%2 A: X& \# n' Q4 N9 U: v
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%
! C$ b6 y0 ^" ngmt psconvert %ps% -A -P -Tf
/ ^- X: K" x% J9 Y/ e$ e' I: o0 s/ v9 _* F% A. w6 k- k: k7 ^+ r
结果 8 S( K0 o' W) S* H& v& n+ Z
图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。
0 Q' C( h/ N5 W8 G$ c/ \9 Y0 o+ [
附赠平面图绘制代码
0 H/ m- M/ u% f" P: o% ^3 {: W$ `0 V3 v8 r# _( a
) B" R7 l" U! x+ ^/ a1 h# y
- A4 w" D. B& w- ^0 j3 {
REM For windows
3 p) i" |, u$ m7 R2 p$ c, _8 Q f5 e ^
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p+ {$ P) M+ O( C' J* m" a; Q
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black 0 \: F/ c4 V. t+ o2 g
/ a8 R- E" m/ O2 C. e7 V6 _set ps=fig1.ps
* ^! E, z. m8 q5 l( Z( mgmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%
: y2 N7 @) a# n, l* u4 w9 ugawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%/ a& [2 T" e2 Q& p, O" p
gawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%
' w* x1 w4 s) w5 z. f2 V* }3 f) |2 _4 N5 A @
gmt psconvert %ps% -A -P -Tg' m3 J# r; Z3 T! }6 M
注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。
! G7 ^! I! D. _$ U# M
样本数据下载 8 O" d$ S" g; B4 h8 p1 S4 ~
https://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
" N# N0 B5 |$ d: u) n& g& h# c0 `( s- u) W' m
|
