血管的三维重建_图文资料.docx

作者：张雄、李宁娟、贾雪娟 血管的三维重建 摘要 随着现代医学的发展，科学对人类病例的研究不再局限在表面现象，在实际研究中利用断面可了解生物组织、器官等的的横截面形态和结构.从而可大大提高人类对某些疾病的预防和治疗 . 针对这一问题，本文由血管的 张连续的平行切片图象计算血管的中轴线与半径，并绘制血管在三个坐标平面上的投影来探讨血管的三维重建 .由于血管的表面是由球心沿着某一曲线（即中轴线）的球滚动而成，由此我们得出结论 :每个切片一定包含滚动球的大圆，并且他一定为切片的最大内切圆，而最大圆对应的半径即为血管的半径，所以求血管半径就转化为求每一个切片内部的点到切片外部轮廓线的所有最短距离中的最大值即为血管半径 . 本文从 张切片图中随机抽取 张切片图，运用 MATLAB 软件，得到 其最大内切圆的圆心及半径，求取平均值，再用圆心拟合求出中轴线 . 最后根据中轴线求出它在 、 、 平面的投影图 . 关键字 MATLAB 软件 中轴线 半径 平均法 一、 问题重述 断面可用于了解生物组织，器官等的形态 .例如，将样本染色后切 成厚约 的切片，在显微镜下观察该横断面的组织形态结构 .如果用切片 机连续不断地将样本切成数十、成百的平行切片，可依次逐片观察 .根 据拍照并采样得到的平行切片数字图像，运用计算机可重建组织、器官等准确的三位形态 . 假设某些血管可视为一类特殊的管道，该管道的表面是由球心沿着 某一曲线的球滚动包络而成 .现有某管道的相继 张平行切片图像，记录 了管道与切片的交 .图像文件名依次为 、 、 ? ， 格式均为 ，宽，高均为 个象素 .为简化起见，假设：管道中轴线 与每张切片有且只有一个交点 ;球半径固定；切片间距以及图像象素的 尺寸均为 .试计算管道的中轴线与半径，给出具体的算法，并绘制中轴 线在 、 、 平面的投影图 . 二、 模型假设 1.假设管道中轴线与每张切片有且只有一个交点； 2.假设球半径固定； 3.假设切片间距以及图像象素的尺寸均为 ； 4.假设血管无严重扭曲； 5.假设切片拍摄不存在误差，数据误差仅与切片数字图像的分辨率有关 . 三、 符号说明 内点的 X 坐 内点的 y 坐 切片 廓 上的点的 X 坐 切片 廓 上的点的 y 坐 坐标为 的内点到轮廓线的距离 第 张切片图的最大内切圆半径 四、 模型分析 于 个血管的三 重建模型 ，由于血管的表面是由球心沿着某一曲 （即中 ）的球 而成，我 此得出 : 若切片与中轴线有交点，且管道的法向横断面是圆，则该 切片必含有半径与球体相同的最大圆，即 切片的最大内切 ，而最大 的半径即 血管的半径， 圆心则在交点处 .所以求血管半径就 化 求每一个切片内部的点到切片外 利用计算机， 运用 MATLAB 件，搜索出 切片 的最大内切 部 廓 的最大半径 . 的半径，并找到每 切片中 与切片交点的坐 ， 中 坐 ，即 心坐 . 利用 些 坐 ，求出血管的中 .在根据中 求出它在 、 、 平面的投影 . 五、 模型的建立与求解 （1 ）半径和圆心的求取（见附录 1 ） a : 运用 MATLAB  软件将每张切片的 文件转化为  矩阵，  代表黑 色， 代表白色  . 同时将切片的轮廓线也存为 矩阵 . b : 在 ? 小距离  张图片中随机抽取了 张切片的图片（ ），做出它们的轮廓线，找出每个内点距离轮廓线的最 ，即为以这个内点为圆心的最小内切圆的半径；  、 在以内点为圆心的最小内切圆中找出距离最大的那个内切圆，即为这幅图的最大内切圆，该内点的坐标即为圆心的坐标，该距离即为最大内切圆的半径 （见表一） . 表一 最大内切圆 切片号 最大内切圆的圆心坐标的半径 X轴 Y轴 Z轴 0 29.0689 96 257 0 9 29.5367 96 259 9 19 29.9672 96 268 19 29 29.6142 98 290 29 39 29.9362 115 338 39 49 29.6873 146 377 49 59 29.8526 202 411 59 69 30.0134 268 423 69 79 29.7302 361 396 79 89 29.6974 396 369 89 99 30.0000 446 257 99 用算数平均法求取半径 . 即 （ 2）求解拟合曲线的方程及平面投影图 通过表 的数据 , 运用 MATLAB 软件先进行 次线性拟合得 面的投影图，再进行 次线性 拟合得 及 面的投影图和中轴线的空间分布图及拟合方程 . 图依次如下： （附录 2和3） 中轴线在 面的拟合方程： + - 中轴线在 面的拟合方程： + + 中轴线在 面的拟合方