血管的三维重建论文.pdfVIP

血管的三维重建 摘要 本文对于血管的三维重建问题，通过分析给定图片的像素数据，通过以下方 法完成了血管的三维重建。 首先我们假定血管为等径管道，将血管看作半径不变的球沿着一条中轴线滚 动形成。根据题目中提供的100 张BMP 图片的像素点坐标数据，通过分析其几 何特性，给出寻找其中轴线和半径的方法，由此找到每张图片中最大内切圆的圆 心和半径，利用平均法，得到半径的平均值约为29.25m. ，即为管道半径。 接下来利用100 个圆心的坐标利用最小二乘法拟合得到中轴线的方程： 7 6 5 4 3 x t = 7.9021t −1.6321t + 0.1378t −0.0056t + 0.0001t 7 6 5 4 3 y t = −2.9429t + 0.6330t −0.0546t + 0.0023t −0.0001t z t = t 通过对圆心坐标点绘制出的曲线与拟合后的曲线进行对比，发现拟合的曲线 与得到的圆心坐标点非常吻合，从而更好的反映了中轴线上各点。然后再用枚举 法以得到的圆心坐标点为球心，画出100 个球面上的点组成的还原图与100 张切 片的边界叠加成的还原图进行比对，形状相符，成功的还原了血管的三维图像。 关键词：血管；三维重建；像素；还原；最小二乘拟合；滚动法；枚举法 1.问题的重述 断面可用于了解生物组织、器官等的形态。例如，将样本染色后切成厚约1 m 的切片，在显微镜下可以观察该横断面的组织形态结构。如果用切片机连续不 断地将样本切成数十、成百的平行切片， 可依次逐片观察。根据拍照并采样得 到的平行切片数字图象，运用计算机可重建组织、器官等准确的三维形态。 假设某些血管可视为一类特殊的管道，该管道的表面是由球心沿着某一曲线 （称为中轴线）的球滚动包络而成。例如圆柱就是这样一种管道，其中轴线为直 线，由半径固定的球滚动包络形成。 现有某管道的相继100 张平行切片图像，记录了管道与切片的交。图像文件 名依次为0.bmp、1.bmp、…、 99.bmp，格式均为BMP，宽、高均为512 个像素 （pixel）。为简化起见，假设：管道中轴线与每张切片有且只有一个交点；球 半径固定；切片间距以及图像像素的尺寸均为1。 取坐标系的Z 轴垂直于切片，第1 张切片为平面Z=0，第100 张切片为平面 Z=99。Z=i 切片图像中象素的坐标依它们在文件中出现的前后次序为 （-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）。 图1 我们需要解决的问题是通过相应的算法，计算出中轴线的坐标以及球心半径， 2 并作出中轴线在XY 、XZ、YZ 面的投影图。 2.问题的分析与假设 2.1 问题的分析 由题意可知，整个管道的表面是由球心沿着某一曲线（中轴线）的球滚动包 络而成，且管道中轴线与每张切片有且只有一个交点，这样经过这个交点必切下 球的最大圆面且这个焦点为球滚动时的球心，因此每张管道切面的最大内切圆就 是经过球心的球切面。 因此如果要求出中轴线的方程，我们需要求出每张图片的最大内切圆的圆心， 内切圆的半径也就是球的半径。另外我们可以根据中轴线绘制在三个平面上的投 影。 现在我们要解决的问题就是找到最大圆的圆心和半径，通过圆心的坐标作出 中轴线，然后绘制中轴线在XY、YZ、ZX 平面的投影图。 2.2 问题的假设 1.血管可视为一类特殊的管道，该管道的表面是由球心沿着某一曲线（称为 中轴线）的球滚动包络而成。 2. 管道中轴线与每张切片有且只有一个交点。