CN110123451A - 应用于无标记点光学手术导航系统的病人表面注册方法 - Google Patents
应用于无标记点光学手术导航系统的病人表面注册方法 Download PDFInfo
- Publication number
- CN110123451A CN110123451A CN201910306492.XA CN201910306492A CN110123451A CN 110123451 A CN110123451 A CN 110123451A CN 201910306492 A CN201910306492 A CN 201910306492A CN 110123451 A CN110123451 A CN 110123451A
- Authority
- CN
- China
- Prior art keywords
- point set
- point
- matching
- matrix
- matched
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 22
- 230000003287 optical effect Effects 0.000 title claims abstract description 15
- 239000011159 matrix material Substances 0.000 claims abstract description 44
- 238000013507 mapping Methods 0.000 claims abstract description 11
- 230000009466 transformation Effects 0.000 claims description 23
- 238000013519 translation Methods 0.000 claims description 6
- 238000004321 preservation Methods 0.000 claims description 2
- 208000027418 Wounds and injury Diseases 0.000 description 3
- 210000000988 bone and bone Anatomy 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 239000007943 implant Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 210000003484 anatomy Anatomy 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000002513 implantation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000002560 therapeutic procedure Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/20—Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
Landscapes
- Health & Medical Sciences (AREA)
- Surgery (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Robotics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种应用于无标记点光学手术导航系统的病人表面注册方法,包括步骤:1)获取病人图像空间表面数据:原始点集与简化点集;2)获取病人实际空间表面数据:激光点集与简化点集;3)通过网格化和八方位匹配的方法筛选出最优匹配的粗匹配变换矩阵集合;4)在上一步的基础上求取精匹配的匹配矩阵集合和对应的匹配误差。根据误差筛选出精匹配步骤的最佳变换矩阵。5)根据3)和4)可求取总的最佳变换矩阵,完成空间注册。本发明使用红外激光灯照射病人表面获取实际空间的数据空间注册,采用网格化进行寻优,并将待匹配点集以八种不同的方位与目标点集进行匹配,这两步操作避免了匹配结果陷入局部最优的情况,实现了自动注册的功能。
Description
技术领域
本发明涉及光学手术导航的技术领域,尤其是指一种应用于无标记点光学手术导航系统的病人表面注册方法。
背景技术
手术导航系统在临床上越来越广泛的使用,成为外科手术治疗精准化和微创化的重要发展方向。在手术精度、手术消耗时间、手术创伤、术后疗效等方面与对应的传统外科手术相比有长足的改进。在术中导航时,系统可以通过可视化手术工具和患者特征结构的在实际空间中的位置来帮助医生进行手术,从而提高手术精度和减少对病人的伤害。导航系统的关键步骤是实际空间与图像空间的变换关系,这一步被称为注册。注册方法的优劣直接影响手术导航精度和导航的方便性。
目前手术导航中注册方法多采用标记点注册方法。标记点的种类主要包括骨植入螺钉标记点、解剖标记点和粘贴于皮肤表面的标记点。骨植入螺钉的标记点精度最高,但是需要术前植入器械,会对患者造成额外的创伤和疼痛。解剖标记点利用人体比较显著的解剖特征作为标记点,需要医生在术前提取图像空间的解剖标记点,以及在术中提取患者空间的解剖标记点,存在一定的人为操作误差。粘贴于皮肤表面的标记点精度较高并且操作简单,是最常用的一种标记点。但是在使用过程中存在标记点光线易被遮挡、标记点易掉落以及会由于皮肤的位移形变产生移动等问题,从而影响到导航精度。
发明内容
本发明的目的在于克服现有技术的缺点与不足,提出了一种应用于无标记点光学手术导航系统的病人表面注册方法,该方法有效解决标记点光线易被遮挡,标记点易掉落以及标记点粘贴部位受限等问题,避免了手动提取标记点产生误差的情况,同时,注册过程操作简单,注册时间短并且注册精度可达到临床要求。
为实现上述目的,本发明所提供的技术方案为:应用于无标记点光学手术导航系统的病人表面注册方法,包括以下步骤:
1)导入病人CT图像序列,对CT图像序列提取外表面轮廓,记为点集P,将点集P简化为点集P';
2)用红外激光笔扫描病人表面,光学导航仪实时记录病人表面反射的红外光点的三维坐标,将这些坐标保存为原始激光点集Q,将点集Q简化得到点集Q';
3)分别将点集P'作为目标点集,点集Q'作为原始点集进行分网格初始化粗匹配,得到的候选粗匹配变换矩阵的集合M={Mi|i=1,…,k},其中k为网格的个数;
4)结合步骤3)中粗匹配求得的k种粗匹配变换矩阵结果,将点集P作为匹配的目标点集,将点集Q作为匹配的原始点集进行精匹配,求得精匹配中的最佳匹配矩阵M1s和对应的粗匹配变换矩阵Ms,0<s<k+1;
5)根据步骤4)中的结果得到表面注册的实际空间和图像空间的变换矩阵为Mat=M1s*Ms。
在步骤1)、步骤2)中,点集简化方法是:设原始点集为S,简化后的点集为S',简化边长为l,将点集S用边长为l的正方体进行划分,对于每个正方体,求取其中心点坐标,若在该正方体内部有属于点集S的点,则找到这些点中距离正方体中心点坐标最近的点并保存得到简化后的点集S'。
在步骤2)中,使用的红外激光笔包含两束激光,一束是波长为650nm的可见光,用于指示红外光的位置;另一束是波长为980nm的不可见红外光,将其照射到人体表面,红外光学定位仪获取实际空间人体表面信息。
在步骤3)中,将点集P'作为目标点集和点集Q'作为原始点集进行分网格初始化粗匹配求取候选初始位置对应的变换矩阵的集合,包括以下步骤:
3.1)将点集P'用边长为50mm的正方体网格进行划分,得到n个正方体网格,网格中心点保存为点集C={ci|i=1,…,n},求取点集Q'的中心点cq,对点集C中每一个点ci,分别设ci为点集C的中心点,将点集Q'和点集C进行八方位匹配,以ci为点集C的中心点时点集C和点集Q'的最佳变换矩阵Mi,将矩阵Mi保存到集合M={Mi|i=1,…,k};
3.2)对于步骤3.1)中提到的设置两点集的八方位匹配,设A为原始点集,B为目标点集,其方式是:平移点集A和点集B,使平移后的点集A'和点集B'的中心点都与坐标系原点重合;旋转点集A'和点集B',使旋转后的点集A1”和点集B”的拟平面法向量与坐标系z轴平行,点集B”为目标点集B的匹配初始位置,点集A1”为原始点集A的第一种匹配初始位置;分别将点集A1”绕z轴旋转90°、180°和270°得到点集A2”、A3”、A4”作为第二、三、四种初始位置;将点集A1”绕y轴旋转180°得到点集A5”作为的第五种匹配初始位置;分别将点集A5”绕z轴旋转90°、180°和270°得到点集A6”、A7”、A8”作为第六、七、八种初始位置;分别将点集A1”~A8”与点集B”进行ICP匹配,得到八组变换矩阵和对应的误差,比较误差值,最小的误差对应的变换矩阵则为该八方位匹配步骤的最佳变换矩阵;
在步骤4)中,分别用步骤3.1)中求出来的变换矩阵集合M={Mi|i=1,…,k}中的每一个矩阵Mi对点集P中的坐标进行变换得到点集Pi,将变换后的点集Pi分别与点集Q进行匹配,求得k个匹配误差e={ei|i=1,…,k}和变换矩阵M1={M1i|i=1,…,k},找到最小误差值对应的序号s,0<s<k+1,则精匹配结果求得的变换矩阵为M1s。
本发明与现有技术相比,具有如下优点与有益效果:
1、本发明首次提出了通过将目标点集网格化初始位置匹配和八方位匹配的方式进行匹配寻优,避免在匹配过程中陷入局部最优的情况,从而得到全局最优匹配的结果。
2、本发明通过将患者在图像空间和实际空间的表面信息匹配起来,从而完成空间注册,这一过程中,不需要通过医生手动确定图像空间和实际空间的特征点来完成粗匹配,简化空间注册的操作。
附图说明
图1为本发明方法的流程示意图。
图2为简化的患者图像空间表面点集。
图3为患者实际空间表面激光点集。
图4为将患者图像空间表面点集进行正方体网格化。
图5为八方位匹配中激光点集的八种位置情况。
图6为图像空间点集和实际空间激光点集的最终匹配结果。
具体实施方式
下面结合具体实施例对本发明作进一步说明。
如图1所示,本实施例所提供的应用于无标记点光学手术导航系统的病人表面注册方法,包括以下步骤:
步骤S101,患者术前三维医学图像表面数据的获取。首先是获得CT扫描医学影像序列,然后根据影响序列重建出患者三维解剖结构。用移动立方体算法(MC)从三维解剖结构中提取患者表面数据得到患者原始表面点集P。然后再将P均匀网格简化为点集P'。
步骤S102,患者实际空间表面数据的获取。首先用红外激光笔扫描人体表面,光学定位仪记录患者表面反射的激光点坐标保存为原始激光点集。对记录到的原始激光点集进行去噪得到患者表面激光点集Q。此处使用的红外激光笔包含两束激光,一束是波长为650nm的可见光,用于指示红外光的位置;另一束是波长为980nm的不可见红外光,将其照射到人体表面,红外光学定位仪获取实际空间人体表面信息。然后再将点集Q进行均匀网格简化得到点集Q'。
步骤S103,粗匹配。首先以特定边长的正方体将点集P'进行网格化(参见图4所示),保存这些网格的中心点C={ci|i=1,…,n},作为下一步粗匹配中点集P'的初始位置中心点。根据这n种初始位置中心点,设置点集P'的n种匹配初始位置与点集Q'进行八方位匹配(参见图5所示),从而求得候选的最佳变换矩阵的集合。其中,将点集P'作为目标点集和点集Q'作为原始点集进行分网格初始化粗匹配求取候选初始位置对应的变换矩阵的集合,包括以下步骤:
步骤S1031,将点集P'用边长为50mm的正方体网格进行划分,得到n个正方体网格,网格中心点保存为点集C={ci|i=1,…,n},求取点集Q'的中心点cq,对点集C中每一个点ci,分别设ci为点集C的中心点,将点集Q'和点集C进行八方位匹配,以ci为点集C的中心点时点集C和点集Q'的最佳变换矩阵Mi,将矩阵Mi保存到集合M={Mi|i=1,…,k};
步骤S1032,对于步骤S1031中提到的设置两点集的八方位匹配,设A为原始点集,B为目标点集,其方式是:平移点集A和点集B,使平移后的点集A'和点集B'的中心点都与坐标系原点重合;旋转点集A'和点集B',使旋转后的点集A1”和点集B”的拟平面法向量与坐标系z轴平行,点集B”为目标点集B的匹配初始位置,点集A1”为原始点集A的第一种匹配初始位置;分别将点集A1”绕z轴旋转90°、180°和270°得到点集A2”、A3”、A4”作为第二、三、四种初始位置;将点集A1”绕y轴旋转180°得到点集A5”作为的第五种匹配初始位置;分别将点集A5”绕z轴旋转90°、180°和270°得到点集A6”、A7”、A8”作为第六、七、八种初始位置;分别将点集A1”~A8”与点集B”进行ICP匹配,得到八组变换矩阵和对应的误差,比较误差值,最小的误差对应的变换矩阵则为该八方位匹配步骤的最佳变换矩阵。
步骤S104,精匹配。分别用步骤S1031中求出来的变换矩阵集合M={Mi|i=1,…,k}中的每一个矩阵Mi对点集P中的坐标进行变换得到点集Pi,将变换后的点集Pi分别与点集Q进行匹配,求得k个匹配误差e={ei|i=1,…,k}和变换矩阵M1={M1i|i=1,…,k},找到最小误差值对应的序号s,0<s<k+1,则精匹配结果求得的变换矩阵为M1s。
步骤S105,求取总的变换矩阵。根据S103和S104求得的最佳匹配结果对应的粗匹配步骤和精匹配步骤的变换矩阵,可以求得最终总的最佳变换矩阵为Mat=M1s*Ms。将激光点集P用最佳变换矩阵Mat进行变换,可以得到点集P对应到图像空间中的位置。如图6所示,为最终匹配结果,其中灰度较暗,范围较大的是头部在图像空间中形成的点集,在该点集表面的灰度较亮,范围较小的是激光点集。
以上所述实施例只为本发明之较佳实施例,并非以此限制本发明的实施范围,故凡依本发明之形状、原理所作的变化,均应涵盖在本发明的保护范围内。
Claims (4)
1.应用于无标记点光学手术导航系统的病人表面注册方法,其特征在于,包括以下步骤:
1)导入病人CT图像序列,对CT图像序列提取外表面轮廓,记为点集P,将点集P简化为点集P';
2)用红外激光笔扫描病人表面,光学导航仪实时记录病人表面反射的红外光点的三维坐标,将这些坐标保存为原始激光点集Q,将点集Q简化得到点集Q';
3)分别将点集P'作为目标点集,点集Q'作为原始点集进行分网格初始化粗匹配,得到的候选粗匹配变换矩阵的集合M={Mi|i=1,…,k},其中k为网格的个数;
4)结合步骤3)中粗匹配求得的k种粗匹配变换矩阵结果,将点集P作为匹配的目标点集,将点集Q作为匹配的原始点集进行精匹配,求得精匹配中的最佳匹配矩阵M1s和对应的粗匹配变换矩阵Ms,0<s<k+1;
5)根据步骤4)中的结果得到表面注册的实际空间和图像空间的变换矩阵为Mat=M1s*Ms。
2.根据权利要求1所述的应用于无标记点光学手术导航系统的病人表面注册方法,其特征在于:在步骤1)、步骤2)中,点集简化方法是:设原始点集为S,简化后的点集为S',简化边长为l,将点集S用边长为l的正方体进行划分,对于每个正方体,求取其中心点坐标,若在该正方体内部有属于点集S的点,则找到这些点中距离正方体中心点坐标最近的点并保存得到简化后的点集S'。
3.根据权利要求1所述的应用于无标记点光学手术导航系统的病人表面注册方法,其特征在于:在步骤2)中,使用的红外激光笔包含两束激光,一束是波长为650nm的可见光,用于指示红外光的位置;另一束是波长为980nm的不可见红外光,将其照射到人体表面,红外光学定位仪获取实际空间人体表面信息。
4.根据权利要求1所述的应用于无标记点光学手术导航系统的病人表面注册方法,其特征在于:在步骤3)中,将点集P'作为目标点集和点集Q'作为原始点集进行分网格初始化粗匹配求取候选初始位置对应的变换矩阵的集合,包括以下步骤:
3.1)将点集P'用边长为50mm的正方体网格进行划分,得到n个正方体网格,网格中心点保存为点集C={ci|i=1,…,n},求取点集Q'的中心点cq,对点集C中每一个点ci,分别设ci为点集C的中心点,将点集Q'和点集C进行八方位匹配,以ci为点集C的中心点时点集C和点集Q'的最佳变换矩阵Mi,将矩阵Mi保存到集合M={Mi|i=1,…,k};
3.2)对于步骤3.1)中提到的设置两点集的八方位匹配,设A为原始点集,B为目标点集,其方式是:平移点集A和点集B,使平移后的点集A'和点集B'的中心点都与坐标系原点重合;旋转点集A'和点集B',使旋转后的点集A1”和点集B”的拟平面法向量与坐标系z轴平行,点集B”为目标点集B的匹配初始位置,点集A1”为原始点集A的第一种匹配初始位置;分别将点集A1”绕z轴旋转90°、180°和270°得到点集A2”、A3”、A4”作为第二、三、四种初始位置;将点集A1”绕y轴旋转180°得到点集A5”作为的第五种匹配初始位置;分别将点集A5”绕z轴旋转90°、180°和270°得到点集A6”、A7”、A8”作为第六、七、八种初始位置;分别将点集A1”~A8”与点集B”进行ICP匹配,得到八组变换矩阵和对应的误差,比较误差值,最小的误差对应的变换矩阵则为该八方位匹配步骤的最佳变换矩阵;
在步骤4)中,分别用步骤3.1)中求出来的变换矩阵集合M={Mi|i=1,…,k}中的每一个矩阵Mi对点集P中的坐标进行变换得到点集Pi,将变换后的点集Pi分别与点集Q进行匹配,求得k个匹配误差e={ei|i=1,…,k}和变换矩阵M1={M1i|i=1,…,k},找到最小误差值对应的序号s,0<s<k+1,则精匹配步骤求得的变换矩阵为M1s。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910306492.XA CN110123451B (zh) | 2019-04-17 | 2019-04-17 | 应用于无标记点光学手术导航系统的病人表面注册方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910306492.XA CN110123451B (zh) | 2019-04-17 | 2019-04-17 | 应用于无标记点光学手术导航系统的病人表面注册方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110123451A true CN110123451A (zh) | 2019-08-16 |
CN110123451B CN110123451B (zh) | 2020-07-28 |
Family
ID=67570039
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910306492.XA Expired - Fee Related CN110123451B (zh) | 2019-04-17 | 2019-04-17 | 应用于无标记点光学手术导航系统的病人表面注册方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110123451B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112837391A (zh) * | 2021-03-04 | 2021-05-25 | 北京柏惠维康科技有限公司 | 坐标转换关系获取方法、装置、电子设备及存储介质 |
CN113274130A (zh) * | 2021-05-14 | 2021-08-20 | 上海大学 | 用于光学手术导航系统的无标记手术注册方法 |
CN114081626A (zh) * | 2021-10-21 | 2022-02-25 | 华南理工大学 | 基于光学定位仪的兼容少量标记点丢失的工具识别方法 |
CN117372661A (zh) * | 2023-12-07 | 2024-01-09 | 华科精准(北京)医疗科技有限公司 | 一种手术导航系统、手术机器人系统及注册方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102651145A (zh) * | 2012-04-06 | 2012-08-29 | 哈尔滨工业大学 | 股骨三维模型可视化方法 |
-
2019
- 2019-04-17 CN CN201910306492.XA patent/CN110123451B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102651145A (zh) * | 2012-04-06 | 2012-08-29 | 哈尔滨工业大学 | 股骨三维模型可视化方法 |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112837391A (zh) * | 2021-03-04 | 2021-05-25 | 北京柏惠维康科技有限公司 | 坐标转换关系获取方法、装置、电子设备及存储介质 |
CN113274130A (zh) * | 2021-05-14 | 2021-08-20 | 上海大学 | 用于光学手术导航系统的无标记手术注册方法 |
CN114081626A (zh) * | 2021-10-21 | 2022-02-25 | 华南理工大学 | 基于光学定位仪的兼容少量标记点丢失的工具识别方法 |
CN114081626B (zh) * | 2021-10-21 | 2023-06-20 | 华南理工大学 | 基于光学定位仪的兼容少量标记点丢失的工具识别方法 |
CN117372661A (zh) * | 2023-12-07 | 2024-01-09 | 华科精准(北京)医疗科技有限公司 | 一种手术导航系统、手术机器人系统及注册方法 |
CN117372661B (zh) * | 2023-12-07 | 2024-03-12 | 华科精准(北京)医疗科技有限公司 | 一种手术导航系统、手术机器人系统及注册方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110123451B (zh) | 2020-07-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110123451A (zh) | 应用于无标记点光学手术导航系统的病人表面注册方法 | |
CN110946654B (zh) | 一种基于多模影像融合的骨科手术导航系统 | |
CN103325143B (zh) | 基于模型匹配的标记点自动注册方法 | |
CN107392995B (zh) | 机械轴导航系统中的人体下肢配准系统 | |
Gunz et al. | Semilandmarks in three dimensions | |
CN106340015B (zh) | 一种关键点的定位方法和装置 | |
JP2966089B2 (ja) | 不均質組織内部における局部手術用対話型装置 | |
CN101474075B (zh) | 微创手术导航系统 | |
CN110464459A (zh) | 基于ct-mri融合的介入计划导航系统及其导航方法 | |
CN110738701B (zh) | 一种肿瘤三维定位系统 | |
CN107714082A (zh) | 一种基于光学定位的超声探头标定装置及方法 | |
CN107392922B (zh) | 基于对称性匹配滤波器组和区域生长的三维血管分割方法 | |
CN104867104B (zh) | 基于xct图像非刚度配准的目标鼠解剖结构图谱获取方法 | |
CN101849813A (zh) | 三维心脏超声虚拟内窥镜系统 | |
US20220304751A1 (en) | Optical scale and method for coordinate system registration | |
CN100418489C (zh) | 手术导航中基于基准面膜的多模式医学图像配准系统 | |
CN112907642B (zh) | 一种配准重合的方法、系统、存储介质及设备 | |
US20240115322A1 (en) | Ultrasonic puncture guidance planning system based on multi-modal medical image registration | |
CN103942772A (zh) | 一种多模态多维度的血管融合方法及系统 | |
CN109994188A (zh) | 基于ndi的神经外科导航配准测试方法和系统 | |
CN103908346A (zh) | 一种高精度自动神经导航空间配准方法 | |
CN109146931A (zh) | 一种三维图像处理方法、系统、装置及存储介质 | |
CN106137395B (zh) | 应用于无标记点光学手术导航系统的全自动病人注册方法 | |
CN116883471B (zh) | 面向胸腹部经皮穿刺的线结构光无接触点云配准方法 | |
CN113648060B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200728 |
|
CF01 | Termination of patent right due to non-payment of annual fee |