读取STL文件后，将三角面片的顶点保存在osg::kdTree中，如何利用osg::kdTree判断是否存在非闭合的三角面片

在OpenSceneGraph (osg) 中，处理STL文件并将其内容存储在`osg::Geometry`对象里通常是第一步。然后，你可以创建一个`osg::KDTree`，这是一种空间划分的数据结构，用于快速查找三维空间内的点是否存在于已知的几何体范围内。

要检查三角面片是否闭合，你需要遍历每个面（即三角形），计算其三个顶点形成的向量围成的角度。如果某个面的任意两个顶点之间的角度小于180度，那么这个面就是非闭合的。这是因为在一个封闭的三角形中，所有边的方向应该指向内部。

以下是大致步骤：

1. **加载STL文件**：
- 通过`osgDB::readNodeFile`函数读取STL文件，并获取包含三角面片的`osg::Geode`或`osg::Geometry`。

2. **提取顶点信息**：
- 使用`osg::Geometry::getVertexArray()`和`getNormalArray()`获取顶点和法线数据。

3. **构建KDTree**：
- 创建一个`osg::ref_ptr<osg::KDTree>`实例，并设置其范围以包括所有几何体的顶点。

4. **遍历三角面片**：
- 对于每个三角形的三个顶点，使用`KDTree::containsPoint`来检查它是否在树内。若不在，可能是非闭合的。

5. **检查闭合性**：
- 计算每个三角形的三个顶点形成的角度，如果最大角小于180度，则面是开放的。

// 假设我们已经有了一个名为"geometry"的osg::Geometry对象
for (int i = 0; i < geometry.getNumTriangles(); ++i) {
    const osg::Vec3* v0 = geometry.getVertexArray()[i * 3];
    const osg::Vec3* v1 = geometry.getVertexArray()[i * 3 + 1];
    const osg::Vec3* v2 = geometry.getVertexArray()[i * 3 + 2];

    // 检查v0-v1、v1-v2 和 v2-v0 的角度
    float angle1 = osg::Math::degrees(osg::calculateAngle(v1, v0, v2));
    float angle2 = osg::Math::degrees(osg::calculateAngle(v2, v1, v0));

    if (angle1 > 180 || angle2 > 180) { // 如果有一个角度大于180，说明面是闭合的
        continue;
    }

    // 使用KDTree检查是否非凸
    if (!kdTree.containsPoint(osg::Vector3(v0.x, v0.y, v0.z))) {
        std::cout << "Non-closed triangle found at (" << v0.x << ", " << v0.y << ", " << v0.z << ")" << std::endl;
    }
}