本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。
( w I1 ^- h8 q( T8 ~7 s2 QWoS数据 6 r- u5 f8 a4 a" k1 H
根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,5057 H6 D: N M: s. h3 @) _
NATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747! D$ T b( Z/ x: B" A* b4 X t
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735
6 o1 V3 G$ d; N. w* G7 n2 {INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.313. ~; \7 K2 t6 _. S5 F$ b$ z
CHINESE ACADEMY OF SCIENCES 629 7.396: {8 h$ c6 v: T I' `2 ^: }0 Z* ~
UNIVERSITE DE TOULOUSE 517 6.0793 F7 t; N8 ^2 I7 m& E7 ?- Z
UNIVERSITE TOULOUSE III PAUL SABATIER 498 5.855* ` o3 b! h- K+ H& |
NASA GODDARD SPACE FLIGHT CENTER 483 5.6795 g4 g6 Z9 G+ B c0 U# o
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162
" s9 z B. ]' i% FNASA JET PROPULSION LABORATORY 431 5.068
! N H2 z: t5 b/ @...国家数据： USA 3118
: d2 F. g' m! ZCHINA 1493
3 F$ Q; b7 |% C# ]) ]" cFRANCE 1385
8 c5 U: {- r( x: b& k) pENGLAND 804
" p4 ~0 b7 t s- K) j; g* U+ kGERMANY 748, J+ I6 o* ~5 ^
AUSTRALIA 421
- _2 B$ M x/ ?6 p. {: A3 R) Y9 |, hITALY 421: _* t# E1 S1 o; F5 I: H w
SPAIN 371+ J! m3 E0 g1 T! f8 I: y
NETHERLANDS 340
. [: l u6 `, W2 u' v( ^ s" j BCANADA 323
. W, A( H3 [! d1 H4 T- ` A0 \7 b- [1 s& s+ Y7 u9 M/ m0 t( [, M
从地名到坐标
, v; ~/ @1 r$ B! w2 C8 ~+ k K使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati
4 i; q- I! [+ E; ?' I, V% K& uNATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.883056& }6 O* D7 d, ]( }# f: S
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037; a1 Y0 ]& @1 |) t C
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.9 ]5 x& I* X) c1 m: B: }' h6 i
CHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327 {! ^6 j( T3 g. z
UNIVERSITE DE TOULOUSE 1.4 43.66 y# P, c7 N+ Z8 J2 U
UNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.56195867 U- }8 {: H2 W) q
NASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495$ c6 z( E9 ]7 I1 V
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65
- V# E! Z; M: C6 C( Q7 KNASA JET PROPULSION LABORATORY -118.1713944 34.2013081
1 e1 @$ F- x1 N' A( @% t! J4 S! F
* w* @. }2 x! b# W! ]
绘图 " ^4 \4 r1 C( d/ F5 v. k
数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows
, l8 c# }# m0 g5 n
1 i: [* z8 w0 K) D( a3 pgmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
. @" O% I1 j" o, vgmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
% b" i) z; S$ K" N! A2 l0 Q q) N- I" {) s9 x. M
set ps=3d.ps: h& F( h, v) B3 J2 d& T! c+ C* r
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%
1 x9 w" ]6 w1 g ? h$ J3 xgawk "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%
" ?4 D# }5 I# }) W) U
8 o0 [# a' O7 l( u" ^& x/ [ }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%# G1 \7 q0 W$ O* l6 M
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%
! _: Z$ ^ y! u& K1 M, b7 vgmt psconvert %ps% -A -P -Tf
; W5 i2 r/ u+ T% g8 D' E' D, M2 T
结果 : m+ u1 c* j. X( {, `' }6 \3 [; `
图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。 / l3 q1 E- L' s7 \
附赠平面图绘制代码 9 ~' q1 }. i3 B+ p0 q" q+ K
; @# ]/ r! [4 y. k1 n' o
, H- ` ?7 O4 Z, s' w. m; k5 E! H% Z) H7 e
REM For windows
; u% {# s( i! N) x
1 |; ?( e; @4 pgmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p7 R) ?+ ?$ b, x2 D! A+ X# V
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
2 C% M# P2 L7 ]3 M7 p0 ]: v- S# y: Y/ f- m; R' F& P
set ps=fig1.ps( \3 b9 U8 j) Y
gmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%# y& j+ [. i* ]6 [
gawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%5 ^! Z# ?; Q$ _
gawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%. e+ q8 y( E4 ^0 j4 J2 A9 Q
5 q: Y; P" ?9 D# H; U
gmt psconvert %ps% -A -P -Tg0 v) y! @9 L+ O* Q) P( Z. k" K
注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。
9 N# p' V8 Z" T
样本数据下载
% q2 Q4 }4 O, r* @5 B& |/ ^https://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
% ]: H3 s5 _0 Y, W! g
' a' \- k1 Y2 q& @$ y |
