我有一个文件，其中第一列是x坐标，第二列是y坐标，第三列是z坐标，第四列是与每个点相关的值。

我想画出这些点，每个点都应该根据第四栏的颜色。

我会用python来做这个。我在Windows上使用anaconda和vtk_可视化工具。

我有几百万分。我发现更快的方法是使用pythonvtk。

这是我现在掌握的代码：import vtk

import numpy as np

## DATA

# Generate random points w/ random RGB colors

n = 10**5

xyz = 100*np.random.rand(n, 3)

color = 10*np.random.rand(n, 1)

# Point size

point_size = 10

## COLORMAP

cmax = np.max(color)

cmin = np.min(color)

cmed = (cmax+cmin)/2

normalizzato = color / np.max( np.absolute(cmax), np.absolute(cmin) )

rgb = np.zeros((len(color), 3))

for i in range(0, len(color) ):

if color[i] >= cmed:

# Red

rgb[i][0] = 255*normalizzato[i]

if color[i] < cmed:

# Blue

rgb[i][2] = 255*normalizzato[i]

## VTK

# Create the geometry of a point (the coordinate)

points = vtk.vtkPoints()

# Create the topology of the point (a vertex)

vertices = vtk.vtkCellArray()

# Setup colors

Colors = vtk.vtkUnsignedCharArray()

Colors.SetNumberOfComponents(3)

Colors.SetName("Colors")

# Add points

for i in range(0, len(xyz)):

p = xyz[i]

id = points.InsertNextPoint(p)

vertices.InsertNextCell(1)

vertices.InsertCellPoint(id)

Colors.InsertNextTuple3(rgb[i][0], rgb[i][1], rgb[i][2])

# Create a polydata object

point = vtk.vtkPolyData()

# Set the points and vertices we created as the geometry and topology of the polydata

point.SetPoints(points)

point.SetVerts(vertices)

point.GetPointData().SetScalars(Colors)

point.Modified()

# Visualize

mapper = vtk.vtkPolyDataMapper()

if vtk.VTK_MAJOR_VERSION <= 5:

mapper.SetInput(point)

else:

mapper.SetInputData(point)

## ACTOR

# Create an actor

actor = vtk.vtkActor()

actor.SetMapper(mapper)

actor.GetProperty().SetPointSize(point_size)

axes = vtk.vtkAxesActor()

## RENDER

renderer = vtk.vtkRenderer()

# Add actor to the scene

renderer.AddActor(actor)

# Background

renderer.SetBackground(0.1, 0.2, 0.3)

# Reset camera

renderer.ResetCamera()

# Render window

renderWindow = vtk.vtkRenderWindow()

renderWindow.AddRenderer(renderer)

# Interactor

renderWindowInteractor = vtk.vtkRenderWindowInteractor()

renderWindowInteractor.SetRenderWindow(renderWindow)

# Begin interaction

renderWindow.Render()

renderWindowInteractor.Start()

这很快。

正如你所看到的，有一个颜色栏，但我不能得到适当的范围和颜色。有什么想法吗？

对于如何替换“COLORMAP”部分并使用引用真正COLORMAP的内容，您有什么建议吗？

我也试着用mayavi.mlab.point3d但它非常慢，而且vtk\u visualizer的代码如下：

^{pr2}$

但它比vtk慢，我不能改变点的大小，有一个颜色条和轴。在

谢谢