本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。
7 g8 z& G9 O, j. k5 y, H/ oWoS数据 / N' f9 X4 `( ]
根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,505
4 |0 ]- Z% r* n; ENATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747
. V0 }: r8 C. M. Q* o1 `4 Y. wCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735% e# p8 Y) C9 w8 K6 v
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.3136 w" o( m+ a! ]
CHINESE ACADEMY OF SCIENCES 629 7.3968 _/ U; g) u' }8 I, D6 m
UNIVERSITE DE TOULOUSE 517 6.079' p! E+ X# h# v2 K& z) T
UNIVERSITE TOULOUSE III PAUL SABATIER 498 5.855
5 F- f) U- `4 S) L$ Y/ ~ p5 sNASA GODDARD SPACE FLIGHT CENTER 483 5.6796 ]* p* v' @+ |. m$ b' U$ x
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162
$ s- K2 B/ q1 SNASA JET PROPULSION LABORATORY 431 5.068
) q- L( v) }8 {# [! o N% O M...国家数据： USA 3118. W {9 }9 k5 m; x( z# [
CHINA 1493# K: }" U" t; S% j
FRANCE 1385
% i# _ A1 H @8 `/ ^1 ZENGLAND 804/ W% e z9 A- B) Y
GERMANY 7482 C% ]* R, r6 f7 {' q) x1 u
AUSTRALIA 421
4 G+ z( u7 @' R1 R3 L& WITALY 421) g8 |; B8 M" @* w E4 l t* Y
SPAIN 371
4 u1 l* d, D1 |7 e) r' jNETHERLANDS 3402 ^* ^* O# h; v8 c, W' U
CANADA 323
: {' A" V6 _% g1 z; ^0 x3 [9 T" p0 H3 t, A
从地名到坐标
! r& ?, g7 s Z; v6 Q' k+ A: k使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati. Z" [' I; B9 X e
NATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.8830561 g* f+ R% p* ?! K
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037
; |# I8 k( q. |$ \" Y/ z8 r; i, GINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.90 F, r% {; e, o$ i- k: Y( W ~/ L
CHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327
6 p& ]# Z! \0 W7 r3 RUNIVERSITE DE TOULOUSE 1.4 43.6
) _1 F1 h2 L% F2 YUNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586
4 |/ {' H3 s6 @8 @ dNASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495
7 A1 L7 B" L2 g1 kLABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65
% U, K5 O6 S! x- J! Q& r4 ]0 ]NASA JET PROPULSION LABORATORY -118.1713944 34.2013081
; v7 J, I" P3 l- u" L
# d; t2 k, b" V& T9 ~0 h. E
- _$ M @6 _" N4 ~. f绘图
9 R7 a7 o! c. ?7 f( d数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows$ T: F6 o/ N, q% D8 H
; F" v7 D1 \$ Egmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p! O2 O) G1 i- i1 n# U0 U
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
8 t) A4 F/ A8 V; o
& g* _ W; P4 Z5 S0 Xset ps=3d.ps( Y; A; j; D0 p* S
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%2 A: w4 }% R4 `& ?% E
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%
" U0 p' k: N. t5 N" u# r K
+ F- y6 a5 C$ B% \: ^4 q( ?, 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%: i' P" O+ c6 J/ 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%, _0 `. ?: @7 M
gmt psconvert %ps% -A -P -Tf
" N' M# L1 ^% m
0 @, Q6 T# V& K结果
. W8 ^% P! k& o4 o6 C图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。
; ?/ n. i! g! W9 `0 H
附赠平面图绘制代码 1 V0 G W; X9 R- e7 K* C
$ z( W' [( W6 k, Q1 _8 ^
6 B2 j x" y; C! F) }7 I
2 C( [% C* ]8 T& x# O! d
REM For windows3 o; P3 u- K3 f
- {! _, \' k. `* n1 Igmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
# y* R$ Y7 `7 x' y. I1 F0 E+ Mgmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black # Y& o' d6 F+ m- \/ Y# X
& J* e o( Y9 g3 f2 ]/ x, k
set ps=fig1.ps+ Z' O% y+ ~8 m2 _; p$ W+ f' A
gmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%4 A, a: P$ o- j' O# l/ ^# W+ Y
gawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%
; I0 h4 W! ]( S+ |4 ]gawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%4 p1 i% r' E8 R% M) j/ d3 y, ]/ [
+ n5 t8 J2 l9 Y0 @9 T5 m* R. J* bgmt psconvert %ps% -A -P -Tg
' P. q1 H. u9 L+ B$ G7 h注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。
% R2 p& [" s. |
样本数据下载
' k; W$ x! G7 [# c" \' K* ~; O5 |+ y) {https://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
, I5 e: [& d6 O" U( R# t
2 h: _( \3 W* F" _% a6 k0 O1 G7 | |
