表格模板-医学图象三维表面重建技术研究博士生秦绪佳导师欧宗瑛教授 精品.ppt

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.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 三维图象二值化 其余若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闭合操作闭合操作 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),(是常数cczyxfzyx,),(| ),(等值面是三次曲面 3.3 移动立方体（Marching Cubes）算 法抽取等值面 1 MC1 MC算法的基本原理算法的基本原理 a a 体素中等值面剖分方式的确定体素中等值面剖分方式的确定 1 1） 如立方体顶点的数据值如立方体顶点的数据值等值面的值，则定义该顶点等值面的值，则定义该顶点位于等值面之外，记为位于等值面之外，记为“0”0”； 2 2） 如立方体顶点的数据值等值面的值，则定义该顶点如立方体顶点的数据值等值面的值，则定义该顶点位于等值之内，记为位于等值之内，记为“1”1”。 8 8个顶点，每个顶点共有个顶点，每个顶点共有2 2个状态，因此共个状态，因此共256256种组合状态种组合状态 根据互补对称性根据互补对称性, 256 128, 256 128 根据旋转对称性，根据旋转对称性， 256 15256 15 (1)体素中由三角片逼近的等值面计算 (2)三角片各顶点法向量计算 2 等值面连接方式上的二义性 (a)(a)连接方式二义性的二维表示连接方式二义性的二维表示 (b) (b) 连接方式二义性的三维表示连接方式二义性的三维表示 图图3.5 3.5 拓扑不一致造成孔隙拓扑不一致造成孔隙 图图3.4 MC3.4 MC方法的二义性方法的二义性 3.3.3 渐近线判别法消除二义性 3.3.4 MC算法的重建结果及分析 (a)(a) 256256256256109MRI109MRI表皮重建表皮重建 （b b）12812812812893CT93CT颅骨重建颅骨重建 （c c）12812812812893CT93CT表皮重建表皮重建 三角面片：三角面片：696889顶点：顶点：347322三角面片：三角面片：187559顶点：顶点：94015三角面片：三角面片：137799顶点：顶点：69331 图图3.83.8MCMC算法重建的表面模型算法重建的表面模型 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(a) 128128128113CT113CT颅骨重建颅骨重建 （b b）104104185185220CT220CT脚骨骼重建脚骨骼重建 (c) 128(c) 128128128113CT113CT表皮重建表皮重建 三角面片：三角面片：423998顶点：顶点：211905三角面片：三角面片：365858顶点：顶点：183056三角面片：三角面片：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 基于边收缩的网格简化算法Hoppe采用显式能量函数E（M）来度量简化网格与原始网格的逼近度Hoppe96：其中Edist(M)为M的距离能量，它定义为点集到网格的距离平方：Espring(M)为弹性能量，这相当于在的每条边上均放置一条弹性系数为k的弹簧，即： Escalar(M)度量M的标量属性的精度，而Edisc(M)则度量了M上视觉不连续的特征线（如边界线、侧影轮廓线等）的几何精度。 )()()()()(MEMEMEMEMEdiscscalarspringdistnividistKxdME12|)(|,()(KjijispringvvkME,2|)()()(MEMEE 边收缩过程示意图边收缩过程示意图 (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 = aX + bY + cZD = aX + bY + cZ + d + 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.64.6切面与三角面片的交切面与三角面片的交 3 三角面片与剖切面的切割运算 4 表面模型的剖切 (1)边表和顶点表均为动态链表结构 (2)表面模型的剖切计算 class Cedgeclass Cedge class CedgeVertexclass CedgeVertex 4.4 剖切截面的生成 4.4.1 边界多边形包含关系检测与确定 1 封闭环的检出 2 封闭轮廓的包含性检测 夹角之和检验法夹角之和检验法: : 0n1iinii12 4.4.2 剖切面区域的三角剖分 1.任意平面多边形Delaunay三角剖分示意图图4.11 图4.7对应轮廓的三角剖分 4.5 手术开窗操作 开窗操作一般是用立方体或棱柱对重建模型进行切割，模型处于剖切体之内的部分被切割掉，之外的部分被保留下来。 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 边界模型的数据结构 图图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)的表面模型说明构建实体的主要步骤： ( (a)a) ( (b)b) （b b）光照图）光照图 图图5.11 5.11 实体模型实体模型 （a a）线框图）线框图 实体几何模型的构建结果 图象输入图象输入 二维图象预处理二维图象预处理 图象分割与提取图象分割与提取 三维重建三维重建 （CT/MRICT/MRI） （滤波、插值）（滤波、插值） （自动分割、手工勾画）（自动分割、手工勾画） （MCMC、MTMT表面重建）表面重建） 病变体投影轮廓照射射束设置病变体投影轮廓照射射束设置 效果显示效果显示 几何操作几何操作（由此计算光栅廓线）（由此计算光栅廓线） （放射治疗规划）（放射治疗规划） （ 颜色、半透明）颜色、半透明） （剖切、手术开窗）（剖切、手术开窗） 图图6.16.1三维重建过程示意图三维重建过程示意图l 重建系统结构重建系统结构 图图6.26.2系统结构系统结构 l 系统程序流程 图图 6.3 6.3 系统程序流程系统程序流程 l 系统数据结构 图图6.4 6.4 系统数据结构系统数据结构 系统界面系统界面 治疗射束安排与光栅轮廓线计算 7.1 工作总结工作总结(1) 对输入图象进行了滤波、断层插值并封装成规则体数据。定义了体数据