【-数学建模竞赛】血管的三维重建模型g.docVIP

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