Question

官方文档

triangulate

简介

笛卡尔数据的 Delaunay 三角剖分或 Voronoi 图和网格化

triangulate 读取表数据并进行 Delaunay 三角化，该三角化方法能最大程度地避免狭窄的三角形。如果使用了

gmt triangulate [ table ] [

table

一个或多个ASCII或二进制表数据。若不提供表数据，则会从标准输入中读取。

可选选项

-Cslpfile

读取坡度网格 slpfile （单位为度），并使用 CURVE 算法计算水深测量中传播的不确定性 (Zambo et al, 2016)。需要

为了确定点在多边形内，外或在边界上，GMT 会平衡数据类型以及多边形形状等因素来确定算法。对于笛卡尔坐标，GMT 使用 non-zero winding 算法，该算法非常快。对于地理坐标，如果多边形不包括两极点且多边形的经度范围不超过 360 度，同样使用该算法。否则，GMT 会采用 full spherical ray-shooting 方法。

示例

对 samples.xyz 进行三角化，在二进制文件中保存三角网结果，并创建网格

gmt triangulate samples.xyz -bo -R0/30/0/30 -I2 -Gsurf.nc > samples.ijk

使用 15 cm 的墨卡托图绘制基于上述同样数据生成的三角网

gmt triangulate samples.xyz -M -R-100/-90/30/34 -JM15c |\
    gmt plot -R-100/-90/30/34 -JM15c -W0.5p -B1 -pdf network

同样使用上述数据，但绘制 Voronoi 图

gmt triangulate samples.xyz -M -Q -R-100/-90/30/34 -JM15c |\
    gmt plot -R-100/-90/30/34 -JM15c -W0.5p -B1 -pdf cells

将 Voronoi 的线段组合成多边形并根据他们的 ID 填充

gmt triangulate samples.xyz -M -Qn -R-100/-90/30/34 -JM15c | \
    gmt plot -R-100/-90/30/34 -JM15c -W0.5p+cf -L -B1 -Ccolors.cpt -L -pdf polygons

使用自然最邻近算法进行网格化

gmt triangulate samples.xyz -Gnnn.nc -Qn -R-100/-90/30/34 -I0.5

参考文献

Shewchuk, J. R., 1996, Triangle: Engineering a 2D Quality Mesh Generator and Delaunay Triangulator, First Workshop on Applied Computational Geometry (Philadelphia, PA), 124-133, ACM, May 1996.

Watson, D. F., 1982, Acord: Automatic contouring of raw data, Comp. & Geosci., 8, 97-101.

Zambo, S., Elmore, P. A., Bourgeois, B. S., and Perkins, A. L., 2016, Uncertainty estimation for sparse data gridding algorithms, Proceedings of the U.S. Hydro Conference,National Harbor, MD, 16-19 March 2015.

Zhou, Q., and Liu, X., 2004, Error analysis on grid-based slope and aspect algorithms, Photogrammetric Eng. & Remote Sensing, 70 (8), 957-962.

Shewchuk’s Homepage

相关模块

greenspline, nearneighbor, contour, sphdistance, sphinterpolate, sphtriangulate, surface