血管重建.docVIP

血管的三维重建 摘要： 为了解决这个问题，我们提出运用“腐蚀的思想”来解决此类问题。 利用现有图片资料重建管道的三维形态，文中首先将图像数字化并提取了相关信息，应用数学形态学中图像腐蚀处理的思想，提出利用逐层去除边缘的算法求解每张图片中图像最大内切圆的圆心(即中轴点)，对这些中轴线上已知的中轴点进行拟合，由于一个只对应于一个和一个,故可分别对其投影分别在、平面上进行多项式拟合，求出和。则中轴线的空间方程即为上两式的联立便得到了血管管道的中轴线函数。 分别将中轴点坐标在、、平面投影拟合得到了在对应平面的投影图，其中最大的内切圆半径的平均值为球的半径r=。 算法中采用了改进的边缘检测算子，该算子能够适合于各种形状图像边缘检测。这一算法在实际求解过程中提高了运算效率和计算精度，具有较强的推广性。 关 键 词：最大内切圆；腐蚀处理；管道重建；边缘检测 1．问题重述 1.1背景 为了认识和研究人体或生物体的有关部分的内部结构和变化，人们利用的是分解和合成的办法，分解的过程是用一簇等间距平行平面将人体和生物体有关部分切成一片片薄片，即分解成切片的序列。切片就是用一组等间距的平行平面将生物体中需要研究的部位切成簿薄的一片，每一片就是生物体某一横断面的图象。按顺序排列起来就形成切片图象序列，或称序列图象。切片的制作过程实际上是一个分解的过程，即将一个空间中的生物体的有关部分，分解为一系列的平面图象。序列图象的计算机三维重建是切片制作的逆过程，它建立了反映生物体内部结构的几何模型或数学模型，利用这种模型，人们可以进行各种分析、处理以及信息的传输等操作。医学序列的三维重建是近些年来医学图像和科学计算可视化领域热门的研究课题。 1.2问题 假设某些血管可视为一类特殊的管道，该管道的表面是由球心沿着某一曲线（称为中轴线）的球滚动包络而成。例如圆柱就是这样一种管道，其中轴线为直线，由半径固定的球滚动包络形成。 现有某管道的相继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平面的投影图。（附100张图片） 2.问题分析 题目中已知某管道的相继100张平行切片图像，记录了管道与切片的交面。根据给定的等间距的平行平面的血管横断面的图象数据,重建血管空间结构的数学模型,、、平面的投影图 对于问题1，调取图像后每张图像为离散的像素点，在各个方向上以相同的速度收缩血管的边缘点，用连续去除边界的方法，将最大内切圆圆心保留在最后的封闭图像内，基于这种形态学中的腐蚀思想可得到最大内切圆圆心所在的点。从而记录下所有图片最大内切圆的半径，求平均值后可得出半径。 对于问题2，在问题1的基础上得出最大内切圆圆心的坐标，建立三维坐标系取z轴垂直于切片，按切片出现的前后顺序将圆心点进行拟合，得出血管中轴线的函数。 对于问题3，利用matlab将得出的中轴线描述的点分别在，，坐标系上拟合，便可得到所要的投影图。 3．模型假设 3.1.为简化模型而做的假设 a.管道中轴线与每张切片有且只有一个交点 b.球半径固定 c.切片间距以及图象象素的尺寸均为1，切片足够薄，其厚度对计算的影响可以忽略不计。 （如果不满足以上的任何一个假设，将有可能使本模型变的相当复杂）。 3.2.为保证模型的合理性而做的假设 d．血管无断裂无严重扭曲，即管道表面光滑且连续； e．对切片拍照的过程中不存在误差，图片分辨率大小对数据造成的误差可忽略不记； f. 切片时，将组织块包埋于石蜡或火棉胶中或以低温冰冻，用切片机切片，外界温度对其误差影响不大。 4．术语及符号说明 中轴点：切片与血管中轴线相交的点。 切面：切片图象上血管边缘所包围的面。 切面围线：图上血管切片的边缘所围成的封闭曲线。 ：原始图像 ：结构图像 ：血管管道的半径。 ：血管管道半径（含管壁厚度） ：图像倒入生成的0—1矩阵 ：矩阵边界点的位置 ：有效点的坐标。 ：中轴点二维的坐标 ：中轴点三维的坐标