数学建模论文-血管的三维重建.docVIP

数学建模论文 题目：血管的三维重建 组号：19组 学院：理学院 姓名 ：############## 2012年7月13日血管的三维重建 摘要 本论文通过假设血管为一类特殊的管道，该管道的表面是由球心沿着某一曲线（称为中轴线）的球滚动包络而成。就像一个弯曲的圆柱在此基础上，建立管道中轴线参数方程，并综合考虑实际情况中由于切片厚度及数字图像离散化带来的偏差，通过建立模型求解各个切片的最大内切圆，从而找出球半径。问题第一部分需要求出管道的中轴线方程和半径，第二部分需要绘制中轴线在各个平面上的投影。 首先，要将血管切片图片转化为二维矩阵，我们把每一个切片图的BMP文件转化为一个512×512的0,1矩阵（0代表无象素，1代表有）。运用四邻域法找出切片的边界线，存为0，1矩阵。在此基础上建立坐标系，也即得出了血管切片边界线上各个边界点的坐标。 其次，在题设条件下我们可以得出以下几个结论： 能够被切片包含的半径最大的圆的半径等于原始球（形成包络的球）的半径； 不能被包含于切片的圆的半径一定大于原始球的半径； 可以被切片包含的圆的半径一定小于等于原始球的半径； 根据上述结论，对每一个切片，只要找到能够被切片包含的半径最大的圆就能求出原始球的半径。我们设计了一个界心距离法求该半径以及中轴线与各切片的交点坐标，算法概要： 利用 matlab中的edge函数找出一系列的边界点并转化成坐标； 根据界心法以及各边界点坐标找出最大内切圆半径以及圆心； 利用多项式拟合的方法对中轴线上已知的点进行拟合。 由于一个z只对应于一个x和一个y,故可分别对其投影分别在YZ、ZX平面上进行多项式拟合，求出y=f1(z)和y= x=f2(z)。则中轴线的空间方程即为上两式的联立。 用本模型可以对每一个切片求出一个r0(原始球的半径)，共求出100个r0。这些r0的平均值为30.1706，方差为0.0176。可见r0的精度很高。 利用本模型可解决简单的切片三维重建问题，如应用于在医学、地质、地理等领域进行粗略的分析和三维重建。 关键字：matalab软件 四邻域法 数学建模 中轴线 三维重建 问题重述 据拍照并采样得到的平行切片数字图像，运用计算机可重建组织、器官等准确的三维形态。 本题将血管作为一类特殊的管道看待，或者说，将血管理解为一个半径为常数的滚动球随着其球心沿着一条称之为中轴线的曲线运动时所形成的包络。它要求根据给定的100张等间距的平行平面的血管横断面的图象数据，重建血管空间结构的数学模型，即确定滚动球半径和中轴线方程。 现有某管道的相继100张平行切片图象，记录了管道与切片的交。图象文件名依次为0.bmp、1.bmp、…、 99.bmp，格式均为BMP，宽、高均为512个象素（pixel）。为简化起见，假设：管道中轴线与每张切片有且只有一个交点；球半径固定；切片间距以及图象象素的尺寸均为1。 取坐标系的Z轴垂直于切片，第1张切片为平面Z=0，第100张切片为平面Z=99。Z=z切片图象中象素的坐标依它们在文件中出现的前后次序为 （-256，-256，z），（-256，-255，z），…（-256，255，z）， （-255，-256，z），（-255，-255，z），…（-255，255，z）， …… （ 255，-256，z），（ 255，-255，z），…（255，255，z）。 所以我们的问题是：计算管道的中轴线与半径，给出具体的算法，并绘制中轴线在XY、YZ、ZX平面的投影图。 二．基本假设 球半径固定； 切片间距以及图象象素的尺寸均为1； 切片足够薄，其厚度对计算的影响，可忽略不计； 管道中轴线与每张切片有且只有一个交点； 假设中轴线的曲率相对于标准圆的曲率较小； 血管无断裂无突变，即管道表面光滑且连续； 对切片拍照的过程中不存在误差，数据误差仅与切片数字图象的分辨率有关； 术语及符号说明 3.1术语 切片： 指包络表面与切平面相交而成的曲线围成的区域。 原始球： 指沿中轴线滚动形成包络的球。 3.2符号说明 原始球的半径。（由基本假设中2.可知半径r0固定，为常数） 第i片切片最大内切圆的圆心 第i片切片的最大内切圆的半径 第i片切片的第m个边界点的横坐标