3D-форма с использованием VTK:

Question

У меня есть трехмерный массив (с именем V). он содержит информацию о вокселе. Где бы то ни было присутствующего воксела, значение V [i] [j] [k] равно 1. Если когда-либо нет воксела, значение V [i] [j] [k] равно 0 Я хочу визуализировать эту форму с помощью VTK.

я написал этот код:

MyVTKPointPlotter pointPlotter;

for(int i=0;i<x_count;i++) 
{ 
    for(int j=0;j<y_count;j++) 
    { 
     for(int k=0;k<z_count;k++) 
     { 
     if(V[i][j][k] != 0) 
     { 
      pointPlotter.PlotPoint(i,j,k,128,128,128); 
     } 
     } 
    } 
} 

Примечание: функции MyVTKPointPlotter все полученные по этой ссылке: http://nawigacjarobota.googlecode.com/svn-history/r10/trunk/wykObMAT/myVTKPointPlotter.cpp

Проблема в этом коде является то, что все точки становятся построены, но рендеринг крайне медленно. Кроме того, точки являются 2 размерными точками и, следовательно, они не имеют толщины. Поэтому, когда я поворачиваю свой объект, я могу видеть точки без толщины (точки видны как диски).

Может ли кто-нибудь сказать мне, как визуализировать эту трехмерную форму?

Answer 1

Если ваши воксели вместе образуют большие концентрации, вы можете использовать визуализацию изоповерхности: http://www.evl.uic.edu/aspale/cs526/final/3-5-2-0.htm. Если предыдущий сложный, google it, существует немало примеров того, как это сделать.

ИЛИ

Если белые вокселей разбросаны, я думаю, что вы можете создать облако точек.

У меня был этот код в Python, который я изменяю, чтобы соответствовать вашему коду вопросов. Перевести на свой Languaje (это не тестировался, но поможет вам идею):

points = vtk.vtkPoints() 
colors = vtk.vtkUnsignedCharArray() 
colors.SetNumberOfComponents(3) 

#Create the point cloud 
for i in xrange(0, x_count): 
    for j in xrange(0, y_count): 
     for k in xrange(0, z_count): 
     if V[i][j][k] != 0: 
      #If they are many points better use the faster combination of 
      #SetPointCount and SetPoint methods 
      points.InsertNextPoint(i, j, k) 
      colors.InsertNextTuple3(128, 128, 128) 



pointsPolyData = vtk.vtkPolyData() 
vertexFilter = vtk.vtkVertexGlyphFilter() 
polyData = vtk.vtkPolyData() 

mapper = vtk.vtkPolyDataMapper() 
actor = vtk.vtkActor() 
actor.GetProperty().SetPointSize(5) 

pointsPolyData.SetPoints(points) 
vertexFilter.SetInputConnection(pointsPolyData.GetProducerPort()) 
polyData.ShallowCopy(vertexFilter.GetOutput()) 
polyData.GetPointData().SetScalars(colors) 
mapper.SetInputConnection(polyData.GetProducerPort()) 
actor.SetMapper(mapper) 


//Then you'll have to add the actor to your renderer, of course 
myRenderer.AddViewProp(actor)