CN112472294B - 不同超声波设备探头在电磁导航系统中空间位置获取方法 - Google Patents
不同超声波设备探头在电磁导航系统中空间位置获取方法 Download PDFInfo
- Publication number
- CN112472294B CN112472294B CN202011471546.7A CN202011471546A CN112472294B CN 112472294 B CN112472294 B CN 112472294B CN 202011471546 A CN202011471546 A CN 202011471546A CN 112472294 B CN112472294 B CN 112472294B
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- probe
- ultrasonic image
- calibration template
- ultrasonic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000000523 sample Substances 0.000 title claims abstract description 61
- 238000000034 method Methods 0.000 title claims abstract description 13
- 230000009466 transformation Effects 0.000 claims abstract description 25
- 239000011159 matrix material Substances 0.000 claims abstract description 20
- 238000005457 optimization Methods 0.000 claims abstract description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 abstract description 4
- 238000002604 ultrasonography Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000003902 lesion Effects 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 230000005865 ionizing radiation Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/20—Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
- A61B2034/107—Visualisation of planned trajectories or target regions
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
- A61B2034/108—Computer aided selection or customisation of medical implants or cutting guides
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/20—Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
- A61B2034/2046—Tracking techniques
- A61B2034/2051—Electromagnetic tracking systems
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/20—Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
- A61B2034/2046—Tracking techniques
- A61B2034/2063—Acoustic tracking systems, e.g. using ultrasound
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/20—Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
- A61B2034/2046—Tracking techniques
- A61B2034/2065—Tracking using image or pattern recognition
Landscapes
- Health & Medical Sciences (AREA)
- Surgery (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biomedical Technology (AREA)
- Robotics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明提出一种不同超声波设备探头在电磁导航系统中空间位置获取方法,该方法根据标定需要确定四个坐标系,分别为超声波图像坐标系,记作I,标定模板坐标系,记作M,黑色磁片坐标系,记作B,以及探针坐标系,记作A,利用最小二乘法优化得到最优的超声波图像坐标系到黑色磁片坐标系的空间变换矩阵;计算出该空间变换矩阵即完成了超声波探头的标定,也就实现了不同超声波设备的探头在电磁导航系统中空间位置的获取。上述不同超声波设备探头在电磁导航系统中空间位置获取方法将二维超声波图像与电磁导航系统相结合,可以实现二维超声波图像信息到三维空间的高精度转换,为基于超声波图像的手术导航系统提供依据。
Description
技术领域
本发明涉及标定技术领域,尤其涉及一种不同超声波设备探头在电磁导航系统中空间位置获取方法,该方法的实质就是建立超声波图像坐标系到黑色磁片坐标系的转换关系,即二维图像和三维空间之间的转换关系。
背景技术
与CT、MRI图像相比,超声波图像具有采集时间短、无电离辐射、经济方便、可实时成像等优点,成为临床必不可少的诊断和治疗方法。然而,如果想更准确地了解病灶结构时,传统的二维超声就不能很好地满足临床要求。与传统的二维超声波图像相比,三维超声能显示病灶的三维形态及其与周围组织器官的空间位置关系,其成为二维超声的重要辅助手段。但是目前市场上的三维超声基本都是基于二维面阵探头的三维彩超,其只能为一些复杂声像结构的判断提供辅助信息,并对某些病变的诊断起到二维超声无法替代的作用。但是它并不能像CT或MRI那样提供准确的3D器官分割信息,因此尚不可能完成与3D CT/MRI数据的配准导航工作。
发明内容
为了解决现有技术中存在的问题,本申请提出了一种不同超声波设备探头在电磁导航系统中空间位置获取方法,以便将二维超声波图像与电磁导航系统相结合,可以实现二维超声波图像信息到三维空间的高精度转换,为基于超声波图像的手术导航系统提供依据。
为了实现上述目的,本申请提出了一种不同超声波设备探头在电磁导航系统中空间位置获取方法,包括以下步骤:
步骤1、将磁场发射器和电磁导航系统相连接;
步骤2、根据标定需要确定四个坐标系,分别为超声波图像坐标系,记作I,标定模板坐标系,记作M,黑色磁片坐标系,记作B,以及探针坐标系,记作A,这四个坐标系对应着三个空间变换矩阵,分别记作TIB,TBA,TAM,其中TIB是超声波图像坐标系到黑色磁片坐标系的空间变换矩阵,TBA是黑色磁片坐标系到探针坐标系的空间变换矩阵,TAM是探针坐标系到标定模板坐标系的空间变换矩阵;标定模板包括主体,在所述主体上设有至少八个锥形孔,各锥形孔在标定模板坐标系下的坐标已知,即PM(i)(i=1,2,3…n)是已知的;已知标定模板中的一点,记作PM,用超声波探头扫描该点,能够得到超声波图像中的一点,记作PI,则该点在标定模板坐标系和超声波图像坐标系中的坐标转换关系能够表示为:
PM=TAM*TBA*TIB*PI (1);
步骤3、将所述标定模板放置在水槽中,将黑色磁片固定在超声波探头上;
步骤4、将探针插入标定模板的第一个锥形孔中,能够由电磁导航系统得到探针尖端点的位置信息,即得到该锥形孔在探针坐标系下的坐标,记为PA(1);由计算机捕捉超声波图像上该锥形孔的位置信息,该锥形孔在超声波图像上显示的是一个圆形亮点,即得到该锥形孔在超声波图像坐标系下的坐标,记作PI(1);
步骤5、重复步骤4,得到标定模板上的其余锥形孔在探针坐标系下的坐标以及在超声波图像坐标系下的坐标,将标定模板上的各锥形孔在探针坐标系下的坐标记为PA(i),将标定模板上的各锥形孔在超声波图像坐标系下的坐标记作PI(i),(i=1,2,3…n);
步骤6、根据公式(2),利用最小二乘法优化得到最优的探针坐标系到标定模板坐标系的空间变换矩阵TAM;
步骤7、根据公式(3),利用最小二乘法优化得到最优的超声波图像坐标系到黑色磁片坐标系的空间变换矩阵TIB;
其中探针和黑色磁片都属于电磁导航系统,所以黑色磁片坐标系到探针坐标系的空间变换矩阵TBA是已知的,计算出TIB即完成了超声波探头的标定,也就实现了不同超声波设备的探头在电磁导航系统中空间位置的获取。
在一些实施例中,在所述步骤3中,所述水槽中的水的温度在10摄氏度至50摄氏度。
本申请的该方案的有益效果在于上述不同超声波设备探头在电磁导航系统中空间位置获取方法将二维超声波图像与电磁导航系统相结合,可以实现二维超声波图像信息到三维空间的高精度转换,为基于超声波图像的手术导航系统提供依据。
附图说明
图1示出了实施例中标定模板的结构示意图。
附图标记:1-主体,2-锥形孔。
具体实施方式
下面结合附图对本申请的具体实施方式作进一步的说明。
本申请所涉及的不同超声波设备探头在电磁导航系统中空间位置获取方法包括以下步骤:
步骤1、将磁场发射器和电磁导航系统相连接。
由于电磁导航系统下的黑色磁片能够固定在不同超声波设备的探头上,故可以忽略不同超声波探头的差异性,从而建立起数学模型。
步骤2、根据标定需要确定四个坐标系,分别为超声波图像坐标系,记作I,标定模板坐标系,记作M,黑色磁片坐标系,记作B,以及探针坐标系,记作A,这四个坐标系对应着三个空间变换矩阵,分别记作TIB,TBA,TAM,其中TIB是超声波图像坐标系I到黑色磁片坐标系B的空间变换矩阵,TBA是黑色磁片坐标系B到探针坐标系A的空间变换矩阵,TAM是探针坐标系A到标定模板坐标系M的空间变换矩阵;其中标定模板的结构示意图如图1所示,包括主体1,在所述主体1上设有至少八个锥形孔2,各锥形孔在标定模板坐标系下的坐标已知,即PM(i)(i=1,2,3…n)是已知的;所述标定模板是为了兼容各种超声波设备的超声波图像与电磁导航系统的标定工作而设计的,该标定模板具备标定速度快,体积小,步骤简单,对图像成像质量有较明显提高等优点。
已知标定模板中的一点PM,用超声波探头扫描该点,能够得到超声波图像中的一点PI,则该点在标定模板坐标系M和超声波图像坐标系I中的坐标转换关系能够表示为:
PM=TAM*TBA*TIB*PI (1)。
步骤3、将所述标定模板放置在水槽中,为了尽可能降低水温对声速的影响,在水槽中灌入10摄氏度至50摄氏度的水作为介质,将黑色磁片固定在超声波探头上。
步骤4、将探针插入标定模板的第一个锥形孔中,能够由电磁导航系统得到探针尖端点的位置信息,即得到该锥形孔在探针坐标系下的坐标,记为PA(1);由计算机捕捉超声波图像上该锥形孔的位置信息,该锥形孔在超声波图像上显示的是一个圆形亮点,即得到该锥形孔在超声波图像坐标系下的坐标,记作PI(1)。
步骤5、重复步骤4,得到标定模板上的其余锥形孔在探针坐标系下的坐标以及在超声波图像坐标系下的坐标,将标定模板上的各锥形孔在探针坐标系下的坐标记为PA(i),将标定模板上的各锥形孔在超声波图像坐标系下的坐标记作PI(i),(i=1,2,3…n)。
步骤6、根据公式(2),利用最小二乘法优化得到最优的探针坐标系A到标定模板坐标系M的空间变换矩阵TAM;
步骤7、根据公式(3),利用最小二乘法优化得到最优的超声波图像坐标系I到黑色磁片坐标系B的空间变换矩阵TIB;
其中探针和黑色磁片都属于电磁导航系统,所以黑色磁片坐标系B到探针坐标系A的空间变换矩阵TBA是已知的,计算出TIB即完成了超声波探头的标定,也就实现了不同超声波设备的探头在电磁导航系统中空间位置的获取。
本申请所涉及的不同超声波设备探头在电磁导航系统中空间位置获取方法将二维超声波图像与电磁导航系统相结合,可以实现二维超声波图像信息到三维空间的高精度转换,为基于超声波图像的手术导航系统提供依据。
Claims (1)
1.一种不同超声波设备探头在电磁导航系统中空间位置获取方法,其特征在于:包括以下步骤:
步骤1、将磁场发射器和电磁导航系统相连接;
步骤2、根据标定需要确定四个坐标系,分别为超声波图像坐标系,记作I,标定模板坐标系,记作M,黑色磁片坐标系,记作B,以及探针坐标系,记作A,这四个坐标系对应着三个空间变换矩阵,分别记作TIB,TBA,TAM,其中-T IB 是超声波图像坐标系到黑色磁片坐标系的空间变换矩阵,TBA是黑色磁片坐标系到探针坐标系的空间变换矩阵,TAM是探针坐标系到标定模板坐标系的空间变换矩阵;标定模板包括主体,在所述主体上设有至少八个锥形孔,各锥形孔在标定模板坐标系下的坐标已知,即PM(i)(i=1,2,3…n)是已知的;已知标定模板中的一点,记作PM,用超声波探头扫描该点,能够得到超声波图像中的一点,记作PI,则该点在标定模板坐标系和超声波图像坐标系中的坐标转换关系能够表示为:PM=TAM*TBA*TIB*PI(1);
步骤3、将所述标定模板放置在水槽中,所述水槽中的水的温度在10摄氏度至50摄氏度,将黑色磁片固定在超声波探头上;
步骤4、将探针插入标定模板的第一个锥形孔中,能够由电磁导航系统得到探针尖端点的位置信息,即得到该锥形孔在探针坐标系下的坐标,记为PA(1);由计算机捕捉超声波图像上该锥形孔的位置信息,该锥形孔在超声波图像上显示的是一个圆形亮点,即得到该锥形孔在超声波图像坐标系下的坐标,记作PI(1);
步骤5、重复步骤4,得到标定模板上的其余锥形孔在探针坐标系下的坐标以及在超声波图像坐标系下的坐标,将标定模板上的各锥形孔在探针坐标系下的坐标记为PA(i),将标定模板上的各锥形孔在超声波图像坐标系下的坐标记作PI(i),(i=1,2,3…n);
步骤6、根据公式(2),利用最小二乘法优化得到最优的探针坐标系到标定模板坐标系的空间变换矩阵TAM;
步骤7、根据公式(3),利用最小二乘法优化得到最优的超声波图像坐标系到黑色磁片坐标系的空间变换矩阵TIB;
其中探针和黑色磁片都属于电磁导航系统,所以黑色磁片坐标系到探针坐标系的空间变换矩阵TBA是已知的,计算出TIB即完成了超声波探头的标定,也就实现了不同超声波设备的探头在电磁导航系统中空间位置的获取。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011471546.7A CN112472294B (zh) | 2020-12-15 | 2020-12-15 | 不同超声波设备探头在电磁导航系统中空间位置获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011471546.7A CN112472294B (zh) | 2020-12-15 | 2020-12-15 | 不同超声波设备探头在电磁导航系统中空间位置获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112472294A CN112472294A (zh) | 2021-03-12 |
CN112472294B true CN112472294B (zh) | 2023-01-20 |
Family
ID=74917088
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011471546.7A Active CN112472294B (zh) | 2020-12-15 | 2020-12-15 | 不同超声波设备探头在电磁导航系统中空间位置获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112472294B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114176779B (zh) * | 2021-12-31 | 2023-12-26 | 杭州三坛医疗科技有限公司 | 一种手术机器人导航定位方法及装置 |
CN117224208B (zh) * | 2023-09-11 | 2024-03-12 | 北京索诺普科技有限公司 | 超声引导下的介入磁导航系统及导航方法 |
CN117137628B (zh) * | 2023-09-26 | 2024-08-30 | 磅客策(上海)智能医疗科技有限公司 | 一种基于磁导航的活检针及活检方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103110429A (zh) * | 2012-06-11 | 2013-05-22 | 大连理工大学 | 超声波探头的光学标定方法 |
CN105193445A (zh) * | 2015-09-01 | 2015-12-30 | 中国科学院深圳先进技术研究院 | 一种超声探头标定体模、超声探头标定系统及其标定方法 |
CN106510760A (zh) * | 2016-11-29 | 2017-03-22 | 苏州国科康成医疗科技有限公司 | 基于双截面超声探头的标定系统及标定方法 |
CN107714082A (zh) * | 2017-09-04 | 2018-02-23 | 北京航空航天大学 | 一种基于光学定位的超声探头标定装置及方法 |
WO2018035942A1 (zh) * | 2016-08-23 | 2018-03-01 | 深圳先进技术研究院 | 一种柔性穿刺针针尖自动跟踪装置及方法 |
CN107928705A (zh) * | 2017-12-14 | 2018-04-20 | 暨南大学 | 一种基于电磁定位技术的超声探头标定方法与标定装置 |
CN208447639U (zh) * | 2017-12-14 | 2019-02-01 | 暨南大学 | 一种基于电磁定位技术的超声探头标定装置 |
CN110169823A (zh) * | 2019-04-24 | 2019-08-27 | 艾瑞迈迪科技石家庄有限公司 | 超声探头标定方法、装置、终端及存储介质 |
-
2020
- 2020-12-15 CN CN202011471546.7A patent/CN112472294B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103110429A (zh) * | 2012-06-11 | 2013-05-22 | 大连理工大学 | 超声波探头的光学标定方法 |
CN105193445A (zh) * | 2015-09-01 | 2015-12-30 | 中国科学院深圳先进技术研究院 | 一种超声探头标定体模、超声探头标定系统及其标定方法 |
WO2018035942A1 (zh) * | 2016-08-23 | 2018-03-01 | 深圳先进技术研究院 | 一种柔性穿刺针针尖自动跟踪装置及方法 |
CN106510760A (zh) * | 2016-11-29 | 2017-03-22 | 苏州国科康成医疗科技有限公司 | 基于双截面超声探头的标定系统及标定方法 |
CN107714082A (zh) * | 2017-09-04 | 2018-02-23 | 北京航空航天大学 | 一种基于光学定位的超声探头标定装置及方法 |
CN107928705A (zh) * | 2017-12-14 | 2018-04-20 | 暨南大学 | 一种基于电磁定位技术的超声探头标定方法与标定装置 |
CN208447639U (zh) * | 2017-12-14 | 2019-02-01 | 暨南大学 | 一种基于电磁定位技术的超声探头标定装置 |
CN110169823A (zh) * | 2019-04-24 | 2019-08-27 | 艾瑞迈迪科技石家庄有限公司 | 超声探头标定方法、装置、终端及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN112472294A (zh) | 2021-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112472294B (zh) | 不同超声波设备探头在电磁导航系统中空间位置获取方法 | |
CN110464459B (zh) | 基于ct-mri融合的介入计划导航系统及其导航方法 | |
US9521994B2 (en) | System and method for image guided prostate cancer needle biopsy | |
CN110087550B (zh) | 一种超声图像显示方法、设备及存储介质 | |
Blackall et al. | An image registration approach to automated calibration for freehand 3D ultrasound | |
JP6097452B2 (ja) | 超音波撮像システム及び超音波撮像方法 | |
US9713458B2 (en) | Ultrasonic imaging system with body marker annotations | |
US20100208963A1 (en) | System and method for fusing real-time ultrasound images with pre-acquired medical images | |
CN106137249A (zh) | 在窄视场情况下进行配准用于多模态医学成像融合 | |
CN103635143A (zh) | 超声波诊断装置以及超声波图像处理方法 | |
CN106725595B (zh) | 一种电磁定位与b超一体化探头标定装置及其标定方法 | |
WO2021018101A1 (zh) | 数据处理方法、装置、设备及存储介质 | |
CN113974830A (zh) | 一种用于超声引导甲状腺肿瘤热消融的手术导航系统 | |
WO2023179756A1 (zh) | 一种双平面超声图像规划方法和装置 | |
EP3934539B1 (en) | Methods and systems for acquiring composite 3d ultrasound images | |
CN113194837B (zh) | 用于帧索引和图像复查的系统和方法 | |
JP2017506537A (ja) | 超音波ガイド処置に対するゾーン視覚化 | |
US20110054319A1 (en) | Ultrasound system and method for providing a plurality of slice plane images | |
JP7442718B2 (ja) | 画像アーチファクト除去のためのシステム、方法及び装置 | |
JP7261870B2 (ja) | 超音波画像内のツールを追跡するためのシステム及び方法 | |
CN115005864B (zh) | 基于电磁定位的超声探头标定方法 | |
EP4243695B1 (en) | Methods and systems for tracking a motion of a probe in an ultrasound system | |
CN112472132A (zh) | 一种定位成像区域的装置、方法及医学成像装置 | |
Rong et al. | Geometric calibration of freehand ultrasound system with electromagnetic tracking | |
Punithakumar et al. | Cardiac ultrasound multiview fusion using a multicamera tracking system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |