vtkPolyDataConnectivityFilter提取模型的连通分区和遍历访问连通分量

864306337

已于 2023-09-11 11:46:09 修改

阅读量421

点赞数

分类专栏： VTK进阶算法文章标签： c++ 算法图论图形渲染 3d

于 2023-04-24 12:51:00 首次发布

本文链接：https://blog.csdn.net/qq_37366618/article/details/130340278

版权

VTK进阶算法专栏收录该内容

18 篇文章 12 订阅

订阅专栏

vtkPolyDataConnectivityFilter提取模型的连通分区和遍历访问连通分量

实现遍历访问模型中所有连通分区。

运行结果

在这里插入图片描述

方法介绍

重点是多次使用vtkPolyDataConnectivityFilter。先第一次使用vtkPolyDataConnectivityFilter提取所有连通分区出来，再在for循环使用vtkPolyDataConnectivityFilter提取指定的分区。

源代码

完整代码可以直接编译运行

//找出连通分区数量并遍历连通分区
#include <iostream>
#include <vtkNew.h>
#include <vtkSphereSource.h>
#include <vtkCleanPolyData.h>
#include <vtkAppendPolyData.h>
#include <vtkPolyDataConnectivityFilter.h>
using namespace std;

int main()
{
	//创建三个位置不同的球体
	vtkNew<vtkSphereSource> sphere1;
	sphere1->SetCenter(0, 0, 10);
	sphere1->SetRadius(2);

	vtkNew<vtkSphereSource> sphere2;
	sphere2->SetCenter(0, 10, 0);
	sphere2->SetRadius(2);

	vtkNew<vtkSphereSource> sphere3;
	sphere3->SetCenter(10, 0, 0);
	sphere3->SetRadius(2);
	//将三个球体合并为一个vtkPolyData数据
	//这一个vtkPolyData数据里就有三个球体，即三个连通分量
	vtkNew<vtkAppendPolyData> append;
	append->AddInputConnection(sphere1->GetOutputPort());
	append->AddInputConnection(sphere2->GetOutputPort());
	append->AddInputConnection(sphere3->GetOutputPort());
	//提取所有连通分量，获取连通分量数目
	vtkNew<vtkPolyDataConnectivityFilter> connectivity;
	connectivity->SetInputConnection(append->GetOutputPort());
	connectivity->SetExtractionModeToAllRegions();
	connectivity->Update();
	int num = connectivity->GetNumberOfExtractedRegions();
	cout << "连通分区的总顶点数量: " << connectivity->GetOutput()->GetNumberOfPoints() << endl;
	//遍历所有连通分区并访问数据
	for (int i = 0; i < num; i++)
	{
		//提取指定区域
		vtkNew<vtkPolyDataConnectivityFilter> region;
		region->SetInputConnection(connectivity->GetOutputPort());
		region->SetExtractionModeToSpecifiedRegions();
		region->AddSpecifiedRegion(i);
		//需要清除孤立点
		vtkNew<vtkCleanPolyData> clean;
		clean->SetInputConnection(region->GetOutputPort());
		clean->Update();
		cout << i + 1 << "分区的顶点数量: " << clean->GetOutput()->GetNumberOfPoints() << endl;
	}
	return 0;
}