本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。4 a+ q( x% Z! h0 E+ m
WoS数据 + t7 w9 R, Y" H- ?
根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,505% a B! W6 }( R! c. ~7 T3 d3 T! x. c5 h
NATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747
& w) l' V2 Y; ^8 qCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735
4 g- a4 X) j# o; ]2 Y9 ^! dINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.313 ? s7 U: p4 a' d. P/ z/ Q
CHINESE ACADEMY OF SCIENCES 629 7.396
, C! D3 @* i ?! q9 K6 FUNIVERSITE DE TOULOUSE 517 6.079( |9 F5 H/ c- q y. Z; I9 a1 k
UNIVERSITE TOULOUSE III PAUL SABATIER 498 5.855
. M k* i0 I5 {& B3 PNASA GODDARD SPACE FLIGHT CENTER 483 5.6795 g0 i* F3 z' d5 }3 [8 V% O
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162
! q+ X6 r6 y+ `9 ^/ ^3 wNASA JET PROPULSION LABORATORY 431 5.0688 R- Y& j* q; d: Q$ o
...国家数据： USA 3118# y9 X4 c7 {, C2 e+ a
CHINA 1493% c1 v6 x4 D- ], a, Y& X: U/ y. m
FRANCE 1385, V7 ~ c9 q2 \0 S! y( Q. ? @" [& a
ENGLAND 8040 u' j l. O' K$ L" p3 b
GERMANY 748. o( N: `& Q" }
AUSTRALIA 421
! A6 m* ~, @; ~2 L/ A* gITALY 421/ H! ]# l8 g7 c- O/ b6 ~3 ^
SPAIN 371
( U0 J+ M: Q8 j; A3 e, I6 ?2 p; F* pNETHERLANDS 340
1 Y% g) ^$ h. |& L2 ECANADA 323
1 \, T. P6 T9 ]3 ^( O7 X% B- w0 @ Z6 K/ @" r H5 t
从地名到坐标 : H z; P, v& f6 S% \" v, i/ y/ N
使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati
; N5 i$ A) A: s2 C$ ANATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.8830569 a7 x8 z6 B5 n/ w, s0 s" p1 r2 M
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037
1 l( |8 A7 |7 }9 B+ dINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.9
e5 h$ Q$ X. A0 UCHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327
3 u [$ y2 U, U& U9 n* D; z0 qUNIVERSITE DE TOULOUSE 1.4 43.6, k" t' O g+ j9 @
UNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586* N/ q+ h7 k5 _5 ^2 k
NASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495
& B/ A! U: x4 k" f0 L; TLABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.652 _: a( ?1 ]4 j5 B
NASA JET PROPULSION LABORATORY -118.1713944 34.2013081: e9 S! X- V9 _
0 v1 a. Q2 Q1 j8 ^6 e7 _; c9 w' E9 _8 ^9 z
绘图
3 T7 \2 B" g/ D# E# u( J% K7 f数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows
/ r5 C3 P) f( r. S$ x2 k2 M! ^3 z6 k: ]8 F, ?* p! E
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
{5 g. Q. b" S# q- v/ Y% kgmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
( A& b1 Q+ \4 ^
# N1 `& E$ P) Y( O& pset ps=3d.ps, p/ Y6 L" n; }8 n! r
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%
: w) `: Q- V) x& S/ M$ u0 P( igawk "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%: ^6 u# W- A; U+ l+ M
% u# p! G+ p! _5 b3 c
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%
: i1 D `9 F3 x$ w Y) I" ugawk "{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%
( |/ U/ U5 @8 ygmt psconvert %ps% -A -P -Tf
# Y4 d9 J) n& B6 @3 f3 C$ S
4 D1 f; u0 x$ c; n0 L7 M结果
( ~ U: V; ^$ j) b1 F图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。
8 N( T9 d+ E3 U# @
附赠平面图绘制代码
3 _. G+ E, u) i. d6 J9 U7 r4 p" [) n- Y! ~. A- Y! `( A. Q
3 ~6 o2 J' B5 w) a6 Q) Y& C! t5 h! R* C; i. `7 ?
REM For windows4 N) O! Y# n8 B* r. k: Y
# y# o5 y; n/ Y5 ~$ U
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
$ B, o; X9 v+ w" Q% agmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black 4 T8 S, c$ M7 C. A$ C: F
, Q6 O6 b' W/ \# z. yset ps=fig1.ps4 B* l! e$ T! ?. l
gmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%. `& s$ M }) Z
gawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%
8 F) n3 l9 v9 T4 Agawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%3 g8 [9 q8 D3 ~6 X4 Z, J/ f
; \$ Z. a* k. m: F
gmt psconvert %ps% -A -P -Tg7 ^1 k* j, Y4 T, i
注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。 , q2 k/ a4 y9 j* v' N& r
样本数据下载 + r. X' N4 l" u. h7 P
https://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
) b3 S+ R) P/ m
' O, j" F9 ?! Z, v% k0 u8 H |
