血管管道的三维重建(范文).pdfVIP

第 19 卷　建模专辑 工 　程 　数 　学 　学 　报 Vol. 19 Supp. 2002 年 02 月 Feb. 2002 JOURNAL OF EN GIN EERIN G MA THEMA TICS 文章编号 (2002) 血管管道的三维重建 丁峰平 , 　周立丰 , 　李孝朋 指导教师 : 　数模组 (浙江工业大学 98 提高班 ,杭州 310032) 编者按 :该文作者能正确运用几何方法重建了等径管道序列切片的中轴线模型及管道形状。其特色是能用所求得的半径的 球沿着该中轴线滚动重新生成管道 ,并再作截面图与原切片图形作比较 ,验证了模型的可信性。 摘　要 :文章对血管管道的三维重建进行了讨论。根据题目所给信息 ,首先读取 100 张血管切面图 ,把它们转换成数据矩 阵 ,然后分三步进行处理 :第一步 ,通过搜索切面最大内切圆求出管道的半径 ,提出两种方案 ,分别是切线法和最大 覆盖法 ;第二步 ,轨迹的搜索 ,本文提出了三种方法 ,分别为网格法、蒙特卡罗法和非线性规划法 ;第三步 , 中轴线在 三平面上投影的精确定位 ,分别用最小二乘和分段最小二乘进行了曲线的拟合。最后又对三维重建的血管管道进 行了检验和误差分析。利用以上算法较好地进行了管道的重建 ,从而得出所求半径为 29. 529 , 中轴线上 100 点的坐 标见表 1 ,其在 X Y , YZ 和 X Z 平面上的投影分别为图 8 到图 15 。 关键词 : 三维图象重建 ;轨迹 ;最大覆盖 ( ) 分类号 : AMS 2000 65D17 　　　　中图分类号: O242. 1 　　　　文献标识码 : A 1 　问题的重述 1. 1 　问题 μ 断面可用于了解生物组织、器官等的形态。例如 ,将样本染色后切成厚约 1 m 的切片 ,在 显微镜下观察该横断面的组织形态结构。如果用切片机连续不断地将样本切成数十、成百的 平行切片 ,可依次逐片观察。根据拍照并采样得到的平行切片数字图象 ,运用计算机可重建组 织、器官等准确的三维形态。 ( 假设某些血管可视为一类特殊的管道 ,该管道的表面是由球心沿着某一曲线 称为中轴 ) 线 的球滚动包络而成。例如圆柱就是这样一种管道 ,其中轴线为直线 , 由半径固定的球滚动 包络形成。 现有某管道的相继 100 张平行切片图象 ,记录了管道与切片的交。图象文件名依次为 0. ( ) bmp 、1. bmp 、…、99. bmp , 宽、高均为 512 个象素 pixel 。为简化起见 ,假设 :管道中轴线与每 张切片有且只有一个交点 ;球半径固定 ;切片间距以及图象象素的尺寸均为 1 。 取坐标系的 Z 轴垂直于切片, 第 1 张切片为平面 Z = 0 , 第 100 张切片为平面 Z = 99 。 试计算管道的中轴线与半径, 给出具体的算法, 并绘制中轴线在 X Y 、Y Z 、ZX 平面的投影 图。 1. 2 　切片的重组 为了便于直观的分析 ,我们对切片图象提取轮廓后进行三维重组 ,得出图 1 ( 为了防止轮 © 1994-2007 China Academic Journal Electronic Publishing House. All rights reserved. 48 工 　程 　数 　学 　学 　报 　　　　　　　　　　　　　　第 19 卷