本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。
$ h0 r; B- N- ~* ]WoS数据 * B6 N/ B9 ?5 h7 O5 F
根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,505
0 _; u4 ?. m5 ^, n, p" W0 QNATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747
' S5 b. N0 a1 L. b1 }; c; m' ACENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.7359 R m0 ?" w- x
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.313
' E% s2 r3 O3 z. o) |CHINESE ACADEMY OF SCIENCES 629 7.396
( L3 J) m# ? N* F: @UNIVERSITE DE TOULOUSE 517 6.079
- [, L! k, g! L8 pUNIVERSITE TOULOUSE III PAUL SABATIER 498 5.8559 K' }$ n; K* v! x! F. u: D% d
NASA GODDARD SPACE FLIGHT CENTER 483 5.679) t: t. I( V5 M' ~1 [* s
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162
/ c# j( a% n9 L' M2 |4 x$ bNASA JET PROPULSION LABORATORY 431 5.068% @8 ^8 _+ X( x
...国家数据： USA 31180 C& ?% k. B' l1 d. I" L
CHINA 1493
2 i8 D3 k3 r7 E, W6 j$ FFRANCE 1385
. m- b! ^' y) k; T& WENGLAND 8044 E1 d! X# q& R2 }$ V
GERMANY 7480 S# R: M+ T2 s3 z; L4 T: m
AUSTRALIA 421
' @4 p2 s& e; v! I* C) Q. F% e( @ITALY 421
- K9 d2 P' k# z* c& Q; ASPAIN 371
' P) x" I$ C- z+ Y$ d; FNETHERLANDS 340
6 q- _' }$ U* o( RCANADA 323" d: F" U: B# X' F; o
) e) `0 n6 I' f: G$ F/ W5 V; C从地名到坐标
% \8 O. d% P6 M. Q% c' p: G使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati
( o- H3 T% G' @* B5 TNATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.883056
' u. J6 W* m' b/ t2 ]% `1 @- ZCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037* `; O7 W$ o' }' V G/ ]) O
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.93 Q- e b; e" v; q+ p
CHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327
" h& j; u% {! r9 C3 ?( `0 L$ NUNIVERSITE DE TOULOUSE 1.4 43.6
; ] n3 V8 J( h+ }" eUNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586
" r/ H" D2 Q$ `. yNASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495
; Q+ C- R# z. GLABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65; ]6 f. G3 R- l: C
NASA JET PROPULSION LABORATORY -118.1713944 34.2013081
. t9 O7 o( D8 r$ S1 Z8 f
1 X+ w) {3 I* k8 P- x! ?1 f% K. P" ?: r9 h( l3 `
绘图
3 _% u Q5 u, o数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows* n ?. G1 }2 }2 s
" L3 O' M; `. X r' [$ J, u3 m6 Wgmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p/ k' U" k! U9 U
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
" K5 x( z, U0 ^% K
: m1 b' A3 e8 t8 oset ps=3d.ps. ?4 L1 I4 C+ |5 E% u
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%
6 f2 H5 p, ]# }( u, Ygawk "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%, N3 M3 z# Q* a. ^+ l2 t3 B" c
+ p- k. z# n5 v1 @+ Z0 `! q3 ^
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%" z" c+ \! w: x) ]* c" w) s
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%
# ?. W$ t! P T: A( @gmt psconvert %ps% -A -P -Tf
) P' _/ g! j: S& o# p- y9 |& _$ L% w' \: }2 s& f- o
结果 . G y5 P# G: O: G/ Q! s; e
图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。 & y# d: [ m7 \$ p9 [
附赠平面图绘制代码 & b3 P8 `8 D% r# c& Q: a
% f9 i. D9 ^; n; f; {
/ y2 S* ~$ a) R% H0 u, E% t
9 @$ h _8 G! G" W" @* SREM For windows
0 M% u7 A$ r" } s5 T$ J# O6 p
) M* m+ L6 q2 x; agmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
7 w/ n6 k' g1 |9 M% bgmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black 0 C, J1 ?+ g# @) X9 @( C# N) }
- j9 ^9 k- x- K. Uset ps=fig1.ps) f8 B! U, |1 S
gmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%% L6 ?6 S% [5 x F+ S
gawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%/ ]& Y" ^' J0 s1 T l3 i4 L: }+ j
gawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%* U+ `/ [2 Y$ F% n0 h# t! _2 n3 Q: E$ z
& A! ?# F2 v/ T# r4 H+ t
gmt psconvert %ps% -A -P -Tg
% h* I0 }, v8 I% q( M9 j* G注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。 1 W% s n8 b) w1 F" f
样本数据下载
" O" j" s; g' Hhttps://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
+ @$ y' ~; }# Z) y7 {
& q5 y" Q+ @5 y' D/ b/ w! f
|
