CN108961413A - 一种用于心脏标测系统中的心内膜逐点插入动态表面重建方法 - Google Patents

一种用于心脏标测系统中的心内膜逐点插入动态表面重建方法 Download PDF

Info

Publication number
CN108961413A
CN108961413A CN201810778217.3A CN201810778217A CN108961413A CN 108961413 A CN108961413 A CN 108961413A CN 201810778217 A CN201810778217 A CN 201810778217A CN 108961413 A CN108961413 A CN 108961413A
Authority
CN
China
Prior art keywords
point
replaced
current
tri patch
delaunay
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
CN201810778217.3A
Other languages
English (en)
Other versions
CN108961413B (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.)
Shenzhen Graduate School Tsinghua University
Original Assignee
Shenzhen Graduate School Tsinghua University
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 Shenzhen Graduate School Tsinghua University filed Critical Shenzhen Graduate School Tsinghua University
Priority to CN201810778217.3A priority Critical patent/CN108961413B/zh
Publication of CN108961413A publication Critical patent/CN108961413A/zh
Application granted granted Critical
Publication of CN108961413B publication Critical patent/CN108961413B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/30Polynomial surface description
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • A61B5/004Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part
    • A61B5/0044Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part for the heart
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/20Image enhancement or restoration using local operators
    • G06T5/30Erosion or dilatation, e.g. thinning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/90Dynamic range modification of images or parts thereof
    • G06T5/94Dynamic range modification of images or parts thereof based on local image properties, e.g. for local contrast enhancement
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10028Range image; Depth image; 3D point clouds
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30048Heart; Cardiac

Landscapes

  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Mathematical Physics (AREA)
  • Biomedical Technology (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Cardiology (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Algebra (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Computer Graphics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Processing Or Creating Images (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

Abstract

本发明公开了一种用于心脏标测系统中的心内膜逐点插入动态表面重建方法,包括:A1:导管采集四个仿射独立的点,并通过上述四个仿射独立的点初始化为一个四面体网格;A2:插入新采样点,并且根据新插入的采样点计算当前重建的表面网格中需要被替换的表面三角面片集合;A3:更新Delaunay三角剖分结构,同时也更新当前的表面网格结构;A4:利用局部膨胀和雕刻的方法优化当前重建的表面网格;A5:插入已经在Delaunay剖分结构但不在表面网格中的孤立点,提高表面网格的质量。上述表面重建方法更贴合手术中导管的采点步骤,更友好地服务医生掌握当前采样情况,显著提高心脏标测系统中心内膜表面的重建效率,更好地为医生提高病灶定位和病灶分析的服务。

Description

一种用于心脏标测系统中的心内膜逐点插入动态表面重建 方法
技术领域
本发明涉及计算机图形学与计算几何领域,特别涉及一种用于心脏标测系统中的心内膜逐点插入动态表面重建方法。
背景技术
近些年来,导管消融技术可防止房颤复发在临床上已经被证明是一个相对高效和安全的方法,并且在一些情况下可以治愈房颤。因此,研发高效廉价的导管射频消融仪器具有重要意义。
对于导管消融手术而言,医生首先需要把导管从肱骨动脉处插入病人体内,并移动至主动脉瓣进入心室,接着医生操作导管开始采集心内膜表面上的点云,并在显示终端上显示通过点云重建的心内膜表面。同时根据重建的三维表面,医生能够确定病灶点的位置,并进行相应的消融手术。因此,通过采集的散乱点云可视化一个三维心脏表面是导管消融中至关重要的环节,通过散乱点云重构一个合理的心内膜表面对消融手术靶点定位起着举足轻重的作用。
对于当前流行的心脏标测系统而言,在心内膜上采样的过程是逐点进行的,然而当前系统中配套的表面重建算法多为一次重建,即一次性输入点云数据,然后输出心内膜表面网格。但是这种重建策略不符合导管采样的过程,重建效率低下,并且不方便医生直观地感觉心脏重建的局部变化范围,从而影响到病灶点的定位。
发明内容
为了能体现心脏标测系统的逐点采样过程,重建表面策略也应该被设计为渐进过程,即重建的表面可以随着采样点的新增进行局部更新。为解决上述技术问题,本发明提出一种动态心内膜表面重建方法,可以有效地提高手术过程中的表面重建效率。
为此,本发明采用了如下技术方案,其包括如下步骤:
一种用于心脏标测系统中的心内膜逐点插入动态表面重建方法,包括如下步骤:
A1:导管采集四个仿射独立的点,并通过上述四个仿射独立的点初始化为一个四面体网格;
A2:插入新采样点,并且根据新插入的采样点计算当前重建的表面网格中需要被替换的表面三角面片集合;
A3:更新Delaunay三角剖分结构,同时也更新当前的表面网格结构,即把新采样点插入到当前的表面网格中;
A4:利用局部膨胀和雕刻的方法优化当前重建的表面网格;
A5:插入已经在Delaunay剖分结构但不在表面网格中的孤立点,以再次提高表面网格的质量。
进一步,步骤A2包括:当新插入一个采样点到当前重建的心内膜表面网格中,通过一个种子三角面片,以广度优先的方式进行增长,从而求得需要被替换的表面三角面片集合。
进一步,种子三角面片为选择距离新待插入点最近的表面三角面片,其计算步骤如下:
A21:通过点定位计算当前插入点的四面体单元;
A22:以当前四面体单元为种子,以区域增长的方式求得该种子领域内距离新的待插入点最近的表面三角面片作为种子三角面片。
进一步,计算得到种子表面三角面片后,通过广度优先的方式计算当前表面网格中被替换的表面面片集合,其计算步骤如下:
以区域增长的方式检测种子面片相邻的表面网格上的面片,检测这些面片是否满足条件:在待插入点的Delaunay影响域中,或者不在Delaunay影响域中但不满足Gabriel准则;对于满足条件的表面面片一律被标记为需要被替换的表面面片;如果被替换的表面面片集合不是一个正则的二维单纯形,则不插入当前采样点并重新采样;如果被替换的表面面片集合是一个正则的二维单纯形,则计算其边界。
进一步,步骤A3包括:如果当前网格中被替换的表面面片集合是正则的二维单纯形,在计算被替换面片集合的边界后,以深度优先的方式将当前待插入点用增量式算法插入到Delaunay剖分结构中,在计算剖分的过程中,如果发现有新生成的四面体有一条边落在被替换表面集合的边界中,则生成一个表面面片。
进一步,步骤A4括:在满足不改变二维流形的情况下运用局部膨胀和雕刻,用其滤除表面网格中的新生成的非Gabriel三角面片,对于每个新生成的四面体,局部膨胀是把非Gabriel表面面片关联的Delaunay四面体移到当前重建表面网格的内部,而局部雕刻是把非Gabriel表面面片关联的Delaunay四面体移到当前表面网格的外部。
进一步,步骤A5包括:插入孤立点,就是在当前表面网格中,通过距离孤立点最近的表面面片的领域,计算不满足Gabriel准则的表面面片的集合,并且这个集合应该为正则二维单纯形,然后,将孤立点与这个集合的边界生成表面面片。
相对于现有技术,本发明的有益效果是:重建表面策略不仅更贴合手术过程中导管的采点步骤,而且以局部更新的方式重建心内膜表面能够更友好地服务医生掌握当前采样情况,以便医生可以快速地定位病灶点。
附图说明
图1是本发明的心内膜逐点插入动态表面重建方法流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
结合图1,本发明所提出的心内膜逐点插入动态表面重建方法,包括以下步骤:
A1:选取仿射独立的四个采样点初始化为四面体网格;
心脏标测在术中初始阶段,医生首先通过导管采集四个仿射独立的点,然后通过这四个仿射独立的点构建一个含有无穷远点的Delaunay剖分结构,并将这四个点初始化为一个四面体网格。
A2:插入新采样点,并且根据新插入的采样点计算当前重建的表面网格中需要被替换的表面三角面片集合;
当新插入一个采样点到当前重建的心内膜表面网格中,首先通过一个种子三角面片,以广度优先的方式进行增长,从而求得需要被替换的表面三角面片集合。
选择距离新待插入点最近的表面三角面片为种子三角面片,其计算步骤如下:
A21:通过点定位计算当前插入点的四面体单元;
A22:以当前四面体单元为种子,以区域增长的方式求得该种子领域内距离新的待插入点最近的表面三角面片作为种子三角面片。
计算得到种子表面三角面片后,通过广度优先的方式计算当前表面网格中被替换的表面面片集合,其计算步骤如下:
以区域增长的方式检测种子面片相邻的表面网格上的面片,检测这些面片是否满足条件:在待插入点的Delaunay影响域中,或者不在Delaunay影响域中但不满足Gabriel准则。对于满足条件的表面面片一律被标记为需要被替换的表面面片。如果被替换的表面面片集合不是一个正则的二维单纯形,则不插入当前采样点并重新采样。如果被替换的表面面片集合是一个正则的二维单纯形,则计算其边界。
A3:更新Delaunay三角剖分结构,同时也更新当前的表面网格结构,即把新采样点插入到当前的表面网格中;
如果当前网格中被替换的表面面片集合是正则的二维单纯形,在计算被替换面片集合的边界后,以深度优先的方式将当前待插入点用增量式算法插入到Delaunay剖分结构中,在计算剖分的过程中,如果发现有新生成的四面体有一条边落在被替换表面集合的边界中,则生成一个表面面片。
A4:利用局部膨胀和雕刻的方法优化当前重建的表面网格;
由于A3的替换策略,会产生违背Gabriel准则的三角面片。因此,在满足不改变二维流形的情况下运用局部膨胀和雕刻,用其滤除表面网格中的新生成的非Gabriel三角面片。对于每个新生成的四面体,局部膨胀是把非Gabriel表面面片关联的Delaunay四面体移到当前重建表面网格的内部,而局部雕刻是把非Gabriel表面面片关联的Delaunay四面体移到当前表面网格的外部。
A5:插入已经在Delaunay剖分结构但不在表面网格中的孤立点,以再次提高表面网格的质量。
孤立点是指在Delaunay剖分中,但不在表面网格之中的采样点。插入孤立点,就是在当前表面网格中,通过距离孤立点最近的表面面片的领域,计算不满足Gabriel准则的表面面片的集合,并且这个集合应该为正则二维单纯形。然后,将孤立点与这个集合的边界生成表面面片。
以上内容是结合具体/优选的实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,其还可以对这些已描述的实施方式做出若干替代或变型,而这些替代或变型方式都应当视为属于本发明的保护范围。

Claims (7)

1.一种用于心脏标测系统中的心内膜逐点插入动态表面重建方法,其特征在于,所述方法包括如下步骤:
A1:导管采集四个仿射独立的点,并通过上述四个仿射独立的点初始化为一个四面体网格;
A2:插入新采样点,并且根据新插入的采样点计算当前重建的表面网格中需要被替换的表面三角面片集合;
A3:更新Delaunay三角剖分结构,同时也更新当前的表面网格结构,即把新采样点插入到当前的表面网格中;
A4:利用局部膨胀和雕刻的方法优化当前重建的表面网格;
A5:插入已经在Delaunay剖分结构但不在表面网格中的孤立点,以再次提高表面网格的质量。
2.根据权利要求1所述的方法,其特征在于,步骤A2包括:当新插入一个采样点到当前重建的心内膜表面网格中,通过一个种子三角面片,以广度优先的方式进行增长,从而求得需要被替换的表面三角面片集合。
3.根据权利要求2所述的方法,其特征在于,种子三角面片为选择距离新待插入点最近的表面三角面片,其计算步骤如下:
A21:通过点定位计算当前插入点的四面体单元;
A22:以当前四面体单元为种子,以区域增长的方式求得该种子领域内距离新的待插入点最近的表面三角面片作为种子三角面片。
4.根据权利要求3所述的方法,其特征在于,计算得到种子表面三角面片后,通过广度优先的方式计算当前表面网格中被替换的表面面片集合,其计算步骤如下:
以区域增长的方式检测种子面片相邻的表面网格上的面片,检测这些面片是否满足条件:在待插入点的Delaunay影响域中,或者不在Delaunay影响域中但不满足Gabriel准则;对于满足条件的表面面片一律被标记为需要被替换的表面面片;如果被替换的表面面片集合不是一个正则的二维单纯形,则不插入当前采样点并重新采样;如果被替换的表面面片集合是一个正则的二维单纯形,则计算其边界。
5.根据权利要求1所述的方法,其特征在于,步骤A3包括:如果当前网格中被替换的表面面片集合是正则的二维单纯形,在计算被替换面片集合的边界后,以深度优先的方式将当前待插入点用增量式算法插入到Delaunay剖分结构中,在计算剖分的过程中,如果发现有新生成的四面体有一条边落在被替换表面集合的边界中,则生成一个表面面片。
6.根据权利要求1所述的方法,其特征在于,步骤A4括:在满足不改变二维流形的情况下运用局部膨胀和雕刻,用其滤除表面网格中的新生成的非Gabriel三角面片,对于每个新生成的四面体,局部膨胀是把非Gabriel表面面片关联的Delaunay四面体移到当前重建表面网格的内部,而局部雕刻是把非Gabriel表面面片关联的Delaunay四面体移到当前表面网格的外部。
7.根据权利要求1所述的方法,其特征在于,步骤A5包括:插入孤立点,就是在当前表面网格中,通过距离孤立点最近的表面面片的领域,计算不满足Gabriel准则的表面面片的集合,并且这个集合应该为正则二维单纯形,然后,将孤立点与这个集合的边界生成表面面片。
CN201810778217.3A 2018-07-16 2018-07-16 一种用于心脏标测系统中的心内膜逐点插入动态表面重建方法 Active CN108961413B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810778217.3A CN108961413B (zh) 2018-07-16 2018-07-16 一种用于心脏标测系统中的心内膜逐点插入动态表面重建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810778217.3A CN108961413B (zh) 2018-07-16 2018-07-16 一种用于心脏标测系统中的心内膜逐点插入动态表面重建方法

Publications (2)

Publication Number Publication Date
CN108961413A true CN108961413A (zh) 2018-12-07
CN108961413B CN108961413B (zh) 2023-07-04

Family

ID=64481350

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810778217.3A Active CN108961413B (zh) 2018-07-16 2018-07-16 一种用于心脏标测系统中的心内膜逐点插入动态表面重建方法

Country Status (1)

Country Link
CN (1) CN108961413B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110153582A (zh) * 2019-05-09 2019-08-23 清华大学 焊接方案生成方法、装置以及焊接系统
CN110428386A (zh) * 2019-06-25 2019-11-08 口口相传(北京)网络技术有限公司 地图网格合并方法、装置、存储介质、电子装置

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2278380A1 (en) * 1998-07-24 2000-01-24 Biosense, Inc. Three-dimensional reconstruction of intrabody organs
US6278457B1 (en) * 1998-01-15 2001-08-21 International Business Machines Corporation Methods and apparatus for performing sampling based synthesis of three-dimensional geometric models
US20030067461A1 (en) * 2001-09-24 2003-04-10 Fletcher G. Yates Methods, apparatus and computer program products that reconstruct surfaces from data point sets
US20080218509A1 (en) * 2007-03-09 2008-09-11 Voth Eric J Method and system for repairing triangulated surface meshes
CN102999937A (zh) * 2011-09-08 2013-03-27 上海翰攀信息科技有限公司 心脏散乱点云数据曲面重建的方法
US20140257083A1 (en) * 2011-10-12 2014-09-11 The Johns Hopkins University Methods for Evaluating Regional Cardiac Function and Dyssynchrony from a Dynamic Imaging Modality Using Endocardial Motion
CN104063903A (zh) * 2014-07-08 2014-09-24 清华大学 三维实体模型的四面体网格生成方法和装置
CN104574472A (zh) * 2014-12-31 2015-04-29 北京大学 基于嵌入网格的固体碎裂模拟和动画方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6278457B1 (en) * 1998-01-15 2001-08-21 International Business Machines Corporation Methods and apparatus for performing sampling based synthesis of three-dimensional geometric models
CA2278380A1 (en) * 1998-07-24 2000-01-24 Biosense, Inc. Three-dimensional reconstruction of intrabody organs
US20010009974A1 (en) * 1998-07-24 2001-07-26 Daniel Reisfeld Three-dimensional reconstruction of intrabody organs
US20030067461A1 (en) * 2001-09-24 2003-04-10 Fletcher G. Yates Methods, apparatus and computer program products that reconstruct surfaces from data point sets
US20080218509A1 (en) * 2007-03-09 2008-09-11 Voth Eric J Method and system for repairing triangulated surface meshes
CN102999937A (zh) * 2011-09-08 2013-03-27 上海翰攀信息科技有限公司 心脏散乱点云数据曲面重建的方法
US20140257083A1 (en) * 2011-10-12 2014-09-11 The Johns Hopkins University Methods for Evaluating Regional Cardiac Function and Dyssynchrony from a Dynamic Imaging Modality Using Endocardial Motion
CN104063903A (zh) * 2014-07-08 2014-09-24 清华大学 三维实体模型的四面体网格生成方法和装置
CN104574472A (zh) * 2014-12-31 2015-04-29 北京大学 基于嵌入网格的固体碎裂模拟和动画方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
肖鹏飞等: "基于快速建立四面体网格的有限元心脏建模", 《计算机仿真》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110153582A (zh) * 2019-05-09 2019-08-23 清华大学 焊接方案生成方法、装置以及焊接系统
CN110153582B (zh) * 2019-05-09 2020-05-19 清华大学 焊接方案生成方法、装置以及焊接系统
CN110428386A (zh) * 2019-06-25 2019-11-08 口口相传(北京)网络技术有限公司 地图网格合并方法、装置、存储介质、电子装置

Also Published As

Publication number Publication date
CN108961413B (zh) 2023-07-04

Similar Documents

Publication Publication Date Title
EP2886047B1 (en) Dynamic feature rich anatomical reconstruction from a point cloud
CN108961413A (zh) 一种用于心脏标测系统中的心内膜逐点插入动态表面重建方法
Hutchinson et al. An organized approach to the localization, mapping, and ablation of outflow tract ventricular arrhythmias
CN104809730B (zh) 从胸部ct图像提取气管的方法和装置
CN102422335B (zh) 用于交互式术前评估的系统、方法和装置
CN100589775C (zh) 一种重建人体器官三维表面的方法及系统
Sarhadi et al. Nerve supply of the breast with special reference to the nipple and areola: Sir Astley Cooper revisited
CN109715055A (zh) 具有效率算法的心脏标测系统
CN110464459A (zh) 基于ct-mri融合的介入计划导航系统及其导航方法
US10285611B1 (en) Cardiac mapping using a 3D grid
CA3015087A1 (en) Mesh fitting algorithm
CN110383345A (zh) 用于内腔导航的扁平化视图
TW201109001A (en) Method of fabricating artificial implant
EP4046574B1 (en) Automatic mesh reshaping of an anatomical map to expose internal points of interest
CN105935314B (zh) 一种具有多点标测功能的射频消融导管
CA3055923A1 (en) Iterative coherent mapping of cardiac electrophysiological (ep) activation including scar effects
Jiang et al. An anatomical study of the full-length phrenic nerve and its blood supply: clinical implications for endoscopic dissection
EP2728552A2 (en) Combining three-dimensional surfaces
CN105796094A (zh) 一种基于ecgi的室性早搏异常激动点定位方法
Yang et al. Chengia laxispicata gen. et sp. nov., a new ephedroid plant from the Early Cretaceous Yixian Formation of western Liaoning, Northeast China: evolutionary, taxonomic, and biogeographic implications
CN101443792A (zh) 用于将电生理学信息映射到复杂几何形状上的系统和方法
Bames Reduction of massive breast hypertrophy
CN112043378A (zh) 对人进行导航支持以关于切除部分进行导航的方法和系统
IL274726B2 (en) Volumetric lat map
CN108338785A (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