CUMCM优秀论文-利用切片的二维空间相关操作实现血管的三维重建【数学建模】.pdfVIP

维普资讯 第19卷 建模专辑 工 程 数 学 学 报 19supp 。 JOURNALOFENGINEERING MATHEMATICS Feb2002 文章编号 ：1005—3085(2002)05—0029—06 利用切片的二维空问相关操作实现 血管的三维重建 胡亦斌 ， 向 杰 ， 程 翔 指导老师 ： 王若鹏 (北京大学物理 系 ，北京 100871) j胄者按 ：车文利用相关作为两 图形相似程度 的度量，通过快速傅立叶变换厦反变换对切片进行空阃相关操作，具有一定特 色。这种方法能够快速确定出中轴线与切片的交点，所得到的计算结果较为准确。请注意文巾x轴的方向与鼯意 要求是相反的。 摘 要 ：相关(C~ tafion)作为两叶、图形相似程度 的度量 ．筏广泛的用于图形 图像 自动识尉 中。为对血管的二维切片圈像进 行分析并重构出血管 以丑血管中轴线的三维空问形貌 ，我仍利用快速傅立叶变换 (FFT)厦反变换 对切 片进行空间 相关操作，几乎 步即可确定出中轴线与切片的交点，从而蛤 出中轴甥的空间坐标 。我们求出了血管的半径 ，利用 这些结槊，路 出了血管 中轴线 的三摊 曲线皿其投影拨 并且和用计算机软件画 出了血管 的三堆造型 ．在诙造型中柞 血管切片 ．结果与初始的切片数据一致 文中分析了相关法进行 图像处理的优点与局限性 ，对利用近代光学信息处 理的手段进行切片三维重建的思路进行了讨论 差挺宇 ：傅立叶变换 ，相关操作 ．三维重建 分类号 {AMS(2000)65D17 中圄分类号 ：02421 文碱标识码 ：A 1 问题分析与半径的获得 血管可以看成无穷多个等径并且 圆心相距无穷小 的球包络面组成 。因此 ，切片上的二维 图形就应该是 由无穷多个球被截 的圆叠加而成 。这些 圆都是被截球 的大 圆或者小 圆，其半径 有一极大值 R，R同时也是球 的半径 。这样一个半径 R的圆是球心在切片平面 内的球被截而 成的，其圆心为 中轴线与切片平面的交点 。假设 ，中轴线与每张切片有且 只有一个交点，所 以 每一张切片 图上包含且 只包含一个半径为 R的圆。我们 只要找到这个 圆 ，就可 以定 出中轴线 与切片平面交点的坐标 ，用这些交点坐标我们可 以建立起中轴线 的空间形态 。 现在我们把 问题分解为两部分 ：第一步 ，找到半径 R；第二步 ，在每一张切片 图中找到半径 为 R的圆的圆心 。第二个 问题我们使用 了空间相关操作方法米解决 ，这将在下一部分 中详细 描述 。我们先来解决第一个问题 图1是一张叠加 图，是 由0bmp到99bmp的相应像素进行或运算所得。这张图代表了 待求血管在 x—Y平面上的投影的形状 。我们不难看出，除去这一 “弯月”的两端之外 ，这一图 形的宽度是保持不变的，这一宽度就是我们所要求 的球的直径 2R。我们找到直径 的方法是这 样 的：在 图形的中段也就是宽度很均匀部分的中段 (Y=260-280)左侧固定一个在轮廓上的点 维普资讯 工 程 数 学 学 报 第 19卷 A．在 A 点对应的另～侧轮廓线上从与 A 同一 Y坐标 的 一 点 向附近寻找一点 B，使得 lAEl在这些点中最短 ，可 以证明lABl就是我们要求的直径。这样 ．我们就可 以求 出2R=60，R=30(像素)。 2 血管中轴线的确定 在 已获得球半径 的基础上 ，我们对切片 图样进行相 关操作来确定半径为 R的圆的圆心的位置 。 所谓 相关操 作是指利 用两 函数 f(z，)、g(z，) 作相关积分时的特性，在样品图样 中定位某种模版图样 的位置方法 。我