CN110352447A - 临床图像中的基准检测 - Google Patents
临床图像中的基准检测 Download PDFInfo
- Publication number
- CN110352447A CN110352447A CN201880014983.0A CN201880014983A CN110352447A CN 110352447 A CN110352447 A CN 110352447A CN 201880014983 A CN201880014983 A CN 201880014983A CN 110352447 A CN110352447 A CN 110352447A
- Authority
- CN
- China
- Prior art keywords
- filter
- image
- clinical image
- reference mark
- clinical
- 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
Links
- 238000000034 method Methods 0.000 claims abstract description 92
- 230000002708 enhancing effect Effects 0.000 claims description 59
- 238000003384 imaging method Methods 0.000 claims description 55
- 238000001914 filtration Methods 0.000 claims description 47
- 238000012545 processing Methods 0.000 claims description 26
- 238000009792 diffusion process Methods 0.000 claims description 22
- 239000011159 matrix material Substances 0.000 claims description 22
- 230000009466 transformation Effects 0.000 claims description 19
- 238000006243 chemical reaction Methods 0.000 claims description 16
- 230000004044 response Effects 0.000 claims description 15
- 230000003044 adaptive effect Effects 0.000 claims description 11
- 238000004458 analytical method Methods 0.000 claims description 9
- 230000002526 effect on cardiovascular system Effects 0.000 claims description 6
- 238000001514 detection method Methods 0.000 claims description 5
- 238000012163 sequencing technique Methods 0.000 claims description 5
- 238000000605 extraction Methods 0.000 claims description 3
- 210000000038 chest Anatomy 0.000 description 23
- 230000000007 visual effect Effects 0.000 description 14
- 238000010276 construction Methods 0.000 description 11
- 238000003745 diagnosis Methods 0.000 description 11
- 230000008859 change Effects 0.000 description 8
- 239000000463 material Substances 0.000 description 8
- 230000005855 radiation Effects 0.000 description 8
- 230000008569 process Effects 0.000 description 7
- 230000000712 assembly Effects 0.000 description 6
- 238000000429 assembly Methods 0.000 description 6
- 238000002001 electrophysiology Methods 0.000 description 6
- 230000007831 electrophysiology Effects 0.000 description 6
- 230000004048 modification Effects 0.000 description 6
- 238000012986 modification Methods 0.000 description 6
- 210000004165 myocardium Anatomy 0.000 description 6
- 229910002056 binary alloy Inorganic materials 0.000 description 5
- 239000011248 coating agent Substances 0.000 description 4
- 238000000576 coating method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000002594 fluoroscopy Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000001737 promoting effect Effects 0.000 description 4
- 238000002679 ablation Methods 0.000 description 3
- 210000000988 bone and bone Anatomy 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000000284 extract Substances 0.000 description 3
- 238000001727 in vivo Methods 0.000 description 3
- 238000013507 mapping Methods 0.000 description 3
- 238000006467 substitution reaction Methods 0.000 description 3
- 238000013519 translation Methods 0.000 description 3
- 230000000747 cardiac effect Effects 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000005684 electric field Effects 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000002595 magnetic resonance imaging Methods 0.000 description 2
- 230000000877 morphologic effect Effects 0.000 description 2
- 238000009206 nuclear medicine Methods 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 238000002600 positron emission tomography Methods 0.000 description 2
- 238000009877 rendering Methods 0.000 description 2
- 238000001931 thermography Methods 0.000 description 2
- 230000002792 vascular Effects 0.000 description 2
- DIDGPCDGNMIUNX-UUOKFMHZSA-N 2-amino-9-[(2r,3r,4s,5r)-5-(dihydroxyphosphinothioyloxymethyl)-3,4-dihydroxyoxolan-2-yl]-3h-purin-6-one Chemical compound C1=2NC(N)=NC(=O)C=2N=CN1[C@@H]1O[C@H](COP(O)(O)=S)[C@@H](O)[C@H]1O DIDGPCDGNMIUNX-UUOKFMHZSA-N 0.000 description 1
- 238000012935 Averaging Methods 0.000 description 1
- RTAQQCXQSZGOHL-UHFFFAOYSA-N Titanium Chemical compound [Ti] RTAQQCXQSZGOHL-UHFFFAOYSA-N 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 239000005030 aluminium foil Substances 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 238000002583 angiography Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000002146 bilateral effect Effects 0.000 description 1
- 210000004413 cardiac myocyte Anatomy 0.000 description 1
- 210000000748 cardiovascular system Anatomy 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000002224 dissection Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000001839 endoscopy Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 229910001092 metal group alloy Inorganic materials 0.000 description 1
- 238000007500 overflow downdraw method Methods 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000012925 reference material Substances 0.000 description 1
- 235000002020 sage Nutrition 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000001225 therapeutic effect Effects 0.000 description 1
- 239000010936 titanium Substances 0.000 description 1
- 229910052719 titanium Inorganic materials 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
- G06T7/73—Determining position or orientation of objects or cameras using feature-based methods
- G06T7/75—Determining position or orientation of objects or cameras using feature-based methods involving models
-
- 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
- A61B90/00—Instruments, 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/39—Markers, e.g. radio-opaque or breast lesions markers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration by the use of local operators
-
- G06T5/70—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B18/00—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body
- A61B18/04—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body by heating
- A61B18/12—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body by heating by passing a current through the tissue to be heated, e.g. high-frequency current
- A61B18/14—Probes or electrodes therefor
- A61B18/1492—Probes or electrodes therefor having a flexible, catheter-like structure, e.g. for heart ablation
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B90/00—Instruments, 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/36—Image-producing devices or illumination devices not otherwise provided for
- A61B90/37—Surgical systems with images on a monitor during operation
- A61B2090/376—Surgical systems with images on a monitor during operation using X-rays, e.g. fluoroscopy
- A61B2090/3762—Surgical systems with images on a monitor during operation using X-rays, e.g. fluoroscopy using computed tomography systems [CT]
- A61B2090/3764—Surgical systems with images on a monitor during operation using X-rays, e.g. fluoroscopy using computed tomography systems [CT] with a rotating C-arm having a cone beam emitting source
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B90/00—Instruments, 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/39—Markers, e.g. radio-opaque or breast lesions markers
- A61B2090/3966—Radiopaque markers visible in an X-ray image
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B90/00—Instruments, 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/39—Markers, e.g. radio-opaque or breast lesions markers
- A61B2090/3983—Reference marker arrangements for use with image guided surgery
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/25—Bioelectric electrodes therefor
- A61B5/279—Bioelectric electrodes therefor specially adapted for particular uses
- A61B5/28—Bioelectric electrodes therefor specially adapted for particular uses for electrocardiography [ECG]
- A61B5/283—Invasive
- A61B5/287—Holders for multiple electrodes, e.g. electrode catheters for electrophysiological study [EPS]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61M—DEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
- A61M25/00—Catheters; Hollow probes
- A61M25/01—Introducing, guiding, advancing, emplacing or holding catheters
- A61M25/0105—Steering means as part of the catheter or advancing means; Markers for positioning
- A61M2025/0166—Sensors, electrodes or the like for guiding the catheter to a target zone, e.g. image guided or magnetically guided
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20004—Adaptive image processing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20016—Hierarchical, coarse-to-fine, multiscale or multiresolution image processing; Pyramid transform
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30021—Catheter; Guide wire
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30204—Marker
Abstract
本发明涉及用于识别患者的解剖区域的一个或多个图像上的基准标记的系统和方法。在识别基准标记之前,识别并去除其它对象。特别地,一些实施例涉及在荧光透视图像中存在导管和其它医疗设备的情况下检测基准标记。在这种实施例中,从临床图像中识别和去除医疗设备有助于正确识别基准标记。
Description
相关申请的交叉引用
本申请要求2017年3月9日提交的美国临时申请号62/469,423的优先权,其通过引用并入本文,如同在本文完全阐述一样。
技术领域
本发明涉及医疗系统,诸如医疗设备导航系统。具体地,本发明涉及一种用于在诸如医疗设备导航系统的医疗系统的坐标系中配准患者的解剖区域的一个或多个图像的系统。
背景技术
期望当医疗设备在体内移动时跟踪医疗设备的位置,使得例如在适当的位置施用药物和其它形式的治疗,并且可以更有效和安全地完成医疗过程。跟踪医疗设备在体内的位置的一种传统手段是荧光透视成像。然而,在一些应用中,荧光透视可能是不受欢迎的,因为它使患者和临床医生受到高水平的电磁辐射。因此,已经开发了替代的医疗设备导航系统来跟踪医疗设备在体内的位置。这些系统通常依赖于电场或磁场的生成以及连接到医疗设备和/或身体外部的位置传感器上的感应电压和电流的检测。然后,通过例如视觉显示器将从这些系统获得的信息提供给临床医生。通常,医疗设备的表现形式相对于正在被操作的设备所在的解剖区域的计算机模型或者一个或多个图像(包括但不限于荧光透视图像)显示。为了将医疗设备显示在相对于模型或图像的正确位置,可以在导航系统的坐标系内配准模型或图像。
可以以多种方式将图像配准在医疗设备导航系统的坐标系中。用于捕获图像的成像系统可以与导航系统物理集成,诸如由圣犹达医疗公司销售的MediGuideTM系统,或者如在公开为US 2008/0183071的共同转让的美国专利申请序列号11/971,004中所述,其全部公开内容通过引用并入本文。在用于捕获图像的成像系统与导航系统物理集成的情况下,成像系统可以在安装期间与导航系统配准,并且导航系统与成像系统的空间关系此后是恒定且已知的,从而避免了需要在每个新过程期间配准。然而,在导航系统和成像系统在物理上分离的情况下,系统的变化的空间关系使得配准更加复杂。一种解决方案是在成像系统的视野中在相对于导航系统的坐标系已知的位置处放置多个不透射线的基准标记。因为标记在由成像系统产生的图像中是可见的,所以可以通过将图像中的每个标记的可见位置与导航坐标系中的标记的已知位置进行协调来使图像与坐标系配准。然而,该解决方案的一个缺点是基准标记通常在执行该过程的临床医生看到的图像中可见,可能干扰临床医生对被成像的解剖区域的视图。
发明人在此认识到需要一种用于在医疗设备导航系统的坐标系中配准患者的解剖区域的一组图像的系统和方法,其将最小化和/或消除一个或多个以上确定的缺陷。前述讨论仅旨在说明本领域,而不应视为对权利要求范围的否定。
发明内容
本发明涉及用于在诸如医疗设备导航系统的医疗系统的坐标系中配准患者的解剖区域的一个或多个图像的系统和方法;特别地,在一个或多个图像中存在导管和其它医疗设备的情况下检测基准标记。
本发明的各方面涉及使用诸如基准标记的对象来促进在物理上分离的医疗系统(诸如医疗设备导航系统)的坐标系中的图像配准的系统和方法。结果,向临床医生呈现解剖区域和相关关注项目的更准确和详细的图像。
本发明的一些实施例涉及一种用于在临床图像上识别和定位基准标记的方法。该方法包括从成像系统接收临床图像,处理临床图像以去除假基准标记,以及识别临床图像中的基准标记。在一些更具体的实施例中,该方法进一步包括将增强滤波器应用于临床图像以增强基准标记的外观。
本发明的各方面涉及一种检测临床图像内的医疗设备电极的方法。该方法包括从成像系统接收临床图像,将增强滤波器应用于临床图像以增强电极的外观,去除临床图像的背景,以及识别临床图像内的电极。在更详细的实施例中,识别临床图像内的电极的步骤包括自适应阈值处理、基于大小和形状的分量滤波,以及电极的提取。
本发明的其它实施例涉及一种用于将医疗设备的位置数据从导航系统叠加到临床图像上的方法。该方法包括:从临床成像系统接收临床图像;处理临床图像以去除假基准标记;识别临床图像内的基准标记;创建变换模型,该变换模型将协调导航系统的第一坐标系与临床成像系统的第二坐标系;确定第一坐标系中医疗设备的位置;将变换模型应用于第一坐标系中的医疗设备的位置,以确定第二坐标系中医疗设备的位置;以及基于第二坐标系中医疗设备的已知位置,将医疗设备的图像叠加到临床图像上。
又一些实施例涉及用于心血管导管的导航系统。该导航系统包括心血管导管、光磁配准板、临床成像系统、导管定位系统和控制器电路。心血管导管包括位于导管远侧末端附近的一个或多个电极。电极有助于导管远侧末端的定位。光磁配准板包括多个基准标记。临床成像系统在第一坐标系中显示包括患者、基准标记和心血管导管的临床图像。导管定位系统在第二坐标系中检测一个或多个电极的位置和取向。控制器电路通信地耦合到临床成像系统和导管定位系统。基于第二坐标系中的基准标记的已知位置和临床图像内的基准标记的位置,控制器电路确定第一和第二坐标系之间的变换模型。基于变换模型,控制器电路确定第一坐标系内电极的位置。在更具体的实施例中,导航系统可以包括显示器,其中控制器电路通信地耦合到显示器并且进一步基于所确定的第一坐标系中的电极的位置,生成用于显示器的信号,该信号包括叠加有导管的代表性图像的临床图像。
通过阅读以下描述和权利要求以及阅读附图,本发明的前述和其它方面、特征、细节、效用和优点将变得显而易见。
附图说明
考虑到以下结合附图的详细描述,可以更全面地理解多种示例实施例,在附图中:
图1A是与本发明的多种实施例一致的用于在医疗系统的坐标系中配准一组图像的系统的一个实施例的示意图。
图1B是与本发明的多种实施例一致的用于在医疗系统的坐标系中配准一组图像的系统的另一实施例的示意图。
图1C是示出与本发明的多种实施例一致的用于在医疗系统的坐标系中配准一组图像的方法的另一实施例的流程图。
图1D是与本发明的多种实施例一致的用于基准检测和背景过滤的流程图。
图2A是与本发明的多种实施例一致的患者胸部的荧光透视图像,其具有3个微创诊断/治疗导管。
图2B是与本发明的多种实施例一致的图2A的荧光透视图像的滤波版本,在图2A中图像中的基准已经增强。
图3是与本发明的多种实施例一致的其中具有微创诊断/治疗导管的患者胸部的荧光透视图像,该荧光透视图像采用非线性扩散反应滤波器过滤以去除背景图像。
图4是与本发明的多种实施例一致的其中具有微创诊断/治疗导管的患者胸部的荧光透视图像,该荧光透视图像用海森(Hessian)滤波器过滤以增强图像中的基准。
图5是与本发明的多种实施例一致的图4的荧光透视图像,其中突出显示了聚类基准元素并且使用测地距离抑制背景。
图6A是与本发明的多种实施例一致的患者胸部的荧光透视图像,在该胸部中具有微创诊断/治疗导管以及患者胸部上的磁定位系统和患者参考传感器的线束。
图6B是与本发明的多种实施例一致的在自适应阈值处理之后应用海森滤波器的图6A的荧光透视图像。
图7A是与本发明的多种实施例一致的患者胸部的荧光透视图像,在该胸部中具有微创诊断/治疗导管,其中通过算法识别真和假基准。
图7B是与本发明的多种实施例一致的在过滤之后具有仅识别的真实基准的图7A的荧光透视图像。
虽然本文讨论的各种实施例适合于修改和替代形式,但是其各方面已经通过附图中的示例示出并且将被详细描述。然而,应该理解,意图不是将本发明限制于所描述的特定实施例。相反,意图是覆盖落入本发明范围内的所有修改、等同物和替代物,包括权利要求中限定的方面。另外,在整个本申请中使用的术语“示例”仅是为了说明而非限制。
具体实施方式
本发明涉及用于在医疗系统(诸如医疗设备导航系统)的坐标系中配准患者的解剖区域的一个或多个图像的系统和方法。特别地,在一个或多个图像中存在导管和其它医疗设备的情况下检测基准标记。下面具体参考附图描述本发明的多种实施例的细节。
在多种实施例中,诸如电磁跟踪和导航系统的医疗导航系统可提供在心脏导管工作站(Cathlab)的X射线图像(在本文中也称为荧光透视图像)上呈现从磁传感器获得的磁信息的能力。为了促进这一点,可以计算从三维(3D)世界到二维(2D)X射线图像的投影矩阵(本文也称为变换矩阵)。投影矩阵可以使用在2D X射线图像上可见并且还与例如磁坐标系内的已知位置相关联的一组不透射线的元素(基准)来计算。不透射线的基准可以位于光磁配准板(OMRP)上,并且作为与本发明的多种实施例一致的电磁导航系统的一部分。
在临床干预和诊断探查期间捕获的荧光透视图像通常包括医疗设备,诸如导管和传感器-其可以包含视觉上看起来类似于OMRP基准的不透射线元素。在荧光透视图像中存在这些“假”类似基准元素可能导致基准检测准确性降低。基准检测准确性对计算的投影矩阵的准确性具有显著影响,并且因此对整个医疗导航系统的准确性具有显著影响。
一种示出患者的解剖区域中的医疗设备的方法是将医疗设备的实时渲染的表现形式叠加到解剖区域的一组图像上。该方法有利地允许临床医生操作医疗设备以查看设备相对于解剖区域的位置,而不会使患者和临床医生暴露于过度的成像辐射。图像由成像系统捕获,并且导航系统基本上实时地跟踪医疗设备的位置。编程的电子控制单元创建用于显示的合成图像,该合成图像包括叠加在一组图像上的跟踪医疗设备的表现形式。为了创建合成图像,可以首先将图像配准在导航系统的坐标系中-即,可以使组中的每个图像的坐标系与导航系统的坐标系协调。
本发明的各方面涉及用于在存在由出现在荧光透视图像中的医疗设备(诸如导管和传感器)引入的其它不相关的不透射线元素的情况下稳健检测位于OMRP上的不透射线的基准球的算法。本发明检测OMRP基准,同时使用多种方法排除假类似基准元素,以从荧光透视图像中提取和过滤球形不透射线元素。OMRP板上的基准有助于将医疗设备的实时渲染的表现形式叠加到医疗设备正在其中操作的解剖区域的一组图像上。
虽然本发明的多种实施例涉及诸如在荧光透视图像上可见的不透射线元素的基准标记,但是本发明的各方面可以应用于非荧光透视成像并且完全不限于OMRP类型的基准。相反,本发明的各方面可以应用于多种医学成像技术,包括超声波检查、磁共振成像、内窥镜检查、弹性成像、触觉成像、热成像、核医学功能成像技术,包括正电子发射断层成像和单光子发射计算机断层成像等(以下称为成像系统,其产生用于与本发明一致的图像处理的临床图像)。
应当理解,虽然本发明的多种实施例涉及包括OMRP的基准标记的临床图像,但是根据本发明的图像处理可以被实施为没有基准标记(诸如在OMRP内的基准标记)的临床图像。在这种实施例中,基准可以是例如标测导管或消融导管的电极。此外,除了临床医生可能关注的图像内的任何其它对象之外,本文公开的图像处理算法和过滤技术可用于检测其它基准标记,诸如导管、线、其它电极、引线、静脉口、病变、导管末端、外部钛球。在一个特定实施例中,例如,本文公开的图像处理系统可识别两个导管在临床图像内彼此交叉的位置。
日期为2011年12月25日的美国专利公开号2013/0272592阐述了用于在医疗系统的坐标系中配准荧光透视图像的系统和方法,并且通过引用并入本文,如同在本文完全公开一样。
图1A示出了成像和导航系统10,其用于对患者12(诸如心脏14)的解剖区域进行成像,并用于确定解剖区域内的医疗设备16的位置和导航。设备16可以包括例如电生理学(EP)标测导管、心腔内超声心动图(ICE)导管、消融导管等。然而,应该理解,所公开的系统可以用于成像心脏14以外的各种解剖区域,并且可以取决于关注的解剖区域与多种诊断和治疗设备相结合。系统10包括成像系统18和医疗设备导航系统20。根据本发明,系统10还包括用于在导航系统的坐标系22中配准患者12的解剖区域的一组图像的配准系统。配准系统可包括一个或多个对象24和电子控制单元(ECU)26。
提供成像系统18以获取心脏14或其他相关解剖区域的图像,并且在所示实施例中包括临床成像系统。尽管在该实施例中描述了临床成像系统,但是本文描述的发明可以用于配置成捕获一组图像的其它类型的成像系统,该图像包括例如但不限于计算机断层扫描(CT)成像系统和三维旋转(radio)血管造影(3DRA)系统等。系统18可包括C形臂支撑结构28、辐射发射器30和辐射检测器32。发射器30和检测器32设置在支撑结构28的相对端上,并当患者12放置在手术台34上时设置在患者12的相对侧上。发射器30和检测器32限定视野36,并且定位成使得当患者12放置在手术台34上时,视野36包括关注的解剖区域。成像系统18被配置为捕获视野36内的解剖特征和其它对象的图像。支撑结构28可以具有围绕患者旋转的自由度,如由线38、40所示。支撑结构28也可以具有沿线42、44(即,沿患者12的颅尾轴线)和/或沿线46、48(即垂直于患者12的颅尾轴线)滑动的自由度。支撑结构28的旋转和平移运动产生对应的视野36的旋转和平移运动。在另外的实施例中,包括成像和导航系统10的电生理学实验室(EP实验室)可以具有多种其它运动自由度;例如,手术台34还能够平移、源-增强器距离(“SID”)、颅尾旋转和检测器32绕其轴线旋转。可以使用本文公开的成像和导航系统来解决这些多种其它运动自由度。
成像系统18可以通过首先沿着线42、44、46、48移动,将关注的解剖区域放置在视野36内,以获取患者12的解剖区域的一组图像。其次,支撑结构28可以围绕患者12旋转辐射发射器30和辐射检测器32,将解剖区域保持在视野36内。当支撑结构28旋转时,成像系统18可以捕获解剖区域的图像,从各种角度提供一组解剖区域的二维图像。该组图像可以被传送到ECU 26以进行图像处理并在显示器58上显示。该组图像可以包括在预定的时间段内拍摄的一系列图像。
提供导航系统20以确定医疗设备16在患者12的身体内的位置,并允许临床医生在身体内导航设备16。在所示实施例中,系统20包括磁导航系统,其在解剖区域中生成磁场,并且与设备16相关联的位置传感器生成响应于传感器在磁场内的位置而改变的输出。系统20可以包括例如由圣犹达医疗公司以商标“GMPS”出售并且通常在例如美国专利号6,233,476、7,197,354和7,386,339中示出和描述的系统,其全部公开内容通过引用并入本文,如同在本文完全公开一样;或者由生物传感韦伯斯特公司(Biosense Webster,Inc)以商标“CARTO XP”出售并且通常例如在美国专利号5,391,199、5,443,489、5,558,091、6,498,944、6,788,967和6,690,963中示出和描述的系统,其全部公开内容通过引用并入本文,如同在本文完全公开一样。尽管在所示实施例中示出了磁导航系统,但是应该理解,本发明可以用于多种导航系统,包括基于轴特定电场的创建和检测的导航系统,诸如用于由圣犹达医疗公司以商标“ENSITE NAVX”销售所提供的系统。导航系统20可包括发射器组件50和一个或多个位置传感器52以及电子控制单元(ECU)54。
发射器组件50在本领域中以其传统意义,并且可包括彼此正交布置以在关注的解剖区域中和/或周围产生磁场的多个线圈。应该注意的是,尽管发射器组件50在图1A中在患者12的身体下方并且在桌子34下方示出,但是发射器组件50可以放置在其它位置中,诸如附接到辐射发射器30或辐射检测器32,磁场发生器可以从该位置将磁场投射到关注的解剖区域中。根据本发明的某些实施例,发射器组件50位于视野36内。在一些特定实施例中,辐射发射器30定位在患者12下方,并且对应的辐射检测器32位于患者上方。
位置传感器52被配置成取决于由发射器组件50生成的场内的传感器52的相对位置而产生输出。传感器52与设备16处于已知的位置关系,并且可以附接到医疗设备16。在图1A中,位置传感器52和医疗设备16示为设置在心脏14内。位置传感器52可以附接到设备16的远端或近端,或其间的任何点。当医疗设备16被引导并穿过关注的解剖区域时,导航系统20确定位置传感器52在所生成的区域中的位置,并因此确定医疗设备16的位置。
提供ECU 54以控制发射器组件50磁场的生成,并处理从传感器52接收的信息。ECU54可包括可编程微处理器或微控制器,或者可包括专用集成电路(ASIC)。ECU 54可以包括中央处理单元(CPU)和输入/输出(I/O)接口。ECU 54经由I/O接口可以接收包括由传感器52生成的信号的多个输入信号,并生成多个输出信号,包括用于控制发射器50的输出信号。尽管在所示实施例中ECU 54被示为单独的部件,但是应该理解,ECU 54和ECU 26可以集成到单个单元中。
提供对象24以允许将由成像系统18捕获的图像配准在导航系统20的坐标系22中。每个对象24可包括例如一个或多个基准标记。对象24定位于成像系统18的视野36内。对象24也位于坐标系22内的已知位置,或者连接到可以在坐标系22内提供位置的传感器。在所示的实施例中,例如,对象24被放置在导航系统20的位于视野36内的部件(诸如发射器50)上。可替代地,对象24可以放置在距导航系统20的部件(例如发射器50)固定的距离处,但仍然在视野36内。对象24也可以固定到设置在视野36内的导管的远端部分,或固定到视野36内的患者的身体部分。可替代地,对象24可以包括类似于可以在导航系统20中定位的传感器52的传感器(未示出)。对象24可以随时间推移保持在坐标系22内的相同已知位置中,或者可以在坐标系22中的不同已知位置之间移动(例如,后续位置与坐标系22内的先前已知位置具有已知关系)。在对象24包括多个基准标记的情况下,标记可以以预定图案布置或以其中标记彼此具有已知关系的方式布置。标记也可具有不同程度的射线不透性。
本发明的各方面涉及在过程开始时、在各个阶段,或者取决于临床医生的需要接近连续地将图像系统18的图像配准在导航系统20的坐标系22中。
参考图1B,在本发明的另一个实施例中,配准系统包括用于在设置在视野36内的位置(使得基准标记25处于第一状态并且在由成像系统18捕获的图像中可见)和设置在视野36外部的另一个位置(使得基准标记25处于第二状态并且在由成像系统18捕获的图像中不可见)之间移动基准标记25的装置56。移动设备56可以包括任何传统的机械设备,该设备通过在ECU 26的控制下操作的线性、旋转或其它运动(例如,直线驱动器或电机驱动平台)在位置之间移动基准标记25。在多种实施例中,基准标记25可以包括一个或多个对象24。在这种实施例中,对象24可以经由OMRP在空间上彼此耦合,从而限制对象24相对于彼此的独立运动。
在本发明的又一个实施例中,对象24配置有相对低强度和/或相对渐变的强度变化,使得对象24在由成像系统18生成的图像中对人眼基本上不可见,但是通过如下所述的ECU 26的图像处理可检测。对象24可以例如相对较薄或者可以具有低但存在的射线不透性。在一个实施例中,对象24由铝箔或铝片制成。然而,应该理解的是,对象24可以由多种金属、金属合金或具有至少最小水平的射线不透性的其它材料制成。在另一个实施例中,利用人眼难以检测相对小的变化的事实,从对象24上的一点移动到对象24上的另一点,对象24的射线不透性逐渐增大或减小(例如,射线不透性在对象24的中心附近较大,但当远离对象24的中心时逐渐减少)。对象24可以呈现多种大小和形状中的任何一种。
提供ECU 26用于处理由成像系统18生成的图像,并基于图像中的对象24的外观将图像配准在坐标系22中。ECU 26还可以控制医疗设备16、成像系统18、导航系统20和/或显示器58的操作。ECU 26可以包括可编程微处理器或微控制器,或者可以包括专用集成电路(ASIC)。ECU 26可以包括中央处理单元(CPU)和输入/输出(I/O)接口,ECU 26可以通过该接口接收多个输入信号,包括由医疗设备16、成像系统18和导航系统20生成的信号。此外,ECU26可以生成多个输出信号,包括用于控制医疗设备16、成像系统18、导航系统20和显示器58的输出信号。ECU 26还可以从诸如ECG监视器的器官监视器(未示出)接收输入信号,并基于被监视器官的定时信号对来自成像系统18的图像进行分类或分离。例如,ECU 26可以基于所收集每个图像的患者心动周期的相位对图像进行分类,如在例如美国专利号7,697,973中更全面描述的,其全部内容通过引用并入本文,如同在本文完全公开一样。
根据本发明的一个实施例,ECU 26配置有适当的编程指令或代码(即,软件)以执行用于在导航系统20的坐标系22中配准患者心脏14或患者12的另一解剖区域的一组图像的方法中的若干步骤。
现在参考图1C,本发明的各方面涉及将医疗设备的图像叠加到临床图像上的方法。在这种实施例中,第一步骤88包括从由成像系统18生成的一组图像中收集多个图像。该组收集的图像包括包含处于坐标系中的已知位置的一个或多个对象24的区域。取决于对象的强度和ECU的处理能力,由ECU收集的多个图像可以包括整个一组图像或该组图像的子集。在可以确定实际图像位置之前,可以进行基准检测和背景过滤89以识别和去除可能导致错误位置确定的假基准。基准检测和背景过滤89都滤除了假基准,同时还增强了该组图像内的真实基准。一旦识别出真实基准,该方法可以继续确定每个图像中的基准(也称为对象和基准标记)的实际图像位置90。响应于来自每个图像的图像数据的总和,ECU可以识别每个图像中的基准的潜在图像位置。ECU可以使用传统算法来对图像数据进行求和,包括对数据求和和/或求平均,以便消除数据中的噪声,并确定多个图像中的每个图像中的基准的实际图像位置。
该方法可以继续响应于多个图像中的对象的实际图像位置以及对象在坐标系中的已知位置来创建变换模型,步骤94。步骤96使用变换模型将该组图像配准在坐标系22中。最后,步骤98响应于医疗设备在坐标系内的位置,将医疗设备的图像叠加到该组图像中的每个图像上。
图1D是与本发明的多种实施例一致的用于基准检测和背景过滤的流程图100,其可以在相对于图1C的步骤89中进行。下面更详细描述的算法(并且其结果在图2-7中呈现)是参考球形基准来讨论。然而,应该理解,基准检测和背景过滤过程以及特定算法通常易于适用于多种其它形状和基准配置。例如,其它基准布置可包括细长的基准,诸如线和杆。由于导管和其它医疗工具可在临床图像上呈现杆状形状,因此本发明的方面还涉及滤除与这种导管和医疗工具相关的假基准(参见例如图6A-6B及其讨论)。
本文公开的多种基准检测和背景过滤算法可以应用于输入X射线图像(也称为荧光透视图像或临床图像)、一组X射线图像和/或从一系列X射线图像获得的平均X射线图像。基准检测和背景过滤算法还可以应用于超声波检查、磁共振成像、内窥镜检查、弹性成像、触觉成像、热成像、核医学功能成像技术(包括正电子发射断层成像和单光子发射计算机断层成像)等。
在处理器电路(也称为ECU)(参见图2A)接收输入图像101之后,可以进行基准增强102。基准增强102可以利用增强滤波器或其它类似的滤波技术/算法来改善基准元素在X射线图像上的外观(参见图2B)。在多种实施例中,磁定位系统可以利用基准,该基准是在X射线图像上表现为暗点的小金属球。一种基准实施方式利用具有已知基准图案的OMRP。在一些实施例中,基准图案可以是均匀或不均匀的图案。在基准处于非均匀图案的情况下,基准的低错误率识别(如本文所公开的)可能是有益的。基准增强102可以利用一个或多个算法/滤波器来增强基准元素相对于图像中的多种其它特征的外观。可用于促进基准增强102的算法的一些示例包括高斯拉普拉斯算子(“LoG”)、高斯差分(“DoG”)、海森滤波器、方向可控滤波器和非线性扩散反应滤波器(例如,基于贝尔特拉米(Beltrami)的滤波器)。
OMRP还可以利用包括不透射的线和/或小杆的基准。类似于球形基准,例如,由于形状与导管的相似性,在现有技术系统中检测杆状基准也可能导致假基准检测。这些假的杆状基准也可以利用图1D中的基准检测和背景过滤方法来去除这种假基准。
以下公开了临床图像中的基准增强102的若干具体实施例。
在一个特定实施例中,使用LoG算法。LoG是一种基于算子的算法,其中算子基于二阶导数-ΔI=Ixx+Ixy。虽然LoG算法可用于“斑点(blob)”检测,LoG算法通常对噪声敏感。LoG可以通过不同尺度的DoG来近似,从而使得该算子的实施更有效。
在临床图像内的基准增强102的另一特定实施例中,海森滤波器可用于自动检测斑点状结构。基于特征值,对海森矩阵的多尺度分析,可以识别局部图案并将其分类为属于板状、线状或斑点状结构。图像中用于每个像素的图像的海森矩阵是:
其中是图像的平滑二阶导数。图像由高斯西格玛参数 来平滑。然后排序特征值和 然后可以将局部结构(通过对象尺度)和形状辨别(通过分析海森矩阵的特征值)整合到共同的响应测量中。增强斑点多尺度响应是:
增强线结构响应是:
在临床图像中的基准增强102的一些实施例中,可以使用贝尔特拉米滤波器。在一个特定实施例中,贝尔特拉米滤波器是非线性扩散反应滤波器,其被设计在保留图像边缘时过滤图像噪声。如果将反应项加到扩散项上,可以实现图像内基准球的对比度增强。由基准增强102输出的处理图像可以被视为基于进化模型的偏微分方程的渐近状态。一个示例进化模型是:
It=ΔgI+f(I)
其中Δg是拉普拉斯-贝尔特拉米(Laplace-Beltrami)算子(例如,拉普拉斯算子从平坦空间到流形的扩展),并且函数f是反应项。反应项功能选择可以基于多种因素(例如,图像处理应用和后处理要求)而变化。
返回参考图1D,背景去除103可以应用于图像以减少患者的心脏和临床图像内的其它背景阴影的影响。临床图像通常可包括大图像区域,包括暗区域(即阴影)。引起这种较暗区域的特征的常见示例包括骨骼(例如,在患者胸部的临床图像中,肋骨和脊柱通常可见为暗梯度),以及患者的心肌。这些暗区域降低了图像中基准的对比度。因此,由于真实基准、与患者的身体特征相关联的暗区域以及临床图像内的医疗设备之间的最小对比度变化,现有系统更可能识别假基准。背景去除103增加了基准标记和非基准特征之间的图像内的对比度。背景去除103可以使用多种算法,包括例如形态滤波器和/或中值滤波器。中值和形态滤波器可以消除图像上的大对象(诸如阴影),同时保留较小的对象,包括基准标记。
在完成背景去除103之后,执行三个并行过滤路径。第一过滤路径120专用于识别临床图像内的每个基准(也称为质心、对象和基准标记)。第二过滤路径130识别、增强和过滤图像内的任何导管。第三过滤路径140识别、增强和滤除电极(例如,用于电生理学标测和/或消融、磁定位等的导管电极-通常位于导管轴的远端)。
三个并行过滤路径120、130和140中的每一个可以基于它们相应的大小/形状来过滤基准候选分量。例如,第一滤波器路径120消除了在大小上大于或小于典型基准的分量。虽然第一滤波器路径120将基于大小/形状选择真实的基准,但是第一过滤路径120可以忽略较小的分量,诸如由图像噪声生成的分量。
在背景去除103之后,第一过滤路径120可以包括自适应阈值处理104(也称为二进制图计算)。可以通过对来自背景去除103的输出图像进行阈值处理来计算二进制图,以生成二进制的连通分量图像。该分量对应于图像中捕获的暗区域的位置和大小(参见例如图7A-7021-2和7251-N)。一些暗区域由OMRP基准生成,而其它暗区域源自假基准。假基准的一些示例可包括导管、导管电极、用于磁定位系统的患者参考传感器(“PRS”)、医师工具、导丝、可能在图像上显示的其它医疗设备、图像噪声、其它图像处理伪像以及可能导致假基准的出现的造影剂(例如,由于在图像的特定位置处造成对比度)。
第一过滤路径120进一步包括分量滤波器105。分量滤波器105消除在大小上大于或小于典型基准大小的分量。在一些情况下,图像噪声可生成小的分量。然而,这种图像噪声生成的分量小于典型的基准分量并因此被过滤。分量滤波器105进一步消除了具有与典型的基准分量显著偏离的形状的分量。例如,在OMRP由球形基准组成的情况下,细长分量将被分量滤波器105丢弃。
一旦所有假基准已被分量滤波器105过滤,则每个剩余基准与质心(例如,质量中心)相关联,该质心表征每个基准的精确坐标,在本文中称为质心提取106。质心可以被提取为分量的平均值,或者作为分量的加权平均值(即质量中心)。还可以通过在应用一些平滑之前或之后定位分量中的最暗像素来提取质心以克服图像噪声。
第二过滤路径130包括导管候选物107的增强。具体地,导管增强107利用一种或多种算法来增强导管候选物在图像中的多种其它各种特征上的外观。通过增强导管候选物,导管候选物、图像背景和各种其它特征之间的对比度变化显著增加。导管增强107可以利用多种算法,包括增强滤波器或其它类似的滤波技术/算法,以改善导管元件在X射线图像上的外观。可用于促进导管增强107的算法的一些示例包括高斯拉普拉斯算子(“LoG”)、高斯差分(“DoG”)、海森滤波器、方向可控滤波器和非线性扩散反应滤波器。在更具体的实施例中,导管增强107可以使用扩散反应滤波器完成。导管、图像背景和各种其它特征之间增加的对比度有助于从图像中过滤导管候选物108-从图像中消除这些导管候选物进一步消除了在医疗设备(诸如导管和传感器)上或附近的假基准的可能性。
例如,可以在图2A中看到,导管220A-B(即,与周围背景具有不同对比度的双边对象)在图像200上显示为脊。可以使用海森滤波器来应用导管增强107,其对细长结构(诸如导管220A-B)以及其它细长结构(诸如导丝和导线)产生强烈响应。发现可用于导管增强107以及产生与海森滤波器基本上相同的结果的其它脊增强滤波器包括:Frangi滤波器、Oriented滤波器和任何其它脊增强滤波器。
在一些特定实施例中,过滤导管候选物108包括对来自先前步骤107的图像进行阈值处理,以创建二进制图,其中细长结构作为大分量存在。这些大的分量对应于细长工具,诸如图像内的导管和其它医疗设备。虽然本实施例实现自适应阈值,但是其它实施例可以实现恒定阈值以过滤导管候选物。
在第二过滤路径130中,在图像内的导管附近检测到的质心可以被提取为分量的平均值,或者作为分量的加权平均值(即质量中心)。还可以通过在应用一些平滑之前或之后定位分量中的最暗像素来提取质心以克服图像噪声。然后可以排除在检测到的导管附近或检测到的导管上检测到的质心作为基准候选物。
第三过滤路径140包括电极候选物109的增强。具体地,电极增强109利用一种或多种算法来增强电极候选物在图像中的多种其它各种特征上的外观。假基准候选物可以源自在导管(例如,诊断和治疗导管)和传感器(例如,患者参考传感器)上发现的导管电极。通过增强电极候选物,电极候选物、图像背景和各种其它特征之间的对比度显著增加。电极增强109可以利用多种算法,包括增强滤波器或其它类似的滤波技术/算法,以改善电极元件在X射线图像中的对比度。可用于促进电极增强109的算法的一些示例包括高斯拉普拉斯算子(“LoG”)、高斯差分(“DoG”)、海森滤波器、方向可控滤波器和非线性扩散反应滤波器。在更具体的实施例中,可以使用扩散反应滤波器完成电极增强109。电极、图像背景和各种其它特征之间对比度的增加有助于从图像中过滤电极候选物110-从图像中消除这些电极候选物进一步消除了在医疗设备(诸如导管和传感器上发现的电极)上或附近的假基准的可能性。
在一些特定实施例中,滤波器电极候选物110包括对来自前一步骤109的图像进行阈值处理,以创建二进制图,其中细长结构作为大分量存在。这些大分量对应于细长工具,诸如图像内的导管和其它医疗设备。虽然本实施例实现自适应阈值,但是其它实施例可以实现恒定阈值以过滤电极候选物。
在与本发明一致的多种实施例中,电极增强109可以利用一个或多个尺度,该尺度被用于增强多种大小的球形基准。作为示例,可以使用具有多种尺度的海森滤波器来实现这种电极增强。海森滤波器的两种不同变体可用于增强球形电极的外边缘。
在滤波器电极候选物110的步骤期间,可以产生包含源自背景去除103和电极增强109的分量的二进制图。二进制图可以由“AND”算子或通过融合方法应用。滤波器电极候选步骤110产生输出图像,其中剩余的连通分量限定基准分量。
滤波器电极候选物110的步骤可以基于大小阈值排除小的和大的分量。然后可以基于阈值距离对来自电极增强步骤109的其余分量进行聚类。大的聚类主要包含导管电极,而(稀疏的)基准保持隔离。最后,排除位于前一步骤中形成的在来自大的聚类的阈值距离以下的基准候选。
一旦完成图像处理,包括由图像的三个并行过滤路径120、130和140中的每一个过滤路径进行的并行图像处理,就可以完成基准111的最终选择。由于三个并行过滤路径120、130和140已经去除了与导管和电极相关联的假基准,因此唯一剩余的基准将是由第一滤波器路径120增强的真实基准。因此,选择剩余的基准,并且系统利用剩余的基准来确定将临床图像的第一笛卡尔坐标系与导航系统的第二笛卡尔坐标系相关联的变换模型。
图2A是与本发明的多种实施例一致的患者胸部的临床图像200,在该胸部中具有3个微创诊断/治疗导管。如图2A中所示,探索性导管210A-B通过患者的脉管系统向心肌延伸,并且电生理学环导管205已经延伸通过患者的脉管系统并进入心肌。临床图像的背景包含表示患者身体部分(包括肋骨、椎骨和心肌本身)的多种阴影元素。肋骨、椎骨和心脏的阴影限制了背景图像与多种导管205和210A-B以及导管轴220A-B之间的对比度。
在利用图像处理来识别临床图像中的基准2021-2的情况下,这种对比度的缺乏是有问题的。已经发现使用包括基准的OMRP,其中OMRP静态地定位在手术室内(通常在患者下方,例如在手术台内/下面),便于临床图像的第一笛卡尔坐标系与导航系统的第二笛卡尔坐标系的关联。许多临床成像系统包括能够进行多轴旋转的C臂,以便于患者的多种图像。为协调临床图像的第一笛卡尔坐标系,其取决于C臂取向和图像的平移(以及其它因素,诸如手术台平移/旋转、探测器围绕其轴的旋转和SID)而改变,采用导航系统的第二笛卡尔坐标系,OMRP可以在显示期间放置在临床图像内。OMRP上的基准有助于协调唯一坐标系,因为OMRP的位置(出现在临床图像上)相对于导航系统是已知的。基于变换模型和来自导航系统的跟踪数据,跟踪对象(例如,导管)的图像可以覆盖在临床图像上。
为了促进唯一坐标系的自动协调(例如,临床成像系统的第一笛卡尔坐标系和导航系统的第二笛卡尔坐标系),图像处理系统可以能够在没有检测到假基准的情况下准确地识别临床图像内的基准元素。一旦检测到基准元素,就可以确定变换模型以将由导航系统跟踪的对象的确定位置适当地叠加到临床图像上。如果假基准由图像处理系统检测到并且据此确定变换模型,则坐标系的协调将是不准确的-从而在临床图像上的导航系统信息的超定位中产生错误。
如图2A中所示,临床图像200的背景与多种对象(例如,多种导管205和210A-B以及导管轴220A-B)之间缺乏对比度使得图像处理系统难以将多种对象和基准2021-2与背景区分开-导致临床图像200内的基准识别自动化中的错误率增加。
图2B是与本发明的多种实施例一致的图2A的临床图像200的过滤版本,在该图2A中图像201中的基准2021-2已经增强。虽然导管轴220A-B与假基准(例如,多种导管205和210A-B上的电极2151-N)之间的对比率保持相同,但是改善了基准2021-2与临床图像201的背景对象之间的对比率。可用于促进基准增强的算法的一些示例包括高斯拉普拉斯算子(“LoG”)、高斯差分(“DoG”)、海森滤波器、方向可控滤波器和非线性扩散反应滤波器(例如,基于贝尔特拉米的滤波器)。
图3是其中具有微创诊断/治疗导管的患者胸部的临床图像300,临床图像300用两个滤波器过滤。第一滤波器是非线性滤波器(贝尔特拉米),其增强球形基准,以及在多种导管305和310A-B上的诸如电极315 1-N的假基准。与本发明的多种实施例一致,应用第二滤波器以去除背景图像(例如,图像的可包括诸如肋骨、阴影、骨骼、心肌等的大尺度元素同时排除诸如基准、导管和电极的小尺度元素的部分)。通过从图像中去除背景细节,图像处理系统进一步增加了基准3021-2和电极3151-N以及背景图像(其中几乎全部被去除)之间的对比度。背景去除步骤主要从图像300中去除诸如肋骨、阴影、骨骼、心肌等的大尺度细节。
图4是其中具有微创诊断/治疗导管的患者胸部的临床图像400。采用图3的图像300的输入,图像处理系统进一步采用海森滤波器过滤图像,以增强图像400中的基准4021-2。海森滤波器从图像300中的导管轴320A-B去除剩余的伪像。在将海森滤波器应用于图像300之后,图像400内的剩余对象是基准标记4021-2和假基准(例如,多种导管405和410A-B上的电极4151-N)。
在图5中,使用测地距离进一步处理图4的临床图像400以突出图像500内的聚类的假基准元素,并进一步抑制背景。测地距离利用导管上发现的假基准的紧密接近关系。具体地,图像500内的多种导管505和510A-B上的电极5151-N通常聚集在导管的远侧末端附近,以便于例如在导管的远侧末端处的治疗和/或定位。由于图像400中的真实基准具有大得多的间距,因此假基准、基准标记4021-2不出现在处理图像500中。因此,图像处理系统可以从图像400和图像500的比较推断出真实基准标记4021-2的位置。如上面更详细讨论的,基于真实基准的位置,可以使用变换模型来协调导航系统和临床图像500的坐标系。
图6A是其中具有微创诊断/治疗导管610A-B和导管轴620A-B的患者胸部的临床图像600。图像600进一步包括线束622和位于患者胸部上的患者参考传感器621(“PRS”)。如在图2A的图像200中,图像600的背景包括指示患者身体部分(包括肋骨、椎骨和心肌本身)的多种阴影元素。肋骨、椎骨和心脏的阴影限制了图像背景、基准标记6021-2和导管电极6151-N之间的对比度,以及图像600中的其它特征。图像处理系统还可能易于错误地识别PRS 621和/或线束622内的线的交叉点作为基准标记。因此,图像处理系统可以对图像600实施过滤以去除这种假基准诱导对象。
图6B示出了在自适应阈值处理和海森滤波器已经应用于原始临床图像600(如图6A中所示)之后的临床图像601。自适应阈值处理和海森滤波器的组合使得与原始图像600内的导管610A-B和线束622相关联的暗段呈现为连通分量。只留下PRS 621的一小部分伪像。然后,图像处理系统可以从图像600和601的比较推断出真实基准标记6021-2的位置。
图7A是患者胸部的临床图像700,具有延伸通过患者的心血管系统的微创诊断/治疗导管705和710A-B。在图7A的图像700中,分别由图像处理系统识别真假基准7021-2和7251-N,但是根据本发明未对图像700进行过滤。因此,所显示的图像700已被增强以覆盖与导管705的电极相关联的所有真实基准7021-2和假基准7251-N之上的高光。这种结果将显著地损害用于协调导航系统和临床成像系统的坐标系的变换模型的准确性。
图7B示出了与本发明的多个方面一致的在过滤之后的图7A的临床图像700。由图像处理系统进行的滤波,质心提取,识别与导管705相关联的假基准7251-N,并且增强显示以覆盖仅在真实基准7021-2之上的高光。
虽然图7A-B中所示的显示器旨在通过图像处理系统在视觉上示出基准标记的自动识别,但是这种增强显示器可以在一些实施例中用于进一步验证图像处理系统对基准选择的准确性。当需要新的变换模型来协调导航系统和临床成像系统的坐标系时,图像处理系统可以实现本文公开的多种过滤方案并向临床医生显示类似于图7A-B的增强图像,以便于自动化过程中的人工干预。在这种实施例中,用户将验证图像处理系统对基准的识别的准确性,或者修改选择以校正错误。
在一些特定实施例中,为了促进改进的准确性,图像处理系统可以将每个临床图像分成一个或多个子帧。在这种实施例中,可以基于其中潜在的基准标记的定位将临床图像划分为一个或多个子帧。在特定实施例中,每个子帧将仅具有一个潜在的基准标记或潜在的基准标记的一部分。然后可以使用本文公开的多种图像处理、滤波和增强技术来单独处理每个子帧,以识别真实的基准标记并滤除子帧内的任何错误基准标记。一旦处理了图像的所有子帧,就可以重新编译图像以便显示给例如临床医生。
本发明的一些特定过滤方法包括从临床图像中去除覆盖层。覆盖层是不在原始临床图像中的附加标记,但是由临床成像系统添加。例如,覆盖层可以包括固定和/或变化的标记。本文公开的过滤方法可以在识别基准标记之前从图像中去除这种覆盖层,这进一步降低了识别假基准的可能性。
尽管上面已经以一定程度的特殊性描述了若干实施例,但是本领域技术人员可以在不脱离本发明的精神的情况下对所公开的实施例进行多种改变。旨在将以上描述中包含的或附图中示出的所有内容解释为仅是说明性的而非限制性的。在不脱离本说明的情况下,可以进行细节或结构的改变。前面的描述和所附权利要求旨在涵盖所有这些修改和变化。
基于以上讨论和说明,本领域技术人员将容易认识到,可以对多种实施例进行多种修改和改变,而不严格遵循本文示出和描述的示例性实施例和应用。例如,上面公开的一个或多个成像处理滤波器可以串行和/或并行使用,以提供附加的基准标记结果,以交叉检查利用不同滤波器和/或滤波器组合的另一结果。这些修改不脱离本发明各个方面的真实精神和范围,包括权利要求中阐述的方面。
可以实现多种模块或其它电路以执行本文描述和/或图中所示的一个或多个操作和活动。在这些上下文中,“模块”是执行这些或相关操作/活动中的一个或多个的电路(例如,图像处理系统)。例如,在某些上述实施例中,一个或多个模块是离散逻辑电路或可编程逻辑电路,其被配置和布置用于实现这些操作/活动。在某些实施例中,这种可编程电路是被编程为执行一组(或多组)指令(和/或配置数据)的一个或多个计算机电路。指令(和/或配置数据)可以是存储在存储器(电路)中并可从存储器(电路)访问的固件或软件的形式。作为示例,第一和第二模块包括基于CPU硬件的电路和固件形式的一组指令的组合,其中第一模块包括具有一组指令的第一CPU硬件电路,并且第二模块包括具有另一组指令的第二CPU硬件电路。
某些实施例涉及一种计算机程序产品(例如,非易失性存储器设备),其包括其上存储有指令的机器或计算机可读介质,该指令可由计算机(或其他电子设备)执行,以执行这些操作/活动。
本文描述了多种设备、系统和方法的多种实施例。阐述了许多具体细节以提供对说明书中描述的和附图中示出的实施例的整体结构、功能、制造和使用的透彻理解。然而,本领域技术人员将理解,可以在没有这些具体细节的情况下实践这些实施例。在其它情况下,没有详细描述公知的操作、组件和元件,以免模糊说明书中描述的实施例。本领域普通技术人员将理解,本文描述和示出的实施例是非限制性实施例,并且因此可以理解,本文公开的具体结构和功能细节可以是代表性的,并不一定限制实施例的范围,其范围仅由所附权利要求限定。
在整个说明书中对“多种实施例”、“一些实施例”、“一个实施例”、“实施例”等的引用意味着结合该实施例描述的特定特征、结构或特性包括在至少一个实施例中。因此,在整个说明书中的地方出现的短语“在多种实施例中”、“在一些实施例中”、“在一个实施例中”、“在实施例中”等等不一定都指代相同的实施例。此外,特定特征、结构或特性可以在一个或多个实施例中以任何合适的方式组合。因此,结合一个实施例示出或描述的特定特征、结构或特性可以整体或部分地与一个或多个其它实施例的特征结构或特性组合而没有限制。
应当理解,术语“近端”和“远端”可以在整个说明书中参考临床医生操作用于治疗患者的器械的一端来使用。术语“近端”是指器械中最靠近临床医生的部分,并且术语“远端”是指距离临床医生最远的部分。将进一步理解,为了简明和清楚起见,这里可以相对于所示实施例使用诸如“垂直”、“水平”、“向上”和“向下”的空间术语。然而,手术器械可以在许多取向和位置中使用,并且这些术语不是限制性的和绝对的。
声明通过引用并入本文的全部或部分的任何专利、出版物或其它公开材料仅在所并入的材料不与现有定义、陈述或本发明中阐述的其它公开材料冲突的程度内并入本文。因此,并且在必要的程度上,本文明确阐述的公开内容取代通过引用并入本文的任何冲突材料。声明通过引用并入本文但与现有定义、陈述或本文阐述的其它公开材料相冲突的任何材料或其部分仅在所并入的材料与现有公开材料之间不发生冲突的程度上被并入。
Claims (49)
1.一种用于在临床图像上识别和定位基准标记的方法,所述方法包括:
从成像系统接收所述临床图像;
处理所述临床图像以去除假基准标记;以及
识别所述临床图像中的基准标记。
2.根据权利要求1所述的方法,进一步包括将增强滤波器应用于所述临床图像以增强基准标记的外观。
3.根据权利要求2所述的方法,其中,所述增强滤波器能够是以下中的一个或多个:高斯拉普拉斯算子、高斯差分、海森滤波器、方向可控滤波器和非线性扩散反应滤波器。
4.根据权利要求2所述的方法,其中,所述增强滤波器是基于贝尔特拉米的滤波器。
5.根据权利要求2所述的方法,其中,所述增强滤波器是具有基于所述图像的二阶导数算子-ΔI=Ixx+Ixy的高斯拉普拉斯滤波器。
6.根据权利要求5所述的方法,其中,所述增强滤波器包括所述高斯拉普拉斯滤波器和高斯差分滤波器。
7.根据权利要求2所述的方法,其中,所述增强滤波器是海森滤波器,所述海森滤波器包括:
提取所述临床图像内的每个像素的海森矩阵,每个像素的海森矩阵是:
其中是图像的平滑二阶导数;
应用高斯西格玛参数来平滑所述图像;
对所述海森矩阵执行特征值多尺度分析;
排序特征值和
通过对象尺度对局部结构进行排序,并通过分析所述海森矩阵的所述特征值来进行形状辨别;
用以下方式增强对象的多尺度响应:
以及
用以下方式增强线结构响应:
8.根据权利要求2所述的方法,其中,所述增强滤波器是包括非线性扩散反应滤波器的基于贝尔特拉米的滤波器,所述非线性扩散反应滤波器被配置为在保留边界线的同时过滤图像噪声。
9.根据权利要求8所述的方法,其中,所述非线性扩散反应滤波器包括扩散项和反应项,以促进所述图像内的基准标记的对比度增强。
10.根据权利要求2所述的方法,其中,在所述增强滤波器之后,基于进化模型,采用偏微分方程的渐近状态来解读所述临床图像。
11.根据权利要求10所述的方法,其中,所述进化模型是:
It=ΔgI+f(I)
其中Δg是拉普拉斯-贝尔特拉米算子,并且所述函数f是反应项。
12.一种检测临床图像内的医疗设备电极的方法,所述方法包括:
从成像系统接收所述临床图像;
将增强滤波器应用于所述临床图像以增强电极的外观;
去除所述临床图像的背景;以及
识别所述临床图像内的电极。
13.根据权利要求12所述的方法,其中,识别所述临床图像内的电极包括自适应阈值处理、基于大小和形状的分量滤波,以及电极的提取。
14.根据权利要求12所述的方法,其中,识别所述临床图像内的电极包括自适应阈值处理。
15.根据权利要求12所述的方法,进一步包括过滤所述临床图像以去除假基准。
16.根据权利要求15所述的方法,其中,过滤所述临床图像以去除假基准的步骤包括使用以下中的一个或多个:基于贝尔特拉米的增强滤波器、具有基于所述图像的二阶导数的算子的高斯拉普拉斯增强滤波器、以及高斯拉普拉斯滤波器和高斯差分滤波器。
17.根据权利要求15所述的方法,其中,过滤所述临床图像以去除假基准的步骤包括用海森滤波器过滤所述临床图像,用海森滤波器过滤包括:
提取所述临床图像内的每个像素的海森矩阵,用于每个像素的海森矩阵为:
其中是图像的平滑二阶导数;
应用高斯西格玛参数来平滑所述图像;
对所述海森矩阵执行特征值多尺度分析;
排序特征值和
通过对象尺度对局部结构进行排序,并通过分析所述海森矩阵的所述特征值来进行形状辨别;
用以下方式增强对象的多尺度响应:
以及
用以下方式增强线结构响应:
18.一种从基准标记中检测和区分导管和导管部件的方法,所述方法包括:
从成像系统接收临床图像;
处理所述临床图像以从临床图像中去除导管和导管部件;以及
识别所述临床图像中的基准标记。
19.根据权利要求18所述的方法,进一步包括将增强滤波器应用于所述临床图像以增强基准标记的外观。
20.根据权利要求19所述的方法,其中,所述增强滤波器包括以下中的一个或多个:高斯拉普拉斯算子、高斯差分、海森滤波器、方向可控滤波器和非线性扩散反应滤波器。
21.一种用于将医疗设备的位置数据从导航系统叠加到临床图像上的方法,所述方法包括:
从临床成像系统接收所述临床图像;
处理所述临床图像以去除假基准标记;
识别所述临床图像内的基准标记;
创建变换模型,所述变换模型将协调所述导航系统的第一坐标系与所述临床成像系统的第二坐标系;
确定所述第一坐标系中医疗设备的位置;
将所述变换模型应用于所述第一坐标系中医疗设备的位置,以确定所述第二坐标系中医疗设备的位置;以及
基于所述第二坐标系中医疗设备的已知位置,将所述医疗设备的图像叠加到临床图像上。
22.根据权利要求21所述的方法,其中处理所述临床图像以去除假基准标记的步骤包括将增强滤波器应用于临床图像以增强基准标记的外观。
23.根据权利要求22所述的方法,其中,所述增强滤波器能够是以下中的一个或多个:高斯拉普拉斯算子、高斯差分、海森滤波器、方向可控滤波器和非线性扩散反应滤波器。
24.根据权利要求22所述的方法,其中,所述增强滤波器是基于贝尔特拉米的滤波器。
25.根据权利要求22所述的方法,其中,所述增强滤波器是具有基于图像的二阶导数算子-ΔI=Ixx+Ixy的高斯拉普拉斯滤波器。
26.根据权利要求25所述的方法,其中,所述增强滤波器包括所述高斯拉普拉斯滤波器和高斯差分滤波器。
27.根据权利要求22所述的方法,其中,所述增强滤波器是海森滤波器,所述海森滤波器包括:
提取所述临床图像内的每个像素的海森矩阵,用于每个像素的海森矩阵是:
其中是所述图像的平滑二阶导数;
应用高斯西格玛参数来平滑所述图像;
对所述海森矩阵执行特征值多尺度分析;
排序特征值和
通过对象尺度对局部结构进行排序,并通过分析所述海森矩阵的所述特征值来进行形状辨别;
用以下方式增强对象的多尺度响应:
以及
用以下方式增强线结构响应:
28.根据权利要求22所述的方法,其中,所述增强滤波器是包括非线性扩散反应滤波器的基于贝尔特拉米的滤波器,所述非线性扩散反应滤波器被配置为在保留边界线的同时过滤图像噪声。
29.根据权利要求28所述的方法,其中,所述扩散反应滤波器包括扩散项和反应项,以促进所述图像内的基准标记的对比度增强。
30.根据权利要求22所述的方法,其中,在所述增强滤波器之后,基于进化模型,采用偏微分方程的渐近状态来解读所述临床图像。
31.根据权利要求30所述的方法,其中,所述进化模型是:
It=ΔgI+f(I)
其中Δg是拉普拉斯-贝尔特拉米算子,并且所述函数f是反应项。
32.一种用于识别临床图像内的基准标记的方法,所述方法包括:
处理所述临床图像以去除假基准标记;以及
识别所述临床图像内的基准标记。
33.根据权利要求32所述的方法,进一步包括增强所述临床图像内的基准标记。
34.根据权利要求32所述的方法,进一步包括去除所述临床图像内的背景。
35.根据权利要求32所述的方法,其中,处理所述临床图像以去除假基准标记的步骤包括增强和过滤临床图像内的多个导管。
36.根据权利要求32所述的方法,其中,处理所述临床图像以去除假基准标记的步骤包括增强和过滤临床图像内的多个电极。
37.根据权利要求32所述的方法,其中,识别所述临床图像内的基准标记包括自适应阈值处理、基于大小和形状的分量滤波、以及基准标记的提取。
38.根据权利要求32所述的方法,其中,识别所述临床图像内的基准标记包括自适应阈值处理。
39.根据权利要求32所述的方法,其中,处理所述临床图像以去除假基准标记的步骤包括使用基于贝尔特拉米的增强滤波器过滤临床图像。
40.根据权利要求32所述的方法,其中,处理所述临床图像以去除假基准标记的步骤包括采用具有基于图像的二阶导数的算子的高斯拉普拉斯增强滤波器过滤所述临床图像。
41.根据权利要求32所述的方法,其中,处理所述临床图像以去除假基准标记的步骤包括用高斯拉普拉斯滤波器和高斯差分滤波器过滤所述临床图像。
42.根据权利要求32所述的方法,其中,处理所述临床图像以去除假基准标记的步骤包括用海森滤波器过滤临床图像,用所述海森滤波器过滤包括:
提取所述临床图像内的每个像素的海森矩阵,用于每个像素的海森矩阵是:
其中是所述图像的平滑二阶导数;
应用高斯西格玛参数来平滑所述图像;
对所述海森矩阵执行特征值多尺度分析;
排序特征值和
通过对象尺度对局部结构进行排序,并通过分析所述海森矩阵的所述特征值来进行形状辨别;
用以下方式增强对象的多尺度响应:
以及
用以下方式增强线结构响应:
43.根据权利要求39所述的方法,其中,所述基于贝尔特拉米的滤波器包括非线性扩散反应滤波器,以在保留边界线的同时过滤图像噪声。
44.一种用于心血管导管的导航系统,所述导航系统包括:
心血管导管,其包括位于所述导管的远侧末端附近的一个或多个电极,所述电极被配置和布置成便于导管的远侧末端的定位;
光磁配准板,其包括多个基准标记;
临床成像系统,其被配置和布置成在第一坐标系中显示包括患者、基准标记和心血管导管的临床图像;
导管定位系统,其被配置和布置成在第二坐标系中检测所述一个或多个电极的位置和方向;以及
控制器电路,其通信地耦合到所述临床成像系统和导管定位系统,所述控制器电路被配置和布置成
基于所述第二坐标系中的基准标记的已知位置和所述临床图像内基准标记的位置,确定第一和第二坐标系之间的变换模型,以及
基于所述变换模型,确定所述第一坐标系内电极的位置。
45.根据权利要求44所述的导航系统,进一步包括显示器;所述控制器电路通信地耦合到所述显示器,并且进一步被配置和布置成基于所确定的所述第一坐标系中的电极的位置,生成用于所述显示器的信号,所述信号包括叠加有导管的代表性图像的临床图像。
46.根据权利要求44所述的导航系统,其中,所述控制器电路进一步被配置和布置成
过滤所述临床图像以去除假基准标记,以及
识别所述临床图像内的基准标记。
47.根据权利要求46所述的导航系统,其中,所述控制器电路进一步被配置和布置成将增强滤波器应用于所述临床图像,以增强基准标记和假基准标记相对于背景的外观。
48.根据权利要求47所述的导航系统,其中,所述增强滤波器包括以下中的一个或多个:高斯拉普拉斯算子、高斯差分、海森滤波器、方向可控滤波器和非线性扩散反应滤波器。
49.根据权利要求1所述的方法,其中,识别所述临床图像中的基准标记的步骤包括以非均匀图案识别基准标记。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201762469423P | 2017-03-09 | 2017-03-09 | |
US62/469423 | 2017-03-09 | ||
PCT/IB2018/051527 WO2018163105A2 (en) | 2017-03-09 | 2018-03-08 | Detection of fiducials in a clinical image |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110352447A true CN110352447A (zh) | 2019-10-18 |
Family
ID=62067744
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201880014983.0A Pending CN110352447A (zh) | 2017-03-09 | 2018-03-08 | 临床图像中的基准检测 |
Country Status (5)
Country | Link |
---|---|
US (1) | US20210142504A1 (zh) |
EP (1) | EP3563340A2 (zh) |
JP (1) | JP2020509826A (zh) |
CN (1) | CN110352447A (zh) |
WO (1) | WO2018163105A2 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111539354A (zh) * | 2020-04-27 | 2020-08-14 | 易普森智慧健康科技(深圳)有限公司 | 一种液基细胞玻片扫描区域识别方法 |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11334998B2 (en) * | 2019-04-22 | 2022-05-17 | Canon Medical Systems Corporation | Medical image processing apparatus, X-ray diagnostic apparatus, and computer-implemented method |
WO2021061960A1 (en) * | 2019-09-24 | 2021-04-01 | Nuvasive, Inc. | Systems and methods for three-dimensional navigation of objects |
US20220370146A1 (en) * | 2021-05-19 | 2022-11-24 | Globus Medical, Inc. | Intraoperative alignment assessment system and method |
US11756178B1 (en) * | 2022-09-30 | 2023-09-12 | Illuscio, Inc. | Systems and methods for employing invisible fiducials for aligning scans or images to form composite images, models, or datasets |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070286336A1 (en) * | 2006-06-07 | 2007-12-13 | Sylvain Bernard | Method for the processing of radiological images for a detection of radiological signs |
US20090324087A1 (en) * | 2008-06-27 | 2009-12-31 | Palo Alto Research Center Incorporated | System and method for finding stable keypoints in a picture image using localized scale space properties |
CN102203826A (zh) * | 2008-12-25 | 2011-09-28 | 梅迪奇视觉成像解决方案有限公司 | 医学图像的降噪 |
US20120301053A1 (en) * | 2011-05-26 | 2012-11-29 | The University Of Tokyo | Image processing method and image processing apparatus |
US20130272592A1 (en) * | 2010-12-30 | 2013-10-17 | Uzi Eichler | System and Method for Registration of Fluoroscopic Images in a Coordinate System of a Medical System |
WO2014155203A1 (en) * | 2013-03-28 | 2014-10-02 | Innomed Srl | Method of localization of electrodes for monitoring bioelectric signals from the patient's brain |
US20140376796A1 (en) * | 2012-03-13 | 2014-12-25 | Fujifilm Corporation | Image Processing Apparatus, Method, and Program |
US20150087965A1 (en) * | 2013-09-20 | 2015-03-26 | Junichi Tokuda | System and method for automatic detection and registration of medical images |
CN104519793A (zh) * | 2012-05-08 | 2015-04-15 | 安吉奥梅特里克斯公司 | 用于内腔的线性映射的系统 |
CN105188550A (zh) * | 2013-03-12 | 2015-12-23 | 光学实验室成像公司 | 血管数据处理和图像配准系统、方法及装置 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5799099A (en) * | 1993-02-12 | 1998-08-25 | George S. Allen | Automatic technique for localizing externally attached fiducial markers in volume images of the head |
US5391199A (en) | 1993-07-20 | 1995-02-21 | Biosense, Inc. | Apparatus and method for treating cardiac arrhythmias |
US5558091A (en) | 1993-10-06 | 1996-09-24 | Biosense, Inc. | Magnetic determination of position and orientation |
US6690963B2 (en) | 1995-01-24 | 2004-02-10 | Biosense, Inc. | System for determining the location and orientation of an invasive medical instrument |
DE69835422T2 (de) | 1998-01-22 | 2006-12-21 | Biosense Webster, Inc., Diamond Bar | Messung im körperinneren |
US6233476B1 (en) | 1999-05-18 | 2001-05-15 | Mediguide Ltd. | Medical positioning system |
US7386339B2 (en) | 1999-05-18 | 2008-06-10 | Mediguide Ltd. | Medical imaging and navigation system |
US20030228044A1 (en) * | 2002-06-05 | 2003-12-11 | Canon Kabushiki Kaisha | Radiographic marker location |
US7197354B2 (en) | 2004-06-21 | 2007-03-27 | Mediguide Ltd. | System for determining the position and orientation of a catheter |
JP2008009549A (ja) * | 2006-06-27 | 2008-01-17 | Yamaguchi Univ | 画像処理方法、画像処理装置及び画像処理プログラム |
IL188262A (en) | 2007-01-10 | 2011-10-31 | Mediguide Ltd | System and method for superimposing a representation of the tip of a catheter on an image acquired by a moving imager |
-
2018
- 2018-03-08 JP JP2019548474A patent/JP2020509826A/ja active Pending
- 2018-03-08 EP EP18720771.7A patent/EP3563340A2/en not_active Withdrawn
- 2018-03-08 US US16/491,564 patent/US20210142504A1/en not_active Abandoned
- 2018-03-08 WO PCT/IB2018/051527 patent/WO2018163105A2/en unknown
- 2018-03-08 CN CN201880014983.0A patent/CN110352447A/zh active Pending
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070286336A1 (en) * | 2006-06-07 | 2007-12-13 | Sylvain Bernard | Method for the processing of radiological images for a detection of radiological signs |
US20090324087A1 (en) * | 2008-06-27 | 2009-12-31 | Palo Alto Research Center Incorporated | System and method for finding stable keypoints in a picture image using localized scale space properties |
CN102203826A (zh) * | 2008-12-25 | 2011-09-28 | 梅迪奇视觉成像解决方案有限公司 | 医学图像的降噪 |
US20130272592A1 (en) * | 2010-12-30 | 2013-10-17 | Uzi Eichler | System and Method for Registration of Fluoroscopic Images in a Coordinate System of a Medical System |
US20120301053A1 (en) * | 2011-05-26 | 2012-11-29 | The University Of Tokyo | Image processing method and image processing apparatus |
US20140376796A1 (en) * | 2012-03-13 | 2014-12-25 | Fujifilm Corporation | Image Processing Apparatus, Method, and Program |
CN104519793A (zh) * | 2012-05-08 | 2015-04-15 | 安吉奥梅特里克斯公司 | 用于内腔的线性映射的系统 |
CN105188550A (zh) * | 2013-03-12 | 2015-12-23 | 光学实验室成像公司 | 血管数据处理和图像配准系统、方法及装置 |
WO2014155203A1 (en) * | 2013-03-28 | 2014-10-02 | Innomed Srl | Method of localization of electrodes for monitoring bioelectric signals from the patient's brain |
US20150087965A1 (en) * | 2013-09-20 | 2015-03-26 | Junichi Tokuda | System and method for automatic detection and registration of medical images |
Non-Patent Citations (1)
Title |
---|
R. MALLADI: "Fast Difference Schemes for Edge Enhancing Beltrami Flow", 《COMPUTER VISION》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111539354A (zh) * | 2020-04-27 | 2020-08-14 | 易普森智慧健康科技(深圳)有限公司 | 一种液基细胞玻片扫描区域识别方法 |
Also Published As
Publication number | Publication date |
---|---|
US20210142504A1 (en) | 2021-05-13 |
WO2018163105A3 (en) | 2018-10-25 |
WO2018163105A2 (en) | 2018-09-13 |
EP3563340A2 (en) | 2019-11-06 |
JP2020509826A (ja) | 2020-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP7093801B2 (ja) | 手術中の位置調整および誘導を容易にするシステム | |
JP7190145B2 (ja) | 患者に貼付される光学コードを用いた患者の画像データと患者の実景との位置合わせ | |
CN110352447A (zh) | 临床图像中的基准检测 | |
US10546396B2 (en) | System and method for registration of fluoroscopic images in a coordinate system of a medical system | |
CN102713976B (zh) | 对介入装置进行导航 | |
JP2022507622A (ja) | 拡張現実ディスプレイでの光学コードの使用 | |
EP2680755B1 (en) | Visualization for navigation guidance | |
US20090012382A1 (en) | Method and system for detection of obstructions in vasculature | |
CN107809955B (zh) | 经由感兴趣界标的自动检测在x射线成像中进行实时准直和roi过滤器定位 | |
EP3224804B1 (en) | Apparatus for determining positions of an interventional instrument in a projection image | |
JP5676269B2 (ja) | 脳画像データの画像解析 | |
WO2012117321A1 (en) | Visualization for navigation guidance | |
CN101563708A (zh) | 用于指示医学成像数据中可能经计算机检测的假阳性的装置和方法 | |
JP7463625B2 (ja) | ナビゲーションサポート | |
EP2572333A1 (en) | Handling a specimen image | |
US11816821B2 (en) | Method and system for generating an enriched image of a target object and corresponding computer program and computer-readable storage medium | |
CN108135552A (zh) | 投影x射线图像的改进的可视化 | |
JP7303144B2 (ja) | 放射線撮影装置、放射線画像の対象物検出プログラムおよび放射線画像における対象物検出方法 | |
US20220230328A1 (en) | Information processing apparatus, information processing method, and non-transitory storage medium | |
WO2022006586A1 (en) | Extended-reality visualization of endovascular navigation |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20191018 |
|
WD01 | Invention patent application deemed withdrawn after publication |