解析OCC几何体，并在OSG中渲染显示

最新推荐文章于 2024-06-17 17:31:16 发布
最新推荐文章于 2024-06-17 17:31:16 发布
阅读量3.4k 收藏 23
点赞数 3
分类专栏： OpenCASCADE OpenSceneGraph
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_43712770/article/details/100972907
版权 
osg::Node* BuildShapeMesh(const TopoDS_Shape& aShape)
{
	{
		osg::ref_ptr<osg::Group> root = new osg::Group();
		osg::ref_ptr<osg::Geode> geode = new osg::Geode();

		//1几何体//
		osg::ref_ptr<osg::Geometry> triGeom = new osg::Geometry();
		osg::ref_ptr<osg::Vec3Array> vertices = new osg::Vec3Array();
		osg::ref_ptr<osg::Vec3Array> normals = new osg::Vec3Array();
		osg::ref_ptr<osg::Vec4Array> colors = new osg::Vec4Array;
		osg::ref_ptr<osg::UIntArray> arrayIndexs2 = new osg::UIntArray;//创建三角形顶点索引数据
		osg::ref_ptr<osg::DrawElementsUInt> arrayIndexs = new osg::DrawElementsUInt(osg::PrimitiveSet::TRIANGLES);

		//2法线///
		osg::ref_ptr<osg::Geometry> lineGeom = new osg::Geometry();
		osg::ref_ptr<osg::Vec3Array> vertices2 = new osg::Vec3Array();
		osg::ref_ptr<osg::Vec4Array> colors2 = new osg::Vec4Array;

		//3关键点//
		osg::ref_ptr<osg::Geometry> keyPointGeom = new osg::Geometry();
		osg::ref_ptr<osg::Vec3Array> vertices3 = new osg::Vec3Array();
		osg::ref_ptr<osg::Vec4Array> colors3 = new osg::Vec4Array;

		BRepMesh_IncrementalMesh(aShape, 1);
		//BRepMesh::Mesh(aShape, 3);//6.8.0后删除了
		cout << std::endl;
		gp_Pnt vertex1;
		gp_Pnt vertex2;
		gp_Pnt vertex3;

		Standard_Integer nVertexIndex1 = 0;
		Standard_Integer nVertexIndex2 = 0;
		Standard_Integer nVertexIndex3 = 0;

		TopExp_Explorer faceExplorer;
		for (faceExplorer.Init(aShape, TopAbs_FACE); faceExplorer.More(); faceExplorer.Next())
		{
			TopLoc_Location loc;
			TopoDS_Face aFace = TopoDS::Face(faceExplorer.Current());
			Handle_Poly_Triangulation triFace = BRep_Tool::Triangulation(aFace, loc);
			if (triFace.IsNull())
			{
				continue;
			}

			Standard_Boolean hasNormal = triFace->HasNormals();
			Standard_Boolean hasuvNormal = triFace->HasUVNodes();
			Standard_Integer l = triFace->Nodes().Length();
			Standard_Integer nTriangles = triFace->NbTriangles();

			TColgp_Array1OfPnt nodes(1, triFace->NbNodes());
			Poly_Array1OfTriangle triangles(1, triFace->NbTriangles());
			nodes = triFace->Nodes();
			triangles = triFace->Triangles();

			for (Standard_Integer i = 1; i <= nTriangles; i++)
			{
				vertex1.SetCoord(0.0, 0.0, 0.0);
				vertex2.SetCoord(0.0, 0.0, 0.0);
				vertex3.SetCoord(0.0, 0.0, 0.0);

				nVertexIndex1 = 0;
				nVertexIndex2 = 0;
				nVertexIndex3 = 0;

				Poly_Triangle aTriangle = triangles.Value(i);
				aTriangle.Get(nVertexIndex1, nVertexIndex2, nVertexIndex3);
				cout << "\tnVertexIndex1 = " << nVertexIndex1 << "\tnVertexIndex2 = " << nVertexIndex2 << "\tnVertexIndex3 = " << nVertexIndex3 << std::endl;

				//顶点
				//vertex1 = nodes.Value(nVertexIndex1).Transformed(loc.Transformation());
				//vertex1 = nodes.Value(nVertexIndex2).Transformed(loc.Transformation());
				//vertex1 = nodes.Value(nVertexIndex3).Transformed(loc.Transformation());
				vertex1 = nodes.Value(nVertexIndex1);
				//vertex1.SetCoord(vertex1.X() + 10.0, vertex1.Y() + 10.0, vertex1.Z() + 10.0);
				vertex2 = nodes.Value(nVertexIndex2);
				//vertex2.SetCoord(vertex2.X() + 10.0, vertex2.Y() + 10.0, vertex2.Z() + 10.0);
				vertex3 = nodes.Value(nVertexIndex3);
				//vertex3.SetCoord(vertex3.X() + 10.0, vertex3.Y() + 10.0, vertex3.Z() + 10.0);

				cout << "******" << "\tvertex1 = " << vertex1.X() << "," << vertex1.Y() << "," << vertex1.Z()
					<< "\tvertex2 = " << vertex2.X() << "," << vertex2.Y() << "," << vertex2.Z()
					<< "\tvertex3 = " << vertex3.X() << "," << vertex3.Y() << "," << vertex3.Z() << std::endl;
				cout << std::endl;
				//gp_Dir dir(0.0,1.0,0.0);
				//gp_Dir dir = triFace->Normal(i);

				//1几何体//
				{
					vertices->push_back(osg::Vec3(vertex1.X(), vertex1.Y(), vertex1.Z()));
					vertices->push_back(osg::Vec3(vertex2.X(), vertex2.Y(), vertex2.Z()));
					vertices->push_back(osg::Vec3(vertex3.X(), vertex3.Y(), vertex3.Z()));

					colors->push_back(osg::Vec4(1.0, 1.0, 0.0, 1.0));

					//arrayIndexs->push_back(nVertexIndex1);
					//arrayIndexs->push_back(nVertexIndex2);
					//arrayIndexs->push_back(nVertexIndex3);
				}

				///2法线//
				double ratio = 0.1;
				{
					gp_XYZ vector12(vertex2.XYZ() - vertex1.XYZ());
					gp_XYZ vector13(vertex3.XYZ() - vertex1.XYZ());
					gp_XYZ normal123 = vector13.Crossed(vector12);
					Standard_Real rModulus = normal123.Modulus();
					if (rModulus > gp::Resolution())
					{
						normal123.Normalize();
					}
					else {
						normal123.SetCoord(0., 0., 0.);
					}

					if (aFace.Orientation() != TopAbs_REVERSED)
					{
						normal123.Reverse();
					}


					//第1个顶点
					{
						gp_XYZ vector12(vertex2.XYZ() - vertex1.XYZ());
						gp_XYZ vector13(vertex3.XYZ() - vertex1.XYZ());
						gp_XYZ normal10 = vector13.Crossed(vector12);
						Standard_Real rModulus = normal10.Modulus();
						if (rModulus > gp::Resolution())
						{
							normal10.Normalize();
						}
						else {
							normal10.SetCoord(0., 0., 0.);
						}

						if (aFace.Orientation() != TopAbs_REVERSED)
						{
							normal10.Reverse();
						}

						//normals->push_back(osg::Vec3(dir.X(), dir.Y(), dir.Z()));
						normals->push_back(osg::Vec3(normal123.X(), normal123.Y(), normal123.Z()));

						vertices2->push_back(osg::Vec3(vertex1.X(), vertex1.Y(), vertex1.Z()));
						vertices2->push_back(osg::Vec3(vertex1.X(), vertex1.Y(), vertex1.Z()) + osg::Vec3(ratio * normal10.X(), ratio * normal10.Y(), ratio * normal10.Z()));
						colors2->push_back(osg::Vec4(1.0, 0.0, 0.0, 1.0));
						colors2->push_back(osg::Vec4(1.0, 0.0, 0.0, 1.0));
					}
					//第2个顶点
					{
						gp_XYZ vector21(vertex1.XYZ() - vertex2.XYZ());
						gp_XYZ vector23(vertex3.XYZ() - vertex2.XYZ());
						gp_XYZ normal20 = vector21.Crossed(vector23);
						Standard_Real rModulus = normal20.Modulus();
						if (rModulus > gp::Resolution())
						{
							normal20.Normalize();
						}
						else {
							normal20.SetCoord(0., 0., 0.);
						}

						if (aFace.Orientation() != TopAbs_REVERSED)
						{
							normal20.Reverse();
						}

						//normals->push_back(osg::Vec3(dir.X(), dir.Y(), dir.Z()));
						normals->push_back(osg::Vec3(normal123.X(), normal123.Y(), normal123.Z()));

						vertices2->push_back(osg::Vec3(vertex2.X(), vertex2.Y(), vertex2.Z()));
						vertices2->push_back(osg::Vec3(vertex2.X(), vertex2.Y(), vertex2.Z()) + osg::Vec3(ratio * normal20.X(), ratio * normal20.Y(), ratio * normal20.Z()));
						colors2->push_back(osg::Vec4(0.0, 1.0, 0.0, 1.0));
						colors2->push_back(osg::Vec4(0.0, 1.0, 0.0, 1.0));
					}
					//第3个顶点
					{
						gp_XYZ vector31(vertex1.XYZ() - vertex3.XYZ());
						gp_XYZ vector32(vertex2.XYZ() - vertex3.XYZ());
						gp_XYZ normal30 = vector32.Crossed(vector31);
						Standard_Real rModulus = normal30.Modulus();
						if (rModulus > gp::Resolution())
						{
							normal30.Normalize();
						}
						else {
							normal30.SetCoord(0., 0., 0.);
						}

						if (aFace.Orientation() != TopAbs_REVERSED)
						{
							normal30.Reverse();
						}

						//normals->push_back(osg::Vec3(dir.X(), dir.Y(), dir.Z()));
						normals->push_back(osg::Vec3(normal123.X(), normal123.Y(), normal123.Z()));

						vertices2->push_back(osg::Vec3(vertex3.X(), vertex3.Y(), vertex3.Z()));
						vertices2->push_back(osg::Vec3(vertex3.X(), vertex3.Y(), vertex3.Z()) + osg::Vec3(ratio * normal30.X(), ratio * normal30.Y(), ratio * normal30.Z()));
						colors2->push_back(osg::Vec4(0.0, 0.0, 1.0, 1.0));
						colors2->push_back(osg::Vec4(0.0, 0.0, 1.0, 1.0));
					}
				}

				3关键点//
				{
					vertices3->push_back(osg::Vec3(vertex1.X(), vertex1.Y(), vertex1.Z()));
					vertices3->push_back(osg::Vec3(vertex2.X(), vertex2.Y(), vertex2.Z()));
					vertices3->push_back(osg::Vec3(vertex3.X(), vertex3.Y(), vertex3.Z()));

					colors3->push_back(osg::Vec4(1.0, 1.0, 1.0, 1.0));
			}
		}
	}

		//1三角面片（几何体）//
		if (vertices->size() > 0)
		{
			//triGeom->setUseDisplayList(false);
			triGeom->setVertexArray(vertices.get());
			triGeom->addPrimitiveSet(new osg::DrawArrays(osg::PrimitiveSet::TRIANGLES, 0, vertices->size()));
			//triGeom->addPrimitiveSet(arrayIndexs);
			std::cout << vertices->size() << "\t" << arrayIndexs->size() << std::endl;
#if 1
			triGeom->setNormalArray(normals);
			triGeom->setNormalBinding(osg::Geometry::BIND_PER_VERTEX);
#else
			//自动生成法线，效果不好
			//NodePrimitiveSetNormalVisitor mv;
			//triGeom->accept(mv);
#endif
			triGeom->setColorArray(colors, osg::Array::BIND_OVERALL);

			//triGeom->getOrCreateStateSet()->setMode(GL_BLEND, osg::StateAttribute::ON);
			//triGeom->getOrCreateStateSet()->setRenderingHint(osg::StateSet::TRANSPARENT_BIN);
			//triGeom->getOrCreateStateSet()->setAttributeAndModes(new osg::CullFace(osg::CullFace::BACK), osg::StateAttribute::ON | osg::StateAttribute::PROTECTED);

			geode->addDrawable(triGeom);
		}


		//2法线//
		if (vertices2->size())
		{
			lineGeom->setVertexArray(vertices2.get());
			lineGeom->addPrimitiveSet(new osg::DrawArrays(osg::PrimitiveSet::LINES, 0, vertices2->size()));

			lineGeom->setColorArray(colors2, osg::Array::BIND_PER_VERTEX);
			lineGeom->getOrCreateStateSet()->setMode(GL_LIGHTING, osg::StateAttribute::OFF | osg::StateAttribute::PROTECTED);

			geode->addDrawable(lineGeom);
		}

		//3点//
		if (vertices3->size() > 0)
		{
			keyPointGeom->setVertexArray(vertices3.get());
			keyPointGeom->addPrimitiveSet(new osg::DrawArrays(osg::PrimitiveSet::POINTS, 0, vertices3->size()));
			keyPointGeom->setColorArray(colors3, osg::Array::BIND_OVERALL);
			osg::ref_ptr<osg::Point> pointSize = new osg::Point;
			pointSize->setSize(7.5);
			keyPointGeom->getOrCreateStateSet()->setAttribute(pointSize);
			keyPointGeom->getOrCreateStateSet()->setMode(GL_LIGHTING, osg::StateAttribute::OFF | osg::StateAttribute::PROTECTED);

			geode->addDrawable(keyPointGeom);

			//osg::ref_ptr<osgFX::Outline> outline = new osgFX::Outline;
			//outline->setWidth(8);
			//outline->setColor(osg::Vec4(1, 0, 0, 1));
			//outline->addChild(triGeom.get());
			//root->addChild(outline.get());

			//osg::PolygonMode* pm = new osg::PolygonMode(osg::PolygonMode::FRONT_AND_BACK, /*osg::PolygonMode::FILL | */osg::PolygonMode::LINE);
			//triGeom->getOrCreateStateSet()->setAttributeAndModes(pm, osg::StateAttribute::ON | osg::StateAttribute::PROTECTED);
			//triGeom->getOrCreateStateSet()->setAttributeAndModes(new osg::CullFace(osg::CullFace::BACK), osg::StateAttribute::ON | osg::StateAttribute::PROTECTED);
		}

		Use SmoothingVisitor to find the vertex average normals.
		//osgUtil::SmoothingVisitor sv;
		//sv.apply(*geode);

		root->addChild(geode);

		osg::MatrixTransform* mtAxes = new osg::MatrixTransform;
		mtAxes->setMatrix(osg::Matrix::scale(2.0, 2.0, 2.0));
		mtAxes->addChild(osgDB::readNodeFile("axes.osgt"));
		root->addChild(mtAxes);

		return root.release();
}
}

 

小胖七少爷
关注
  • 3
    点赞
  • 踩
  • 23
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 4
    评论
专栏目录
python OCC获取模型三角剖分的坐标点
qq_35342292的博客
04-07 2003
通过获取三角剖分的坐标点可以大概提取三维模型分布范围，通过浏览OCC C++的API网站，可以查找和阅读三角剖分的相关API。 OCC C++ API 我们可以在BRep_Tool这个包里可以找到有关三角坐标点的函数Triangulation，他的第一个参数是拓扑形状中的面类型，我们可以通过拓扑形状探索类TopologyExplorer来遍历一个shape的所有面。然后他的第二个参数是locati...
OCC把TopoDS_Edge转换为一系列点（二维）
qq_40330801的博客
05-27 904
OCC TopoDS_Edge细分
初探OSG+OpenCascade（简称：OCC）在QT上的实践
热门推荐
XXM的博客
05-19 1万+
因本人属于在OCC以及OSG领域是小白水平，因此文章可能存在错误的地方，恳请指正。 因为对几何建模比较有兴趣，而且受身边的人的影响，所以利用晚上和周末的时间写了一个简单的基于QT5+OCC7.0+OSG3.4的初级3D CAD(准确来说不算CAD,算是3D模型查看器)。 主要功能有： 读取3D文件，如图： 视角切换、按住滚轮拖动模型、放大缩小、视觉样式显示、导出文件功能。其他操作在接下来...
基于OCC+OSG的读取IGS模型显示其装配以及模型颜色
最新发布
XXM的博客
06-17 364
OSG+OCC读取igs模型，显示模型结构和颜色
初探基于OSG+OCC的CAD之Netgen体网格划分与显示
XXM的博客
07-28 1863
利用netgen网格划分程序的nglib接口划分体网格
OCC笔记：访问拓扑边TopoDS_Edge的起末点
憧憬，思考，奋斗，飞翔
04-12 667
使用TopExp取出拓扑边TopoDS_Edge的起末拓扑顶点TopoDS_Vertex； 使用BRep_Tool从拓扑顶点取出几何点gp_Pnt
集成OCC+Gmsh+OSG的CAD测试程序
12-24
在集成OCC和Gmsh的数据后，OSG可以将这些几何模型和网格数据转化为可视化的3D场景，让用户能够直观地查看和交互模型。 在"TestProject"这个压缩包中，包含了集成这些组件的测试程序。这个程序可能实现了读取Gmsh的...
用OCC结合QT创建并显示一个几何.docx
03-13
在本篇文章中，我们将介绍如何使用OpenCASCADE（OCC）结合Qt创建并显示一个几何体。OpenCASCADE是一个功能强大的三维几何模型库，提供了广泛的几何计算和处理功能。Qt是一个功能强大的跨平台应用程序开发框架，提供...
opencascade OCC中文帮助文档
01-11
4. **图形渲染**：了解如何使用OCC的图形渲染模块将几何模型显示在屏幕上，以及自定义视图和光照设置。 5. **数据交换**：学习OCC提供的数据导入和导出接口，如IGES、STEP等标准格式，以与其他CAD系统进行交互。 6...
关于OCC中操作和计算总结.doc
08-07
在Python OCC库中，OCC（Open CASCADE Technology）是一个强大的3D建模和几何处理框架，主要用于CAD、CAE和CAM应用。这个文档主要涵盖了OCC中的一些关键操作和计算，包括STP文件的导入导出、AIS（Application ...
基本体理论之OCC体系结构1
08-03
Open CASCADE（以下简称 OCC）是一种基于边界表示（B-Rep）的三维形体模型，在 OCC 中，系统包括拓扑（Topologic）和几何（Geometric）两种主要结构。 OCC 模型的几何拓扑结构 在 OCC 中，三维形体采用了基于 B-...
Opencascade 判断TopoDS_Edge
HNU_slacker的博客
05-14 264
TopoDS_Face类型通过BRep_Tool::Surface判断。其他继承Geom_Curve的子类也可判断。
OCCT中的裁剪平面（Geom_RectangularTrimmedSurface）
万里归来少年心
09-04 1305
    裁剪平面（Geom_RectangularTrimmedSurface）是基础平面的一部分，它在基础平面的u方向上用两个值限制范围，在基础平面的v方向上也用两个值限制范围。裁剪平面的区域一定位于基础平面的区域内，它的定义如下：    基础面    u方向和v方向的参数值(umin, umax) 和(vmin, vmax)。裁剪平面是在基础面的一份拷贝上做裁剪，因此，当基础面改变时，裁剪...
第二章：Pythonocc官方demo 案例59（通过继承TopoDS_Edge类构造新的类）
loujiand的博客
11-27 181
pythonocc 通过继承TopoDS_Edge类构造新的类。
关于OpenCasCade 曲线均分的测试，显示效果用OSG
XXM的博客
05-20 1997
抽空写了个曲线均分测试，把任意曲线平成分成长度相等的N段。仅此记录, 下面测试了三种，前2种是测试方案： 第一种是通过参数进行分割:分别获得曲线的 FirstParameter 和 LastParameter ，然后对参数进行分割，获得n个ui，并对每个ui调用D0（获得这个点的坐标值）或D1（获得这个点的坐标值和切向量）。这个方法的优点是，简单易行，好操作，缺点均分参数轴获得的曲线的分割并不是均匀的，在OSG中显示如下： 测试代码如下： /////////////////////////...
基于OCC与OSG的CAD平台构建值gmsh划分效果测试
XXM的博客
03-19 823
大半年没更新博客了，目前处在项目的重构中，一时没有增加特别的功能模块，所以本文主要是对gmsh网格划分模块的一个效果展示，更多功能将持续改进中...........四、gmsh的parallelograms效果。后记：后续将对网格模块作进一步优化。一、gmsh的波前法划分效果。二、gmsh的三角法划分效果。三、gmsh的网格合并效果。
初探OSG+OpenCascade（简称：OCC）在QT上的实践之尝试集成Netgen6.2网格划分程序
XXM的博客
10-29 3792
Netgen是一款优秀的开源四面体网格划分程序，如果是准备做CAE方向，网格时其中的一个重点，本文就Netgen6.1集成做个简要流程，流程也是比较简单，具体如下： 第一.编译netgen with OCC,这个网上资料很多，因为我们是准备集成netgen的库(nglib),所以我在cmake时候，在UESER选项中只勾选了USC_OCC,截图如下： 基本的cmake配置是：camke3.12版本，VS2015release 64，选择OCC的时候我直接把OCC7.4的cmake对应就行了。截图如.
osg加载occ模型并显示
05-13
要加载并显示 OCC 模型，您需要使用 OpenSceneGraph（OSG）和 OpenCASCADE（OCC），并将它们集成在一起。以下是加载 OCC 模型并在 OSG 中显示的基本步骤： 1. 安装 OSG 和 OCC 并将其设置为您的项目依赖项。 2. 使用 OCC API 加载模型并将其转换为 OSG 几何体。 ```cpp #include <OpenCASCADE/TopoDS.hxx> #include <OpenCASCADE/BRepTools.hxx> #include <OpenCASCADE/IGESControl_Reader.hxx> #include <osg/Geometry> #include <osg/Geode> osg::ref_ptr<osg::Node> loadOCCModel(const std::string& fileName) { // Step 1: Load the model using OCC API TopoDS_Shape shape; IGESControl_Reader reader; reader.ReadFile(fileName.c_str()); reader.TransferRoots(); shape = reader.OneShape(); // Step 2: Convert the OCC geometry to OSG geometry osg::ref_ptr<osg::Geode> geode = new osg::Geode; BRepTools::Clean(shape); BRepTools::Update(shape); TopExp_Explorer exp(shape, TopAbs_FACE); while (exp.More()) { TopoDS_Face face = TopoDS::Face(exp.Current()); Handle(Geom_Surface) surface = BRep_Tool::Surface(face); if (!surface.IsNull()) { // Create OSG geometry from OCC surface osg::ref_ptr<osg::Geometry> geom = new osg::Geometry; // Add vertices to the geometry // ... // Add primitives to the geometry // ... geode->addDrawable(geom); } exp.Next(); } return geode; } ``` 3. 将 OSG 几何体添加到场景图中并显示。 ```cpp #include <osgViewer/Viewer> int main(int argc, char** argv) { // Create the viewer osgViewer::Viewer viewer; // Load the model and add it to the scene osg::ref_ptr<osg::Node> model = loadOCCModel("model.igs"); viewer.setSceneData(model); // Start the viewer return viewer.run(); } ``` 这只是一个基本示例，您可能需要根据您的要求进行适当的修改。
评论 4
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

小胖七少爷

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值