Question

我想计算具有由三角形构成的表面的3D网格物体的体积。

Answer 1

Reading this paper，这实际上是一个非常简单的计算。

诀窍是计算四面体的有符号体积 - 基于您的三角形并在原点处加满。音量的符号来自您的三角形是否指向原点方向。（三角形的法线本身取决于顶点的顺序，这就是为什么你没有看到它在下面明确引用的原因。）

这一切归结为以下简单功能：

public float SignedVolumeOfTriangle(Vector p1, Vector p2, Vector p3) {
    var v321 = p3.X*p2.Y*p1.Z;
    var v231 = p2.X*p3.Y*p1.Z;
    var v312 = p3.X*p1.Y*p2.Z;
    var v132 = p1.X*p3.Y*p2.Z;
    var v213 = p2.X*p1.Y*p3.Z;
    var v123 = p1.X*p2.Y*p3.Z;
    return (1.0f/6.0f)*(-v321 + v231 + v312 - v132 - v213 + v123);
}

然后是一个计算网格体积的驱动程序：

public float VolumeOfMesh(Mesh mesh) {
    var vols = from t in mesh.Triangles
               select SignedVolumeOfTriangle(t.P1, t.P2, t.P3);
    return Math.Abs(vols.Sum());
}

Answer 2

Yip Frank Kruegers的回答很有效。如果你有可用的向量函数，你也可以使用它：

    public static float SignedVolumeOfTriangle(Vector p1, Vector p2, Vector p3)
    {
        return p1.Dot(p2.Cross(p3)) / 6.0f;
    }

编辑..添加了impl。如果您不确定，请使用Dot（）和Cross（）。大多数数学库都会有这些。如果您使用的是WPF，则它们将作为Vector3D类的静态方法实现。

    public class Vector
    {
        ... 

        public float Dot(Vector a)
        {
            return this.X * a.X + this.Y * a.Y + this.Z * a.Z;
        }

        public Vector Cross(Vector a)
        {
            return new Vector(
              this.Y * a.Z - this.Z * a.Y,
              this.Z * a.X - this.X * a.Z,
              this.X * a.Y - this.Y * a.X
            );
        }
        ...
    }

Answer 3

GNU Triangulated Surface Library can do this for you。请记住，表面必须关闭。对于不少3D模型来说情况并非如此。

如果你想自己实现它，你可以先看看他们的代码。

Answer 4

上述方法适用于“简单”对象（没有交叉/重叠三角形），如球体四面体等。对于更复杂的形状，一个好主意可能是分割网格（关闭它）并分别计算每个网段的体积。希望这会有所帮助。

如何计算表面由三角形组成的3D网格对象的体积

4 个答案: