【FastCAE源码阅读3】几何模型显示：从OCC对象到VTK对象

已于 2023-11-06 15:48:57 修改
阅读量945 收藏 12
点赞数 7
分类专栏： FastCAE源码阅读 文章标签： 算法 c++ 图形渲染
于 2023-11-03 14:59:38 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/loveoobaby/article/details/134203489
版权

从几何到显示还是比较麻烦的，需要将几何对象转换成渲染对象，涉及几何建模、面的三角化、图形渲染等学科，阅读本文需了解一些基本的OCC、VTK编程

一、几何体显示基本流程

FastCAE几何内核使用的是OCC，显示渲染用的VTK，那么就存在将OCC建模后的几何对象变成VTK支持类型的数据这个过程。

FastCAE在用OCC创建完模型之后，会生成一个Geometry::GeometrySet对象，这个对象包含了最原始的OCC中TopoDS_Shape实例对象。创建完成后都会发出信号emit showSet(set);，这个信号会触发将TopoDS_Shape变成VTK显示对象。这个段代码是在void GeometryViewProvider::showGeoSet(Geometry::GeometrySet *set, bool render)函数中实现。

二、几何对象变成显示对象过程

几何体是由点、线、面构成的，进行显示的时候是分别提取这些数据，进行单独显示的。也就是说我们看到的一个立方体等几何对象是由三个VTK显示对象Actor拼起来的。注意一个体是没有直接的数据对象与其对应的，体只是逻辑上的概念，所以很多软件中隐藏一个体的面会发现内部是空的。这个函数代码如下：

void GeometryViewProvider::showGeoSet(Geometry::GeometrySet *set, bool render)
{
	QList<vtkPolyData *> viewPolys = _viewData->transferToPoly(set); // 对几何体进行点、线、面的拆分
	vtkPolyData *facePoly = viewPolys.at(0); // 面的多边形数据
	vtkPolyData *edgePoly = viewPolys.at(1); // 边的数据集
	vtkPolyData *pointPoly = viewPolys.at(2); // 点数据集
	GeoViewObj viewObj;
	if (facePoly != nullptr) // 创建显示面的Actor
	{
		vtkSmartPointer<vtkActor> actor = vtkSmartPointer<vtkActor>::New();
		vtkSmartPointer<vtkPolyDataMapper> mapper = vtkSmartPointer<vtkPolyDataMapper>::New();
		mapper->SetInputData(facePoly);
		actor->SetMapper(mapper);
		bool vis = set->isVisible();
		bool show = Setting::BusAPI::instance()->getGraphOption()->isShowGeoSurface();
		actor->SetVisibility(show && vis); // 设置是否显示面
		actor->SetPickable(false);   // 不可鼠标拾取？哪里打开呢？
		actor->GetProperty()->SetRepresentationToSurface(); // 这个函数不调用也没看到啥影响
		_preWindow->AppendActor(actor, ModuleBase::D3, false); // 将面actor添加到场景中
		viewObj._faceObj = QPair<vtkActor *, vtkPolyData *>(actor, facePoly);
	}

	// 创建显示边的Actor
	if (edgePoly != nullptr) 
	{
		vtkSmartPointer<vtkActor> actor = vtkSmartPointer<vtkActor>::New();
		vtkSmartPointer<vtkPolyDataMapper> mapper = vtkSmartPointer<vtkPolyDataMapper>::New();
		mapper->SetInputData(edgePoly);
		actor->SetMapper(mapper);
		bool vis = set->isVisible();
		bool show = Setting::BusAPI::instance()->getGraphOption()->isShowGeoEdge();
		actor->SetVisibility(show && vis);
		actor->SetPickable(false);
		actor->GetProperty()->SetRepresentationToWireframe();
		//			actor->GetProperty()->EdgeVisibilityOn();
		float width = Setting::BusAPI::instance()->getGraphOption()->getGeoCurveWidth();
		actor->GetProperty()->SetLineWidth(width);
		_preWindow->AppendActor(actor, ModuleBase::D3, false);
		viewObj._edgeObj = QPair<vtkActor *, vtkPolyData *>(actor, edgePoly);
	}

	// 创建显示点的Actor
	if (pointPoly != nullptr)
	{
		vtkSmartPointer<vtkActor> actor = vtkSmartPointer<vtkActor>::New();
		vtkSmartPointer<vtkPolyDataMapper> mapper = vtkSmartPointer<vtkPolyDataMapper>::New();
		float size = Setting::BusAPI::instance()->getGraphOption()->getGeoPointSize();
		mapper->SetInputData(pointPoly);
		actor->SetMapper(mapper);
		bool vis = set->isVisible();
		bool show = Setting::BusAPI::instance()->getGraphOption()->isShowGeoPoint();
		actor->SetVisibility(show && vis);
		actor->SetPickable(false);
		actor->GetProperty()->SetRepresentationToPoints();
		actor->GetProperty()->SetPointSize(size);
		_preWindow->AppendActor(actor, ModuleBase::D3, false);
		viewObj._pointObj = QPair<vtkActor *, vtkPolyData *>(actor, pointPoly);
	}
	_geoViewHash.insert(set, viewObj);
	if (render)
		_preWindow->resetCamera();
}

三、几何面数据的提取

几何体的面、边、点转换为vtkPolyData对象的方式类似，都在函数GeometryViewData::transferToPoly()中，这里只分析面数据的提取。

/* 提取几何表面的多边形数据集 */
vtkPolyData *GeometryViewData::transferFace(Geometry::GeometrySet *gset)
{
	TopoDS_Shape *shape = gset->getShape(); // 获取OCC的Shape
	TopExp_Explorer faceExp(*shape, TopAbs_FACE); // 这个类可访问Shape的拓扑关系
	QList<Handle(TopoDS_TShape)> tshapelist; // 放置已访问的Face对象，防止重复访问
	vtkSmartPointer<vtkAppendPolyData> appendFilter = vtkSmartPointer<vtkAppendPolyData>::New();
	vtkPolyData *polyData = vtkPolyData::New();
	int beg = 0;
	for (int index = 0; faceExp.More(); faceExp.Next(), ++index) // 开始遍历面
	{
		const TopoDS_Shape &s = faceExp.Current();
		Handle(TopoDS_TShape) ts = s.TShape();
		if (tshapelist.contains(ts)) // 防止重复访问，会存在这种情况吗？
			continue;
		tshapelist.append(ts);
		IVtkOCC_Shape::Handle aShapeImpl = new IVtkOCC_Shape(s); // OCC提供IVtkOCC_Shape类
		vtkSmartPointer<IVtkTools_ShapeDataSource> DS = vtkSmartPointer<IVtkTools_ShapeDataSource>::New(); // OCC提供的VTK数据源
		DS->SetShape(aShapeImpl); 
		vtkSmartPointer<vtkCleanPolyData> cleanFilter = vtkSmartPointer<vtkCleanPolyData>::New();
		cleanFilter->SetInputConnection(DS->GetOutputPort());
		cleanFilter->Update();
		vtkSmartPointer<vtkPolyData> tpolys = vtkSmartPointer<vtkPolyData>::New();
		vtkPolyData *tpolydata = cleanFilter->GetOutput();
		const int np = tpolydata->GetNumberOfPoints(); // 点的数量
		const int nc = tpolydata->GetNumberOfCells();  // cell的数量
		vtkPoints *points = vtkPoints::New();
		for (int i = 0; i < np; i++) //  提取几何点数据
		{
			double *coor = tpolydata->GetPoint(i);
			points->InsertNextPoint(coor);
		}
		tpolys->SetPoints(points); // 设置几何点数据
		vtkCellArray *cells = vtkCellArray::New();
		for (int i = 0; i < nc; ++i)
		{
			vtkCell *cell = tpolydata->GetCell(i);
			vtkIdList *ceid = cell->GetPointIds();
			if (ceid->GetNumberOfIds() == 3) // 只提取三角形，这里获取的cell包含点、线、三角形，一个立方体最后应该有12个三角形
			{
				vtkTriangle *triangle = vtkTriangle::New();
				triangle->DeepCopy(cell); // 有DeepCopy接口
				cells->InsertNextCell(triangle);
			}
		}
		tpolys->SetPolys(cells); // 设置拓扑多边形数据集的cell

		// 法线数据？打开与否貌似不影响显示效果
		vtkSmartPointer<vtkPolyDataNormals> normals = vtkSmartPointer<vtkPolyDataNormals>::New();
		normals->SetInputData(tpolys);
		normals->FlipNormalsOn();
		normals->Update();
		vtkPolyData *facePoly = normals->GetOutput();
		const int ncell = facePoly->GetNumberOfCells();
		if (ncell < 1)
			continue;
		GeometryViewObject *obj = new GeometryViewObject(GeometryViewObject::Face, beg, beg + ncell - 1, ts);
		beg += ncell;
		appendFilter->AddInputData(facePoly);  // 追加到appendFilter中
		auto setViewObj = this->getGeosetObj(gset);
		setViewObj->appendFaceViewObj(index, ts, obj);
	}
	appendFilter->Update();
	polyData->DeepCopy(appendFilter->GetOutput());
	auto setViewObj = this->getGeosetObj(gset);
	setViewObj->setFacePoly(polyData); // 保存面多边形数据
	const int npc = polyData->GetNumberOfCells();
	if (npc < 1)
		return nullptr;
	return polyData;
}

这里要注意几个类：

  1. TopExp_Explorer: OCC提供的遍历几何体拓扑结构的类；
  2. IVtkOCC_Shape、IVtkTools_ShapeDataSource：这俩类可以将面三角化，并变成VTK的cell，注意cell包含面上的点、边、三角形，所以提取面的时候只提取了三角形，这俩类同样是OCC提供的；
  3. vtkPolyDataNormals： 看类名猜测是生成法线的，但是不用这类显示也没啥问题；
  4. vtkAppendPolyData：可以将各个面的数据追加在一起。

四、面数据的最终数据结构

OCC对象在转换过程中，会生成成一系列对象，用于后续的拾取高亮等功能。FastCAE中有个组件叫GeometryViewProvider，其中的字段_viewData（类型 GeometryViewData）保存这些数据。GeometryViewData类中有哈希表，key是几何对象，value是GeoSetViewObject。GeoSetViewObject又有hash表_faceObjHash、_faceViewObjs。以立方体为例，_faceViewObjs保存6个面对应的GeometryViewObject对象，每个GeometryViewObject对象记录这这个面对应的cell范围及所有面的多边形数据集vtkPloyData。这个数据集在GeoSetViewObject也有引用（_facePoly字段），也是面Actor中的数据集，一个体的所有面组成一个Actor。最终形成如下的一种数据结构：
在这里插入图片描述

边和点的数据与面类似，这里就不再分析了。

最后来张图吧：
在这里插入图片描述

loveoobaby
关注
  • 7
    点赞
  • 踩
  • 12
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
vtk解析转换OCC中的数据
天冷请你吃火锅
03-01 799
vtk转化OCC数据
pythonocc-core：用于3D CADBIMPLMCAM的Python软件包
01-28
pythonocc-core 关于 pythonocc提供3D建模和数据交换功能。 它旨在与CAD / PDM / PLM和BIM相关的开发。 最新版本： 特征 pythonocc提供以下功能： 从Python到几乎所有数千个OpenCascade C ++类的完整访问。 类和方法/函数共享相同的名称，并尽可能使用相同的签名 最著名的Python Gui（pyQt，PySide1和2，wxPython）的3D可视化 使用WebGL和/或x3dom渲染器的Web浏览器中的3D可视化 在Jupyter笔记本电脑中进行3D可视化和工作 用于数据交换，拓扑运算，惯性计算等的各种实用Python类/方法。 在mybinder在线尝试 请点击 打开运行最新pythonocc-core 7.4.1，gmesh 4.5.3（ ）和最新IfcOpenshell（ ）的jupyter笔记本开发分支。 下载/安装Linux / OSX / Windows二进制文件 pythonocc为python 3.6、3.7和3.8提供了预编译的conda（它们依赖于dlr-sc和conda-forge co
OCCT 建模，VTK显示
pengzeyu0524的博客
06-21 3448
2019/6/21 在OCCT中visualization中有现成的VTK接口 从Topo::DS_Shape到Mapper的接口如下所示： IVtkOCC_Shape::Handle aShapeImpl = new IVtkOCC_Shape(MakeBottle(20,30,8)); vtkSmartPointer<IVtkTools_ShapeDataSource>...
FastCAE:FastCAE，是一套开源国产CAE软件集成开发平台。面向求解器开发者，平台定义了规范的数据接口，支持插件开发模式，可视化配置手段，基于平台可轻松集成用户自研求解程序和第三方求解器。 FastCAE的本质是一套CAE软件集成接口标准和数据标准框架，基于该框架集成任意的求解器、前处理、后处理等功能模块。平台默认提供完整的几何交互、网格划分、结果可视化、边界条件、材料库等模块，为自主CAE软件产品化提供一体化解决方案，让用户享受集合前处理、求解计算、后处理于一体的平台化集成体验
05-24
FastCAE源码下载与编译说明（源码下载编译必读） 特别注意： 源码绝对路径中不要出现中文字符; 源码文件夹中不要出现空格(){}*/?|\等特殊字符 要求Windows7以上的64位系统 要求自行安装Visual Studio 2017和Qt5.14.2 源码下载： 两种方式可以下载源码，直接下载和使用git克隆 直接下载，从GitHub直接一次性打包下载全部源码。 git克隆，使用git克隆代码仓库，后续可以不断拉取更新。 #文件夹结构 Code FastCAE源码 OCC OpenCASCADE7.4.0编译的第三方库，适用于windows平台 VTK VTK9.0.1编译的第三方库，适用于windows平台 output 程序的编译输出 cgns cgns第三方库，适用于windows平台 quazip 压缩文件的第三方库 编译过程 源码编译支
pythonocc-demos：pythonocc CAD软件包的一组示例和演示
02-04
一些与pythonocc相关的代码段，示例，jupter笔记本等。 要求 ： pythonocc-core 7.4.1 conda install -c conda-forge pythonocc-core=7.4.1 jupyter，如果要测试jupyter笔记本以及pythreejs。 仓库结构： 资源：各种格式的2D图像，3D模块。 由python脚本使用 示例：每个都描述pythonocc功能的小型python脚本 $ cd examples $ python core_helloworld.py jupyter_notebook：在jupyter笔记本中
用OCC结合QT创建并显示一个几何.docx
03-13
OCC结合Qt创建并显示一个几何体
FastCAE 2.0版本产品介绍
FastCAE的博客
12-01 730
FastCAE，是一套开源国产CAE软件集成开发平台，采用BSD3开源协议，用户可免费使用该框架并鼓励衍生软件商业化。 FastCAE本质是一套CAE软件集成接口标准和数据标准框架，基于该框架集成任意的求解器、前处理和后处理模块。平台面向求解器开发者提供CAE软件界面与通用功能模块，快速研发集成框架，提高自研求解程序产品化效率。平台提供了完善的前后处理功能与丰富的可视交互界面，支持API级二次开发和插件模式功能扩展。
【OCC学习3】OCC与VTK整合
loveoobaby的博客
07-02 2057
3. Visual Studio配置 在Visual Studio中配置OCC的头文件、VTK的头文件； 配置附加库目录，包括OCC、VTK的库目录； 配置附加依赖项，如果不知道怎么根据链接错误找lib库，可使用Far Manager查找； 4. 配置调试PATH环境变量。PATH变量中需增加OCC、VTK、TBB、Openvr、freeimage、ffmpeg库的dll目录。官网编译好的包包含这些库的，可直接下载。......
【FastCAE源码阅读5】使用VTK实现鼠标拾取对象并高亮
loveoobaby的博客
11-06 1124
其中当拾取模式是点的时候，使用vtkPointPicker，如果拾取成功返回点的id及对应的Actor（一个Actor可能包含多个点）。注意当拾取线的时候，需要提高拾取Tolerance，不然线很难拾取到。另外要注意的是，当拾取线的时候，面Actor是拾取不到的，同样当拾取面的时候，线Actor也是拾取不到。拾取对象是在PropPickerInteractionStyle类实现的，该类是vtkInteractorStyleRubberBandPick的子类，重写原来的鼠标处理函数。拾取对象之后，会发出。
FastCAE版本更新——以“工匠”精神对待每一个版本
FastCAE的博客
12-02 179
FastCAE开发人员不断对产品进行完善，以回馈用户更好的使用体验，在此期间，FastCAE在原有基础上进行了更新优化，更新后的安装包以及相应源码现已发布，欢迎大家在 FastCAE官网下载使用。
【FastCAE源码阅读2】FastCAE如何调用OpenCasCade
loveoobaby的博客
08-28 921
FastCAE使用的是OCCT（OpenCascade）几何内核，FastCAE几何模块菜单上每个功能对应着一个创建几何的类，本文简单分析这部分源码。打开FastCAE几何模块，在功能菜单中可以看到很多几何操作功能，包括创建立方体、创建圆柱，创建点，布尔操作，物体进行变换等。这每一个功能被封装成一个类，底层会调用OCCT的API，代码在GeometryCommand模块中。其中，GeoCommandBase是所有几何创建命令的基类，也是QObject对象，定义了通用接口，基础的QT槽函数，通用字段。...
node-occ：使用OpenCascade和NodeJS构建BREP实体-3D建模
01-28
节点occ 用于实体建模的OpenCascade nodejs扩展。 这个nodejs扩展为nodejs提供了坚实的基础。 它提供了一个简单而功能强大的javascript api来构建3D几何模型。 该项目带有一组围绕OpenCascade API的V8包装器和一个示例Web应用程序。 简单的例子 var occ = require ( "node-occ" ) . occ ; // construct a box var box = occ . makeBox ( [ 0 , 0 , 0 ] , [ 100 , 100 , 50 ] ) ; // construct a cylinder var cyl = occ . makeCylinder ( [ 50 , 50 , - 10 ] , [ 50 , 50 , 60 ] , 40 ) ; // cut the box with cylinder box = occ . cut ( box , cyl ) ; // save result to a STEP file occ . writeSTEP ( "somefile.
JellyCAD:基于OCC的简单CAD软件
05-11
JellyCAD 一个简单的CAD软件。 发展 在Windows 7平台上开发基于：Qt 5.9.0 MSVC 14 OpenCASCADE-7.3.0 建造 1.从以下下载最新版本的OpenCASCADE： //www.opencascade.com/content/latest-release 2.使用二进制文件安装OpenCASCADE，或从源代码构建OCC。 3.使用Qt打开JellyCAD.pro。 4.配置OpenCASCADE环境。 影响 反馈 以色列（ ）
【FastCAE源码阅读8】调用gmsh生成网格
最新发布
loveoobaby的博客
11-13 614
FastCAE使用gmsh进行网格划分，划分的时候直接启动一个新的gmsh进程，个人猜测这么设计是为了规避gmsh的GPL协议风险。
FastCAE集成FEMAP软件的网格剖分功能
FastCAE的博客
12-02 609
Siemens FEMAP是一款独立于CAD且基于Windows技术的前/后处理高级工程有限元分析（FEA）应用程序，它可以帮助工程师和分析师利用FEA建模解决方案以合理的成本轻松、准确地完成复杂的仿真任务。 FastCAE提供了统一的交互和显示功能，通过接口可集成自研、开源和商业的几何、网格、后处理内核模块。本案例采用OEM方式集成西门子FEMAP网格划分组件，支持复杂工程级的网格划分功能，在功能强大的同时，界面上可实现自定义、个性化的风格。
FastCAE工业仿真软件开源学院
FastCAE的博客
12-01 528
FastCAE工业仿真软件开源学院，简称FastCAE开源学院，由FastCAE开源团队牵头组建的非营利性组织。面向社会提供工业仿真软件研发及应用相关的技术培训、资源发布、产业对接等服务。
验证occ和vtk整合工作的demo
weixin_30784141的博客
08-14 487
在编译occ通过过后，我需要验证occ是否能够正常结合vtk进行开发工作 使用CMake进行环境变量设置： CMakeList.txt PROJECT (IGESReader) #VTK Part: FIND_PACKAGE(VTK) IF (VTK_FOUND) INCLUDE(${VTK_USE_FILE}) ELSE(VTK_FOUND) MESSAGE(FA...
基于FastCAE的电磁仿真求解器集成中个性化功能与技术实现方法分享
FastCAE的博客
02-06 313
基于FastCAE源码的个性化定制开发是一种灵活度较高的定制开发方式，能够实现复杂的个性化功能，也能够对FastCAE的基础功能进行修改。本文以中算电磁仿真软件作为项目介绍案例，主要涉及到前后处理显示以及GMSH网格划分的相关个性化优化操作。
Windows occ(OpenCascade)编译 (occ+vtk+qt+cmake示例)
abc的博客
08-21 4573
OCC在windows平台使用vs2019+qt5.15+vtk9.1进行编译。并创建cmake项目及测试。
使用gmsh::model::occ需要引入的头文件
05-28
使用 `gmsh::model::occ` 需要引入以下头文件： ```cpp #include <gmsh.h> #include <gmsh/GModel.h> #include <gmsh/GModelOCC.h> ``` 其中，`<gmsh.h>` 是 Gmsh 的主头文件，包含了大部分的 Gmsh 类和函数的声明；`<gmsh/GModel.h>` 是 Gmsh 中的模型类头文件，包含了 `GModel` 类的声明；`<gmsh/GModelOCC.h>` 是 Gmsh 中的 OCC 模型类头文件，包含了 `GModelOCC` 类的声明。 在使用 `gmsh::model::occ` 时，还需要在代码中添加以下命名空间： ```cpp using namespace gmsh; using namespace gmsh::model; using namespace gmsh::model::occ; ``` 这样才能够使用 `gmsh::model::occ` 中定义的类和函数。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值