1. 思路

这个问题其实涉及到OSG中的两个问题:多边形分格化和几何图元遍历。

1) 多边形分格化

在OpenGL/OSG中,由于效率的原因,默认是直接显示的简单的凸多边形。如果直接强行显示凹多边形,渲染结果是不确定的。所以对于复杂的凹多边形,需要将其分解成简单的凸多边形,这个过程就是多边形分格化。在OSG中是通过osgUtil::Tessellator类来实现多边形分格化的。

2) 几何图元遍历

对于二维的凹多边形,可以有办法计算其面积。但是对于三维空间的凹多边形,计算其面积却很困难。这是因为三维空间凹多边形甚至都有可能不是共面的。而我们知道,任何复杂的图形都是通过分解成三角形进行绘制的,只要获取分解成的三角形,计算其面积并相加(空间三角形的面积计算比较简单),就可以得到凹多边形的总面积。

在OSG中提供了一个用来访问图元的类:osg::PrimitiveFunctor,其继承类osg::TriangleFunctor可以获取其三角面图元。几何体类osg::Geometry提供了遍历几何图元的访问器接口。

2. 实现

其具体实现如下。注意在查找多边形分格化的资料的时候,提到了环绕数和环绕规则的概念。在OSG里面也有相应的参数设置。可惜这一段没有看明白,只能根据仿照例子来设置了。

#include <iostream>
#include <Windows.h> #include <osgViewer/Viewer>
#include <osgDB/ReadFile> #include <osgUtil/Tessellator>
#include <osg/TriangleFunctor> using namespace std;
using namespace osg; osg::ref_ptr<osg::Geometry> redPolygon; //计算空间三角形的面积
double CalTriangleArea(const osg::Vec3& a, const osg::Vec3& b, const osg::Vec3& c)
{
double area = 0; double side[3];//存储三条边的长度; side[0] = sqrt(pow(a.x() - b.x(), 2) + pow(a.y() - b.y(), 2) + pow(a.z() - b.z(), 2));
side[1] = sqrt(pow(a.x() - c.x(), 2) + pow(a.y() - c.y(), 2) + pow(a.z() - c.z(), 2));
side[2] = sqrt(pow(c.x() - b.x(), 2) + pow(c.y() - b.y(), 2) + pow(c.z() - b.z(), 2)); //不能构成三角形;
if (side[0] + side[1] <= side[2] || side[0] + side[2] <= side[1] || side[1] + side[2] <= side[0]) return area; //利用海伦公式。s=sqr(p*(p-a)(p-b)(p-c));
double p = (side[0] + side[1] + side[2]) / 2; //半周长;
area = sqrt(p*(p - side[0])*(p - side[1])*(p - side[2])); return area;
} //三角面片访问器
struct TriangleAreaFunctor
{
TriangleAreaFunctor()
{
sumArea = new double;
} ~TriangleAreaFunctor()
{
if (sumArea)
{
delete sumArea;
sumArea = nullptr;
}
} void operator() (const osg::Vec3& v1, const osg::Vec3& v2, const osg::Vec3& v3) const
{
*sumArea = *sumArea + CalTriangleArea(v1, v2, v3);
} double GetSumArea()
{
return *sumArea;
} protected:
double *sumArea = nullptr;
}; //
ref_ptr<Geode> createPolygon()
{
const float wall[6][3] =
{
{ -115.54f, 70.873f, -118.952f},
{ -111.516f, 70.7189f, -71.8492f },
{ -88.5345f, 70.8667f, -86.3565f },
{ -64.9495f, 71.8231f, -53.6525f },
{ -52.9755f, 69.028f, -129.093f },
{ -89.2272f, 71.1478f, -105.434f }
}; //
ref_ptr<Geode> geode = new Geode();
redPolygon = new osg::Geometry; //
osg::ref_ptr<osg::Vec3Array> redVex = new osg::Vec3Array;
redPolygon->setVertexArray(redVex); for (int i = 0; i< 6; i++)
{
redVex->push_back(osg::Vec3(wall[i][0], wall[i][1], wall[i][2]));
} redPolygon->addPrimitiveSet(new osg::DrawArrays(osg::PrimitiveSet::POLYGON, 0, 6)); //设置颜色数组
osg::ref_ptr<osg::Vec4Array> redColors = new osg::Vec4Array;
redColors->push_back(osg::Vec4(1.0, 0.0, 0.0, 0.5));
redPolygon->setColorArray(redColors);
redColors->setBinding(osg::Array::BIND_PER_PRIMITIVE_SET); //如果需要透明,则加入这个
redPolygon->getOrCreateStateSet()->setMode(GL_BLEND, osg::StateAttribute::ON);
redPolygon->getOrCreateStateSet()->setRenderingHint(osg::StateSet::TRANSPARENT_BIN); //创建分格化对象(支持凹多边形)
osg::ref_ptr<osgUtil::Tessellator> tscx = new osgUtil::Tessellator;
//设置分格类型为几何体
tscx->setTessellationType(osgUtil::Tessellator::TESS_TYPE_GEOMETRY);
//设置只显示轮廓线为false。设置环绕规则,这里不太懂
tscx->setWindingType(osgUtil::Tessellator::TESS_WINDING_ODD);
//使用分格化
tscx->retessellatePolygons(*(redPolygon.get())); geode->addDrawable(redPolygon); return geode;
} int main()
{
//
ref_ptr<Group> root = new Group();
root->getOrCreateStateSet()->setMode(GL_LIGHTING, osg::StateAttribute::OFF | osg::StateAttribute::OVERRIDE); //关闭默认光照
root->addChild(createPolygon()); //
osgViewer::Viewer viewer;
viewer.setSceneData(root);
viewer.setUpViewInWindow(100, 100, 800, 600);
viewer.run(); osg::TriangleFunctor<TriangleAreaFunctor> tf;
redPolygon->accept(tf);
cout << "面积:" << tf.GetSumArea() << endl; return 0;
}

渲染结果如下:

3. 参考

  1. OSG学习笔记(三)之如何将非三角面转换为三角面
  2. osg几何体的图元的遍历
  3. OSG计算并绘制模型中每一个三角面片的法向量
  4. OSG(OpenSceneGraph)基础学习9:OSG多边形分格化

OSG绘制空间凹多边形并计算其面积的更多相关文章

  1. opencv学习之路(27)、轮廓查找与绘制(六)——外接圆、椭圆拟合、逼近多边形曲线、计算轮廓面积及长度、提取不规则轮廓

    一.最小外接圆 #include "opencv2/opencv.hpp" #include<iostream> using namespace std; using ...

  2. google map 计算地图面积方法

    花了几个小时把js的google计算地图面积的算法改成了c# 的. class Program { static void Main(string[] args) { // a = new qq.ma ...

  3. ArcGIS应用——四种计算图斑面积的方法

    ArcGIS中有多种方法可计算出图斑面积,本文总结了四种方法,是否可堪称史上最全? 1.计算几何 本人认为这是最适合非专业人士的方法,直接利用ArcGIS中的计算几何功能进行计算. a.首先添加一do ...

  4. OpenJudge计算概论-计算三角形面积【海伦公式】

    /*============================================== 计算三角形面积 总时间限制: 1000ms 内存限制: 65536kB 描述 平面上有一个三角形,它的 ...

  5. 有一台机器,上面有m个储存空间。然后有n个请求,第i个请求计算时需要占 R[i]个空间,储存计算结果则需要占据O[i]个空间(据O[i]个空间(其中O[i]<R[i])。问怎么安排这n个请求的顺序,使

    有一台机器,上面有m个储存空间.然后有n个请求,第i个请求计算时需要占 R[i]个空间,储存计算结果则需要占据O[i]个空间(据O[i]个空间(其中O[i]<R[i]).问怎么安排这n个请求的顺 ...

  6. [math] 绘制空间几何体的直观图

    这么多年,一直凭着从天而降的神来之灵感画着立体图. 而今才知道在二维平面上绘制空间几何体的直观图也是有方法的.叫做“画法几何” 1. 斜二测图 就是倾斜y轴,使y轴与x轴成45度的夹角.见: http ...

  7. 计算平面面积和斜面面积-ArcGIS案例学习笔记

    计算平面面积和斜面面积-ArcGIS案例学习笔记 联系方式:谢老师,135_4855_4328,xiexiaokui#139.com 数据:实验数据\Chp8\Ex5\demTif.tif 平面面积= ...

  8. 计算概论(A)/基础编程练习2(8题)/3:计算三角形面积

    #include<stdio.h> #include<math.h> int main() { // 声明三角形的三个顶点坐标和面积 float x1, y1, x2, y2, ...

  9. 创建一个接口Shape,其中有抽象方法area,类Circle 、Rectangle实现area方法计算其面积并返回。又有Star实现Shape的area方法,其返回值是0,Star类另有一返回值boolean型方法isStar;在main方法里创建一个Vector,根据随机数的不同向其中加入Shape的不同子类对象(如是1,生成Circle对象;如是2,生成Rectangle对象;如是3,生成S

    题目补充: 创建一个接口Shape,其中有抽象方法area,类Circle .Rectangle实现area方法计算其面积并返回. 又有Star实现Shape的area方法,其返回值是0,Star类另 ...

随机推荐

  1. 【开发者portal在线开发插件系列三】字符串 及 可变长度字符串

    基础篇 基础场景见上面两个帖子,这里单独说明字符串和可变长度字符串的用法. 话不多说,开始今天的演(表)示(演) Profile和插件开发 添加一个string类型的属性: 在插件里添加一条数据上报消 ...

  2. Python之HTTP协议

    HTTP协议,又称超文本传输协议,主要用于浏览器与服务器之间的通信. HTTP 协议的制作者是蒂姆·伯纳斯-李,1991年设计出来的,HTTP 协议设计之前目的是传输网页数据的,现在允许传输任意类型的 ...

  3. 求局域网内所有在线主机的ip

    在一个局域网内,已知其中一台主机的ip为192.,子网掩码为255.,求所有其他在线主机的ip. shell 编码实现 #!/bin/bash netWorkIP=. ) do $netWorkIP$ ...

  4. es6 babel 安装以及使用

    1,安装好node(需要使用npm包管理工具) 2,在本地项目路径下npm init,格式化成功后会在项目下生成一个配置文件package.json 3,本地安装bable npm install - ...

  5. JsonModel的使用

    本人最近在开发一款医疗类的APP 发现接口返回的数据很复杂 手动解析的话对新手来说就是一场灾难 在分解成所需要的model类型时工作量非常的大,于是从网上查阅相关资料,发现JsonModel这个第三方 ...

  6. 数据挖掘算法(三)--logistic回归

    数据挖掘算法学习笔记汇总 数据挖掘算法(一)–K近邻算法 (KNN) 数据挖掘算法(二)–决策树 数据挖掘算法(三)–logistic回归 在介绍logistic回归之前先复习几个基础知识点,有助于后 ...

  7. 了解一下Mysql分布式事务及优缺点、使用案例(php+mysql)

    在开发中,为了降低单点压力,通常会根据业务情况进行分表分库,将表分布在不同的库中(库可能分布在不同的机器上),但是一个业务场景可能会同时处理两个表的操作.在这种场景下,事务的提交会变得相对复杂,因为多 ...

  8. PlantUML Viewer Chrome 插件 画时序图

    PlantUML通过简单直观的语言来定义示意图 使用 Chrome+ PlantUML Viewer的插件画图 1,打开chrome网上应用店 2,搜索plantuml viewer 并添加 3,扩展 ...

  9. numpy的基本API(三)——索引

    numpy的基本索引API iwehdio的博客园:https://www.cnblogs.com/iwehdio/ 1.单个元素的索引 对于一维数组,索引方式与内置的List相同.正索引从0开始,负 ...

  10. 一文讲透Dubbo负载均衡之最小活跃数算法

    本文是对于Dubbo负载均衡策略之一的最小活跃数算法的详细分析.文中所示源码,没有特别标注的地方均为2.6.0版本. 为什么没有用截止目前的最新的版本号2.7.4.1呢?因为2.6.0这个版本里面有两 ...