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

1、资料仅供参考,不当之处，请联系改正。1.1 引言引言 1.2 基于三维数据的建模与可视化基于三维数据的建模与可视化 1.三维数据的来源与分类 2.三维数据建模及可视化研究内容 数据预处理 建模 绘制与显示 资料仅供参考,不当之处，请联系改正。3.数据建模技术综述 三维几何模型 基于三维数据的建模方法 1)基于断层轮廓的表面重建 2)基于体素的等值面重建 3)几何变形模型 4)体素建模 资料仅供参考,不当之处，请联系改正。1.3 1.3 医学图象三维重建技术综述医学图象三维重建技术综述 1 医学图象的预处理 2 医学图象的分割 3 三维重建方法 4 模型的网格简化 资料仅供参考,不当之处，请联系

2、改正。1.4 1.4 医学图象三维重建在医疗中的应用医学图象三维重建在医疗中的应用 1 在医疗诊断中的应用 2 在手术规划及放射治疗规划中的应用 3 在整形与假肢外科中的应用 4 在虚拟手术及解剖教育中的应用资料仅供参考,不当之处，请联系改正。1.5 1.5论文背景及主要工作论文背景及主要工作 1 论文背景及研究意义 2 本文的主要工作 1)图象预处理，组织器官分割与提取 2)MC、MT算法构建表面几何模型 3)模型表面网格简化，剖切与开窗 4)由表面几何模型转换成实体几何模型 5)适用于适形放射治疗规划的医学图象 三维重建系统的开发资料仅供参考,不当之处，请联系改正。医学图象预处理医学图象预

3、处理 分割流程分割流程 图图2.1 2.1 三维医学图象分割流程三维医学图象分割流程 资料仅供参考,不当之处，请联系改正。2.2 CT 2.2 CT、MRIMRI图象的获取与输入图象的获取与输入 2.2 2.2 二维图象处理与规则体数据封装二维图象处理与规则体数据封装 1 二维图象滤波 (1)邻域平均法 (2)中值滤波法 (3)保持边缘滤波法 资料仅供参考,不当之处，请联系改正。2 断层图象间插值 3 三维规则体数据封装 (1)内存记录方式 (2)体数据文件格式 图图2.4 2.4 体数据内存记录方式体数据内存记录方式 资料仅供参考,不当之处，请联系改正。2.3 交互分割过程交互分割过程 1

4、三维图象二值化 其余若0),(1),(21qzyxfqzyxf资料仅供参考,不当之处，请联系改正。二值化结果 图图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

5、闭合操作闭合操作 资料仅供参考,不当之处，请联系改正。3 种子填充法进行组织提取 图图2.12 2.12 交互分割结果交互分割结果 图图2.13 2.13 对分割区域的重建对分割区域的重建 资料仅供参考,不当之处，请联系改正。3.2.1 体素模型(a)(a)方向无关的三线性插值模型方向无关的三线性插值模型 (b)(b)方向有关的三线性插值模型方向有关的三线性插值模型图图3.1 3.1 体素模型体素模型 资料仅供参考,不当之处，请联系改正。3.2.2 等值面（IsoSurface）定义 1三线性插值结果 2等值面定义xyzazxayzaxyazayaxaazyxf76543210),(是常数cc

6、zyxfzyx,),(|),(等值面是三次曲面资料仅供参考,不当之处，请联系改正。3.3 移动立方体（Marching Cubes）算 法抽取等值面 1 MC1 MC算法的基本原理算法的基本原理 a a 体素中等值面剖分方式的确定体素中等值面剖分方式的确定 1 1）如立方体顶点的数据值如立方体顶点的数据值等值面的值，则定义该顶点等值面的值，则定义该顶点位于等值面之外，记为位于等值面之外，记为“0”“0”；2 2）如立方体顶点的数据值等值面的值，则定义该顶点如立方体顶点的数据值等值面的值，则定义该顶点位于等值之内，记为位于等值之内，记为“1”“1”。8 8个顶点，每个顶点共有个顶点，每个顶点共有

7、2 2个状态，因此共个状态，因此共256256种组合状态种组合状态 根据互补对称性根据互补对称性,256 128 256 128 根据旋转对称性，根据旋转对称性，256 15 256 15 (1)体素中由三角片逼近的等值面计算 (2)三角片各顶点法向量计算资料仅供参考,不当之处，请联系改正。资料仅供参考,不当之处，请联系改正。2 等值面连接方式上的二义性(a)(a)连接方式二义性的二维表示连接方式二义性的二维表示 (b)(b)连接方式二义性的三维表示连接方式二义性的三维表示 图图3.5 3.5 拓扑不一致造成孔隙拓扑不一致造成孔隙 图图3.4 MC3.4 MC方法的二义性方法的二义性 资料仅供

8、参考,不当之处，请联系改正。3.3.3 渐近线判别法消除二义性资料仅供参考,不当之处，请联系改正。3.3.4 MC算法的重建结果及分析(a)(a)256256256256109109MRIMRI表皮重建表皮重建（b b）12812812812893CT93CT颅骨重建颅骨重建（c c）12812812812893CT93CT表皮重建表皮重建 三角面片：三角面片：696889顶点：顶点：347322三角面片：三角面片：187559顶点：顶点：94015三角面片：三角面片：137799顶点：顶点：69331 图图3.83.8MCMC算法重建的表面模型算法重建的表面模型 资料仅供参考,不当之处，请联

9、系改正。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相邻立方体公共面上的剖分一致性相邻立方体公共面上的剖分一致性 资料仅供参考,不当之处，请

10、联系改正。3.4.3 相关性处理加速MT重建速度 1 体素内的相关性处理 2 体素间的相关性处理 图图3.14 3.14 剖分后立方体的顶点及棱边编号剖分后立方体的顶点及棱边编号 资料仅供参考,不当之处，请联系改正。3.4.4 MT算法的重建结果及分析(a)128(a)128128128113CT113CT颅骨重建颅骨重建（b b）104104185185220CT220CT脚骨骼重建脚骨骼重建(c)128(c)128128128113CT113CT表皮重建表皮重建 三角面片：三角面片：423998顶点：顶点：211905三角面片：三角面片：365858顶点：顶点：183056三角面片：三角面

11、片：331290顶点：顶点：165808图图3.153.15MTMT算法重建的表面模型算法重建的表面模型资料仅供参考,不当之处，请联系改正。4.2 基于边收缩的网格简化算法 1 网格简化算法简述 （1 1）抽样（）抽样（SamplingSampling）（2 2）自适应细分）自适应细分(Adaptive subdivision)(Adaptive subdivision)（3 3）删除）删除(Decimation)(Decimation)（4 4）顶点合并）顶点合并 （Vertex mergingVertex merging）资料仅供参考,不当之处，请联系改正。2 基于边收缩的网格简化算法Ho

12、ppe采用显式能量函数E（M）来度量简化网格与原始网格的逼近度Hoppe96：其中Edist(M)为M的距离能量，它定义为点集到网格的距离平方：Espring(M)为弹性能量，这相当于在的每条边上均放置一条弹性系数为k的弹簧，即：Escalar(M)度量M的标量属性的精度，而Edisc(M)则度量了M上视觉不连续的特征线（如边界线、侧影轮廓线等）的几何精度。)()()()()(MEMEMEMEMEdiscscalarspringdistnividistKxdME12|)(|,()(KjijispringvvkME,2|)()()(MEMEE资料仅供参考,不当之处，请联系改正。边收缩过程示意图边

13、收缩过程示意图 (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

14、,Y,Z）,定义定义“距离距离”D”D：D=aX+bY+cZ+d D=aX+bY+cZ+d 则有：则有：1 1）若）若 D D0 0，P P点处在点处在A A半空间；半空间；2 2）若）若 D D0 0，P P点处在平点处在平面上；面上；3 3）若）若 D D0 0，P P点处在点处在B B半空间。半空间。图图 4.6 4.6切面与三角面片的交切面与三角面片的交 资料仅供参考,不当之处，请联系改正。3 三角面片与剖切面的切割运算 4 表面模型的剖切(1)边表和顶点表均为动态链表结构 (2)表面模型的剖切计算 class Cedgeclass Cedge class CedgeVertexcla

15、ss CedgeVertex资料仅供参考,不当之处，请联系改正。4.4 剖切截面的生成 4.4.1 边界多边形包含关系检测与确定 1 封闭环的检出 2 封闭轮廓的包含性检测 资料仅供参考,不当之处，请联系改正。夹角之和检验法夹角之和检验法:0n1iinii12资料仅供参考,不当之处，请联系改正。4.4.2 剖切面区域的三角剖分 1.任意平面多边形Delaunay三角剖分示意图图4.11 图4.7对应轮廓的三角剖分 资料仅供参考,不当之处，请联系改正。4.5 手术开窗操作 开窗操作一般是用立方体或棱柱对重建模型进行切割，模型处于剖切体之内的部分被切割掉，之外的部分被保留下来。资料仅供参考,不当之

16、处，请联系改正。4.6 实验结果分析 图图4.12模型的剖切与开窗模型的剖切与开窗(a)剖切剖切(b)开窗开窗 (c)开窗开窗资料仅供参考,不当之处，请联系改正。5.1 引言 5.2 相关工作 （1）提取边界轮廓线 （2）提取轮廓线上的特征点 （3）轮廓对应 （4）三维表面重建 资料仅供参考,不当之处，请联系改正。（a a）轮廓线）轮廓线 （b b）基础轮廓表面模型）基础轮廓表面模型（c c）左分支表面模型）左分支表面模型 （d d）右分支表面模型）右分支表面模型 图图5.3 5.3 轮廓及表面子模型轮廓及表面子模型 资料仅供参考,不当之处，请联系改正。5.3 实体几何模型的构建 1 边界模型

17、的数据结构 图图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分别表示体、边、顶点、面、环、孔。)(2hsfevrhsfev)(2资料仅供参考,不当之处，请联系改正。3 由轮廓重建的表面模型重建实体几何模型的方法 用图5.10(a)

18、的表面模型说明构建实体的主要步骤：(a)(a)(b)(b)资料仅供参考,不当之处，请联系改正。（b b）光照图）光照图 图图5.11 5.11 实体模型实体模型（a a）线框图）线框图 实体几何模型的构建结果资料仅供参考,不当之处，请联系改正。图象输入图象输入 二维图象预处理二维图象预处理 图象分割与提取图象分割与提取 三维重建三维重建 （CT/MRICT/MRI）（滤波、插值）（滤波、插值）（自动分割、手工勾画）（自动分割、手工勾画）（MCMC、MTMT表面重建）表面重建）病变体投影轮廓照射射束设置病变体投影轮廓照射射束设置 效果显示效果显示 几何操作几何操作（由此计算光栅廓线）（由此计算光

19、栅廓线）（放射治疗规划）（放射治疗规划）（颜色、半透明）颜色、半透明）（剖切、手术开窗）（剖切、手术开窗）图图6.16.1三维重建过程示意图三维重建过程示意图资料仅供参考,不当之处，请联系改正。l 重建系统结构重建系统结构 图图6.26.2系统结构系统结构 资料仅供参考,不当之处，请联系改正。l 系统程序流程 图图 6.3 6.3 系统程序流程系统程序流程 资料仅供参考,不当之处，请联系改正。l 系统数据结构 图图6.4 6.4 系统数据结构系统数据结构 资料仅供参考,不当之处，请联系改正。系统界面系统界面 资料仅供参考,不当之处，请联系改正。治疗射束安排与光栅轮廓线计算 资料仅供参考,不当之

20、处，请联系改正。7.1 工作总结工作总结(1)对输入图象进行了滤波、断层插值并封装成规则体数据。定义了体数据的内存记录方式及外存文件格式，压缩存储空间。(2)提出并实现了三维医学图象交互分割的方法，交互分割的技术路线是：先分析断层图象，交互给定分割阈值，对图象二值化，然后选择适当的形态学操作进行区域修整，最后用种子填充的方法填充出所要区域。(3)实现了MC算法和MT算法构造表面模型。针对MT算法，为避免体元棱边与等值面交点的重复性插值计算，提出了相关性处理方法。采用相关性处理，加快了MT算法的重建速度。资料仅供参考,不当之处，请联系改正。(4)实现了Hoppe的边收缩算法，重建模型经简化90，

21、依然能较好地保持原模型特征，基本不影响视觉效果。模型经网格简化，绘制时间大大缩短，提高了交互时绘制能力。(5)提出了对重建模型实施剖切及手术开窗的一种方法。(6)提出了基于轮廓重建的表面模型构建实体几何模型的方法实现步骤。可作为造型系统的一种造型方式。(7)开发了一个适用于多叶光栅适形调强放射治疗的医学图象三维重建系统。提出并建立了系统数据结构，提出了自动分割与手工勾画轮廓的方法。对分割出的组织重建其三维几何模型，并对模型实现了网格简化。实现了治疗射束的设置及多叶光栅轮廓的计算。资料仅供参考,不当之处，请联系改正。7.2 展望展望（1）基于知识模型的医学图象的分割（2）多模态图象的匹配和融合（3）表面绘制与直接体绘制的结合（4）虚拟手术工作环境资料仅供参考,不当之处，请联系改正。谢谢大家光临！谢谢大家光临！