CN108697402B - 在三维图像中确定深度脑刺激电极的旋转方位 - Google Patents
在三维图像中确定深度脑刺激电极的旋转方位 Download PDFInfo
- Publication number
- CN108697402B CN108697402B CN201780003514.4A CN201780003514A CN108697402B CN 108697402 B CN108697402 B CN 108697402B CN 201780003514 A CN201780003514 A CN 201780003514A CN 108697402 B CN108697402 B CN 108697402B
- Authority
- CN
- China
- Prior art keywords
- data
- image
- image data
- computer
- orientation
- 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
- 230000000638 stimulation Effects 0.000 title description 12
- 210000004556 brain Anatomy 0.000 title description 7
- 238000000034 method Methods 0.000 claims abstract description 93
- 210000003484 anatomy Anatomy 0.000 claims abstract description 72
- 238000003384 imaging method Methods 0.000 claims abstract description 62
- 238000002059 diagnostic imaging Methods 0.000 claims abstract description 23
- 238000003672 processing method Methods 0.000 claims abstract description 5
- 230000009466 transformation Effects 0.000 claims description 42
- 238000004422 calculation algorithm Methods 0.000 claims description 31
- 230000004927 fusion Effects 0.000 claims description 24
- 238000013500 data storage Methods 0.000 claims description 20
- 238000003325 tomography Methods 0.000 claims description 14
- 238000004590 computer program Methods 0.000 claims description 11
- 238000002595 magnetic resonance imaging Methods 0.000 claims description 10
- 238000002601 radiography Methods 0.000 claims description 5
- 238000003708 edge detection Methods 0.000 claims description 2
- 238000004070 electrodeposition Methods 0.000 claims description 2
- 239000003550 marker Substances 0.000 description 30
- 238000002591 computed tomography Methods 0.000 description 17
- 238000012545 processing Methods 0.000 description 15
- 238000005457 optimization Methods 0.000 description 13
- 206010028980 Neoplasm Diseases 0.000 description 8
- 239000000463 material Substances 0.000 description 7
- 239000013598 vector Substances 0.000 description 7
- BASFCYQUMIYNBI-UHFFFAOYSA-N platinum Chemical compound [Pt] BASFCYQUMIYNBI-UHFFFAOYSA-N 0.000 description 6
- 238000011524 similarity measure Methods 0.000 description 5
- 210000001519 tissue Anatomy 0.000 description 5
- 238000007408 cone-beam computed tomography Methods 0.000 description 4
- 230000005684 electric field Effects 0.000 description 4
- 239000011521 glass Substances 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 238000002600 positron emission tomography Methods 0.000 description 4
- 239000004065 semiconductor Substances 0.000 description 4
- 208000003174 Brain Neoplasms Diseases 0.000 description 3
- 241001465754 Metazoa Species 0.000 description 3
- 230000003190 augmentative effect Effects 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000001225 therapeutic effect Effects 0.000 description 3
- 238000002604 ultrasonography Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000002583 angiography Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 239000002872 contrast media Substances 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000002598 diffusion tensor imaging Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000000968 medical method and process Methods 0.000 description 2
- 231100000915 pathological change Toxicity 0.000 description 2
- 230000036285 pathological change Effects 0.000 description 2
- 230000001575 pathological effect Effects 0.000 description 2
- 229910052697 platinum Inorganic materials 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 230000001131 transforming effect Effects 0.000 description 2
- 238000012285 ultrasound imaging Methods 0.000 description 2
- JBRZTFJDHDCESZ-UHFFFAOYSA-N AsGa Chemical compound [As]#[Ga] JBRZTFJDHDCESZ-UHFFFAOYSA-N 0.000 description 1
- 229910001218 Gallium arsenide Inorganic materials 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000004931 aggregating effect Effects 0.000 description 1
- 229910045601 alloy Inorganic materials 0.000 description 1
- 239000000956 alloy Substances 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 210000001638 cerebellum Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 210000002451 diencephalon Anatomy 0.000 description 1
- 238000002091 elastography Methods 0.000 description 1
- 238000001839 endoscopy Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 210000000527 greater trochanter Anatomy 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 210000000528 lesser trochanter Anatomy 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 210000003141 lower extremity Anatomy 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 210000001259 mesencephalon Anatomy 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 230000004118 muscle contraction Effects 0.000 description 1
- 210000004126 nerve fiber Anatomy 0.000 description 1
- 239000012811 non-conductive material Substances 0.000 description 1
- 238000009206 nuclear medicine Methods 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 210000002975 pon Anatomy 0.000 description 1
- 230000002980 postoperative effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
- 238000002922 simulated annealing Methods 0.000 description 1
- 238000002603 single-photon emission computed tomography Methods 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 230000004936 stimulating effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 210000001587 telencephalon Anatomy 0.000 description 1
- 238000002560 therapeutic procedure Methods 0.000 description 1
- 238000001931 thermography Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 210000000689 upper leg Anatomy 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/12—Arrangements for detecting or locating foreign bodies
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/44—Constructional features of apparatus for radiation diagnosis
- A61B6/4417—Constructional features of apparatus for radiation diagnosis related to combined acquisition of different diagnostic modalities
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
- A61B6/5247—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from an ionising-radiation diagnostic technique and a non-ionising radiation diagnostic technique, e.g. X-ray and ultrasound
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N1/00—Electrotherapy; Circuits therefor
- A61N1/02—Details
- A61N1/04—Electrodes
- A61N1/05—Electrodes for implantation or insertion into the body, e.g. heart electrode
- A61N1/0526—Head electrodes
- A61N1/0529—Electrodes for brain stimulation
- A61N1/0534—Electrodes for deep brain stimulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/168—Segmentation; Edge detection involving transform domain methods
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
- G16H20/30—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to physical therapies or activities, e.g. physiotherapy, acupressure or exercising
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H40/00—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices
- G16H40/60—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices
- G16H40/63—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices for local operation
-
- 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/08—Accessories or related features not otherwise provided for
- A61B2090/0807—Indication means
-
- 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
- A61B2090/364—Correlation of different images or relation of image positions in respect to the body
-
- 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]
-
- 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/3954—Markers, e.g. radio-opaque or breast lesions markers magnetic, e.g. NMR or MRI
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/4808—Multimodal MR, e.g. MR combined with positron emission tomography [PET], MR combined with ultrasound or MR combined with computed tomography [CT]
- G01R33/4812—MR combined with X-ray or computed tomography [CT]
-
- 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/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- 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/10—Image acquisition modality
- G06T2207/10116—X-ray image
-
- 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/10—Image acquisition modality
- G06T2207/10116—X-ray image
- G06T2207/10124—Digitally reconstructed radiograph [DRR]
-
- 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
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Veterinary Medicine (AREA)
- Animal Behavior & Ethology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Optics & Photonics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- General Physics & Mathematics (AREA)
- Neurology (AREA)
- Neurosurgery (AREA)
- Psychology (AREA)
- Pulmonology (AREA)
- Physical Education & Sports Medicine (AREA)
- General Business, Economics & Management (AREA)
- Business, Economics & Management (AREA)
- Cardiology (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Electrotherapy Devices (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Nuclear Medicine (AREA)
Abstract
公开了一种计算机实现的用于确定电极的方位的医学数据处理方法,电极被配置为电刺激患者的解剖结构并且包括旋转方位标记,所述方法包括在至少一个计算机的至少一个处理器上执行步骤:a)在所述至少一个处理器处获取(S1.1)描述解剖结构和电极的二维医学图像的旋转图像数据,二维医学图像是在医学成像装置相对于解剖结构旋转期间用二维医学成像装置拍摄的,旋转图像数据针对每个二维医学图像进一步描述相对于与相应二维医学图像相关联的解剖结构的成像视角;b)由所述至少一个处理器并且基于旋转图像数据来确定(S1.2)描述二维医学图像的参考系中的电极的旋转方位的旋转方位数据;c)在所述至少一个处理器处获取(S1.3)描述解剖结构的一组断层医学图像的断层图像数据;d)由所述至少一个处理器并且基于旋转图像数据和断层图像数据以及旋转方位数据,来确定(S1.4)描述断层医学图像数据的参考系中的电极的旋转方位的电极方位数据。
Description
技术领域
本发明涉及一种计算机实现的用于确定电极方位的医学方法。电极例如被配置为电(例如,电磁地)刺激患者的解剖结构并且该电极包括旋转方位标记。本发明还涉及一种被配置为执行对应于该方法的程序的计算机和用于确定电极的方位的、包括这种计算机的医学系统。
背景技术
定向深度脑刺激(DBS)引导体(也称为电极)允许将刺激集中到引导体(即电极)的特定角度。为了规划刺激或评估刺激效果,需要知道引导体相对于患者解剖结构的正确的定向方位。引导体配备有旋转非对称标记,该标记设计为通过术后成像来被检测和评估。很遗憾,引导体中的金属导致MRI不合格,并且CT太模糊以至于无法足够准确地识别标记。
所有之前的方法都不允许对限定电极围绕其纵向轴线的旋转角方位的旋转角度进行自动化、可追踪和准确的测量。
本发明被设计为提供一种用于确定三维医学图像的参考系中的电极的方位的可靠方法。
以下公开了本发明的各方面、示例和示例性步骤及其实施方式。只要技术上适宜且可行,本发明的不同示例性特征可根据本发明进行组合。
本发明的示例性简述
在下文中,将给出本发明的具体特征的简短描述,其不应被理解为将本发明仅限于本部分中描述的特征或特征的组合。
所公开的方法使用旋转x射线(例如,来自旋转血管造影术,锥形束CT,C形臂x射线)。在检测器围绕头部旋转时拍摄的一系列x射线照片中,引导体(电极)及其标记在不同角度的图像上清晰可见。这些图像允许确定标记相对于图像平面和/或图像的参考系的角度。为了将这一信息传送至3D患者系统中,借助于自动算法将2D x射线图像配准到解剖三维图像体(volume)中。通过将3D图像数据的数字重建射线照片融合到2D x射线图像上,用户可以立即验证2D和3D图像的匹配并评估计算出的方位角度的准确性。而且,能够确定在3D图像体的参考系中的引导体的旋转方位。
本发明概述
在本节中,例如通过参考本发明的可能实施方式来给出本发明的一般特征的描述。
总体上,本发明通过在第一方面提供一种计算机实现的用于确定电极方位的医学方法来实现上述目标。电极例如被配置为电(例如,电磁地)刺激患者的解剖结构,并且该电极包括旋转方位标记。
电极可以是用于深度脑刺激的定向电极(DBS电极),其包括沿着电极的纵向轴线位于离散段中的不同的且明显的触点,并且电极可以采取不完全围绕电极的周边延伸而是仅沿着围绕电极的纵向轴线的角扇形(环形段)延伸的切片形式。电极通常采用细长形,基本上圆柱形,并且在一个示例中为刚性形状,其具有指向患者身体内部的远端和指向患者身体外部的近端。由于电极的延伸,其被限定为具有例如从其近端延伸到远端的纵向轴线。电极通常是定向DBS电极,其允许仅在相应接触的方向上选择性地刺激三维角扇区。为了实现对作为待刺激的目标区域或包括待刺激的目标区域的解剖结构的期望刺激,必须确定定向电极相对于解剖结构(目标区域)的位置的旋转角方位,以使得解剖结构(目标区域)位于由电极的至少一个定向触点的定向特性(发射特性)覆盖的扇区中。
根据第一方面的方法涉及确定电极的角方位,例如相对于其在二维医学图像中描绘的平面的角方位,以便基于如下所述的进一步的信息来确定例如在描述解剖结构和电极的三维图像数据集的参考系中电极相对于解剖结构的旋转方位。角方位可被定义为在围绕由电极的纵向轴线的延伸限定的旋转轴线的医学二维图像中的电极(例如,方位标记的)的旋转相位,例如被定义为通过二维图像的图像平面中的电极的所描绘的平面部分相对于电极圆周上的方位标记的周向位置的旋转相位。换句话说,角方位被定义为相对于例如圆周上的方位标记的位置沿着电极的圆周的旋转角。甚至换句话说,角方位由一方面方位标记(例如,其第一部分)的预定(例如,已知的或固定的中的至少一种)位置或与该位置相交的电极的半径中的至少一个所处的、以及电极的纵向轴线所在的平面,与另一方面二维医学图像的图像平面之间的角度(例如,可以在如球坐标的三维参考系中定义的空间角度)来定义。这两个平面的方位通常彼此平行和/或两者的方位位于一个共同的平面中。为了实现根据第一方面的方法的效果,可以忽略这两个方位之间的这种空间关系的较小的三维偏差。例如从包括描述(例如,定义或表示)电极的几何结构的结构化数据的电极模板数据中预先确定(例如,已知的或固定的中的至少一种)方位标记的预定位置和电极的其余部分之间的空间关系。其中定义角方位的参考系与其中定义二维医学图像中的电极的剩余部分或图像位置中的至少一个的坐标系之间的坐标变换是预定的(例如,已知的或固定的中的至少一种)。在本公开的含义内,术语参考系涵盖并且通常等同于术语坐标系。因此,可以通过分析二维医学图像来确定二维医学图像(其中定义了定义二维医学图像的图像内容的位置)的参考系中的电极的旋转方位。
注意,将这种电极放置在人体或动物体内并不属于所公开的方法的一部分。此外,执行所公开的方法不需要在执行所公开的方法的同时放置电极。换句话说,根据第一方面的所公开的方法是操作或控制如计算机的数据处理设备的方法中的至少一种,例如通过在该数据处理设备上运行相应的程序以使该数据处理设备执行根据第一方面的方法的步骤。
根据第一方面的方法是例如数据处理方法。根据第一方面的方法包括在至少一台计算机(例如,至少一台计算机是深度脑刺激校准或调谐系统的一部分)的至少一个处理器上执行由该至少一个处理器执行的以下示例性步骤。
在(例如第一)示例性步骤中,获取旋转图像数据,其描述(例如,定义或表示中的至少一种)解剖结构和电极的二维医学图像(即,描绘或示出)。具体而言,二维图像(例如每个)既描绘解剖结构又描绘电极。通常,解剖结构可以是包含至少一个神经纤维的任何解剖学身体部位,该至少一个神经纤维可以被电(电磁地)刺激(例如,解剖结构包括,例如是至少部分脑),或者该解剖结构可以是可以由电刺激而引起例如肌肉收缩或伸展的肌肉组织。
通过旋转血管造影术或锥形束计算机x射线断层扫描成像(锥形束CT),已经或者正在使用二维医学成像装置或例如利用C形臂(C形臂x射线,也称为C形臂放射线摄影术)成像的方法来采集二维医学图像。例如,在医学成像装置相对于解剖结构旋转期间已经或正在采集二维图像,即,在医学成像装置相对于解剖结构和电极的位置的不同旋转位置,已经拍摄了每一个二维图像,例如这使得医学成像装置的紧接着的后续位置与拍摄的紧接其后的二维图像按其生成顺序相关联。例如,每个二维医学图像与不同的成像视角(例如,相对于解剖结构的位置)相关联。旋转图像数据针对每个二维医学图像进一步描述相对于与相应二维医学图像相关联的解剖结构的成像视角。成像视角优选地由医学成像装置的位置来定义,该位置为当相应的二维医学图像生成时,医学成像装置相对于解剖结构以及例如相对于电极的位置。生成旋转图像数据(即,进行相关成像过程)的步骤通常不是所要求保护的方法的一部分,并且在一个示例中,该步骤在根据第一方面的方法被执行之前进行。然后通过该方法在稍后的时间点获取(读取)旋转图像数据。然而,在该方法的另一个示例中,旋转图像数据的生成可以是根据第一方面的方法的一部分。
在又一(例如第二)示例性步骤中,确定旋转方位数据,其描述(例如定义或表示)二维医学图像的参考系(其中定义了定义旋转图像数据的图像内容的位置)中的电极的旋转方位。二维图像的参考系是例如用于生成旋转图像数据的参考系,并且可以由二维医学成像装置(如C形臂x射线设备)预先确定。例如基于如来自相应二维图像中的电极的图像描绘的旋转图像数据来确定旋转方位数据。为此,分析关于二维图像中的方位标记的图像外观的至少一个二维医学图像。在一个示例中,相应地分析所有的二维医学图像。这可以例如通过基于旋转图像数据而确定方位标记的图像外观来完成,例如通过以下至少一项:
-分割至少一个或每一个二维医学图像中的电极的图像外观;
-对至少一个或每一个二维医学图像的组成部分进行边缘检测;
-对至少一个或每一个二维医学图像中的电极的图像外观与描述电极的结构化数据(如其几何形状和至少一个定向触点和方位标记之间的空间关系中的至少一个,其中空间关系为位置和方向中的至少一个)的先前获取并预定的电极模板数据进行比较。
然后从关于与二维医学图像的参考系中的相应二维图像相关联的成像视角的信息以及例如先前获取的电极模板数据来确定电极的旋转方位,如必要,将电极的部分的位置或方位或由电极模板数据描述的空间关系的任何重新计算考虑到二维图像的参考系中。
在进一步的示例中,例如第二示例性步骤包括,获取描述(例如,定义或表示)二维医学图像的参考系中的电极的预定旋转方位的预定方位数据。方位标记相对于电极的位置和方向例如是预定的(例如,已知的或固定的中的至少一种),并且包括在电极模板数据中或者可以从电极模板数据导出,为此,电极模板数据可以通过根据第一个方面的方法获取。预定方位可以是例如方位标记在二维图像中具有预定(例如,已知的或固定的中的至少一种)外观的方位。例如,方位标记具有第一部分,该第一部分小于周向180°地覆盖周向表面,并沿着电极的纵向轴线具有预定(例如,已知的或固定的中的至少一种)长度(在更具体的示例中为小于总长度),针对其余圆周,在一个示例中,方位标记在沿着纵向轴线的另一位置处还具有第二部分,该第二部分至少在圆周的其余部分覆盖电极的表面。通常通过从二维医学图像中的至少一个提取方位标记的图像外观并将其与方位标记的预定的(例如,已知的或固定的中的至少一种)图像外观相比较来确定电极的方位。例如,如果第一部分在特定的成像平面中二维成像,例如在与电极的纵向轴线相交并且表示方位标记的对称镜面的平面中二维成像,则第一部分具有预定的厚度。第一部分的图像外观的厚度可以定义方位标记和电极的预定方位。然后该方法基于旋转方位数据和预定方位数据确定最佳方位数据。
最佳方位数据描述(例如,定义或表示)与满足关于预定旋转方位的预定条件的电极的旋转方位(在二维医学图像的图像平面中)相关联的二维医学图像。在本公开内容中,该二维医学图像也被称为最佳二维医学图像,具体而言,因为当与预定方位相比时,其与方位标记的所谓的最佳旋转方位相关联。例如,通过将二维医学图像中至少一个中(例如,在每一个二维医学图像中)的方位标记的图像外观与方位标记的预定图像外观进行比较来确定最佳二维图像,以便确定方位标记的图像外观是否满足相对于预定图像外观的预定条件,例如其在至少一个二维医学图像中的厚度是否至少在等于其预定厚度的预定(例如,已知的或固定的中的至少一种)极限内。如果确定二维医学图像中的至少一个中的方位标记的图像外观满足预定条件,则相关联的二维医学图像被确定为最佳二维医学图像。如果以这种方式分析多于一个的二维医学图像,并且如果确定二维医学图像中的另一个能够更好地满足预定条件,则该另一个二维医学图像被确定为最佳二维医学图像。
在(例如第三)示例性步骤中,获取描述(例如,定义或表示)解剖结构的一组断层医学图像的断层图像数据。定义断层图像数据的图像内容的位置被定义在三维中(即,在三维参考系中),并且在执行根据第一方面的方法之前例如生成断层图像数据。例如,断层图像数据可以以用于规划设想的医学过程的计划图像数据集的形式存在。断层图像数据允许从包括与最优二维医学图像相关联的视角的不同视角生成解剖结构的切片图像。可以通过计算机x射线断层摄影术、医学共振成像、超声成像、正电子发射断层摄影术、医学共振扩散张量成像或从旋转图像数据重建三维图像体中的至少一个来生成断层图像数据。
在(例如第四)示例性步骤中,确定电极方位数据,其描述在定义断层医学图像数据的参考系中的电极的旋转方位。基于旋转图像数据、断层图像数据以及旋转方位数据来确定电极方位数据。例如,将最佳二维医学图像的参考系中(其中定义了定义旋转图像数据的图像内容的位置)的电极的旋转方位变换到断层图像数据的参考系(其中定义了定义断层图像数据的图像内容的位置)中。这种变换是预定的(例如,已知的或固定的中的至少一种),并且该变换基于两个参考系之间的预定的(已知的或固定的中的至少一种)空间关系。该变换例如是线性变换且可以通过矩阵乘法来体现,并且该变换可以是例如两个参考系之间的基底变换。
在一个示例中,该方法可以包括以下步骤:确定断层图像数据的参考系中的电极的位置和方位,以及确定解剖结构是否被电极的预定触点的方向特性所覆盖,即被用于电刺激解剖结构的电场覆盖,该电场具有例如预定的场强度以及相对于定向触点的位置和解剖结构的位置中的至少一个的空间关系。为此,还可以通过该方法获取包括描述电极的预定的(例如,已知的或固定的中的至少一种)几何结构(具体地,方位标记与期望的定向触点之间的空间关系)的结构化数据的电极模板数据,并且将该电极模板数据用作确定定向触点与解剖结构的位置之间的空间关系的基础。解剖结构的位置可以由解剖结构的几何或物理重心定义。该方法随后还可以确定由电场覆盖的解剖结构的整个体积的百分比,以及该百分比是否在用于执行解剖结构的电刺激的期望范围内。
根据第一方面的方法例如包括确定投影图像数据的步骤。基于(例如从)断层医学图像数据来确定投影图像数据,并且投影图像数据描述(例如,定义或表示)从断层图像数据以与最佳二维医学图像相关联的成像视角(例如相对于解剖结构,具体为解剖结构的位置)生成(例如,合成)的投影断层医学图像。换句话说,投影断层医学图像至少在与最佳二维医学图像相同的(即,至少在预定程度上相似的)外观的预定限制(例如精度限制)内描述(例如,示出或描绘)解剖结构。投影的医学图像可以是数字重建的射线照片。投影图像数据例如通过确定变换来确定,例如,定义的旋转图像数据的参考系(其中定义了定义旋转图像数据的图像内容的位置)与断层图像数据的参考系(其中定义了定义断层图像数据的图像内容的位置)之间的位置变换,然后通过将该变换应用于由旋转方位数据描述的旋转方位来确定电极方位数据。这种变换是预定的(例如,已知的或固定的中的至少一种),并且该变换基于两个参考系之间的预定的(已知的或固定的中的至少一种)空间关系。该变换例如是线性变换且可以通过矩阵乘法来体现,并且该变换例可以是例如两个参考系之间的基底变换。
然后,该方法可以包括确定投影方位数据的步骤。投影方位数据描述(例如,定义或表示)投影医学图像的参考系中的电极的旋转方位。基于投影图像数据、旋转图像数据以及旋转方位数据,例如通过将先前确定的电极的旋转方位从二维医学图像的参考系(其中定义了定义旋转图像数据的图像内容的位置)变换到投影断层医学图像的参考系(其中定义了定义投影断层医学图像的图像内容的位置)中,来确定投影方位数据。这种变换是预定的(例如,已知的或固定的中的至少一种),并且该变换基于两个参考系之间的预定的(已知的或固定的中的至少一种)空间关系。该变换例如是线性变换且可以通过矩阵乘法来体现,并且该变换例可以是例如两个参考系之间的基底变换。因为与投影断层图像相关联的成像视角是已知的,所以可以将投影断层摄影图像中的电极的旋转方位变换成整组断层医学图像中的电极的旋转方位,使得例如结合关于电极的几何结构的已知信息,也可以确定在断层图像数据的参考系中的电极的方位以及例如电极的位置。因此,根据第一方面的方法例如基于投影方位数据和投影图像数据来确定电极方位数据。
在一个示例中,针对与一系列二维医学图像相关联的所有的成像视角来确定投影图像数据,并且投影图像数据描述对于成像视角中的每一个的投影断层医学图像,其中基于与最佳二维图像最匹配的投影断层摄影医学图像来确定投影方位数据,该匹配例如由以下至少一个确定:
-比较最佳二维医学图像和投影断层医学图像之间的颜色值(例如,灰度值);然后,例如,可以在两个图像之间建立颜色值的相关关系,可能将相邻图像元素(如像素或体素)的空间关系和它们的相关颜色值考虑在内;
-运用融合算法(图像融合算法,例如刚性或弹性融合算法)以将最佳二维医学图像配准到投影断层医学图像中的每一个中,并且从投影断层医学图像中选择仅与最佳融合结果相关联的投影断层医学图像,用于确定投影方位数据。最佳融合结果可以通过融合算法计算的度量来定义,该度量指示数据集之间的最佳对应关系(例如,诸如图像内容的几何特征之间的对应关系,如通过比较两个方面的对应关系,即一方面是一个图像的图像内容的组成部分相对于另一个图像的位置关系,另一方面是再一个图像的图像内容的组成部分相对于又一图像的位置关系)。
在根据第一方面的方法的示例中,已经通过应用基于x射线的成像模态生成了旋转图像数据和断层图像数据,例如已经通过应用旋转常规x射线成像,如旋转放射线摄影术生成了旋转图像数据,并且通过应用计算机x射线断层摄影术已经生成了断层图像数据。具体而言,定义旋转图像数据的图像内容的颜色值(例如,灰度值)与定义断层图像数据的图像内容的颜色值(例如,灰度值)之间的对应关系是预定的(例如,已知的或固定的中的至少一种),并且可以用作根据第一方面的方法的输入。这种对应关系通常从特定材料(如形成例如解剖结构的生物组织)的x射线吸收的物理学中得知。
在根据第一方面的方法的另一示例中,通过应用基于x射线的成像模态已经生成了旋转图像数据,并且通过应用不涉及应用x射线的成像模态已经生成了断层图像数据,例如,应用旋转常规x射线成像(如在前面的示例中)已经生成了旋转图像数据,并且通过应用磁共振成像或正电子发射成像(PET)或超声成像(超声波扫描术)或磁共振扩散张量成像(MR-DTI)已经生成了断层图像数据。为了使定义旋转图像数据的图像内容的颜色值(例如,灰度值)与定义断层图像数据的图像内容的颜色值(例如,灰度值)具有可比性,在一个示例中可以将分别与旋转图像数据和断层图像数据相关联的两个不同颜色值标度之间的映射输入到该方法中。这种映射可以被包括在描述解剖结构的多模式图谱的图谱数据中。
在本示例中,该方法因此可以包括以下步骤:
在不涉及应用x射线的成像模态中以及在基于x射线的断层摄影成像模态中(例如由定义与这两种模态中的解剖结构对应的图像组成部分的如灰度值的颜色值来定义),获取描述解剖结构的模型(例如,基于图像的模型)的图谱数据和关于解剖结构的图像外观的信息,其中基于x射线的成像模态是例如计算机x射线断层摄影术;
基于断层图像数据和图谱数据确定描述在基于x射线的断层摄影成像模态中的解剖结构的图像外观的变换的外观数据(例如由定义与基于x射线的断层摄影成像模态中的解剖结构对应的图像组成部分的如灰度值的颜色值来定义)。
然后可以基于变换的外观数据进一步确定电极方位数据,并且例如可以基于(例如,从)变换的外观数据进一步确定投影图像数据。
图谱数据可以描述解剖结构的多模态图谱,其中存储了用每个不同成像模态生成的解剖学身体部位的模型以及模型的解剖学对应部分之间的变换规则。这允许将断层图像数据变换成与用于生成该断层图像数据的成像模态不同的成像模态,例如使医学图像数据与该不同的成像模态相比较。变换规则可以基于针对每个模型存储的组织类别信息,其描述相应各个成像模态中解剖学身体部位的组成部分的图像外观(例如,如多色值或灰度值的颜色值)。
此外,可能已经从多个患者的解剖结构的医学图像生成了图谱数据。可选地,可能已经仅从生成旋转图像数据和断层图像数据的特定患者的解剖结构的至少一个医学图像生成图谱数据,即模型可以是个别患者图谱的一部分。
图谱数据包括例如在三维坐标系(代表用于定义包含在图谱数据中的位置信息的参考系)中定义的位置信息。例如,已经从解剖结构的断层图像中生成了图谱数据。
变换的外观数据可以通过将图谱数据与断层图像数据匹配来确定,例如通过确定变换,例如断层图像数据的参考系(其中定义了定义断层图像数据的图像内容的位置)与图谱数据的参考系(其中定义了定义图谱数据的图像内容的位置)之间的位置变换,或者通过例如执行断层图像数据与图谱数据之间的颜色值匹配。位置变换可以是预定的(例如,已知的或固定的中的至少一种)或者作为根据第一方面的方法的一部分而建立,例如通过将融合算法(如刚性或弹性融合算法)应用于断层图像数据和图谱谱数据。
在根据第一方面的方法的一个示例中,电极包括至少两个彼此间隔开的定向触点,并且其中使用二维医学图像中的至少两个定向触点之间的至少两个间隔中的每一个的至少一部分的图像外观来验证由旋转方位数据描述的旋转方位。
第二方面,本发明涉及一种计算机程序,其在至少一台计算机(例如,计算机)的至少一个处理器(例如,处理器)上运行时或者当被加载到至少一台计算机的至少一个存储器(例如,存储器)时,使得至少一台计算机执行根据第一方面的上述方法。
第三方面,本发明涉及一种存储有根据第二方面的程序的非暂时性计算机可读程序存储介质。
第四方面,本发明涉及至少一台计算机(例如,计算机),包括至少一个处理器(例如,处理器)和至少一个存储器(例如,存储器),其中根据第二方面的程序在处理器上运行或者被加载到存储器中,或者其中至少一台计算机可操作地耦合到根据第三方面的程序存储介质以执行存储在程序存储介质上的程序。
在第五方面,本发明涉及一种承载表示根据第二方面的程序的信息的(物理的,例如电的,例如技术上生成的)信号波,例如数字信号波。
第六方面,本发明涉及一种用于确定电极的方位的医学系统,电极被配置为电刺激患者的解剖结构并且包括旋转方位标记,该系统包括:
a)根据前述权利要求的至少一台计算机;以及
b)至少一个电子数据存储设备,其存储旋转图像数据或断层图像数据或者在适用的情况下的图谱数据中的至少一个,
其中至少一台计算机可操作地耦合到至少一个电子数据存储设备,以用于从至少一个数据存储设备获取旋转图像数据或断层图像数据或在适用的情况下的图谱数据中的至少一个。
在一个示例中,根据第六方面的系统还包括二维医学成像装置,其中二维医学成像装置可操作地耦合到至少一台计算机,以允许至少一台计算机从二维医学成像装置接收对应于旋转图像数据的至少一个电子信号。
一般而言,本发明不涉及或例如包含或涵盖侵入式步骤(invasive steps),该侵入式步骤表现出要利用专业医学知识来执行对身体的实质性物理干扰,并且即使在运用所需的专业护理及专业知识执行时也具有实质性健康风险。例如,本发明不包括用紫外线照射解剖学身体部分和/或患者身体的步骤,这使得其不包括人体或动物体的任何治疗步骤,例如其不包括任何治疗或手术步骤。更具体地,本发明不涉及或特别包含或涵盖任何手术或治疗活动。相反,本发明适用于读取和处理数据以及操作或控制计算机来执行使计算机执行根据第一方面的数据处理方法的程序。仅凭此原因,通过实施本发明,不需要或未暗示手术或治疗活动,尤其是手术或治疗步骤。
为了在技术上适宜和/或可行的情况下形成新的实施例,将本发明的一个或多个实施例或方面的一个或多个特征进行组合,这在本发明的范围内。具体地,具有与另一实施例的另一特征相同或相似的功能的一个实施例的特征可与该另一特征交换,并且例如将附加功能添加到另一实施例的一个实施例的特征可以被添加到该另一实施例。
定义
在本节中,提供了本公开中使用的特定术语的定义,其也形成本公开的一部分。
根据本发明的方法是例如一种计算机实现的方法。例如,根据本发明的方法的所有步骤或仅仅一些步骤(即,少于步骤总数)可以由计算机(例如,至少一台计算机)执行。计算机实现的方法的一个实施例是使用计算机来执行数据处理方法。计算机实现的方法的实施例是关于计算机操作的方法,以操作计算机来执行该方法的一个、多个或所有步骤。
计算机例如包括至少一个处理器和例如至少一个存储器,以便例如电子地和/或光学地(在技术上)处理数据。处理器例如由作为半导体的物质或组合物制成,例如至少部分地n型和/或p型掺杂半导体,例如II-,III-,IV-,V-,VI-半导体材料中的至少一个,例如(掺杂的)硅和/或砷化镓。所描述的计算步骤例如由计算机执行。确定步骤或计算步骤是例如在技术方法的框架内,例如在程序的框架内确定数据的步骤。计算机是例如任何种类的数据处理设备,例如电子数据处理设备。计算机可以是通常被认为是这样的设备,例如台式电脑、笔记本电脑、上网本等,但也可以是任何可编程装置,例如移动电话或嵌入式处理器。计算机可以例如包括“子计算机”的系统(网络),其中每个子计算机本身代表计算机。术语“计算机”包括云计算机,例如云服务器。术语“云计算机”包括云计算机系统,其例如包括至少一个云计算机的系统,以及例如多个可操作地互连的云计算机,如服务器农场。这样的云计算机优选连接到如万维网(WWW)的广域网,并且位于所有连接到万维网的计算机云中。这样的基础设施被用于“云计算”,其描述了计算、软件、数据访问和存储服务,其不要求最终用户知道提供特定服务的计算机的物理位置和/或配置。例如,术语“云”在这方面被用作互联网(万维网)的隐喻。例如,云提供计算的基础设施即服务(IaaS)。云计算机可以用作用于执行本发明的方法的操作系统和/或数据处理应用程序的虚拟主机。云计算机是例如由亚马逊网络服务TM提供的弹性计算云(EC2)。计算机例如包括接口以便接收或输出数据和/或执行模数变换。数据例如是代表物理性质的数据和/或从技术信号中产生的数据。技术信号例如通过(技术)检测设备(例如用于检测标记设备的设备)和/或(技术)分析设备(例如用于执行(医学)成像方法的设备)来生成,其中技术信号是例如电信号或光信号。技术信号例如表示由计算机接收或输出的数据。计算机优选可操作地耦合到显示设备,该显示设备允许将由计算机输出的信息显示给例如用户。显示设备的一个示例是可以用作用于导航的“护目镜”的增强现实设备(也被称为增强现实眼镜)。这种增强现实眼镜的具体示例是谷歌眼镜(Google Glass,Google Inc.的商标)。增强现实设备既可以用于通过用户交互将信息输入到计算机中,也可以用于显示计算机输出的信息。显示设备的另一个示例是标准计算机监视器,包括例如可操作地耦合到计算机的液晶显示器,用于接收来自计算机的显示控制数据,以产生用于在显示设备上显示图像信息内容的信号。这种计算机监视器的具体实施例是数字灯箱。监视器也可以是便携式设备的监视器,例如手持设备,手持设备例如是智能电话或个人数字助理或数字媒体播放器。
在本发明的框架内,计算机程序单元可以通过硬件和/或软件(这包括固件、常驻软件、微代码等)来体现。在本发明的框架内,计算机程序单元可以采用计算机程序产品的形式,该计算机程序产品可以通过计算机可用的,例如计算机可读的数据存储介质来实现,该计算机可读的数据存储介质包括体现在该存储介质中以在指令执行系统中使用或与指令执行系统结合使用的计算机可用的(例如计算机可读的)程序指令、“代码”或“计算机程序”。这样的系统可以是计算机;计算机可以是包括用于执行根据本发明的计算机程序单元和/或程序的装置的数据处理设备,例如数据处理设备包括执行计算机程序单元的数字处理器(中央处理单元或CPU)以及可选地包括易失性存储器(例如随机存取存储器或RAM)的数据处理设备,该易失性存储器用于存储用于执行计算机程序单元元件和/或通过执行计算机程序单元产生的数据。在本发明的框架内,计算机可用的,例如计算机可读的数据存储介质可以是任何数据存储介质,其可以包括、存储、通信、传播或传输程序以供在指令执行系统、装置或设备中使用或与指令执行系统、装置或设备结合使用。计算机可用的,例如计算机可读的数据存储介质可以例如是但不限于电子、磁性、光学、电磁、红外或半导体系统、装置或设备或例如互联网的传播介质。因为程序可以例如通过光学扫描纸张或其他合适的介质来电子地捕获,然后编译,解释或以适当的方式处理,因此计算机可用或计算机可读的数据存储介质甚至可以例如是在其上打印程序的纸张或另一合适介质。数据存储介质优选为非易失性数据存储介质。这里描述的计算机程序产品和任何软件和/或硬件形成了用于在示例实施例中执行本发明的功能的各种装置。计算机和/或数据处理设备可以例如包括指导信息设备,其包括用于输出指导信息的装置。指导信息可以例如通过视觉指示装置(例如,监视器和/或灯)视觉地输出给用户,和/或通过声学指示装置(例如,扬声器和/或数字语音输出设备)声学地输出给用户,和/或通过触觉指示装置(例如,安装在仪器中的振动元件)触觉地输出给例如用户。为了实现本文的目的,计算机是技术性计算机,其例如包括技术性的,例如有形的部件,例如机械和/或电子部件。本文中提到的任何设备都是技术性的,以及例如有形的设备。
表述“获取数据”例如包括(在计算机实现的方法的框架内)其中数据由计算机实现的方法或程序确定的情况。确定数据例如包括测量物理量并将测量值转换成数据,例如转换成数字数据,和/或利用计算机以及例如在根据本发明的方法的框架内计算数据。“获取数据”的含义例如还包括其中数据被计算机实现的方法或程序例如从另一程序、先前的方法步骤或数据存储介质接收或读取,例如用于被计算机实现的方法或程序进一步处理的情况。待获取的数据的生成可以但不必是根据本发明的方法的一部分。因此,表述“获取数据”也可以例如意为等待接收数据和/或正在接收数据。接收到的数据可以例如通过接口来输入。表述“获取数据”还可以表示计算机实现的方法或程序执行步骤以(主动地)从数据源,例如数据存储介质(例如,ROM、RAM、数据库、硬盘驱动器等)或通过接口(例如,来自于另一台计算机或网络)接收或恢复数据。分别由所公开的方法或设备获取的数据可以从位于数据存储设备中的数据库获得,该数据存储设备可操作地连接到用于数据库和计算机之间的(例如从数据库到计算机的)数据传输的计算机。计算机获取数据以用作确定数据的步骤的输入。所确定的数据可以再次输出到相同或另一个数据库中被存储以供以后使用。数据库或用于实现所公开的方法的数据库可以位于网络数据存储设备或网络服务器(例如,云数据存储设备或云服务器)或本地数据存储设备(例如可操作地连接到执行所公开的方法的至少一台计算机的大容量存储设备)。在获取步骤之前,通过执行附加步骤,可以使数据“随时可用”。根据该附加步骤,生成数据以便其被获取。(例如通过分析装置)检测或捕获数据。可选地或附加地,根据附加步骤,例如经由接口输入数据。能够将所生成的数据输入(例如输入至计算机中)。根据附加步骤(其在获取步骤之前),还可以通过执行将数据存储在数据存储介质(例如,ROM、RAM、CD和/或硬盘驱动器)中的附加步骤来提供数据,以使得这些数据在根据本发明的方法或程序的框架内随时可用。因此,“获取数据”的步骤也可以涉及命令设备获取和/或提供待获取的数据。特别地,获取步骤不涉及侵入式步骤,该侵入式步骤会表现出对身体的实质性物理干扰,需要用专业医学知识来执行,并且即使在使用所需的专业护理及专业知识执行时也具有实质性健康风险。具体而言,获取数据的步骤,例如确定数据,不涉及手术步骤,并且尤其不涉及使用手术或疗法来治疗人体或动物体的步骤。为了区分本方法使用的不同数据,将数据表示为(即,称为)“XY数据”等,并且根据它们描述的信息来定义,然后优选地将其称为“XY信息”等。
在医学领域中,使用成像方法(也称为成像模态和/或医学成像模态)来生成人体的解剖结构(如软组织、骨骼、器官等)的图像数据(例如,二维或三维图像数据)。术语“医学成像方法”被理解为意指(有利地基于设备的)成像方法(例如所谓的医学成像模态和/或放射成像方法),如计算机断层摄影术(CT)和锥形束计算机断层摄影术(CBCT,例如容积CBCT)、x射线断层摄影术、磁共振断层摄影术(MRT或MRI)、常规x射线、超音波和/或超声波检查以及正电子发射断层摄影术。例如,医学成像方法由分析设备执行。医学成像方法应用的医学成像模态的实例是:x射线成像、磁共振成像、医学超声波检查或超声波、内窥镜检查、弹性成像、触觉成像、热成像、医学摄影以及核医学功能成像技术,如维基百科所述的如正电子发射断层摄影术(PET)和单光子发射计算机断层摄影术(SPECT)。这样生成的图像数据也被称为“医学成像数据”。分析设备例如用于在基于装置的成像方法中生成图像数据。成像方法例如用于医学诊断,以分析解剖体以便生成由图像数据描述的图像。成像方法也例如用于检测人体内的病理变化。然而,解剖结构中的一些变化,例如结构(组织)中的病理变化可能无法检测到,并且例如在由成像方法生成的图像中可能是不可见的。肿瘤代表了解剖结构中的变化的一个例子。如果肿瘤生长,那么可以说它代表了扩张的解剖结构。这种扩张的解剖结构可能无法检测到,例如,只有扩张解剖结构的一部分可以被检测到。当造影剂用于浸润肿瘤时,原发性/高级别的脑肿瘤例如通常在MRI扫描中可见。MRI扫描代表成像方法的示例。在这种脑肿瘤的MRI扫描的情况下,MRI图像中的信号增强(由于造影剂浸润肿瘤)被认为代表实体肿瘤块。因此,肿瘤是可检测的,并且例如在由成像方法生成的图像中可辨别。除了被称为“增强”肿瘤的这些肿瘤之外,据认为大约10%的脑肿瘤在扫描中不可辨别,并且例如对于查看由成像方法生成的图像的用户不可见。
图像融合可以是弹性图像融合或刚性图像融合。在刚性图像融合的情况下,2D图像的像素之间和/或3D图像的体素之间的相对位置是固定的,而在弹性图像融合的情况下,允许相对位置改变。
在本申请中,术语“图像变形”也被用作术语“弹性图像融合”的替代,但是具有相同的含义。
弹性融合变换(例如弹性图像融合变换)例如被设计成使得实现从一个数据集(例如第一数据集,例如第一图像)到另一个数据集(例如第二数据集,例如第二图像)的无缝过渡。例如,变换被设计为使第一和第二数据集(图像)中的一个变形,这样使得相应的结构(例如,对应的图像元素)被设置在与第一和第二图像中的另一个相同的位置处。从第一和第二图像之一变换的变形图像(变换图像)例如与第一和第二图像中的另一个尽可能相似。优选地,应用(数值)优化算法以便找到导致最佳相似度的变换。相似度优选通过测量相似性来测量(在下文中也称为“相似性量度”)。优化算法的参数例如是变形场的向量。这些向量由优化算法确定,以便产生最佳的相似度。因此,最佳相似度表示优化算法的条件,例如约束条件。向量的基部例如位于待变换的第一和第二图像之一的体素位置处,并且向量的顶点位于变换的图像中的对应体素位置处。优选提供多个这些向量,例如多于二十或一百或一千或一万等等。优选地,对变换(变形)有(其他)约束,例如为了避免病态变形(例如,通过变换将所有体素移位到相同的位置)。这些约束包括例如变换是规则的约束,例如这意味着从变形场(例如,向量场)的矩阵计算的雅可比行列式大于零,并且还包括变换的(变形的)图像不是自相交的约束,并且例如变换的(变形的)图像不包含缺陷和/或断裂。约束包括例如这样的约束,即如果规则网格与图像同时变换并且以相应的方式变换,则不允许网格在其任何位置处进行交叠。例如,迭代地解决该优化问题,例如通过优化算法,其例如是一阶优化算法,例如梯度下降算法。优化算法的其他示例包括不使用导数的优化算法,如下坡单纯形算法,或者使用高阶导数的算法,如牛顿类算法。优化算法优选地执行局部优化。如果存在多个局部最优值,则可以使用如模拟退火或遗传算法的全局算法。在线性优化问题的情况下,例如可以使用单纯形法。
在优化算法的步骤中,体素例如在一定方向上移位一个量值,以使得相似度增加。这个量值优选地小于预定限度,例如小于图像直径的十分之一或百分之一或千分之一,并且例如大约等于或小于相邻体素之间的距离。例如由于大量的(迭代)步骤,可以实现大变形。
所确定的弹性融合变换可以例如用于确定第一和第二数据集(第一和第二图像)之间的相似度(或相似性量度,参见上文)。为此,确定了弹性融合变换与恒等变换之间的偏差。例如可以通过确定弹性融合变换与恒等变换的行列式之间的差异来计算偏差度。偏差越高,相似性越低,因此可以使用偏差度来确定相似性量度。
例如可以基于第一和第二数据集之间的确定的相关性来确定相似性量度。
优选地,获取描述(例如,定义,更具体地表示和/或是)解剖学身体部位的一般三维形状的图谱数据。因此该图谱数据代表解剖学身体部位的图谱。图谱通常由多个对象的通用模型(generic model)组成,其中对象的通用模型一起形成复杂的结构。例如,图谱构成患者的身体(例如身体的一部分)的统计模型,该统计模型从集合了来自多个人体的解剖信息生成,例如从包含这些人体的图像的医用图像数据来生成。原则上,图谱数据因此表示对于多个人体的这种医学图像数据的统计分析的结果。该结果可以作为图像输出-因此图谱数据包含医学图像数据或可与其相比较。这种比较可以例如通过应用在图谱数据和医学图像数据之间进行图像融合的图像融合算法来执行。比较的结果可以是图谱数据和医学图像数据之间的相似性量度。图谱数据包括可以与例如医学图像数据中包含的位置信息匹配(例如通过应用弹性或刚性图像融合算法)的位置信息,从而例如将图谱数据与医学图像数据进行比较,以便确定医学图像数据中与由图谱数据定义的解剖结构相对应的解剖结构的位置。
其中解剖结构用作生成图谱数据的输入的人体优选地共享如性别、年龄、种族、身体尺度(例如尺寸和/或质量)和病理状态中的至少一个的共同特征。解剖信息描述例如人体的解剖结构,并且例如从关于人体的医学图像信息中提取。例如,股骨的图谱可以包括头部、颈部、身体、大转子、小转子和下肢作为一起构成完整结构的对象。例如,大脑的图谱可以包括端脑、小脑、间脑、脑桥、中脑和髓质作为一起构成复杂结构的对象。这种图谱的一个应用是医学图像的分割,其中图谱与医学图像数据相匹配,并且将图像数据与匹配的图谱进行比较,以便将图像数据的点(像素或体素)分配给匹配的图谱的对象,由此将图像数据分割成各个对象。
附图描述
在下文中,参照描述本发明的具体实施例的附图来描述本发明。然而,本发明的范围不限于在附图的上下文中公开的具体特征,其中:
图1是示出根据第一方面的公开的方法的基本步骤的流程图;
图2示出了具有定向标记的定向DBS电极;
图3示出了根据第一方面的公开的方法的实施例;以及
图4示出了从二维医学图像确定电极的旋转方位的具体方法。
图1是示出了根据第一方面的公开的方法的基本步骤的流程图,在图1的说明性示例中,该方法开始于获取旋转图像数据的步骤S1.1。在随后的步骤S1.2中,确定旋转方位数据。然后步骤S1.3继续获取断层图像数据,接着步骤S1.4,其包括基于在前的步骤S1.1到步骤S1.3中获取和确定的数据来确定电极方位数据。
图2示出了定向DBS电极1(也称为引导体),包括其上设置有触点2的圆柱形且细长的主体。触点2包括定向触点8,定向触点8不沿着基体的整个圆周延伸,而是设计成仅在跨度小于圆周的360°的三维扇区内发射电场。在图2的示例中,每个定向触点8都至少覆盖大约圆周的60°并且通过在定向触点8之间形成空间的狭缝7彼此分开。该空间可以是空的或者填充有非导电材料,例如制造电极1的主体的材料。狭缝7平行于电极1的主体的纵向轴线延伸并因此位于形成电极1的基本形状的圆柱体的母线上。电极1还包括方位标记3,该方位标记3包括具有第一部分5和第二部分4的标记带以及标记窗5。方位标记的旋转方位被定义为指向电极1的主体的外部的方向并垂直于第一部分5的中心的表面的向量,即,其基部(base)位于处于第一部分1的对称镜面内的点中,电极的主体的纵向轴线也位于该平面内。标记带由铂(Pt)或至少由包含铂(Pt)的材料(例如合金)制成,并且比标记窗口5更不透射线。在一个示例中,标记带也比电极1的主体更不透射线。最佳的二维医学图像是其中在二维医学图像中第一部分5的描绘(图像外观)最薄的放射线摄影术,即其由最小数量的图像元素(如像素)限定。
如图3所示,可以通过使用C型臂x射线装置对患者体内的解剖结构进行成像来获取旋转图像数据。然后,例如通过对二维医学图像应用分割算法,从二维医学图像(2D x射线图像/射线摄像术)中提取电极1(引导体)的图像描绘。在图3左下角的框中,最右边的二维医学图像被选为最佳的二维医学图像,因为它包含了标记带的第一部分5的最薄的描绘。如框中的箭头所示,电极1(以及因此方位标记3)的旋转方位被确定为位于最右边的二维医学图像的图像平面中并指向该图像的右下方。
又如图3所示,获取与断层图像数据或变换的外观数据对应的3D(三维)图像。可以通过CT或者从2D x射线图像重建3D图像来生成断层图像数据。可选地或附加地,可以基于已经变换为CT外观的MR断层摄影来获取3D图像,从而生成也是三维图像数据集的变换的外观数据。MR断层摄影的变换体现了断层图像数据,并通过获取多模式图谱数据并在图谱数据和MR断层摄影之间进行匹配而被变换为CT灰度值。该匹配在图3中被称为通用图谱匹配。基于组织依赖性MR灰度值标度与图谱中包含的组织依赖性CT灰度值标度之间的对应关系,通过用相应CT图像的灰度值代替MR断层摄影的灰度值(例如体素取向),从MR断层摄影术模拟CT(计算机x射线断层摄影)。
从3D图像中产生DRR。在一个示例中,从朝向对应于与最佳二维医学图像相关联的视角的解剖结构的角度仅产生一个DRR。在另一示例中,针对已经采集或生成二维医学图像的所有视角生成DRR。然后,选择与最佳二维医学图像最佳匹配的DRR用于进一步处理。例如通过将图像数据集彼此融合,建立2D图像的参考系与一个或多个DRR之间的变换(在图3中称为2D-3D配准)。基于该变换,可以基于相对于与生成2D x射线相关联(即,使用时)的解剖结构的成像视角的信息,在2D x射线图像的参考系中确定电极的旋转方位。然后可以通过应用2D-3D配准将旋转方位从2D x射线图像的参考系变换到一个或多个DRR的参考系中。然后,基于用于生成一个或多个DRR的成像视角的信息,如图4的上半部分所示,可以在3D图像的参考系中确定电极的旋转方位。然后,例如在3D图像的参考系中,这样确定的旋转方位例如可以与电极的预定(例如,期望的)旋转方位进行比较。如果所确定的旋转方位至少在基本上等于预定旋转方位的预定距离内,则可以根据需要开始所设想的医学过程。如果所确定的旋转方位不是至少在基本上等于预定旋转方位的预定距离内,则可以调节电极的旋转方位直到与预定旋转方位一致。
如图4所示,方向触点8之间的狭缝7可以用于支持旋转方位的确定,例如用于验证从关于方位标记的图像外观分析二维医学图像接收到的结果。图4的右侧部分示出用于生成二维医学图像的x射线是否可能以有限的概率穿过两个狭缝7并因此在二维图像中生成相应的标记,因为制造定向触点8的材料与制造狭缝7的材料或填充它们的材料的不透射线性不同。例如,穿过两个狭缝7的x射线可以产生指示与定向触点8的图像外观相比较的半透明性的标记。从而,在图1和4的示例中,可以以大约60°的精度接收从方位标记的图像外观确定的旋转方位是否有效的附加指示。
所公开的方法可以进一步总结如下。
该方法使用由x射线系统在如患者头部的解剖结构周围旋转时拍摄的图像。可以使用来自任何断层扫描仪的头部的3D(体)图像。此外,使用算法将2D x射线图像配准到3D图像。另一种算法检测2D x射线图像中的引导体并且定义引导体的方位标记相对于图像平面(即相对于2D x射线图像的参考系)的角度。角度信息通过配准从2D图像空间映射回3D图像(解剖结构)空间。计算机系统能够显示混合在任何以所计算的配准为模的2D x射线图像上的3D图像的模拟投影,以允许用户快速检查配准精度。
可以考虑的几个实施方式:
·使用一个x射线,两个x射线或多个(旋转的)x射线。使用的图像越多,角度的确定就越鲁棒。
·使用标准CT作为3D图像。它们被校准到Hounsfield值。数字重建摄像照片(DRR)可以从头开始计算。
·将x射线扫描仪直接重建的CT用作3D图像。那么不需要额外的扫描。可能有必要重新校准到HU值,这可以在通用图谱的帮助下完成。
·将手术前MRI用作3D图像。然后,需要对其进行预处理并将其变换为相应CT表示的Hounsfield值。这可以在用图谱的帮助下完成。
·数字重建射线照片(DRR)的计算可以在图形卡上完成。
·定向引导体在分段触点之间有狭缝,这些夹缝在x射线中在某些角度可见。这些可以额外用于提高角度精度。
·引导体可以通过x射线或3D图像中的算法进行分割。
·标记方向可以通过模板匹配或机器学习找到。
·可以从存储在一起/与x射线图像信息相关联的DICOM信息中获取2D x射线图像的照相机参数(如成像视角)。
·如果不能从DICOM获得,则例如从2D图像中的引导体几何结构估计X射线的照相机参数。
除了旋转方位之外,还可以例如基于对电极的至少一部分(例如方位标记)的图像外观与电极模板数据进行比较,来在x射线(制造商,型号)中检测引导体类型。
所公开的方法旨在快速且可追踪地传输从2D图像到3D解剖结构的关于DBS引导体方位的信息。此传输很快,因为系统包含允许配准2D和3D坐标的自动算法,否则这将是需要用户手动完成的繁琐步骤(即为了匹配2D和3D图像而手动采用若干自由度,如移位/旋转/变焦/焦距)。此外,2D-3D配准使得用户可以追踪最终的角度分配。用户可以立即检查计算出的角度与2D和3D图像的关系。
Claims (17)
1.一种计算机实现的用于确定电极的方位的医学数据处理方法,所述电极被配置为电刺激患者的解剖结构并且包括旋转方位标记,所述方法包括在至少一个计算机的至少一个处理器上执行步骤:
a)在所述至少一个处理器处获取描述所述解剖结构和所述电极的二维医学图像的旋转图像数据,所述二维医学图像是在医学成像装置相对于解剖结构旋转期间用二维医学成像装置拍摄的,所述旋转图像数据针对每个二维医学图像进一步描述相对于与相应二维医学图像相关联的解剖结构的成像视角;
b)由所述至少一个处理器并且基于最佳二维医学图像及其方向标记的图像外观来确定描述二维医学图像的参考系中的电极的旋转方位的旋转方位数据,所述二维医学图像的参考系中的电极的旋转方位满足关于预定旋转方位的预定条件,此处的二维医学图像被称为所述最佳二维医学图像;
c)在所述至少一个处理器处获取描述所述解剖结构的一组断层医学图像的断层图像数据;以及
d)由所述至少一个处理器并且基于所述旋转图像数据和所述断层图像数据以及所述旋转方位数据,来确定描述断层医学图像的参考系中的电极的旋转方位的电极方位数据。
2.根据前述权利要求所述的方法,包括:
由所述至少一个处理器获取描述二维医学图像的参考系中的电极的预定旋转方位的预定方位数据;
基于所述旋转方位数据和所述预定方位数据,来确定描述与二维医学图像的参考系中的电极的旋转方位相关联的二维医学图像的最佳方位数据;
由所述至少一个处理器并且基于断层图像数据来确定投影图像数据,该投影图像数据描述以与所述最佳二维医学图像相关联的成像视角从断层图像数据生成的投影断层医学图像;
由所述至少一个处理器并且基于所述投影图像数据和所述旋转图像数据来确定投影方位数据,该投影方位数据描述投影医学图像的参考系中的电极的旋转方位。
3.根据权利要求2所述的方法,其中,电极方位数据由所述至少一个处理器基于所述投影方位数据和所述投影图像数据来确定。
4.根据权利要求2或3所述的方法,其中,针对与一系列二维医学图像相关联的所有成像视角确定投影图像数据,并且所述投影图像数据针对每个成像视角描述投影断层医学图像,并且其中基于与最佳二维图像最匹配的投影断层医学图像来确定投影方位数据,该匹配由以下至少一项确定:
-比较最佳二维医学图像和投影断层医学图像的颜色值;
-运用融合算法以将所述最佳二维医学图像配准到每个投影断层医学图像中,并且从所述投影断层医学图像中选择仅与最佳融合结果相关联的投影断层医学图像,以用于确定投影方位数据。
5.根据权利要求1所述的方法,其中,所述旋转方位数据通过以下至少一项确定:
-分割每个二维医学图像中的电极的图像外观;
-对二维医学图像的组成部分进行边缘检测;
-对二维医学图像中的电极的图像外观与描述电极的结构化数据的先前获取并预定的电极模板数据进行比较。
6.根据权利要求1所述的方法,其中,每个二维医学图像与不同的成像视角相关联。
7.根据权利要求2所述的方法,其中,确定投影图像数据包括:确定定义了旋转图像数据的参考系与定义了断层图像数据的参考系之间的位置变换,并且其中确定电极方位数据包括:将所述变换应用于由旋转方位数据描述的旋转方位。
8.根据权利要求1所述的方法,其中,旋转图像数据和断层图像数据已经通过应用基于x射线的成像模态而生成,旋转图像数据已经应用旋转射线摄影术的旋转常规x射线成像而生成,并且断层图像数据已经应用计算机x射线断层摄影术而生成。
9.根据权利要求2所述的方法,其中,旋转图像数据已经通过应用基于x射线的成像模态而生成,并且断层图像数据已经通过应用不涉及x射线的成像模态而生成,旋转图像数据已经通过应用旋转常规x射线成像而生成,并且断层图像数据已经通过应用磁共振成像而生成。
10.根据权利要求9所述的方法,包括:
在所述至少一个处理器处获取描述解剖结构的模型的图谱数据以及关于不涉及应用x射线的成像模态中的以及基于x射线的断层摄影成像模态中的解剖结构的图像外观的信息,其中所述基于x射线的成像模态是计算机x射线断层摄影术;
由所述至少一个处理器并且基于断层图像数据和图谱数据,来确定描述基于x射线的断层成像模态中的解剖结构的图像外观的变换的外观数据,
其中进一步基于变换的外观数据来确定电极方位数据,并且,进一步基于变换的外观数据来确定投影图像数据。
11.根据权利要求10所述的方法,其中,确定变换的外观数据包括:通过确定定义了断层图像数据的参考系与定义了图谱数据的参考系之间的位置变换,或者通过执行断层图像数据与图谱数据之间的颜色值匹配,来匹配图谱数据和断层图像数据。
12.根据权利要求2所述的方法,其中,所述投影医学图像是数字重建的射线照片。
13.根据权利要求1所述的方法,其中,所述电极包括至少两个彼此间隔开的定向触点,并且其中使用二维医学图像中的至少两个定向触点之间的至少两个间隔中的每一个的至少一部分的图像外观,来验证由旋转方位数据描述的旋转方位。
14.一种存储有计算机程序的非暂时性计算机可读程序存储介质,当所述计算机程序在至少一个计算机的至少一个处理器上运行时或当被加载到至少一个计算机的存储器中时,所述计算机程序使所述至少一个计算机执行根据权利要求1至13中任一项所述的方法或者产生承载代表所述程序的信息的数字信号波。
15.至少一个计算机,包括至少一个处理器和存储器,其中,根据权利要求14所述的非暂时性计算机可读程序存储介质中存储的程序在所述至少一个处理器上运行或者被加载到所述存储器中,或者其中所述至少一个计算机包括所述非暂时性计算机可读程序存储介质。
16.一种用于确定电极的方位的医学系统,电极被配置为电刺激患者的解剖结构并且包括旋转方位标记,所述系统包括:
a)根据权利要求15的所述至少一个计算机;以及
b)至少一个电子数据存储设备,其存储有旋转图像数据或断层图像数据或者只要所述程序使至少一个计算机执行根据权利要求10或11所述的方法的步骤时的图谱数据中的至少一个,
其中,所述至少一个计算机可操作地耦合到至少一个电子数据存储设备,以用于从至少一个数据存储设备获取旋转图像数据或断层图像数据或只要程序使至少一个计算机执行根据权利要求9或10所述的方法的步骤时的图谱数据中的至少一个。
17.根据权利要求16所述的系统,进一步包括二维医学成像装置,其中所述二维医学成像装置可操作地耦合到至少一个计算机,以允许所述至少一个计算机从所述二维医学成像装置接收对应于旋转图像数据的至少一个电子信号。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/EP2017/051426 WO2018137759A1 (en) | 2017-01-24 | 2017-01-24 | Determining rotational orientation of a deep brain stimulation electrode in a three-dimensional image |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108697402A CN108697402A (zh) | 2018-10-23 |
CN108697402B true CN108697402B (zh) | 2022-05-31 |
Family
ID=57956255
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201780003514.4A Active CN108697402B (zh) | 2017-01-24 | 2017-01-24 | 在三维图像中确定深度脑刺激电极的旋转方位 |
Country Status (5)
Country | Link |
---|---|
US (2) | US11317880B2 (zh) |
EP (1) | EP3376960B1 (zh) |
JP (1) | JP6675495B2 (zh) |
CN (1) | CN108697402B (zh) |
WO (1) | WO2018137759A1 (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3532866A1 (en) | 2016-10-28 | 2019-09-04 | PPG Industries Ohio, Inc. | Coatings for increasing near-infrared detection distances |
EP3881236A1 (en) | 2018-11-13 | 2021-09-22 | PPG Industries Ohio Inc. | Method of detecting a concealed pattern |
US11561329B2 (en) | 2019-01-07 | 2023-01-24 | Ppg Industries Ohio, Inc. | Near infrared control coating, articles formed therefrom, and methods of making the same |
JP7359860B2 (ja) * | 2019-03-12 | 2023-10-11 | エーエムディーティー ホールディングス インコーポレイテッド | モノスコピック放射線画像と3次元モデルの登録方法およびシステム |
CN110742652A (zh) * | 2019-09-18 | 2020-02-04 | 中国科学院西安光学精密机械研究所 | 一种对消化道肿瘤磁辅助超声图像的三维重建设备及方法 |
WO2021121631A1 (en) * | 2019-12-20 | 2021-06-24 | Brainlab Ag | Correcting segmentation of medical images using a statistical analysis of historic corrections |
CN111544020B (zh) * | 2020-04-17 | 2023-08-01 | 北京东软医疗设备有限公司 | X射线成像设备的几何校正方法及装置 |
US11238591B1 (en) | 2020-07-15 | 2022-02-01 | Taipei Medical University (Tmu) | Medical image processing system and method thereof |
TWI737404B (zh) * | 2020-07-15 | 2021-08-21 | 臺北醫學大學 | 醫療影像處理系統及方法 |
DE102020130079A1 (de) | 2020-11-13 | 2022-05-19 | Boston Scientific Neuromodulation Corporation | Verfahren und System zur Ermittlung der Rotationsorientierung einer in einem Patienten implantierten Elektrode |
EP4173590A1 (en) * | 2021-10-26 | 2023-05-03 | Ecential Robotics | Method for registering a 3d medical image with a registration phantom |
WO2024002476A1 (en) | 2022-06-29 | 2024-01-04 | Brainlab Ag | Determining electrode orientation using optimized imaging parameters |
CN116524011B (zh) * | 2023-04-26 | 2024-04-30 | 北京航空航天大学 | 一种居家环境下机器人对目标位姿的精化方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007092316A2 (en) * | 2006-02-03 | 2007-08-16 | The University Of Florida Research Foundation, Inc. | Image guidance system for deep brain stimulation |
CN102369529A (zh) * | 2009-04-02 | 2012-03-07 | 皇家飞利浦电子股份有限公司 | 用于图像引导治疗规划的自动解剖结构描绘 |
CN103959345A (zh) * | 2011-09-30 | 2014-07-30 | 博医来股份公司 | 使用颜色的剂量分布显示方法 |
Family Cites Families (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102005028746B4 (de) * | 2005-06-21 | 2018-02-22 | Siemens Healthcare Gmbh | Verfahren zum Ermitteln der Position und Orientierung eines Objekts, insbesondere eines Katheters, aus zweidimonsionalen Röntgenbildern |
US7394889B2 (en) * | 2006-05-18 | 2008-07-01 | Varian Medical Systems Technologies, Inc. | Contrast-enhanced cone beam X-ray imaging, evaluation, monitoring and treatment delivery |
US20080147086A1 (en) * | 2006-10-05 | 2008-06-19 | Marcus Pfister | Integrating 3D images into interventional procedures |
EP2594315A1 (en) * | 2007-07-10 | 2013-05-22 | Sapiens Steering Brain Stimulation B.V. | Neurostimulation system |
US9220889B2 (en) * | 2008-02-11 | 2015-12-29 | Intelect Medical, Inc. | Directional electrode devices with locating features |
EP2321002B1 (en) * | 2008-05-15 | 2014-04-23 | Intelect Medical Inc. | Clinician programmer system and method for calculating volumes of activation |
US8995731B2 (en) * | 2008-11-26 | 2015-03-31 | Medtronic, Inc. | Image-based characterization of implanted medical leads |
US8983577B2 (en) * | 2009-02-12 | 2015-03-17 | Koninklijke Philips N.V. | System for determining the orientation of a catheter |
CN102147919B (zh) * | 2010-02-10 | 2013-09-25 | 昆明医学院第一附属医院 | 一种校正术前三维图像的术中配准方法和装置 |
EP2580710B1 (en) * | 2010-06-14 | 2016-11-09 | Boston Scientific Neuromodulation Corporation | Programming interface for spinal cord neuromodulation |
WO2012009181A2 (en) * | 2010-07-16 | 2012-01-19 | Boston Scientific Neuromodulation Corporation | Systems and methods for radial steering of electrode arrays |
EP2595542A1 (en) * | 2010-07-19 | 2013-05-29 | Koninklijke Philips Electronics N.V. | 3d-originated cardiac roadmapping |
DE102010027692A1 (de) * | 2010-07-20 | 2012-01-26 | Siemens Aktiengesellschaft | Verfahren zur Bildüberwachung bei der Implantation eines Chochlea-Implantats und Röntgeneinrichtung |
US8849375B2 (en) | 2011-11-01 | 2014-09-30 | Siemens Medical Solutions Usa, Inc. | System for detecting rotation angle of a catheter in an X-ray image |
AU2013327731B2 (en) * | 2012-10-03 | 2016-09-15 | Boston Scientific Neuromodulation Corporation | Micro-electrode recording-guided implantation of directional leads |
US9406130B2 (en) * | 2012-10-26 | 2016-08-02 | Brainlab Ag | Determining an anatomical atlas |
EP2934378A4 (en) * | 2012-12-21 | 2016-08-10 | Jason Spencer | CATHETER ORIENTATION MARKER |
US10799145B2 (en) * | 2013-03-15 | 2020-10-13 | The Cleveland Clinic Foundation | Method and system to facilitate intraoperative positioning and guidance |
WO2015003748A1 (en) * | 2013-07-11 | 2015-01-15 | Brainlab Ag | Generation of a patient-specific anatomical atlas |
CN105408939B (zh) * | 2013-07-23 | 2021-08-24 | 皇家飞利浦有限公司 | 用于配准成像设备与跟踪设备的配准系统 |
JP2015181773A (ja) * | 2014-03-25 | 2015-10-22 | オリンパス株式会社 | 神経刺激用リード |
US9999772B2 (en) | 2014-04-03 | 2018-06-19 | Pacesetter, Inc. | Systems and method for deep brain stimulation therapy |
JP6376873B2 (ja) * | 2014-07-16 | 2018-08-22 | キヤノン株式会社 | 画像処理装置、画像処理方法及びプログラム |
JP6457238B2 (ja) * | 2014-10-24 | 2019-01-23 | オリンパス株式会社 | 電気刺激カテーテル |
WO2017040573A1 (en) * | 2015-09-01 | 2017-03-09 | Boston Scientific Neuromodulation Corporation | Detection of lead orientation |
-
2017
- 2017-01-24 WO PCT/EP2017/051426 patent/WO2018137759A1/en unknown
- 2017-01-24 JP JP2018544352A patent/JP6675495B2/ja active Active
- 2017-01-24 US US15/761,752 patent/US11317880B2/en active Active
- 2017-01-24 CN CN201780003514.4A patent/CN108697402B/zh active Active
- 2017-01-24 EP EP17702561.6A patent/EP3376960B1/en active Active
-
2022
- 2022-05-02 US US17/734,939 patent/US12070343B2/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007092316A2 (en) * | 2006-02-03 | 2007-08-16 | The University Of Florida Research Foundation, Inc. | Image guidance system for deep brain stimulation |
CN102369529A (zh) * | 2009-04-02 | 2012-03-07 | 皇家飞利浦电子股份有限公司 | 用于图像引导治疗规划的自动解剖结构描绘 |
CN103959345A (zh) * | 2011-09-30 | 2014-07-30 | 博医来股份公司 | 使用颜色的剂量分布显示方法 |
Also Published As
Publication number | Publication date |
---|---|
EP3376960A1 (en) | 2018-09-26 |
EP3376960B1 (en) | 2019-03-27 |
WO2018137759A1 (en) | 2018-08-02 |
JP2019511268A (ja) | 2019-04-25 |
US20220249038A1 (en) | 2022-08-11 |
US20200237326A1 (en) | 2020-07-30 |
US12070343B2 (en) | 2024-08-27 |
JP6675495B2 (ja) | 2020-04-01 |
US11317880B2 (en) | 2022-05-03 |
CN108697402A (zh) | 2018-10-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US12070343B2 (en) | Determining rotational orientation of a deep brain stimulation electrode in a three-dimensional image | |
CN110770792B (zh) | 确定临床靶体积 | |
EP3589355B1 (en) | Optimal deep brain stimulation electrode selection and placement on the basis of stimulation field modelling | |
EP4343707A2 (en) | Indication-dependent display of a medical image | |
US20190012805A1 (en) | Automatic detection of an artifact in patient image data | |
JP7561819B2 (ja) | 人体部分の撮像方法、コンピュータ、コンピュータ読み取り可能記憶媒体、コンピュータプログラム、および医療システム | |
US20180315204A1 (en) | Determination of Registration Accuracy | |
WO2024002476A1 (en) | Determining electrode orientation using optimized imaging parameters | |
WO2023110116A1 (en) | Ct-less free breathing image contour for planning radiation treatment | |
EP3794550B1 (en) | Comparison of a region of interest along a time series of images | |
US10832423B1 (en) | Optimizing an atlas |
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 |