本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。
: J. z. W; @% Z6 AWoS数据 ! K9 W6 `4 N5 V. w+ u5 w) K: R& V. s
根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,505
. ?% k" N$ A7 aNATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747* `/ y" a0 r. k& O
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735( `: J' n6 _! I) {8 t
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.313
+ x9 t# t& \# f8 SCHINESE ACADEMY OF SCIENCES 629 7.396$ _7 A k+ p. M0 o4 w$ p, s
UNIVERSITE DE TOULOUSE 517 6.079
3 d/ G( H1 d5 P3 L& G! nUNIVERSITE TOULOUSE III PAUL SABATIER 498 5.855
( G: r1 c5 H- B, u# TNASA GODDARD SPACE FLIGHT CENTER 483 5.679' ]' I5 f' z9 _6 @0 C
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162/ e$ ]8 d2 L" t3 b8 h# q9 B
NASA JET PROPULSION LABORATORY 431 5.068, W \2 {4 X8 Z
...国家数据： USA 3118
4 q' V* S1 g, V6 J8 `CHINA 14931 ]! b, j( u$ v* q2 `# t2 ~3 C
FRANCE 1385- l u) h$ x% O9 c
ENGLAND 804$ c' a, f9 I. T. q* U* u9 G
GERMANY 748
5 |9 |! P# T8 p+ m4 ~* y uAUSTRALIA 421
! Z6 L0 _# R; z- |8 eITALY 421 N: y# C9 Z- i4 }2 l/ {3 o
SPAIN 3719 W- l6 E) o0 Q& X5 l
NETHERLANDS 3401 T/ C" \& L! q. z# V1 E
CANADA 323
- l! c! x* n# S. B3 C- L
4 F9 e3 K# f* C# l5 {4 K从地名到坐标 6 z9 p: p9 n' [8 f, w4 p
使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati
2 d* X9 {3 C% s, q4 Y0 XNATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.883056
3 |* O$ p# j5 y% ^9 tCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037
Z8 q. W3 q: Z( yINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.9, m, F# `# l" C* Y' ~' n |$ @
CHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327
! H# q) g9 V+ sUNIVERSITE DE TOULOUSE 1.4 43.6
( n2 ]) d3 a. }; AUNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586
+ e# p4 n. x" U; K- ZNASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495
5 ^" {8 B. G. [6 w# N' fLABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65% `( [" W3 k9 x7 _
NASA JET PROPULSION LABORATORY -118.1713944 34.2013081
+ ?& |3 T$ d$ [9 L9 a4 W/ E3 b+ B- o/ z) p+ O5 f
6 {# t5 @+ [- b. j5 W% W- U3 @. a绘图
+ G8 M$ U Z! B, o0 ]2 |数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows1 t& M3 S6 l# r: a
# v$ r4 _: s- v2 r6 B! K" ?: g y- q4 a/ Vgmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p2 N: ?7 ~, R! g- S) s' j
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black ) G* ~+ c6 q" m) [
" c7 q! X; Y: @. Y) P' c. {set ps=3d.ps$ \! {7 U2 m2 n4 X
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%$ G; D/ T. ]( W: r+ `
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%
2 p5 a1 B9 K% O& ~- r& H) |) ]/ j* C$ }5 ~
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%6 ]* Z/ X; u+ o- X" P5 D& E
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%
# H" V2 F0 ^, E, V7 M1 ?! V4 |5 b6 X+ fgmt psconvert %ps% -A -P -Tf
2 S' j- G! G9 y/ T" [+ D: j F; T( `+ ?& K. Z" t- d
结果 - _' M( ]) }. h: C
图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。
1 Y2 H" u4 ~. Z) X3 D
附赠平面图绘制代码
3 F$ S4 }' F; n* ]' Z* V& I1 M, [+ u% V$ o
( j3 N/ u" _9 a
. a$ B$ I" n6 e# b3 E9 `8 \
REM For windows e1 L, J: B' R- M
. v1 `; j) }6 X v4 q
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p9 y C& d: Y% V/ s7 U
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black . ^% d2 o1 N" u# D
4 g$ A1 \& ?* S3 R% x# Lset ps=fig1.ps
- S q( @/ j% Zgmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%
$ e) t7 }3 r9 G5 r3 Cgawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%4 N; ]: d0 Q: c. V9 ~$ f
gawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%% y" t; m O% l: P
; D4 r( J+ ^% I
gmt psconvert %ps% -A -P -Tg9 J7 I% o- s, a% u4 D; Z
注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。 : z1 o5 W( C N4 K& D& k
样本数据下载
( @" a- _1 w) \4 q1 e9 z1 g1 _https://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
+ z: W, ]1 j+ a
# S: t( I2 ^. S: ~4 w$ w% i# R0 h) N2 N
|
