CN115363717A - 一种脑部穿刺路径精确规划方法及3d打印的穿刺定位装置 - Google Patents

一种脑部穿刺路径精确规划方法及3d打印的穿刺定位装置 Download PDF

Info

Publication number
CN115363717A
CN115363717A CN202211077113.2A CN202211077113A CN115363717A CN 115363717 A CN115363717 A CN 115363717A CN 202211077113 A CN202211077113 A CN 202211077113A CN 115363717 A CN115363717 A CN 115363717A
Authority
CN
China
Prior art keywords
puncture
brain
tumor
dimensional matrix
point
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
CN202211077113.2A
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.)
Shantou University Medical College
Original Assignee
Shantou University Medical College
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 Shantou University Medical College filed Critical Shantou University Medical College
Priority to CN202211077113.2A priority Critical patent/CN115363717A/zh
Publication of CN115363717A publication Critical patent/CN115363717A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B17/00Surgical instruments, devices or methods, e.g. tourniquets
    • A61B17/34Trocars; Puncturing needles
    • A61B17/3403Needle locating or guiding means
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B10/00Other methods or instruments for diagnosis, e.g. instruments for taking a cell sample, for biopsy, for vaccination diagnosis; Sex determination; Ovulation-period determination; Throat striking implements
    • A61B10/02Instruments for taking cell samples or for biopsy
    • A61B10/0233Pointed or sharp biopsy instruments
    • 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
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/501Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the head, e.g. neuroimaging or craniography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5217Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5258Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5294Devices using data or image processing specially adapted for radiation diagnosis involving using additional data, e.g. patient information, image labeling, acquisition parameters
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B90/00Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges
    • A61B90/10Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges for stereotaxic surgery, e.g. frame-based stereotaxis
    • A61B90/11Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges for stereotaxic surgery, e.g. frame-based stereotaxis with guides for needles or instruments, e.g. arcuate slides or ball joints
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B90/00Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges
    • A61B90/50Supports for surgical instruments, e.g. articulated arms
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B17/00Surgical instruments, devices or methods, e.g. tourniquets
    • A61B17/34Trocars; Puncturing needles
    • A61B17/3403Needle locating or guiding means
    • A61B2017/3405Needle locating or guiding means using mechanical guide means
    • A61B2017/3407Needle locating or guiding means using mechanical guide means including a base for support on the body
    • 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
    • A61B90/00Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges
    • A61B90/50Supports for surgical instruments, e.g. articulated arms
    • A61B2090/502Headgear, e.g. helmet, spectacles

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Molecular Biology (AREA)
  • Biomedical Technology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Optics & Photonics (AREA)
  • Biophysics (AREA)
  • Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Robotics (AREA)
  • Neurology (AREA)
  • Neurosurgery (AREA)
  • Dentistry (AREA)
  • Physiology (AREA)
  • Pulmonology (AREA)
  • Theoretical Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本发明公开了一种脑部穿刺路径精确规划方法,为病灶定位并自动选取特征点算法、以特征点为三角形参考平面自动生成穿刺路径算法,病灶定位并自动选取特征点算法通过创建三维矩阵、数据处理、图像处理、设置阈值扫描肿瘤,获得肿瘤坐标,并根据肿瘤位置自动选取特征点并提示。以特征点为三角形参考平面自动生成穿刺路径算法通过CT扫描获得特征点距离,结合肿瘤位置、进针点位置,通过数学关系计算得出穿刺角度和穿刺路径。本发明还公开了3D打印的脑部穿刺定位装置。采用本发明,事先形成穿刺路径,可实现穿刺路径可视化和精准化,辅助穿刺设备穿戴简便,便于医生穿刺操作,且可根据医生和患者的需求进行调整,减轻患者医疗负担。

Description

一种脑部穿刺路径精确规划方法及3D打印的穿刺定位装置
技术领域
本发明涉及医疗辅助类器械技术领域,尤其涉及一种3D打印脑部穿刺路径精确规划方法及穿刺定位装置。
背景技术
脑部肿瘤压迫血管、神经、脑组织;脑部出血形成的颅内血肿造成颅内压增高,压迫脑组织进而造成神经损伤。脑部出血、积水、囊肿、肿瘤都严重影响人类的健康和生活,对人类造成巨大的医疗负担。脑部穿刺活检是诊断和治疗脑部肿瘤和颅内血肿的一种方法,具有手术风险低、创伤小、临床效果显著等优点。穿刺活检通过穿刺针取出一部分病变组织进行病理检查,通过病理检查诊断肿瘤。穿刺的精确度是穿刺手术的重要影响因素。
目前脑部穿刺技术分为两种:一是传统的脑部穿刺方法,外科医生通过CT扫描确定病灶和头皮进针点后,根据经验自行把握进针角度及深度。但对医生穿刺手术经验要求高;二是在立体定向、神经导航等设备辅助下进行穿刺。但设备昂贵,且会增加创面。
一种改进的颅内穿刺导向方法(201510507773.3),包括如下步骤:1)穿刺规划导航,1.1)加载头颅扫描图像序列,建立脑部三维模型,同时加载测量角度修正数据;1.2)设置基准点,在三维模型中建立相互垂直的基准参考平面;1.3)设置靶点和进入点,计算出靶点和进入点组成的穿刺路径的长度及其与基准参考平面之间的夹角;通过角度修正数据进行修正,得到最终计划穿刺角度和深度。
一种基于3D打印个性化的颅脑粒子植入导向系统(2017207405506),该发明为一种3D打印个性化颅脑粒子植入导向系统,包括头部固定真空垫、粒子植入计划软件、个性化版骨打孔模板、个性化粒子植入穿刺校板、3D建模模板制作和打印软件、3D打印机。本发明顽骨打孔模板与粒子植入穿刺模板配套使用,使每根粒子植入针体交叉点在颅骨打孔处,穿刺针顺利通过颅骨,模板的粒子植入针柱可以精确引导穿刺针的穿刺方向,避开血管和重要中枢神经,精准穿刺到病灶,减小手术难度,缩短手术时间,降低手术风险。具体实施方式如下:患者佩戴头部固定真空垫进行增强CT扫描,根据扫描结果计算穿刺针道;将CT图像导入3D建模软件并输入术前计划针道坐标信息,重建出个性化植入穿刺模板和与之配套使用的个性化颅骨打孔模板;3D打印出模板;低温等离子消毒,手术中将个性化精确颅骨打孔模板的定位孔对准头皮对应的标记点,准确复位模板,在骨钻引导柱对应的头皮处局部麻醉头皮并切开,骨钻插入骨钴引导柱中精确颅骨打孔,打孔后麻醉硬脑膜,更换粒子植入穿刺模板,定位孔对准头皮标记点准确复位模板,将粒子植入针插入到粒子植入针柱,按照术前计划插入到颅内肿瘤中,到达位置后植入粒子,取出粒子植入针,更换植入针插入到另一个设计针道,重复上述操作直到全部植入。
中国专利公开号为CN105054999A所公开的一种改进的颅内穿刺导向方法,包括如下步骤:1)穿刺规划导航,1.1)加载头颅扫描图像序列,建立脑部三维模型,同时加载测量角度修正数据;1.2)设置基准点,在三维模型中建立相互垂直的基准参考平面;1.3)设置靶点和进入点,计算出靶点和进入点组成的穿刺路径的长度及其与基准参考平面之间的夹角;通过角度修正数据进行修正,得到最终计划穿刺角度和深度。其虽提供一种通过补偿头颅坐标系与地面坐标系间的角度差完成两个坐标系的匹配、实施简单、便于操作的改进的颅内穿刺导向方法及装置,但仍需要陀螺仪辅助操作,步骤较为繁杂。
中国专利公开号为CN207804782U公开了一种基于3D打印个性化的颅脑粒子植入导向系统(2017207405506),该发明为一种3D打印个性化颅脑粒子植入导向系统,包括头部固定真空垫、粒子植入计划软件、个性化版骨打孔模板、个性化粒子植入穿刺校板、3D建模模板制作和打印软件、3D打印机。本发明顽骨打孔模板与粒子植入穿刺模板配套使用,使每根粒子植入针体交叉点在颅骨打孔处,穿刺针顺利通过颅骨,模板的粒子植入针柱可以精确引导穿刺针的穿刺方向,避开血管和重要中枢神经,精准穿刺到病灶,减小手术难度,缩短手术时间,降低手术风险。具体实施方式如下:患者佩戴头部固定真空垫进行增强CT扫描,根据扫描结果计算穿刺针道;将CT图像导入3D建模软件并输入术前计划针道坐标信息,重建出个性化植入穿刺模板和与之配套使用的个性化颅骨打孔模板;3D打印出模板;低温等离子消毒,手术中将个性化精确颅骨打孔模板的定位孔对准头皮对应的标记点,准确复位模板,在骨钻引导柱对应的头皮处局部麻醉头皮并切开,骨钻插入骨钴引导柱中精确颅骨打孔,打孔后麻醉硬脑膜,更换粒子植入穿刺模板,定位孔对准头皮标记点准确复位模板,将粒子植入针插入到粒子植入针柱,按照术前计划插入到颅内肿瘤中,到达位置后植入粒子,取出粒子植入针,更换植入针插入到另一个设计针道,重复上述操作直到全部植入。其不同模板之间的替换依靠头皮标记点及模板上的定位孔进行定位以便复位模板,然而,头皮活动度较大,头皮上标记点容易发生移位,多个模板之间多次复位,难以保证个性化颅骨打孔模板与前一步骤的头部固定真空垫处于同一位置。同时,多次更换模板难以保证各个模板在头皮上固定的稳定度;CT扫描后需要打印三个模板,耗时较长,且需要多次替换复位模板,操作复杂;3D建模软件操作复杂;粒子植入针柱缺乏植入深度信息指示。
中国专利公开号为CN111195148A提供了一种简易脑部穿刺手术定位支架、设计辅助支架及方法,包括第一支架、高度测量器、固定设置于所述第一支架下方的形成等边三角形关系的三个第一圆柱体,所述高度测量器垂直固定设置于所述第一支架上且测量端可伸出底部,所述高度测量器的设置位置对应所述三个圆柱体的重心位置。首先CT扫描确定患者头皮的进针点;将设计辅助支架固定放置于所述患者头部,并使所述高度测量器垂直于所述进针点;旋转所述高度测量器测量所述第一支架与所述患者头皮的高度差;再对所述患者头部的位置进行第二次CT扫描;测量所述患者头部内的目标点与所述三个第一圆柱体之间的距离;计算所述目标点与所述进针点之间的连线在设定的三维空间中的空间角度;以所述第一支架、三个第一圆柱体位置关系重新设计位置关系一致的第二支架与三个第二圆柱体,并于所述第二支架上设计与所述空间角度一致且与所述目标点与所述进针点之间连线重合的进针管。其中限制为等边三角形,在某些特殊体表位置可能存在固定不稳定的问题;影像特征点非体表特征点,医生难以辨别定位;要求两次CT扫描,不仅增加患者经济压力,而且延长了手术时间可能耽误最佳治疗时机;唯一特定角度的进针管难以实现多种穿刺路径的选择。
发明内容
本发明实施例所要解决的技术问题在于,提供一种脑部穿刺路径精确规划方法及3D打印的穿刺定位装置,无需额外仪器辅助测量,简化操作步骤,3D打印方便装置制作,3D打印材料成本低、环保,降低脑部穿刺成本,提高穿刺定位精确度,可用于解决脑部穿刺缺乏相应辅助器械的问题,达到穿刺路径的可视化和精准化。
为了解决上述技术问题,本发明实施例提供了一种脑部穿刺路径精确规划方法,包括以下步骤:
S1:输入CT图像序列,并创建三维矩阵;
S2:扫描所述三维矩阵并定位脑部肿瘤的位置;
S3:根据肿瘤位置选取三个特征点并提示;
S4:根据所述三个特征点形成三角形参考平面并建立空间直角坐标系;
S5:寻找特征点与靶点的关系,通过直线与平面距离建立数学关系,计算穿刺角度、穿刺距离,形成穿刺路径。
其中,所述S1还包括以下步骤:
S11:获取所述三维矩阵的其中一个平面;
S12:使用降噪、二值化排除无效信息的处理步骤对所述平面的对应的二维矩阵进行处理;
S13:将处理后的二维矩阵更新所述三维矩阵;
S14:重复上述步骤直到所述三维矩阵中所有二维平面均处理完毕。
其中,所述S2还包括以下步骤:
为扫描检测脑部肿瘤创建3*3*3检测矩阵,根据肿瘤特征设置阈值,历遍三维矩阵进行扫描,当符合阈值时停止扫描并返回肿瘤坐标。
其中,所述S3还包括以下步骤:
S31:预先获取CT序列中中间的CT图像;
S32:对所述CT图像进行二值化处理,获取脑部图像上下边缘;
S33:根据灰度值特点,进行脑部图像上下边缘检索,获得上边缘凸点所在行数记为top,下边缘凸点所在行数记为bollow;
S34:判断获得的肿瘤坐标位于脑部前半部分或是后半部分。
其中,所述S34还包括以下步骤:
肿瘤行坐标row在上边缘凸点所在行数top与中间行数(x-bollow+top)/2范围内,判断为肿瘤位于脑部前半部分;
肿瘤行坐标row在中间行数(x-bollow+top)/2)与下边缘凸点所在行数bollow范围内,判断为肿瘤位于脑部后半部分;
若在脑部前半部分,选取两侧耳洞和一侧眉骨作为特征点;
若在脑部后半部分,选取两侧耳洞和第七颈椎作为特征点。
其中,所述S5还包括步骤:预先通过DICOM软件三维重建并获取肿瘤与入针点、三个特征点间的距离,计算入穿刺角度的垂直旋转角、水平旋转角,形成穿刺路径。
相应地,本发明实施例还提供了一种3D打印的脑部穿刺定位装置,包括支撑框架、半圆挂式固定部件、挂耳式固定部件、穿刺引导管,所述半圆挂式固定部件设置于所述支撑框架的其中一个面的内侧,所述支撑框架的顶面或相临所述半圆挂式固定部件的侧面设置有参考平面板,所述穿刺引导管以上述所得的穿刺角度设置于所述参考平面板上,所述挂耳式固定部件设置于所述半圆挂式固定部件的两个端部,用于与耳部连接。
其中,所述支撑框架于所述半圆挂式固定部件对向的下边上设置有眉骨吻合滑块。
实施本发明实施例,具有如下有益效果:本发明只需进行一次CT扫描获取所需数据对CT序列进行分析形成穿刺路径,无需其他装置辅助测量;辅助穿刺设备穿戴便捷,不影响医生操作,不会产生额外创面,设备成本低,本发明通过软件实现算法辅助计算,实现穿刺路径可视化和精准化,辅助穿刺设备穿戴简便,便于医生穿刺操作,且可根据医生和患者的需求进行调整,减轻患者医疗负担。
附图说明
图1表示本发明方法实现流程图;
图2表示本发明以特征点为三角形参考平面自动生成穿刺路径算法中的空间直角坐标系图解;
图3是本发明的脑部穿刺定位装置的简化结构示意图;
图4是脑部穿刺定位装置的结构示意图;
图5是脑部穿刺定位装置另一结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步地详细描述。
如图1所示,本发明实施例的一种脑部穿刺路径精确规划方法,可以解决多种问题,比如脑部出血、积水、囊肿、肿瘤;以脑部肿瘤为例,通过以下步骤进行实施。
S1:病灶定位并智能选取特征点算法:通过MATLAB实现输入CT图像序列后创建三维矩阵,扫描定位脑部肿瘤位置,并且根据脑部肿瘤位置自动选取三个特征点并提示。
具体通过下述步骤进行:
S1.1读取CT序列并创建三维矩阵;
S1.2数据处理,将CT值转化为灰度值;
S1.3图像进行降噪、二值化、排除无效信息处理;
S1.4对三维灰度矩阵进行扫描,定位肿瘤;
S1.5根据肿瘤位置自动选取三个特征点并提示。
其中,S1.1通过以下步骤进行:
S1.1.1读取单个DICOM文件:CT1=dicomread('文件名');
S1.1.2为进行三维矩阵的创建,重复上述步骤至完成CT序列读取得到CT2、......、CTn;
S1.1.3创建三维矩阵:三维矩阵=cat(3,CT1,CT2,......,CTn)。
S1.2通过以下步骤进行:
S1.2.1预先读取CT数据元素:
CTdata=dicominfo('文件名');
S1.2.2获取缩放截距:
缩放截距=CTdata.RescaleIntercept;
S1.2.3预先获取缩放斜度:
缩放斜度=CTdata.RescaleSlope;
S1.2.4根据公式将CT值转为灰度值(pixelvalue):
(公式:HU=pixelvalue*RescaleSlope+RescaleIntercept;)
注:(公式:CT值=灰度值*缩放斜坡+缩放截距;)
三维矩阵(灰度值)=(三维矩阵(CT值)-缩放截距)/缩放斜度;
S1.2.5预先将灰度值范围转化成0-1:
最大灰度值=max(max(max(三维矩阵)));
最小灰度值=min(min(min(三维矩阵)));
三维矩阵=(原三维矩阵-最小灰度值)/(最大灰度值-最小灰度值)
S1.3通过以下步骤进行:
S1.3.1预先获取三维矩阵的第一个平面:二维矩阵=三维矩阵(:,:,1);
S1.3.2为排除图像噪点,避免噪点对肿瘤检测造成干扰,通过均值滤波器对数据进行降噪:
滤波掩膜=fspecial('average',3);
二维矩阵=imfilter(原二维矩阵,滤波掩膜,'replicate');
S1.3.3预先将图像二值化,获取脑部轮廓:
二值化二维矩阵=imbinarize(二维矩阵);
S1.3.4预先进行填充图像,获取脑部主体:填充后二维矩阵=imfill(二值化二维矩阵,'holes');
S1.3.5为将脑部主体外无效信息进行排除,避免影响扫描三维矩阵并定位肿瘤,预先排除无效信息:二维矩阵(填充后二维矩阵==0)=1;
S1.3.6为更新数据,二维矩阵处理后需重新输入三维矩阵,将处理后的二维矩阵输入三维矩阵:原三维矩阵(:,:,i)=新二维矩阵;
S1.3.7重复S1.3.1—S1.3.6步骤直至三维矩阵中所有二维平面均处理完毕。
S1.4通过以下步骤进行:
S1.4.1为进行步骤S1.4.2对三维矩阵进行历遍扫描,需要预先获取三维矩阵大小数据进行历遍操作:[行,列,层]=size(三维矩阵);
S1.4.2为扫描检测脑部肿瘤创建3*3*3检测矩阵,根据肿瘤特征设置阈值,历遍三维矩阵进行扫描,当符合阈值时停止扫描并返回肿瘤坐标。
S1.4.3对于S1.4.2步骤中设置阈值时,阈值根据不同肿瘤类型的特征而改变;如血肿,CT值范围一般为60-90,根据S1.2.1-S1.2.5步骤计算得灰度值(0-1范围)为0.515-0.5225,所以检测阈值设置为0.515-0.5225;
(脑部平均CT值一般为30-40,灰度值(0-1范围)为0.5075-0.51)。
S1.5通过以下步骤进行:
S1.5.1为进行步骤S1.5.5对肿瘤位置属于脑部前半部分或是脑部后半部分的判断,需要预先获取CT序列中中间的CT图像,即三维矩阵中第1/2层:
中间层二维矩阵=CT_3Dmatrix(:,:,z/2);(z为总层数);
S1.5.2二值化处理,方便获取脑部图像上下边缘(上边缘位于额骨,下边缘位于枕骨);
S1.5.3脑部轮廓灰度值为1,根据灰度值特点,进行脑部图像上下边缘检索,获得上边缘凸点所在行数记为top,下边缘凸点所在行数记为bollow,方便步骤S1.5.5对肿瘤位置的判断;
S1.5.4为选取三个特征点,判断步骤S1.4.3中获得的肿瘤坐标位于脑部前半部分或是后半部分;
S1.5.5步骤S1.5.4判断方法:通过肿瘤行坐标与二维图像上下边缘凸点所在行和中间行数的关系进行判断;
肿瘤行坐标row在上边缘凸点所在行数top与中间行数(x-bollow+top)/2范围内,即肿瘤位于脑部前半部分;
(公式:row>top&row<(x-bollow+top)/2);
肿瘤行坐标row在中间行数(x-bollow+top)/2)与下边缘凸点所在行数bollow范围内,即肿瘤位于脑部后半部分;
(公式:row<bollow&row>(x-bollow+top)/2);
其中x为二维矩阵行数;
S1.5.6若在脑部前半部分,即返回“肿瘤位于脑部前半部分,选取两侧耳洞和一侧眉骨作为特征点”;
S1.5.7若在脑部后半部分,即返回“肿瘤位于脑部后半部分,选取两侧耳洞和第七颈椎作为特征点”。
S2:根据特征点为三角形参考平面自动生成穿刺路径算法:根据所选取的三个特征点形成三角形参考平面建立空间直角坐标系,寻找特征点与靶点的关系,通过直线与平面距离建立数学关系,计算穿刺角度、穿刺距离,形成穿刺路径。
步骤S2具体通过以下步骤进行实施。
S2.1 根据提示选取特征点,通过CT扫描获取三个点的距离;
S2.2 三个特征点作为三角形参考平面建立空间直角坐标系,通过三角函数、三角形面积建立数学关系,通过直线与平面角度关系建立数学关系计算得出穿刺角度,形成穿刺路径,具体方法步骤为S2.1.1-S2.2.11;
步骤S2.1通过以下过程实施:
S2.1.1 为计算穿刺路径的距离、垂直旋转角和水平旋转角,需要根据步骤S1.5.6-S1.5.7提示,预先通过DICOM软件三维重建并获取肿瘤与入针点、三个特征点间的距离;
S2.1.2 DICOM软件为Radiant DICOM Viewer。
S2.2.1 为求得垂直旋转角(向量EF与三角形平面ABC夹角θ)和水平旋转角(向量E’F’与Y轴的夹角α)(参考图2),明确已知条件并结合图2进行合理设定:
三角形ABC为三个特征点,点A、B、C分别代表三角形ABC参考平面的三个对应的顶点,以A为眉骨,B、C为两侧耳洞为例:设置E点为肿瘤点,E'为肿瘤点E在平面ABC上的投影点。点F代表头皮进针点,点F’代表在三角形ABC平面上的投影点,点E、F在同一条直线上。由影像重建过程中得到AB长度为n1,AC长度为n2,BC长度为n3;EA长度为a,EB长度为b,EC长度为c;FA长度为a1,FB长度为b1,FC长度为c1,EF长度l。
S2.2.2为建立数学关系,需要预先建立空间直角坐标系并对未知坐标点进行坐标设定:以点B为原点,作EE’⊥平面ABC,FF’⊥平面ABC,以
Figure 640859DEST_PATH_IMAGE001
Figure 107236DEST_PATH_IMAGE002
Figure 226370DEST_PATH_IMAGE003
为方向建立空间直角坐标系;
设E点坐标为
Figure 931283DEST_PATH_IMAGE004
,设F点坐标为
Figure 553894DEST_PATH_IMAGE005
已知AB长度为n1,AC长度为n2,BC长度为n3;
EA长度为a,EB长度为b,EC长度为c;
FA长度为a1,FB长度为b1,FC长度为c1;
∴ B(0,0,0) C(0,n3,0) A(-h,yA,0)
S2.2.3 为步骤S2.2.4求得AB、AC间正弦值,需要预先根据余弦定理求AB、AC余弦值:
Figure 683131DEST_PATH_IMAGE006
S2.2.4 为步骤S2.2.5利用面积公式求A点在x轴上坐标h,需要预先根据三角函数求AB、AC求正弦值:
Figure 961666DEST_PATH_IMAGE007
S2.2.5为步骤S2.2.6利用勾股定理求A点在y轴上坐标,需要预先利用面积公式求A点在x轴上坐标h:
Figure 968805DEST_PATH_IMAGE008
Figure 947387DEST_PATH_IMAGE009
Figure 889936DEST_PATH_IMAGE010
S2.2.6为求得点A坐标并建立步骤S2.2.7中方程组求得点E坐标,需要预先利用勾股定理求A点在y轴上坐标:
yA
Figure 452504DEST_PATH_IMAGE011
Figure 433973DEST_PATH_IMAGE012
Figure 31177DEST_PATH_IMAGE011
Figure 941364DEST_PATH_IMAGE013
S2.2.7 通过EF与三角形参考平面的夹角关系求得穿刺路径的垂直旋转角和水平旋转角,需求E点坐标;为求得E点坐标:
根据向量的模长计算公式可得
Figure 492693DEST_PATH_IMAGE014
=
Figure 779318DEST_PATH_IMAGE015
┉┉┉┉┉①
Figure 729563DEST_PATH_IMAGE016
=
Figure 13914DEST_PATH_IMAGE017
┉┉┉┉┉┉┉②
Figure 488757DEST_PATH_IMAGE018
=
Figure 313494DEST_PATH_IMAGE019
┉┉┉┉┉┉③
联立式子①,②,③可得
Figure 557393DEST_PATH_IMAGE020
Figure 12646DEST_PATH_IMAGE021
Figure 538567DEST_PATH_IMAGE022
此时E点有两组解,用EF长度l确定其中一组解,即得E
Figure 104678DEST_PATH_IMAGE004
S2.2.8 为求得F点坐标,根据步骤S2.2.7同理可求得F
Figure 203084DEST_PATH_IMAGE023
,其中:
Figure 891554DEST_PATH_IMAGE024
Figure 839525DEST_PATH_IMAGE025
Figure 6064DEST_PATH_IMAGE026
为求得穿刺路径距离,求得EF模长:
Figure 224556DEST_PATH_IMAGE027
S2.2.9 求垂直旋转角:
Figure 83928DEST_PATH_IMAGE028
与三角形平面ABC夹角θ,平面ABC法向量
Figure 692763DEST_PATH_IMAGE029
,通过直线与平面数学关系求得:
Figure 430037DEST_PATH_IMAGE030
Figure 503036DEST_PATH_IMAGE031
Figure 533309DEST_PATH_IMAGE011
Figure 691758DEST_PATH_IMAGE032
Figure 698634DEST_PATH_IMAGE033
S2.2.10 求水平旋转角:
Figure 891718DEST_PATH_IMAGE034
与Y轴的夹角α,
E在平面ABC上投影E'
Figure 30576DEST_PATH_IMAGE035
,F在平面ABC上投影点F'
Figure 738638DEST_PATH_IMAGE036
,通过直线与平面数学关系求得:
Figure 552135DEST_PATH_IMAGE037
Figure 599725DEST_PATH_IMAGE011
Figure 971801DEST_PATH_IMAGE038
Figure 42525DEST_PATH_IMAGE039
Figure 158249DEST_PATH_IMAGE040
S2.2.11 步骤S2.2.9中获得的垂直旋转角、步骤S2.2.10中获得的水平旋转角为穿刺角度,结合步骤S2.1.1中通过CT扫描获得的肿瘤与入针点的距离,即可形成穿刺路径。
本发明实施例还提供了一种3D打印的脑部穿刺定位装置,包括支撑框架、半圆挂式固定部件、挂耳式固定部件、穿刺引导管,如图3所示。
支撑框架由多面体支架构成,包括了四条等长的位于竖直方向上的直边椭圆立体支架6、7、8、9,以及两条等长的位于底部左右边的直边椭圆立体支架13、14,以及位于顶面四边的螺栓形支架10、11。
直边椭圆立体支架的两端带圆孔,直边椭圆立体支架为两端带圆孔且一端为U形结构。
半圆挂式固定部件16设置于支撑框架的一侧上,支撑框架的顶面或相临所述半圆挂式固定部件的侧面设置有参考平面板1。
当选取两侧耳洞和一侧眉骨作为特征点时,参考平面板1位于支撑框架的顶面。
穿刺引导管以上述方法中所确定的穿刺角度设置于参考平面板上,挂耳式固定部件3设置于所述半圆挂式固定部件的两个端部,用于与耳部连接。
挂耳式固定部件3通过螺栓与半圆挂式固定部件16的端部固定连接。入耳式固定部件4具有螺栓,通过螺帽配合固定于半圆挂式固定部件16的端部,用于置入耳道固定。
作为另一种实施方式,挂耳式固定部件3、入耳式固定部件4的形状与人耳部形状匹配,并可与半圆挂式固定部件16通过3D打印一体成型。
支撑框架于半圆挂式固定部件16对向的下边12上设置有眉骨吻合滑块15,用于与眉骨支撑配合固定本发明的支撑框架。
参考平面板1为四角带圆孔的扁平长方体板,以便通过双头螺栓形支架将参考平面板固定于固定部件,对穿刺部件起到固定作用。
挂耳式固定部件3具有可挂于耳部的弯部。
形成3D打印文件时,需要确定穿刺引导管2与参考平面板1的位置关系,包括长度、角度与交点D的位置,以点坐标作为3D打印参数;因此,在图2建立的坐标系的基础上进行拓展说明;
图3说明:
a.平面ABC、xBy面、BCPQ面为同一平面的不同说法;
b.图2中以A为眉骨,B、C为两侧耳洞为例:以点B为原点,作EE’⊥平面ABC,FF’⊥平面ABC,以
Figure 558881DEST_PATH_IMAGE001
Figure 101858DEST_PATH_IMAGE002
Figure 456616DEST_PATH_IMAGE003
为方向建立空间直角坐标系;设置E点为肿瘤点,E'为肿瘤点E在平面ABC上的投影点。点F代表头皮进针点,点F’代表在三角形ABC平面上的投影点,点E、F在同一条直线上;
c.参考平面B’C’P’Q’平行于xBy面,两平面间距离处处相等,投影为BCPQ;穿刺引导管FG与参考平面B’C’P’Q’相交于倾斜圆柱形孔,设其底面中心点为D;此时EFDG共线,即D、G在EF延长线上,其在xBy面(BCPQ面)投影为E’、F’、D’、G’;
d.BB’为竖直方向上的两端带圆孔的直边椭圆立体支架,长度规定为H;根据c中参考平面B’C’P’Q’平行于xBy面,两平面间距离处处相等,即为BB’长度H;规定点G与参考平面B’C’P’Q’距离为h;
e.EE’、FF’、DD'、GG’相互平行、D’’为DD’上一点,使得FD’’平行于E’G’;
计算:
f.步骤2.2.1-2.2.10已求得A、B、C、E、F、E’、F’各点坐标、穿刺路径的垂直旋转角
Figure 313713DEST_PATH_IMAGE041
、水平旋转角
Figure 70316DEST_PATH_IMAGE042
g.设D坐标为(
Figure 285659DEST_PATH_IMAGE043
Figure 393293DEST_PATH_IMAGE044
Figure 850819DEST_PATH_IMAGE045
),设D’’坐标为(
Figure 461929DEST_PATH_IMAGE046
Figure 549971DEST_PATH_IMAGE047
Figure 377856DEST_PATH_IMAGE048
),设G点坐标为(
Figure 373494DEST_PATH_IMAGE049
Figure 104689DEST_PATH_IMAGE050
Figure 425949DEST_PATH_IMAGE051
);DD’’垂直于xBy面,所以
Figure 914699DEST_PATH_IMAGE043
=
Figure 481072DEST_PATH_IMAGE046
Figure 66774DEST_PATH_IMAGE044
=
Figure 293356DEST_PATH_IMAGE047
;根据e中FD’’平行于E’G’,
Figure 597299DEST_PATH_IMAGE048
=
Figure 722291DEST_PATH_IMAGE052
(步骤2.2.8中求得F在z轴上坐标);即D’’坐标为(
Figure 100183DEST_PATH_IMAGE043
Figure 763245DEST_PATH_IMAGE044
Figure 288904DEST_PATH_IMAGE053
);
h.DD’ = BB’ = H,
Figure 695615DEST_PATH_IMAGE045
= H,可求得DD’’ = DD’ - D’D’’ =
Figure 757374DEST_PATH_IMAGE054
-
Figure 263442DEST_PATH_IMAGE052
④;
i.在三角形FDD’’中,∠DFD’’ =
Figure 541977DEST_PATH_IMAGE041
,由正切定理求得FD’’ = DD’’ / tan
Figure 486799DEST_PATH_IMAGE041
⑤;
j.根据勾股定理,FD’’ =
Figure 901600DEST_PATH_IMAGE055
⑥;
k.联立④⑤⑥式即可求得
Figure 139421DEST_PATH_IMAGE056
,即求得D坐标为(
Figure 577355DEST_PATH_IMAGE043
Figure 325868DEST_PATH_IMAGE044
Figure 595176DEST_PATH_IMAGE054
);
l.根据步骤h、i、j、k同理可求得G点坐标为(
Figure 770942DEST_PATH_IMAGE049
Figure 259955DEST_PATH_IMAGE050
Figure 546579DEST_PATH_IMAGE057
);
m.将所得D、F、G坐标作为参数形成3D打印文件,解决确定穿刺引导管2与参考平面板1的位置关系,包括长度、角度与交点D的位置的问题;其中通过D点坐标确定穿刺引导管2与参考平面板1的位置关系,穿刺引导管2的F端即为穿刺入针点,从而确定穿刺入针点与参考平面板1的位置关系;
n.进行穿刺操作时,医生由G点通过穿刺引导管GF进行穿刺,F点为穿刺入针点;
3D打印时将穿刺引导管2与参考平面板1固定为一体,解决穿刺引导管2与参考平面板1固定问题。
以上所揭露的仅为本发明一种较佳实施例而已,当然不能以此来限定本发明之权利范围,因此依本发明权利要求所作的等同变化,仍属本发明所涵盖的范围。

Claims (8)

1.一种脑部穿刺路径精确规划方法,其特征在于,包括以下步骤:
S1:输入CT图像序列,并创建三维矩阵;
S2:扫描所述三维矩阵并定位脑部肿瘤的位置;
S3:根据肿瘤位置选取三个特征点并提示;
S4:根据所述三个特征点形成三角形参考平面并建立空间直角坐标系;
S5:寻找特征点与靶点的关系,通过直线与平面距离建立数学关系,计算穿刺角度、穿刺距离,形成穿刺路径。
2.根据权利要求1所述的脑部穿刺路径精确规划方法,其特征在于,所述S1还包括以下步骤:
S11:获取所述三维矩阵的其中一个平面;
S12:使用降噪、二值化排除无效信息的处理步骤对所述平面的对应的二维矩阵进行处理;
S13:将处理后的二维矩阵更新所述三维矩阵;
S14:重复上述步骤直到所述三维矩阵中所有二维平面均处理完毕。
3.根据权利要求2所述的脑部穿刺路径精确规划方法,其特征在于,所述S2还包括以下步骤:
为扫描检测脑部肿瘤创建3*3*3检测矩阵,根据肿瘤特征设置阈值,历遍三维矩阵进行扫描,当符合阈值时停止扫描并返回肿瘤坐标。
4.根据权利要求3所述的脑部穿刺路径精确规划方法,其特征在于,所述S3还包括以下步骤:
S31:预先获取CT序列中中间的CT图像;
S32:对所述CT图像进行二值化处理,获取脑部图像上下边缘;
S33:根据灰度值特点,进行脑部图像上下边缘检索,获得上边缘凸点所在行数记为top,下边缘凸点所在行数记为bollow;
S34:判断获得的肿瘤坐标位于脑部前半部分或是后半部分。
5.根据权利要求4所述的脑部穿刺路径精确规划方法,其特征在于,所述S34还包括以下步骤:
肿瘤行坐标row在上边缘凸点所在行数top与中间行数(x-bollow+top)/2范围内,判断为肿瘤位于脑部前半部分;
肿瘤行坐标row在中间行数(x-bollow+top)/2与下边缘凸点所在行数bollow范围内,判断为肿瘤位于脑部后半部分;
若在脑部前半部分,选取两侧耳洞和一侧眉骨作为特征点;
若在脑部后半部分,选取两侧耳洞和第七颈椎作为特征点。
6.根据权利要求5所述的脑部穿刺路径精确规划方法,其特征在于,所述S5还包括步骤:预先通过DICOM软件三维重建并获取肿瘤与入针点、三个特征点间的距离,计算入穿刺角度的垂直旋转角、水平旋转角,形成穿刺路径。
7.一种3D打印的脑部穿刺定位装置,其特征在于,包括支撑框架、半圆挂式固定部件、挂耳式固定部件、穿刺引导管,所述半圆挂式固定部件设置于所述支撑框架的其中一个面的内侧,所述支撑框架的顶面或相临所述半圆挂式固定部件的侧面设置有参考平面板,所述穿刺引导管以权利要求1~6任一项所得的穿刺角度设置于所述参考平面板上,所述挂耳式固定部件设置于所述半圆挂式固定部件的两个端部,用于与耳部连接。
8.根据权利要求7所述的3D打印的脑部穿刺定位装置,其特征在于,所述支撑框架于所述半圆挂式固定部件对向的下边上设置有眉骨吻合滑块。
CN202211077113.2A 2022-09-05 2022-09-05 一种脑部穿刺路径精确规划方法及3d打印的穿刺定位装置 Pending CN115363717A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211077113.2A CN115363717A (zh) 2022-09-05 2022-09-05 一种脑部穿刺路径精确规划方法及3d打印的穿刺定位装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211077113.2A CN115363717A (zh) 2022-09-05 2022-09-05 一种脑部穿刺路径精确规划方法及3d打印的穿刺定位装置

Publications (1)

Publication Number Publication Date
CN115363717A true CN115363717A (zh) 2022-11-22

Family

ID=84069929

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211077113.2A Pending CN115363717A (zh) 2022-09-05 2022-09-05 一种脑部穿刺路径精确规划方法及3d打印的穿刺定位装置

Country Status (1)

Country Link
CN (1) CN115363717A (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2528401Y (zh) * 2001-12-31 2003-01-01 李娜 确定头颅穿刺点的立体定位架
CN102592283A (zh) * 2012-02-21 2012-07-18 华南理工大学 一种脑部肿瘤头皮定位图像处理方法
CN103800051A (zh) * 2014-02-17 2014-05-21 浙江工业大学 一种颅内穿刺辅助导向系统
CN104091347A (zh) * 2014-07-26 2014-10-08 刘宇清 一种基于3d打印技术的颅内肿瘤手术规划模拟方法
WO2019029359A1 (zh) * 2016-10-28 2019-02-14 崔建忠 一种可实现精准定位的脑内血肿穿刺导向仪
CN209018877U (zh) * 2018-04-17 2019-06-25 广东省医疗器械研究所 一种脑支架立体定位仪
CN111407376A (zh) * 2020-04-16 2020-07-14 赵华 一种头颅表面立体定位穿刺颅内脑空间结构的技术
CN113208730A (zh) * 2021-05-18 2021-08-06 赵雨萱 个体化脑空间穿刺导航定位器的制作方法
US20210282743A1 (en) * 2020-02-04 2021-09-16 Tianli Zhao Puncture needle positioning system and method

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2528401Y (zh) * 2001-12-31 2003-01-01 李娜 确定头颅穿刺点的立体定位架
CN102592283A (zh) * 2012-02-21 2012-07-18 华南理工大学 一种脑部肿瘤头皮定位图像处理方法
CN103800051A (zh) * 2014-02-17 2014-05-21 浙江工业大学 一种颅内穿刺辅助导向系统
CN104091347A (zh) * 2014-07-26 2014-10-08 刘宇清 一种基于3d打印技术的颅内肿瘤手术规划模拟方法
WO2019029359A1 (zh) * 2016-10-28 2019-02-14 崔建忠 一种可实现精准定位的脑内血肿穿刺导向仪
CN209018877U (zh) * 2018-04-17 2019-06-25 广东省医疗器械研究所 一种脑支架立体定位仪
US20210282743A1 (en) * 2020-02-04 2021-09-16 Tianli Zhao Puncture needle positioning system and method
CN111407376A (zh) * 2020-04-16 2020-07-14 赵华 一种头颅表面立体定位穿刺颅内脑空间结构的技术
CN113208730A (zh) * 2021-05-18 2021-08-06 赵雨萱 个体化脑空间穿刺导航定位器的制作方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
侯东奥等: "MRI的脑部肿瘤分割及其三维重建", 应用科学学报, vol. 36, no. 5, 30 September 2018 (2018-09-30), pages 808 - 817 *
吴海龙;胥亦龙;罗树芳;毛林花;祝金香;: "颅内血肿微创清除术中穿刺路径立体定位标记线的优化设计研究", 中国实用神经疾病杂志, vol. 2009, no. 3, 15 February 2009 (2009-02-15), pages 7 - 9 *

Similar Documents

Publication Publication Date Title
US20220047244A1 (en) Three dimensional mapping display system for diagnostic ultrasound
Montúfar et al. Automatic 3-dimensional cephalometric landmarking based on active shape models in related projections
CN109419524B (zh) 医学成像系统的控制
CN105963002A (zh) 3d打印微创导向模板及其制造方法
CN202751447U (zh) 一种基于结构光扫描的椎弓根内固定手术导航系统
US20090220136A1 (en) Image Guidance System for Deep Brain Stimulation
CN107126619A (zh) 一种基于3d打印个性化的颅脑粒子植入导向系统
US20230177696A1 (en) Method and system for generating guide information, and computer-readable storage medium
CN105381534A (zh) 粒子植入用导板及其制造方法、装置
JP2008126063A (ja) フルオロスコピィ画像投影への用具及び/又は埋植物の統合を備えた医用ナビゲーション・システム並びに方法
CN102949240A (zh) 一种影像导航肺部介入手术系统
JP2008119472A (ja) 埋植物距離測定のシステム及び方法
CN106308919B (zh) 辅助胸椎椎弓根螺钉置入的导航模板制作方法及导航模板
Hauser et al. Noninvasive tracking of patient's head movements during computer‐assisted intranasal microscopic surgery
CN113425411B (zh) 一种病变定位导航的装置
Tebo et al. An optical 3D digitizer for frameless stereotactic surgery
CN105894508A (zh) 一种医学图像自动定位质量的评估方法
Richey et al. Impact of deformation on a supine-positioned image-guided breast surgery approach
CN108113629A (zh) 硬管内窥镜旋转角度测量方法和装置
Pu et al. A novel ultrasound probe calibration method for multimodal image guidance of needle placement in cervical cancer brachytherapy
Zhang et al. A novel device for preoperative registration and automatic tracking in cranio-maxillofacial image guided surgery
CN115363717A (zh) 一种脑部穿刺路径精确规划方法及3d打印的穿刺定位装置
CN113244516B (zh) 一种基于骨性多点定位的非共面穿刺模板制作方法及模板
CN205658974U (zh) 一种基于3d打印的放射性粒子植入进针定位辅助系统
CN104490415A (zh) 用于腹盆部pet、ct及mr图像融合的立体定向框架

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