[发明专利]基于直接线性变换和奇异值分解的磁共振坐标定位方法

说明书

技术领域

本发明涉及一种超声探头坐标定位方法，尤其是涉及一种基于直接线性变换和奇异值分解的磁共振坐标定位方法。

背景技术

治疗超声是利用超声波的穿透性将具有一定能量的超声波传输至人体组织，在超声的机械效应、热效应以及空化效应等多种物理效应作用下，进行肿瘤治疗、溶栓、微泡载药、血脑屏障打开等多种非侵入式治疗，众多医生和科技工作者认为治疗超声技术的广泛推广和应用，将会大力推进无创治疗技术的进展，是人类探究无创治疗技术过程中的重要里程碑。但是单一的治疗超声技术由于存在无法进行精确定位、无法精确监控治疗靶区温度、无法进行实时的疗效评估等难题，使该项技术的发展一度遇到了很大的瓶颈难题。

随着磁共振(Magnetic Resonance Imaging，MRI)技术的迅猛发展，将磁共振技术与治疗超声技术紧密融合，突破了上述一些制约传统治疗超声技术发展的瓶颈，使治疗超声技术在近十年又得到了长足的发展。MRI和治疗超声技术的融合，MRI至少在以下三个方面有重大优势：1、MRI能够精确的定位治疗目标区域，为超声治疗提供实时精确的引导；2、MRI利用特定的对温度敏感的序列，能够准确的测量治疗目标区域的温度，可以实时的监控治疗目标区域的治疗情况，从而实现在超声治疗过程中实时准确的反馈控制热量的加载，既可治疗目标组织，又可以不过度治疗目标组织周围的正常组织，实现真正意义上的适形治疗，保证治疗的安全性和有效性；3、MRI术后疗效评估也更准确。因此MRI和治疗超声技术融合的治疗方案在临床上有重大的应用前景。

在磁共振引导的超声治疗系统中，超声探头，尤其是相控型的超声探头，是否准确对准了治疗目标区域，其定位精度，将直接影响治疗的安全性和有效性。因此，需要研发一种能够对超声探头进行精确定位的方法，以提高超声治疗的有效性。中国专利CN201210541512X公开一种基于磁共振图像的超声探头三维坐标定位方法，该方法采用磁共振成像获取超声探头上定位标记物的图像，并记录该图像的图像信息，根据图像信息计算定位标记物上的定位标记点在MRI坐标系中的三维坐标值，根据定位标记点在MRI坐标系中的三维坐标值和其在超声治疗系统坐标系中的三维坐标值计算MRI坐标系与超声治疗系统坐标系的转换关系矩阵，根据转换关系矩阵和磁共振成像获取的超声探头图像对超声探头进行三维坐标定位。但该方法在转换关系矩阵的求解上还存在不足。

发明内容

本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种定位精度高的基于直接线性变换和奇异值分解的磁共振坐标定位方法。

本发明的目的可以通过以下技术方案来实现：

一种基于直接线性变换和奇异值分解的磁共振坐标定位方法，该方法在超声探头的边沿分布有六条用于定位的定位标记物，每个定位标记物的上下两个端点作为定位点，根据各定位点的位置计算MRI坐标系与超声治疗系统坐标系的转换关系矩阵，实现超声探头三维坐标定位，所述MRI坐标系与超声治疗系统坐标系的转换关系矩阵的计算过程具体为：

1)通过超声探头的横向扫描获取各定位标记物的线方程，并剔除误差最大的一个线方程；

2)通过超声探头的纵向扫描获得超声探头上下表面的面方程；

3)将步骤1)得到的5个线方程分别与步骤2)获得的2个面方程求交，获得10个定位点的坐标；

4)根据获得的定位点坐标信息，利用直接线性变换和奇异值分解算法计算从MRI坐标系到超声治疗系统坐标系的中间转换矩阵H；

5)在步骤3)获得的10个定位点中，依次去掉第i个点，用剩下的9个定位点的坐标信息重复步骤4)获得对应的中间转换矩阵Hi，i＝1,2,…,10；

6)比较步骤4)获得的H和步骤5)获得的Hi剔除其中误差最大的两个定位点，基于剩余的8个定位点的坐标信息，利用直接线性变换和奇异值分解算法计算从MRI坐标系到超声治疗系统坐标系的转换矩阵T。

所述定位标记物为分布在超声探头边沿的定位槽或设于超声探头边沿且垂直于超声探头表面的定位圆柱体。

所述步骤1)具体为：

通过两次横向扫描获得两组横截面上各6个点的坐标，由每个定位标记物的两个点确定各定位标记物的线方程，在获得的6个线方程中，剔除由上下两个横截面中距离最近的两个点所确定的线方程。

所述步骤2)具体为：

通过两次纵向扫描获得两组纵截面，获得每组纵截面上左上角、左下角、右上角、右下角上的四个点的坐标，利用这八个点的坐标计算超声探头上下表面的面方程：