第 31 卷 第 3 期
2014 年 9 月
上海第二工业大学学报
JOURNAL OF SHANGHAI SECOND POLYTECHNIC UNIVERSITY
Vol. 31 No. 3
Sep. 2014
文章编号: 1001-4543(2014)03-0211-04
基基基于于于 VTK 的的的足足足部部部骨骨骨骼骼骼三三三维维维可可可视视视化化化研研研究究究
戴 虹
(上海第二工业大学电子与电气工程学院, 上海 201209)
摘 要: 基于 CT 图像序列利用 VTK(Visualization Toolkit, 可视化工具箱) 实现了人体足部骨骼的三维表面重建。采
用最大熵阈值算法进行足部骨骼图像分割, 将该阈值作为骨骼等值面抽取值; 在 VC++ 中利用 VTK 编程实现了基
于移动立方体法的足部骨骼表面重建, 并对所得三维表面进行了优化; 将所得结果存为三维模型常用的 STL 文件格
式。实验结果表明, 此方法获得了理想的可视化效果, 可在足部疾病诊断和治疗、手术仿真等方面发挥重要作用。
关键词: VTK; 足部骨骼; 三维可视化; 面绘制; 最大熵阈值算法
中图分类号: TP391.4 文献标志码: A
0 引言
图像的三维可视化技术是指利用图像的二维投
影恢复物体的三维信息, 目前广泛应用于医学和工
程领域。在临床医学上利用可视化技术, 将一系列
二维断层图像, 如 CT 与 MRI(磁共振) 图像序列重
建成三维模型, 医生根据模型准确地获得病变组织
的位置、大小和体积等信息, 并进行手术仿真, 可提
高手术的成功率
[1]
。目前常用 VC++ 软件进行医学
图像编程工作, 但其复杂度较高, 开发周期长且三维
显示速度较慢。在使用通用三维图形接口 OpenGL
(Open Graphics Library, 开放的图形接口) 后, 三维绘
制的速度和质量有了较大提高, 但编程困难问题并
没有根本解决。美国 Kitware 公司于 1998 年推出的
VTK (Visualization Toolkit, 可视化工具箱), 是一种
基于 OpenGL 的用于三维图形学、图像处理及三维
可视化的开发包
[2]
, 可用简单代码实现图像的三维
重建。
本文首先基于最大熵阈值法求得足部 CT 图像
中骨骼的分割阈值, 将其作为等值面抽取值, 再利用
VTK 编程实现足部骨骼的三维表面重建。实验结
果表明, 本文提出的方法获得了较为理想的重建结
果。
1 基于最大熵阈值法的足部 CT 图像
骨骼分割
在对足部骨骼进行三维表面绘制之前须求得
骨骼等值面抽取值, 也就是足部骨骼图像的分割阈
值, 大于该阈值的图像像素被认定为骨骼。由于凭
经验手动设置阈值准确率不高且效率低下, 所以采
用一种自动的最大熵阈值分割算法: 由信息论的定
义, 熵是平均信息量的表征, 根据最大熵原则, 用灰
度的一维熵求取阈值即选取一个阈值 T, 把图像分
成背景和目标两类区域, 使图像用这个阈值分割出
的两部分的一维熵最大, 可实现背景和目标的最佳
分类
[3]
。
该算法的步骤如下: 设一幅足部图像 f 的灰度
级为 0 ∼ 255, 灰度级为 i 的像素个数为 n
i
, 总像素
个数为 N, 则各灰度级出现的概率 (即归一化直方
图) 为 P
i
=
n
i
N
(i = 0, 1, · · · , 255)。设分割骨骼与
其他组织的阈值为 T , 则获得目标 A 和背景 B 两
类, 它们的概率分布分别为 A:
P
0
P
,
P
1
P
, · · · ,
P
T
P
, B:
P
T +1
1 − P
,
P
T +2
1 − P
, · · · ,
P
255
1 − P
, 其中 P =
T
X
i=0
P
i
, 则 A 和
B 的熵分别为:
收稿日期: 2014-03-02
通讯作者: 戴 虹 (1977–), 女, 上海人, 副教授, 博士, 主要研究方向为信号与信息处理、医学图像处理。
电子邮箱 daihong@sspu.edu.cn。
基金项目: 上海第二工业大学校级重点学科 (通信与信息系统) 第四期建设项目 (No. XXKZD1302) 资助