本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。, Q, w- u+ M; c/ r ]: j0 c
WoS数据 1 E6 F. b2 t3 t0 }+ N' U
根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,5059 k5 I: O Y" W5 W' U- b' k5 S- E) h
NATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747
1 Y" B4 S5 @. E2 q" t# G! o9 e) vCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735
" E t6 f+ W n' LINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.3137 D/ o) q+ T/ @$ c. F$ d
CHINESE ACADEMY OF SCIENCES 629 7.396, [8 X9 V$ D3 Z: u: p/ {6 M
UNIVERSITE DE TOULOUSE 517 6.079
* N: ?' d' K$ ~9 ~- E! ]5 eUNIVERSITE TOULOUSE III PAUL SABATIER 498 5.855, ~! M! ]* |* R3 C; Z% `5 z
NASA GODDARD SPACE FLIGHT CENTER 483 5.679
, x+ j; s' R' G. FLABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162
/ |* B7 K) ?- ]$ I2 \3 _: _NASA JET PROPULSION LABORATORY 431 5.068
3 w4 }/ g3 b; Z4 F# Q" P9 B...国家数据： USA 31187 H3 z5 D. _0 }, R1 ]
CHINA 1493, t$ b5 Z" ?! ?" b7 V- o) U
FRANCE 1385' n0 \, n. w# s: @
ENGLAND 804
9 [1 T4 j7 ^. u- f0 hGERMANY 748
; H: A W1 N; f+ nAUSTRALIA 421
! d! E5 k# I1 g/ h. g, IITALY 421
$ M( L% V8 @- f5 h4 aSPAIN 371
$ |# E0 K! Y f1 XNETHERLANDS 340
* }' |4 M; @8 Q L0 C6 ^6 @* r: dCANADA 3230 O, Q8 K. b) _2 A7 e
; _1 b# H- r8 ~6 P
从地名到坐标 7 E" x8 t' n, q. [: U# Y9 _$ S" d
使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati g& W; X+ [4 e) \: t$ d0 l
NATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.883056
, J- B4 W, S# `' ACENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.84760370 B/ }' L* O6 d* u% b
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.9
" j7 N4 k5 i" D2 F0 XCHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327+ [. u3 K' a: b
UNIVERSITE DE TOULOUSE 1.4 43.6; V% r5 H+ C5 t# }+ U
UNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586- {0 m3 u: w* d1 y2 U
NASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495# U/ t$ B$ a4 c s: c
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.655 f$ g/ F- [% e# X: M4 u5 L
NASA JET PROPULSION LABORATORY -118.1713944 34.2013081
6 c! v2 T3 v) m) _% _) x! E
7 `# d' V5 u; m- V" d
, p7 _8 ~+ s1 J+ _% u绘图
' r+ Z* o, I) e数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows
: D& Z) g2 M! P" G: O$ S5 i/ X3 H3 M
: K( ]* D# \4 \gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
! w1 L& r: j( B" v: agmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
7 g; S) |3 X- t& z: k
3 y. U, N$ w& s% N( ?$ mset ps=3d.ps
' }0 Z3 c# G) K0 j9 t. [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%
! P, l# }( w$ P- _: ^- z) l& Hgawk "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%
7 r/ S: F( C8 h t2 H; _, p+ {* s' O4 Y {
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%
4 z* H, j4 @, m M2 dgawk "{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 g' I8 h+ P( o% A8 Dgmt psconvert %ps% -A -P -Tf" ~7 U$ g7 V0 D( e- h
: g0 u2 _* x$ C2 q/ |$ d
结果
% y- ` `# H. S3 ?图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。 4 T8 ^- u1 e* \2 y
附赠平面图绘制代码
/ j- ^: D+ z3 _ V7 G. v) T5 p% Y) y" \' f2 r! a
/ k( t& O7 W' k$ Q
( W7 O. H* H9 ?) y4 B9 Y
REM For windows
* a$ I0 d% k3 u$ [+ C/ T1 b, Q6 ^1 R. X5 [7 ]
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
B& n$ p" W0 \* r) J8 ?gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
& _* O% C# ~3 K: k+ P
4 ^# b1 m+ e/ B9 ?4 a3 zset ps=fig1.ps
8 J+ @4 Z' h s) O, V4 d' |gmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%
7 E/ o! T* H+ w9 @1 I+ Ygawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%+ L& ? s% T9 A- F
gawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%
. i1 e6 p( O; \- t3 [8 R
- W. W5 _6 n6 `gmt psconvert %ps% -A -P -Tg5 M+ t: s& q( C( T# Q6 G3 y
注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。 & @7 k% K- |( I. [4 _: n2 m8 B
样本数据下载
: H1 X" Y8 R* P0 H0 Fhttps://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
1 Z0 w& B7 k, M5 z- B/ d) ?" L4 V3 { @& p" C2 }
|
