空间八叉树---vtkCellLocator

关于空间点，VTK有一堆我还没弄明白的东西：vtkKdTree , vtkOBBTree ,vtkOctreePointLocator , vtkModifiedBSPTree

自己写了半天，发现VTK一个函数搞定的事情，悲催
vtkSmartPointer<vtkPolyData> pd = point0;
//--------------------------------------------
// Create the tree
vtkSmartPointer<vtkCellLocator> cellLocator =
vtkSmartPointer<vtkCellLocator>::New();
cellLocator->SetDataSet(pd);
cellLocator->BuildLocator();
double testPoint[3] = {146,128,111};
cellLocator->FindClosestPoint(testPoint, closestPoint, cellId, subId, closestPointDist2);
std::cout << "Coordinates of closest point: " << closestPoint[0] << " " << closestPoint[1] << " " << closestPoint[2] << std::endl;
std::cout << "Squared distance to closest point: " << closestPointDist2 << std::endl;
std::cout << "CellId: " << cellId << std::endl;

double point_temp[3] = {146,128,111};
double squaredDistance = 1000;
double first_point[3];
for (vtkIdType i = 0; i< pPos->GetNumberOfPoints()-1 ; i++)
{
double point_temp2[3];
point_temp2[0] = pPos->GetPoint(i)[0];
point_temp2[1] = pPos->GetPoint(i)[1];
point_temp2[2] = pPos->GetPoint(i)[2];
double squaredDistance2 = vtkMath::Distance2BetweenPoints(point_temp, point_temp2);
if (squaredDistance > squaredDistance2)
{
squaredDistance = squaredDistance2;
first_point[0] = point_temp2[0];
first_point[1] = point_temp2[1];
first_point[2] = point_temp2[2];
}
}
std::cout << "squaredDistance:"<< squaredDistance<< std::endl;
std::cout << "first_point:"<< first_point[0]<<","<< first_point[1]<<","<< first_point[2]<< std::endl;

继续~KD-Tree：搜索最近的n个点
//Create the tree
vtkSmartPointer<vtkKdTreePointLocator> pointTree =
vtkSmartPointer<vtkKdTreePointLocator>::New();
pointTree->SetDataSet(pd);
pointTree->BuildLocator();
// Find the k closest points to (0,0,0)
unsigned int k = 50;
//double testPoint[3] = {0.0, 0.0, 0.0};
vtkSmartPointer<vtkIdList> result =
vtkSmartPointer<vtkIdList>::New();
pointTree->FindClosestNPoints(k, testPoint, result);
for(vtkIdType i = 0; i < k; i++)
{
vtkIdType point_ind = result->GetId(i);
double p[3];
pd->GetPoint(point_ind, p);
std::cout << "Closest point " << i << ": Point "
<< point_ind << ": (" << p[0] << ", " << p[1] << ", " << p[2] << ")" << std::endl;
}

吧主 你好 请问可以提供一下您上面所看的八叉树细分方面的资料/程序吗？最近一直在做，比较纠结 自己是在搞不懂了

楼主 可以再仔细讲讲这个八叉树的？？

这些类都是空间查询的，各个类实现方法略有不同，用过kdtree，该算法是先构建k维空间点集，然后可以快速查询到某一个位置附近的一个或者几个点