数学建模常用模型与案例整理.pdfVIP

数学建模常用模型与案例整理 薛力源 谭楚婧 夏亮军 建模基本方法：机理分析，测试分析 最常见的建模类型：优化问题，或结合模拟建模 案例1：血管三维重建 本题要求计算 100 张平行血管切片的中轴线与半径并给出具体算法，最后用 MatLab 重构整个血管模型，本题的第一个突破点是一定要意识到每一个平行切痕 的最大内切圆半径均相等，且等于管道半径。这主要是由于在假设中假定血管的 每一个正切切片是标准正圆（或看做小球沿一定曲线运动包络形成，求新的轨迹 即是中心轴线），因此本问题转换成为如何内在一张二值图像中快速找到最大内切 圆的算法设计，并最后通过对各个切片圆心在空间的曲线进行插值确定其中轴线。 仔细观察第一张图可以发现其近似正切，因此可以依据此基本确定半径的取 值范围，其大约在28-32 之间，而之后找寻圆心和半径算法可如下设计： 1.确定每一张切片二值图像的外接四边形，将四边形向内缩小约 26 个像素， 计算边界像素到小四边形内各个点之间的距离() ，记第一次迭代 (1) = 28，淘汰 () (1) ̅ 掉对于小四边形内像素0 ，存在 的像素点，此时剩下的 即是可能的圆心。 ̅ (2) 2.在所有 中迭代 = 29，重复进行以上过程，再次淘汰部分像素点，剩下 的像素点集合即是可能的圆心。 { } 3.依次迭代直至 = ∅，此时，上一轮迭代所剩余的像素点即是可能的中心。 最后可以计算出100 张图像的半径约为30，至此已经完全得到了关于血管中 轴线的坐标数据(, ) ，基于此在空间中绘制出中轴线的曲线，将其分别投影到xoy ， yoz ，xoz 平面运用三次样条法分别进行插值即可最终得到平滑的血管切片中轴线。 然而此种方法的缺点在于速度缓慢且效率不高，且投影到三个平面分别做三次样 条插值的做法极其不可取，因此在这里我们给出用B 样条曲线逼近的方法计算出 其中轴曲线。 B 样条曲线是广泛应用的形状数学方法，在三维医学图像学、人体解剖学、机 械设计等学科均有应用。 p 次B 样条曲线的方程为 n ( ) C u = ∑ ( ) ， 0 ≤ ≤ 1 ， , i=0 其中 为控制顶点； ( ) 为 p 次样条基函数，由节点矢量 = { = ⋯ = , 0 , , …, u , = ⋯ = }确定。先计算数据点的参数值u̅ 和节点矢量就 n +1 ++1 i 可以建立并求解位置控制顶点的线性方程组。 在模型和算法的改进部分，应当注意到本快速算法的两个问题，首先，将原二 值图像边界外切四边形向内缩小约 26 个像素缩小圆心的所有范围的办法对于具 有近似对称关系的图像满足需要，但若图像不具有对称关系，此方法将有可能完 全「躲开」最大内切圆圆心，如下图 检索四边形 最大内切圆 外接四边形 一个好的办法是设定一个对称性指标作为外接四边形内陷像素点的个数，对