毕业论文--血管三维重建.docVIP

血管三维重建 作者：汤志高 等 ［摘要］ 本文通过对样本断面的数字化处理，运用计算机重建血管三维形态，得出血管轴线方程和管体半径。 对断面特征图像文件数值化：视位图图像黑色象素点为0，白色为1，在MATLAB中建立图像矩阵A。对A进行转化，建立图像轮廓矩阵，通过绘制出血管三维重建图形。 在求解管道半径时，结合三维图形与断面图形特征，发现管道半径即为断面图形最大内切圆半径，轴线与断面的交点即为最大内切圆圆心，并在6.1.1中给出了严格证明。就此，将问题转化为求解断面的最大内切圆的半径。 基于转化的目标，首先加入内圆约束（内切圆半径不大于圆心到断面边界上任意一点的距离）；然后采用约束算法“四向扫描法”（见6.4.1）使约束加细，限制圆心坐标。主要原理是对图象矩阵A再次转化，将每个断面的边界都分割成四部分，分离出各断面图形的上、下、左、右轮廓，建立四向轮廓矩阵，对圆心坐标从四个方向约束，当圆心坐标同时满足四个约束时，圆心一定位于断面内部。最后，以为决策变量，内切圆半径最大为目标建立非线性规划模型，采用LINGO软件编程求解，最终管道半径为29.7069；坐标集（见附录10.1）。 在求解中轴线的过程中，首先建立向量（即各内切圆圆心的z坐标），运用100张断面与中轴线的交点在和平面的投影，建立面的投影坐标集，采用高次非线性多项式拟合得到平面曲线方程，同理对面的投影坐标集采用同样方法得到平面曲线方程，以z为参数，由经过上述两条曲线的空间曲面方程的联立，建立了中轴线方程。 最后根据空间中轴线的方程，可以分别画出中轴线在面，面，面的轴线投影图和中轴线的空间图象。 问题重述 断面可用于了解生物组织、器官等的形态。例如，将样本染色切成厚约1? m的切片，在显微镜下观察该横断面的组织形态结构。如果用切片机连续不断地将样本切成数十、成百的平行切片， 可依次逐片观察。根据拍照并采样得到的平行切片数字图象，运用计算机可重建组织、器官等准确的三维形态。 假设某些血管可视为一类特殊的管道，该管道的表面是由球心沿着某一曲线（称为中轴线）的球滚动包络而成。例如圆柱就是这样一种管道，其中轴线为直线，由半径固定的球滚动包络形成。 现有某管道的相继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平面的投影图。 模型假设 血管管道粗细均匀； 切片相继平行； 切片厚度相同； 管道壁表面光滑； 样本染色之后对拍照采样获取数据、切片厚度没有影响； 管道中轴线与每张切片有且只有一个交点； 切片图像真实反映切片原貌及性质； 符号说明 管道半径 切片的平面区域 第i张切片的最大内切圆的半径 G 中轴线在XY平面投影方程的系数矩阵 E 中轴线在XZ平面投影方程的系数矩阵 C 中轴线在YZ平面的投影方程系数矩阵 问题分析 用切片机连续不断地将血管样本切成数十、成百的平行切片， 可依次逐片观察。根据拍照并采样得到的平行切片数字图象，运用计算机可重建血管准确的三维形态，克服了直接通过显微镜等仪器不能直观的总体观察血管形态的缺点。 提取重建血管的关键为：血管中轴线和固定半径的确定。主要有下面几个要点： 血管的重建 运用计算机重建血管的三维形态的过程，是利用从大量的切片的数字图像中获取数据，把切片有序排列起来，本质上是在计算机中将血管数字化。 对管道半径的换位思考 “血管管道”示以人的概念比较复杂，可以将其表面视为由球心沿着某一曲线（称为中轴线）的球滚动包络而成，管道半径就是这个球的半径。设想球沿着我们用若干切片重组的管道顺势滚动，球顺利通过每一张切片的必要条件是：每一张切片的最大内切圆的半径等于球半径，也就是管道的半径。就此，将空间问题转化在平面内解决——在切片所在的平面求解管道半径。 中轴线如何确定 管道中轴线与每张切片有且只有一个交点，平行切片是相继的，相继交点的轨迹即为中轴线的走向。当切片的数量有限时，所有交点亦不足以精确确定中轴线的方程，可以通过拟合空间曲线的方法近似求解。 血