本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。
( K/ n# p _2 C& pWoS数据
( J9 h. n4 U& [4 m0 P* }根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,505, u" b1 y+ z+ e4 n% n: r! i* V
NATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747
2 A% x2 H7 E/ G, ECENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.7355 v- C# D( Q3 ~3 a+ e# U4 J
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.313! O* h, [2 x% H" o$ ?
CHINESE ACADEMY OF SCIENCES 629 7.396
% s5 R. e+ f, L1 A( k, uUNIVERSITE DE TOULOUSE 517 6.079' K0 E' ^2 s. G7 ^9 a
UNIVERSITE TOULOUSE III PAUL SABATIER 498 5.8557 J7 |* p% k2 ?
NASA GODDARD SPACE FLIGHT CENTER 483 5.679
3 f9 G5 e& @) uLABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162* Y* y0 U( }# B
NASA JET PROPULSION LABORATORY 431 5.068
' }, U9 ~7 X9 ]" j" @...国家数据： USA 3118; {$ j& e6 f5 X) Q# n
CHINA 1493
; M- W- V4 t0 V+ z6 K% C8 y% T; gFRANCE 1385/ Z8 d) U, m0 a- i0 T3 k0 o& {
ENGLAND 804. Q5 n9 m0 R6 ^+ g0 v: X% h
GERMANY 748
& ]; @5 Q) Q6 v4 S- hAUSTRALIA 421
. ^& U$ O* J$ T2 sITALY 421
: {4 {2 Z$ L* P" u2 HSPAIN 371+ a, ^% w. }$ _2 L& e0 W5 {3 p" L
NETHERLANDS 340
- q4 ]; ?- O0 }0 s, H5 tCANADA 3231 x, k) U3 c# G! ^* u
4 c7 R$ N; C1 B9 V
从地名到坐标
" q h" B" }" d使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati
' g8 W6 u8 x& P7 y+ t) |+ X) ~NATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.8830566 n, r' s4 `$ h/ p; i6 U D' C# v
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037! P/ P0 Z' ^8 E3 v. r8 T. k
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.9 W+ m% g# e0 f, }+ o% Q7 `
CHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327
+ h! h E* W# S3 B+ A: v& C O/ o, LUNIVERSITE DE TOULOUSE 1.4 43.6
' L+ l$ o7 h+ O# bUNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586
. C4 i9 F$ [/ X: c% X6 gNASA GODDARD SPACE FLIGHT CENTER -76.852344 38.994953 v4 `6 @* g" O2 X' Q5 z
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65
, a! l0 U, r) n Z4 ^7 l% oNASA JET PROPULSION LABORATORY -118.1713944 34.2013081
6 i! t' }9 ~0 L3 `/ `& V7 J/ W- q, C/ F& c
9 o5 J4 M1 t! x0 Q! h. G" R
绘图 T! [" D9 _/ K) M: x0 j
数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows. o: `' J; ^" x3 s/ ]) E
3 W+ X6 W& I5 e+ b6 S7 O8 E% Y$ u/ }gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p. W0 U- r! X: ~0 V( p
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
1 P, y( ~* i" x) L0 }8 }3 j7 f# n
1 F% i" K% @" p( b1 l7 }/ Q0 Aset ps=3d.ps
5 L% L/ L; V) S9 @. v! {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%
; [0 @* o5 S! C9 j+ I3 ?9 Vgawk "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%) w" J& ?6 F5 ]! Q$ G
/ S/ m6 q) C, u+ c6 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%
" m: L* o- ^, Z9 Agawk "{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%3 Y8 j. \8 p; w2 @* i& v5 T: V+ ?' J
gmt psconvert %ps% -A -P -Tf+ K( L6 y) k0 P6 h2 W' ^
5 i3 K8 B0 e5 d" V结果 ( [/ m, s% k+ x) A$ u( {/ t
图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。 " ^$ i9 v: }# Y1 M- A& C
附赠平面图绘制代码
; V) f+ p- J& _
4 ]9 }( g' R u5 U% d& t$ L9 _0 m2 ^; H% y y
4 V# Y* A, ~2 QREM For windows6 V% s" p8 c; i* {
6 M# C1 E! X& c: wgmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p! o1 Y, E) }2 Q6 m
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
D8 c1 @1 f) S# z, O" F6 o4 q9 M$ V; R
set ps=fig1.ps; ]8 I2 k$ \! |" n
gmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%
& S3 G; ~6 b$ H' ogawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%
4 ~" f$ s6 J3 {! t' s* mgawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%8 G3 e" k ]$ ?5 L$ y7 s
& v) y& c. A3 z E8 D4 z* @8 U
gmt psconvert %ps% -A -P -Tg
* f; C3 Y/ b7 w4 p7 m注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。
' Y3 E+ H$ W! V4 I
样本数据下载
/ S8 I$ F: |& D7 l7 ~( ?https://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
+ t" u$ G( ?2 x% t' z4 E( a3 }7 p) b, b( E- ]. a" n" ]+ J
|
