本期再次回到GMT绘图，介绍一种使用GMT绘制3D柱状图的方法，本例子数据采用的是卫星测高领域的研究国家和机构的文章数量。本图件可以直观的展示某一领域研究机构的空间地理分布和强度。对于综述文章或者研究生的文献调研有一定用途。
# R4 A5 T; K- Z3 B7 b; m# aWoS数据
2 P) ~+ G- J, o$ M根据前几期文章介绍的方法，从Web of Science数据库检索相关主题，例如本次我们采用satellite altimet*，得到8000多篇研究结果。下载研究国家/地区和研究机构的文章数量数据。 机构数据： Affiliations Record Count % of 8,5052 N1 t, C7 F% ~7 J5 h
NATIONAL AERONAUTICS SPACE ADMINISTRATION 914 10.7472 n1 I1 ` t4 R- ^2 R
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 828 9.735' S4 i4 y8 z' Z
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 707 8.313
6 q: Q% l6 B3 Q/ }0 K5 eCHINESE ACADEMY OF SCIENCES 629 7.396
8 f+ S" h0 q4 x5 G" R, i7 yUNIVERSITE DE TOULOUSE 517 6.0790 n0 H( E! X0 i$ |2 D$ M
UNIVERSITE TOULOUSE III PAUL SABATIER 498 5.8555 i* J) m+ Z* \7 A0 @# j! m
NASA GODDARD SPACE FLIGHT CENTER 483 5.679
5 d4 z7 j2 f1 b9 K [LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 439 5.162* a) t. X" n, U4 H5 U
NASA JET PROPULSION LABORATORY 431 5.068
* A4 W4 S! r8 @; w2 T...国家数据： USA 3118: l8 ^ D7 p; r# I" ^& f
CHINA 1493
2 e+ [$ E' Y9 `. ~0 l1 P: N! N* `FRANCE 1385
, g7 V, W8 k0 `2 p' UENGLAND 804& z1 ?4 ?1 \, m2 s, y
GERMANY 7489 i+ i5 A% p* V& E( x5 U" s: L
AUSTRALIA 421
9 C/ n- S, s/ u3 ]3 q/ bITALY 421
, K% h; O4 G5 S) M/ lSPAIN 371& q: B; V, ^3 j# [% x* B
NETHERLANDS 340
% `- r8 S( s5 iCANADA 323
- |: Y5 P4 \/ c" P# N$ r6 p6 X9 u8 r% M1 U9 u/ b
从地名到坐标 ) t3 v+ n2 W2 p. ]: x
使用google sheet结合Geocode插件对研究机构和国家名称进行坐标解析。 得到： Affiliations long lati
- u0 _, m+ N3 x) a' f3 B" L# Z! GNATIONAL AERONAUTICS SPACE ADMINISTRATION -77.016389 38.883056( X) Z6 T: k) e: p1 y; G
CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 2.2639934 48.8476037' [: t' d9 m j
INSTITUT DE RECHERCHE POUR LE DEVELOPPEMENT 2.3 48.9
" n* v0 }' ~2 i& }. lCHINESE ACADEMY OF SCIENCES 116.3356321 39.9106327
! K! g( {( h* {' w& T3 dUNIVERSITE DE TOULOUSE 1.4 43.61 L* ^# M5 s3 X
UNIVERSITE TOULOUSE III PAUL SABATIER 1.4699399 43.5619586 {$ G5 b8 G- u' g. `+ j+ r
NASA GODDARD SPACE FLIGHT CENTER -76.852344 38.99495 Z& z$ D: l/ P' s& V9 Q. r
LABORATOIRE D ETUDES EN GEOPHYSIQUE ET OCEANOGRAPHIE SPATIALES 1.45 43.65$ ~2 }$ b5 Z- @$ _0 j: \
NASA JET PROPULSION LABORATORY -118.1713944 34.2013081
* U$ f; h6 J2 |( b. F1 f5 d+ I' D/ e; P1 i8 Z! `
0 R6 |1 Q. V3 ]8 ]8 Q: O/ l+ z7 d绘图 , ]' W/ S, @3 R: r; ~; h2 B7 b
数据整理完毕后，使用GMT绘制3D柱状图。 REM For windows
! C& W; Z( _8 ?6 U/ `2 U! }' ~' }$ z8 Y! x4 A+ M9 ]" ~
gmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p
4 R5 U( g! w ]/ `( z9 ugmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
+ U- F; i9 j+ z& g9 F# x3 n& R8 l' u% m/ p
set ps=3d.ps
# K; W/ _# K0 j H5 T5 K% B8 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%
5 a3 J/ R. m# t& K7 u. Ggawk "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 ?+ Z/ w! i+ V& b
) k& G M/ o, B- y g' s. T! Agmt 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%
' n& y0 k0 J; h6 Z# B G, sgawk "{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%
9 @% S' ^4 Y. T9 ?% X9 j) N- Pgmt psconvert %ps% -A -P -Tf
6 k$ @8 B. k9 e! ?: Y# A) ^8 c$ R1 N8 V3 J6 p/ q4 Y# D
结果
/ f5 ~; e; q) q) c" d' j图中上部分为国家和地区的发表卫星测高相关的文章数量，下部分为机构发表的文章数量。首先，国家和地区较集中在美国、法国、中国、英国和德国，再以中国为例，相关研究机构主要分布在东部地区，最高值为位于北京（中科院大学）。
: I6 p j$ r4 A: n1 A
附赠平面图绘制代码 ; Q! [ n5 r# X0 I& H
/ c2 |; v! h* l8 H: Q% [% ^7 C" L" s
6 J) T: d9 g& |/ H
; g4 M( a' Q7 b: F, xREM For windows3 d0 z7 y9 ^5 Z/ s0 v# t% r4 x
# W; E$ {0 j' |* _0 h0 I- qgmt gmtset FORMAT_GEO_MAP = dddF MAP_FRAME_WIDTH=2p" U( l5 [: }2 e4 b2 p
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 7p,Helvetica,black
3 b/ \/ L0 {; M3 l
8 @' S9 t5 {- X3 Kset ps=fig1.ps
4 E1 f2 G/ G: Y' G. A. q, @gmt pscoast -JN4.5i -Rd -Dc -A10000/0/1 -Ba --FONT_TITLE=10p -K -W0.1p --MAP_ANNOT_OBLIQUE=45 > %ps%
# P7 D+ D! R* L* G ?gawk "{print $1, $2, log($3)/log(10)/10}" loc | gmt psxy -R -J -Sc -Gred -K -O -t55 >>%ps% W# ?% z$ w$ p% k* f/ k6 B! R' k
gawk "{print $2, $3, $4/2000}" loc_c | gmt psxy -R -J -Sc -Gblue -O -t75 >>%ps%
: y, f, y7 Y5 O5 X$ x
q/ S! Y8 ^. U/ y" igmt psconvert %ps% -A -P -Tg
! E4 _- O& y* w注意国家发文数量做了处理（log）。 图件缺少比例尺，有兴趣的同学可以修改代码加上。 至于哪一种表达方式效果更佳，小编也难以下结论，这或许和数据特点和个人喜好有关。 u% S. g( ^7 Z. {* G
样本数据下载 % d# A9 r$ ^; \1 i% y
https://www.jianguoyun.com/p/DSoojxEQ2PCQBxjioswEIAA
" t/ C- e6 F2 }( j7 |* G: D8 g6 z" |- y( K" o1 d& t$ Y
|
