本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。( H' E' `* y9 l* S6 J
WoS数据 4 E: } d+ P% |+ w6 X
根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,505
7 W, N1 c9 d7 ]0 y. F% }; }* HNATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747
5 s$ T8 c( [, T. [CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.7357 i) q0 k# E3 h5 f
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.313" d; y6 W1 j8 X0 {$ \
CHINESE ACADEMY OF SCIENCES 629 7.396
% p7 W, V+ L1 U3 `+ \, y5 YUNIVERSITE DE TOULOUSE 517 6.079
0 y! W% k& P! j }. A% VUNIVERSITE TOULOUSE III PAUL SABATIER 498 5.8553 m! @$ `6 ^7 k. i C8 P: i3 Z
NASA GODDARD SPACE FLIGHT CENTER 483 5.6790 f8 [+ G z* Q6 N5 v0 @; V
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162
( d' p7 I# C4 I8 `" MNASA JET PROPULSION LABORATORY 431 5.068
w. G; h! o ] g3 Y- A! |...国家数据： USA 3118; t8 b* o3 z; I2 i* s6 Z
CHINA 1493
8 y% a" N0 X: x3 y0 r2 Y6 U) q% }FRANCE 1385
" n3 u/ l& S. T- N( j# p( L! qENGLAND 804
% r) r% A! w! d. fGERMANY 748" M0 L4 I7 N8 r( ]1 ?6 C5 l
AUSTRALIA 421
2 f( D6 Q* r$ L8 HITALY 421
; s c! Z( m7 W; b7 B6 {0 qSPAIN 3711 \* |% ?# _! x9 P/ v6 T$ W
NETHERLANDS 340
" X" j* j3 G% \4 b5 X0 |. w- ?5 OCANADA 323
9 h6 l- y3 `/ m) t
' ^: j% h7 P- I8 o从地名到坐标 " j2 ]- V# b8 w' n+ p( D4 x( J
使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati9 M- F3 A! D# h
NATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.883056
' t$ q$ x' W. ]8 Y/ I& I2 h0 A: C) Z2 ACENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037; F( ?7 @3 L5 p, m$ a" t1 {
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.95 S& K$ ?' E! t5 i# m, Y
CHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327$ C% c) I/ ^! e+ t1 B5 x+ e0 A% p
UNIVERSITE DE TOULOUSE 1.4 43.6
6 \& }2 a, m N" @& z3 JUNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586
3 p" F4 n: k+ N0 INASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495& \" g. f& \$ N
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65
+ o/ l! F O+ v0 }# xNASA JET PROPULSION LABORATORY -118.1713944 34.2013081
' J. u9 t/ ] W2 H; [0 P. C3 T) M+ S0 R. r1 o) l% P9 N: a
. P S5 [0 g, S0 ~: x$ t% [
绘图 $ W, c" m, }7 s
数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows
! T( u& t! m9 F4 ?" l
+ f9 t+ `0 y3 s! s. t9 L' Zgmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p* d0 C. M5 P; e( K) I
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
" o: I" U A/ h {: O- x
- r$ D1 i% \: `+ ` mset ps=3d.ps+ H# c/ D; {# i: I
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%
+ U3 D3 q2 b/ P( 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%4 g7 J5 }) ?8 ]/ z' j% Z- i
' N* o- f; p1 s# h3 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%
3 S# B/ P/ S* t3 hgawk "{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%$ T% L: r5 U# Q2 F2 ?- z) v
gmt psconvert %ps% -A -P -Tf& V7 j& P ^9 b+ X) k
5 Q( } X0 | Y7 E
结果 ' j4 r+ U$ J0 G& h! m: u: W1 h
图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。 % ^- F2 F, f% b
附赠平面图绘制代码
5 Q& z) S6 f( {- F, W: @5 J
7 l+ J8 u" [0 j! Q+ {' j2 o& F3 f$ Q7 N
. d$ |2 O+ f6 c( d$ k3 lREM For windows
) X: Q6 j* l9 X( b5 J; M
1 ^1 S4 r( Y# P; }3 Z* Xgmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
3 ~7 H5 x+ w& s3 ugmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black " D0 q# N8 J7 x" e0 J% T, G0 r
# v% e; Q/ m2 k8 ?' X3 o% W, X0 O9 uset ps=fig1.ps
! V) ?! b# T" C5 q5 R) ~+ {* w. R/ Rgmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%
8 B+ z1 Y- `0 Q6 c5 A9 Vgawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%
) [; b F7 W, E& C' S( e0 j* Z6 E, A0 qgawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%
# p: z0 i+ G/ V7 h$ [
: f! L- w2 Q7 h3 O5 J# ^+ g7 ygmt psconvert %ps% -A -P -Tg
6 o& U5 q. @/ b! n+ `注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。
) E$ \, B( H: I" M: [
样本数据下载
0 O0 t2 `& K$ ^6 S3 Khttps://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
: u) ~0 ~; t* i7 o0 [5 X, d3 y7 o/ a
|
