C#如何实现3D数组中三点定义的任意2D平面计算

核心实现步骤 #

1. 基于输入三点推导平面方程 #

给定三个不共线的三维点P0/P1/P2，先计算平面内两个边向量v1 = P1-P0、v2 = P2-P0，对两个向量做叉乘得到平面法向量n = Cross(v1, v2)，最终整理为平面标准式A*x + B*y + C*z + D = 0，其中A/B/C对应法向量n的三个分量，D = -Dot(n, P0)。

所有几何计算统一使用double类型，避免float精度不足导致的截面空洞、条纹问题。

2. 计算平面与三维数组包围盒的相交区域 #

三维数组本质是轴对齐立方体，边界范围固定：

• x维度范围：[0, xMax]，xMax为单张图像宽度-1
• y维度范围：[0, yMax]，yMax为单张图像高度-1
• z维度范围：[0, zMax]，zMax为图像栈总层数-1
  逐次计算平面与立方体12条棱的交点，过滤掉超出棱范围的无效点，最终得到的有效交点会构成一个凸多边形，这就是平面在三维数组内的全部有效取值范围——不管用户输入的三个点是否落在数组边界上，都能覆盖平面穿过数组的所有区域。

3. 构建平面局部二维坐标系，确定输出截面尺寸 #

平面本身是二维结构，需要建立局部正交坐标系来确定输出图像的宽高：

• 取向量v1单位化后作为平面局部U轴（对应输出图像宽度方向）
• 用法向量n叉乘U轴得到V轴方向，单位化后作为输出图像高度方向
• 把第2步得到的所有凸多边形交点投影到U-V坐标系，分别计算U方向最小值uMin、最大值uMax，V方向最小值vMin、最大值vMax
• 输出图像的像素宽度为(int)Math.Ceiling(uMax - uMin)，像素高度为(int)Math.Ceiling(vMax - vMin)，如果z轴和xy轴的物理分辨率不一致，需要按实际间距缩放尺寸避免拉伸。

4. 逐点采样取值 #

遍历输出图像的每个像素坐标(i,j)，通过U/V轴基向量反算对应的三维空间坐标：

pt3D.X = P0.X + (uMin + i) * uAxis.X + (vMin + j) * vAxis.X
pt3D.Y = P0.Y + (uMin + i) * uAxis.Y + (vMin + j) * vAxis.Y
pt3D.Z = P0.Z + (uMin + i) * uAxis.Z + (vMin + j) * vAxis.Z

由于斜平面对应的三维坐标大多不是整数索引，需要对该点做三线性插值：取坐标周围8个相邻整数索引的体素值，按距离加权计算得到当前点的像素值。如果点坐标因浮点误差超出数组边界，直接赋值为背景值即可。

优化点 #

• 加快速分支：如果法向量某两个分量的绝对值小于1e-6，说明平面平行于某一坐标平面（比如法向量只有Z分量就是平行于xy面的普通切片），直接取对应层/行/列的数值即可，不需要走插值逻辑，性能提升明显。
• 大体积栈优化：不需要遍历整个三维数组，可以预计算平面和每个z层切片的交线，仅遍历交线覆盖的像素，能把计算量降低一个数量级。

核心代码参考 #

public struct Point3D
{
    public double X, Y, Z;
    public Point3D(double x, double y, double z)
    {
        X = x;
        Y = y;
        Z = z;
    }
    public static Point3D operator -(Point3D a, Point3D b) 
        => new Point3D(a.X - b.X, a.Y - b.Y, a.Z - b.Z);
    
    public static Point3D Cross(Point3D a, Point3D b) 
        => new Point3D(
            a.Y * b.Z - a.Z * b.Y,
            a.Z * b.X - a.X * b.Z,
            a.X * b.Y - a.Y * b.X
        );
    
    public static double Dot(Point3D a, Point3D b) 
        => a.X * b.X + a.Y * b.Y + a.Z * b.Z;
}

public (double A, double B, double C, double D) GetPlaneEquation(Point3D p0, Point3D p1, Point3D p2)
{
    Point3D v1 = p1 - p0;
    Point3D v2 = p2 - p0;
    Point3D normal = Point3D.Cross(v1, v2);
    double D = -Point3D.Dot(normal, p0);
    return (normal.X, normal.Y, normal.Z, D);
}

内容的提问来源于stack exchange，提问作者cubes247