血管的.docVIP

-------------------------------------------------------------------------------------------------------血管切片的三维重建和分析 摘要： 虽然名称是“三维重建”，但我们考虑问题重点不在于处理三维图形的能力，所以更侧重模型的建立和分析： 管道是等径球的包络，可看作是球沿中轴线滚动形成，则知任意切片可看作是球在滚动过程中渐渐“穿过”平面Zi而留下的“穿孔”（洞），而球穿过此平面的任一时刻与平面相交成一圆，由此知切片边缘Ei是球“穿过”过程中与Zi相交而得一簇圆的包络。 球在穿过此平面时必有某一时刻满足球心恰穿过，即切片上必会留有一个大圆，此大圆属于切片边缘所包络的圆簇，且半径r=球半径R=血管半径R. 球在穿过此平面某时刻（球心不在平面上时），球心到所截得圆的圆心的连线与截平面垂直，即切片边缘所包络的一簇圆的圆心的轨迹是球在穿过此平面过程中球心轨迹在此平面上的射影。 从以上所考虑的两点出发,我们首先给出了一个简易可行的实验算法(模型一),接着又给出了一改进了的强理论高精度的几何算法（模型二）,最后,给出了一个将模型一和模型二较好结合的改进模型(构想三)。其中： 模型一:简单,形象.通过找出边界点后逐点扫描，求得各截面区域内与边界距离最大的内点及其对应的距离(给了两种求法).则知此内点是中轴线与界面交点,所对应的距离必等于球半径R.如此在100个截面上求得了中轴线与之交点,并得到了100个在如题干精确度下较为一致的R值,进而得到球半径R=29.后又利用此100个点,分别在xz,yz平面上拟合了两条x,y关于z的光滑函数,由此知轴线的拟合函数 模型一中关于结果的检验最具特色,得出拟合中轴线后, 让球再次沿此中轴线滚动一次，穿过第Zi平面, 与原第i切片图形比较.若沿所求得的轴心滚动球（即做如上算法）穿截面所得图形与实际图形相差较大，则说明中轴线取法有误（或由于算法不好而导致误差过大）。在我们的检验过程中，每张接受检验的图形都得到了一个较好的“回映”这说明，我们的结果比较好。 另一方面值得指出的是通过此方法我们不但能检验实验所得数据的可靠性，也提供一种求解此问题的新颖的思路。 模型二:可靠,精确.是通过几何理论改进的模型.推导出了中轴线上点与每一切片的中心线（包络圆所在线） (i=0,1……99)为第i张截面 ②N=(0,0,1)为截面单位法向量 ③di=i为切片距Z0距离 ④Di : 截面上切片图形所在区域 ⑤C0(s):血管中轴线 ⑥s:血管中轴线弧长参数 ⑦Ri:第I张截面图上求得球的半径 ⑧R:包络球半径 ⑨M:血管曲面上的点 三 模型的建立及求解 模型一：简易模拟实验性模型 1，模型简介 本模型从“血管可看成等半径的球沿中轴线运动”出发，模拟式的求出球半径R出各截面与中轴线交点，用光滑函数精确的7次拟合中轴线 。 后又从此拟合出的中轴线出发，模拟，让半径为R的球沿此中轴线运动，“穿过”一部分截面Zi ，比较得出的穿孔图形与原Zi平面上切片图形，用二者不同的象素数与原切片图形做比，而确定Zi上所找到的中轴线上点的误差。 2，模型的建立与求解 设球半径R，截面Zi 截血管所得区域为Di ,中轴线与截面Z=i交点为, i=(0,……99) ( 1 ) 求中轴线与截面交点,进而求包络球半径： 解法1 ①将各截面数字化（转化为0，1数据）。 ②圈定Di所在的最小方形区域：求出截面图形的上下左右边界，以此范围圈定为截面图形的最小区域。在此区域中对数字化图形作计算，减少运算量。 ③找出Di 边界点（用离散的象素组成的x-y坐标表示） ④将Di的内点按y值由大到小，x值由小到大的次序逐一扫描，求其与边界的距离。在这所有距离中求最大值，则此最大值及其所对应内点即为R和Ci. 算法为：计算图形内每点到边界点的距离，在此过程中与之前得到的内点与边界点距离的最大值作比较，若小于此值，则退出循环，否则继续。若最后得到它与边界点的距离大于此前所有距离，则将此距离及点坐标记录，并以此为最大距离继续，直至穷尽图形内所有点，如此，最后得到的距离即为轴半径。而最后记录下来的点坐标即为中心轴线与此切面的交点。 ⑤按如上所述算法可编程序，借助计算机求解Ci及Ri，并可对Ri(i=0,……99)取平均值，精确到1个单位，对此典型精确度受象素控制问题，即可求出较精确的R。 （程序详见附录1） 解法2 ①②③同1 ④为提高效率，采取3*3方块中只取中心的方法，而事实上，当此9个点全在图形内部时，与边界距离最大的点只能在中心达到，所以基本上没有损失精度 ⑤改用4中所取小方块中心代替解法Ⅰ中各内点参加如（*）所述循环计算。 ⑥此法缺点：不十分精确 ，但优点不容忽视：节省