本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。
$ K t" s7 b& e! O' w5 oWoS数据 . p: `8 Z2 K; s
根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,505
4 v3 N2 G; r+ T1 c G; ]NATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747
# V7 ?% b5 y3 ^) I: L$ VCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735
* s- e* s' |, rINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.313
* ~) v. s* s: Q, v1 `* jCHINESE ACADEMY OF SCIENCES 629 7.3968 y$ {. F6 B5 Y5 _, c
UNIVERSITE DE TOULOUSE 517 6.079 j1 u/ n' `0 \% b
UNIVERSITE TOULOUSE III PAUL SABATIER 498 5.8551 l6 a* a/ ^) ?! J& E. g
NASA GODDARD SPACE FLIGHT CENTER 483 5.679% w1 ]4 h9 y B8 g1 I* U
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162! U" ]. F$ _4 t! m9 v
NASA JET PROPULSION LABORATORY 431 5.068
/ r4 z) \/ ~! @; S, ?2 |...国家数据： USA 3118% [; V- d* o) [, d2 V1 F
CHINA 1493
8 s1 U# \) e' y; L+ s1 W# ~8 b2 ^( |FRANCE 13852 P# s( }( `9 w- `0 m* d7 N( o
ENGLAND 804
& c; R5 v0 m8 ?9 q( `2 @( ]GERMANY 7484 T% v$ ?5 B8 @! u
AUSTRALIA 421
}- h6 R# Z4 u# J) uITALY 421
$ \% r: w ^, [5 N: V6 uSPAIN 371
% g9 R& m& \0 g, S: eNETHERLANDS 340' n0 G) D* c% w' v: ^ `( ^5 ?
CANADA 3237 ?6 H! m/ K8 M! T: n- n& E
2 s& L8 \4 U) B d/ u }
从地名到坐标 ' T9 p7 d* v! C1 |' J+ p
使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati
Y$ e0 i+ n' k; X5 u3 ]3 b: D6 `NATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.883056
7 Y% O+ I7 U& a3 H& ?# F5 CCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037
, [# C- T) b% a* p3 @- zINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.9
: H! f$ b8 T' h/ s+ oCHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327
/ A: u b) j: o4 @6 M8 @3 lUNIVERSITE DE TOULOUSE 1.4 43.6
/ Q8 C+ n: i/ r8 i* R) F& AUNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586
( A. `' B/ `$ c* DNASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495
$ x* d- W h" D- H3 K& ^LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65; G1 y( t2 u" E* z$ }5 F' X
NASA JET PROPULSION LABORATORY -118.1713944 34.2013081- j' v2 e" o d) u- j# g
1 z g2 V4 T$ a# s: ^
# z! y, E+ n0 P绘图 . s+ M- J, D& J: I; h3 `0 k
数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows Y% d& @1 a8 r$ {4 W) k+ `
6 g& x. B2 P; Ngmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p7 P' [4 e, r4 N. `+ p
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
) F* F2 ?% [" ?' T: ]4 G
5 D: G3 v; y0 \1 H }& A6 ?set ps=3d.ps
; @/ l2 [, y0 s$ Q( Sgmt 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%
5 I9 J' d' ~! L: G; O* x& Bgawk "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%: ?( G5 J& a+ d' N5 G+ F
; P- v* A9 L, m: v* S( {8 n' M1 W
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%
- a# R# n7 }. ]: Tgawk "{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%/ X* G# u. l9 E3 D! e; a1 U8 h5 z
gmt psconvert %ps% -A -P -Tf% a6 g- `* n; ~0 X. ~
" q& F! P# e; S3 j& y/ k/ p5 U结果
3 R4 C, E, ?% k ~图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。 6 f1 o& |* Z1 [! ]0 z
附赠平面图绘制代码
! @. T" n: s( z1 I- E2 U# U1 _2 ~! J6 X0 y4 N/ |5 }
/ T7 q/ \1 a P+ F/ Y5 G
! S1 J: J- K( t: Y
REM For windows* `* o+ r8 K& U& F) k" |8 G/ [
4 G# s; z0 M+ t2 C3 F% v2 Cgmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p) O! X- f+ @# Z+ k) _) s* g
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black ! x! V% [+ d, F1 R) u8 J! e7 j
6 b4 p% P4 D" m/ x* ?
set ps=fig1.ps& ]; ]0 b) S4 o
gmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%
# `( M0 s# r- H5 E4 bgawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%
6 O3 x1 y: Y* Y3 Dgawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%
- @% B [' l0 d. @8 N
" V! E4 ]! ?. R% v1 o1 Wgmt psconvert %ps% -A -P -Tg
" E: V# J* E) q注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。
9 `, x2 l# R( T5 }, D1 h
样本数据下载
" n& b! a( B6 v/ }https://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
! d) i$ _* I6 d& B. N! W3 j
4 H4 M/ O4 @. n9 g
|
