Question

我使用VTK工具包加载OBJ文件和vtkCutter以通过播放切割数据集，然后绘制剪切的轮廓。对于大型对象，这可能会变得非常慢，正如另一位用户在VTK Users Forum中指出的那样。

有没有办法让切割机使用分层数据结构来获得更好的性能？

这是代码：

#include <vtkSmartPointer.h>
#include <vtkCubeSource.h>
#include <vtkPolyDataMapper.h>
#include <vtkPlane.h>
#include <vtkCutter.h>
#include <vtkProperty.h>
#include <vtkActor.h>
#include <vtkRenderer.h>
#include <vtkRenderWindow.h>
#include <vtkRenderWindowInteractor.h>
#include <vtkOBJReader.h>

int main(int argc, char *argv[])
{
    // Parse command line arguments
    if (argc != 2) {
        std::cout << "Usage: " << argv[0] << " Filename(.obj)" << std::endl;
        return EXIT_FAILURE;
    }

    std::string filename = argv[1];
    vtkSmartPointer<vtkOBJReader> obj = vtkSmartPointer<vtkOBJReader>::New();
    obj->SetFileName(filename.c_str());
    obj->Update();

    vtkSmartPointer<vtkPolyDataMapper> mapper = vtkSmartPointer<vtkPolyDataMapper>::New();
    mapper->SetInputConnection(obj->GetOutputPort());

    // Create a plane to cut,here it cuts in the XZ direction (xz normal=(1,0,0);XY =(0,0,1),YZ =(0,1,0)
    vtkSmartPointer<vtkPlane> plane = vtkSmartPointer<vtkPlane>::New();
    plane->SetOrigin(0, 0, 0);
    plane->SetNormal(1, 0, 0);

    // Create cutter
    vtkSmartPointer<vtkCutter> cutter = vtkSmartPointer<vtkCutter>::New();
    cutter->SetCutFunction(plane);
    cutter->SetInputConnection(obj->GetOutputPort());
    cutter->Update();

    vtkSmartPointer<vtkPolyDataMapper> cutterMapper = vtkSmartPointer<vtkPolyDataMapper>::New();
    cutterMapper->SetInputConnection(cutter->GetOutputPort());

    // Create plane actor
    vtkSmartPointer<vtkActor> planeActor = vtkSmartPointer<vtkActor>::New();
    planeActor->GetProperty()->SetColor(1.0, 1, 0);
    planeActor->GetProperty()->SetLineWidth(2);
    planeActor->SetMapper(cutterMapper);

    // Create cube actor
    vtkSmartPointer<vtkActor> cubeActor = vtkSmartPointer<vtkActor>::New();
    cubeActor->GetProperty()->SetColor(0.5, 1, 0.5);
    cubeActor->GetProperty()->SetOpacity(0.5);
    cubeActor->SetMapper(mapper);

    // Create renderers and add actors of plane and cube
    vtkSmartPointer<vtkRenderer> renderer = vtkSmartPointer<vtkRenderer>::New();
    renderer->AddActor(planeActor); //display the rectangle resulting from the cut
    renderer->AddActor(cubeActor); //display the cube

    // Add renderer to renderwindow and render
    vtkSmartPointer<vtkRenderWindow> renderWindow = vtkSmartPointer<vtkRenderWindow>::New();
    renderWindow->AddRenderer(renderer);
    renderWindow->SetSize(600, 600);

    vtkSmartPointer<vtkRenderWindowInteractor> interactor = vtkSmartPointer<
            vtkRenderWindowInteractor>::New();
    interactor->SetRenderWindow(renderWindow);
    renderer->SetBackground(0, 0, 0);
    renderWindow->Render();

    interactor->Start();

    return EXIT_SUCCESS;
}

Answer 1

vtkCutter使用任意复杂的func(x,y,z)对网格进行切片，这里使用一个简单的平面来描述该函数，这是一个常见且覆盖良好的特殊情况，因为切割计数位于一个简单的平面，因此将是一个简单的（平面）多边形。

这些通用实现通常会耗费大量的CPU时间，因为预期在vtkCutter的情况下会出现所有特殊的多边形切割情况。
在VTK的庞大类层次结构中调用虚函数也会导致速度减慢。没有特殊的黑客攻击，它只依赖于编译器来优化循环中的虚函数指针查找，而VTK在一个或多个嵌套循环中多次调用虚函数（例如过滤函数）。有关相关信息，请参阅此处：about the cost of virtual function
VTK几乎随处可见double，即使可以使用float s。转换和高精度也会增加一点计算和内存开销。
VTK（5.8）没有明确涉及像SSE，afaik这样的SIMD操作。
...

搜索以下主题：

尽管在CPU上执行此操作，也可以在变换反馈过程中使用OpenGL几何着色器来提取由切割平面确定的切割轮廓。在OpenCL中执行此操作也是可能的，但是，如果没有基于GPU的计算设备可用，它可能比C或C ++实现慢。

要渲染网格，可以使用任何支持OpenGL 3+的渲染器：

Ogre3D
Unity3D
Irrlicht
OSG
一个简单的自制OpenGL 3渲染器。
...

更多：What is the best way to have realtime 3D rendering in an engineering application?

VTK工具包 - vtkCutter性能

1 个答案: