数学建模竞赛血管的三维重建模型g.docxVIP

血管的三维重建模型 摘要： 本文对血管三维重建中，中轴线及球的半径确定问题进行了讨论。首先，根据问题及图象处理提取有效数据，给出两种可行算法，利用上述数据建立了最大最小方法和二次规划方法。 搜索中心点，并给出全局和局部搜索，得到各切片中心点坐标 (见表 1) ，并通过插值方式得到中轴线图象及其各投影。最后对模型给出检验方式。 一 、问题的重述 假设某些血管可视为一类特殊的管道 ,该管道的表面是由球心沿着某一曲线 ( 称为中轴线 )的球（命名为包络球）滚动包络而成。 现有某管道的相继 100 张平行切片图象 ,记录了管道与切片的交。 假设：管道中轴线与每 张图片有且只有一个交点；球半径固定；切片间距以及图象象素的尺寸均为 1. 取坐标的 Z 轴垂直于切片，第 1 张切片为平面 Z 0 ，第 100 张切片为平面 Z 99 . 计算管道的中轴线与半径，给出具体的算法，并绘制中轴线在 XY 、YZ 、 ZX 平面的投 影图。 二、模型假设与符号说明 1、 基本假设： （1） 该管道的表面为一定长半径的球沿一固定的曲线运动所得曲 面族包络的光滑表面。 （2） 该管道的中轴线连续而且光滑。 （3） 该管道的中轴线与每个切面有且只有一个交点。 （4） 图象象素的尺寸为 1. （5） 切片的间距尺寸为 1. 2、 符号说明： Oi  L R x, y, z  中轴线 包络球的半径 中轴线与第 i 个切片的交点（定为此切片的中心） Si  第 i 个切片切得的图形 D i  第 i 个切片的图象数据矩阵 三、问题分析及建模准备 问题分析： 通常血管的表面可认为是连续且光滑的曲面，断面可用于了解其形态等特性。本问题给 出的是一些离散的切面， 要求重建出原图中轴线和求出包络球半径。因为每一个切面与 中轴线 L 有且只有一个交点 Oi ，如果找出所有 Oi , 就可以用插值或拟合的方式作出  L 的 近似图象，其在坐标平面上的投影就很容易画出。 问题的关健转变为求每个平面上的 Oi . 建模准备： 1、 图象的读取 由于切片图象中只有黑、白两种颜色的象素 , 而且所给的 BMP格式图象文 件是 512× 512 象素的 . 因此 , 把图象读取为一个 512× 512 的数字矩阵； 用数字 1 表示黑 色的象素，用数字 0 表示白色的象素。如图 1 把 (a) 的图象读取为 (b) 的矩阵 : 1 1 1 1 1 1 1 1 1 (a) 图 1  (b) 2、 边界提取 512×512 的数字矩阵所占的空间较大 , 不便处理 . 为节省空间 , 只用提取图象的边界进行处理 , 就可达到目的 . 下面给出两种提取边界的方法： (1) Roberts 交叉算子方法 如图 2, a1 , a2 , , a8 为相邻的象素 , 设黑色象素的值为 1, 白色象素的值为 0. a0 a1 a2 对于中间象素 a8 a7 a8 a3 当 K a0 a4 2 2 0 时, a6 a5 a4 a2 a6 可判断其为边界点 . 用此方法对 0.bmp 判别边界 , 的出的边界图如图 3. 图 2 200 190 180 170 160 150 140 130 120 -50-40 -30 -20-10 010 203040 50 图 3 可见，边界点重复 , 多达 512 个, 对求轴心不利 . 自行设计一种判别法 : 为使边界细化给出如下判别方法 : 对 3×3 膜, 当中心点 a8 为黑色象素且白象素≥ 4 时, 方判其为边界点 , 如图 4. a0 a1 a2 a0 a1 a2 a7 a8 a3 a7 a8 a3 a6 a5 a4 a6 a5 a4 （1） （2） 图 4 而对白象素 =2 或 3 时, 则要用 5× 5 膜类似判定 此时主要处理非凹图形的凹边提取 . 对 0.bmp 实施此种边界提取方法后 , 边界点降至 82 个.如图 5 190 180 170 160 150 140 130 120 -20 -10 0 10 20 30 40 -30 图 5 四 模型建立与求解： 1 模型的理论基础： 命题 1：由球心沿着某一曲线（称为中轴线）的半径不变的球（即包络球）包络而成的曲面，若被某一平面横截，且截面与中轴线有且仅有一个交点，则此截面中包含的 圆中半径最大者 r 即为包络球半径。 O，球O以 证明：如图 6，截面与中轴线的交点为 O 为球心，以 R 为半径。则经过 O 的所有截面必包含 球 O 的大圆，所以包络球半径 R r . 若 R r , 则球 O 与切面边界点必无交点，由中轴 O 线的连续性，容易知道，曲面必不是由球 O 包络而成， 产生矛盾。所以， R r