MedicalDemo4 体素绘制与展示

一：主要的知识点

1、说明

本文只是教程内容的一小段，因博客字数限制，故进行拆分。主教程链接：vtk教程——逐行解析官网所有Python示例-CSDN博客

2、知识点纪要

本段代码主要涉及的有①vtkFixedPointVolumeRayCastMapper体素渲染，②体渲染的其它内容

二：代码及注释

import vtkmodules.vtkInteractionStyle import vtkmodules.vtkRenderingOpenGL2 from vtkmodules.vtkCommonColor import vtkNamedColors from vtkmodules.vtkIOImage import vtkMetaImageReader from vtkmodules.vtkRenderingCore import vtkColorTransferFunction, vtkRenderer, vtkRenderWindow, \ vtkRenderWindowInteractor, vtkVolume, vtkVolumeProperty from vtkmodules.vtkRenderingVolume import vtkFixedPointVolumeRayCastMapper from vtkmodules.vtkCommonCore import vtkLookupTable from vtkmodules.vtkCommonDataModel import vtkPiecewiseFunction def main(): colors = vtkNamedColors() colors.SetColor('BkgColor', [51, 77, 102, 255]) file_name = "Data/FullHead.mhd" reader = vtkMetaImageReader() reader.SetFileName(file_name) reader.Update() """ vtkFixedPointVolumeRayCastMapper TK 体绘制（Volume Rendering）中非常经典的体渲染算法类， 尤其在需要 高精度但又相对高效的 CPU 光线投射渲染 时使用。 它通过计算光线穿过体积数据时的累积颜色与透明度，实现 半透明体渲染 效果，通俗的讲，它不是绘制表面，而是“看透体积”的渲染器 """ volume_mapper = vtkFixedPointVolumeRayCastMapper() volume_mapper.SetInputConnection(reader.GetOutputPort()) # 建立标量与颜色之间的对应关系 volume_color = vtkColorTransferFunction() volume_color.AddRGBPoint(0, 0.0, 0.0, 0.0) volume_color.AddRGBPoint(500, 240.0 / 255.0, 184.0 / 255.0, 160.0 / 255.0) volume_color.AddRGBPoint(1000, 240.0 / 255.0, 184.0 / 255.0, 160.0 / 255.0) volume_color.AddRGBPoint(1150, 1.0, 1.0, 240.0 / 255.0) """ vtkPiecewiseFunction 定义标量值（scalar）到“单个实数值”（通常是透明度或强度）的映射关系， 并通过分段线性插值（piecewise linear interpolation）实现平滑过渡 """ # 标量透明度 volume_scalar_opacity = vtkPiecewiseFunction() volume_scalar_opacity.AddPoint(0, 0.00) volume_scalar_opacity.AddPoint(500, 0.15) volume_scalar_opacity.AddPoint(1000, 0.15) volume_scalar_opacity.AddPoint(1150, 0.85) """ 将体素强度的梯度（变化率）作为输入，将不透明度（Opacity）作为输出 在3D体积数据中，梯度测量的是体素强度在单位距离内变化的速度和方向 低梯度(0~90): 发生在均匀的区域内,例如纯净的空气，软组织内部等 高梯度(90, 100) 发生在组织边界，例如软组织过度到谷歌表面，组织过度到空气的表面 如果没有梯度不透明度函数，体积渲染的结果看起来就像一团“云雾” 标量不透明度函数 (volume_scalar_opacity) 已经为骨骼设置了较高的不透明度，但这仍可能导致骨骼看起来模糊且缺乏锐利的表面 解决方案: 通过这个梯度函数，我们明确告诉渲染器：“只有在体素强度发生急剧变化的地方（即物体的表面边界），才将其设置为高不透明度。” 有效地抑制了体素内部的“噪声”或细节，并突出了体素间的过渡边界 """ # 梯度透明度 volume_gradient_opacity = vtkPiecewiseFunction() volume_gradient_opacity.AddPoint(0, 0.0) volume_gradient_opacity.AddPoint(90, 0.5) volume_gradient_opacity.AddPoint(100, 1.0) """ vtkVolumeProperty 是 VTK 体绘制（Volume Rendering）管线中最核心的类之一， 它定义了 体数据在渲染时的视觉属性（颜色、透明度、光照、插值方式等） 可以把它理解为三维体的“材质(Material)” + “视觉属性(Visual Property)” 控制器 """ volume_property = vtkVolumeProperty() volume_property.SetColor(volume_color) volume_property.SetScalarOpacity(volume_scalar_opacity) volume_property.SetGradientOpacity(volume_gradient_opacity) """ SetInterpolationTypeToLinear 设置体素差值方式为线性插值 在体积渲染的光线投射过程中，光线会沿着一条直线穿过体素（Voxel）。 Mapper 需要在光线上的许多非整数体素坐标点进行采样（即在体素的中心之间取值） 如果不设置或设置为最近邻插值，采样点的值将直接采用其最近的那个体素中心的值。 这会导致渲染结果出现明显的方块状、锯齿状的伪影，因为体素的边界被清晰地显示出来 """ volume_property.SetInterpolationTypeToLinear() volume_property.ShadeOn() # 开启着色 volume_property.SetAmbient(0.4) # 设置环境光， 控制在没有直接光源照射到的地方，物体表面仍能看到的基本亮度 volume_property.SetDiffuse(0.6) # 控制物体表面对光的普遍反射。这是形成物体主要颜色和形状的关键因素 0.6 使物体的基本颜色和纹理细节清晰可见 volume_property.SetSpecular(0.2) # 设置镜面反射，意味着物体表面有一定的光泽，但不会像镜面那样反射出强烈的高光 """ vtkVolume 体对象类(Volume Actor) 是体渲染中最顶层、最关键的“场景对象”之一，它是你最终放进 Renderer 里看到的“3D体对象” 一个 vtkVolume 对象由两个主要部分组成 Property： vtkVolumeProperty 控制体渲染的外观（颜色、透明度、光照等) 决定渲染的算法 Mapper: vtkVolumeMapper 控制体渲染算法（如何从数据生成影像） 决定渲染的外观 可搭配的Mapper vtkSmartVolumeMapper 自动选择最佳渲染器 推荐首选 vtkGPUVolumeRayCastMapper GPU 光线投射 快速、现代GPU加速 vtkFixedPointVolumeRayCastMapper CPU 光线投射 CPU版，可用于无GPU环境 """ volume = vtkVolume() volume.SetMapper(volume_mapper) volume.SetProperty(volume_property) render = vtkRenderer() render.AddViewProp(volume) renderWindow = vtkRenderWindow() renderWindow.AddRenderer(render) iren = vtkRenderWindowInteractor() iren.SetRenderWindow(renderWindow) camera = render.GetActiveCamera() c = volume.GetCenter() camera.SetViewUp(0, 0, -1) camera.SetPosition(c[0], c[1] - 400, c[2]) camera.SetFocalPoint(c[0], c[1], c[2]) camera.Azimuth(30.0) camera.Elevation(30.0) render.SetBackground(colors.GetColor3d('BkgColor')) # Increase the size of the render window renderWindow.SetSize(640, 480) renderWindow.SetWindowName('MedicalDemo4') # Interact with the data. iren.Start() if __name__ == '__main__': main()