CN111150461A - 一种穿刺针定位系统及方法 - Google Patents

一种穿刺针定位系统及方法 Download PDF

Info

Publication number
CN111150461A
CN111150461A CN202010079873.1A CN202010079873A CN111150461A CN 111150461 A CN111150461 A CN 111150461A CN 202010079873 A CN202010079873 A CN 202010079873A CN 111150461 A CN111150461 A CN 111150461A
Authority
CN
China
Prior art keywords
virtual
puncture needle
probe
positioning
coordinate system
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.)
Pending
Application number
CN202010079873.1A
Other languages
English (en)
Inventor
赵天力
罗衡
肖锦涛
胡世军
柯建源
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to CN202010079873.1A priority Critical patent/CN111150461A/zh
Publication of CN111150461A publication Critical patent/CN111150461A/zh
Priority to EP20917786.4A priority patent/EP3973896A4/en
Priority to PCT/CN2020/091532 priority patent/WO2021155649A1/zh
Priority to US17/334,951 priority patent/US11980496B2/en
Pending legal-status Critical Current

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B17/00—Surgical instruments, devices or methods, e.g. tourniquets
    • A61B17/34—Trocars; Puncturing needles
    • A61B17/3403—Needle locating or guiding means
    • 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
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/42—Details of probe positioning or probe attachment to the patient
    • A61B8/4245—Details of probe positioning or probe attachment to the patient involving determining the position of the probe, e.g. with respect to an external reference frame or to the patient
    • A61B8/4254—Details of probe positioning or probe attachment to the patient involving determining the position of the probe, e.g. with respect to an external reference frame or to the patient using sensors mounted on the probe
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B17/00—Surgical instruments, devices or methods, e.g. tourniquets
    • A61B17/34—Trocars; Puncturing needles
    • A61B17/3403—Needle locating or guiding means
    • A61B2017/3413—Needle locating or guiding means guided by 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/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/2072—Reference field transducer attached to an instrument or patient

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Surgery (AREA)
  • Engineering & Computer Science (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • Pathology (AREA)
  • Robotics (AREA)
  • Biophysics (AREA)
  • Physics & Mathematics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明涉及手术定位技术领域，特别是涉及一种穿刺针定位系统及方法。所述系统包括：超声单元，包括用于提供病灶处的超声切面的第一探头，所述第一探头上设有用于提供第一探头的坐标信息的多个第一定位装置；穿刺针单元，包括穿刺针，所述穿刺针上设有用于提供穿刺针的坐标信息的多个第二定位装置；处理及显示单元，分别与超声单元、各所述第一定位装置及各所述第二定位装置通信连接。本发明提供的穿刺针定位系统可大大提高经胸壁穿刺的精准度，不仅能有效避免反复穿刺带来的并发症，也大大缩短了手术时穿刺的时间。

Description

一种穿刺针定位系统及方法
技术领域
本发明涉及手术定位技术领域，特别是涉及一种穿刺针定位系统及方法。
背景技术
随着医疗科技的高度发展，介入技术已经成为继药物治疗和外科手术之后的第三种有效临床治疗手段。介入路径是成功实施治疗的关键。理想的介入路径应该距离病变近，便于操控器械，同时创伤小。为了便于控制出血，经股静脉和股动脉是结构性心脏病的主要介入路径。但这种介入路径远离心脏，而且行程曲折，因此不利于心内介入的精确操作，同时也限制了介入治疗的适应症。近十年来，经心表路径心内介入技术已广泛应用于临床。由于介入距离心脏更近和超声引导下精准的可视化操作，此项技术不仅拓宽了手术的适应症，而且在某些疾病例如室间隔缺损治疗上取得了比传统介入方法更好的效果。因此，经心表路径应该是结构性心脏病最理想的介入通道。经心表路径实施介入治疗的方法有切开胸壁和经胸壁穿刺两种。无疑，经胸壁穿刺更具优势，不仅创伤更小，而且没有手术切口疤痕。但此项技术关键之处在于选择精准的穿刺位点和路径，确保一次性成功。若反复多次穿刺心脏，会导致因穿刺点出血而产生的心包填塞等严重并发症。为此，术前需要采用经胸超声规划穿刺路径。但在临床实际操作过程中，术者只能凭借个人经验在脑海里重构穿刺针的空间位置，在穿刺点及角度上，存在很大的错选率。因此，此项技术难以在临床上推广应用。如果想要实现经胸壁穿刺介入治疗，一种高度精确的导航系统成为了迫切的需求。
发明内容
鉴于以上所述现有技术的缺点，本发明要解决的技术问题在于提供一种穿刺针定位系统及方法，用于提供精准的穿刺位点和路径，确保一次性穿刺成功。
为实现上述目的及其他相关目的，本发明一方面提供一种穿刺针定位系统，所述系统包括：
超声单元，包括用于提供病灶处的超声切面的第一探头，所述第一探头上设有用于提供第一探头的坐标信息的多个第一定位装置；
穿刺针单元，包括穿刺针，所述穿刺针上设有用于提供穿刺针的坐标信息的多个第二定位装置；
处理及显示单元，分别与超声单元、各所述第一定位装置及各所述第二定位装置通信连接，用于：获取病灶处的超声切面并在虚拟的坐标系统中显示；获取第一探头在获取病灶处的超声切面时的坐标信息并在虚拟的坐标系统中显示；在虚拟坐标系统中确定规划路径；获取各第二定位装置提供的穿刺针的坐标信息并在虚拟的坐标系统中显示虚拟穿刺针中轴线和虚拟穿刺针顶点；在虚拟的坐标系统中，比较虚拟穿刺针中轴线与规划路径是否重合；
本发明另一方面提供一种穿刺针定位方法，包括如下步骤：
S100获取病灶处的超声切面并在虚拟的坐标系统中显示；获取第一探头在获取病灶处的超声切面时的坐标信息并在虚拟的坐标系统中显示；
S200在虚拟的坐标系统中确定规划路径；
S300获取穿刺针上的第二定位装置在虚拟坐标系统中的坐标信息并显示虚拟穿刺针中轴线和虚拟穿刺针顶点；
S400在虚拟的坐标系统中，比较虚拟穿刺针中轴线和所规划路径是否重合。
本发明另一方面提供一种设备，所述设备包括：存储器及处理器；所述存储器，其上存储有计算机程序；所述处理器，用于执行所述存储器存储的计算机程序，该程序被执行时实现本发明所述的穿刺针定位方法。
本发明另一方面提供一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现本发明所述的穿刺针定位方法。
本发明达到了以下有益效果：
本发明提供的穿刺针定位系统可大大提高了经胸壁穿刺的精准度，不仅能有效避免反复穿刺带来的并发症，也大大缩短了手术时穿刺的时间。引入第一定位装置和第二定位装置，将以往单一的超声定位导航予以了极大的丰富，并将所有位置信息予以数字化解析，所有磁场中的位置信息通过第一定位装置和第二定位装置均可转化为三维坐标系。穿刺定位系统依据每个人不同的病情提供了精确的穿刺路径规划，真正实现了对病人的个体化精准医疗，将使手术相关并发症大幅度减少，使经胸穿刺介入治疗变成常规治疗手段成为可能。
附图说明
图1为本发明的穿刺针定位系统应用场景；
图2为本发明第一探头的结构示意图。
图3为本发明超声探头卡扣结构示意图。
图4为本发明穿刺针及穿刺针卡扣件结构示意图。
图5为本发明穿刺针定位方法的流程示意图。
图6为本发明第一探头中轴线投射至超声切面上的坐标转换方法的流程示意图。
图7为本发明穿刺针定位系统可视化界面示意图。
图8为本发明虚拟穿刺针和规划路径重合判定流程示意图。
图中元件标号
1 超声单元
11 第一探头
12 第二探头
2 穿刺针单元
3 处理及显示单元
4 第一定位装置
5 第二定位装置
6 探头卡扣件
7 穿刺针卡扣件
具体实施方式
以下通过特定的具体实例说明本发明的实施方式，本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用，本说明书中的各项细节也可以基于不同观点与应用，在没有背离本发明的精神下进行各种修饰或改变。需说明的是，在不冲突的情况下，以下实施例及实施例中的特征可以相互组合。
需要说明的是，以下实施例中所提供的图示仅以示意方式说明本发明的基本构想，遂图式中仅显示与本发明中有关的组件而非按照实际实施时的组件数目、形状及尺寸绘制，其实际实施时各组件的型态、数量及比例可为一种随意的改变，且其组件布局型态也可能更为复杂。
如图1所示，本发明实施例提供一种穿刺针定位系统。如图1所示场景中，所述系统包括：超声单元1、穿刺针单元2、处理及显示单元3。所述超声单元1包括用于提供病灶处超声切面的第一探头11，所述第一探头11上如图2所示，设有用于提供第一探头11的坐标信息的多个第一定位装置4；所述穿刺针单元2包括穿刺针，所述穿刺针上如图4所示，设有用于提供穿刺针的坐标信息的多个第二定位装置5；所述处理及显示单元3分别与超声单元1、各所述第一定位装置4及各所述第二定位装置5通信连接，用于：获取病灶处的超声切面并在虚拟的坐标系统中显示；获取第一探头11在获取病灶处的超声切面时的坐标信息并在虚拟的坐标系统中显示，进一步在虚拟坐标系统中确定规划路径L1；获取穿刺针上的第二定位装置5在虚拟坐标系统中的坐标信息并显示虚拟穿刺针中轴线L2和虚拟穿刺针顶点C2；在虚拟的坐标系统中，比较虚拟穿刺针中轴线L2与规划路径L1是否重合。当所述虚拟穿刺针和中轴线L2与规划路径L1重合时采用穿刺针进行穿刺。
在所述超声单元1中，超声单元1可以是带多个探头的超声机。如图2所示，超声单元1包括用于提供病灶处的超声切面的第一探头11，所述第一探头11为经胸心脏超声探头。所述第一探头可采用现有常规的超声探头，例如飞利浦的EPIQ7C。使用时，将第一探头11在人体的肋间隙扫查病灶处，例如室间隔缺损，进而获得病灶处的超声切面，通常情况下，当病灶处出现在超声切面上时则为病灶处的超声切面。
所述第一探头11上设有用于提供第一探头11坐标信息的多个第一定位装置4。所述第一探头11的坐标信息包括第一探头11上位若干位置点的坐标信息。所述位置点应选择便于确定第一探头11中轴线的位置点。在一具体实施例中，如图2所示，所述第一探头11上设有3个第一定位装置4，3个所述第一定位装置4位于同一水平面且与所述第一探头11的中轴面重合，其中两个第一定位装置4位于第一探头11的同一横截面上，所述横截面是垂直中轴面的截面，且3个所述第一定位装置4呈直角三角形分布。该设置是为了通过各第一定位装置4进一步确定第一探头11的中轴线。所述第一定位装置4选自传感器。通常情况下，传感器与磁源和磁定位仪配套使用，传感器与磁定位仪通信连接，磁源与磁定位仪通信连接。使用时，磁源设于手术台附近。在一具体实施例中，可以采用加拿大NDI公司的3DGuidancetrakSTAR仪器，所述3DGuidance trakSTAR仪器包括本申请使用的传感器、磁源和磁定位仪。进一步地，为了第一探头11和第一定位装置4的安装稳定性，在所述第一探头11上设有与第一探头11相配合的探头卡扣件6，所述探头卡扣件6如图3所示，其上设有用于安装第一定位装置4的3个第一插孔，3个第一插孔位置与传感器Q1、Q2和Q3匹配。由3个传感器Q1、Q2和Q3，可以获取并记录下此时传感器Q1、Q2和Q3的坐标信息，作为第一探头11的坐标信息。由于三个传感器是设在第一探头11的多个位置点上，即可以知道第一探头11的坐标信息。
所述超声单元1还包括用于监测整个穿刺过程的第二探头12。第二探头12与处理及显示单元3通信连接，可将第二探头扫描到的图像信息在处理及显示单元3中显示。所述第二探头12为经食道超声，所述第二探头12可采用现有常规的超声探头，例如，所述超声探头可以是飞利浦的X7-2T。第二探头12可以为穿刺过程提供实时监测和评估，进一步确保了心内介入操作的安全性。
在所述穿刺针单元2中包括穿刺针，所述穿刺针上设有用于提供穿刺针的坐标信息的多个第二定位装置5。所述穿刺针的坐标信息包括穿刺针上若干个位置点的坐标信息。所述位置点应选择便于确定穿刺针中轴线的位置点。在一具体实施例中，所述穿刺针上设有2个第二定位装置5，所述第二定位装置5分别设于穿刺针的中轴线上，2个第二定位装置5设于穿刺针中轴线的不同位置，该设置是为了通过第二定位装置5进一步确定穿刺针的中轴线。所述第二定位装置5选自传感器。通常情况下，传感器与磁源和磁定位仪配套使用，磁源和传感器分别与磁定位仪通信连接，使用时，磁源设于手术台附近。在一具体实施例中，可以采用加拿大NDI公司的3DGuidance trakSTAR仪器，所述3DGuidance trakSTAR仪器包括本申请使用的传感器、磁源和磁定位仪。进一步地，为了穿刺针和第二定位装置5的安装稳定性，所述穿刺针上设有与穿刺针配合的穿刺针卡扣件7，所述穿刺针卡扣件7如图4所示，其上设有用于安装第二定位装置5的2个第二插孔，2个第二插孔位置与传感器B1和B2匹配。上述可知，由2个传感器B1和B2，可以获取并记录下此时B1和B2的坐标信息，作为穿刺针的坐标信息。由于2个传感器B1和B2是设在穿刺针的多个位置点上，即可以知道穿刺针的坐标信息。
所述处理及显示单元3包括处理器及显示器，所述处理器可以是服务器，还可以是通用处理器，包括中央处理器(Central Processing Unit，简称CPU)、网络处理器(NetworkProcessor，简称NP)等；还可以是数字信号处理器(Digital Signal Processing，简称DSP)、专用集成电路(Application Specific Integrated Circuit，简称ASIC)、现场可编程门阵列(Field－Programmable Gate Array，简称FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。所述处理器与显示器通信连接。所述处理器分别与超声单元1、各所述第一定位装置4及各所述第二定位装置5通信连接。所述处理器与显示器通信连接，用于显示。更具体的，所述处理器与第一探头和第二探头通信连接。
进一步的，处理及显示单元3包括下列模块：
第一数据获取模块：用于获取病灶处的超声切面并在虚拟的坐标系统中显示；获取第一探头11在获取病灶处的超声切面时的坐标信息并在虚拟的坐标系统中显示；
规划路径确定模块：用于在虚拟的坐标系统中确定规划路径L1；
第二数据获取模块：用于获取穿刺针上的第二定位装置5在虚拟坐标系统中的坐标信息并显示虚拟穿刺针中轴线L2和虚拟穿刺针顶点C2；
比较模块：用于在虚拟的坐标系统中，比较虚拟穿刺针中轴线L2与规划路径L1是否重合。
具体的，第一数据获取模块中，病灶处的超声切面通过第一探头11扫查病灶处获得，一般认为，当病灶处出现在超声切面上时则为病灶处的超声切面。
病灶处的超声切面在虚拟坐标系统中的显示位置并无特殊要求。虚拟的坐标系统的原点并无特殊要求，为了便于数据处理，在一较佳的实施方式中，以病灶处超声切面的扇形切面顶点C1作为虚拟坐标系统的原点。
第一探头11可通过第一定位装置4提供其坐标信息。各第一定位装置4的坐标信息即为匹配虚拟坐标系统的坐标信息。可以用各第一定位装置4在虚拟的坐标系统的位置表示第一探头11。为了更方便获得规划路径，需要确定第一探头11中轴线的位置。在一优选的实施方式中，将第一探头11的第一定位装置4均设于第一探头11中轴面上以便于确定第一探头11中轴线的位置。在如图2所示的第一探头11中，可以将3个传感器Q1、Q2和Q3在虚拟的坐标系统的位置表示第一探头11，Q1、Q2和Q3所在的中轴面中，Q1和Q2的垂直平分线即为第一探头11的中轴线。
为了便于后续步骤确定规划路径，在优选的实施方式中，可将预进行穿刺取样的病灶点置于病灶处超声切面的中轴线上。
规划路径确定模块中，可根据虚拟的坐标系统中，第一探头11的空间坐标信息，经平面矩阵转换，得出第一探头11中轴线在病灶处的超声切面上的位置，从而确定规划路径L1。
病灶处超声切面为扇形切面。一般情况下，可将过病灶处超声切面的扇形切面顶点及一病灶穿刺取样点的直线作为实际规划路径。当预进行穿刺取样的病灶点置于病灶处超声切面的中轴线上时，病灶处超声切面的中轴线即可作为实际规划路径。在真实情况下，第一探头11的中轴线应位于病灶处的超声切面上，第一探头11的顶点与病灶处的超声切面顶点重合，第一探头11中轴线与病灶处超声切面的中轴线重合。因此，确定了第一探头11中轴线的位置即可确定规划路径。
但在虚拟坐标系统中，病灶处的超声切面与第一探头11中轴线的相对位置与两者的实际相对位置并不匹配，因此需要进行转换，将第一探头11的中轴线投影至病灶处的超声切面上。
具体的，可采用下列转换公式Ι进行平面矩阵转换，获得第一探头11的中轴线投影至病灶处的超声切面上的转换坐标信息。
(x′，y′，z′，1)＝T(x，y，z，1)T(转换公式Ι)，其中(x′，y′，z′)T转换后的坐标信息，(x，y，z)T为转换前的坐标信息。
其中，T＝T4T3T2T1，T为虚拟第一探头中轴面到其在虚拟超声切面的投影平面的变换矩阵。
Figure BDA0002379914930000071
Figure BDA0002379914930000072
Figure BDA0002379914930000073
Figure BDA0002379914930000074
Q′1＝T3T2T1Q1；
Figure BDA0002379914930000075
Figure BDA0002379914930000076
Figure BDA0002379914930000077
Figure BDA0002379914930000078
Figure BDA0002379914930000079
Figure BDA0002379914930000081
Figure BDA0002379914930000082
Figure BDA0002379914930000083
Figure BDA0002379914930000084
其中，P1、P2、P3为病灶处的虚拟超声切面F1上的三个不共线取样点，Q1、Q2、Q3为虚拟第一探头中轴面F2上的三个不共线的点，
Figure BDA0002379914930000085
Figure BDA0002379914930000086
代表各点的坐标；
Figure BDA0002379914930000087
为虚拟超声切面F1的单位法向量，
Figure BDA0002379914930000088
为虚拟第一探头中轴面F2的单位法向量；
Figure BDA0002379914930000089
及
Figure BDA00023799149300000810
均为将虚拟第一探头中轴面F2平移至Q1与虚拟坐标系统原点重合，获得平面F2'后，再将平面F2'的单位法向量围绕原点旋转至与单位法向量
Figure BDA00023799149300000811
重合时，平面F2'的单位法向量绕X轴逆时针旋转的角度。
Figure BDA00023799149300000812
及
Figure BDA00023799149300000813
均为将虚拟第一探头中轴面F2平移至Q1与虚拟坐标系统原点重合，获得平面F2'后，再将平面F2'的单位法向量围绕原点旋转至与单位法向量
Figure BDA00023799149300000814
重合时，平面F2'的单位法向量绕Y轴逆时针旋转的角度。
将上述第一探头11中轴线投射至超声切面上的平面矩阵转换装置包括：
取样模块，用于获取虚拟坐标系统OXYZ中虚拟超声切面F1上三个不共线取样点坐标
Figure BDA00023799149300000815
以及虚拟第一探头中轴面F2上三个不共线取样点坐标
Figure BDA00023799149300000816
第一平移模块，用于获得平面F2'与虚拟第一探头中轴面F2之间的变换矩阵T1，具体地，平面F2'为平移所述虚拟第一探头中轴面F2直至所述虚拟第一探头中轴面F2上任一取样点与坐标系原点重合获得；
旋转模块，用于获得平面F2”与平面F2'之间的变换矩阵T2，具体地，平面F2”为将所述平面F2'绕坐标原点旋转至其单位法向量与所述虚拟超声切面F1的单位法向量
Figure BDA0002379914930000091
重合获得；
第二平移模块，用于获得平面F2”'与所述平面F2”之间的变换矩阵T3，其中平面F2”'为平移所述平面F2”直至第一平移模块中与坐标系原点重合的取样点回到其初始位置获得。
第三平移模块，用于获得平面F2””与所述平面F2”'之间的变换矩阵T4，其中平面F2””为获取所述平面F2”'与所述虚拟超声切面F1之间的距离，将所述平面F2”'沿其单位法向量(此处所述单位法向量等于F1的单位法向量
Figure BDA0002379914930000092
)平移向量
Figure BDA0002379914930000093
至平面F2”'与虚拟超声切面F1重合获得，利用坐标转换公式及平面F2中对应超声探头中轴线的坐标获取超声切面上超声探头中轴线的坐标。
具体的：
在取样模块中，虚拟超声切面F1是通过第一探头11扫查病灶处获得的病灶处超声切面在虚拟的坐标系统中显示的虚拟切面。虚拟第一探头中轴面F2是第一探头11扫查病灶获得病灶处超声切面时，在虚拟坐标系统中显示的虚拟第一探头中轴面。虚拟第一探头中轴面可依据第一探头11上的传感器获得。
在一优选的实施方式中，虚拟坐标系统还显示病灶处超声切面的扇形切面顶点C1，以便于辅助确定穿刺点。具体的，虚拟超声切面F1为扇形切面，此时真实情况下的第一探头11顶点与胸壁接触的位置在虚拟坐标系统中显示即为病灶处超声切面的扇形切面顶点C1。
在一种实施方式中，可根据第一定位装置4在第一探头11上所设的位置，确定第一探头11的中轴面及中轴线。在一较佳实施方式中，为了便于确定第一探头11的中轴线，第一定位装置4均位于第一探头11的中轴面上。
第一定位装置4可以是传感器，在第一探头11的中轴面上设有如图2所示的三个不共线的传感器Q1、Q2和Q3，用于确定第一探头11的中轴面及中轴线。
为了便于获得中轴线的坐标位置，Q1Q2Q3呈直角三角形，且超声探头的中轴线为线段Q1Q2的中垂线。
在一优选的实施例中，虚拟三维坐标系OXYZ以虚拟超声切面F1的扇形切面顶点C1作为坐标系的原点。
在另一优选的实施例中，第一平移模块中，将虚拟第一探头中轴面F2平移至点Q1与原点重合。
以虚拟三维坐标系OXYZ以虚拟超声切面F1的扇形切面顶点C1作为坐标系的原点，第一平移模块中，将虚拟第一探头中轴面F2平移至点Q1与坐标系原点重合，
Figure BDA0002379914930000101
其中：
Figure BDA0002379914930000102
旋转模块：用于将平面F2'的单位法向量绕坐标系原点旋转至其单位法向量与虚拟超声切面F1的单位法向量
Figure BDA0002379914930000103
重合，获得平面F2”；对应变换矩阵T2为：
Figure BDA0002379914930000104
其中，
Figure BDA0002379914930000105
满足：
Figure BDA0002379914930000106
Figure BDA0002379914930000107
满足：
Figure BDA0002379914930000108
Figure BDA0002379914930000109
满足：
Figure BDA00023799149300001010
Figure BDA00023799149300001011
满足：
Figure BDA00023799149300001012
上述公式解释为：
Figure BDA0002379914930000111
及
Figure BDA0002379914930000112
均为将虚拟第一探头中轴面F2平移至Q1与虚拟坐标系统原点重合，获得平面F2'后，再将平面F2'的单位法向量围绕坐标系原点旋转至与单位法向量
Figure BDA0002379914930000113
重合时，平面F2'的单位法向量绕X轴逆时针旋转的角度。
Figure BDA0002379914930000114
及
Figure BDA0002379914930000115
均为将虚拟第一探头中轴面F2平移至Q1与虚拟坐标系统原点重合，获得平面F2'后，再将平面F2'的单位法向量围绕坐标系原点旋转至与单位法向量
Figure BDA0002379914930000116
重合时，平面F2'的单位法向量绕Y轴逆时针旋转的角度。
通过旋转模块可将平面F2'旋转变换后得到平面F2”。
第二平移模块中，用于平移平面F2”直至第一平移模块中与坐标系原点重合的取样点回到其初始位置获得平面F2”'。以取样点为Q1为例，对应变换矩阵T3为：
Figure BDA0002379914930000117
其中变换矩阵T3:
Figure BDA0002379914930000118
第三平移模块中，用于获取平面F2”'与虚拟超声切面F1之间的距离，设平面F2”'需要平移的向量为
Figure BDA0002379914930000119
其中Q＇1＝T3T2T1Q1。变换矩阵T4为：
Figure BDA0002379914930000121
平面F2”'平移向量
Figure BDA0002379914930000122
后与平面F1重合，获得平面F2””。
综上，得到总变换矩阵为：T＝T4T3T2T1，得到如下结论：对于虚拟第一探头中轴面F2上任一点(x，y，z)T其对应齐次坐标满足：
(x′，y′，z′，1)＝T(x，y，z，1)T，其中(x′，y′，z′)T转换后的坐标信息，(x，y，z)T为转换前的坐标信息。因此可以知道虚拟第一探头中轴面F2与虚拟超声切面F1变换重合后的对应点。
根据上述公式，可知虚拟第一探头中轴面F2上Q1、Q2、Q3与虚拟超声切面F1变换重合后的对应点，进一步可将虚拟第一探头中轴面F2中Q1和Q2的底边垂线，虚拟第一探头中轴面F2上与Q1和Q2的底边中垂线重合的路径即为虚拟超声切面F1中的规划路径L1。
所述第二数据获取模块中，根据虚拟的坐标系统中，第二定位装置5的坐标信息，经平面矩阵转换，获得穿刺针中轴线及穿刺针顶点投射至病灶处的超声切面上的虚拟穿刺针中轴线L2及虚拟穿刺针顶点C2的转换坐标信息，并根据转换坐标信息，在虚拟坐标系统中显示虚拟穿刺针中轴线及其顶点。
具体的，所述虚拟的坐标系统中，以病灶处超声切面的扇形切面顶点C1作为原点时，采用穿刺针定位系统中的转换公式I进行平面矩阵转换，获得对应的虚拟穿刺针中轴线L2及虚拟穿刺针顶点C2的转换坐标信息。
比较模块中，用于在虚拟的坐标系统中，比较虚拟穿刺针顶点与超声切面顶点是否重合，以及虚拟穿刺针中轴线L2与规划路径L1的走向是否重合，当两者均重合时判断为重合。
更具体的，虚拟穿刺针中轴线L2与规划路径L1的重合判断流程如图8所示，在虚拟的坐标系统中，比较虚拟穿刺针顶点C2与病灶处超声切面的扇形切面顶点C1是否重合，若不重合，则判断为与规划路径L1不重合，若重合，则进一步判断虚拟穿刺针中轴线L2与规划路径L1是否重合，若不重合，则判断为与规划路径L1不重合，若重合，则判断为与规划路径L1重合。一般情况下，重合时给出重合提示信号，不重合时给出不重合提示信号。在一具体实施例中，以虚拟穿刺针顶点C2为轴心变换穿刺针的空间位姿，如图7所示，初始状态下，系统指示灯首先显示为红色，当虚拟穿刺针顶点C2和病灶处超声切面的扇形切面顶点C1位置信息重合时系统指示灯为蓝色，当穿刺针的中轴线L2和规划路径L1重合时，系统指示灯显示为绿色，并实时监测显示穿刺针的中轴线L2和规划路径L1的重合信息，在穿刺针行进中根据虚拟界面的标示颜色调整穿刺针，以辅助穿刺顺利进行。以系统指示灯为绿色时的穿刺路径进行穿刺。
基于转换公式，一坐标信息仅能对应唯一的转换坐标信息。
同理，空间中只存在唯一的直线，其坐标经转换后能对应超声切面的中轴线，只有当穿刺针的中轴线与该直线在空间上重合时，其对应的虚拟穿刺针中轴线才能与超声切面的中轴线(即规划路径L1)重合，当穿刺针的中轴线与该直线在空间上不重合时，其对应的虚拟穿刺针中轴线不会与规划路径L1重合。
本发明的穿刺针定位系统的工作过程：
将设有第一定位装置4的第一探头11在患者病灶处扫查，获得包括病灶处的超声切面，并在处理及显示单元3的虚拟坐标系统中显示病灶处的超声切面及第一探头11。病灶处的超声切面为第一探头11扫查到的扇形区域，病灶处超声切面的扇形切面顶点C1对应超声探头与肌肤的接触点，将该接触点作为胸壁穿刺点位置。通过坐标转换公式将第一探头11的中轴面转换成与病灶处的超声切面重合，获得规划路径L1。将穿刺针顶点置于胸壁穿刺点位置，记录各第二定位装置5提供的穿刺针的坐标信息并在虚拟的坐标系统中显示穿刺针的中轴线L2，以胸壁穿刺点位置为轴心来变换穿刺针的空间位姿直至系统提示虚拟穿刺针中轴线L2与规划路径L1位置信息重合，保持穿刺针此时的位姿进行穿刺。在优选的实施方式中，在穿刺针穿刺行进过程中，系统可实时提示穿刺针的位姿是否与规划路径L1重合，并在第二探头12的监测下辅助穿刺顺利进行。
本发明实施例还提供一种穿刺针的定位方法，包括：
S100获取病灶处的超声切面并在虚拟的坐标系统中显示；获取第一探头11在获取病灶处的超声切面时的坐标信息并在虚拟的坐标系统中显示；
S200在虚拟的坐标系统中确定规划路径L1；
S300获取穿刺针上的第二定位装置5在虚拟坐标系统中的坐标信息并显示虚拟穿刺针中轴线L2和虚拟穿刺针顶点C2；
S400在虚拟的坐标系统中，比较虚拟穿刺针中轴线L2和所规划路径L1是否重合。
具体的，步骤S100中，病灶处的超声切面通过第一探头11扫查病灶处获得，一般认为，当病灶处出现在超声切面上时则为病灶处的超声切面。
病灶处的超声切面在虚拟坐标系统中的显示位置并无特殊要求。虚拟的坐标系统的原点并无特殊要求，为了便于数据处理，在一较佳的实施方式中，以病灶处超声切面的扇形切面顶点C1作为虚拟坐标系统的原点。
第一探头11可通过第一定位装置4提供其坐标信息。各第一定位装置4的坐标信息即为匹配虚拟坐标系统的坐标信息。可以用各第一定位装置4在虚拟的坐标系统的位置表示第一探头11。为了更方便获得规划路径，需要确定第一探头11中轴线的位置。在一优选的实施方式中，将第一探头11的第一定位装置4均设于第一探头11中轴面上以便于确定第一探头11中轴线的位置。在如图2所示的第一探头11中，可以将3个传感器Q1、Q2和Q3在虚拟的坐标系统的位置表示第一探头11，Q1、Q2和Q3所在的中轴面中，Q1和Q2的垂直平分线即为第一探头11的中轴线。
为了便于后续步骤确定规划路径，在优选的实施方式中，可将预进行穿刺取样的病灶点置于病灶处超声切面的中轴线上。
步骤S200中，可根据虚拟的坐标系统中，第一探头11的空间坐标信息，经平面矩阵转换，得出第一探头11中轴线在病灶处的超声切面上的位置，从而确定规划路径L1。
病灶处超声切面为扇形切面。一般情况下，可将过病灶处超声切面的扇形切面顶点及一病灶穿刺取样点的直线作为实际规划路径。当预进行穿刺取样的病灶点置于病灶处超声切面的中轴线上时，病灶处超声切面的中轴线即可作为实际规划路径。在真实情况下，第一探头11的中轴线应位于病灶处的超声切面上，第一探头11的顶点与病灶处的超声切面顶点重合，第一探头11中轴线与病灶处超声切面的中轴线重合。因此，确定了第一探头11中轴线的位置即可确定规划路径。
但在虚拟坐标系统中，病灶处的超声切面与第一探头11中轴线的相对位置与两者的实际相对位置并不匹配，因此需要进行转换，将第一探头11的中轴线投影至病灶处的超声切面上。
具体的，可采用下列转换公式Ι进行平面矩阵转换，获得第一探头11的中轴线投影至病灶处的超声切面上的转换坐标信息。
(x′，y′，z′，1)＝T(x，y，z，1)T(转换公式Ι)，其中(x′，y′，z′)T转换后的坐标信息，(x，y，z)T为转换前的坐标信息。
其中，T＝T4T3T2T1，T为虚拟第一探头中轴面到其在虚拟超声切面的投影平面的变换矩阵。
Figure BDA0002379914930000151
Figure BDA0002379914930000152
Figure BDA0002379914930000153
Q′1＝T3T2T1Q1
Figure BDA0002379914930000154
Figure BDA0002379914930000155
其中，
Figure BDA0002379914930000156
满足：
Figure BDA0002379914930000157
Figure BDA0002379914930000158
满足：
Figure BDA0002379914930000159
Figure BDA00023799149300001510
满足：
Figure BDA00023799149300001511
Figure BDA00023799149300001512
满足：
Figure BDA00023799149300001513
Figure BDA00023799149300001514
Figure BDA0002379914930000161
Figure BDA0002379914930000162
其中，P1、P2、P3为病灶处的虚拟超声切面F1上的三个不共线取样点，Q1、Q2、Q3为虚拟第一探头中轴面F2上的三个不共线的点，
Figure BDA0002379914930000163
Figure BDA0002379914930000164
代表各点的坐标；
Figure BDA0002379914930000165
为虚拟超声切面F1的单位法向量，
Figure BDA0002379914930000166
为虚拟第一探头中轴面F2的单位法向量；
Figure BDA0002379914930000167
及
Figure BDA0002379914930000168
均为将虚拟第一探头中轴面F2平移至Q1与虚拟坐标系统原点重合，获得平面F2'后，再将平面F2'的单位法向量围绕坐标系原点旋转至与单位法向量
Figure BDA0002379914930000169
重合时，平面F2'的单位法向量绕X轴逆时针旋转的角度。
Figure BDA00023799149300001610
及
Figure BDA00023799149300001611
均为将虚拟第一探头中轴面F2平移至Q1与虚拟坐标系统原点重合，获得平面F2'后，再将平面F2'的单位法向量围绕坐标系原点旋转至与单位法向量
Figure BDA00023799149300001612
重合时，平面F2'的单位法向量绕Y轴逆时针旋转的角度。
如图6所示，将上述第一探头11中轴线投射至超声切面上的坐标转换公式获得过程如下：
S201获取虚拟坐标系统OXYZ中虚拟超声切面F1上三个不共线取样点坐标
Figure BDA00023799149300001613
以及虚拟第一探头中轴面F2上三个不共线取样点坐标
Figure BDA00023799149300001614
S202获得平面F2'与虚拟第一探头中轴面F2之间的变换矩阵T1，具体地，平面F2'为平移所述虚拟第一探头中轴面F2直至所述虚拟第一探头中轴面F2上任一取样点与坐标系原点重合获得；
S203获得平面F2”与平面F2'之间的变换矩阵T2，具体地，平面F2”为将所述平面F2'的绕坐标系原点旋转至其单位法向量与所述虚拟超声切面F1的单位法向量
Figure BDA00023799149300001615
重合获得；
S204获得平面F2”'与所述平面F2”之间的变换矩阵T3，具体地，平面F2”'为平移所述平面F2”直至步骤S202中与坐标系原点重合的取样点回到其初始位置获得；
S205获得平面F2””与所述平面F2”'之间的变换矩阵T4，其中平面F2””为获取所述平面F2”'与所述虚拟超声切面F1之间的距离，将所述平面F2”'沿其单位法向量(此处所述平面F2”'单位法向量等于F1的单位法向量
Figure BDA0002379914930000171
)平移向量
Figure BDA0002379914930000172
至平面F2”'与虚拟超声切面F1重合获得，利用坐标转换公式及虚拟第一探头中轴面F2中对应第一探头中轴线的坐标获取超声切面上第一探头中轴线的坐标。
具体的，步骤S201中，虚拟超声切面F1是通过第一探头11扫查病灶处获得的病灶处超声切面在虚拟的坐标系统中显示的虚拟切面。虚拟第一探头中轴面F2是第一探头11扫查病灶获得病灶处超声切面时，在虚拟坐标系统中显示的虚拟第一探头中轴面。虚拟第一探头中轴面可依据第一探头11上的传感器获得。
在一优选的实施方式中，虚拟坐标系统还显示病灶处超声切面的扇形切面顶点C1，以便于辅助确定穿刺点。具体的，虚拟超声切面F1为扇形切面，此时真实情况下的第一探头11顶点与胸壁接触的位置在虚拟坐标系统中显示即为病灶处超声切面的扇形切面顶点C1。
在一种实施方式中，可根据第一定位装置4在第一探头11上所设的位置，确定第一探头11的中轴面及中轴线L。在一较佳实施方式中，为了便于确定第一探头11的中轴线，第一定位装置4均位于第一探头11的中轴面上。
第一定位装置4可以是传感器，在第一探头11的中轴面上设有如图2所示的三个不共线的传感器Q1、Q2和Q3，用于确定第一探头11的中轴面及中轴线。
为了便于获得中轴线的坐标位置，Q1Q2Q3呈直角三角形，且第一探头的中轴线为线段Q1Q2的中垂线。
在一优选的实施例中，虚拟三维坐标系OXYZ以虚拟超声切面F1的扇形切面顶点C1作为坐标系的原点。
在另一优选的实施例中，步骤S202中，将虚拟第一探头中轴面F2平移至点Q1与原点重合。
在一具体的实施方式中，以虚拟三维坐标系OXYZ以虚拟超声切面F1的扇形切面顶点C1作为坐标系的原点，步骤S202中，将虚拟第一探头中轴面F2平移至点Q1与坐标系原点重合。
对应坐标变换为：
Figure BDA0002379914930000181
其中变换矩阵T1:
Figure BDA0002379914930000182
步骤S203中，将平面F2'的单位法向量绕坐标系原点旋转至其单位法向量与虚拟超声切面F1的单位法向量
Figure BDA0002379914930000183
重合，获得平面F2”；对应变换为：
Figure BDA0002379914930000184
其中，
Figure BDA0002379914930000185
满足：
Figure BDA0002379914930000186
Figure BDA0002379914930000187
满足：
Figure BDA0002379914930000188
Figure BDA0002379914930000189
满足：
Figure BDA00023799149300001810
Figure BDA00023799149300001811
满足：
Figure BDA00023799149300001812
上述公式解释为：
Figure BDA0002379914930000191
及
Figure BDA0002379914930000192
均为将虚拟第一探头中轴面F2平移至Q1与虚拟坐标系统原点重合，获得平面F2'后，再将平面F2'的单位法向量围绕坐标系原点旋转至与单位法向量
Figure BDA0002379914930000193
重合时，平面F2'的单位法向量绕X轴逆时针旋转的角度。
Figure BDA0002379914930000194
及
Figure BDA0002379914930000195
均为将虚拟第一探头中轴面F2平移至Q1与虚拟坐标系统原点重合，获得平面F2'后，再将平面F2'的单位法向量围绕坐标系原点旋转至与单位法向量
Figure BDA0002379914930000196
重合时，平面F2'的单位法向量绕Y轴逆时针旋转的角度。
通过上述步骤可将平面F2'旋转变换后得到平面F2”。
步骤S204中，平移平面F2”直至步骤S202中与坐标系原点重合的取样点回到其初始位置获得平面F2”'。以取样点为Q1为例，对应变换为：
Figure BDA0002379914930000197
其中变换矩阵T3：
Figure BDA0002379914930000198
步骤S205中，获取平面F2”'与虚拟超声切面F1之间的距离，设平面F2”'需要平移的向量为
Figure BDA0002379914930000199
其中Q′1＝T3T2T1Q1。变换矩阵
Figure BDA00023799149300001910
平面F2”'平移向量
Figure BDA00023799149300001911
后与平面F1重合，获得平面F2””。
综上，得到总变换矩阵为：T＝T4T3T2T1，得到如下结论：对于虚拟第一探头中轴面F2上任一点(x，y，z)T其对应齐次坐标满足：
(x′，y′，z′，1)＝T(x，y，z，1)T，其中(x′，y′，z′)T转换后的坐标信息，(x，y，z)T为转换前的坐标信息。因此可以知道虚拟第一探头中轴面F2与虚拟超声切面F1变换重合后的对应点。
根据上述公式，可知虚拟第一探头中轴面F2上Q1、Q2、Q3与虚拟超声切面F1变换重合后的对应点，进一步可将虚拟第一探头中轴面F2中Q1和Q2的底边垂线，虚拟第一探头中轴面F2上与Q1和Q2的底边中垂线重合的路径即为虚拟超声切面F1中的规划路径L1。
步骤S300中，根据根据虚拟的坐标系统中，第二定位装置5的坐标信息，经平面矩阵转换，获得穿刺针中轴线及穿刺针顶点投射至病灶处的超声切面上的虚拟穿刺针中轴线L2及虚拟穿刺针顶点C2的转换坐标信息，并根据转换坐标信息，在虚拟坐标系统中显示虚拟穿刺针中轴线及其顶点。
具体的，所述虚拟的坐标系统中，以病灶处超声切面的扇形切面顶点C1作为原点时，采用本发明穿刺针定位方法中的转换公式I进行平面矩阵转换。获得对应的虚拟穿刺针中轴线L2及虚拟穿刺针顶点C2的转换坐标信息。
步骤S400中，在虚拟的坐标系统中，比较虚拟穿刺针顶点与超声切面顶点是否重合，以及虚拟穿刺针中轴线L2与规划路径L1的走向是否重合，当两者均重合时判断为重合。
更具体的，虚拟穿刺针中轴线L2与规划路径L1的重合判断流程如图8所示，在虚拟的坐标系统中，比较虚拟穿刺针顶点C2与病灶处超声切面的扇形切面顶点C1是否重合，若不重合，则判断为与规划路径L1不重合，若重合，则进一步判断虚拟穿刺针中轴线L2与规划路径L1是否重合，若不重合，则判断为与规划路径L1不重合，若重合，则判断为与规划路径L1重合。一般情况下，重合时给出重合提示信号，不重合时给出不重合提示信号。在一具体实施例中，以虚拟穿刺针顶点C2为轴心变换穿刺针的空间位姿，如图7所示，初始状态下，系统指示灯首先显示为红色，当虚拟穿刺针顶点C2和病灶处超声切面的扇形切面顶点C1位置信息重合时系统指示灯为蓝色，当穿刺针的中轴线L2和规划路径L1重合时，系统指示灯显示为绿色，并实时监测显示穿刺针的中轴线L2和规划路径L1的重合信息，在穿刺针行进中根据虚拟界面的标示颜色调整穿刺针，以辅助穿刺顺利进行。以系统指示灯为绿色时的穿刺路径进行穿刺。
基于转换公式，一坐标信息仅能对应唯一的转换坐标信息。
同理，空间中只存在唯一的直线，其坐标经转换后能对应超声切面的中轴线，只有当穿刺针的中轴线与该直线在空间上重合时，其对应的虚拟穿刺针中轴线才能与超声切面的中轴线(即规划路径L1)重合，当穿刺针的中轴线与该直线在空间上不重合时，其对应的虚拟穿刺针中轴线不会与规划路径L1重合。
本发明实施例还提供一种设备，所述设备包括：存储器、处理器；所述存储器，其上存储有计算机程序；所述处理器，用于执行所述存储器存储的计算机程序，该程序被执行时实现本发明所述的穿刺针定位方法。
所述存储器可能包含随机存取存储器(Random Access Memory，简称RAM)，也可能还包括非易失性存储器(non-volatile memory)，例如至少一个磁盘存储器。
所述处理器可以是通用处理器，包括中央处理器(Central Processing Unit，简称CPU)、网络处理器(Network Processor，简称NP)等；还可以是数字信号处理器(DigitalSignal Processing，简称DSP)、专用集成电路(Application Specific IntegratedCircuit，简称ASIC)、现场可编程门阵列(Field－Programmable Gate Array，简称FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。
本发明实施例还提供一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现本发明所述的穿刺针定位方法。
所述计算机可读存储介质，本领域普通技术人员可以理解：实现上述各方法实施例的全部或部分步骤可以通过计算机程序相关的硬件来完成。前述的计算机程序可以存储于一计算机可读存储介质中。该程序在执行时，执行包括上述各方法实施例的步骤；而前述的存储介质包括：ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
综上所述，本发明提供的穿刺针定位系统及方法可大大提高了经胸壁穿刺的精准度，不仅能有效避免反复穿刺带来的并发症，也大大缩短了手术时穿刺的时间。引入第一定位装置和第二定位装置，将以往单一的超声定位导航予以了极大的丰富，并将所有位置信息予以数字化解析，所有磁场中的位置信息通过第一定位装置和第二定位装置均可转化为三维坐标系。穿刺定位系统依据每个人不同的病情提供了精确的穿刺路径规划，真正实现了对病人的个体化精准医疗，将使手术相关并发症大幅度减少，使经胸穿刺介入治疗变成常规治疗手段成为可能。
综上所述，本发明有效克服了现有技术中的种种缺点而具高度产业利用价值。
上述实施例仅例示性说明本发明的原理及其功效，而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下，对上述实施例进行修饰或改变。因此，举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变，仍应由本发明的权利要求所涵盖。

Claims (21)

1.一种穿刺针定位系统，其特征在于，所述系统包括：
超声单元(1)，包括用于提供病灶处的超声切面的第一探头(11)，所述第一探头(11)上设有用于提供第一探头(11)的坐标信息的多个第一定位装置(4)；
穿刺针单元(2)，包括穿刺针，所述穿刺针上设有用于提供穿刺针的坐标信息的多个第二定位装置(5)；
处理及显示单元(3)，分别与所述超声单元(1)、各所述第一定位装置(4)及各所述第二定位装置(5)通信连接，用于：获取病灶处的超声切面并在虚拟的坐标系统中显示；获取第一探头(11)在获取病灶处的超声切面时的坐标信息并在虚拟的坐标系统中显示；
在虚拟坐标系统中确定规划路径L1；获取各第二定位装置(5)提供的穿刺针的坐标信息并在虚拟的坐标系统中显示虚拟穿刺针中轴线L2和虚拟穿刺针顶点C2；在虚拟的坐标系统中，比较虚拟穿刺针中轴线L2与规划路径L1是否重合。
2.如权利要求1所述的穿刺针定位系统，其特征在于，所述第一探头(11)上设有3个第一定位装置(4)，3个所述第一定位装置(4)位于同一水平面且与所述第一探头(11)的中轴面重合，其中两个第一定位装置(4)位于第一探头(11)的同一横截面上，且3个所述第一定位装置(4)呈直角三角形分布。
3.如权利要求2所述的穿刺针定位系统，其特征在于，所述第一探头(11)上设有与第一探头(11)相配合的探头卡扣件(6)，所述探头卡扣件(6)上设有用于安装第一定位装置(4)的3个第一插孔，所述3个第一插孔位置与3个第一定位装置(4)匹配。
4.如权利要求1所述的穿刺针定位系统，其特征在于，所述穿刺针上设有2个第二定位装置(5)，所述第二定位装置(5)分别设于穿刺针的中轴线上。
5.如权利要求4所述的穿刺针定位系统，其特征在于，所述穿刺针上设有与穿刺针配合的穿刺针卡扣件(7)，所述穿刺针卡扣件(7)上设有用于安装第二定位装置(5)的2个第二插孔，所述2个第二插孔位置与2个第二定位装置(5)匹配。
6.如权利要求1所述的穿刺针定位系统，其特征在于，所述第一定位装置(4)和所述第二定位装置(5)均选自传感器。
7.如权利要求1所述的穿刺针定位系统，其特征在于，所述超声单元(1)还包括用于监测整个穿刺过程的第二探头(12)。
8.如权利要求1所述的穿刺针定位系统，其特征在于，所述处理及显示单元(3)包括下列模块：
第一数据获取模块：用于获取病灶处的超声切面并在虚拟的坐标系统中显示；获取第一探头(11)在获取病灶处的超声切面时的坐标信息并在虚拟的坐标系统中显示；
规划路径确定模块：用于在虚拟的坐标系统中确定规划路径L1；
第二数据获取模块：用于获取穿刺针上的第二定位装置(5)在虚拟坐标系统中的坐标信息并显示虚拟穿刺针中轴线L2和虚拟穿刺针顶点C2；
比较模块：用于在虚拟的坐标系统中，比较所述虚拟穿刺针中轴线L2与所述规划路径L1是否重合。
9.如权利要求8所述的穿刺针定位系统，其特征在于，所述规划路径确定模块中，根据虚拟的坐标系统中，第一探头(11)的空间坐标信息，经平面矩阵转换，得出第一探头(11)中轴线在病灶处的超声切面上的位置，从而确定规划路径L1。
10.如权利要求9所述的穿刺针定位系统，其特征在于，所述虚拟的坐标系统中，以病灶处超声切面的扇形切面顶点C1作为原点，采用下列转换公式I进行平面矩阵转换，获得第一探头(11)的中轴线投影至病灶处的超声切面上的转换坐标信息；
所述转换公式I为：
(x′，y′，z′，1)＝T(x，y，z，1)T
其中，(x′，y′，z′)T转换后的坐标信息，(x，y，z)T为转换前的坐标信息；
T＝T4T3T2T1
Figure FDA0002379914920000021
Figure FDA0002379914920000022
Figure FDA0002379914920000023
Q′1＝T3T2T1Q1；
Figure FDA0002379914920000031
Figure FDA0002379914920000032
Figure FDA0002379914920000033
Figure FDA0002379914920000034
Figure FDA0002379914920000035
Figure FDA0002379914920000036
Figure FDA0002379914920000037
Figure FDA0002379914920000038
Figure FDA0002379914920000039
其中，P1、P2、P3为病灶处的虚拟超声切面F1上的三个不共线取样点，Q1、Q2、Q3为虚拟第一探头中轴面F2上的三个不共线的点，
Figure FDA00023799149200000310
Figure FDA00023799149200000311
代表各点的坐标；
Figure FDA0002379914920000041
为虚拟超声切面F1的单位法向量，
Figure FDA0002379914920000042
为虚拟第一探头中轴面F2的单位法向量。
11.如权利要求8所述的穿刺针定位系统，其特征在于，所述第二数据获取模块中，根据虚拟的坐标系统中，第二定位装置(5)的坐标信息，经平面矩阵转换，获得穿刺针中轴线及穿刺针顶点投射至病灶处的超声切面上的虚拟穿刺针中轴线L2及虚拟穿刺针顶点C2的转换坐标信息，并根据所述转换坐标信息，在虚拟坐标系统中显示虚拟穿刺针中轴线及其顶点。
12.如权利要求11所述的穿刺针定位系统，其特征在于，所述虚拟的坐标系统中，以病灶处超声切面的扇形切面顶点C1作为原点时，采用权利要求10所述的穿刺针定位系统中的转换公式I进行平面矩阵转换。
13.如权利要求8所述的穿刺针定位系统，其特征在于，所述比较模块中，在虚拟的坐标系统中，比较所述虚拟穿刺针顶点C2与所述超声切面顶点C1是否重合，以及所述虚拟穿刺针中轴线L2与所述规划路径L1的走向是否重合，当两者均重合时判断为重合。
14.一种穿刺针定位方法，包括如下步骤：S100获取病灶处的超声切面并在虚拟的坐标系统中显示；获取第一探头(11)在获取病灶处的超声切面时的坐标信息并在虚拟的坐标系统中显示；S200在虚拟的坐标系统中确定规划路径L1；S300获取穿刺针上的第二定位装置(5)在虚拟坐标系统中的坐标信息并显示虚拟穿刺针中轴线L2和虚拟穿刺针顶点C2；S400在虚拟的坐标系统中，比较所述虚拟穿刺针中轴线L2与所述规划路径L1是否重合。
15.如权利要求14所述的穿刺针定位方法，其特征在于，所述步骤S200中，根据虚拟的坐标系统中，第一探头(11)的空间坐标信息，经平面矩阵转换，得出第一探头(11)中轴线在病灶处的超声切面上的位置，从而确定规划路径L1。
16.如权利要求15所述的穿刺针定位方法，其特征在于，所述虚拟的坐标系统中，以病灶处超声切面的扇形切面顶点C1作为原点，采用下列转换公式I进行平面矩阵转换，获得第一探头(11)的中轴线投影至病灶处的超声切面上的转换坐标信息；所述转换公式I为：
(x′，y′，z′，1)＝T(x，y，z，1)T
其中，(x′，y′，z′)T转换后的坐标信息，(x，y，z)T为转换前的坐标信息；
T＝T4T3T2T1
Figure FDA0002379914920000051
Figure FDA0002379914920000052
Figure FDA0002379914920000053
Q′1＝T3T2T1Q1；
Figure FDA0002379914920000054
Figure FDA0002379914920000055
Figure FDA0002379914920000056
Figure FDA0002379914920000057
Figure FDA0002379914920000058
Figure FDA0002379914920000059
Figure FDA00023799149200000510
Figure FDA0002379914920000061
Figure FDA0002379914920000062
其中，P1、P2、P3为病灶处的虚拟超声切面F1上的三个不共线取样点，Q1、Q2、Q3为虚拟第一探头中轴面F2上的三个不共线的点，
Figure FDA0002379914920000063
Figure FDA0002379914920000064
代表各点的坐标；
Figure FDA0002379914920000065
为虚拟超声切面F1的单位法向量，
Figure FDA0002379914920000066
为虚拟第一探头中轴面F2的单位法向量。
17.如权利要求14所述的穿刺针定位方法，其特征在于，所述步骤S300中，根据根据虚拟的坐标系统中，第二定位装置(5)的坐标信息，经平面矩阵转换，获得穿刺针中轴线及穿刺针顶点投射至病灶处的超声切面上的虚拟穿刺针中轴线L2及虚拟穿刺针顶点C2的转换坐标信息，并根据转换坐标信息，在虚拟坐标系统中显示虚拟穿刺针中轴线及其顶点。
18.如权利要求17所述的穿刺针定位方法，其特征在于，所述虚拟的坐标系统中，以病灶处超声切面的扇形切面顶点C1作为原点时，采用权利要求16所述的穿刺针定位系统中的转换公式I进行平面矩阵转换。
19.如权利要求14所述的穿刺针定位方法，其特征在于，所述步骤S400中，在虚拟的坐标系统中，比较所述虚拟穿刺针顶点C2与所述超声切面顶点C1是否重合，以及所述虚拟穿刺针中轴线L2与所述规划路径L1的走向是否重合，当两者均重合时判断为重合。
20.一种设备，其特征在于，所述设备包括：存储器及处理器；所述存储器，其上存储有计算机程序；所述处理器，用于执行所述存储器存储的计算机程序，所述程序被执行时实现如权利要求14～19任一项权利要求所述的方法。
21.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述程序被处理器执行时实现如权利要求14～19任一项权利要求所述的方法。
CN202010079873.1A 2020-02-04 2020-02-04 一种穿刺针定位系统及方法 Pending CN111150461A (zh)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN202010079873.1A CN111150461A (zh) 2020-02-04 2020-02-04 一种穿刺针定位系统及方法
EP20917786.4A EP3973896A4 (en) 2020-02-04 2020-05-21 PUNCTION NEEDLE POSITIONING SYSTEM AND METHOD
PCT/CN2020/091532 WO2021155649A1 (zh) 2020-02-04 2020-05-21 一种穿刺针定位系统及方法
US17/334,951 US11980496B2 (en) 2020-02-04 2021-05-31 Puncture needle positioning system and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010079873.1A CN111150461A (zh) 2020-02-04 2020-02-04 一种穿刺针定位系统及方法

Publications (1)

Publication Number Publication Date
CN111150461A true CN111150461A (zh) 2020-05-15

Family

ID=70565201

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010079873.1A Pending CN111150461A (zh) 2020-02-04 2020-02-04 一种穿刺针定位系统及方法

Country Status (1)

Country Link
CN (1) CN111150461A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111026269A (zh) * 2019-12-04 2020-04-17 上海褚信医学科技有限公司 基于力反馈的生物组织结构的触觉反馈方法、装置和设备
WO2021155649A1 (zh) * 2020-02-04 2021-08-12 赵天力 一种穿刺针定位系统及方法
WO2023116185A1 (zh) * 2021-12-21 2023-06-29 诺创智能医疗科技(杭州)有限公司 路径确定方法、电子装置及计算机可读存储介质
CN117860381B (zh) * 2024-03-13 2024-05-17 福建省立医院 一种穿刺路径动态规划的系统

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111026269A (zh) * 2019-12-04 2020-04-17 上海褚信医学科技有限公司 基于力反馈的生物组织结构的触觉反馈方法、装置和设备
CN111026269B (zh) * 2019-12-04 2023-10-10 上海褚信医学科技有限公司 基于力反馈的生物组织结构的触觉反馈方法、装置和设备
WO2021155649A1 (zh) * 2020-02-04 2021-08-12 赵天力 一种穿刺针定位系统及方法
WO2023116185A1 (zh) * 2021-12-21 2023-06-29 诺创智能医疗科技(杭州)有限公司 路径确定方法、电子装置及计算机可读存储介质
CN117860381B (zh) * 2024-03-13 2024-05-17 福建省立医院 一种穿刺路径动态规划的系统

Similar Documents

Publication Publication Date Title
US11769292B2 (en) Treatment procedure planning system and method
CN111150461A (zh) 一种穿刺针定位系统及方法
US10258413B2 (en) Human organ movement monitoring method, surgical navigation system and computer readable medium
EP3422297B1 (en) System and method for glass state view in real-time three-dimensional (3d) cardiac imaging
JP2019517291A (ja) 内視鏡画像及び超音波画像の画像ベースの融合
JP7111680B2 (ja) デバイス対画像レジストレーションアルゴリズムからの結果の可視化及び操作
US20070244369A1 (en) Medical Imaging System for Mapping a Structure in a Patient's Body
WO2021155649A1 (zh) 一种穿刺针定位系统及方法
US20200113636A1 (en) Robotically-assisted surgical device, robotically-assisted surgery method, and system
JP7221190B2 (ja) 最適化されたデバイス対画像の位置合わせのための構造マスキングまたはアンマスキング
CN113679470A (zh) 一种用于颅脑穿刺手术的计算机辅助穿刺路径规划方法、装置及存储介质
US11771508B2 (en) Robotically-assisted surgical device, robotically-assisted surgery method, and system
CN111134843B (zh) 一种在超声切面上获得超声探头中轴线的方法及装置
CN212281548U (zh) 一种定位设备
CN116077152A (zh) 穿刺路径规划方法及相关产品
CN209996471U (zh) 一种基于ct三维模型映射混合定位导航系统
US20240090866A1 (en) System and method for displaying ablation zone progression
JP7495216B2 (ja) 鏡視下手術支援装置、鏡視下手術支援方法、及びプログラム
CN111658141B (zh) 一种胃切除术端口位置导航系统、装置及存储介质
JP2021045341A (ja) 鏡視下手術支援装置、鏡視下手術支援方法、及びプログラム
TW202018300A (zh) 超音波探頭的定位方法及超音波系統
CN111148473A (zh) 用于消融可视化的系统和方法
Adebar A system for intraoperative transrectal ultrasound imaging in robotic-assisted laparoscopic radical prostatectomy
Sadeghi Neshat Real-Time Segmentation of Curved Needles For 3D Ultrasound Image Guidance

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