CN116687569B - 一种编码标识手术导航方法、系统和存储介质 - Google Patents
一种编码标识手术导航方法、系统和存储介质 Download PDFInfo
- Publication number
- CN116687569B CN116687569B CN202310935763.4A CN202310935763A CN116687569B CN 116687569 B CN116687569 B CN 116687569B CN 202310935763 A CN202310935763 A CN 202310935763A CN 116687569 B CN116687569 B CN 116687569B
- Authority
- CN
- China
- Prior art keywords
- coding
- information
- reference plate
- coding mark
- mark
- 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
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000012545 processing Methods 0.000 claims abstract description 48
- 238000004458 analytical method Methods 0.000 claims abstract description 22
- 230000003287 optical effect Effects 0.000 claims description 41
- 239000011159 matrix material Substances 0.000 claims description 37
- 238000007781 pre-processing Methods 0.000 claims description 16
- 230000009466 transformation Effects 0.000 claims description 15
- 238000012216 screening Methods 0.000 claims description 11
- 230000005484 gravity Effects 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000003706 image smoothing Methods 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 4
- 238000013461 design Methods 0.000 description 7
- 230000009286 beneficial effect Effects 0.000 description 5
- 239000003550 marker Substances 0.000 description 5
- 229910002056 binary alloy Inorganic materials 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 239000004053 dental implant Substances 0.000 description 2
- 239000007943 implant Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T9/00—Image coding
-
- 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/20—Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
- A61B2034/2046—Tracking techniques
- A61B2034/2059—Mechanical position encoders
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Surgery (AREA)
- Theoretical Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Robotics (AREA)
- Physics & Mathematics (AREA)
- Multimedia (AREA)
- Biomedical Technology (AREA)
- General Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明提供了一种编码标识手术导航方法、系统和存储介质,所述方法包括:获取编码标识参考板上带有编码信息的板载二维图像;对板载二维图像进行图像处理,识别各个编码标志点;对各个编码标志点进行解码,得到对应的编码信息;通过编码信息,确定各个编码标志点所在参考板上的对应位置;实时扫描板载二维图像,分析得到各个编码标志点的三维坐标;通过多个编码标志点的三维坐标分析获取编码标识参考板的位姿信息,从而对参考板进行实时定位。编码标识参考板上的四个编码标志点通过各自的中心圆区域和信息环带区域携带编码信息,在导航过程中,经过图像处理、解码和实时定位即可确定编码标识参考板的实时位置,从而实现精确导航。
Description
技术领域
本发明涉及图像处理和医疗技术领域,涉及一种用于在实施手术过程中进行导航的编码标识手术导航方法、系统和存储介质。
背景技术
手术导航系统至少包括光学定位仪、手术仪器、人体参考板和仪器参考板。实施手术过程中,人体参考板与患者手术部件保持相互静止,仪器参考板固定设置于手术仪器上,人体参考板和仪器参考板发出或者折射红外光分别被光学定位仪获取,光学定位仪通过标定配准,将人体参考板坐标系与仪器参考板坐标系统一实时导航图像中。如此即可在实施手术过程中为手术仪器的方向、角度等提供参考,从而实现精准导航。
在种植牙手术导航领域,目前市面上口腔种植导航系统主要是通过识别参考板上标记点组成的不同几何形状来实现定位的,其过程主要分为两个步骤:一是识别标记点,二是多个标记点组合的几何形状与模型遍历匹配。这样的方式存在以下问题:标记点均为相同的几何结构,在光学定位仪识别并构建三维坐标时,需通过复杂的算法来确定标记点在光学定位仪两个图像传感器上的对应关系;同时,在设计参考板时,需考虑标记点组合的唯一几何形状,对参考板的设计带来很大限制。业界亟需提供一种新方案,能够更好兼顾导航精准度和设计科学性。
发明内容
本发明要解决的技术问题在于,针对现有技术的上述缺陷,提供一种编码标识手术导航方法、系统和存储介质,在实施手术的过程中基于编码标识参考板进行动态导航,具有导航精准度高和编码标识参考板设计科学合理等优点。
本发明解决技术问题所采用的技术方案如下:
一种编码标识手术导航方法,包括以下步骤:
S1.获取编码标识参考板上带有编码信息的板载二维图像,所述编码标识参考板上设置有四个编码标志点,所述编码标志点包括中心圆区域和信息环带区域;
S2.对板载二维图像进行图像处理,从经过图像处理后的板载二维图像识别出各个编码标志点;
S3.对各个编码标志点进行解码,得到各个编码标志点对应的编码信息;
S4.通过各个编码标志点对应的编码信息,确定各个编码标志点所在参考板上的对应位置;
S5.实时扫描编码标识参考板上各个编码标志点的板载二维图像,分析得到各个编码标志点的三维坐标;
S6.通过多个编码标志点的三维坐标分析获取编码标识参考板的位姿信息,从而对参考板进行实时定位,进而实现手术导航。
与现有技术相比,本技术方案的有益效果是:编码标识参考板上的四个编码标志点通过各自的中心圆区域和信息环带区域携带编码信息,在导航过程中,经过图像处理、解码和实时定位即可确定编码标识参考板的实时位置,从而实现精确导航。
进一步地,所述步骤S2中对板载二维图像进行图像处理具体包括:
对板载二维图像进行预处理;所述预处理包括图像平滑处理、删除噪声处理、灰度化处理、高斯滤波处理;
根据板载二维图像中各个编码标志点的中心圆,对各个编码标志点进行定位;
对定位后的编码标志点进行二值化处理,提取编码标志点的轮廓并进行椭圆拟合;
通过预设的长短轴比例、椭圆尺寸和白色单位环段的面积占信息环带总面积的比例,对拟合得到的椭圆进行筛选,从而识别出各个编码标志点。
采用上述方案的有益效果是:对板载二维图像进行图像处理的目的在于识别出光学定位仪视场范围内所有的编码标志点,通过预处理去除信号干扰、再利用中心圆进行定位、并基于椭圆拟合和椭圆筛选,即可在众多疑似识别点中精确识别出携带编码信息的编码标志点。
进一步地,所述步骤S3中对各个编码标志点进行解码具体包括:
根据中心圆定位得到编码标志点,选取编码标志点上的信息环带;所述信息环带等分为10个首尾相连的单位环段,所述单位环段的颜色为白色或者黑色;
分析计算确定白色单位环段的面积占信息环带总面积的比例,确定白色单位环段的数量;
通过图像识别确定黑色单位环段的角度和位置,从而确定白色单位环段的间隔角度;
从预设位置上沿顺时针方向扫描信息环带,通过白色单位环段的数量和间隔角度,解码得到编码标志点对应的编码信息。
采用上述方案的有益效果是:在编码时将编码信息记载在白色单位环段和黑色单位环段中,利用若干个白色单位环段和若干个黑色单位环段记载编码信息,识别过程中仅需要识别白色单位环段在信息环带区域内的数量、排布方式,即可进行译码,从而区分四个编码标志点。
进一步地,所述步骤S5中分析得到各个编码标志点的三维坐标具体包括:
获取光学定位仪的内参矩阵、畸变矩阵和位置变换矩阵,所述光学定位仪上设置有第一相机和第二相机;
分别通过第一相机和第二相机获取同一编码标志点的坐标,结合位置变换矩阵和两个坐标值分析得到编码标志点的三维坐标。
采用上述方案的有益效果是:基于三角定位测量,通过光学定位仪的第一相机和第二相机,利用内参矩阵、畸变矩阵和位置变换矩阵等对各个编码标志点进行分析,即可分析得到各个编码标志点的三维坐标。
进一步地,所述步骤S6中通过多个编码标志点的三维坐标分析获取编码标识参考板的位姿信息具体包括:
结合四个编码标志点的三维坐标分析出编码标识参考板的重心坐标,从而得到编码标识参考板在光学定位仪的坐标系中的三维坐标;
通过四个编码标志点的三维坐标拟合得到编码标识参考平面;
确定码标识参考平面在光学定位仪的坐标系中的X轴、Y轴和Z轴,进而得到获取编码标识参考板的位姿信息。
采用上述方案的有益效果是:综合四个编码标志点的三维坐标即可分析得到编码标识参考板的重心坐标,拟合出编码标识参考平面并确定该平面在光学定位仪的坐标系中的X轴、Y轴和Z轴后,即可实时获取编码标识参考板的位姿信息。
一种编码标识手术导航系统,包括:
图像获取模板,用于获取编码标识参考板上带有编码信息的板载二维图像,所述编码标识参考板上设置有四个编码标志点,所述编码标志点包括中心圆区域和信息环带区域;
标志点识别模块,用于对板载二维图像进行图像处理,从经过图像处理后的板载二维图像识别出各个编码标志点;
信息解码模块,用于对各个编码标志点进行解码,得到各个编码标志点对应的编码信息;
关联定位模块,用于通过各个编码标志点对应的编码信息,确定各个编码标志点所在参考板上的对应位置;
坐标采集模块,用于实时扫描参考板上各个编码标志点的板载二维图像,分析得到各个编码标志点的三维坐标;
位姿获取模块,用于通过多个编码标志点的三维坐标分析获取编码标识参考板的位姿信息,从而对参考板进行实时定位,进而实现手术导航。
进一步地,所述标志点识别模块具体包括:
图像预处理单元,用于对板载二维图像进行预处理;所述预处理包括图像平滑处理、删除噪声处理、灰度化处理、高斯滤波处理;
标志点定位单元,用于根据板载二维图像中各个编码标志点的中心圆,对各个编码标志点进行定位;
椭圆拟合单元,用于对定位后的编码标志点进行二值化处理,提取编码标志点的轮廓并进行椭圆拟合;
椭圆筛选单元,用于通过预设的长短轴比例、椭圆尺寸和白色单位环段的面积占信息环带总面积的比例,对拟合得到的椭圆进行筛选,从而识别出各个编码标志点。
进一步地,所述信息解码模块具体包括:
环带选取单位,用于根据中心圆定位得到编码标志点,选取编码标志点上的信息环带,所述信息环带等分为10个首尾相连的单位环段,所述单位环段的颜色为白色或者黑色;
数量识别单位,用于分析计算确定白色单位环段的面积占信息环带总面积的比例,确定白色单位环段的数量;
间隔识别单元,用于通过图像识别确定黑色单位环段的角度和位置,从而确定白色单位环段的间隔角度;
解码单元,用于从预设位置上沿顺时针方向扫描信息环带,通过白色单位环段的数量和间隔角度,解码得到编码标志点对应的编码信息。
进一步地,所述坐标采集模块具体包括:
矩阵获取单元,用于获取光学定位仪的内参矩阵、畸变矩阵和位置变换矩阵,所述光学定位仪上设置有第一相机和第二相机;
坐标计算单元,用于分别通过第一相机和第二相机获取同一编码标志点的坐标,结合位置变换矩阵和两个坐标值分析得到编码标志点的三维坐标;
所述位姿获取模块具体包括:
重心分析单元,用于结合四个编码标志点的三维坐标分析出编码标识参考板的重心坐标,从而得到编码标识参考板在光学定位仪的坐标系中的三维坐标;
平面拟合单元,用于通过四个编码标志点的三维坐标拟合得到编码标识参考平面;
位姿获取单元,用于确定码标识参考平面在光学定位仪的坐标系中的X轴、Y轴和Z轴,进而得到获取编码标识参考板的位姿信息。
一种存储介质,所述存储介质存储有计算机程序,所述的计算机程序包括程序指令,当所述程序指令被处理器执行时,处理器执行如上所述的编码标识手术导航方法。
附图说明
图1是本发明编码标识手术导航方法的流程示意图。
图2是本发明编码标识手术导航方法中处理板载二维图像的流程示意图。
图3是本发明编码标识手术导航方法中解码编码标志点的流程示意图。
图4是本发明编码标识手术导航方法中识别编码标志点三维坐标的流程示意图。
图5是本发明编码标识手术导航方法中编码标识参考板位姿信息获取的流程示意图。
图6是本发明编码标识手术导航系统的示意图。
图7是本发明编码标识手术导航方法中编码标识参考板的示意图。
图中,各标号所代表的部件列表如下:
图像获取模板1、标志点识别模块2、信息解码模块3、关联定位模块4、坐标采集模块5、位姿获取模块6、编码标识参考板7。
实施方式
为使本发明的目的、技术方案及优点更加清楚、明确,以下参照附图并举实施例对本发明进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
在本发明的描述中,需要理解的是,术语中“中心”、“上”、“下”、“前”、“后”、“左”、“右”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或组件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性。
在本发明的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“连接”、“相连”应做广义理解,例如,可以是固定连接,也可以是拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以是通过中间媒介间接相连,可以是两个组件内部的连通。当组件被称为“固定于”或“设置于”另一个元件,它可以直接在另一个组件上或者也可以存在居中的组件。当一个组件被认为是“连接”另一个元件,它可以是直接连接到另一个元件或者可能同时存在居中元件。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明的具体含义。
手术导航系统至少包括光学定位仪、手术仪器、人体参考板和仪器参考板。实施手术过程中,人体参考板与患者手术部件保持相互静止,仪器参考板固定设置于手术仪器上,人体参考板和仪器参考板发出红外光分别被光学定位仪获取,光学定位仪通过标定配准,将人体参考板坐标系与仪器参考板坐标系统一实时导航图像中。如此一来,即可在实施手术过程中为手术仪器的方向、角度等提供参考,从而实现精准导航。
在种植牙手术导航领域,目前市面上口腔种植导航系统主要是通过识别参考板上标记点组成的不同几何形状来实现定位的,其过程主要分为两个步骤:一是识别标记点,二是多个标记点组合的几何形状与模型遍历匹配。通常情况下,标记点分为主动式和被动式两种,其中主动式标记点主要为LED发光点或发光球,被动式标记点主要为反光圆盘、反光球、或者黑白格。如下图所示。标记点所组成的几何形状必须是唯一几何,即该几何图形自身不能够产生对称,其位置和方向均唯一。
然而,不管是主动式标记点还是被动式标记点,这样的方式存在以下问题:标记点均为相同的几何结构,在光学定位仪识别并构建三维坐标时,需通过复杂的算法来确定标记点在光学定位仪两个图像传感器上的对应关系;同时,在设计参考板时,需考虑标记点组合的唯一几何形状,对参考板的设计带来很大限制。因此,业界亟需提供一种新方案,能够更好兼顾导航精准度和设计科学性。为了解决上述问题,本发明提供一种编码标识手术导航方法、系统和存储介质。
如图1所示,一种编码标识手术导航方法,所述方法包括以下步骤:
S1.获取编码标识参考板上带有编码信息的板载二维图像,所述编码标识参考板上设置有四个编码标志点,所述编码标志点包括中心圆区域和信息环带区域;
S2.对板载二维图像进行图像处理,从经过图像处理后的板载二维图像识别出各个编码标志点;
S3.对各个编码标志点进行解码,得到对应的编码信息;
S4.通过各个编码标志点对应的编码信息,确定各个编码标志点所在参考板上的对应位置;
S5.实时扫描编码标识参考板上各个编码标志点的板载二维图像,分析得到各个编码标志点的三维坐标;
S6.通过多个编码标志点的三维坐标分析获取编码标识参考板的位姿信息,从而对参考板进行实时定位,进而实现手术导航。
与现有技术中的通过不同形状、不同位置和不同角度来对标志点进行区分为,本发明技术方案的核心是通过编码赋予各个编码标志点不同的数据特征:在编码标识参考板上设置有带有编码信息的板载二维图像,四个编码标志点记录不同特征信息,并在导航过程中通过译码出不同特征信息而识别不同的编码标志点,来对多个编码标志点进行区分,进而获取编码标识参考板的位姿信息,以实现手术过程中实时导航。
具体地,为了在板载二维图像中能区分不同的编码标志点,本发明设计了一种带有编码信息的编码标志点,相机拍摄图片后,可根据设计规则反编译出信息。以图7为例,在编码标识参考板7的板载二维图像中,共有四个编码标志点,每个编码标志点的中间为中心圆区域,在图像识别的过程中通过圆心的位置可精确定位编码标志点的位置。编码标志点周围一圈为信息环带区域,均分为10等份,每一份颜色为黑色或者白色,分别表示二进制的“0”和“1”。根据信息环带区域上白色单位环段的面积比例以及角度,可组合成不同的二进制数。图7中,左上角的编码标志点对应的二进制数为1,右上角的编码标志点对应的二进制数为101,左下角的编码标志点对应的二进制数为111,右下角的编码标志点对应的二进制数为1011。因此,虽然四个编码标志点从结构上看完全一模一样,但通过在信息环带区域内排布不同的白色单位环段,即可区别不同的编码标志点,使得编码标识参考板的结构设置更加科学、合理。总的而言,编码标识参考板上的四个编码标志点通过各自的中心圆区域和信息环带区域携带编码信息,经过图像处理、解码和实时定位等过程,确定编码标识参考板的实时位置,从而实现精确导航。
如图2所示,所述步骤S2中对板载二维图像进行图像处理具体包括以下步骤:
S201.对板载二维图像进行预处理;所述预处理包括图像平滑处理、删除噪声处理、灰度化处理、高斯滤波处理;
S202.根据板载二维图像中各个编码标志点的中心圆,对各个编码标志点进行定位;
S203.对定位后的编码标志点进行二值化处理,提取编码标志点的轮廓并进行椭圆拟合;
S204.通过预设的长短轴比例、椭圆尺寸和白色单位环段的面积占信息环带总面积的比例,对拟合得到的椭圆进行筛选,从而识别出各个编码标志点。
在成像过程中,图像数据会因受到干扰而产生噪声,这些噪声会对后续处理产生负面影响,而这些噪声叠加可用高斯噪声来近似,将输入图像灰度化后,采用高斯滤波去噪。因此,基于上述步骤,通过预处理去除信号干扰、再利用中心圆进行定位、并基于椭圆拟合和椭圆筛选,即可在众多疑似识别点中精确识别出携带编码信息的编码标志点,即可在对板载二维图像进行图像处理后,识别出光学定位仪视场范围内所有的编码标志点。
如图3所示,所述步骤S3中对各个编码标志点进行解码具体包括以下步骤:
S301.根据中心圆定位得到编码标志点,选取编码标志点上的信息环带,所述信息环带等分为10个首尾相连的单位环段,所述单位环段的颜色为白色或者黑色;在此步骤中,由于此前已经确定中心圆位置,同时也可通过图像识别的方式根据获知信息环带的分布方向,由此即可选取编码标志点上的信息环带;
S302.分析计算确定白色单位环段的面积占信息环带总面积的比例,确定白色单位环段的数量;
S303.通过图像识别确定黑色单位环段的角度和位置,从而确定白色单位环段的间隔角度;
S304.从预设位置上沿顺时针方向扫描信息环带,通过白色单位环段的数量和间隔角度,解码得到编码信息。
对照图7,4个编码标志点上的信息环带均等分为10个首尾相连的单位环段,区别不同在于白色单位环段的位置和位置不同,左上角有1个,右上角有2个,左下角有3个,右下角有3个,在图像识别的过程中,可通过计算白色单位环段的面积占信息环带总面积的比例来确定白色单位环段的数量。接着,通过图像识别确定黑色单位环段的角度和位置,从而确定白色单位环段的间隔角度:右上角中,两个白色单位环段之间的间隔角度为72°,由此可确定二者之间存在一个黑色单位环段;左下角中,三个白色单位环段首尾相连;右下角中,一个白色单位环段为一个区域,两个白色单位环段相连,由此可确定二者之间存在一个黑色单位环段。最后,从预设位置上沿顺时针方向扫描信息环带,在本技术方案中,均将编码标志点中线左上侧的区域作为起始位置,如上所述,白色代表二进制中的“1”,黑色代表二进制中的“1”,因此可得,左上角的编码标志点对应的二进制数为1,右上角的编码标志点对应的二进制数为101,左下角的编码标志点对应的二进制数为111,右下角的编码标志点对应的二进制数为1011。至此,通过白色单位环段的数量和间隔角度,即解码得到编码信息。
基于上述技术方案,在编码时将编码信息记载在白色单位环段和黑色单位环段中,利用若干个白色单位环段和若干个黑色单位环段记载编码信息,识别过程中仅需要识别白色单位环段在信息环带区域内的数量、排布方式,即可进行译码,从而区分四个编码标志点。
如图4所示,所述步骤S5中分析得到各个编码标志点的三维坐标具体包括以下步骤:
S501.获取光学定位仪的内参矩阵、畸变矩阵和位置变换矩阵,所述光学定位仪上设置有第一相机和第二相机;使用前,先对光学定位仪进行标定,即可获取内参矩阵、畸变矩阵和位置变换矩阵;
S502.分别通过第一相机和第二相机获取同一编码标志点的坐标,结合位置变换矩阵和两个坐标值分析得到编码标志点的三维坐标。
S503.基于三角定位测量,通过光学定位仪的第一相机和第二相机,利用内参矩阵、畸变矩阵和位置变换矩阵等对各个编码标志点进行分析,即可分析得到各个编码标志点的三维坐标。
如图5所示,所述步骤S6中通过多个编码标志点的三维坐标分析获取编码标识参考板的位姿信息具体包括以下步骤:
S601.结合四个编码标志点的三维坐标分析出编码标识参考板的重心坐标,从而得到编码标识参考板在光学定位仪的坐标系中的三维坐标;
S602.通过四个编码标志点的三维坐标拟合得到编码标识参考平面;
S603.确定码标识参考平面在光学定位仪的坐标系中的X轴、Y轴和Z轴,进而得到获取编码标识参考板的位姿信息。
基于上述技术方案,综合四个编码标志点的三维坐标即可分析得到编码标识参考板的重心坐标,拟合出编码标识参考平面并确定该平面在光学定位仪的坐标系中的X轴、Y轴和Z轴后,即可实时获取编码标识参考板的位姿信息。在实施手术的过程中,基于编码标识参考板的位姿信息,即可确定与之相对位置固定的物体的实时位姿信息,如人体或者手术仪器,从而进行实时导航。
如图6所示,对应地,一种编码标识手术导航系统,包括图像获取模板、标志点识别模块、信息解码模块、关联定位模块、坐标采集模块和位姿获取模块;所述图像获取模板用于获取编码标识参考板上带有编码信息的板载二维图像,所述编码标识参考板上设置有四个编码标志点,所述编码标志点包括中心圆区域和信息环带区域;所述标志点识别模块用于对板载二维图像进行图像处理,从经过图像处理后的板载二维图像识别出各个编码标志点;所述信息解码模块用于对各个编码标志点进行解码,得到对应的编码信息;所述关联定位模块用于通过各个编码标志点对应的编码信息,确定各个编码标志点所在参考板上的对应位置;所述坐标采集模块用于实时扫描参考板上各个编码标志点的板载二维图像,分析得到各个编码标志点的三维坐标;所述位姿获取模块用于通过多个编码标志点的三维坐标分析获取编码标识参考板的位姿信息,从而对参考板进行实时定位,进而实现手术导航。
优选地,所述标志点识别模块具体包括图像预处理单元、标志点定位单元、椭圆拟合单元和椭圆筛选单元;图像预处理单元用于对板载二维图像进行预处理,所述预处理包括图像平滑处理、删除噪声处理、灰度化处理、高斯滤波处理;标志点定位单元用于根据板载二维图像中各个编码标志点的中心圆,对各个编码标志点进行定位;椭圆拟合单元用于对定位后的编码标志点进行二值化处理,提取编码标志点的轮廓并进行椭圆拟合;椭圆筛选单元用于通过预设的长短轴比例、椭圆尺寸和白色单位环段的面积占信息环带总面积的比例,对拟合得到的椭圆进行筛选,从而识别出各个编码标志点。
优选地,所述信息解码模块具体包括环带选取单位、数量识别单位、间隔识别单元和解码单元:所述环带选取单位用于根据中心圆定位得到编码标志点,选取编码标志点上的信息环带,所述信息环带等分为10个首尾相连的单位环段,所述单位环段的颜色为白色或者黑色;所述数量识别单位用于分析计算确定白色单位环段的面积占信息环带总面积的比例,确定白色单位环段的数量;所述间隔识别单元用于通过图像识别确定黑色单位环段的角度和位置,从而确定白色单位环段的间隔角度;所述解码单元用于从预设位置上沿顺时针方向扫描信息环带,通过白色单位环段的数量和间隔角度,解码得到编码信息。
优选地,所述坐标采集模块具体包括矩阵获取单元和坐标计算单元;所述矩阵获取单元用于获取光学定位仪的内参矩阵、畸变矩阵和位置变换矩阵,所述光学定位仪上设置有第一相机和第二相机;所述坐标计算单元用于分别通过第一相机和第二相机获取同一编码标志点的坐标,结合位置变换矩阵和两个坐标值分析得到编码标志点的三维坐标。
优选地,所述位姿获取模块具体包括重心分析单元、平面拟合单元和位姿获取单元;所述重心分析单元用于结合四个编码标志点的三维坐标分析出编码标识参考板的重心坐标,从而得到编码标识参考板在光学定位仪的坐标系中的三维坐标;所述平面拟合单元用于通过四个编码标志点的三维坐标拟合得到编码标识参考平面;所述位姿获取单元用于确定码标识参考平面在光学定位仪的坐标系中的X轴、Y轴和Z轴,进而得到获取编码标识参考板的位姿信息。
对应地,一种存储介质,所述存储介质存储有计算机程序,所述的计算机程序包括程序指令,当所述程序指令被处理器执行时,处理器执行如上所述的编码标识手术导航方法。
应当理解的是,本发明的应用不限于上述的举例,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (3)
1.一种编码标识手术导航系统,其特征在于,包括:
图像获取模板,用于获取编码标识参考板上带有编码信息的板载二维图像,所述编码标识参考板上设置有四个编码标志点,所述编码标志点包括中心圆区域和信息环带区域;
标志点识别模块,用于对板载二维图像进行图像处理,从经过图像处理后的板载二维图像识别出各个编码标志点;
信息解码模块,用于对各个编码标志点进行解码,得到各个编码标志点对应的编码信息;
关联定位模块,用于通过各个编码标志点对应的编码信息,确定各个编码标志点所在参考板上的对应位置;
坐标采集模块,用于实时扫描参考板上各个编码标志点的板载二维图像,分析得到各个编码标志点的三维坐标;
位姿获取模块,用于通过多个编码标志点的三维坐标分析获取编码标识参考板的位姿信息,从而对参考板进行实时定位,进而实现手术导航;
所述标志点识别模块具体包括:
图像预处理单元,用于对板载二维图像进行预处理;所述预处理包括图像平滑处理、删除噪声处理、灰度化处理、高斯滤波处理;
标志点定位单元,用于根据板载二维图像中各个编码标志点的中心圆,对各个编码标志点进行定位;
椭圆拟合单元,用于对定位后的编码标志点进行二值化处理,提取编码标志点的轮廓并进行椭圆拟合;
椭圆筛选单元,用于通过预设的长短轴比例、椭圆尺寸和白色单位环段的面积占信息环带总面积的比例,对拟合得到的椭圆进行筛选,从而识别出各个编码标志点;
所述信息解码模块具体包括:
环带选取单位,用于根据中心圆定位得到编码标志点,选取编码标志点上的信息环带,所述信息环带等分为10个首尾相连的单位环段,所述单位环段的颜色为白色或者黑色;
数量识别单位,用于分析计算确定白色单位环段的面积占信息环带总面积的比例,确定白色单位环段的数量;
间隔识别单元,用于通过图像识别确定黑色单位环段的角度和位置,从而确定白色单位环段的间隔角度;
解码单元,用于从预设位置上沿顺时针方向扫描信息环带,通过白色单位环段的数量和间隔角度,解码得到编码标志点对应的编码信息。
2.根据权利要求1所述的一种编码标识手术导航系统,其特征在于,所述坐标采集模块具体包括:
矩阵获取单元,用于获取光学定位仪的内参矩阵、畸变矩阵和位置变换矩阵,所述光学定位仪上设置有第一相机和第二相机;
坐标计算单元,用于分别通过第一相机和第二相机获取同一编码标志点的坐标,结合位置变换矩阵和两个坐标值分析得到编码标志点的三维坐标;
所述位姿获取模块具体包括:
重心分析单元,用于结合四个编码标志点的三维坐标分析出编码标识参考板的重心坐标,从而得到编码标识参考板在光学定位仪的坐标系中的三维坐标;
平面拟合单元,用于通过四个编码标志点的三维坐标拟合得到编码标识参考平面;
位姿获取单元,用于确定编码标识参考平面在光学定位仪的坐标系中的X轴、Y轴和Z轴,进而得到获取编码标识参考板的位姿信息。
3.一种存储介质,其特征在于,所述存储介质存储有计算机程序,所述的计算机程序包括程序指令,当所述程序指令被处理器执行时,处理器执行以下步骤:
S1.获取编码标识参考板上带有编码信息的板载二维图像,所述编码标识参考板上设置有四个编码标志点,所述编码标志点包括中心圆区域和信息环带区域;
S2.对板载二维图像进行图像处理,从经过图像处理后的板载二维图像识别出各个编码标志点;
S3.对各个编码标志点进行解码,得到各个编码标志点对应的编码信息;
S4.通过各个编码标志点对应的编码信息,确定各个编码标志点所在参考板上的对应位置;
S5.实时扫描编码标识参考板上各个编码标志点的板载二维图像,分析得到各个编码标志点的三维坐标;
S6.通过多个编码标志点的三维坐标分析获取编码标识参考板的位姿信息,从而对参考板进行实时定位,进而实现手术导航;
所述步骤S5中分析得到各个编码标志点的三维坐标具体包括:
获取光学定位仪的内参矩阵、畸变矩阵和位置变换矩阵,所述光学定位仪上设置有第一相机和第二相机;
分别通过第一相机和第二相机获取同一编码标志点的坐标,结合位置变换矩阵和两个坐标值分析得到编码标志点的三维坐标;
所述步骤S6中通过多个编码标志点的三维坐标分析获取编码标识参考板的位姿信息具体包括:
结合四个编码标志点的三维坐标分析出编码标识参考板的重心坐标,从而得到编码标识参考板在光学定位仪的坐标系中的三维坐标;
通过四个编码标志点的三维坐标拟合得到编码标识参考平面;
确定编码标识参考平面在光学定位仪的坐标系中的X轴、Y轴和Z轴,进而得到获取编码标识参考板的位姿信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310935763.4A CN116687569B (zh) | 2023-07-28 | 2023-07-28 | 一种编码标识手术导航方法、系统和存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310935763.4A CN116687569B (zh) | 2023-07-28 | 2023-07-28 | 一种编码标识手术导航方法、系统和存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116687569A CN116687569A (zh) | 2023-09-05 |
CN116687569B true CN116687569B (zh) | 2023-10-03 |
Family
ID=87843536
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310935763.4A Active CN116687569B (zh) | 2023-07-28 | 2023-07-28 | 一种编码标识手术导航方法、系统和存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116687569B (zh) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016026349A1 (zh) * | 2014-08-20 | 2016-02-25 | 深圳大学 | 高鲁棒性的标志点解码方法及系统 |
KR101802918B1 (ko) * | 2016-07-08 | 2017-12-28 | 주식회사 메디트 | 3차원 영상 획득 장치 및 그 동작방법 |
DE102017105158A1 (de) * | 2017-03-10 | 2018-09-13 | Carl Zeiss Industrielle Messtechnik Gmbh | Verfahren und system zum bestimmen einer 3d-position eines objekts im raum |
CN109215016A (zh) * | 2018-08-03 | 2019-01-15 | 湖南科技大学 | 一种编码标志的识别定位方法 |
CN112184826A (zh) * | 2019-07-05 | 2021-01-05 | 杭州海康机器人技术有限公司 | 一种标定板和标定方法 |
CN113712665A (zh) * | 2021-11-01 | 2021-11-30 | 北京柏惠维康科技有限公司 | 基于定位标志物的定位方法、装置及计算机存储介质 |
CN115049744A (zh) * | 2022-07-11 | 2022-09-13 | 深圳市易尚展示股份有限公司 | 机器人手眼坐标转换方法、装置、计算机设备和存储介质 |
CN115345870A (zh) * | 2022-08-29 | 2022-11-15 | 贵州宏图伟业科技发展有限公司 | 一种基于单目相机和带编码自发光靶标实现大场景精密形变监测的方法 |
CN115880373A (zh) * | 2022-12-28 | 2023-03-31 | 常熟理工学院 | 基于新型编码特征的立体视觉系统的标定板及标定方法 |
CN116030240A (zh) * | 2023-01-10 | 2023-04-28 | 深圳市格灵精睿视觉有限公司 | 标志点识别方法、标志点识别装置、电子设备及存储介质 |
CN116077184A (zh) * | 2023-02-07 | 2023-05-09 | 上海诠视传感技术有限公司 | 一种ar手术导航系统的定位装置及方法 |
-
2023
- 2023-07-28 CN CN202310935763.4A patent/CN116687569B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016026349A1 (zh) * | 2014-08-20 | 2016-02-25 | 深圳大学 | 高鲁棒性的标志点解码方法及系统 |
KR101802918B1 (ko) * | 2016-07-08 | 2017-12-28 | 주식회사 메디트 | 3차원 영상 획득 장치 및 그 동작방법 |
DE102017105158A1 (de) * | 2017-03-10 | 2018-09-13 | Carl Zeiss Industrielle Messtechnik Gmbh | Verfahren und system zum bestimmen einer 3d-position eines objekts im raum |
CN109215016A (zh) * | 2018-08-03 | 2019-01-15 | 湖南科技大学 | 一种编码标志的识别定位方法 |
CN112184826A (zh) * | 2019-07-05 | 2021-01-05 | 杭州海康机器人技术有限公司 | 一种标定板和标定方法 |
CN113712665A (zh) * | 2021-11-01 | 2021-11-30 | 北京柏惠维康科技有限公司 | 基于定位标志物的定位方法、装置及计算机存储介质 |
CN115049744A (zh) * | 2022-07-11 | 2022-09-13 | 深圳市易尚展示股份有限公司 | 机器人手眼坐标转换方法、装置、计算机设备和存储介质 |
CN115345870A (zh) * | 2022-08-29 | 2022-11-15 | 贵州宏图伟业科技发展有限公司 | 一种基于单目相机和带编码自发光靶标实现大场景精密形变监测的方法 |
CN115880373A (zh) * | 2022-12-28 | 2023-03-31 | 常熟理工学院 | 基于新型编码特征的立体视觉系统的标定板及标定方法 |
CN116030240A (zh) * | 2023-01-10 | 2023-04-28 | 深圳市格灵精睿视觉有限公司 | 标志点识别方法、标志点识别装置、电子设备及存储介质 |
CN116077184A (zh) * | 2023-02-07 | 2023-05-09 | 上海诠视传感技术有限公司 | 一种ar手术导航系统的定位装置及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN116687569A (zh) | 2023-09-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7927014B2 (en) | X-ray diagnostic imaging system with a plurality of coded markers | |
US6724930B1 (en) | Three-dimensional position and orientation sensing system | |
US8180426B2 (en) | Method for transforming head surface coordinates to brain surface coordinates and transcranial brain function measuring method using the transformation data | |
CN110689579A (zh) | 基于合作目标的快速单目视觉位姿测量方法及测量系统 | |
CN110763204B (zh) | 一种平面编码靶标及其位姿测量方法 | |
CN115880373B (zh) | 基于新型编码特征的立体视觉系统的标定板及标定方法 | |
CN107678153B (zh) | 光纤束图像处理方法和装置 | |
EP4067817A1 (en) | System and method for spatial positioning of magnetometers | |
JP5354426B2 (ja) | 測定装置、測定方法、およびプログラム | |
CN109410282B (zh) | 一种输液杆的检测跟踪方法和装置 | |
CN116309829B (zh) | 一种基于多目视觉的长方体扫描体群解码和位姿测量方法 | |
Robson et al. | Optimised multi-camera systems for dimensional control in factory environments | |
CN113963067B (zh) | 一种采用小靶标对大视场视觉传感器进行标定的标定方法 | |
Schaller et al. | Time-of-flight sensor for patient positioning | |
CN116687569B (zh) | 一种编码标识手术导航方法、系统和存储介质 | |
CN107622491B (zh) | 光纤束图像分析方法和装置 | |
JP3423676B2 (ja) | 仮想対象物操作装置および仮想対象物操作方法 | |
JP3301031B2 (ja) | 物体の自動認識方法および自動認識装置 | |
CN116503486A (zh) | 一种标定系统、标定方法、标定装置及可读存储介质 | |
CN113456228B (zh) | 辅助手术的组织辨识方法及系统 | |
CN109872353A (zh) | 基于改进迭代最近点算法的白光数据与ct数据配准方法 | |
CN116125489A (zh) | 一种室内物体三维检测方法、计算机设备及存储介质 | |
WO2015136499A1 (en) | Three-dimensional target for a vehicle wheel alignment system, and related method | |
CN115546783A (zh) | 一种三维影像下的示踪球定位方法 | |
CN116672082B (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 |