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

维普资讯 第 19卷 建模专辑 工 程 数 学 学 报 Vo[19Supp 20 年 02月 Feb 2∞ 2 JOURNAL OF ENGINEERING MATHEMATICS 文章编号 ：1005．3085(2002)05—0035—06 血管的三维重建 徐 晋 (电子工程与信息科学系)， 刘雪峰 (数学系)． 柏容刚(电子工程与信 ．科学系) 指 导老师 ：窦 斗 (数 学 系) (中国科 学技术大学 ．合肥 230026) 编者按：奉文分析了每张切片与营道 曲面的交线是一族 圆的包络缱 ．且这族 圆中半径最大者即为最 大 内切 圆这一几何 特性． 建立 了相应的算法，将 中轴鳍拟合成 B岛；er曲线 ．并用两种算法耐模垂进行了检验 ．方法有一定特色 。 擅 要 !对血管的三维重建问题 ．我们假定血管为等径管道 ．通过分析其几何 特性．培 出了确定其管道 中轴缱和半径的数学 模童——搜索每十切片截面 ．求最大 内切 圆，谈 内切 圆圆心即为切片截面与管道 中轴线 的交点，谈 内切匿半径即为 管道半径 ，再通过拟台各十交点求 出轴心缱 本模型 中，我们确立 了两种有效 的误差分析方法；井 由此发现 由于 中轴线与切 片交角过小会使结果产生较大偏 差 。为解挟此 问题 ．我们从其它方 向重新对血管进行切割 ．再进行处理求解 ，得剜更加精确的结果 美t词：血管 ；等径管道 ；旋转切面 分类号 ；AMS(2000)65D17 中围分类 号 ：0242 1 文鼍标识码 ：A 1 问题重述 (略) 条件假设 1)．血管的表面是 由半径 固定 、圆心连续变化的一族球滚动形成的包络面。 2)医学上 ，血管不存在严重扭 曲。 3)．管道 中轴线与每张切片有且 只有一个交点。 2 问题分析与模型建立 根据假设 ．血管可视为表面是 由球心沿着某一 曲线 (称为 中轴线)的球滚动包络而成 的管 道 。 根据所查文献 (参考文献_6)，可知这种管道有如下几何特性 ： 定理 等径管道每个切片的轮廓线是一族半径 、圆心连续变化的圆的包络线，而这族圆中 半径最大 的圆的圆心即为管道 的中轴线与切片的交点 ．半径 即为管道半径 。 利用这个定理 ．我们建立模型如下 ： 1)对第 个切片 ．求其轮廓线最大内切圆的半径 Rf与 圆心Ci，i=0，1，。一，99； 1 99 2)．由R= ∑R 确定管道半径R； 3)．将各 圆心 G 拟合成 B6zier曲线 ，此 曲线 即为管道 的中轴线 。 4)．确立误差分析方法 ，提 出调整算法 对 2)、3)步结果进行调整 。 维普资讯 工 程 数 学 学 报 第 19卷 3 模型求解 1)．导人数据 、转换存储方式 为方便计算 ，先将 Bmp文件用程序 (本文使用 Mathematica)导人计算机转换成三维矩 阵 存 储 。 2)求截面轮廓线上各点的坐标 由 “管道 中轴线与每张切片有且只有一个交点”，知截面为单连通的区域，轮廓线为一闭合 曲线 。本文使用了 Matlab的edge()函数求 出此轮廓线。由于数据量较大 ，以稀疏矩 阵存取数 据 。 3)．求轮廓线最大 内切 圆的半径与圆心的算法 总体思想 ： 在截面 内选取一定数量 的点 ，对每一个点分别计算它到所有边界点的距离，取距离中最小 的一个 ，这个距离就是以该点为圆Ii,的、与轮廓线 内切的若干个圆中最小的圆的半径 。在这些 由一定数量的点确定 的最