基于ICP的脑部血管配准系统技术方案

本发明专利技术公开一种基于ICP的脑部血管配准系统，包括点集获取单元、点集对应单元、点对筛选单元、向量配准单元以及误差估计单元，点集获取单元用于获取术前的MRA影像以及术中的多普勒超声影像，获取静态点集和动态点集；点集对应单元根据动态点集在静态点集中匹配距离最近的点，两个点构成一个点对；所述点对筛选单元使用最小截平方和算法除去离群点对；所述向量配准单元计算最小平方配准向量并应用于原点集得新的点集；所述误差估计单元计算所得新的点集的与原点集距离的均方误差，并与设定的容忍度τ进行比较，多次迭代后所得到的新的点集接近于大脑移位后的实际血管形态，解决大脑移位对手术过程的影响，提高神经导航系统影像信息的准确性。

Brain blood vessel registration system based on ICP

The invention discloses a brain vessel registration system based on ICP, including the point set acquisition unit, set the corresponding unit, for screening unit, the registration unit vector and the error estimation unit set acquisition unit is used for acquiring MRA images of preoperative and intraoperative ultrasound Doppler shadow image, obtain the static and dynamic point point set set point set; the corresponding unit according to the dynamic point set matching the nearest point in the static point, two points form a pair; the point to remove outliers on the screening unit using the minimum cut and square algorithm; the vector registration unit calculates the least square registration vector is applied to a new set of origin the set of points; the error estimation unit calculated new point set MSE distance and the origin, and compare with the set tolerance by new points obtained after several iterations of the close set The effect of the brain displacement on the operation process and the accuracy of the image information of the neuronavigation system are improved after the displacement of the brain.

【技术实现步骤摘要】
基于ICP的脑部血管配准系统
本专利技术涉及医疗
，特别是指一种基于ICP的脑部血管配准系统。
技术介绍
在脑部手术过程中，医生往往需要参考术前获取的头部MRI影像进行手术。然而由于脑脊液的丢失、患者体位乃至药物的影响，都会导致大脑发生移位和变型。大脑移位使得神经导航的准确性下降，会对手术过程造成相当程度的困难，对病灶的定位造成影响。目前针对大脑移位的配准方法主要分为两大类，分别是基于术中影像的配准方法以及基于生物力学模拟的配准方法。前者主要利用MRI、CBCT以及超声的多种影响技术，使用基准点对脑组织进行配准。后者主要通过构造生物力学模型模拟大脑的各种生物力学特征，将计算所得的变力作用于大脑模型上，获得配准后的脑部结构。
技术实现思路
针对
技术介绍
中存在的问题，本专利技术的目的是提供一种基于ICP的脑部血管配准系统，根据术中获得的多普勒超声影像对术前获得的MRA影像进行校正，使得MRA影像与患者脑部移位情况匹配，解决大脑移位对手术过程的影响，提高神经导航系统影像信息的准确性。本专利技术的技术方案是这样实现的：一种基于ICP的脑部血管配准系统，包括点集获取单元、点集对应单元、点对筛选单元、向量配准单元以及误差估计单元，其中，所述点集获取单元：用于获取术前的MRA影像以及术中的多普勒超声影像，然后使用ITK对血管中心线进行提取，获取的两组中心线便是两组点集，术前MRA获得的点集是静态点集，术中多普勒超声获得的是动态点集；所述点集对应单元：根据动态点集，在静态点集中匹配距离最近的点，两个点构成一个点对；所述点对筛选单元：使用最小截平方和算法除去离群点对；所述向量配准单元：计算最小平方配准向量并应用于原点集得新的点集；所述误差估计单元：以上三步为一个迭代流程，经过一次迭代后，计算所得新的点集的与原点集距离的均方误差，并与设定的容忍度τ进行比较，如果均方误差小于容忍度τ，或者迭代次数达到预设上线，则结束迭代；否则重新进行迭代；经过多次迭代后，所得到的新的点集接近于大脑移位后的实际血管形态。在上述技术方案中，所述术前MRA影像需要使用区域生长法进行图像分割得到MRA血管树，术中获取的多普勒超声影像直接使用阈值分割，将图像中灰度小于一预设阈值的像素移除。在上述技术方案中，所述区域生长法以MRA图像中的一个像素点为基点，检查该基点相邻的像素，将其周围的像素与其进行比较，如果差值小于预设的阈值，则将其合并；再以新像素点为基点，进行上述检测，直到周围没有相似点为止。在上述技术方案中，所述最小截平方和算法计算点对距离的最小协方差矩阵来估计稳健距离，将稳健距离超过预设阈值的点对删除。本专利技术基于ICP的脑部血管配准系统，包括点集获取单元、点集对应单元、点对筛选单元、向量配准单元以及误差估计单元，点集获取单元用于获取术前的MRA影像以及术中的多普勒超声影像，获取静态点集和动态点集；点集对应单元根据动态点集在静态点集中匹配距离最近的点，两个点构成一个点对；所述点对筛选单元使用最小截平方和算法除去离群点对；所述向量配准单元计算最小平方配准向量并应用于原点集得新的点集；所述误差估计单元计算所得新的点集的与原点集距离的均方误差，并与设定的容忍度τ进行比较，多次迭代后所得到的新的点集接近于大脑移位后的实际血管形态，解决大脑移位对手术过程的影响，提高神经导航系统影像信息的准确性。附图说明图1为本专利技术基于ICP的脑部血管配准系统流程图；图2为本专利技术中迭代流程图；图3为本专利技术一具体实例迭代前的点集图；图4为本专利技术一具体实例迭代后的配准结果图。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。本专利技术所述的一种基于ICP的脑部血管配准系统，首先获取多普勒超声影像中的特征点集，随后在MRA血管树中寻找特征点集的最近点，构成点对；使用最小截平方和除去离群点对后，计算各个点对的配准向量，再将配准向量应用于特征点集中。将以上步骤进行多次迭代，便可完成多普勒影像对MRA血管树的配准。本专利技术基于ICP的脑部血管配准系统包括点集获取单元、点集对应单元、点对筛选单元、向量配准单元以及误差估计单元，以下是结合图1所示的本专利技术流程图对上述各单元的详细说明。(1)点集获取单元：需要两种图像作为输入，分别是术前的MRA影像以及术中的多普勒超声影像。术前MRA影像需要进行图像分割，所使用的方法是区域生长法。区域生长法以MRA图像中的一个像素点为基点，检查该基点相邻的像素，将其周围的像素与其进行比较，如果差值小于预设的阈值，则将其合并；再以新像素点为基点，进行上述检测，直到周围没有相似点为止。经过区域生长，得到MRA血管树。而术中获取的多普勒超声影像直接使用阈值分割，将图像中灰度小于一预设阈值的像素移除。获得术前与术中的血管影像后，使用ITK对血管中心线进行提取，获取的两组中心线便是两组点集，术前MRA获得的点集是静态点集，而术中多普勒超声获得的是动态点集。(2)点集对应单元：根据动态点集，在静态点集中匹配距离最近的点，两个点构成一个点对。设动态点集中的一个点为整个静态点集设为点集A，其中包含NA个点。则点到点集A的距离为：为欧几里得距离。最近的一个点满足条件：将匹配为一个点对。记Y为A的最近点的点集，记为最近点算子，则有：(3)点对筛选单元：使用最小截平方和算法除去离群点对。最小截平方和算法主要通过计算点对距离的最小协方差矩阵来估计稳健距离。最小协方差矩阵为：其中：通过MCD求得稳健距离后，将稳健距离超过预设阈值的点对删除。(4)向量配准单元：最小平方四分数操作记为：其中dms是平均平方点匹配误差。计算最小平方配准向量：将配准向量应用到点集P上，可得新的点集P：(5)误差估计单元：以上三步成为一个迭代流程。经过一次迭代后，求一次所得新的点集P的与原点集距离的均方误差，将这一所得的均方误差与设定的容忍度τ进行比较，如果均方误差小于容忍度τ，或者迭代次数达到预设上线，则结束迭代；否则重新进行迭代。均方误差的计算公式为：经过多次迭代后，所得到的新的点集越来越接近于大脑移位后的实际血管形态。其中，迭代流程图如图2所示。以下是结合一具体实例对本专利技术的进一步说明，其中，迭代前的点集如图3所示：如下表所示，是动态点集中某点到静态点集中各点的距离：1234561.42.31.74.33.13.2可见该点到静态点集中的点1距离最近。经过点对匹配后，得到一系列点对，这些点对之间的距离如下：11.421.732.142.558.9....通过最小截平方和的计算，除去当中的离群点(如第五对点对)。随后根据这些点对计算配准向量：将这些配准向量应用于原来的静态点集上，得到新的静态点集Pk+1。设定的迭代次数上限为100，容忍度为0.001，如果新的静态点集Pk+1与上一个静态点集之间的MSE小于容忍度0.001，或者迭代次数达到100，则结束迭代，得到的新的静态点集也就符合目前的患者的脑部血管情况。迭代后的配准结果如图4所示。本专利技术基于ICP的脑部血管配准系统，相较于现有技术，具本文档来自技高网...

【技术保护点】
一种基于ICP的脑部血管配准系统，其特征在于：包括点集获取单元、点集对应单元、点对筛选单元、向量配准单元以及误差估计单元，其中，所述点集获取单元：用于获取术前的MRA影像以及术中的多普勒超声影像，然后使用ITK对血管中心线进行提取，获取的两组中心线便是两组点集，术前MRA获得的点集是静态点集，术中多普勒超声获得的是动态点集；所述点集对应单元：根据动态点集，在静态点集中匹配距离最近的点，两个点构成一个点对；所述点对筛选单元：使用最小截平方和算法除去离群点对；所述向量配准单元：计算最小平方配准向量并应用于原点集得新的点集；所述误差估计单元：以上三步为一个迭代流程，经过一次迭代后，计算所得新的点集的与原点集距离的均方误差，并与设定的容忍度τ进行比较，如果均方误差小于容忍度τ，或者迭代次数达到预设上线，则结束迭代；否则重新进行迭代；经过多次迭代后，所得到的新的点集接近于大脑移位后的实际血管形态。

【技术特征摘要】
1.一种基于ICP的脑部血管配准系统，其特征在于：包括点集获取单元、点集对应单元、点对筛选单元、向量配准单元以及误差估计单元，其中，所述点集获取单元：用于获取术前的MRA影像以及术中的多普勒超声影像，然后使用ITK对血管中心线进行提取，获取的两组中心线便是两组点集，术前MRA获得的点集是静态点集，术中多普勒超声获得的是动态点集；所述点集对应单元：根据动态点集，在静态点集中匹配距离最近的点，两个点构成一个点对；所述点对筛选单元：使用最小截平方和算法除去离群点对；所述向量配准单元：计算最小平方配准向量并应用于原点集得新的点集；所述误差估计单元：以上三步为一个迭代流程，经过一次迭代后，计算所得新的点集的与原点集距离的均方误差，并与设定的容忍度τ进行比较，如果均方误差小于容忍度τ，或者迭代次数达到预设上线，则结束迭代...

【专利技术属性】
技术研发人员：彭逢安，陈卫单，李卓妍，
申请(专利权)人：广州慧扬健康科技有限公司，
类型：发明
国别省市：广东,44