CN115568942A - 多椎体联合约束的椎弓根螺钉最优轨迹规划方法及装置 - Google Patents

多椎体联合约束的椎弓根螺钉最优轨迹规划方法及装置 Download PDF

Info

Publication number
CN115568942A
CN115568942A CN202211259373.1A CN202211259373A CN115568942A CN 115568942 A CN115568942 A CN 115568942A CN 202211259373 A CN202211259373 A CN 202211259373A CN 115568942 A CN115568942 A CN 115568942A
Authority
CN
China
Prior art keywords
screw
optimal
trajectory
vertebral body
nut
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.)
Granted
Application number
CN202211259373.1A
Other languages
English (en)
Other versions
CN115568942B (zh
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.)
Capital Medical University
Beijing Jishuitan Hospital
Original Assignee
Capital Medical University
Beijing Jishuitan Hospital
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 Capital Medical University, Beijing Jishuitan Hospital filed Critical Capital Medical University
Priority to CN202211259373.1A priority Critical patent/CN115568942B/zh
Publication of CN115568942A publication Critical patent/CN115568942A/zh
Application granted granted Critical
Publication of CN115568942B publication Critical patent/CN115568942B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/10Computer-aided planning, simulation or modelling of surgical operations
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B17/00Surgical instruments, devices or methods
    • A61B17/56Surgical instruments or methods for treatment of bones or joints; Devices specially adapted therefor
    • A61B17/58Surgical instruments or methods for treatment of bones or joints; Devices specially adapted therefor for osteosynthesis, e.g. bone plates, screws or setting implements
    • A61B17/68Internal fixation devices, including fasteners and spinal fixators, even if a part thereof projects from the skin
    • A61B17/70Spinal positioners or stabilisers, e.g. stabilisers comprising fluid filler in an implant
    • A61B17/7074Tools specially adapted for spinal fixation operations other than for bone removal or filler handling
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B17/00Surgical instruments, devices or methods
    • A61B17/56Surgical instruments or methods for treatment of bones or joints; Devices specially adapted therefor
    • A61B17/58Surgical instruments or methods for treatment of bones or joints; Devices specially adapted therefor for osteosynthesis, e.g. bone plates, screws or setting implements
    • A61B17/68Internal fixation devices, including fasteners and spinal fixators, even if a part thereof projects from the skin
    • A61B17/70Spinal positioners or stabilisers, e.g. stabilisers comprising fluid filler in an implant
    • A61B17/7074Tools specially adapted for spinal fixation operations other than for bone removal or filler handling
    • A61B17/7083Tools for guidance or insertion of tethers, rod-to-anchor connectors, rod-to-rod connectors, or longitudinal elements
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B17/00Surgical instruments, devices or methods
    • A61B17/56Surgical instruments or methods for treatment of bones or joints; Devices specially adapted therefor
    • A61B17/58Surgical instruments or methods for treatment of bones or joints; Devices specially adapted therefor for osteosynthesis, e.g. bone plates, screws or setting implements
    • A61B17/68Internal fixation devices, including fasteners and spinal fixators, even if a part thereof projects from the skin
    • A61B17/84Fasteners therefor or fasteners being internal fixation devices
    • A61B17/86Pins or screws or threaded wires; nuts therefor
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B17/00Surgical instruments, devices or methods
    • A61B17/56Surgical instruments or methods for treatment of bones or joints; Devices specially adapted therefor
    • A61B17/58Surgical instruments or methods for treatment of bones or joints; Devices specially adapted therefor for osteosynthesis, e.g. bone plates, screws or setting implements
    • A61B17/88Osteosynthesis instruments; Methods or means for implanting or extracting internal or external fixation devices
    • A61B17/90Guides therefor
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/20Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/10Computer-aided planning, simulation or modelling of surgical operations
    • A61B2034/107Visualisation of planned trajectories or target regions
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/10Computer-aided planning, simulation or modelling of surgical operations
    • A61B2034/108Computer aided selection or customisation of medical implants or cutting guides

Landscapes

  • Health & Medical Sciences (AREA)
  • Orthopedic Medicine & Surgery (AREA)
  • Surgery (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Neurology (AREA)
  • Medical Informatics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biomedical Technology (AREA)
  • Molecular Biology (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Robotics (AREA)
  • Surgical Instruments (AREA)

Abstract

本发明涉及一种多椎体联合约束的椎弓根螺钉最优轨迹规划方法,包括如下步骤:确定每个椎体置钉的安全置钉区域,获取每个椎体上螺钉轨迹的入点范围;利用每个椎体上螺钉的入点范围,计算连接杆的最优形态;确定多轴钉钉帽的活动范围;利用椎体置钉的规避区域、连接杆位置及多轴钉钉帽的活动范围形成临床有效螺钉置入轨迹的约束条件,然后确定出每个椎体在联合约束条件下的椎弓根螺钉最优轨迹。本发明能够准确、快速的计算出椎弓根螺钉的临床有效三维置入路径。

Description

多椎体联合约束的椎弓根螺钉最优轨迹规划方法及装置
技术领域
本发明涉及脊柱椎弓根螺钉置入技术领域,具体是关于一种多椎体联合约束的椎弓根螺钉最优轨迹规划方法及装置。
背景技术
图像引导(机器人辅助)椎弓根螺钉置入可以提高螺钉放置的准确性。与手动螺钉置入相比,图像引导(机器人辅助)螺钉置入具有更高的准确率。然而,目前的商业机器人导航系统普遍缺乏自动化。医生必须在CT或CBCT像上手动规划螺钉轨迹的方向和位置。
目前的螺钉轨迹规划主要是在CBCT图像上使用手动拖放策略来获得螺钉轨迹。此过程会中断工作流程并可能导致错误。因此计算机辅助规划方法成为一种可行的替代方案。例如从不同角度揭示了椎体骨密度和螺钉置入轨迹之间的关系,并根据相关关系规划螺钉轨迹。采用3D规划系统,通过CT图像衍生的骨力学性质或骨密度特性能来寻找和优化螺钉轨迹。根据CT图像和椎体几何特性优化螺钉轨迹,例如通过构建多角度投影并在投影面上标记规避区域的方法寻找最优螺钉轨迹。通过构建椎体和螺钉轨迹图谱然后将图谱配准到目标椎体上进行螺钉轨迹规划。基于椎弓根几何特征的椎体和椎弓根自动分割和轨迹规划方法。通过人工智能方法进行椎体分割和螺钉轨迹规划。
但实践中,通常需要同时处理多个椎体。在这种情况下,同侧螺钉需要与刚性的连接杆连接,以实现最终的目的。因此除了对单个椎体螺钉放置的约束外,还需要很好地理解同侧不同椎体螺钉之间的相互约束关系。但到目前为止,传统的螺钉轨迹规划方法都是针对单椎体的,未考虑多椎体约束的情况。
发明内容
针对上述问题,本发明的目的是提供一种多椎体联合约束的椎弓根螺钉最优轨迹规划方法及装置,该方法为多椎体联合置钉提供确定入钉点范围及连接杆设计所需的有效方法,能够准确、快速的计算出椎弓根螺钉的临床有效三维置入路径
为实现上述目的,本发明采取以下技术方案:
本发明所述的多椎体联合约束的椎弓根螺钉最优轨迹规划方法,包括如下步骤:
确定每个椎体置钉的安全置钉区域,获取每个椎体上螺钉轨迹的入点范围;
利用每个椎体上螺钉的入点范围,计算连接杆的最优形态;
确定多轴钉钉帽的活动范围;
利用椎体置钉的规避区域、连接杆位置及多轴钉钉帽的活动范围形成临床有效螺钉置入轨迹的约束条件,然后确定出每个椎体在联合约束条件下的椎弓根螺钉最优轨迹。
所述的多椎体联合约束的椎弓根螺钉最优轨迹规划方法,优选地,确定每个椎体置钉的安全置钉区域,获取每个椎体上螺钉轨迹的入点范围是通过以下方式实现的:
椎弓根螺钉最优路径为:
Figure BDA0003890807600000021
椎体置钉的安全置钉区域为:
Figure BDA0003890807600000022
其中,M是规避区域上的点,M0为最优螺钉轨迹的入钉点;ρ为螺钉到规避区域的安全距离,且ρ>0;
Figure BDA0003890807600000023
为螺钉轨迹的单位方向向量;
Figure BDA0003890807600000024
为特殊符号表示任意;
对公式(2)进行放射变换,使其简化为一个以最优螺钉轨迹为中心线圆柱体,具体为:
Figure BDA0003890807600000025
z2+Y2=(d-ρ)2 (4)
0≤X≤L
其中,
Figure BDA0003890807600000031
为(α,β,γ)的正交的单位向量;
Figure BDA0003890807600000032
Figure BDA0003890807600000033
为(α,β,γ)的正交的单位向量;
Figure BDA0003890807600000034
为空间变换矩阵;L为螺钉长度;
Figure BDA0003890807600000035
为最优螺钉轨迹到规避区域的最短距离;
Figure BDA0003890807600000036
Figure BDA0003890807600000037
为安全置钉区域;ρ为螺钉到避开区域的安全距离;L为螺钉长度;X、Y、Z为原始坐标;
螺钉轨迹的入钉点范围为:
Figure BDA0003890807600000038
z2+Y2=(d-ρ)2 (6)
其中,
Figure BDA0003890807600000039
为螺钉轨迹的入钉点范围。
所述的多椎体联合约束的椎弓根螺钉最优轨迹规划方法,优选地,所述连接杆的最优形态是符合脊柱自然弯曲形态的曲线,且位于尽可能靠近所有同侧最优螺钉轨迹方向上的平面;所述最优连接杆的形态包括其位置和形状,具体为:
设M0(x,y,z)为所有椎体同侧最优螺钉轨迹入点位置,则最优连接杆位置表示为:
AX+BY+CZ+D=0 (7)
Figure BDA00038908076000000310
因此,在连接杆约束下的最优螺钉轨迹表示为:
Figure BDA00038908076000000311
Figure BDA0003890807600000041
其中,
Figure BDA0003890807600000042
A、B、C、D为连接杆所在位置平面的参数;L′为钉帽长度;(αnut,βnut,γnut)为钉帽的单位方向向量;ω为钉帽能够与钉子形成的最大夹角;
设已知连接杆的最佳平面为P:
AX+BY+CZ+D=0 (11)
以该平面为XOZ面可构建新的空间坐标系{iT,jT,kT},则:
Figure BDA0003890807600000043
Figure BDA0003890807600000044
Figure BDA0003890807600000045
其中,N0、Nx和Nz为最佳平面P上的点,且
Figure BDA0003890807600000046
垂直于
Figure BDA0003890807600000047
则钉帽在最佳平面P上的坐标为:
Figure BDA0003890807600000048
Figure BDA0003890807600000049
其中,
Figure BDA00038908076000000410
为将
Figure BDA00038908076000000411
坐标转换为
Figure BDA00038908076000000412
坐标上的空间变换矩阵;
Figure BDA00038908076000000413
为钉帽在空间坐标系{iT,jT,kT}下的坐标;由于,定义最佳平面P为空间坐标系{iT,jT,kT}的XOZ面,则yT=0,因此,连接杆的形状通过最小二乘法拟合得出:
Figure BDA00038908076000000414
综上,得到连接杆最优位置和形状为:
Figure BDA0003890807600000051
其中,xTl和zTl为钉帽在新坐标系下的坐标;xol、yol和zol为钉帽在原始坐标系下的坐标;
Figure BDA0003890807600000054
为钉帽所在曲线的曲线参数。
所述的多椎体联合约束的椎弓根螺钉最优轨迹规划方法,优选地,确定多轴钉的钉帽的活动范围是通过以下方式实现的:
设螺钉钉帽长度为L’,则钉帽的参数方程表示为:
(x,y,z)=t(αnut,βnut,γnut)+Mnut (19)
s.t.0<t<L′
其中,(αnut,βnut,γnut)为钉帽的单位方向向量,Mnut为钉帽的位置,t表示钉帽长度的未知数;
则钉帽的可活动范围为:
Figure BDA0003890807600000055
Figure BDA0003890807600000052
其中
Figure BDA0003890807600000053
为螺钉的单位方向向量;M0为最优螺钉轨迹的入钉点;ω为钉帽与螺钉之间的最大可活动夹角。
本发明所述的多椎体联合约束的椎弓根螺钉最优轨迹规划装置,包括:
第一处理单元,用于确定每个椎体置钉的安全置钉区域,获取每个椎体上螺钉轨迹的入点范围;
第二处理单元,用于利用每个椎体上螺钉的入点范围,计算连接杆的最优形态;
第三处理单元,用于确定多轴钉钉帽的活动范围;
第四处理单元,用于利用椎体置钉的规避区域、连接杆位置及多轴钉钉帽的活动范围形成临床有效螺钉置入轨迹的约束条件,然后确定出每个椎体在联合约束条件下的椎弓根螺钉最优轨迹。
本发明所述的计算机存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现所述的多椎体联合约束的椎弓根螺钉最优轨迹规划方法步骤。
本发明所述的计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现所述多椎体联合约束的椎弓根螺钉最优轨迹规划方法步骤。
本发明由于采取以上技术方案,其具有以下优点:
本发明为多椎体联合置钉提供确定入钉点范围及连接杆设计所需的有效方法,能够准确、快速的计算出椎弓根螺钉的临床有效三维置入路径。
附图说明
通过阅读下文优选实施方式的详细描述,各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的,而并不认为是对本发明的限制。在整个附图中,用相同的附图标记表示相同的部件。在附图中:
图1是本发明的螺钉的活动范围示意图;
图2是本发明的连接杆的最优形态的示意图。
具体实施方式
下面将参照附图更详细地描述本发明的示例性实施方式。虽然附图中显示了本发明的示例性实施方式,然而应当理解,可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。相反,提供这些实施方式是为了能够更透彻地理解本发明,并且能够将本发明的范围完整的传达给本领域的技术人员。
本发明提供一种多椎体联合约束的椎弓根螺钉最优轨迹规划方法,通过确定每个椎体置钉的安全区域,获取每个椎体上螺钉轨迹的入点范围;利用每个椎体上螺钉的入点范围,计算连接杆的临床有效位置;确定多轴钉钉帽的活动范围;利用椎体置钉的规避区域、连接杆位置及多轴钉钉帽的活动范围形成最优螺钉置入轨迹的约束条件,然后计算出每个椎体在联合约束条件下的有效螺钉轨迹。本发明为多椎体联合置钉提供确定入钉点范围及连接杆设计所需的有效方法,能够准确、快速的计算出椎弓根螺钉的临床有效三维置入路径。
本发明提供的多椎体联合约束的椎弓根螺钉最优轨迹规划方法,包括如下步骤:
S1、确定每个椎体置钉的安全置钉区域,获取每个椎体上螺钉轨迹的入点范围;
S2、利用每个椎体上螺钉的入点范围,计算连接杆的最优形态;
S3、确定多轴钉钉帽的活动范围;
S4、利用椎体置钉的规避区域、连接杆位置及多轴钉钉帽的活动范围形成临床有效螺钉置入轨迹的约束条件,然后确定出每个椎体在联合约束条件下的椎弓根螺钉最优轨迹。
在上述实施例中,优选地,确定每个椎体置钉的安全置钉区域,获取每个椎体上螺钉轨迹的入点范围是通过以下方式实现的:
椎弓根螺钉最优路径为:
Figure BDA0003890807600000071
椎体置钉的安全置钉区域为:
Figure BDA0003890807600000072
其中,M是规避区域上的点,M0为最优螺钉轨迹的入钉点;ρ>0,为螺钉到规避区域的安全距离;
Figure BDA0003890807600000073
为螺钉轨迹的单位方向向量;
Figure BDA0003890807600000074
为特殊符号表示任意;
为置钉的安全性和计算的简便性可将公式(2)进行放射变换,使其简化为一个以最优螺钉轨迹为中心线圆柱体,具体为:
Figure BDA0003890807600000075
z2+Y2=(d-ρ)2 (4)
0≤X≤L
其中,
Figure BDA0003890807600000076
为(α,β,γ)的正交的单位向量;
Figure BDA0003890807600000077
Figure BDA0003890807600000078
为(α,β,γ)的正交的单位向量;
Figure BDA0003890807600000079
为空间变换矩阵;L为螺钉长度;
Figure BDA00038908076000000710
为最优螺钉轨迹到规避区域的最短距离;
Figure BDA00038908076000000711
Figure BDA0003890807600000081
为安全置钉区域;ρ为螺钉到避开区的安全距离;L为螺钉长度;X、Y、Z为原始坐标;
其中,螺钉轨迹的入钉点范围为:
Figure BDA0003890807600000082
z2+Y2=(d-p)2 (6)
其中,
Figure BDA0003890807600000083
为螺钉轨迹的入钉点范围。
在上述实施例中,优选地,所述连接杆的最优形态是符合脊柱自然弯曲形态的曲线,且位于尽可能靠近所有同侧最优螺钉轨迹方向上的平面;所述最优连接杆的形态包括其位置和形状,具体为:
设M0(x,y,z)为所有椎体同侧最优螺钉轨迹入点位置,则最优连接杆位置可表示为:
AX+BY+CZ+D=0 (7)
Figure BDA0003890807600000084
因此,在连接杆约束下的最优螺钉轨迹可表示为:
Figure BDA0003890807600000085
Figure BDA0003890807600000086
其中,
Figure BDA0003890807600000087
A、B、C、D为连接杆所在位置平面的参数;L′为钉帽长度;(αnut,βnut,γnut)为钉帽的单位方向向量;ω钉帽能够与钉子形成的最大夹角;
设已知连接杆的最佳平面为P:
AX+BY+CZ+D=0 (11)
以该平面为XOZ面可构建新的空间坐标系{iT,jT,kT},则:
Figure BDA0003890807600000091
Figure BDA0003890807600000092
Figure BDA0003890807600000093
其中,N0,Nx,Nz为最佳平面P上的点,月
Figure BDA0003890807600000094
垂直于
Figure BDA0003890807600000095
则钉帽在最佳平面P上的坐标为:
Figure BDA0003890807600000096
Figure BDA0003890807600000097
其中,
Figure BDA0003890807600000098
为将
Figure BDA0003890807600000099
坐标转换为
Figure BDA00038908076000000910
坐标上的空间变换矩阵;
Figure BDA00038908076000000911
为钉帽在空间坐标系{iT,jT,kT}下的坐标;由于,定义最佳平面P为空间坐标系{iT,jT,kT}的XOZ面,则yT=0,因此,连接杆的形状可以通过最小二乘法拟合得出:
Figure BDA00038908076000000913
如图2所示,综上可得连接杆最优位置和形状为:
Figure BDA00038908076000000912
其中,xTl和zTl为钉帽在新坐标系下的坐标;xol、yol和zol为钉帽在原始坐标系下的坐标;
Figure BDA0003890807600000103
为钉帽所在曲线的曲线参数。
在上述实施例中,优选地,确定多轴钉的钉帽的活动范围是通过以下方式实现的:
设螺钉钉帽长度为L’,则钉帽的参数方程可表示为:
(x,y,z)=t(αnut,βnut,γnut)+Mnut (19)
s.t.0<t<L′
其中,(αnut,βnut,γnut)为钉帽的单位方向向量,Mnut为钉帽的位置,t表示钉帽长度的未知数;
则,如图1所示,钉帽的可活动范围为:
Figure BDA0003890807600000104
Figure BDA0003890807600000101
其中
Figure BDA0003890807600000102
为螺钉的单位方向向量;M0为最优螺钉轨迹的入钉点;ω为钉帽与螺钉之间的最大可活动夹角。
本发明还提供一种多椎体联合约束的椎弓根螺钉最优轨迹规划装置,包括:
第一处理单元,用于确定每个椎体置钉的安全置钉区域,获取每个椎体上螺钉轨迹的入点范围;
第二处理单元,用于利用每个椎体上螺钉的入点范围,计算连接杆的最优形态;
第三处理单元,用于确定多轴钉钉帽的活动范围;
第四处理单元,用于利用椎体置钉的规避区域、连接杆位置及多轴钉钉帽的活动范围形成临床有效螺钉置入轨迹的约束条件,然后确定出每个椎体在联合约束条件下的椎弓根螺钉最优轨迹。
本发明还提供一种计算机存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现上述的多椎体联合约束的椎弓根螺钉最优轨迹规划方法步骤。
本发明还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述多椎体联合约束的椎弓根螺钉最优轨迹规划方法步骤。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (7)

1.一种多椎体联合约束的椎弓根螺钉最优轨迹规划方法,其特征在于,包括如下步骤:
确定每个椎体置钉的安全置钉区域,获取每个椎体上螺钉轨迹的入点范围;
利用每个椎体上螺钉的入点范围,计算连接杆的最优形态;
确定多轴钉钉帽的活动范围;
利用椎体置钉的规避区域、连接杆位置及多轴钉钉帽的活动范围形成临床有效螺钉置入轨迹的约束条件,然后确定出每个椎体在联合约束条件下的椎弓根螺钉最优轨迹。
2.根据权利要求1所述的多椎体联合约束的椎弓根螺钉最优轨迹规划方法,其特征在于,确定每个椎体置钉的安全置钉区域,获取每个椎体上螺钉轨迹的入点范围是通过以下方式实现的:
椎弓根螺钉最优路径为:
Figure FDA0003890807590000011
椎体置钉的安全置钉区域为:
Figure FDA0003890807590000012
其中,M是规避区域上的点,M0为最优螺钉轨迹的入钉点;ρ为螺钉到规避区域的安全距离,且ρ>0;
Figure FDA0003890807590000013
为螺钉轨迹的单位方向向量;
Figure FDA0003890807590000014
为特殊符号表示任意;
对公式(2)进行放射变换,使其简化为一个以最优螺钉轨迹为中心线圆柱体,具体为:
Figure FDA0003890807590000015
z22=(d-ρ)2 (4)
0≤X≤L
其中,
Figure FDA0003890807590000021
为(α,β,γ)的正交的单位向量;
Figure FDA0003890807590000022
Figure FDA0003890807590000023
为(α,β,γ)的正交的单位向量;
Figure FDA0003890807590000024
为空间变换矩阵;L为螺钉长度;
Figure FDA0003890807590000025
为最优螺钉轨迹到规避区域的最短距离;
Figure FDA0003890807590000026
Figure FDA0003890807590000027
为安全置钉区域;ρ为螺钉到避开区域的安全距离;L为螺钉长度;X、Y、Z为原始坐标;
螺钉轨迹的入钉点范围为:
Figure FDA0003890807590000028
z2+Y2=(d-ρ)2 (6)
其中,
Figure FDA0003890807590000029
为螺钉轨迹的入钉点范围。
3.根据权利要求2所述的多椎体联合约束的椎弓根螺钉最优轨迹规划方法,其特征在于,所述连接杆的最优形态是符合脊柱自然弯曲形态的曲线,且位于尽可能靠近所有同侧最优螺钉轨迹方向上的平面;所述最优连接杆的形态包括其位置和形状,具体为:
设M0(x,y,z)为所有椎体同侧最优螺钉轨迹入点位置,则最优连接杆位置表示为:
AX+BY+CZ+D=0 (7)
Figure FDA00038908075900000210
因此,在连接杆约束下的最优螺钉轨迹表示为:
Figure FDA0003890807590000031
Figure FDA0003890807590000032
其中,
Figure FDA0003890807590000033
A、B、C、D为连接杆所在位置平面的参数;L′为钉帽长度;(αnutnutnut)为钉帽的单位方向向量;ω为钉帽能够与钉子形成的最大夹角;
设已知连接杆的最佳平面为P:
AX+BY+CZ+D=0 (11)
以该平面为XOZ面可构建新的空间坐标系{iT,jT,kT},则:
Figure FDA0003890807590000034
Figure FDA0003890807590000035
Figure FDA0003890807590000036
其中,N0、Nx和Nz为最佳平面P上的点,且
Figure FDA0003890807590000037
垂直于
Figure FDA0003890807590000038
则钉帽在最佳平面P上的坐标为:
Figure FDA0003890807590000039
Figure FDA00038908075900000310
其中,
Figure FDA00038908075900000311
为将
Figure FDA00038908075900000312
坐标转换为
Figure FDA00038908075900000313
坐标上的空间变换矩阵;
Figure FDA00038908075900000314
为钉帽在空间坐标系{iT,jT,kT}下的坐标;由于,定义最佳平面P为空间坐标系{iT,jT,kT}的XOZ面,则yT=0,因此,连接杆的形状通过最小二乘法拟合得出:
Figure FDA00038908075900000315
综上,得到连接杆最优位置和形状为:
Figure FDA0003890807590000041
其中,xTl和zTl为钉帽在新坐标系下的坐标;xol、yol和zol为钉帽在原始坐标系下的坐标;
Figure FDA0003890807590000045
为钉帽所在曲线的曲线参数。
4.根据权利要求1所述的多椎体联合约束的椎弓根螺钉最优轨迹规划方法,其特征在于,确定多轴钉的钉帽的活动范围是通过以下方式实现的:
设螺钉钉帽长度为L’,则钉帽的参数方程表示为:
(x,y,z)=t(αnutnutnut)+Mnut (19)
s.t.0<t<L′
其中,(αnutnutnut)为钉帽的单位方向向量,Mnut为钉帽的位置,t表示钉帽长度的未知数;
则钉帽的可活动范围为:
Figure FDA0003890807590000042
Figure FDA0003890807590000043
其中
Figure FDA0003890807590000044
为螺钉的单位方向向量;M0为最优螺钉轨迹的入钉点;ω为钉帽与螺钉之间的最大可活动夹角。
5.一种多椎体联合约束的椎弓根螺钉最优轨迹规划装置,其特征在于,包括:
第一处理单元,用于确定每个椎体置钉的安全置钉区域,获取每个椎体上螺钉轨迹的入点范围;
第二处理单元,用于利用每个椎体上螺钉的入点范围,计算连接杆的最优形态;
第三处理单元,用于确定多轴钉钉帽的活动范围;
第四处理单元,用于利用椎体置钉的规避区域、连接杆位置及多轴钉钉帽的活动范围形成临床有效螺钉置入轨迹的约束条件,然后确定出每个椎体在联合约束条件下的椎弓根螺钉最优轨迹。
6.一种计算机存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1-4中任一项所述的多椎体联合约束的椎弓根螺钉最优轨迹规划方法步骤。
7.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1-4任意一项所述多椎体联合约束的椎弓根螺钉最优轨迹规划方法步骤。
CN202211259373.1A 2022-10-14 2022-10-14 多椎体联合约束的椎弓根螺钉最优轨迹规划方法及装置 Active CN115568942B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211259373.1A CN115568942B (zh) 2022-10-14 2022-10-14 多椎体联合约束的椎弓根螺钉最优轨迹规划方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211259373.1A CN115568942B (zh) 2022-10-14 2022-10-14 多椎体联合约束的椎弓根螺钉最优轨迹规划方法及装置

Publications (2)

Publication Number Publication Date
CN115568942A true CN115568942A (zh) 2023-01-06
CN115568942B CN115568942B (zh) 2024-10-29

Family

ID=84584764

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211259373.1A Active CN115568942B (zh) 2022-10-14 2022-10-14 多椎体联合约束的椎弓根螺钉最优轨迹规划方法及装置

Country Status (1)

Country Link
CN (1) CN115568942B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116740341A (zh) * 2023-06-01 2023-09-12 上海极睿医疗科技有限公司 一种自动生成椎弓根螺钉最优置入路径的方法与系统
CN119405418A (zh) * 2024-09-14 2025-02-11 中欧智薇(上海)机器人有限公司 椎弓根螺钉路径规划方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106659524A (zh) * 2013-06-07 2017-05-10 乔治·弗雷 用于执行外科手术的患者匹配型器械以及方法
US20180116727A1 (en) * 2016-11-02 2018-05-03 Rochester Institute Of Technology Method and fixture for guided pedicle screw placement
CN109069166A (zh) * 2016-04-28 2018-12-21 皇家飞利浦有限公司 确定椎弓根螺钉的最佳放置
US20190231435A1 (en) * 2016-06-23 2019-08-01 Mazor Robotics Ltd. Minimally invasive intervertebral rod insertion
US20190321193A1 (en) * 2018-04-16 2019-10-24 Carlsmed, Inc. Systems and methods for orthopedic implant fixation
WO2021061877A2 (en) * 2019-09-24 2021-04-01 Nuvasive, Inc. Robot assisted rod bending
CN113940752A (zh) * 2021-11-11 2022-01-18 首都医科大学 一种基于多角度投影的椎弓根螺钉最优路径规划方法
EP4011311A1 (en) * 2020-12-14 2022-06-15 Incremed AG Method for planning the nail path of a pedicle bone screw
US20220265352A1 (en) * 2021-02-25 2022-08-25 Mazor Robotics Ltd. Cross-modality planning using feature detection

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106659524A (zh) * 2013-06-07 2017-05-10 乔治·弗雷 用于执行外科手术的患者匹配型器械以及方法
CN109069166A (zh) * 2016-04-28 2018-12-21 皇家飞利浦有限公司 确定椎弓根螺钉的最佳放置
US20190133690A1 (en) * 2016-04-28 2019-05-09 Koninklijke Philips N.V. Determining an optimal placement of a pedicle screw
US20190231435A1 (en) * 2016-06-23 2019-08-01 Mazor Robotics Ltd. Minimally invasive intervertebral rod insertion
US20180116727A1 (en) * 2016-11-02 2018-05-03 Rochester Institute Of Technology Method and fixture for guided pedicle screw placement
US20190321193A1 (en) * 2018-04-16 2019-10-24 Carlsmed, Inc. Systems and methods for orthopedic implant fixation
WO2021061877A2 (en) * 2019-09-24 2021-04-01 Nuvasive, Inc. Robot assisted rod bending
EP4011311A1 (en) * 2020-12-14 2022-06-15 Incremed AG Method for planning the nail path of a pedicle bone screw
US20220265352A1 (en) * 2021-02-25 2022-08-25 Mazor Robotics Ltd. Cross-modality planning using feature detection
CN113940752A (zh) * 2021-11-11 2022-01-18 首都医科大学 一种基于多角度投影的椎弓根螺钉最优路径规划方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
于海龙;项良碧;刘军;陈语;曹艳;: "计算机辅助颈椎椎弓根螺钉置入参数的获取", 中国骨科临床与基础研究杂志, no. 02, 15 December 2009 (2009-12-15) *
刘嘉斌,夏桂锋,左睿: "基于CT值的脊柱椎弓根螺钉轨迹自动规划在骨质疏松患者中的应用价值研究", 中国医学装备, vol. 19, no. 2, 28 February 2022 (2022-02-28) *
邵振轩;吕青波;车灿文;唐笠;王鉴顺;池永龙;王向阳;吴爱悯;: "腰椎一侧椎弓根螺钉联合对侧关节突螺钉3D打印导板的设计和实验研究", 转化医学电子杂志, no. 03, 31 March 2018 (2018-03-31) *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116740341A (zh) * 2023-06-01 2023-09-12 上海极睿医疗科技有限公司 一种自动生成椎弓根螺钉最优置入路径的方法与系统
CN119405418A (zh) * 2024-09-14 2025-02-11 中欧智薇(上海)机器人有限公司 椎弓根螺钉路径规划方法

Also Published As

Publication number Publication date
CN115568942B (zh) 2024-10-29

Similar Documents

Publication Publication Date Title
CN109567942B (zh) 采用人工智能技术的颅颌面外科手术机器人辅助系统
CN109561919B (zh) 微创椎间杆插入
JP7813793B2 (ja) ナビゲーション気管支鏡検査のための動的変形追跡
CA2929319C (en) System and method for generating partial surface from volumetric data for registration to surface topology image data
CN115568942A (zh) 多椎体联合约束的椎弓根螺钉最优轨迹规划方法及装置
JP2022500148A5 (zh)
CN112006777A (zh) 基于表面跟踪的打钉手术机器人系统及控制方法
US20220361959A1 (en) System and Method for Computation of Coordinate System Transformations
CN113781495A (zh) 面向并联骨折手术机器人的复位轨迹自动式规划方法
CN103251457A (zh) 一种用于神经外科立体定向手术的头面部标记点半自动标定方法
CN109360213B (zh) 一种基于脊柱超声冠状面图像的自动化椎体识别方法
US11481969B2 (en) Systems, methods, and computer-readable media for automatic computed tomography to computed tomography registration
CN113850810A (zh) 用于转正股骨的方法及手术系统、存储介质以及电子设备
US20230245589A1 (en) Computer implemented method, a computing device and a system for assisting bending of a reinforcing rod for attachment to a plurality of chirurgical implants
CN115468491B (zh) 一种舱段对接测量方法
Yang et al. Autonomous laparoscope control for minimally invasive surgery with intuition and RCM constraints
CN114404039B (zh) 三维模型的组织漂移校正方法、装置、电子设备及存储介质
CN116363093B (zh) 用于寻找髋臼的旋转中心的方法与装置、手术规划系统
CN115731285B (zh) 机器人末端控制方法、装置、电子设备及存储介质
CN215778615U (zh) 基于自定义标尺的一体化手术机器人术中配准系统
CN119074216B (zh) 截骨机器人辅助作业轨迹生成方法、系统及存储介质
CN114451992A (zh) 一种术后置钉精度评价方法
CN113313766A (zh) 机器人及其位姿计算方法、位姿计算装置
CN115880469B (zh) 一种表面点云数据与三维影像的配准方法
CN119648674A (zh) 一种椎弓根螺钉规划的方法及系统

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