2012-07-09 160 views
2

我正在使用一个项目的VTK,似乎无法弄清楚它的一部分。我试图遍历几千个点,并找到每个点最接近的5个点。看起来像一个简单的循环操作,但我的问题是,由于某种原因,我被告知相同的5分是我的数据中每个点的最近点...我知道这是不正确的。我将附加代码比特我描述下图:如何在找到每个点的最近邻点时迭代VTK中的点?

//test 
    vtkSmartPointer<vtkPointSource> pointSource = 
    vtkSmartPointer<vtkPointSource>::New(); 
    pointSource->SetNumberOfPoints(Output->GetNumberOfPoints()); 
    pointSource->Update(); 

    vtkSmartPointer<vtkKdTreePointLocator> Tree = 
    vtkSmartPointer<vtkKdTreePointLocator>::New(); 
    Tree->SetDataSet(pointSource->GetOutput()); 
    Tree->BuildLocator(); 

    unsigned int k = 5; 
    double testpoint[3]; 

    vtkSmartPointer<vtkIdList> result = 
    vtkSmartPointer<vtkIdList>::New(); 

    for(vtkIdType n = 0; n < Output->GetNumberOfPoints(); n++) 
    { 

    result->Reset(); 
    Output->GetPoint(n,testpoint); 
    Tree->Update(); 
    std::cout << "Point: " << testpoint[0] << ", " << testpoint[1] << ", " << testpoint[2] << ": " << endl; 

    Tree->FindClosestNPoints(k, testpoint, result); 

    for(vtkIdType i = 0; i < k; i++) 
     { 
     vtkIdType point_ind = result->GetId(i); 
     double p[3]; 
     pointSource->GetOutput()->GetPoint(point_ind, p); 
     std::cout << "Closest point " << i+1 << ": Point " 
     << point_ind << ": (" << p[0] << ", " << p[1] << ", " << p[2] << ")" << std::endl; 
     } 

    } 
    //end test 

这是做什么我试图做...它打印出5个最接近点的兴趣的一个点,但尽管点的利益变化,5个最近的点保持不变。我假设我只是在我的代码中传递一个小细节,但我可能是错的。如果你需要更多的信息来帮助,只需要提问即可。

谢谢您的时间和帮助, 卢克^ h

回答

4

我发现这个问题......我用/ vtkPointSource生成随机点,我正在喂我点到该功能。我不知道为什么,但是这使得获得结果非常困难,但是一旦我传递了正确的信息,并在点Id数组和结果Id数组上放置了Reset(),它就像冠军一样工作。希望这会和其他人一样麻烦,因为我花了一些时间在这方面工作。

卢克ħ

vtkSmartPointer<vtkKdTreePointLocator> Tree = 
vtkSmartPointer<vtkKdTreePointLocator>::New(); 
Tree->SetDataSet(Output); 
Tree->BuildLocator(); 

vtkSmartPointer<vtkIdList> result = 
    vtkSmartPointer<vtkIdList>::New(); 

vtkSmartPointer<vtkIdList> point_ind = 
    vtkSmartPointer<vtkIdList>::New(); 

unsigned int k = 5; 

double testpoint[3]; 

for(vtkIdType n = 0; n < Output->GetNumberOfPoints(); n++) 
    { 

    Output->GetPoint(n,testpoint); 
    std::cout << "Point: " << testpoint[0] << ", " << testpoint[1] 
    << ", " << testpoint[2] << ": " << endl; 

    Tree->FindClosestNPoints(k, testpoint, result); 

    for(vtkIdType i = 0; i < k; i++) 
    { 

vtkIdType point_ind = result->GetId(i); 

    double p[3]; 

    Output->GetPoint(point_ind, p); 
std::cout << "Closest point " << i+1 << ": Point " 
    << point_ind << ": (" << p[0] << ", " << p[1] << ", " 
     << p[2] << ")" << std::endl; 

    } 

    result->Reset(); 
point_ind->Reset(); 

    }