CUMCM优秀论文-血管管道的三维重建【数学建模】.pdfVIP

维普资讯 第19卷 建模专辑 工 程 数 学 学 报 vo1 ． 19 supp 。。 。月 JOURNALOFENGINEERING MATHEMATICS Feb2002 文章编号 ：1005—3085(2002)05．0047—07 血管管道的三维重建 丁峰平 ， 周立丰 ， 李孝朋 指导教师： 数模组 (浙江工业大学 98提高班 ．杭州 310032) 螭者按：该文作者能正确运用几何方法重建了等径管道序列切片的中轴绒模型及曾道形状 其特色是能用所求褥的半径的 球沿着该中轴线蒗 曲重新生成管道 ，井再作裁面 图与原切片囝形作 比较，验证 了模型的可信性 。 摘 要：文章对血管管道 的三维重建进行 了讨 论。根据题 目所 培信息 ．首先读取 100张血管 切面 圉，把它们转换戚数据矩 阵 ．I}!l后分三步进秆处 理：第一步 ．通过搜索切面最大内切圆求出管道 的半径 ．提出两种 方案，分 别是忉缱 {击和摄犬 覆盖法 ；第二步，孰迹 的搜索 ，车文提 出了三种方{击，分别为网格法、蒙特卡罗法和非线性规埘法 ；第三步．中轴线在 三平面上投影的精确定位 ．分*用最小二乘和分段最小二乘进行丁 曲线的拟舍 。最后叉对三堆重建的血管管道避 行了检验和误差分析 。利用 上算法驶好地进行了．管道 的重建，从而得 出所求半径为 29529．中轴线上 lt／0点的坐 标见表 1．其在 ．YZ和髓 平面上的投影分别为囝8到图 l5。 关键词：三维图象重建 ；轨迹 ；最大覆盖 分类号：AMS(2000)65D17 中田分类号 ：02421 文献标识码 ：A 1 问题 的重述 11 问题 断面可用于了解生物组织 、器官等的形态 。例如 ，将样本染色后切成厚约 1．am 的切片 ，在 显微镜下观察该横断面的组织形态结构。如果用切片机连续不断地将样本切成数十 、成百的 平行切片 ，可依次逐片观察 。根据拍照并采样得到 的平行切片数字 图象 ，运用计算机可重建组 织 、器官等准确 的三维形态 。 假设某些血管可视为一类特殊 的管道 ，该管道 的表面是 由球心沿着某一 曲线 (称为 中轴 线)的球滚动包络而成。例如圆柱就是这样一种管道，其 中轴线为直线 ，由半径固定的球滚动 包络形成 。 现有某管道的相继 100张平行切片图象，记录了管道与切片的交。图象文件名依次为 0． bmp、1．bmp、…、99．bmp，宽、高均为 512个象索 (pixe[)。为简化起见，假设 ：管道 中轴线与每 张切片有且 只有一个交点 ；球半径 固定 ；切片 间距 以及 图象象素 的尺寸均为 1。 取坐标系 的 z轴垂直于切片 ，第 1张切片为平面 z=0，第 100张切片为平面 z=99。 试计算管道的中轴线与半径 ，给出具体的算法 ，并绘制中轴线在 XY、YZ、ZX平面的投影 图。 12 切 片的重 组 为了便于直观的分析，我们对切片图象提取轮廓后进行三维重组，得出图 1(为了防止轮 维普资讯 工 程 数 学 学 报 第 19卷 廓过密而不易观察 ，这儿每三幅图片提取一次)。 2 问题的分析 由题 意可 知 ，整个 管道 的表面 是 由 球心沿着某一 曲线 (中轴线)的球滚动包 络而成 ，且 管道 中轴 线 与每 张切片有且 仅有一个交点 ，这样经过这个交点必切 下球 的最大 圆面且这个交点为球滚动时 的球心 ，又 因为每一个切面都 可 以看成 由无数个球切面叠加而成 ，根据两个定