对血管切片做三维重建的一种方法

《对血管切片做三维重建的一种方法》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、对血管切片做三维重建的一种方法　　Abstract:Firstly，takinguseofsearchingalgorithmtocalcu-latecoordinateofpointofintersectionforkindsofbloodvessel’saxisdemarcationlineandradiusofincircle.Andthen，averagingtheradiusofincircleulatingallpointsofintersection.　　0引言　　在医学、天文观测、工业非破坏性试验等一些工业问题中

2、，需要确定某个空间物体的形状，但由于技术上的限制，无法将物体分离出来并直观的显现.解决此类问题的一个可行的方法是用等间隔的平行平面去截取这个物体，得到一组平行截面，通过采样得到平行截面的数字图像，进行数据处理后可运用计算机重建物体的三维形态.在2001年全国大学生数学建模竞赛题中，A题给出了某血管管道的相继N张平行切片图像，图像格式、尺寸及坐标均已给定，这些图片记录了血管管道与切片的交点，并且假定该血管管道是由一个小球沿某路径移动而形成的.取z轴垂直于切片，规定第一张切片平面为z=0，第N张切片平面为z=N-1.要求解决的主要问

3、题是:求解血管管道的中轴线与半径，同时给出具体算法;以及绘制中轴线在XY，YZ，ZX平面的投影图.现在我们试图利用这组平行截面所提供的信息来重建被截血管管道.　　1一些假设　　题中已给出的假设:（1）管道的中轴线与每张切片有且只有一个交点;（2）形成管道的球半径固定;（3）切片间距以及图像像素的尺寸均为1.　　求解所需的假设:（1）假设球是沿光滑曲线运动的;（2）将像素抽象成为一个点，由于像素宽度与切片间距相等，所以由原数据得到的是一个均匀的三维点阵.（3）像素点的坐标位置在像素点的中心.　　2问题分析　　此类问题的关键是:如何

4、利用所已知的一组等间距的平行截面，根据其平行截面的形状确定球心轨迹（此轨由于采用本文所用的方法，参与“2001年全国大学生数学建模竞赛”的建模小组获得了“国家一等奖”.迹是一个工程上可接受其精度的近似值）.由于问题中假设管道中轴线与每张切片有且只有一个交点，则切片图像在XZ，YZ平面的投影是以z为自变量的函数.可以通过求出单个截面与中轴线的交点，再将所有截面与中轴线的交点进行曲线拟合得到所要求的结果.具体分析思路如下.　　2.1对单个截面的分析文中所给的单个截面，我们可以认为是形成管道的球体沿某一空间曲线穿过一截平面，球体被截平

5、面连续截取所得的圆在该截平面上形成的包络即为截面.由截面形成过程可知，截平面所截得的包含球心的圆就是截面的最大内切圆，其圆心坐标为截面与中轴线的交点坐标.由单个截面所给的信息，我们便可确定每个切片与管道中轴线的交点，并通过求解其最大内切圆的平均半径得出球体半径.形成管道的平移球体穿过截面时，所形成的相交图形是一系列半径可变的圆，可设其半径函数为r=r（t），t1　　2.2对所有截面的分析单个截面所提供的是局部坐标所包含的内容，按上述同样思路求得所有切片与中轴线的交点.所有截面的信息则包含了连接各切片的中轴线的内容，将所有的交点拟

6、合后便可得到中轴线的近似轨迹.　　3问题求解　　基于上述分析，对于一个给定的截面形状，就可以得到这个截面与管道中轴线的唯一交点，若求出每个截面与中轴线的交点，则通过空间曲线的拟合就可以得到原管道的中轴线.　　3.2搜索最大内切圆圆心Ci由切片及截面形成的过程可知，在截面上必定存在一个最大内切圆，中轴线与截面的交点即为该内切圆的圆心.　　算法描述:（1）将边界点的坐标按顺时针（逆时针也可）的方向依次写入一个坐标序列（Xi，Yi）（i=0，1，2……N-1，N为边界点总数）中;（2）对任意一个截面内的点（即黑色的点），记录其到边界的

7、最短距离dminp;（3）扫描所有截面内的点，对于每一个点按上述方法求出到边界的最短距离，找出这些最短距离中最大的一个所对应的点，则此点即为最大的内切圆圆心，亦即管道中轴线与此截面的交点Ci.　　为了进一步精确求得圆心位置，我们进行了以下修正:将所取的参考点在其周围的8个方向上都减小1/2步长，再以上述同样的方法作圆，这样所得的半径最大的圆就与最大内切圆更加接近，所确定的该参考点即为圆心.　　3.3拟合中轴线通过MATLAB软件中相应功能，应用最小二乘拟合将所有切片与中轴线交点在空间进行拟合.具体拟合方法如下:对于上述方法求出的

8、C（xi，yi，zi），i=0，1，2，…N-1，将其投影到各个坐标平面上，即得到XY平面上的（xi，yi），YZ平面上的（yi，zi），XZ平面上的（xi，zi），应用最小二乘法将二维点用一个多项式拟合，应用MATLAB的PULYFIT函数选取适当的拟合次数进