鼻咽癌近距离放射治疗手术计划和仿真系统的关键技术研究

摘要

现代科学技术的发展越来越体现多门学科的交叉和渗透。医学虚拟现实和仿真技术作为正在发展起来的研究方向，是集医学、计算机图形学、计算机视觉、数学分析等学科为一体的新型交叉研究领域。近20年来，随着计算机技术和医学影像学的发展，肿瘤放射治疗技术也经历了很大的变革。正像著名的美国医学物理学家J．Purdy所描述的那样，放射治疗已进入一个崭新的令人振奋的时代，即三维放射治疗时代(Three—DimensionalRadiationTherapy，3DRT)。 鼻咽癌是我国的高发恶性肿瘤，由于其解剖部位、病理类型、生物学行为等的特殊性，放射治疗是其治疗的主要手段。近年来国内外开展经颌下咽旁间隙插植组织间近距离放疗技术，是近距离治疗技术一项重要进展。但是由于鼻咽癌靶区附近存在着众多正常的神经、血管等重要解剖结构，因此如果能结合医学虚拟现实和仿真技术精确地明确靶区、合理设计剂量计划，将可以进一步做到手术的精确化，保护正常组织、提高疗效和生活质量。同时通过构建手术仿真系统，使医生可以借助虚拟环境进行手术预演，从而使研究结果可以应用于手术计划、手术训练、医学教学等医学科研、临床和教学培训各个环节，具有非常重要的意义。 本文首先简要介绍了近距离放射治疗计划与手术仿真方面的背景。然后对人体鼻咽部组织CT图像的分割、三维重建、网格化简、近距离放射源的剂量计算、剂量优化模型和方法、基于质点弹簧模型的针管插植仿真、基于OBB包围盒的碰撞检测等关键技术进行了研究。 本文工作的主要创新之处在于： 1．提出了一种基于PCA形变模型的医学序列图像分割算法 一个组织在CT序列图像中其形状变化往往是连续的，因此，我们先对病人某个组织的CT扫描序列图像中的相隔的几幅图像用其他方法做一个分割(由于边界不清晰，一般需要医生的参与)。然后将分割后的组织提取出来，生成其距离图。再对这几幅距离图使用主成分分析算法(PrincipleComponentAnalysis)进行分析，求其平均距离图及其协方差矩阵的特征值和特征矢量。然后利用平均距离图和这些特征值和特征矢量生成一个在参数作用下可以连续变形的模板，仿真组织在序列图像中的形状变化。再将这个模板采用互信息法和主动轮廓模型法(snake)去匹配分割序列中其他未分割的图像，从而达到一个比较理想的分割效果。 2．提出了一种基于超椭圆形变模型的医学图像分割算法 椭圆变形模板分割算法对目标组织的轮廓形状有要求，在医学图像中，只要组织的轮廓比较光滑、结构对称，采用超椭圆模板来进行匹配分割往往能够取得比较好的结果，而且其没有序列图像分割这一限制。算法先对待分割组织的典型形状用超椭圆进行拟合，获得先验信息。然后使用优化算法获得超椭圆与目标组织的最佳匹配。然后再使用主动轮廓模型法进行进一步调整，得到最终分割结果。由于超椭圆的控制参数比较少，因此相对于标准的主动形状模型(ASM)，它的优点在于其无论是变形还是形状参数的统计都比较方便，特别是省去了进行形状参数统计分析时各形状图的大小、角度、位置对准。 3．本文提出了一种基于GENOCOP和POWELL算法并支持线性约束的混合遗传算法 本文对标准的POWELL算法进行了修改，使其支持线性约束，并将其和遗传算法GENOCOP进行结合，构成一种新的支持线性约束的混合遗传算法。POWELL算法经过修改后，并不能找到线性约束下的函数极值点，只能达到最优解附近区域。本文根据对其原因的分析已及对实验结果的统计分析，对GENOCOP算法的变异、交叉算子进行了进一步的改进，降低遗传算子的收敛能力，增强其全局搜索能力。新的算法融合了遗传算法的全局搜索能力和POWELL的局部优化能力，并利用遗传算法中交叉算子的特性有效的解决了POWELL算法修改后所带来的问题，从而将两者有机的结合了起来。 4．本文提出了一种连续驻留位置和驻留时间的近距离放疗剂量优化方法 本文提出了一种连续驻留位置和驻留时间的近距离放疗剂量优化方法。首先将驻留t看成是位置x的一个连续函数，建立各参考剂量点剂量计算的积分模型，计算每个参考剂量点的剂量。然后用计算剂量和目标剂量值之差的加权平方和为目标函数，用优化算法求解最优的曲线参数。在优化的过程中由于实际的情况的复杂性，可以根据不同的情况不同的特点对曲线进行分段。得到优化的曲线后再用积分的数值逼近方法将其离散化，得到最终的插值点源治疗计划。算法不仅解决了负的驻留时间问题，还让相邻驻留位置的驻留时间更加平滑。同时，在最后的离散化过程中，可以得到不同的驻留位置和驻留时间结果，使计划具有更好的灵活性。 5．应用质点弹簧模型实现了二维软组织的针管插植仿真 鼻咽部近距离放射治疗时放射源需经针管插植的路径送到肿瘤组织附近。对针管插植过程进行仿真，其结果可以用于医生的训练，同时可以估计插植引起的软组织变形，使剂量计划更准确。本文使用OBBTREE进行碰撞检测，使用质点弹簧模型对软组织进行了建模，对插植手术进行了仿真，模拟了针管插入二维软组织时的组织变形和受力。由于质点弹簧模型模型简单，实现容易，计算速度比较快，因此能够满足仿真的实时性要求。

目录

文摘

英文文摘

声明

第一章绪论

§1.1虚拟人、虚拟器官及仿真技术的发展

§1.2鼻咽部近距离放射治疗的发展和意义

§1.3近距离放射治疗手术计划与插植手术仿真国内外研究现状

§1.4本文的主要研究内容

第二章 鼻咽部组织的图像分割

§2.1当前比较常见的分割方法概述

§2.2鼻咽部组织图像特征及一些经典方法的分割结果

§2.3主动轮廓模型分割的原理及发展

2.3.1基本Snakes模型

2.3.2轮廓模型的改进和发展

§2.4基于主动形变模型的医学序列图像的分割方法

2.4.1距离图的形状模型生成

2.4.2相关配准

2.4.3主动轮廓分割

2.4.4实验结果

§2.5基于超椭圆的变形模板分割

2.5.1超椭圆曲线方程及其变形参数

2.5.2基于超椭圆模板的变形匹配

2.5.3初始轮廓的进一步修正

§2.6本章小结

第三章组织图像的三维重建

§3.1三维重建算法算法简介

§3.2 VTK开发工具包简介

§3.3鼻咽部及其周围组织图像序列的三维重建

§3.2网格化简

§3.3本章小结

第四章 鼻咽癌近距离放疗剂量的优化

§4.1近距离放射治疗及计划概述

§4.2基于蒙特卡罗方法的剂量计算

4.2.1放射物理原理

4.2.2蒙特卡罗剂量计算原理

4.2.3基于EGS的剂量计算

§4.3剂量计算的经验计算方法

4.3.1与近距离放疗的剂量计算有关的物理量

4.3.2 TG43报告推荐的剂量计算模型

4.3.3 ICWG推荐点源的计算公式

4.3.4结果

§4.4剂量优化

4.4.1以参考点剂量约束为目标的优化原理及常见优化方法

4.4.2连续时间连续驻留位置的优化模型

4.4.3复形调优法优化原理及优化结果

4.4.4混合遗传优化算法

4.3.6混合遗传算法的剂量的优化求解

4.4.6连续时间连续驻留位置优化模型的进一步改进

4.4.6离散过程

§4.5本章小结

第五章 插植手术仿真

§5.1软组织变形仿真概述

§5.2软组织变形计算建模

§5.3基于弹簧质点模型的二维软组织变形仿真

5.3.1质点受力分析

5.3.2变形模型的求解

5.3.3插植过程中外力的模拟

§5.4仿真过程中的碰撞检测

5.4.1包围盒算法的比较

5.4.2基于OBB包围盒的碰撞检测

§5.5本章小节

第六章总结与展望

§6.1本论文内容总结

§6.2需进一步研究的方向

博士期间参加的项目和研究成果

致谢

参考文献