Visualize Surface by Delaunay Triangulator

特意加上了材质和光照效果。离散的曲面效果如下图所示：

Figure 2.1 Shaded Surfaces in OpenSceneGraph

Figure 2.2 Mesh of the Surfaces in OpenSceneGraph

为了显示出光照效果，使用了osgUtil::SmoothingVisitor来自动计算三角网格的平均法向。从上图可知，这种方式是将曲面的参数空间均匀剖分得到的曲面，当剖分得密时，数据量大，但显示得逼真。剖分疏时，数据量小，显示失真。

利用特征敏感网格重剖技术，可以使用较少的三角面片来比较精确地表示几何模型，如下图所示：（图片来源： http://cg.cs.tsinghua.edu.cn/papers/TVCG2007featuresensitive.pdf）

Figure 2.3 Feature sensitive remeshing ( http://cg.cs.tsinghua.edu.cn/)

3. Holes in the Surface

当曲面上有开孔时，开孔的信息可以通过WIRE来获得。对于组成开孔的WIRE的每条EDGE，可以通过曲面上的曲线PCurve来将开孔的参数统一到曲面的UV参数空间。

Figure 3.1 Hole in Parametric UV space

在OpenSceneGraph中实现代码如下所示：

osg::Node* BuildSurface(const Handle_Geom_Surface &theSurface, const Handle_Geom2d_Curve &thePCurve)
{
    osg::ref_ptr<osg::Geode> aGeode = new osg::Geode();
    osg::ref_ptr<osg::Geometry> aGeometry = new osg::Geometry();

    osg::ref_ptr<osg::Vec3Array> aUVPoints = new osg::Vec3Array();
    osg::ref_ptr<osg::Vec3Array> aPoints = new osg::Vec3Array();
    osg::ref_ptr<osg::Vec3Array> aBounds = new osg::Vec3Array();

    osg::ref_ptr<osgUtil::DelaunayTriangulator> dt = new osgUtil::DelaunayTriangulator();
    osg::ref_ptr<osgUtil::DelaunayConstraint> dc = new osgUtil::DelaunayConstraint();

    // build triangulation for the parametric space.
    Standard_Real u1 = 0.0;
    Standard_Real u2 = 0.0;
    Standard_Real v1 = 0.0;
    Standard_Real v2 = 0.0;

    theSurface->Bounds(u1, u2, v1, v2);

    Precision::IsNegativeInfinite(u1) ? u1 = -1.0: u1;
    Precision::IsPositiveInfinite(u2) ? u2 =  1.0: u2;
    Precision::IsNegativeInfinite(v1) ? v1 = -1.0: v1;
    Precision::IsPositiveInfinite(v2) ? v2 =  1.0: v2;

    // tesselate the parametric space.
    Standard_Integer aStep = 30;
    Standard_Real uDelta = (u2 - u1) / aStep;
    Standard_Real vDelta = (v2 - v1) / aStep;

    for (Standard_Integer i = 0; i <= aStep; ++i)
    {
        for (Standard_Integer j = 0; j <= aStep; ++j)
        {
            Standard_Real u = u1 + i * uDelta;
            Standard_Real v = v1 + j * vDelta;

            aUVPoints->push_back(osg::Vec3(u, v, 0.0));
        }
    }

    Standard_Real pDelta = (thePCurve->LastParameter() - thePCurve->FirstParameter()) / aStep;
    for (Standard_Integer c = 0; c <= aStep; ++c)
    {
        gp_Pnt2d p = thePCurve->Value(thePCurve->FirstParameter () + c * pDelta);

        aBounds->push_back(osg::Vec3(p.X(), p.Y(), 0.0));
    }

    dc->setVertexArray(aBounds);
    dc->addPrimitiveSet(new osg::DrawArrays(osg::PrimitiveSet::LINE_LOOP, 0, aBounds->size()));

    // triangulate the parametric space.
    dt->addInputConstraint(dc);
    dt->setInputPointArray(aUVPoints);
    dt->triangulate();
    dt->removeInternalTriangles(dc);

    for (osg::Vec3Array::const_iterator j = aUVPoints->begin(); j != aUVPoints->end(); ++j)
    {
        // evaluate the point on the surface
        gp_Pnt aPoint = theSurface->Value((*j).x(), (*j).y());

        aPoints->push_back(osg::Vec3(aPoint.X(), aPoint.Y(), aPoint.Z()));
    }

    //aGeometry->setVertexArray(aUVPoints);
    aGeometry->setVertexArray(aPoints);
    aGeometry->addPrimitiveSet(dt->getTriangles());

    aGeode->addDrawable(aGeometry);

    // use smoothing visitor to set the average normals
    osgUtil::SmoothingVisitor sv;
    sv.apply(*aGeode);

    // set material for the surface
    osg::ref_ptr<osg::StateSet> aStateSet = aGeode->getOrCreateStateSet();
    osg::ref_ptr<osg::Material> aMaterial = new osg::Material();

    aMaterial->setDiffuse(osg::Material::FRONT, osg::Vec4(1.0f, 1.0f, 0.0f, 1.0f));
    aMaterial->setSpecular(osg::Material::FRONT, osg::Vec4(1.0f, 1.0f, 1.0f, 1.0f));
    aMaterial->setShininess(osg::Material::FRONT, 100.0f);

    aStateSet->setAttribute(aMaterial);

    return aGeode.release();
}

开孔的实现主要也是用osgUtil::DelaunayConstraint，将孔中的三角形去除。如下图所示为在球面和锥面开孔：

Figure 3.2 Holes on Sphere and Cone Surface

Figure 3.3 Mesh for Sphere and Cone with holes

由上图可知，对于有拓朴结构的三维模型数据，可以由WIRE得到组成孔的每条边Edge，根据Edge中PCurve可以找到对应的曲面。通过PCurve将开孔数据统一到参数空间中。剖分完带孔的参数空间，再映射回三维空间就得到开孔的曲面了。

4. Conclusion

原来一直百思不得其解的问题现在已经豁然开朗，对三维模型的可视化有了一定的理解。借助于这些开源库，对相关知识的学习要轻松许多，可以将许多国产教材上断续的理论知识与实践衔接起来。

5. Acknowledgement

感谢OpenCascade和OpenSceneGraph的开放和分享，让学习变得轻松有趣。

6. Reference

1. 孙家广, 胡事民等. 计算机图形学. 清华大学出版社. 2000

2. 赵罡, 穆国旺, 王拉柱译. 非均匀有理B样条. 清华大学出版社. 2010

3. Jonathan R. Shewchuk. Triangle: http://www.cs.cmu.edu/~quake/triangle.html

4. 汪嘉业 王文平 屠长河 杨承磊. 计算几何及应用.  科学出版社. 2011

5. 王成恩. 面向科学计算的网格划分与可视化技术. 科学出版社. 2011

6. 周培德. 计算几何－算法设计与分析. 清华大学出版社. 2008

7. http://cg.cs.tsinghua.edu.cn/

8. Mesh Algorithm in OpenCascade：

http://www.cppblog.com/eryar/archive/2014/04/06/206484.html