【作者】 闵小平；

【导师】 王博亮；

【作者基本信息】 厦门大学 ， 通信与信息系统， 2008， 博士

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

【Abstract】 The development of present-day science and technology was more and more characterized by the crossover of multi-disciplinarity.Virtual reality and simulation as a new area of rearch encompasses a wide spectrum of techniques including medicine、computer graphic technology、computer vision technology、mathematic analyzing and so on.Going with the development of computer technology and studying of medical image,cancer treatment technology go through a great change.Just as americal physical scientist J.Purdy described:radiation therapy has entrance a new era,the era of 3DRT(Three-Dimensional Radiation Therapy).Nasopharyngeal Cancer is one of the most common form of human cancer and the leading cause of human cancer deaths.Because of its characteristic,radiation therapy has become its major treatment option and Brachytherapy:which is a type of internal radiation thrapy given by placing radioactive sources directly into the cancerous tissue,is often used.Because there are many important tissues such as vas、nerve,lie around the nasopharynx,the dosimetric treatment planning must be designed carefully to achieves a desired dose distribution to the target,while sparing the critical structures.The image based treatment planning and treatment simulation for brachytherapy which can improve treatment quality,was studied by many researchers in recent years.The results of these research are of great value to treatment teaching,doctor training,clinical diagnosis and therapy.Firstly,we introduced the background of this area,then the key technology of image segmentation of soft tissues,three-dimensional modeling,calculation and optimization of HDR radiation dose,needle insertion simulation based on mass-springs model was researched.The main contents and contributions of this dissertation are as follows:(1) An algorithm based on the active shape model is proposed for the segmentation of medical image series.The edge of soft tissue in medical images usually are not very clearly,so it’s difficult to segment them only by the computer.But a tissue’s shape transform sequentially in medical image series,we asked doctor to select some images from image series and segment the tissue in these images,then we used the results of the segmentation as training maps,created a active shape model by PCA Algorithm.Next we used this model to segemnt other images.The experimental results show the accuracy of this method.(2) A parametric shape modeling using deformable superellipses for segmentation of medical image was proposed.Superellipses can represent the smooth and symmetrical contour perfectly, Most of anatomical structures having such characters can be approached by a superellipse.Prior shape information been got from a statistical modal analysis of a training set.This information was used to restrict the transformation of the surperellipse parameters.We used Hybrid Genetic Optimization Algorithm to find the optimal superellipse parameters,and snake algorithm embed with these shape information to segment the image.With shape guidance,this algorithm is less sensitive to initial contour placement and more robust even in the presence of large boundary gaps.The experiments show the efficiency of this method.(3) A new hybrid genetic algorithm based on a continuous-parameter GA with linear constraints (GENOCOP) and POWELL algorithm was proposed.The genetic algorithm has been widely used in many fields as an easy robust global search and optimization method.To optimizing the dose distribution,a new hybrid genetic algorithm based on a continuous-parameter GA with linear constraints(GENOCOP) and POWELL algorithm was proposed.The POWELL algorithm was modified to support the linear constraints and was embed in hybrid genetic algorithm as a local search operation,and the operators of GENOCOP were amended to improve its ability of global search.The simulation results show that compared to traditional GENOCOP methods,this hybrid genetic algorithm has obvious potential on many respects,such as convergence speed,solution accuracy,ability of global optimization.(4) An algorithm to optimize the dwell time of the radioactive source was proposed. Brachytherapy is the treatment of cancer by means of radioactive sources that are placed at short distances from the target cells.The computer-based HDR planning system of brachytherapy is becoming more and more widely used.To get the solution that best satisfies the clinical dose constraints,we proposed an algorithm to optimize the dwell time of the radioactive source.Firstly we considered the dwell time as a continuous function of the dwell position,computed the dose values of referenced dose points.The minimal difference between these dose values and demanded dose values was used as the objective of an optimization algorithm to get the optimized function parameters.The function curve can be separated to several sects.At last the optimized curve was dispersed using compound trapezoid formula to get the dwell positions and dwell times.This method not only avoids the emergence of negative dwell time,but also reduces the dwell time gradient.And during the process of dispersing the curve,we can get different dwell times and dwell positions.This makes the brachytherapy treatment plan more flexible.(5) Simulation of needle insertion based on mass-spring model was developed to facilitate surgeon training and planning for brachytherapy.Inserting a needle into soft tissues causes the tissues to displace and deform:ignoring these effects during seed implantation leads to imprecise dose distribution.We construct a 2-D dynamic mass-spring model using 400 springs and simulated the needle insertion.The deformation of the tissue and the force can be calculated.Simulation achieved 100 flames per second on a 2Ghz Pentium PC.更多还原