本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。
8 p( z; z9 T! ^WoS数据
/ C' x" i* U0 z' t2 b: N! `根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,5057 k* K( C4 j8 N Q) e& E3 `6 t
NATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.7473 V, l1 j7 O. R" K' K2 i
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735
* }7 y0 m2 j8 y$ c9 }2 g( x0 qINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.3130 _* Q8 p- w6 Q
CHINESE ACADEMY OF SCIENCES 629 7.396
) Y, b- Y: m" v @8 p: T# R6 {UNIVERSITE DE TOULOUSE 517 6.079" s/ m0 @- J1 k
UNIVERSITE TOULOUSE III PAUL SABATIER 498 5.855) Z S! X, G/ F1 j4 Y# s& e
NASA GODDARD SPACE FLIGHT CENTER 483 5.679* `. M3 s- t* }
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162
$ g8 D c# l- ~: b# Z$ n4 ENASA JET PROPULSION LABORATORY 431 5.068' } d# x+ ~" b) l9 d5 u* v
...国家数据： USA 3118" i; \; I6 b1 E) U6 P) d3 \
CHINA 1493
% J( ]4 x8 Q4 V) e- R; T" C5 UFRANCE 1385
. _# \8 S7 m* u, g6 SENGLAND 804' Z# q; W* w# y/ l: x5 p
GERMANY 748
3 Y9 h$ w, Y' j! U2 w) X3 I) ]) H+ GAUSTRALIA 4218 ^$ {, p6 b( O
ITALY 4211 r l) H$ h4 U; r
SPAIN 371
( Z$ X, }' F. M) u8 aNETHERLANDS 340
0 G" }3 y$ a* u9 w' |" kCANADA 323" y+ T% M' c0 g6 B6 Q
' ~! U7 t m% s* F
从地名到坐标
) {' _, t) \3 x6 L) F使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati
6 ?) A1 z+ h4 C+ VNATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.883056
2 R0 }0 ? W2 d$ n9 _' i0 _+ v. {CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037
1 a& K% F' [; [* aINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.9) R4 ?" r+ Y; J; ]
CHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327 E3 z+ h% |0 `/ T, [, k) o
UNIVERSITE DE TOULOUSE 1.4 43.6
( I/ G- R, _* `: eUNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586
# j- N) K2 N l+ M: K) G, d5 H% Q+ dNASA GODDARD SPACE FLIGHT CENTER -76.852344 38.994951 k) y# G3 m) Y. N9 e
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65
, v+ b; ~! h4 X/ C( g$ aNASA JET PROPULSION LABORATORY -118.1713944 34.2013081+ p8 x" ~0 R, N9 @* e ~# C' m* k
3 S F/ \0 g3 w/ f
9 O3 z6 h7 }8 z9 j" z" h
绘图 - x# y. h2 r* i. v& D' S
数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows
- B+ F# x+ h) z' Q
% v4 E0 g; z- a3 n7 }- g- ?gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p* C! K: l( e' F1 }- @" w" }8 U
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black ) G8 d' q8 M9 c* N' [' {1 h8 K
* i2 Y0 B4 ~* {# E; U, Z0 L+ @9 d2 @set ps=3d.ps# D: e7 @# P: w) S2 l: |4 ^- n7 s
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%
8 ~9 W4 o' c e; X4 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%2 y0 n/ r' [1 m) f, ?
$ o6 U* r" h5 ^* [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%
/ ? B- L, ]! T! [/ a* s9 j- n3 [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% x6 V4 J/ Y0 e9 \5 \# {
gmt psconvert %ps% -A -P -Tf9 Y; S0 S1 [3 i9 Y7 [1 v
/ |: Y+ a% _! O4 n: `' B" A结果
T4 m( l C! Y4 O8 m图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。
+ v: K7 z6 ~( q- Z% I7 R- `$ d9 w
附赠平面图绘制代码 8 U, J- L9 L7 P J# P# l6 @, l4 O
2 j B) U# J" k# m7 w$ Z0 l& q
+ U( }1 c( N9 w$ h; D4 V# o7 J: Q( j8 ~
REM For windows- F) h4 L- ?3 o- U3 ~+ l
6 P: L$ I; |1 e: e' N2 b
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
* b% W" Y( P8 U7 Hgmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black 2 f. c$ ?6 u0 S
; O0 w3 D, b/ ~& d4 [
set ps=fig1.ps
( ^8 c( h( Z5 t4 K7 C! k @% ggmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%( ?- O/ r1 N- _
gawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%
4 z$ f( S1 ]4 E; _, Ygawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%( l6 I/ N0 |, p6 u' H( |
3 s- G; J* m5 c2 P
gmt psconvert %ps% -A -P -Tg. G4 A- l0 Y* r' N, P9 T- d" x
注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。
; W; m3 ?* G; p
样本数据下载 $ X- U p5 s& o1 w$ L5 ^% Y) n
https://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
6 B x9 }+ _" H8 e- z: _9 y3 |7 ]0 P& d) V; `
|
