本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。
! c. {$ r$ A2 u7 o3 QWoS数据
3 h* I. Y B5 H% c; n+ H根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,505
# U4 |2 T" N/ D3 ENATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747
* s/ e+ D2 o: l. p+ W; RCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735
! n& o9 I. X& w V+ \/ oINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.313 `' q* V8 [7 l; D! k7 y1 d( t7 C
CHINESE ACADEMY OF SCIENCES 629 7.396
, T& P: v9 F# i- V7 \* i, l$ Z Y" _UNIVERSITE DE TOULOUSE 517 6.0799 e5 Q$ m3 P( a
UNIVERSITE TOULOUSE III PAUL SABATIER 498 5.855
5 h8 z$ i1 P7 j2 ^# TNASA GODDARD SPACE FLIGHT CENTER 483 5.679
9 x; j( L. i9 p+ q6 q8 k' YLABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162
! u. [5 w! C$ ?/ hNASA JET PROPULSION LABORATORY 431 5.068
) j1 M! J" v8 R( o...国家数据： USA 3118" }& s: [$ F/ I0 h( c7 o7 |+ ]2 e0 P
CHINA 1493! N, r0 ]8 Z; k+ I7 N3 d: [
FRANCE 1385
, O! K' W7 Z0 y: W4 U8 e8 k" Z! C# yENGLAND 804
7 H9 o& n1 ~) R7 v* QGERMANY 748
; _8 U0 ^0 z5 {1 RAUSTRALIA 4219 n' L4 Z5 n4 S0 r* @
ITALY 421. X* @0 B- |# T2 J" Q! q X
SPAIN 371
% n( D1 O! b5 o+ s% g( BNETHERLANDS 340
' o( J% o' H$ P* F' z/ u7 {% m4 P) `CANADA 323; O' L6 o7 n- I! S
3 \ |6 a, Q3 r* p: \, s4 o从地名到坐标 * }( L: C0 J% [8 Q4 t
使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati1 K7 W) D; z1 z7 Y* W
NATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.883056
# q% l! v" Q% g# b0 y# X( @4 MCENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037
, P |8 ~$ F c9 {0 MINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.94 L' l; @% R' e g9 _6 P% x
CHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327
8 `7 {1 S$ B! h+ K- ?1 e1 LUNIVERSITE DE TOULOUSE 1.4 43.6
1 q9 @% ?( p* T5 u* OUNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586" l6 Q g- q$ V8 ]2 \3 x: I
NASA GODDARD SPACE FLIGHT CENTER -76.852344 38.994955 \2 }3 r h) n9 w5 g: s
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65( I, d" G' s1 b6 ~) z. G3 F$ x/ |
NASA JET PROPULSION LABORATORY -118.1713944 34.2013081) s* d% V t) A+ g: S2 R3 e
! \- ?. r0 H# u0 E* i* R, H; s
2 F3 [$ Z1 N1 P1 d; J绘图
' |/ C: {4 ~. W2 u数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows
4 E7 ^0 ^+ R$ Z, a1 A' u& b. Z) k( B7 _! s0 X
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
5 O0 [9 D% {. I, O% Ogmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black ' S, f% Y9 u, M- ?' g
- }- y3 d7 w3 B' F8 o# q
set ps=3d.ps
' X8 h0 v) S5 H B2 P8 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%6 b1 Q: r% R0 D+ ~# |
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 y/ ^$ D/ Z, Z9 ~% I
$ @! Z: {7 b2 O' w s3 W+ Mgmt 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%) T) J! ]" q- [7 w+ L4 x9 k
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%7 r! i0 s9 `/ \* q7 L2 b4 G% `
gmt psconvert %ps% -A -P -Tf
6 a- D9 v! N0 m* J
# V& H5 r; y# v Y9 `/ O- t' j结果
1 I- f/ I f9 M" b- d6 `4 _1 {; o! E图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。
& }. t% I, K, e. b. p ~
附赠平面图绘制代码 % [6 x0 c4 r! a0 c2 m, ^! G2 r
0 n c$ Q0 y3 e4 T5 {* u
$ e. x* ~) e: [
% D. E* }& Q& v* G, \4 ?" gREM For windows+ X$ `) K W# Y
2 w$ ~) `) ^9 D6 }& h# y
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p% w, }5 d8 |' O3 I- ?' i
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
$ }/ p5 d' L5 l, V: e9 R( X$ [, x5 v t! j0 j
set ps=fig1.ps1 k' @# w8 b1 p( j7 O5 Q' p
gmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%
. f; w/ Z* b$ wgawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%
Q2 C& `0 `8 b5 H; t& Cgawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%
9 d# Q8 G* ~* h0 G( C' i& }6 H
. q( J9 E& @( u/ u3 hgmt psconvert %ps% -A -P -Tg, K! F' o8 o* @9 L t7 U7 l
注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。
7 U7 x! `5 }: }# x) q' S
样本数据下载 . [/ j# d+ u, g6 _6 Q# R+ b8 U
https://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
2 a" B: J8 b, N9 L( C
! }' b+ U( n- S/ F! R5 O2 _8 M9 A" I
|
