本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。
; s4 Q1 ~$ V. N0 ]* A. UWoS数据 # h8 b) \, f. ]4 f4 `
根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,505- Q# \* I) n: X& N+ G0 v0 ?
NATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747 b* q5 z" z ?9 M9 {
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735( @8 p" g# M, ?1 q
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.313
/ y: @+ Y3 R4 v/ g0 fCHINESE ACADEMY OF SCIENCES 629 7.3960 _, r1 e/ X+ ?
UNIVERSITE DE TOULOUSE 517 6.079- ]; j/ @( D9 \$ H( o
UNIVERSITE TOULOUSE III PAUL SABATIER 498 5.855; ^- w9 n2 {) W f
NASA GODDARD SPACE FLIGHT CENTER 483 5.679
5 Q+ k N; _6 ^2 H3 P* G8 l2 dLABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162
( J6 W4 z* F! C) E' a$ uNASA JET PROPULSION LABORATORY 431 5.068
( w9 I. G4 C1 a...国家数据： USA 31182 S5 U+ O2 K! t5 n: P0 w
CHINA 1493
6 E2 ?$ _* f, M5 s, M9 |FRANCE 1385
9 {& m: Z. ^2 q) rENGLAND 804* ^7 @& B% b6 D% I
GERMANY 748: u6 b1 J$ g5 F) [- J
AUSTRALIA 421
0 W6 P2 `2 z$ R& {! O* QITALY 4216 H s6 m6 r$ J0 Y, `
SPAIN 3719 y, Y" K8 M! J; k5 {
NETHERLANDS 340- E+ h5 y! {. F
CANADA 323
/ `: c# v1 r& _" L
. O7 J/ l' q' X从地名到坐标 ; ~8 i( J2 S F& h" D& W
使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati% O; R1 V1 r e% V, U+ G, g" ?
NATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.8830566 p g5 t! G1 V7 t- i" j/ M# L+ M
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037
/ j8 x) e4 g8 B8 M+ N- f) rINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.9
$ ^% C3 `8 m7 h9 {9 ~+ |CHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327
( ?% o9 K. [0 H( z& WUNIVERSITE DE TOULOUSE 1.4 43.6. j5 f; ^+ g3 \" x+ z- k' S# b
UNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586; Y P, L9 a; j6 j% K
NASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495; A' M8 |0 |2 l5 e1 }
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65
" Y5 k3 B; G/ r9 Z7 F; X( eNASA JET PROPULSION LABORATORY -118.1713944 34.20130815 w: T3 p; r% ?6 T8 U
$ R( B C/ `9 X z4 D/ k/ l- z
! E5 [2 [5 }, ^绘图 9 U# A N+ e3 E% ]
数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows
6 f/ r+ ^) ]# n$ {$ U, C
# |# c" M* R, Ggmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p& C$ n& [# T( P
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
- c6 J+ g7 o& d; N7 L" i, o5 d0 L* p
set ps=3d.ps
' I. w, y' E h' zgmt 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%
2 B3 X" k% Y0 H s. k: H$ Tgawk "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%
* q" K7 o% [; o& J
" t$ Y+ R" ]" w6 ~7 ?7 Ggmt 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%; {* j! n( j* |+ e! N
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%0 Y# Q7 p. B; P$ [! Z* p+ W% l
gmt psconvert %ps% -A -P -Tf. s; k4 s+ }" ?# p9 E
& t3 y3 _$ n2 N2 M4 B" l0 d结果
. [4 m' v6 f2 a6 i" j. d$ d图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。 / o! U# Z0 z* p1 {
附赠平面图绘制代码
: s8 `5 T! n' F2 M: C8 _' W( ~" ^2 q. z9 l5 y% n7 G
3 H0 {+ ^6 W* H+ M4 D# B$ F$ q# o4 }0 |
REM For windows
! X$ v3 p1 e# s. {
3 D3 W/ F: X/ | U, {7 Ugmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
4 J v# m( G4 b3 bgmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
5 y3 d" q4 Q/ \# ?6 L2 h z: f. P1 H' X$ O
set ps=fig1.ps3 F$ T9 A s! M5 n3 s- z+ r
gmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%6 {! g9 `6 O% d! J# _
gawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%% \ k' W" o% i1 K9 m
gawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%
; o7 U5 |5 @. v/ t% A. d" n4 U
5 ?' ] Y+ S8 Z5 D7 u% Bgmt psconvert %ps% -A -P -Tg, V+ v1 _$ k* p/ q4 q4 _
注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。
' \ E; z* V$ c! f5 o
样本数据下载
% _' D2 ?: O: ]4 b0 P4 Z# Ahttps://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
8 j4 z' A8 S* `7 r( ~" v" Q0 n3 I
: B+ V9 H3 ]) A, V1 F! ?1 K |
