本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。
. H. N7 Q" w/ D7 iWoS数据
# x D0 h2 V( x- |( F9 C/ y' J2 n根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,5050 o$ i! O& e3 b0 S1 G( U
NATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.747( I' \4 h( T' F: ]! a, a
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735
& b+ B4 c0 l( @INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.313
3 D& b. u4 \8 n+ v4 \" t+ _" h2 jCHINESE ACADEMY OF SCIENCES 629 7.396+ e* Z( r) c+ Q* i
UNIVERSITE DE TOULOUSE 517 6.079) R* U/ p! @& v
UNIVERSITE TOULOUSE III PAUL SABATIER 498 5.8558 `$ V, ^& Q1 C& o6 U b* ^
NASA GODDARD SPACE FLIGHT CENTER 483 5.679 t9 b7 d/ ?5 u8 @
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162) ?6 L6 T: w6 t3 o- [) e2 q- G: s
NASA JET PROPULSION LABORATORY 431 5.068
( o8 U* u6 o, {* _, B...国家数据： USA 31186 N% W; |/ A5 c0 B' r
CHINA 1493
2 ?5 {3 t& _5 ^; Z7 kFRANCE 13850 E/ W5 n, h- C/ z# Y3 G7 R
ENGLAND 804
: K8 ^ K0 r/ z6 Y1 rGERMANY 748
1 [; [; z4 J7 s8 u- b! P9 l7 PAUSTRALIA 4210 a# j$ c8 C' o7 }6 X; ]
ITALY 4218 d s w1 ]7 m* r4 P; Z
SPAIN 3714 B' \5 w! p% E7 c/ J8 L5 u# {
NETHERLANDS 3402 M2 A( D& @9 }5 y, @
CANADA 323
8 C' f: o9 w! |3 D6 s$ e9 C1 U' p8 @- x7 P9 z
从地名到坐标
. V0 ^8 E* X# ?; x/ Y& t: u# n使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati- z/ x: D1 K% M/ }" ~
NATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.883056+ ~: G+ T& u7 t) N! y; W
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037
. B P$ x0 A! i( r& }! c1 h* R; W0 TINSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.9) k+ K/ |( V4 W& b
CHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327. ^8 K @0 P/ J7 e
UNIVERSITE DE TOULOUSE 1.4 43.6
: ?( x x* o$ ~" i$ gUNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.56195860 I4 e1 A ~/ W& C
NASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495 p$ M/ j. T* H/ ]# r+ h
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65% U, j( S- T0 O( i* s- ]+ X
NASA JET PROPULSION LABORATORY -118.1713944 34.2013081
4 M/ T! e3 N! `& X3 r
1 n$ y8 X& w) f5 J2 B
1 K0 Y/ n4 k' ^ K绘图
2 W6 d1 k0 b Y7 g6 D数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows' j( g6 w9 }2 Q9 `2 H
/ d5 ?) S* [+ h- _4 C0 A% v; U! [0 ^gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
; Z. F2 g x& V, V" {! rgmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black " b, G* K* g1 b. A$ t
+ w" W4 j/ s. t/ h9 `& _2 dset ps=3d.ps g2 ]/ N' D1 O2 Q, k3 C- z
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%# |' g5 h" U" E+ S" @' i1 I
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%
- F. `5 o9 h9 N5 u8 ?
; S+ M! ?( M: P* x) kgmt 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+ E! K- b9 w# Xgawk "{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%
- C6 m3 I/ ^6 S" m3 pgmt psconvert %ps% -A -P -Tf
0 w9 Y3 a6 r3 [$ N: J
# k& Y; e9 x' t8 g- B3 @结果 O Z/ m# d# \3 \0 T3 g! v% R5 u
图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。 " J1 u+ S1 z z% L. ]
附赠平面图绘制代码 5 K: s, u+ O3 X1 e
; b- N( r$ p% Q' I9 {, L, _& O" ]4 V$ F$ X7 \7 B. T
& k, P, F* |! @2 v
REM For windows
& \; O/ z# R! Z$ z) f3 b) g) o
4 l3 _+ G3 q3 r5 J* t3 igmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p7 a9 K: s( ?% y& ]) v( z
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
4 {" c3 I/ k; L# }3 ?% ]+ l Y1 _
1 c/ ~! m( ~( v* zset ps=fig1.ps
( F* ^. q' q% U: k2 wgmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%
9 L* c8 V" L) o6 d- T2 J/ fgawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps%2 w0 W9 E; J# P9 {
gawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%$ Y0 B1 u, V. _, K5 g- z( ~: g& T
: |8 C, h b; K5 ~, j! |
gmt psconvert %ps% -A -P -Tg
1 {4 E7 w8 l; B9 E4 w注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。 . b4 Y$ y7 O5 K& ~
样本数据下载 " T) t3 s5 Z/ v' E7 P
https://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
: g. }4 o7 ~, }1 z! J L( M- W
% H* C1 \# s! ]$ \3 e2 P1 R |
