本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。
- ?" P. T2 N' UWoS数据
7 P1 v# `3 j4 W6 }& A根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,505' J$ T! E* ], C6 j, @
NATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747% ^; Z @- I; A. e
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.7354 [7 V& Z: D( r! {" F
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.313
u1 N* _' l$ E" B9 s5 t6 LCHINESE ACADEMY OF SCIENCES 629 7.396- w, _9 s F2 o: D+ d* `
UNIVERSITE DE TOULOUSE 517 6.0792 t3 F; V. \9 X' U4 y
UNIVERSITE TOULOUSE III PAUL SABATIER 498 5.855
* Y) p$ p0 r8 |2 C4 `NASA GODDARD SPACE FLIGHT CENTER 483 5.679
7 V; P4 k# X$ ?+ g, _/ G( eLABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162: x! @/ w1 b, R0 n
NASA JET PROPULSION LABORATORY 431 5.0681 o) [* t- B! N2 b# N5 h% q
...国家数据： USA 3118
8 l. n( L$ d9 y: m# J( f4 [% XCHINA 1493
% a/ J5 ^/ J5 XFRANCE 13850 q0 ~( i' c6 H) r
ENGLAND 8041 e+ R- a& w: o2 _
GERMANY 748
7 L" ?5 r# e1 s$ I% D' ] TAUSTRALIA 421( z. s9 W a; E
ITALY 421
) Z4 Q; T+ j1 iSPAIN 3711 Q5 \2 s) S; `- ~8 |1 U1 h$ ?
NETHERLANDS 340
( l" @2 S: M, p4 n" Y6 R# ZCANADA 323
. U$ n: ?: F7 O3 c9 y( d4 v8 _, J
从地名到坐标
: P8 ?% s+ K( a6 h7 M% r使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati3 j3 c. H% b' \
NATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.883056
2 f! y- v s) d, v- r) ~6 \% wCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037
6 @+ D1 p0 M7 K8 K! M4 iINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.9
4 C: w( u+ {6 w% y( C ACHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327
/ d( i' n1 f9 [0 o- p6 AUNIVERSITE DE TOULOUSE 1.4 43.6+ d4 B- ~8 w0 W7 S: \9 ^
UNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586
. m% d* \" } ZNASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495
/ J4 o+ ~' s: O7 y2 dLABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65
5 X3 w# a& \$ k4 `- ]9 |NASA JET PROPULSION LABORATORY -118.1713944 34.20130810 M8 u4 B! _3 `8 z/ q, |+ k# Y/ \( w
& R7 W& r% y: m8 p
# J* u8 U8 y9 ?6 s1 y3 m9 H/ K
绘图 6 G$ t; o; S( e+ V' K" J+ i" k
数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows U- r P& X+ U# A. f8 l) R
, q7 T7 Q2 R, \8 c! m6 ~ r9 Dgmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p) E8 N9 G: G) x9 f! i& ~1 n$ a! [2 h
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black A, O" V5 L5 f% ]
; ~9 I4 y, L) M0 v) t' \% Gset ps=3d.ps
) F* J5 l8 V7 P% X2 R) e0 h) ugmt 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%; t2 Z; ^# ]% ^5 `% x) N, k/ y
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%) m8 g: t& S% ? ~' n7 r
) q7 }2 S6 v; N7 C6 X6 p! z
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%
G" ?4 \& j" y' D* jgawk "{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%: K3 y6 N# V% }' J2 {/ B. z
gmt psconvert %ps% -A -P -Tf
; O) h+ r* @& }* V7 [/ b" Y
* e. Z* ]+ d! k; u3 N结果 $ K# t0 j1 A( j) S0 X
图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。 9 \5 ]1 u& N* S4 j& ~7 p
附赠平面图绘制代码
8 ^+ ~2 [ p% y8 o! x1 S& d. X, T3 w7 u$ S" t( j
3 r2 D ^! u: W2 X$ v9 j
9 i* E: ^6 G# G* I% d# U2 W8 dREM For windows
1 I. N5 t9 y- ^2 r/ q( q
/ H. l6 y: A) f1 wgmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p* x; d! {# b( a3 z0 f, g( c: ?
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
; I5 i5 k9 ?: g9 s
5 B+ F0 o4 S2 U( X/ _, p" ]! Zset ps=fig1.ps
2 U' ?: n0 {* dgmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%3 U$ W& C0 N7 A# Z0 ~7 z
gawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%
! r( W. c9 N+ l8 k$ n$ Y3 dgawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%6 f" y7 N( r$ O; F
% p& |1 S8 ?' j* \$ g
gmt psconvert %ps% -A -P -Tg
( I9 \& o6 k- E. D$ S/ c# n注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。 ) ]' K! V P+ D( e
样本数据下载 , N9 J! [! F) r& |2 @: G7 j+ Z- X
https://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
8 k& i: V0 J9 `( j7 u$ ] ~$ I
# D( X1 ^8 y+ z" y) @) S |
