一种新型高精度鲁棒的胶囊内镜定位方法与流程

1.本发明涉及医疗器械技术领域，具体涉及一种新型高精度鲁棒的胶囊内镜定位方法。


背景技术：

2.胶囊内窥镜系统，是一类具有胶囊形状，内置摄像头、无线信号传输装置等传感器的智能系统。通过口服智能胶囊，借助消化道蠕动，胶囊内镜系统可在消化道完整运动并拍摄图像；医务人员通过外部接收器获取到的图像，诊断消化道病情。胶囊内镜系统是传统电子胃肠镜检查系统的有效补充，尤其是小肠区域，传统电子胃肠镜无法直接检查；并且胶囊内镜系统还具备无痛无创伤、检查方便、不影响患者活动等优点。
3.医生需要通过胶囊内镜在消化道的位置推定判断病灶位置，因此胶囊内镜在消化道的定位问题是胶囊内镜系统需要解决的一大关键技术问题。现有的胶囊内镜定位技术，包括磁定位技术，视觉定位技术等，在实验室理想条件下能取得较高的定位精度。但在实际使用中，一方面，目前的磁定位系统要求体外穿戴的磁传感器设备刚性连接才能保证定位精度，影响人体自由活动，而且刚性连接使体外穿戴的信号接收设备容易发生形变，导致定位精度降低。另一方面目前的磁定位系统，需要事先精确标定计算地磁场、零偏等各种地球环境和磁力计设备相关的固有参数，并且需要保持本地磁场环境稳定，但在实际使用场景中，精确地标定不一定有条件实现，而且本地磁场环境复杂多变，导致定位精度低和定位结果不鲁棒。


技术实现要素：

4.针对上述问题，本发明提供一种新型高精度鲁棒的胶囊内镜定位方法，解决目前胶囊内镜定位系统定位精度低和定位结果不鲁棒的问题。
5.本发明采用下述的技术方案：
6.一种新型高精度鲁棒的胶囊内镜定位方法，其特征在于，包括以下步骤：
7.s1、磁力计参数初始标定：
8.s101、将体外磁力计阵列上电，并作8字形运动，记录磁力计读数，通过极小化∑i||m
i-bm||，其中mi是第i个测力计观测值，计算磁力计零偏bm；
9.s102、将体外磁力计阵列静止放置并上电，体内胶囊含有永磁铁，对永磁铁不同位置，记录磁力计读数，构造非线性最小二乘问题，并使用列文伯格-马夸尔特迭代方法求解，计算不同磁力计之间的相对位置和方向参数。
10.s2、定位过程中，磁力计数据预处理：通过插值将磁力计数据统一到相同时间系统，相同采样频率并且时间戳对齐。
11.s3、定位解算体内胶囊的位置p和方向n：
12.s301、基于磁偶极子模型建模永磁铁磁场，磁偶极子模型：
13.14.δp＝p
s-pm15.ρ＝||δp||
[0016][0017]
式中，m为依赖磁体材料的常数；b为位于参考坐标系{nm,pm}的磁棒在ps处产生的磁场；nm为永磁铁方向；pm为永磁铁位置；
[0018]
通过李群扰动方向矢量nm得到扰动模型：
[0019]
nm ∈3×1＝exp
so3
(t3×2δ2×1)nm[0020]
式中，t3×2为nm的正交平面的正交基；∈3×1为方向矢量对应李群的扰动项，exp
so3
为是so3群的指数变换，δ2×1为方向矢量对应李代数的扰动项；
[0021]
s302、建立磁力计测量模型：
[0022][0023]
式中，为磁力计坐标系s到定位参考坐标系w的旋转；ms为测力计量；bs为磁力计零偏；为地磁场在参考坐标系w的值；为本地缓慢变化的磁干扰在参考坐标系w的值；noise为磁力计观测的白噪声；
[0024]
s303、通过磁力计观察模型差分技术，消除地磁场和本地磁场环境的干扰，磁观测单差方程：
[0025][0026]
式中，分别为磁力计si和sj在参考坐标系w的旋转矩阵；分别为对应的磁力计测量值；分别为对应磁力计的零偏；bi,bj分别为永磁铁产生的磁场在磁力计si和sj处的理论值；noise为白噪声；
[0027]
s304、通过coordinatedescent方法计算胶囊较低精度的3自由度位置和2自由度方向，作为后续联合优化的初始值；
[0028]
s305、使用非线性优化技术，联合优化胶囊位置、方向，以及体外磁力计设备的相对位置和方向，获得鲁棒高精度的胶囊内镜定位结果。
[0029]
本发明的有益效果是：
[0030]
1、通过同时优化定位体内胶囊和体外磁力计设备的位姿，解决目前磁定位系统体外设备需要刚性固定的缺点，提升定位精度和鲁棒性；
[0031]
2、设计了一种新的数学模型和求解方法，在磁力计观测模型中消除地磁场和本地局部变化的磁环境的干扰，保证人体自由活动时，也能达到胶囊内镜毫米级定位精度。
附图说明
[0032]
为了更清楚地说明本发明实施例的技术方案，下面将对实施例的附图作简单地介绍，显而易见地，下面描述中的附图仅仅涉及本发明的一些实施例，而非对本发明的限制。
[0033]
图1为本发明流程示意图；
[0034]
图2为本发明设备示意图。
具体实施方式
[0035]
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例的附图，对本发明实施例的技术方案进行清楚、完整地描述。显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。基于所描述的本发明的实施例，本领域普通技术人员在无需创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0036]
除非另外定义，本公开使用的技术术语或者科学术语应当为本公开所属领域内具有一般技能的人士所理解的通常意义。本公开中使用的“包括”或者“包含”等类似的词语意指出现该词前面的元件或者物件涵盖出现在该词后面列举的元件或者物件及其等同，而不排除其他元件或者物件。“上”、“下”、“左”、“右”等仅用于表示相对位置关系，当被描述对象的绝对位置改变后，则该相对位置关系也可能相应地改变。
[0037]
下面结合附图和实施例对本发明进一步说明。
[0038]
如图1所示，一种新型高精度鲁棒的胶囊内镜定位方法，其特征在于，包括以下步骤：
[0039]
s1、磁力计参数初始标定：
[0040]
s101、将体外磁力计阵列上电，并作8字形运动，记录磁力计读数，通过极小化∑i||m
i-bm||，其中mi是第i个测力计观测值，计算磁力计零偏bm；
[0041]
s102、将体外磁力计阵列静止放置并上电，体内胶囊含有永磁铁，对永磁铁不同位置，记录磁力计读数，构造非线性最小二乘问题，并使用列文伯格-马夸尔特迭代方法求解，计算不同磁力计之间的相对位置和方向参数。
[0042]
s2、定位过程中，磁力计数据预处理：通过插值将磁力计数据统一到相同时间系统，相同采样频率并且时间戳对齐；
[0043]
s3、定位解算体内胶囊的位置p和方向n：
[0044]
s301、基于磁偶极子模型建模永磁铁磁场，磁偶极子模型：
[0045][0046]
δp＝p
s-pm[0047]
ρ＝||δp||
[0048][0049]
式中，m为依赖磁体材料的常数；b为位于参考坐标系{nm,pm}的磁棒在ps处产生的磁场；nm为永磁铁方向；pm为永磁铁位置；
[0050]
通过李群扰动方向矢量nm得到扰动模型：
[0051]
nm ∈3×1＝exp
so3
(t3×2δ2×1)nm[0052]
式中，t3×2为nm的正交平面的正交基；∈3×1为方向矢量对应李群的扰动项，exp
so3
为是so3群的指数变换，δ2×1为方向矢量对应李代数的扰动项；
[0053]
s302、建立磁力计测量模型：
[0054]
[0055]
式中，为磁力计坐标系s到定位参考坐标系w的旋转；ms为测力计量；bs为磁力计零偏；为地磁场在参考坐标系w的值；为本地缓慢变化的磁干扰在参考坐标系w的值；noise为磁力计观测的白噪声；
[0056]
s303、通过磁力计观察模型差分技术，消除地磁场和本地磁场环境的干扰，磁观测单差方程：
[0057][0058]
式中，分别为磁力计si和sj在参考坐标系w的旋转矩阵；分别为对应的磁力计测量值；分别为对应磁力计的零偏；bi,bj分别为永磁铁产生的磁场在磁力计si和sj处的理论值；noise为白噪声；
[0059]
s304、通过coordinate descent方法，在磁力计阵列内部任取一点作为永磁铁位置的初值，固定位置利用非线性最小二乘方法计算永磁铁朝向，然后固定朝向，再利用非线性最小二乘方法计算永磁铁位置，反复迭代直至收敛；获得胶囊较低精度的3自由度位置和2自由度方向，作为后续联合优化的初始值；
[0060]
s305、使用非线性最小二乘方法，通过列文伯格-马夸尔特迭代求解技术，联合优化胶囊位置、方向，以及体外磁力计设备的相对位置和方向，获得鲁棒高精度的胶囊内镜定位结果。
[0061]
图2为本发明胶囊与磁力计阵列位置关系示意图。
[0062]
以上所述，仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制，虽然本发明已以较佳实施例揭露如上，然而并非用以限定本发明，任何熟悉本专业的技术人员，在不脱离本发明技术方案范围内，当可利用上述揭示的技术内容做出些许更动或修饰为等同变化的等效实施例，但凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所做的任何简单修改、等同变化与修饰，均仍属于本发明技术方案的范围内。