1、.1主要内容主要内容研究背景研究背景12医学图像生成三维医学图像生成三维CAD模型模型应用验证应用验证43双向平滑滤波和分支表面生成算法双向平滑滤波和分支表面生成算法5结论结论.2研究背景研究背景 如今,医学成像技术,如计算机断层扫描、磁共振、PET和B超等,被广泛用于医疗诊断和外科手术规划中。随着计算机成像和扫描技术的发展,现在已经可以由医学图像(如CT)生成医用三维模型。重建后的三维模型可以提供宝贵的医疗信息,可以作为一个强大的诊断工具使外科医生能够了解人体复杂的解剖学结构,还可以用于制造人造器官(人造骨,血管)和各种医学仿真和分析。.3医学医学图像图像生成三维生成三维CADCAD模型模型
2、 对医学图像的三维重建,现在存在着很多方法,但大多数可以归为两类:医学图像生成三维模型基于体素的方法基于轮廓线的方法.4基于体基于体素的素的方法方法优点优点:运用了marching cube 算法,用三角形体素网面来构成高分辨率的三维模型的表面,在数据结构,数据计算处理方面效率高。缺点:缺点:这种运用三角形体素的方法经常产生错误的形状,而且不经过后期处理,这种方法产生的模型很难应用到生产和工程中。基于体素基于体素的方法的方法.5基于轮廓基于轮廓线的线的方法方法 1.从医学图像中提取出横截面轮廓从医学图像中提取出横截面轮廓数据数据2.对横截面轮廓线进行平滑处理对横截面轮廓线进行平滑处理 3.由各
3、个横截面的轮廓线生成三维由各个横截面的轮廓线生成三维模型表面模型表面.6基于轮廓线的方法基于轮廓线的方法轮廓数据轮廓数据提取提取 医学图像可以看作一个三维点阵,扫描方向沿着轴线方向,则每个截面就是一个二维矩阵。这就可以通过人工构描法或者采用图像边缘检测算法,提取出每个截面的轮廓线数据。运用形态学操作和轮廓跟踪算法处理,生成封闭无缺口的轮廓线。完成图像的轮廓线提取后,医学图像就转变为一系列排列紧凑的封闭轮廓线。.7基于轮廓线的方法基于轮廓线的方法轮廓轮廓线的线的平滑处理平滑处理 一旦有了轮廓数据,就可通过连接相邻轮廓线上的数据点创建表面模型。但医学图像的人为因素较多,有时相邻轮廓的形状差异变化很
4、大,这导致拟合后形成的模型表面粗糙,多褶皱。这就需要用平滑算法处理轮廓数据,减缓轮廓数据中的突变,获得较为平滑的轮廓线,进而通过蒙皮算法生成较为平滑的模型表面。.8基于轮廓线的方法基于轮廓线的方法已获取轮廓线已获取轮廓线双向平滑滤波双向平滑滤波分支表面生成分支表面生成平滑的轮廓线组平滑的轮廓线组.9双向双向平滑平滑滤波滤波 轮廓线的平滑处理的第一步,目的是为了在双方向上(平行于横截面方向和垂直于横截面方向)对轮廓线组进行平滑滤波,消去噪声,并确定相邻轮廓线之间的连通性,方便模型表面的生成,避免模型表面的扭曲和褶皱。u方向的轮方向的轮廓线重构廓线重构预设相邻预设相邻轮廓线的轮廓线的连通性连通性垂
5、直方向的垂直方向的轮廓线平轮廓线平滑滤波滑滤波.10u u方向的方向的轮廓轮廓线重构线重构 u方向的轮廓线重构是基于弦长,对单条轮廓线轨迹进行平滑处理的步骤。首先要确定取样比(sampling ratio),然后计算所需弦长。弦长公式:.11u u方向的轮廓线重构方向的轮廓线重构轮廓线的重建:v一条轮廓线:Ck=P1,P2,.,Pm,共m个点;v根据弦长与实际两点之间的距离的比较分为两种情况,通过用新的整形点来取代原有点来实现轮廓线的重建。.12u u方向的轮廓线重构方向的轮廓线重构新的整形点Pnew点在Ci中插入到Pi和Pi+1之间.13u u方向的轮廓线重构方向的轮廓线重构.14u u方向
6、的轮廓线重构方向的轮廓线重构 s是符合下式的最大整数:弦长余量:得到的新整形点Pnew点在Ci中取代被覆盖的Pi+1,.,Pi+s+1.15u u方向的轮廓线重构方向的轮廓线重构 传统的高斯滤波和均值滤波等滤波器同样可以实现对轮廓线的平滑处理,但是同时随着算法迭代次数的增加会引起轮廓线的过度收缩和对轮廓线形状的篡改。相比之下这种轮廓线重构的方法不仅可以消除轮廓线上形状的突变和噪点,还能够防止轮廓线的过度紧缩,防止丢失轮廓线上的形状细节。.16u u方向的轮廓线重构方向的轮廓线重构a.CTa.CT图像上的最初轮廓线图像上的最初轮廓线b.b.高斯滤波处理高斯滤波处理c.c.均值滤波处理均值滤波处理
7、d.d.轮廓线重构处理轮廓线重构处理.17u u方向的轮廓线重构方向的轮廓线重构 横坐标:迭代次数 纵坐标:轮廓区域面积轮廓线轮廓线重构重构高斯滤高斯滤波波均值均值滤波滤波.18相邻轮廓的连通性相邻轮廓的连通性 u方向的轮廓线重构结束后,相邻轮廓之间的连接决定了垂直方向上的形状,因此在垂直方向平滑滤波前需要预设相邻轮廓线间的连通性,确定相邻轮廓线上采样点的一一对应的连接关系。整个预设过程是通过将轮廓线组进行参数化来实现的。主要步骤如下主要步骤如下:轮廓组的轮廓组的参数化参数化确定各个点的确定各个点的参数位置参数位置确定轮廓组垂确定轮廓组垂直方向上的连直方向上的连接接.19相邻轮廓的连通性相邻轮
8、廓的连通性1.1.轮廓轮廓组组的参数的参数化化定义参数:u方向:其中:其中:i=1i=1,2 2,3m3m;k=1k=1,2 2,3n3n;u u1 1=0=0v方向:其中:其中:j=1j=1,2 2,3 3nn;Z Zk k 是第是第k k个轮廓在个轮廓在Z Z方向上的高度方向上的高度 .20相邻轮廓的连通性相邻轮廓的连通性2.2.确定确定各个点的参数各个点的参数位置位置将轮廓线中的点在上一步确定的参数域中进行表述,各个参数点在u、v方向间距均匀地排布:ui=ui-1+uvj=vj-1+vi=1,2,3m;u=1/mj=1,2,3n;v=1/n.21相邻轮廓的连通性相邻轮廓的连通性3.3.确
9、定确定轮廓组垂直方向上的轮廓组垂直方向上的连接连接 参数域以及轮廓上各个点的参数值确定后,每个参数点表示为(ui,vj)的形式,轮廓组可以看作一个矩形二维点阵。点阵中每一行中各个点的vj是相同的,代表各个点在同一个横截面轮廓线中;每一列各点的ui值是相同的,代表各个点在同一个垂直方向的轮廓中。确定轮廓组垂直方向的连接就是将ui值相等的各相邻点进行连接,即将(ui,vj)点与(ui,vj-1)和(ui,vj+1)两点分别连接。.22垂直方向垂直方向的平滑滤波的平滑滤波 相邻轮廓的连通性预设之后,模型在垂直方向的轮廓就确定了,下面就要对垂直方向的轮廓进行平滑滤波,以消除垂直方向上的突变和噪声。垂直
10、方向平滑滤波主要通过计算垂直方向离散曲率寻找最大缺陷点,再运用拉普拉斯滤波方法进行处理以实现平滑处理。.23垂直方向的平滑滤波垂直方向的平滑滤波构建垂直方向轮廓构建垂直方向轮廓线线通过计算离散曲率通过计算离散曲率寻找最大缺陷点寻找最大缺陷点 P Pimpimp在参数化的轮廓组在参数化的轮廓组点阵中寻找点阵中寻找 P Pimpimp的相邻点的相邻点运用拉普拉斯滤波运用拉普拉斯滤波处理处理是否满意是否满意结束结束否否是是.24垂直方向的平滑滤波垂直方向的平滑滤波拉普拉斯算子拉普拉斯算子:松弛因子;迭代次数和松弛因子的值应取适当值,不能过大,过大会引起轮廓线的大幅变化。.25分支表面生成分支表面生成
11、已获取轮廓线已获取轮廓线双向平滑滤波双向平滑滤波分支表面生成分支表面生成平滑的轮廓线组平滑的轮廓线组.26分支表面生成分支表面生成 已经获得了平滑的轮廓,并确定了相邻轮廓的连接,要完成模型表面的生成要考虑两种情况:与横截面一一对应的轮廓 分支区域的轮廓。第一种情况可以采取蒙皮算法直接进行模型表面生产,但是分支区域轮廓拟合成模型表面需要在分支区域构建复合轮廓。传统的复合轮廓生成不能保证分支轮廓与基础轮廓的平滑过渡,这里提出一种新的方法,这种方法引入了边缘桥和二等分线,确保分支区域的平滑过渡。.27分支表面生成分支表面生成1通过在基础轮廓通过在基础轮廓组和分支轮廓组组和分支轮廓组上各做其垂直轮上各
12、做其垂直轮廓在末端点处的廓在末端点处的切向量,得到初切向量,得到初始的复合轮廓。始的复合轮廓。2在初始的复合轮在初始的复合轮廓上构建边缘桥廓上构建边缘桥和二等分线。和二等分线。3通过二分线对初通过二分线对初始的复合轮廓进始的复合轮廓进行修正,并进行行修正,并进行轮廓曲线拟合。轮廓曲线拟合。.28分支表面生成分支表面生成 1.1.初始复合轮廓获取初始复合轮廓获取 Z Zcomcom 平面:平面:垂直轮廓垂直轮廓切向量:切向量:.29分支表面生成分支表面生成 2.2.构建边缘桥和二等分线构建边缘桥和二等分线边缘桥e ei i 是基础轮廓上的点与分支轮廓上点之间的连线,定义公式:其中,D Dbaba
13、是分支轮廓点到基础轮廓的最短连线,D Dbrbr是到另一个分支轮廓的最短连线。二等分线是一系列的边缘桥中点的连线。.30分支表面生成分支表面生成3.3.拟合生成复合轮廓拟合生成复合轮廓 将将基础轮廓和分基础轮廓和分支轮廓上的点的支轮廓上的点的坐标用二等分线坐标用二等分线上的中点坐标进上的中点坐标进行代换行代换基于二等分线对基于二等分线对基础轮廓和两个基础轮廓和两个分支分支轮廓分别用轮廓分别用B样条曲线进行样条曲线进行拟合拟合 最后在峡谷区域最后在峡谷区域将两个复合轮廓将两个复合轮廓再用再用B样条曲线拟样条曲线拟合合,获得,获得分支区分支区域的复合轮廓域的复合轮廓.31分支表面生成分支表面生成.
14、32应用应用验证验证 基于此方法,使用C+和OpenGL图形库,开发了软件包“MediSurf”。其主要模块:(1 1)轮廓提取模块 (2)双向平滑滤波模块 (3)分支表面生成模块 为了验证这种方法和基于其开发出的软件的有效性,使用一系列的股骨的CT进行实验。.33应用验证应用验证1.1.轮廓轮廓数据提取数据提取a.a.原始的股骨原始的股骨CTCT图像图像b.b.阈值图像阈值图像c.c.形态学处理图像形态学处理图像d.d.提取出的轮廓提取出的轮廓.34应用验证应用验证2.2.双向平滑滤波双向平滑滤波a.a.股骨头模型轮廓线;股骨头模型轮廓线;b.b.没有进行滤波时形成的模型表面;没有进行滤波时
15、形成的模型表面;c.c.经过双向平滑经过双向平滑滤波后形成的模型表面滤波后形成的模型表面.35应用验证应用验证 3.3.表面模型生成及应用表面模型生成及应用a.a.从从CTCT图像中提取出的轮图像中提取出的轮廓组;廓组;b.b.最终形成的表面模型;最终形成的表面模型;c.c.有限元模型,可用于医有限元模型,可用于医学分析;学分析;d.d.运用激光烧结成形技术运用激光烧结成形技术制作出的股骨模型,可制作出的股骨模型,可应用于制作医用植入物应用于制作医用植入物和手术仿真和手术仿真。.36结论结论 从医学图像生成三维模型的整个过程包括三个步骤:(1)从医学图像提取轮廓;(2)轮廓的平滑处理;(3)建
16、立表面模型。本文针对第二个步骤提出了两个新的算法:双向平滑滤波算法和分支表面生成算法。进行了应用性的验证,验证表明运用这两种算法的平滑处理方式可以在保证形状细节的同时,有效的去除从医学图像中提取的轮廓线上的噪声和形状突变,进而能够生产较为平滑的医学三维模型表面。目前的研究只是针对骨的外表面进行三维建模,下一步的研究将会对骨的内部表面结构进行三维重建。.37.38后面内容直接删除就行资料可以编辑修改使用资料可以编辑修改使用资料仅供参考,实际情况实际分析.39主要经营:课件设计,文档制作,网络软件设计、图文设计制作、发布广告等秉着以优质的服务对待每一位客户,做到让客户满意!致力于数据挖掘,合同简历、论文写作、PPT设计、计划书、策划案、学习课件、各类模板等方方面面,打造全网一站式需求.40The user can demonstrate on a projector or computer,or print the presentation and make it into a film to be used in a wider field.41
侵权处理QQ:3464097650--上传资料QQ:3464097650
【声明】本站为“文档C2C交易模式”,即用户上传的文档直接卖给(下载)用户,本站只是网络空间服务平台,本站所有原创文档下载所得归上传人所有,如您发现上传作品侵犯了您的版权,请立刻联系我们并提供证据,我们将在3个工作日内予以改正。