当我们逆时针从V1到V2时,如何确定V3是否在V1和V2之间?

crossProds = [V1(1)*V2(2)-V1(2)*V2(1) ...
              V1(1)*V3(2)-V1(2)*V3(1) ...
              V3(1)*V2(2)-V3(2)*V2(1)];
if (all(crossProds >= 0) || ...
    (crossProds(1) < 0) && ~all(crossProds(2:3) < 0)),
  disp("V3 is between V1 and V2");
else
  disp("out of the interval");
end

解释:

第1版 和 。如果以下两者之间的逆时针角度 第1版 和 第2版 第1版 和 第2版 和 也大于或等于零。这解释了第一个逻辑检查:

all(crossProds >= 0)

第1版 和 第2版 第1版 和 第2版 和 第1版 和 第2版

(crossProds(1) < 0) && ~all(crossProds(2:3) < 0)

short circuit operators

a1 <= a2 < a1 + 2*pi
a1 <= a3 < a1 + 2*pi

:

R1 -= R3;
R2 -= R3;

if (R1 < 0) R1 += 2 * PI;
if (R2 <= 0) R2 += 2 * PI;

return (r1 < r2);

逆时针

要测试三角形的缠绕,只需检查顶点的二维叉积的符号即可。

int IsBetween (vector v1, vector v2, vector v3)
{
  float winding1 = (v1.x * v3.y - v1.y * v3.x);
  float winding2 = (v3.x * v2.y - v3.y * v2.x);

  // this test could be exactly the wrong way around. This depends
  // on how you define your coordinate system (e.g. is Y going up or down?)

  if ((winding1 <0) && (winding2 < 0))
  {
    printf ("V3 is between them\n");
  }
  else
  {
    printf ("it's not\n");
  }
}