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

1、医学图象三维重建及可视化技术研究医学图象三维重建及可视化技术研究1 绪论 1.1 引言 1.2 基于三维数据的建模与可视化 1.三维数据的来源与分类 2.三维数据建模及可视化研究内容 数据预处理 建模 绘制与显示 医学图象三维重建及可视化技术研究21 绪论医学图象三维重建及可视化技术研究2 3.数据建模技术综述 三维几何模型 基于三维数据的建模方法 1)基于断层轮廓的表面重建 2)基于体素的等值面重建 3)几何变形模型 4)体素建模 医学图象三维重建及可视化技术研究3 3.数据建模技术综述 医学图象三维重建及可视化技术研 1.3 医学图象三维重建技术综述 1 医学图象的预处理 2 医学图象的分

2、割 3 三维重建方法 4 模型的网格简化 医学图象三维重建及可视化技术研究4 1.3 医学图象三维重建技术综述 医学图象三维重建及 1.4 医学图象三维重建在医疗中的应用 1 在医疗诊断中的应用 2 在手术规划及放射治疗规划中的应用 3 在整形与假肢外科中的应用 4 在虚拟手术及解剖教育中的应用医学图象三维重建及可视化技术研究5 1.4 医学图象三维重建在医疗中的应用 医学图象三维重 1.5论文背景及主要工作 1 论文背景及研究意义 2 本文的主要工作 1) 图象预处理，组织器官分割与提取 2) MC、MT算法构建表面几何模型 3) 模型表面网格简化，剖切与开窗 4) 由表面几何模型转换成实体

3、几何模型 5) 适用于适形放射治疗规划的医学图象 三维重建系统的开发医学图象三维重建及可视化技术研究6 1.5论文背景及主要工作 医学图象三维重建及可视化技术 2医学图象预处理与人体组织的分割医学图象预处理 分割流程 图2.1 三维医学图象分割流程 医学图象三维重建及可视化技术研究7 2医学图象预处理与人体组织的分割医学图象预 2.2 CT、MRI图象的获取与输入 2.2 二维图象处理与规则体数据封装 1 二维图象滤波 (1) 邻域平均法 (2) 中值滤波法 (3) 保持边缘滤波法 医学图象三维重建及可视化技术研究8 2.2 CT、MRI图象的获取与输入医学图象三维重建及 2 断层图象间插值

4、3 三维规则体数据封装 (1) 内存记录方式 (2) 体数据文件格式 图2.4 体数据内存记录方式 医学图象三维重建及可视化技术研究9 2 断层图象间插值 图2.4 体数据内存记录方式 2.3 交互分割过程 1 三维图象二值化 医学图象三维重建及可视化技术研究10 2.3 交互分割过程 医学图象三维重建及可视化技术研究 二值化结果 图2.8 断层图象二值化结果医学图象三维重建及可视化技术研究11 二值化结果 图2.8 断层图象 2 数学形态学操作进行区域修整 (1) 二值形态学操作简述 (a) 原图象 (b)结构元素 (c)对原图象的腐蚀 (d)对腐蚀图象的膨胀 图2.9 开启操作 (a)原图

5、象 (b)结构元素 (c)对原图象的膨胀 (d)对膨胀图象的腐蚀 图2.10闭合操作 医学图象三维重建及可视化技术研究12 2 数学形态学操作进行区域修整 (a) 原图象 3 种子填充法进行组织提取 图2.12 交互分割结果 图2.13 对分割区域的重建 医学图象三维重建及可视化技术研究13 3 种子填充法进行组织提取 图2.12 交互3 基于规则体数据的三维表面模型的构建 3.2.1 体素模型 (a) 方向无关的三线性插值模型 (b) 方向有关的三线性插值模型图3.1 体素模型 医学图象三维重建及可视化技术研究143 基于规则体数据的三维表面模型的构建 3.2.1 体素 3.2.2 等值面（

6、IsoSurface）定义 1三线性插值结果 2等值面定义等值面是三次曲面医学图象三维重建及可视化技术研究15 3.2.2 等值面（IsoSurface）定义 等 3.3 移动立方体（Marching Cubes）算 法抽取等值面 1 MC算法的基本原理 a 体素中等值面剖分方式的确定 1） 如立方体顶点的数据值等值面的值，则定义该顶点位于等值面之外，记为“0”； 2） 如立方体顶点的数据值等值面的值，则定义该顶点位于等值之内，记为“1”。 8个顶点，每个顶点共有2个状态，因此共256种组合状态 根据互补对称性, 256 128 根据旋转对称性， 256 15 (1)体素中由三角片逼近的等值面

7、计算 (2)三角片各顶点法向量计算医学图象三维重建及可视化技术研究16 3.3 移动立方体（Marching Cubes）算 医学图象三维重建及可视化技术研究17 医学图象三维重建及可视化技术研究17 2 等值面连接方式上的二义性 (a)连接方式二义性的二维表示 (b) 连接方式二义性的三维表示 图3.5 拓扑不一致造成孔隙 图3.4 MC方法的二义性 医学图象三维重建及可视化技术研究18 2 等值面连接方式上的二义性 (a)连接方式二义性的二 3.3.3 渐近线判别法消除二义性医学图象三维重建及可视化技术研究19 3.3.3 渐近线判别法消除二义性医学图象三维重建及可 3.3.4 MC算法的

8、重建结果及分析 256256109MRI表皮重建 （b）12812893CT颅骨重建 （c）12812893CT表皮重建 三角面片：696889顶点：347322三角面片：187559顶点：94015三角面片：137799顶点：69331 图3.8MC算法重建的表面模型 医学图象三维重建及可视化技术研究20 3.3.4 MC算法的重建结果及分析 256256 3.4 移动四面体（Marching Tetrahedra）算法抽取等值面 3.4.1 MT算法的基本原理 图3.9 立方体的四面体剖分 图3.10 四面体中的等值面 医学图象三维重建及可视化技术研究21 3.4 移动四面体（Marchi

9、ng 图3.9 立方 3.4.2 四面体剖分的一致性处理 图3.11 立方体剖分为四面体的不同方式 图3.12 两相邻立方体剖分不一致时共有面的剖分情况 图3.13相邻立方体公共面上的剖分一致性 医学图象三维重建及可视化技术研究22 3.4.2 四面体剖分的一致性处理 图3.11 立方体 3.4.3 相关性处理加速MT重建速度 1 体素内的相关性处理 2 体素间的相关性处理 图3.14 剖分后立方体的顶点及棱边编号 医学图象三维重建及可视化技术研究23 3.4.3 相关性处理加速MT重建速度 图3.14 剖 3.4.4 MT算法的重建结果及分析 (a) 128128113CT颅骨重建 （b）1

10、04185220CT脚骨骼重建 (c) 128128113CT表皮重建 三角面片：423998顶点：211905三角面片：365858顶点：183056三角面片：331290顶点：165808图3.15MT算法重建的表面模型医学图象三维重建及可视化技术研究24 3.4.4 MT算法的重建结果及分析 (a) 1284.三维模型的网格简化与模型的剖切 4.2 基于边收缩的网格简化算法 1 网格简化算法简述 （1）抽样（Sampling） （2）自适应细分(Adaptive subdivision) （3）删除 (Decimation) （4）顶点合并 （Vertex merging） 医学图象三维

11、重建及可视化技术研究254.三维模型的网格简化与模型的剖切 （1）抽样（S 2 基于边收缩的网格简化算法Hoppe采用显式能量函数E（M）来度量简化网格与原始网格的逼近度Hoppe96：其中Edist(M)为M的距离能量，它定义为点集到网格的距离平方：Espring(M)为弹性能量，这相当于在的每条边上均放置一条弹性系数为k的弹簧，即： Escalar(M)度量M的标量属性的精度，而Edisc(M)则度量了M上视觉不连续的特征线（如边界线、侧影轮廓线等）的几何精度。 医学图象三维重建及可视化技术研究26 2 基于边收缩的网格简化算法医学图象三维重建及可视化技术 边收缩过程示意图 (a)收缩前

12、(b)收缩后 图4.1边收缩过程医学图象三维重建及可视化技术研究27 边收缩过程示意图 (a)收缩前 4 网格简化结果 MT重建结果 简化50 简化90简化90表面绘制医学图象三维重建及可视化技术研究28 4 网格简化结果 MT重建结果 MC重建结果简化50简化85 简化85表面绘制医学图象三维重建及可视化技术研究29 MC重建结果简化50简化854.3 三维模型的剖切 4.3.1 模型三角面片的剖切处理 1 平面方程的确定 ax + by + cz + d = 0 医学图象三维重建及可视化技术研究304.3 三维模型的剖切 医学图象三维重建及可视化技术研究30 2 三角面片与剖切平面的求交检

13、测 定义空间一点P（X,Y,Z）, 定义“距离”D： D = aX + bY + cZ + d 则有： 1）若 D0，P点处在A半空间； 2）若 D0，P点处在平面上； 3）若 D0，P点处在B半空间。 图 4.6切面与三角面片的交 医学图象三维重建及可视化技术研究31 2 三角面片与剖切平面的求交检测 定义空间一点P（X, 3 三角面片与剖切面的切割运算 4 表面模型的剖切 (1)边表和顶点表均为动态链表结构 (2)表面模型的剖切计算 class Cedge class CedgeVertex医学图象三维重建及可视化技术研究32 3 三角面片与剖切面的切割运算 class Ce 4.4 剖切

14、截面的生成 4.4.1 边界多边形包含关系检测与确定 1 封闭环的检出 2 封闭轮廓的包含性检测 医学图象三维重建及可视化技术研究33 4.4 剖切截面的生成医学图象三维重建及可视化技术研究 夹角之和检验法: 医学图象三维重建及可视化技术研究34 夹角之和检验法: 医学图象三维重建及可视化技术研究 4.4.2 剖切面区域的三角剖分 1.任意平面多边形Delaunay三角剖分示意图图4.11 图4.7对应轮廓的三角剖分 医学图象三维重建及可视化技术研究35 4.4.2 剖切面区域的三角剖分 图4.11 图4医学图象三维重建及可视化技术研究培训课件 4.6 实验结果分析 图4.12模型的剖切与开窗

15、 (a)剖切 (b)开窗 (c)开窗医学图象三维重建及可视化技术研究37 4.6 实验结果分析 图4.12模型的剖切 5 由基于轮廓重建的表面模型构建实体几何模型 5.1 引言 5.2 相关工作 （1）提取边界轮廓线 （2）提取轮廓线上的特征点 （3）轮廓对应 （4）三维表面重建 医学图象三维重建及可视化技术研究38 5 由基于轮廓重建的表面模型构建实体几何模型 5. （a）轮廓线 （b）基础轮廓表面模型 （c）左分支表面模型 （d）右分支表面模型 图5.3 轮廓及表面子模型 医学图象三维重建及可视化技术研究39 （a）轮廓线 （b）基础轮廓表面模型 （c）左分支表面模型5.3 实体几何模型的

16、构建 1 边界模型的数据结构 图5.4 系统B- rep模型的数据结构及半边的结构示意 医学图象三维重建及可视化技术研究405.3 实体几何模型的构建 图5.4 系统B- 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分别表示体、边、顶点、面、环、孔。 医学图象三维重建及可视化技术研究41 2 实体造型的基本操作 医学图象三维重建及可视化技术 3 由轮廓重建

17、的表面模型重建实体几何模型的方法 用图5.10(a)的表面模型说明构建实体的主要步骤： (a) (b) 医学图象三维重建及可视化技术研究42 3 由轮廓重建的表面模型重建实体几何模型的方法 (（b）光照图 图5.11 实体模型 （a）线框图 实体几何模型的构建结果医学图象三维重建及可视化技术研究43（b）光照图 图5.11 实体模型 （a）线框图 实体几何6 适用于适形放射治疗规划的医学图象三维重建系统的开发 图象输入 二维图象预处理 图象分割与提取 三维重建 （CT/MRI） （滤波、插值） （自动分割、手工勾画） （MC、MT表面重建） 病变体投影轮廓照射射束设置 效果显示 几何操作（由此

18、计算光栅廓线） （放射治疗规划） （ 颜色、半透明） （剖切、手术开窗） 图6.1三维重建过程示意图医学图象三维重建及可视化技术研究446 适用于适形放射治疗规划的医学图象三维重建系统的开发 重建系统结构 图6.2系统结构 医学图象三维重建及可视化技术研究45 重建系统结构 图6.2系统结构 医学 系统程序流程 图 6.3 系统程序流程 医学图象三维重建及可视化技术研究46 系统程序流程 图 6.3 系统程序流程 系统数据结构 图6.4 系统数据结构 医学图象三维重建及可视化技术研究47 系统数据结构 图6.4 系统数据结构 医 系统界面 医学图象三维重建及可视化技术研究48 系统界面 医学图象三维重建及 治疗射束安排与光栅轮廓线计算 医学图象三维重建及可视化技术研究49 治疗射束安排与光栅轮廓线计算 医学图象三维重建及可7 结论与展望 7.1 工作总结(1) 对输入图象进行了滤波、断层插值并封装成规则体数据。定义了体数据的内存记录方式及外存文件格式，压缩存储空间。(2) 提出并实现了三维医学图象交互分割的方法，交互分割的技术路线是：先分析断层图象，交互给定分割阈值，对图象二值化，然后选择适当的形态学操作进行区域修整