血管三维重建模型..docVIP

全国一等奖 血管三维重建模型 指导教师：覃思义 学生：朱天飞 丁理杰 孙红芳 血管三维重建模型 摘要 本文对重建血管三维图象的问题建立了数学模型。 通过对问题的分析，我们认为首先求出半径是解决该问题的关键，在此基础上应该用尽量少的搜索次数找出每个切片的轴心。 我们用“选优搜索”法首先求出了血管半径在28～30个象素点之间。每幅切片图上的轴心一定在“有效象素点”（即血管切片边缘内的点）上，我们对图象的象素点进行了分析，为了减少搜索范围，得到了对“有效象素点”的筛选方法，得到了“筛选点集”。由此确定了一个很小的搜索区域。利用半径的平均值，在此范围内通过较少的搜索次数求得各个切片的与中轴线的交点。 我们对所得结果进行了分析。对明显的偏差点进行了纠偏，对纠偏后结果的进行了有效性检验，所有轴心坐标的平均有效率为99.8％。为了得到较精确的中轴线在X-Y、Y-Z、X-Z面的投影，对纠偏后的所有坐标进行分段线性插值。由此得到了较精确的投影和三维散点图。 最后我们给出了该算法的评价，提出了更有效的搜索算法。 问题的重述与分析 在对生物组织、器官等的形态的研究中，需要知道组织、器官的三维形态。其中的一种途径是通过样本切片来重建其三维图象。其大体的步骤为：用切片机对研究的样本进行平行切片；将切片拍照并采样得到平行切片的数字图象；运用计算机建立三维重建的算法，根据二维数字图象的信息重现原样本的三维形态。 我们现在研究的对象是一种血管样本。这种血管的特点是：血管管道的表面是由球心沿着某一曲线的球滚动包络而成。现在已经得到了该血管的100张平行切片的图象，每张图象的象素点为512×512个，格式为bmp。为了简化，已有下面三个重要的假设： 管道中轴线与每张切片有且只有一个交点。 滚动包络成管道的球的半径是固定。 切片间距以及图象象素的尺寸均为1。 需要建立计算机算法，通过图片所包含的信息，对样本进行三维重建，再现血管的三维形态。最后需要得到血管管道的半径和中轴线三维空间的图象以及它分别在XY、YZ?、XZ平面投影的图象。如何建立有效的算法，是我们所要解决的问题。 问题的一个重要要求是求得中轴线的三维和它投影的二维图象。又根据上面的假设1，可知100张图片中包含且仅包含有100个中轴线上的点，它们是图象包含的最重要的信息，我们称之为“中轴点”。所以，我们将建立算法的目的是“提取每幅图象中的中轴点”。 根据假设1、2，通过分析，我们得到一个重要的结论： 结论一：每幅图上有且只有一个最大内切园，该圆的圆心即是这幅图上的中轴点。我们将每幅图上的这个圆称为“关键圆”。根据这个结论，问题转化为在每幅图上查找其切面围线内的最大内切圆。 在上述的分析基础上，我们另外发现，图象中的对查找中轴点有着决定性作用的象素点集中在关键圆的周围。可先将这些最有效点提取出来，再在点集所在的范围内查找关键圆点。 求出每个面上的中轴点，对其进行处理，可以得到所要求的中轴线图象。 定义与符号说明 中轴点：切片与血管中轴线相交的点。 切面：切片图象上血管边缘所包围的面。 切面围线：图上血管切片的边缘所围成的封闭曲线。 关键圆：切面围线内的最大内切圆。 原图：所提供的100幅切片图。 复合图：将100幅原图重叠，实质是管道在xy平面上的投影图。 R：血管管道的半径。 图阵：原图经过转化的矩阵，大小为512×512，包含着原图象信息。 基本假设 ⑴ 管道中轴线与每张切片有且只有一个交点。 ⑵ 滚动包络成管道的球的半径是固定。 ⑶ 切片间距以及图象象素的尺寸均为1。 模型的建立与求解 1）准备工作 (1)获取初始数据：对所给出的100幅bmp格式的图，通过数学软件matlab的imread命令，转化成matlab能处理的图阵。该图阵是一个元素值只有0 和1 的矩阵，每幅图转换后的大小为512×512。矩阵中为0的点对应图上的灰色度高的点（即切面上的点），元素为1的点对应灰色度低的点（即切面外的点）。 （2）我们将100幅原图重叠，得到一个复合图，如下图 该图是管道在X-Y平面上的投影，同时也将它转换成图阵。 我们认为首先应确定出血管管道半径，然后在此基础上进行中轴点的搜索。 （3）对关键圆进行搜索时，必定使用到原图的象素点的信息。每幅原图有效点（即黑色点）点的数量非常大，如果全使用到的话，其工作量将很大。所以，我们考虑先对图阵上的有效点进行筛选，再作进一步的处理。通过对图象的分析，结合结论一，我们发现切面上分布在最大内切圆附近的象素点，（插图进行说明）尤其是其中分布在切面围线上的点，对这种最大内切圆的查找有着重要的作用。所以，