判断指定点在不在polydata内部的两种方式

1 篇文章 0 订阅
  1. 方法一:vtkCellLocator FindClosestPoint
    计算两点之间的距离,不过这个方法,返回的参数是距离的平方
bool CheckPointInsidePolyData(vtkPolyData * poly, double* pt)
{
	bool isInside=false; 
	auto cellLocator = vtkSmartPointer<vtkCellLocator>::New();
	cellLocator->SetDataSet(poly);
	cellLocator->BuildLocator();
	auto assistCell = vtkSmartPointer<vtkGenericCell>::New();
	// 接收最近点的坐标
	double closestPoint[3];
	// 接收最近点的距离
	double closestPointDist2;
	// 最近cell的id
	vtkIdType cellId; 
	int subId=-1;
	cellLocator->FindClosestPoint(point, closestPoint, assistCell, cellId, subId, closestPointDist2);
 
	if (-1!= subId&&closestPointDist2 < 0.001)
	{
		isInside=true;
	}
	return isInside;
}
  1. 方法二:vtkSelectEnclosedPoints
	double point[3] = { aPoint[0],aPoint[1],aPoint[2] };
	// 输入要检验的点集合,此处只设置一个点
	vtkNew<vtkPoints> checkPoint;
	checkPoint->InsertNextPoint(point[0], point[1], point[2]);
	// 将点转换成polydata对象
	vtkNew<vtkPolyData> pointsPolyData;
	pointsPolyData->SetPoints(checkPoint);

	vtkNew<vtkSelectEnclosedPoints> selectEnclosedPoints;
	selectEnclosedPoints->SetSurfaceData(myPolyData);
	selectEnclosedPoints->SetInputData(pointsPolyData);
	// 进行计算
	selectEnclosedPoints->Update();

	// 判断第一个点是否在polydata内部
	bool isInside = selectEnclosedPoints->IsInside(0);
```cpp

  • 8
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
要使用 VTK 计算两个牙齿的咬合,可以按照以下步骤进行: 1. 加载两个牙齿的 STL 模型文件,使用 vtkSTLReader 类读取模型数据,并分别存储为 vtkPolyData 对象。 ```cpp vtkSmartPointer<vtkSTLReader> reader1 = vtkSmartPointer<vtkSTLReader>::New(); reader1->SetFileName("tooth1.stl"); reader1->Update(); vtkSmartPointer<vtkPolyData> polyData1 = reader1->GetOutput(); vtkSmartPointer<vtkSTLReader> reader2 = vtkSmartPointer<vtkSTLReader>::New(); reader2->SetFileName("tooth2.stl"); reader2->Update(); vtkSmartPointer<vtkPolyData> polyData2 = reader2->GetOutput(); ``` 2. 对每个牙齿的 vtkPolyData 对象进行预处理,包括计算表面法向量、扫描线三角形插值等操作。 ```cpp vtkSmartPointer<vtkPolyDataNormals> normals1 = vtkSmartPointer<vtkPolyDataNormals>::New(); normals1->SetInputData(polyData1); normals1->ComputePointNormalsOn(); normals1->ComputeCellNormalsOff(); normals1->Update(); vtkSmartPointer<vtkPolyDataNormals> normals2 = vtkSmartPointer<vtkPolyDataNormals>::New(); normals2->SetInputData(polyData2); normals2->ComputePointNormalsOn(); normals2->ComputeCellNormalsOff(); normals2->Update(); vtkSmartPointer<vtkCleanPolyData> cleaner1 = vtkSmartPointer<vtkCleanPolyData>::New(); cleaner1->SetInputData(normals1->GetOutput()); cleaner1->Update(); vtkSmartPointer<vtkCleanPolyData> cleaner2 = vtkSmartPointer<vtkCleanPolyData>::New(); cleaner2->SetInputData(normals2->GetOutput()); cleaner2->Update(); ``` 3. 计算两个牙齿的交,使用 vtkIntersectionPolyDataFilter 类进行计算。 ```cpp vtkSmartPointer<vtkIntersectionPolyDataFilter> intersectionFilter = vtkSmartPointer<vtkIntersectionPolyDataFilter>::New(); intersectionFilter->SetInputData(0, cleaner1->GetOutput()); intersectionFilter->SetInputData(1, cleaner2->GetOutput()); intersectionFilter->Update(); vtkSmartPointer<vtkPolyData> intersectionPolyData = intersectionFilter->GetOutput(); ``` 4. 最后,将交vtkPolyData 对象保存为 STL 文件或其他格式的文件。 ```cpp vtkSmartPointer<vtkSTLWriter> stlWriter = vtkSmartPointer<vtkSTLWriter>::New(); stlWriter->SetFileName("intersection.stl"); stlWriter->SetInputData(intersectionPolyData); stlWriter->Update(); stlWriter->Write(); ``` 以上是使用 VTK 计算两个牙齿咬合的大致步骤,具体实现可能需要根据实际情况进行调整。
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值