为什么80%的码农都做不了架构师?>>>   

Visulalize Boost Voronoi in OpenSceneGraph

eryar@163.com

Abstract. One of the important features of the boost polygon library is the implementation of the generic sweepline algorithm to construct Voronoi diagrams of points and linear segments in 2D(developed as part of the Google Summer of Code 2010 program). Voronoi diagram data structure has applications in image segmentation, optical character recognition, nearest neighbor queries execution. It is closely related with the other computational geometry conectps: Delaunay triangulation, medial axis, straight skeleton, the largest empty circle. The paper focus on the usage of Boost.Polygon Voronoi Diagram and visualize it in OpenSceneGraph.

Key words. Voronoi, Boost.Polygon, C++, OpenSceneGraph, Visualization

1. Introduction

计算几何(Computational Geometry)作为一门学科,起源于20世纪70年代,经过近四十多年的发展,其研究内容不断扩大,涉及Voronoi图、三角剖分、凸包、直线与多边形求交、可见性、路径规划、多边形剖分等内容。据相关统计,在数以千计的相关文章中,约有15%是关于Voronoi图及其对偶(dual)图Delaunay三角剖分(Delaunay Triangulation)的研究。由于Voronoi图具有最近性、邻接性等众多性质和比较系统的理论体系,如今已经在计算机图形学、机械工程、地理信息系统、机器人、图像处理、大数据分析与处理、生物计算及无线传感网络等领域得到了广泛应用,同时也是解决碰撞检测、路径规划、可见性计算、骨架计算以及凸包计算等计算几何所涉及的其他问题的有效工具。

Voronoi图的起源最早可以追溯到17世纪。1644年,Descartes用类似Voronoi图的结构显示太阳系中物质的分布。数学家G.L. Dirichlet和M.G.Voronoi分别于1850年和1908年在他们的论文中讨论了Voronoi图的概念,所以Voronoi图又叫Dirichlet tessellation。在其他领域,这个概念也曾独立地出现,如生物学和生理学中称之为中轴变换(Medial Axis Transform)或骨架(Skeleton)。化学与物理学中称之为Wigner-Seitz Zones,气象学与地理学中称之为Thiessen多边形。Voronoi图最早由Thiessen应用于气象观测站中随机分布的研究。由于M.G. Voronoi从更通用的n维情况对其进行研究和定义,所以Voronoi图这个名称为大多数人所使用。

在路径规划、机械加工、模式识别、虚拟现实、生物计算等领域,将站点从离散点扩展到线段圆弧等生成Voronoi图的方式也是非常常见的。

目前可用于生成Voronoi图的库有一些,很多是开源库。像CGAL库、boost中也提供了生成Voronoi图的算法。本文根据Boost.Polygon中的Voronoi库,并用OpenSceneGraph显示出剖分结果。

2. Boost.Polygon Voronoi Diagram

Boost.Polygon库提供了构造Voronoi图的接口,可根据点集、线段集来生成Voronoi图,如下图所示:

Figure 2.1 Voronoi Diagram generated by Boost.Polygon Voronoi Algorithms

Boost.Polygon中的基于扫描线算法(sweep-line algorithm)Voronoi库可以实现如下功能:

v 输入数据可以是点集和线段;

v 算法的稳定性高及输出完整的拓朴信息;

v 可以控制输出的几何信息的精度;

计算几何方面以稳定性著称的CGAL库中的Voronoi算法只满足前两个功能。S-Hull库以上功能都没有很好的满足。下面是一些Boost.Polygon,CGAL,S-Hull库的对比数据:

Figure 2.2 Construction time for 10 random points

Figure 2.3 Construction time for 100 random points

Figure 2.4 Construction time for 1000 random points

Figure 2.5 Construction time for 10000 random points

Figure 2.6 Memory usage for 100000 random points

Figure 2.7 Logarithmic Execution Time

结论:

v 在输入上没有限制这点上CGAL要优于Boost.Polygon;

v Boost.Polygon Voronoi的稳定性要高于S-Hull;

v Boost.Polygon Voronoi和S-Hull的时间复杂度为N*log(N),而CGAL的不是;

v Boost.Polygon Voronoi的输出顶点的精度高于CGAL库;

v Boost.Polygon Voronoi的速度快;

v Boost.Polygon Voronoi根据10000个点或1000个线段来构造Voronoi的时间为0.02秒以内,所以可用来处理有实时性要求的场景;

3. Implementation

Boost.Polygon的Voronoi算法使用简单,只需要输入点集或线段集合,就可以直接构造出Voronoi图了。最简单的程序示例代码如下:

/*
*    Copyright (c) 2014 eryar All Rights Reserved.
*
*        File    : Main.cpp
*        Author  : eryar@163.com
*        Date    : 2014-05-06 18:28
*        Version : V1.0
*
*    Description : The Simplest example for boost voronoi library.
*      Key words : boost voronoi, C++
*
*/#include "boost/polygon/voronoi.hpp"using namespace boost::polygon;typedef int coordinate_type;
typedef point_data<coordinate_type> Point;
typedef voronoi_diagram<double> VD;int main(int argc, char* argv[])
{std::vector<Point> points;points.push_back(Point(0, 0));points.push_back(Point(1, 6));points.push_back(Point(-4, 5));points.push_back(Point(5, -1));points.push_back(Point(3, -11));points.push_back(Point(13, -1));VD vd;construct_voronoi(points.begin(), points.end(), &vd);return 0;
}

且Boost.Polygon的Voronoi算法遍历Voronoi边Edges,Voronoi单元cell,Voronoi顶点Vertex也很直接。如下代码所示为遍历所有边,来将剖分结果可视化:

/*
*    Copyright (c) 2014 eryar All Rights Reserved.
*
*        File    : Main.cpp
*        Author  : eryar@163.com
*        Date    : 2014-05-06 18:28
*        Version : V1.0
*
*    Description : VoronoiViewer for boost voronoi library visulization.
*      Key words : boost voronoi, C++, OpenSceneGraph
*
*/#include <osgViewer/Viewer>
#include <osgGA/StateSetManipulator>
#include <osgViewer/ViewerEventHandlers>#pragma comment(lib, "osgd.lib")
#pragma comment(lib, "osgDBd.lib")
#pragma comment(lib, "osgGAd.lib")
#pragma comment(lib, "osgViewerd.lib")#include "boost/polygon/voronoi.hpp"using namespace boost::polygon;typedef double coordinate_type;
typedef point_data<coordinate_type> Point;
typedef voronoi_diagram<coordinate_type> VD;osg::Node* BuildVoronoiDiagram(void)
{srand(static_cast<unsigned int> (time(NULL)));osg::ref_ptr<osg::Geode> theGeode = new osg::Geode();osg::ref_ptr<osg::Geometry> theLines = new osg::Geometry();osg::ref_ptr<osg::Vec3Array> theVertices = new osg::Vec3Array();VD vd;std::vector<Point> thePoints;// Add points for the Voronoi Diagram.for (int i = 0; i < 100; ++i){int x = rand() % 100;int y = rand() % 100;thePoints.push_back(Point(x, y));// Display the site of the Voronoi Diagram.theVertices->push_back(osg::Vec3(x - 1, 0.0, y));theVertices->push_back(osg::Vec3(x + 1, 0.0, y));theVertices->push_back(osg::Vec3(x, 0.0, y - 1));theVertices->push_back(osg::Vec3(x, 0.0, y + 1));}construct_voronoi(thePoints.begin(), thePoints.end(), &vd);// Visualize the edge of the Voronoi Diagram.// Traversing Voronoi edges using edge iterator.for (VD::const_edge_iterator it = vd.edges().begin(); it != vd.edges().end(); ++it){if (it->is_primary()){if (it->is_finite()){theVertices->push_back(osg::Vec3(it->vertex0()->x(), 0.0, it->vertex0()->y()));theVertices->push_back(osg::Vec3(it->vertex1()->x(), 0.0, it->vertex1()->y()));}else{Point p1 = thePoints[it->cell()->source_index()];Point p2 = thePoints[it->twin()->cell()->source_index()];Point origin;Point direction;coordinate_type koef = 1.0;origin.x((p1.x() + p2.x()) * 0.5);origin.y((p1.y() + p2.y()) * 0.5);direction.x(p1.y() - p2.y());direction.y(p2.x() - p1.x());if (it->vertex0() == NULL){theVertices->push_back(osg::Vec3(origin.x() - direction.x() * koef, 0.0,origin.y() - direction.y() * koef));}else{theVertices->push_back(osg::Vec3(it->vertex0()->x(), 0.0, it->vertex0()->y()));}if (it->vertex1() == NULL){theVertices->push_back(osg::Vec3(origin.x() + direction.x() * koef, 0.0,origin.y() + direction.y() * koef));}else{theVertices->push_back(osg::Vec3(it->vertex1()->x(), 0.0, it->vertex1()->y()));}}}}theLines->setVertexArray(theVertices);// Set the colors.osg::ref_ptr<osg::Vec4Array> theColors = new osg::Vec4Array();theColors->push_back(osg::Vec4(1.0f, 1.0f, 0.0f, 1.0f));theLines->setColorArray(theColors);theLines->setColorBinding(osg::Geometry::BIND_OVERALL);// Set the normal.osg::ref_ptr<osg::Vec3Array> theNormals = new osg::Vec3Array();theNormals->push_back(osg::Vec3(0.0f, -1.0f, 0.0f));theLines->setNormalArray(theNormals);theLines->setNormalBinding(osg::Geometry::BIND_OVERALL);theLines->addPrimitiveSet(new osg::DrawArrays(osg::PrimitiveSet::LINES, 0, theVertices->size()));theGeode->addDrawable(theLines);return theGeode.release();
}int main(int argc, char *argv[])
{osgViewer::Viewer theViewer;theViewer.setSceneData(BuildVoronoiDiagram());theViewer.addEventHandler(new osgGA::StateSetManipulator(theViewer.getCamera()->getOrCreateStateSet()));theViewer.addEventHandler(new osgViewer::StatsHandler);theViewer.addEventHandler(new osgViewer::WindowSizeHandler);return theViewer.run();
}

绘制Voronoi的边时,当边是有限的finite时,直接可以画出直线;当边是infinite时,根据定义计算出了无界边的方向。显示结果如下图所示:

Figure 3.1 Construct Voronoi Diagram by 10 random points by Boost.Polygon

Figure 3.2 Construct Voronoi Diagram by 100 random points by Boost.Polygon

4. Conclusion

Boost.Polygon中的Voronoi图算法稳定性及性能较高,且可以根据站点查找相关的拓朴信息,如根据站点查找Voronoi单元等;惟一不足的就是默认只处理整数点集。

当输入有线段时,生成的Voronoi图中有曲线,曲线的绘制可参考相关例子实现。

Boost.Polygon中的Voronoi算法都以模板实现,编译时只需要包含相关的头文件即可,不依赖其他库,使用还是很方便的。

5. References

1. http://www.boost.org/

2. Boost.Polygon,http://www.boost.org/doc/libs/1_55_0/libs/polygon/doc/index.htm

3. Voronoi Basic Tutorial,\boost_1_54_0\libs\polygon\doc\voronoi_basic_tutorial.htm

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

5. 杨承磊, 吕琳, 杨义军, 孟祥旭. Voronoi图及其应用. 清华大学出版社. 2013

PDF Version and Source Code: Visulalization Boost Voronoi in OpenSceneGraph

转载于:https://my.oschina.net/eryar/blog/347806

Visulalize Boost Voronoi in OpenSceneGraph相关推荐

  1. 跨平台C++开发环境搭建——GNU+STL+Boost+wxWidget+OpenSceneGraph

    学了好长时间的C++了,鉴于对标准和设计模型模式的推崇,一直没有找到合适的GUI语言,也一直苦于GUI开发时我还要借助于Delphi,但是现在不一样了,花了一个周末整个跨平台的C++开发环境,编译器采 ...

  2. boost::polygon模块voronoi相关的测试程序

    boost::polygon模块voronoi相关的测试程序 实现功能 C++实现代码 实现功能 boost::polygon模块voronoi相关的测试程序 C++实现代码 #include < ...

  3. boost一些常用库(数据结构,迭代器,算法及字符串) 简单介绍

    原文地址:http://blog.csdn.net/jinzhuojun/article/details/51923337 C++强大的重要原因之一就是它的库.而boost就是这样一个包罗万象的C++ ...

  4. 所有的 Boost 库文档的索引

    入门 工具 网站 新闻 社区 常见问题 更多的信息 按字母顺序列出的库 按类别列出的库 算法 破碎的编译器的解决方法 并发编程 容器 正确性和测试 数据结构 特定于域的 函数对象和高阶编程 泛型编程 ...

  5. PCL(Point Cloud Library)的第三方库简介(boost,eigen,flann,vtk,qhull)

    PCL因为融合了大量的第三方开源库,导致学习成本升高~在学习之前我们不妨了解一下这些库都是干嘛的,以便有的放矢,在之后更好的使用 boost: C++的标准库的备用版,擅长从数学库到智能指针,从模板元 ...

  6. Linux下安装 boost 库

    1. 先去官网下载压缩包: https://www.boost.org/ 2. 解压后cd 进入根目录,然后执行: ./bootstrap.sh 3. 安装(时间会很长): ./b2 install ...

  7. Mysql依赖库Boost的源码安装,linux下boost库的安装

    boost'准标准库'安装过程. 安装的是boost_1_60_0. (1)首先去下载最新的boost代码包,网址www.boost.org. (2)进入到自己的目录,解压: bzip2 -d boo ...

  8. 函数指针amp;绑定: boost::functoin/std::function/bind

    see link: https://isocpp.org/wiki/faq/pointers-to-members function vs template: http://stackoverflow ...

  9. Boost 1.53.0 发布,可移植的C++标准库

    Boost 1.53.0 发布了,包含了 5 个新的库,修复了一些安全漏洞以及 Boost.Locale 组件的 bug . 新增的 5 个库包括: Boost.Atomic Boost.Corout ...

最新文章

  1. TypeScript 从听说到入门(上篇)
  2. 漫画:据说很多搞软件的羡慕硬件工程师
  3. 深入分析SpringBoot源码如何内嵌Tomcat容器?
  4. 层次聚类(cluster.AgglomerativeClustering)并通过TSNE进行可视化分析
  5. 推荐系统——开源代码
  6. Prthon编码和解码/is 和 ==
  7. 零基础入门Python3-set集合详解
  8. 海奥华预言--第一章 神秘邀请
  9. extjs4.2 弹出注册用户界面实例
  10. 2021年中国微粉磨料市场趋势报告、技术动态创新及2027年市场预测
  11. oracle外部结合,浅谈Oracle外部文件
  12. linkbox php,win10 docker-toolsbox 搭建php开发环境
  13. c语言免杀程序源码,ghost源码免杀教程 步
  14. 洛谷1498-谢尔宾斯基三角形-python-(递归)
  15. iTween使用小坑
  16. 0xC0000005:读取位置0x00000000时发生访问冲突
  17. cati服务器授权信息无效,CATI基础知识介绍(四)
  18. 交互媒体专题设计------《The Wiley Handbook of Human Computer Interaction》
  19. Python爬虫:ZzzFun动漫视频网
  20. Android模拟电话,向Android模拟器打电话和发短消息

热门文章

  1. 开发一个属于自己的app
  2. java 判断网络类型_javaexcel判断类型
  3. Bentley 软件公司将以约 10.5 亿美元的价格收购全球领先的地球科学三维建模软件提供商 Seequent
  4. 计算机专业毕设评阅人评语,毕业论文指导教师评语与评阅人评语写法
  5. Multisim调幅和解调电路仿真实验
  6. AD根据Primitives进行敷铜
  7. 【Linux(一)】最新VMware虚拟机下载与安装
  8. install4j的使用
  9. 解决 No converter found capable of converting from type [java.lang.String] to type ... 的问题
  10. python-docx对Word文档的指定位置(批量)插入图片