本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。7 ^" {3 n2 P C! J! V/ L$ N
WoS数据 . r; x3 L4 F$ R$ {# w
根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,505
' ?: I0 f2 y9 ]- r9 lNATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747
. W7 m: J1 |$ ?7 ^1 yCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735
+ E. w1 v3 c6 [ u* I- nINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.313 W. L8 Y8 a$ Q+ ]2 q* b
CHINESE ACADEMY OF SCIENCES 629 7.3968 n9 b1 Q" N# Q- _3 f
UNIVERSITE DE TOULOUSE 517 6.079
& C0 ~- X. |; Q# S1 s3 ZUNIVERSITE TOULOUSE III PAUL SABATIER 498 5.855
6 m, J6 e9 v3 R" ENASA GODDARD SPACE FLIGHT CENTER 483 5.6794 a4 B3 K4 V: e0 S, f0 ]7 p8 W
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162
% L$ e' C0 C' s# oNASA JET PROPULSION LABORATORY 431 5.068
! s' N- ?0 S# t# } ~4 q8 j5 `9 h0 f...国家数据： USA 3118
0 \' s" ?+ Y& [* e0 F2 O4 B! tCHINA 1493
! D! k1 V+ t; C& g7 ~FRANCE 1385
% b- ?% ~1 G$ g! X- t7 YENGLAND 8047 }& \4 [7 K( |+ S4 B, p
GERMANY 748
3 u" Z- v- k1 Y O' P, R" c! AAUSTRALIA 421% r: P# ]3 h) m' Q, L0 Q: C
ITALY 4212 O8 e; r& `8 r9 v. j
SPAIN 371
, z* {/ G6 s' U& }! F# bNETHERLANDS 340$ a& ^, I6 l0 z, `
CANADA 323
( \5 G T# V2 N( P. ?6 P; [, _
6 z. j% K! A+ X' w2 N" @* m7 J+ v从地名到坐标 5 U Y6 L& w! j# S' m7 B: {4 i
使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati9 `- e) t- p/ I
NATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.883056
8 v) p$ V" o! I, nCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037
7 ]: j/ G: E1 ]- V; xINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.9
0 M* D7 v* N# r+ Z7 P9 ?- T5 FCHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327 | u2 I& N4 G$ d* _! G5 z, |
UNIVERSITE DE TOULOUSE 1.4 43.6
+ O2 {9 K$ u1 oUNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586( I" E! V! P* D, I, V1 U
NASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495/ O, k( M8 g2 k, ?
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.656 v7 z# o& ]7 i7 C8 Q6 p
NASA JET PROPULSION LABORATORY -118.1713944 34.2013081 p( A$ O/ I: {1 u9 ~' S. n
! j V( K% o2 ^; R8 ?! |$ U8 \
5 M9 Z5 ?5 E7 K- S& q% W2 x6 J绘图 6 ^8 e3 T! B) s: q+ B
数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows6 s. \' e/ A: E2 D3 C; Y
" [0 w7 }$ Q$ A1 |' t1 ~9 ?gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
6 \$ `2 K& P- Z" P# vgmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black 8 Z" m! t1 v% O
; s: [ I4 x( X7 }set ps=3d.ps
. H" A" M' z. ^# S# Pgmt 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%8 r, n) K$ O2 C7 F4 V+ {5 I
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%/ K7 w5 E; M! z+ P3 U0 b" u' B3 g( f
! z9 d1 D/ `" N- Z8 qgmt 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%: n) Y" L/ B7 o
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%
# ^& Q( A0 a+ i. r2 l- rgmt psconvert %ps% -A -P -Tf
+ N2 y. n) ^" t/ o$ S& f& s
7 _8 v- h; }, V& z8 {! c结果 ; E% ` Q) _/ w2 i: }' y
图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。 ( n1 ~% f: L" v2 `: g
附赠平面图绘制代码 9 _ h0 l( M) j4 h$ V9 y; h* ]; V" P
( N; d0 L8 f8 q; H0 B" I
: w& S$ C6 Q" k6 i, l* r' f
4 P. ?( G+ i2 i4 t. `REM For windows
- U( K/ N2 S6 Y* B. J) @2 U5 \% Z, g+ l7 x) B
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
- A* p3 B% ?% Ogmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
4 N) ~4 r) r) J' l- A3 T0 K1 s: Y* g4 m5 j7 O
set ps=fig1.ps
/ X. g' o+ r# u* z3 v: w3 agmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%
8 a8 E% ]; b6 P& _* P, g1 wgawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%0 \4 G' o5 ?( X$ u0 ?
gawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%
" U( {9 L+ ?% X0 D8 v
& {+ W) o* S* d) J3 tgmt psconvert %ps% -A -P -Tg
& d" |: `$ i* B注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。
& `. e( ?$ i1 K w3 N/ ~
样本数据下载
4 }. \4 @5 A" k! R6 V* g: zhttps://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
b# l% C R1 ]8 b* [& `
0 ~6 S" S7 S6 g" z$ c% z: y' J$ f |
