GMT 入门
GMT 实例
GMT 进阶
- 基础知识
- 表数据
- 网格文件
- 颜色表 CPT
- 标准选项
- 地图投影
- -JX: Linear, logarithmic, power, and time
- -JP: Polar/Cylindrical
- -JA: Lambert azimuthal equal area
- -JB: Albers conic equal area
- -JC: Cassini cylindrical
- -JCyl_stere: Cylindrical stereographic
- -JD: Equidistant conic
- -JE: Azimuthal equidistant
- -JF: Azimuthal gnomonic
- -JG: Azimuthal orthographic/General perspective
- -JH: Hammer equal area
- -JI: Sinusoidal equal area
- -JJ: Miller cylindrical
- -JK: Eckert equal area
- -JL: Lambert conic conformal
- -JM: Mercator cylindrical
- -JN: Robinson
- -JO: Oblique Mercator
- -JPoly: (American) polyconic
- -JQ: Equidistant cylindrical
- -JR: Winkel Tripel
- -JS: General stereographic
- -JT: Transverse Mercator
- -JU: Universal Transverse Mercator (UTM)
- -JV: Van der Grinten
- -JW: Mollweide
- -JY: Cylindrical equal area
- 配置参数
- 地学数据集
- 中国地理空间数据集
- 中文支持
- API
- 模块手册
- basemap
- begin
- clear
- clip
- coast
- colorbar
- contour
- coupe
- end
- figure
- filter1d
- fitcircle
- gmt
- gmt-config
- gmtbinstats
- gmtconnect
- gmtdefaults
- gmtget
- gmtinfo
- gmtlogo
- gmtmath
- gmtselect
- gmtset
- gmtsimplify
- gmtspatial
- gmtsplit
- gmtvector
- gmtwhich
- grd2cpt
- grd2xyz
- grdblend
- grdclip
- grdcontour
- grdconvert
- grdcut
- grdedit
- grdfill
- grdimage
- grdinfo
- grdlandmask
- grdmask
- grdmath
- grdpaste
- grdproject
- grdsample
- grdtrack
- grdtrend
- grdvector
- grdview
- grdvolume
- gshhg
- histogram
- image
- img2google
- img2grd
- kml2gmt
- makecpt
- mapproject
- mask
- meca
- nearneighbor
- plot
- polar
- project
- psconvert
- rose
- sac
- sample1d
- solar
- spectrum1d
- sph2grd
- sphdistance
- sphinterpolate
- sphtriangulate
- surface
- ternary
- text
- triangiulate
- velo
- wiggle
- x2sys_binlist
- x2sys_cross
- x2sys_datalist
- x2sys_get
- x2sys_init
- x2sys_list
- x2sys_merge
- x2sys_put
- x2sys_report
- x2sys_solve
- xyz2grd
- mgd77manage
- mgd77convert
- mgd77header
- mgd77info
- mgd77list
- mgd77magref
- mgd77path
- mgd77sniffer
- mgd77track
- legend
- docs
- subplot
- inset
附录
其他
gshhg
- 官方文档
- 简介
从 GSHHG 或 WDBII 数据文件中提取数据
gshhg 读取海岸线文件(GSHHG)或者政治边界和河流(WDBII),并输出结果。该模块会自动处理不同数据集之间的差异。可以设置是否显示段头部信息,头部信息的格式为
ID npoints hierarchical-level source area f_area west east south north container ancestor
hierarchical-level 海岸线多边形的级别,从 1(海岸线)到 4(湖中岛等)。
source 来源为 W(World Vector Shoreline)或者 C(CIA World Data Bank II),如果湖泊是河流中的湖泊,将使用小写字母。
west east south north 是矩形边界。
area 是多边形区域的面积,单位为 km^2
f_area 为其父多边形的面积。
container 为包含此多边形 ID 的 ID,如果没有,设置为 -1
ancestor 是 full 分辨率数据集中多边形的 ID,通过该分辨率数据产生当前多边形。如果当前多边形也为 full 分辨率数据,则没有祖先,因此设置为 -1
对于线数据,其头部信息格式为
ID npoints hierarchical-level source west east south north
更详细的信息见
gmt gshhg binaryfile.b [
- binaryfile.b
GSHHG 或 WDBII 二进制数据文件,随 GSHHG 数据补充文件一起分发。可以使用 5 种标准分辨率(full, high, intermediate, low, crude)。
可选选项
- -Amin
输出面积大于等于 min 的特征 [默认输出所有多边形]
- -G
在数据段前加 % ,后加 NaN 以便输出数据可以被 Matlab 或和 Octave 导入
- -Iid
只输出 ID 为 id 的特征信息。使用 -Ic 仅输出所有的大陆,见下面的
将整个 intermediate 分辨率的 GSHHG 二进制文件转换为 Matlab 可使用的格式
gmt gshhg gshhs_i.b -G > gshhs_i.txt
在 full 分辨率河流数据集中提取段头部信息
gmt gshhg wdb_rivers_f.b -L > riverlisting.txt
从 high 分辨率数据中提取湖泊,但不包含河流形成的湖泊
gmt gshhg gshhs_h.b -Qe -N2 > all_lakes.txt
特定多边形
全球范围内的较大的多边形:
ID
Landmass
0
Eurasia
1
Africa
2
North America
3
South America
4
Antarctica (AC grounding line)
5
Antarctica (AC ice line)
6
Australia
7
Greenland
8
New Guinea
9
Borneo
10
Madagascar
11
Baffin Island
12
Indonesia
数据文件
gshhg 读取的文件为原始的二进制多边形文件,可通过 NCEI Shoreline Coastline Resources 或 SOEST GSHHG Site 获得。需要注意的是,这里的数据和 gmt 中 gmtselect,grdlandmask,以及 coast 模块用的 NetCDF 格式的数据是不同的
技术细节
希望使用自己程序访问 GSHHG 和 WDBII 数据的用户可以查看 gmt 源码中的
gshhg.c
和gshhg.h
并熟悉他们的数据格式。GMT 并没有用来读取这些文件的 MATLAB 代码,但 Mathworks 和 IDL 都存在此类工具。但是,这些工具一般都不频繁更新,因此,随着数据文件的更新,这些工具可能已经失效。因此,在这里给出一些技术细节。GSHHG
GSHHG (The Global Self-consistent, Hierarchical, High-resolution Geography Database) 最初为 GSHHS,是一个海岸线数据库。其中包含大陆和岛屿(1 级)、湖泊(2 级)、湖中岛(3 级)以及湖中岛中湖(4 级)共 4 个等级的封闭多边形,使用
Douglas, D. H., and T. K. Peucker, 1973, Algorithms for the reduction of the number of points required to represent a digitized line of its caricature, Can. Cartogr., 10, 112–122.
Gorny, A. J., 1977, World Data Bank II General User GuideRep. PB 271869, 10pp, Central Intelligence Agency, Washington, DC.
Soluri, E. A., and V. A. Woodson, 1990, World Vector Shoreline, Int. Hydrograph. Rev., LXVII(1), 27–35.
Wessel, P., and W. H. F. Smith, 1996, A global, self-consistent, hierarchical, high-resolution shoreline database, J. Geophys. Res., 101(B4), 8741–8743.
相关模块
如果你对这篇内容有疑问,欢迎到本站社区发帖提问 参与讨论,获取更多帮助,或者扫码二维码加入 Web 技术交流群。
绑定邮箱获取回复消息
由于您还没有绑定你的真实邮箱,如果其他用户或者作者回复了您的评论,将不能在第一时间通知您!
发布评论