血管的三维重建

1、目录 问题的提出 血管三维重建的背景 问题分析 模型的假设 模型的建立及求解 模型的检验 问题的提出 断面可用于了解生物组织、器官等的形态。例 如，将样本染色后切成厚约1m m的切片，在显 微镜下观察该横断面的组织形态结构。如果用 切片机连续不断地将样本切成数十、成百的平 行切片， 可依次逐片观察。根据拍照并采样 得到的平行切片数字图象，运用计算机可重建 组织、器官等准确的三维形态。 假设某些血管可视为一类特殊的管道，该管道 的表面是由球心沿着某一曲线（称为中轴线） 的球滚动包络而成。例如圆柱就是这样一种管 道，其中轴线为直线，由半径固定的球滚动包 络形成。 问题的提出 现有某管道的相继100

2、张平行切片图象，记录了管道与切片的 交。图象文件名依次为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

3、，-255，z）（255，255，z）。 试计算管道的中轴线与半径，给出具体的算法，并绘制中轴线 在XY、YZ、ZX平面的投影图。 问题的提出 Z=1Z=0 Z=99Z=98 Z=49Z=50 血管三维重建的背景 这问题的来源于序列图象的计算机三维重建。 序列图象的计算机三维重建是应用数学和计算 机技术在医学与生物学领域的重要应用之一； 是医学和生物学的重要研究方法，它帮助人们 由表及里、由浅人深地认识生物体的内部性质 与变化，理解其空间结构和形态。 血管是血液流通的通路，其在生命活动中的重 要性是众所周知，诊断师在临床中经常需要了 解血管的分布、走向等重要信息。理想的血管 可以看成是粗细均匀

4、的管道，如何建立其数学 模型是图象三维重建的重要一环。 血管三维重建的背景 问题分析 我们将说有100张切片按其在空间的位置 叠加起来，已经可以看到大体的血管的 三维结构。 问题分析 将100张图片叠加在XOY平面，形成血管 在XOY平面上的投影，可以看出血管是 由一系列圆移动形成的。 问题分析 对于每一个切片，仔细观察可以发现： 切片也是由一系列半径不同的圆移动形 成的。 问题分析 现在我们可以知道：血管可以看成无穷多个等 径并且圆心相距无穷小的球包络面组成。因此， 切片上的二维图形就应该是由无穷多个球被截 的圆叠加而成。这些圆都是被截球的大圆或者 小圆，其半径有一极大值R，R同时也是球的半

5、 径。这样一个半径R的圆是球心在切片平面内 的球被截而成的，其圆心为中轴线与切片平面 的交点。假设，中轴线与每张切片有且只有一 个交点，所以每一张切片图上包含且只包含一 个半径为R的圆。我们只要找到这个圆，就可 以定出中轴线与切片平面交点的坐标，用这些 交点坐标我们可以建立起中轴线的空间形态。 问题分析 由以上分析，做出如下的判断： 每个切片中 包含一系列滚球在平面上切出的圆， 其中最大的圆为滚球半径，圆心位于血管的中 轴线上。 血管中轴线在XOY上的投影为所有切片在XOY 平面上叠加形成的阴影的中心线。 在求得所有切片的最大圆圆心后，拟和这些点 形成的曲线就是中轴线。 重建后血管形状应和前面

6、100张切片直接形成 的形状一致（用于进行模型的检验）。 模型的假设 血管粗细均匀，其表面由球心沿某一曲线称 为中轴线)且半径固定的球滚动包络而成。 假设管道中轴线与每张切片有且只有一个交点， 球半径固定，切片间距以及图像象素尺寸均为 1。 切点间距尺寸为1，则假设相继切片间连续， 即相继的切片与中轴线的交点是连续的。 模型的假设(更深一层次） 血管粗细均匀的充要条件是各法截面圆之间不相交，这样 可保证各法截面圆周上的点全落在包络面上 (直观地说， 粗细均匀就是过中轴 上的任意点P处用垂直于 在P点切 线方向的刀片切血管得到的截面是以P为圆心以固定的常数 r为半径的圆，称此为过P点的法截面圆)

7、。为此要求 满足 下列条件： 中轴线 上每一点处的曲率半径大于r 中轴线最窄处的宽度d大于2r 中轴线两端点处的法截面圆不相交。 中轴线 上最窄处的宽度d可以这样决定：当 上两点p，q 的连线垂直于 在这两点处的切线时，或仅垂直在其中一 点处的切线而另一点为 的端点时，称p，q为相关点对。 上可以没有相关点对，也可以不止一对相关点对。如果 上无相关点对，则认为d为无穷大，否则取d为相关点对中 的两点间距离的最小值。 模型的建立及求解 通过以上分析，整个模型的建立分为两个部分：滚球 半径的确定和中轴线的确定。 确定滚球半径的方法有很多：最笨也是最容易想到的 方法就是枚举的方法 求每张切片的图象内

8、的最大 内切圆的圆心时，以位于图象内每一个象素为圆心作 圆遍历所有象素点后再作确定。此种方法，思想简 单，程序简单，但计算量大。 1）平均法 求出每张横断面团象内的最大内切圆半径， 再取r为它们的算术平均值。 2) 抽样法 由于已知 滚动球半径是常数，取前几片切片图象内的最大内切 圆半径的平均值为r的值。 3) 极大似然法 在求得 每一片切片图象内的最大内切圆半径后，进行统计， 以出现频率最大的值为r的值。 模型的建立及求解 为了减小搜索的区域，找到切片中血管所在的最小区 域，然后在此区域里逐个象素的搜索。 模型的建立及求解 %先寻找最小区域 M N=size(Im);xmin=1;xmax=

9、N; for j=1:N for i=1:M if(Im(i,j)0) xmax=j; break; end end end for j=N:-1:1 for i=1:M if(Im(i,j)0) xmin=j; break; end end end ymin=1;ymax=M; for i=1:M for j=1:N if(Im(i,j)0) ymax=i; break; end end end for i=M:-1:1 for j=1:N if(Im(i,j)0) ymin=i; break; end end end 模型的建立及求解 %在最小区域里搜索最 大半径 rmin=1; rmax

10、=min(xmax- xmin)/2 (ymax-ymin)/2); flag=1; while(flag0) R=(rmin+rmax)/2; isR,Ox,Oy=isneiqieyuan (Im,xmin,ymin,xmax,yma x,R); if(isR=1) rmin=R;OX=Ox; OY=Oy; else rmax=R; end D=abs(rmax-rmin); if(D0.1) flag=0; R=rmin; end end 模型的建立及求解 模型的建立及求解 确定中轴线的一些方法： 平行切线法 切片的图象边界上的两点的连线如果同 时垂直边界在这两点处的切线，则这两点连线有可

11、能 是最大内切圆的直径。发现所有具有这样性质的点对， 并检验之，以确定最大内切圆的圆心。 外推法 利用中轴线的连续性，采用插值外推方法， 根据前几片已求得的最大内切圆心位置，推断出新的 一片图象包含最大内切圆心的估计位置，然后经过几 次迭代求得较正确的圆心位置。 滚球法 让球在血管内滚动，保证球与血管相内切， 逐个横断面地定出球心的位置。 模型的建立及求解 投影法 将各切片的图象叠加在XY平面上，形成血管 在XY平面上的投影，其中心线是血管中轴线 在XY平 面上的投影。类似地将血管向XZ平面(或YZ平面)投影， 也可以求得中轴线 在XZ平面(或YZ平面)的投影。这 样做的参赛者为数不少。 模型

12、的建立及求解 变换法 切片图象包含的最大内切圆位置可以理解为图象与固定 半径为r的圆的交的面积达到最大值时图的位置，这可以通过几何 方法实现，也可以将其理解为图象的与半径r的圆的卷积达到最大 值时的情况，可以运用傅里叶变换及其逆变换计算卷积，特别可 以运用快速傅里叶变换的方法。 模型的建立及求解 值得提到的是：由于数据图形（即切片图形）为宽、高均为512个 象素的图像，所以数据精度最高为一个象素。上述的算法得到的 结果却为浮点形的数值，本身存在模型误差，同时有浪费计算时 间。 应用离散模型能一定程度的避免模型误差。针对图像的离散性质， 计算的时候也是使用象素精度，同时做圆的时候也使用Brese

13、nham 算法（详见计算机图形学）来做圆，尽可能的消除数字图像 所具有的锯齿现象。 模型的建立及求解 模型的建立及求解 得到的中轴线在各平面的投影： 模型的建立及求解 模型的建立及求解 模型的建立及求解 对数据点（中轴线坐标）的拟和： 在求得中轴线 与各切片的交点的近似位置 后很多参赛者采用多项式的参数曲线进行拟 合逼近，也有的采用参数样条曲线进行拟合， 这些方法都是可行的。 但要注一点，如果直接拟和三维空间曲线，只 能用参数曲线的形式；要用代数曲线的形式， 只能分别对空间曲线在两个正交坐标平面上的 投影进行拟和，然后将空间曲线表示成两个柱 面的交。 模型的建立及求解 F(x,y,z)=0或z

14、=f(x,y)都表示曲面。 才表示空间曲线。 我们可以简单的用两个柱面的交表示空 间曲线，即 ( , , )0 ( , , )0 F x y z G x y z ( , )0 ( , )0 F x y G y z 模型的建立及求解 模型的建立及求解 参数曲线的形式： 例如螺旋线： 三次B样条曲线： 式中V1，V2，V3，V4为特征多边形顶点，亦称控制顶点；t为参数。 ( )01rr tt cos( ) sin( ) xat ybt zct 1 232 3 4 11223344 1331 3630 1 ( )1 30306 1410 ( )( )( )( ) V V P tttt V V f t

15、 Vf t Vf t Vf t V 模型的建立及求解 模型的建立及求解 模型的建立及求解 重建后的血管三维图像： 模型的建立及求解 模型的建立及求解 简易的血管三维图形绘制效果 模型的检验 本题是应该检验的，不检验只能说完成问题的 一半。一方面，无论以何种方式建模，其过程 都是近似计算，几经近似，效果如何，检验很 必要；另一方面各血管的切面数据已知，按指 定的空间位置放置，就能形成一段血管，完全 可以作为检验的标准。 不论用何种方法建立模型，对照给定血管的切 片图像数据，通过检验，发现模型的误差，修 正模型，可以提高模型的正确性。 模型的检验 检验的主要方法有： 逐片比较 运用求得的滚动球的半径r和中轴线 ，用 球心沿 运动的方法产生一段新的血管：用原来的100 张平面截断血管，生成100张新的切片图像，逐一比较 新、旧100张切片图像之间的差别。若差别越小，则认 为满意，否则要修正模型。这样做的检验效果最好， 但缺点是计算速度慢，时间上可能做不到检验全部100 张切片。 法平面法 设P为中轴线 上任意点，过P点作 的 法平面与原血管相交，若该法平面与血管的交与半径 为r的圆差别小，则认为是满意的，否则修正模型。 模型的检验 模型的假设 血管粗细均匀，其表面由球心沿某一曲线称 为中轴线)且半径固定的球滚动包络而成。 假设管道中轴线与每张切片有且