Question

贡献者

周茂

---

官方文档

sphtriangulate

简介

创建球面 Delaunay 三角网或 Voronoi 图

Delaunay 三角网 是 Delaunay 三角化的结果。Delaunay 三角化是一种空间三角剖分算法。该算法的优点在于能避免三角网结果中出现“极瘦” 的三角形。

Voronoi 图 又称泰森多边形，由 Delaunay 三角网中连接两相邻点线段的中垂线组成的连续多边形组成。

sphtriangulate 读取一个或者多个经纬度坐标文件并进行球面 Delaunay 三角化。-Qv 选项可进步一生成 Voronoi 图。通常，多边形输出结果为闭合的并且可填充，使用

gmt sphtriangulate [ table ] [

table

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

可选选项

-A

计算球面三角的面积或多边形（见

从远程数据文件 hotspots.txt 中创建球面三角网，并绘制

gmt sphtriangulate @hotspots.txt -Qd -T | gmt plot -Rg -JG-120/-30/7i -Bafg -W3p -pdf map

三角化 testdata.txt 中的点，创建 Voroni 图并绘制

gmt sphtriangulate testdata.txt -Qv | gmt plot -Rg -JG30/30/6i -L -W1p -Bag -pdf testdata

基于包含多段记录的数据文件 globalnodes.txt 计算 Delaunay 三角网，并计算每个三角形的面积，输出到段信息中

gmt sphtriangulate globalnodes.txt -Qd -A > global_tri.txt

除上述示例外，GMT 官方示例 35 和 示例 48 也使用了该模块。

注意事项

多边形或三角形面积的计算使用了 STRIPACK 算法。该算法需要输入中不包含重复点。 blockmean 等模块可以将多个接近的点合并成单个点。sphtriangulate 的

Renka, R, J., 1997, Algorithm 772: STRIPACK: Delaunay Triangulation and Voronoi Diagram on the Surface of a Sphere, AMC Trans. Math. Software, 23(3), 416-434.

相关模块

triangiulate, sphdistance