论文部分内容阅读
截至2015年,脑血管疾病患病率高达12‰,已严重威胁到我国居民健康,脑血管三维模型的可视化及数值模拟,可为研究学者和临床医生提供重要的诊断和治疗信息。点云数据处理是目前研究的热点,已经应用在VR、自动驾驶等新型技术中,本文将其应用到脑血管领域中。一些未经过分割处理的血管点云数据,难以得到高质量的脑血管模型和进行下一步血流分析,本文拟对具有复杂拓扑结构下的脑血管分割、骨架线提取、光滑表面重建和血流动力学分析展开研究。 针对脑血管点云数据的分割,提出了基于连通性的分割方法。基于连通性的分割算法首先建立覆盖集,对覆盖集进行特征分析,制定覆盖集连通性判断标准,然后确定分割起始点,通过对覆盖集的连通性判断进行点集扩充,当试探集不连通时,则可能存在分叉区域,这时需要根部反向扩充来优化判断,最后对未处理的覆盖集继续进行扩充,直到所有覆盖集被处理,完成脑血管的局部分割。 在分割的基础上,提出了基于空间中值的骨架提取方法。基于空间中值的方法,通过定义局部半径引入排斥力和稠密度权重函数,实现对骨架点的提取;然后将提取的骨架点连接成线,分别对分支上的骨架线连接和分叉处的骨架线连接提出了解决方案;最后完成脑血管骨架点的提取和骨架线的连接。 脑血管表面重建是在骨架线的基础上进行的。首先对骨架线进行分层编码,重组成各个单条分支;再通过相邻等间距采样的骨架点计算中间骨架的切向量和最优截面直径,构造正交截面圆,将轮廓上的点进行三角拼接重构各单分支的表面;然后提出类Delaunay三角剖分方法处理分叉区域的表面,通过将各分支之间的“互吃”关系映射到二维平面,构造类Delaunay三角,再对空洞区域进行填补,完成分叉区域的表面过渡拼接;最后对重构的表面进行Loop细分处理,获得一个光滑准确的脑血管三维表面模型。 从脑血管三维建模中得到相关血管形状与几何参数对血流动力学的研究至关重要,有助于脑血管疾病的临床研究和诊断。通过对关注区域的脑血管局部模型和狭窄情况下进行网格划分、血液流体参数设置和相关边界条件设置,进行了一个心动周期的数值模拟,分别分析了正常情况下和狭窄情况下不同时刻的血流速度、壁面压力和壁面切应力三个动力学参数。