枚举多边形质心应具有的属性是一个好主意。这是我的标准:
(a)这是多边形内部的属性(而不是顶点或边)。因此,通过插入附加顶点将边一分为二不应更改质心的位置。请注意,Jenness对质心的定义在此标准上失败,因为质心的位置将取决于多边形如何划分为三角形。
(b)稍微干扰多边形的形状应将质心稍微移动一点。这里有必要对多边形的整体范围施加限制(例如,限制到一个半球)。如果没有此限制,则很容易构造出质心突然通过顶点的轻微移动而摆动到地球另一侧的情况。此条件不包括要求质心位于多边形内部的方法。
(c)对于小多边形,应减少到质心的平面定义。
以下是满足这些条件的两种方法:
(1)在三个维度上计算椭球面的质心,然后投影回椭球面(沿椭球的法线)。很大的优势:可以通过将多边形分解为更简单的形状来计算质心。
(2)质心是到多边形内部所有点的RMS测地线距离最小的点。见巴斯和菲尔莫尔,“球形平均值和应用球面样条插值”,在图形ACM交易20,95-126(2001)。巨大的优势:所得的点不依赖于R 3中如何嵌入表面。
不幸的是,这些定义都不容易直接实施。 但是,第一种方法可以简单地用于球体。最好使用的“基本”区域是由多边形的边缘,通过该边缘的端点的两个子午线以及赤道所界定的四边形。整个多边形的结果需要对边上的贡献求和。(如果多边形环绕一个极点,则需要采取其他步骤。)
假定边缘的端点为(φ 1,λ 1)和(φ 2,λ 2)。让所述边缘的通过α的方位角和端点1
和α 2。假设球体的半径为1,则四边形的面积为
甲 =α 2 - α 1
= 2黄褐色-1
[黄褐色½(λ 2 - λ 1)罪½(φ 2 +φ 1)/余弦½(φ 2 +φ 1)]
(由于贝塞尔的缘故,该面积公式在数值上要比常用的三角形面积的L'Huilier公式好。)
该四边形的质心分量由下式给出
2 甲 ⟨ X ⟩=φ 2罪(λ 2 - λ 0) - φ 1罪(λ 1 - λ 0)
2 甲 ⟨ Ŷ ⟩= COSα 0(σ 2 - σ 1) - (φ 2 COS(λ 2 - λ 0) - φ 1 COS(λ 1 - λ 0))
2 甲 ⟨ ž ⟩=(λ 2 - λ 1) -罪α 0(σ 2 - σ1)
其中,σ 2 - σ 1是边缘的长度,λ 0和α 0是经度和它穿过赤道边缘的方位,并且
X和ÿ轴被定向为使得所述赤道交叉是在X = 1,y =0。(当然,z是穿过极点的轴。)