医学图象三维重建及可视化技术研究.ppt

医学图象三维重建及可视化医学图象三维重建及可视化技术研究技术研究秦绪佳秦绪佳浙江大学浙江大学浙江大学浙江大学CADCADCADCAD＆＆＆＆CGCGCGCG国家重点实验室国家重点实验室国家重点实验室国家重点实验室2001.9.282001.9.281 绪论绪论　　　　       1.1 引言引言   1.2 基于三维数据的建模与可视化基于三维数据的建模与可视化        1.三维数据的来源与分类    2.三维数据建模及可视化研究内容  数据预处理         建模         绘制与显示    3.数据建模技术综述 三维几何模型 　基于三维数据的建模方法  1)　基于断层轮廓的表面重建  2)　基于体素的等值面重建  3)　几何变形模型  4)　体素建模 　 1.3 1.3 医学图象三维重建技术综述医学图象三维重建技术综述    1 医学图象的预处理     2  医学图象的分割     3 三维重建方法        4  模型的网格简化  1.4 1.4 医学图象三维重建在医疗中的应用医学图象三维重建在医疗中的应用      1 在医疗诊断中的应用   2 在手术规划及放射治疗规划中的应用   3 在整形与假肢外科中的应用   4 在虚拟手术及解剖教育中的应用 1.5 1.5　论文背景及主要工作　论文背景及主要工作 　   1 论文背景及研究意义 　   2 本文的主要工作          1) 图象预处理，组织器官分割与提取           2) MC、MT算法构建表面几何模型           3) 模型表面网格简化，剖切与开窗           4) 由表面几何模型转换成实体几何模型           5) 适用于适形放射治疗规划的医学图象                      三维重建系统的开发 2医学图象预处理与人体组织的分割医学图象预处理与人体组织的分割　　　　 　　医学图象预处理　医学图象预处理　　　  分割流程分割流程  图图2.1 2.1 三维医学图象分割流程三维医学图象分割流程  2.2 2.2 CTCT、、MRIMRI图象的获取与输入图象的获取与输入 2.2 2.2 二维图象处理与规则体数据封装二维图象处理与规则体数据封装   1 二维图象滤波         (1) 邻域平均法 (2) 中值滤波法 (3) 保持边缘滤波法    2 断层图象间插值    3 三维规则体数据封装      (1) 内存记录方式   (2) 体数据文件格式  图图2.4 2.4 体数据内存记录方式体数据内存记录方式   2.3 交互分割过程交互分割过程    1 三维图象二值化         二值化结果 图图2.8 2.8 断层图象二值化结果断层图象二值化结果   2 数学形态学操作进行区域修整           (1) 二值形态学操作简述 ( (a) a) 原图象原图象 ( (b)b)结构元素结构元素 ( (c)c)对原图象的腐蚀对原图象的腐蚀 ( (d)d)对腐蚀图象的膨胀对腐蚀图象的膨胀 图图2.9 2.9 开启操作开启操作 ( (a)a)原图象原图象 ( (b)b)结构元素结构元素 ( (c)c)对原图象的膨胀对原图象的膨胀 ( (d)d)对膨胀图象的腐蚀对膨胀图象的腐蚀 　　图　　图2.102.10　闭合操作　闭合操作    3 种子填充法进行组织提取  图图2.12 2.12 交互分割结果交互分割结果 图图2.13 2.13 对分割区域的重建对分割区域的重建 3 基于规则体数据的三维表面模型的基于规则体数据的三维表面模型的构建构建 3.2.1 体素模型 ( (a) a) 方向无关的三线性插值模型方向无关的三线性插值模型 ( (b) b) 方向有关的三线性插值模型方向有关的三线性插值模型图图3.1 3.1 体素模型体素模型  3.2.2 等值面（IsoSurface）定义       1　三线性插值结果　　　   2　等值面定义等值面是三次曲面 3.3 移动立方体（Marching Cubes）算 法抽取等值面   1 1 MCMC算法的基本原理算法的基本原理 a a 体素中等体素中等值值面剖分方式的确定面剖分方式的确定 1 1））   如立方体顶点的数据值如立方体顶点的数据值≥≥等值面的值，则定义该顶点位等值面的值，则定义该顶点位于等值面之外，记为于等值面之外，记为“0”“0”；； 2 2））   如如立立方方体体顶顶点点的的数数据据值值＜＜等等值值面面的的值值，，则则定定义义该该顶顶点点位于等值之内，记为位于等值之内，记为“1”“1”。

8 8个个顶顶点，每个点，每个顶顶点共有点共有2 2个状个状态态，因此共，因此共256256种种组组合状合状态态     根据互根据互补对补对称性称性, , 256 128 256 128     根据旋根据旋转对转对称性，称性， 256 15 256 15           (1)体素中由三角片逼近的等值面计算          (2)三角片各顶点法向量计算      2 等值面连接方式上的二义性 ( (a)a)连接方式二义性的二维表示连接方式二义性的二维表示  ( (b) b) 连接方式二义性的三维表示连接方式二义性的三维表示 图图3.5 3.5 拓扑不一致造成孔隙拓扑不一致造成孔隙 图图3.4 3.4 MCMC方法的二义性方法的二义性  3.3.3 渐近线判别法消除二义性 3.3.4 MC算法的重建结果及分析 (a)(a)256×256×109256×256×109MRIMRI表皮重建表皮重建 （（b b））128×128×93CT128×128×93CT颅骨重建颅骨重建 （（c c））128×128×93CT128×128×93CT表皮重建表皮重建 三角面片：三角面片：696889顶点：顶点：347322三角面片：三角面片：187559顶点：顶点：94015三角面片：三角面片：137799顶点：顶点：69331 图图3.83.8　　MCMC算法重建的表面模型算法重建的表面模型  3.4 移动四面体（Marching Tetrahedra）算法抽取等值面 　　3.4.1 MT算法的基本原理 图图3.9 3.9 立方体的四面体剖分立方体的四面体剖分  图图3.10 3.10 四面体中的等值面四面体中的等值面  3.4.2 四面体剖分的一致性处理 图图3.11 3.11 立方体剖分为四面体立方体剖分为四面体的不同方式的不同方式 图图3.12 3.12 两相邻立方体剖分不一致时共有两相邻立方体剖分不一致时共有面的剖分情况面的剖分情况 图图3.133.13相邻立方体公共面上的剖分一致性相邻立方体公共面上的剖分一致性  3.4.3 相关性处理加速MT重建速度      1 体素内的相关性处理    2 体素间的相关性处理 图图3.14 3.14 剖分后立方体的顶点及棱边编号剖分后立方体的顶点及棱边编号  3.4.4 MT算法的重建结果及分析 ( (a) 128×128×113CTa) 128×128×113CT颅骨重建颅骨重建 （（b b））104×185×220CT104×185×220CT脚脚骨骼重建骨骼重建 ( (c) 128×128×113CTc) 128×128×113CT表皮重建表皮重建 三角面片：三角面片：423998顶点：顶点：211905三角面片：三角面片：365858顶点：顶点：183056三角面片：三角面片：331290顶点：顶点：165808图图3.153.15　　MTMT算法重建的表面模型算法重建的表面模型4.三维模型的网格简化与模型的剖切三维模型的网格简化与模型的剖切 4.2 基于边收缩的网格简化算法        1 网格简化算法简述  （（1 1）抽样（）抽样（SamplingSampling）） （（2 2））自适应细分自适应细分( (Adaptive subdivision)Adaptive subdivision) （（3 3））删除删除 ( (Decimation)Decimation) （（4 4））顶点合并顶点合并 （（Vertex mergingVertex merging））    2 基于边收缩的网格简化算法Hoppe采用显式能量函数E（M）来度量简化网格与原始网格的逼近度［Hoppe96］：其中Edist(M)为M的距离能量，它定义为点集到网格的距离平方：Espring(M)为弹性能量，这相当于在的每条边上均放置一条弹性系数为k的弹簧，即： Escalar(M)度量M的标量属性的精度，而Edisc(M)则度量了M上视觉不连续的特征线（如边界线、侧影轮廓线等）的几何精度。

 边收缩过程示意图边收缩过程示意图 ( (a)a)收缩前收缩前   ( (b)b)收缩后收缩后 图图4.14.1边收缩过程边收缩过程   4 网格简化结果 MT重建结果　　　重建结果　　　          简化简化50％　　　　　　％　　　　　　   简化简化90％　　　　简化％　　　　简化90％表面绘制％表面绘制  MC重建结果　　　　　简化重建结果　　　　　简化50％　　　　　简化％　　　　　简化85％　　　％　　　   简化简化85％表面绘制％表面绘制4.3 三维模型的剖切    4.3.1 模型三角面片的剖切处理    1 平面方程的确定 ax + by + cz + d = 0    2 三角面片与剖切平面的求交检测 定定义义空空间间一点一点P P（（X,Y,ZX,Y,Z））, , 定义定义““距离距离””D D：： D = D = aX aX + + bY bY + + cZ cZ + d + d 则有：则有： 1 1）若）若 D D＞＞0 0，，P P点处在点处在A A半空间；半空间； 2 2）若）若 D D＝＝0 0，，P P点处在平点处在平面上；面上；   3 3）若）若 D D＜＜0 0，，P P点处在点处在B B半空间。

半空间 图图 4.6 4.6切面与三角面片的交切面与三角面片的交   3 三角面片与剖切面的切割运算    4 表面模型的剖切 　　(1)边表和顶点表均为动态链表结构  (2)表面模型的剖切计算  class class CedgeCedge  class class CedgeVertexCedgeVertex 4.4 剖切截面的生成      4.4.1 边界多边形包含关系检测与确定   1 封闭环的检出              2 封闭轮廓的包含性检测    夹角之和检验法夹角之和检验法: :  4.4.2 剖切面区域的三角剖分    1.任意平面多边形Delaunay三角剖分示意图图4.11 图4.7对应轮廓的三角剖分  4.5 手术开窗操作 开窗操作一般是用立方体或棱柱对重建模型进行切割，模型处于剖切体之内的部分被切割掉，之外的部分被保留下来 4.6 实验结果分析 图图4.12　模型的剖切与开窗　模型的剖切与开窗 　　　　(a)剖切　　　　　　　　　　剖切　　　　　　　　　　 (b)开窗　　　　　开窗　　　　　                (c)开窗开窗 5 由基于轮廓重建的表面模型构建由基于轮廓重建的表面模型构建实体几何模型实体几何模型 5.1 引言 5.2 相关工作        （1）提取边界轮廓线  （2）提取轮廓线上的特征点  （3）轮廓对应  （4）三维表面重建  （（a a））轮廓线轮廓线 （（b b））基础轮廓表面模型基础轮廓表面模型 （（c c））左分支表面模型左分支表面模型 （（d d））右分支表面模型右分支表面模型  图图5.3 5.3 轮廓及表面子模型轮廓及表面子模型 5.3 实体几何模型的构建 　　   1 边界模型的数据结构 图图5.4 5.4 系统系统B- repB- rep模型的数据结构及半边的结构示意模型的数据结构及半边的结构示意 　   2 实体造型的基本操作  欧拉特征关系： 其中v、e、f、s、h分别代表顶点、边、小面、壳和孔。

基本的欧拉操作包括如下互逆的5对：MVFS，MEV，MEF，MEKR，KFMRH；KVFS，KEV，KEF，KEMR，MFKRH其中M表示构造，K表示删除，S、E、V、F、R、H分别表示体、边、顶点、面、环、孔 　   3 由轮廓重建的表面模型重建实体几何模型的方法    用图5.10(a)的表面模型说明构建实体的主要步骤：  ( (a)a) ( (b)b) （（b b））光照图光照图 图图5.11 5.11 实体模型实体模型 （（a a））线框图线框图 实体几何模型的构建结果6 6 适用于适形放射治疗规划的医学图适用于适形放射治疗规划的医学图象三维重建系统的开发象三维重建系统的开发 图象输入图象输入 二维图象预处理　二维图象预处理　 图象分割与提取图象分割与提取 三维重建三维重建 （（CT/MRICT/MRI）　）　 （（滤波、插值）滤波、插值） （自动分割、手工勾画）（自动分割、手工勾画）       （（MCMC、、MTMT表面重建）表面重建） 病病变变体投影体投影轮轮廓　　照射射束廓　　照射射束设设置置 效果显示　　效果显示　　 几何操作几何操作（（由由此此计计算算光光栅栅廓廓线线）） （（放放射射治治疗疗规规划划）） （（ 颜颜色色、、半半透透明明）） （（剖剖切切、、手手术术开开窗窗）） 　　 图　图　6.16.1三维重建过程示意图三维重建过程示意图l 重建系统结构重建系统结构 图图6.26.2　系统结构　系统结构 l 系统程序流程 图图 6.3 6.3 系统程序流程系统程序流程 l 系统数据结构 图图6.4 6.4 系统数据结构系统数据结构  系统界面系统界面  治疗射束安排与光栅轮廓线计算 7 7 结论与展望结论与展望           7.1 工作总结工作总结(1) 对输入图象进行了滤波、断层插值并封装成规则体数据。

定义了体数据的内存记录方式及外存文件格式，压缩存储空间(2) 提出并实现了三维医学图象交互分割的方法，交互分割的技术路线是：先分析断层图象，交互给定分割阈值，对图象二值化，然后选择适当的形态学操作进行区域修整，最后用种子填充的方法填充出所要区域(3)实现了MC算法和MT算法构造表面模型针对MT算法，为避免体元棱边与等值面交点的重复性插值计算，提出了相关性处理方法采用相关性处理，加快了MT算法的重建速度(4) 实现了Hoppe的边收缩算法，重建模型经简化90％，依然能较好地保持原模型特征，基本不影响视觉效果模型经网格简化，绘制时间大大缩短，提高了交互时绘制能力(5)提出了对重建模型实施剖切及手术开窗的一种方法(6)提出了基于轮廓重建的表面模型构建实体几何模型的方法实现步骤可作为造型系统的一种造型方式 (7)开发了一个适用于多叶光栅适形调强放射治疗的医学图象三维重建系统提出并建立了系统数据结构，提出了自动分割与手工勾画轮廓的方法对分割出的组织重建其三维几何模型，并对模型实现了网格简化实现了治疗射束的设置及多叶光栅轮廓的计算 7.2 展望展望 （1）基于知识模型的医学图象的分割 （2）多模态图象的匹配和融合 （3）表面绘制与直接体绘制的结合 （4）虚拟手术工作环境 　　　　谢谢大家光临！谢谢大家光临！　　　。