CN104145295B - 使用不同指示器来确定放射治疗目标的位置变化 - Google Patents
使用不同指示器来确定放射治疗目标的位置变化 Download PDFInfo
- Publication number
- CN104145295B CN104145295B CN201280071164.2A CN201280071164A CN104145295B CN 104145295 B CN104145295 B CN 104145295B CN 201280071164 A CN201280071164 A CN 201280071164A CN 104145295 B CN104145295 B CN 104145295B
- Authority
- CN
- China
- Prior art keywords
- change
- indicator
- data
- target
- location
- 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
- 230000008859 change Effects 0.000 title claims abstract description 231
- 238000001959 radiotherapy Methods 0.000 title claims abstract description 12
- 238000000034 method Methods 0.000 claims abstract description 41
- 238000013507 mapping Methods 0.000 claims abstract description 38
- 150000001875 compounds Chemical class 0.000 claims abstract description 6
- 238000010276 construction Methods 0.000 claims description 140
- 230000033001 locomotion Effects 0.000 claims description 38
- 230000001225 therapeutic effect Effects 0.000 claims description 11
- 230000004044 response Effects 0.000 claims description 4
- 238000003384 imaging method Methods 0.000 description 26
- 230000005855 radiation Effects 0.000 description 21
- 238000005259 measurement Methods 0.000 description 17
- 206010028980 Neoplasm Diseases 0.000 description 15
- 238000004458 analytical method Methods 0.000 description 15
- 239000013598 vector Substances 0.000 description 14
- 230000006870 function Effects 0.000 description 12
- 238000012545 processing Methods 0.000 description 12
- 238000013500 data storage Methods 0.000 description 11
- 238000003672 processing method Methods 0.000 description 11
- 230000004888 barrier function Effects 0.000 description 10
- 210000000038 chest Anatomy 0.000 description 7
- 238000002591 computed tomography Methods 0.000 description 7
- 238000004590 computer program Methods 0.000 description 7
- 238000007408 cone-beam computed tomography Methods 0.000 description 7
- 210000004072 lung Anatomy 0.000 description 7
- 230000029058 respiratory gaseous exchange Effects 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- 238000002595 magnetic resonance imaging Methods 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 6
- 210000003205 muscle Anatomy 0.000 description 6
- 238000001514 detection method Methods 0.000 description 5
- 230000004927 fusion Effects 0.000 description 5
- 238000002372 labelling Methods 0.000 description 5
- 210000000056 organ Anatomy 0.000 description 5
- 239000002245 particle Substances 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000002224 dissection Methods 0.000 description 4
- 230000005670 electromagnetic radiation Effects 0.000 description 3
- 230000002349 favourable effect Effects 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 238000001356 surgical procedure Methods 0.000 description 3
- 208000003174 Brain Neoplasms Diseases 0.000 description 2
- 206010058467 Lung neoplasm malignant Diseases 0.000 description 2
- 208000003443 Unconsciousness Diseases 0.000 description 2
- 230000003187 abdominal effect Effects 0.000 description 2
- 210000003484 anatomy Anatomy 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 230000000968 intestinal effect Effects 0.000 description 2
- 210000000936 intestine Anatomy 0.000 description 2
- 210000004185 liver Anatomy 0.000 description 2
- 201000005202 lung cancer Diseases 0.000 description 2
- 208000020816 lung neoplasm Diseases 0.000 description 2
- 238000010801 machine learning Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 210000002307 prostate Anatomy 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 239000013589 supplement Substances 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 238000012935 Averaging Methods 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 206010027336 Menstruation delayed Diseases 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 210000000133 brain stem Anatomy 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000013329 compounding Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 238000004141 dimensional analysis Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000002675 image-guided surgery Methods 0.000 description 1
- 210000000876 intercostal muscle Anatomy 0.000 description 1
- 210000005075 mammary gland Anatomy 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 235000013372 meat Nutrition 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 230000000474 nursing effect Effects 0.000 description 1
- 238000011275 oncology therapy Methods 0.000 description 1
- 230000000399 orthopedic effect Effects 0.000 description 1
- 210000001672 ovary Anatomy 0.000 description 1
- 210000000496 pancreas Anatomy 0.000 description 1
- 230000036285 pathological change Effects 0.000 description 1
- 231100000915 pathological change Toxicity 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 210000000989 pectoralis minor Anatomy 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000035479 physiological effects, processes and functions Effects 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 210000001562 sternum Anatomy 0.000 description 1
- 238000012916 structural analysis Methods 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
- 210000000115 thoracic cavity Anatomy 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000004614 tumor growth Effects 0.000 description 1
- 210000004291 uterus Anatomy 0.000 description 1
- 210000001835 viscera Anatomy 0.000 description 1
- 230000000007 visual effect Effects 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/20—Analysis of motion
- G06T7/246—Analysis of motion using feature-based methods, e.g. the tracking of corners or segments
-
- 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/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1048—Monitoring, verifying, controlling systems and methods
- A61N5/1049—Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam
-
- 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
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
-
- 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
- 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
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1048—Monitoring, verifying, controlling systems and methods
- A61N5/1064—Monitoring, verifying, controlling systems and methods for adjusting radiation treatment in response to monitoring
-
- 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/30096—Tumor; Lesion
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/24—Aligning, centring, orientation detection or correction of the image
- G06V10/245—Aligning, centring, orientation detection or correction of the image by locating a pattern; Special marks for positioning
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biomedical Technology (AREA)
- Public Health (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Veterinary Medicine (AREA)
- Animal Behavior & Ethology (AREA)
- Radiology & Medical Imaging (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Optics & Photonics (AREA)
- Theoretical Computer Science (AREA)
- Physiology (AREA)
- Multimedia (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- Radiation-Therapy Devices (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明涉及一种用于确定称为目标变化数据并能用于执行放射治疗的数据的方法,所述目标变化数据描述了包含在患者体内的目标的位置变化的信息,其中该方法通过计算机来执行并且包括以下步骤:·获取描述映射规则的确定规则数据,其中所述映射规则用于将指示器变化集合映射到所述目标变化数据上,所述指示器变化集合是包括多于一个元素的集合,所述指示器变化集合的元素分别表示指示器位置的各个变化的信息,至少两个指示器分别指示称为指示器结构的不同身体结构的位置变化;·获取描述所获取的指示器变化集合的集合数据;以及·通过将映射规则应用于所获取的指示器变化集合上,基于确定规则数据和集合数据来确定目标变化数据。
Description
技术领域
本发明涉及放射治疗技术领域,具体说涉及确定待治疗的身体部位(称为目标)的位置(也由此称为“治疗患体部位”)的领域。指示器用以指示目标位置的变化,具体说以确定目标的实际位置。指示器的位置变化与目标的位置变化相关联。指示器可以是可检测的身体部位(称为“指示性身体部位”)或标记装置(如处于身体内部或附着于身体外部),具体说是单个标记物。
背景技术
可以参考专利申请EP 08169422.6和EP 09160153.4及其对应的美国专利申请US12/621,881和US 12/622,002以及WO 2011/107145。这些申请的公开内容通过引用并入本文。此外,还涉及Tong Lin等人于2008年12月11日在美国新泽西州Piscataway由IEEE举办的2008年MACHINE LEARNING AND APPLICATIONS(ICMLA’08)第七届国际会议上发表的“Tumor Targeting for Lung Cancer Radiotherapy Using Machine LearningTechniques”,其刊载于ISBN:978-0-7695-4395-4,XP31379459,第533–538页。
发明内容
本发明涉及医学领域,具体涉及为了治疗患体部位的辐射利用。
本发明涉及控制治疗束领域。治疗束治疗待治疗的身体部位,该身体部位在下文中还称为“目标”或“治疗患体部位”。这些身体部位是患者身体的特定部位,即解剖学上的患体部位。
本发明涉及医学领域,尤其涉及使用束具体说幅射束(也称为治疗束)来治疗患者的身体部位。电离辐射特别地用于治疗目的。具体说,治疗束包括或由电离辐射构成。电离辐射包括或由粒子(例如,亚原子粒子或离子)或电磁波构成,所述粒子或电磁波的能量足以使电子从原子或分子中分离,由此对其进行电离。这种电离辐射的示例是X射线、高能粒子(高能粒子束)和/或从放射性元件发出的电离辐射。治疗辐射,具体来说是治疗束,具体说用于放射治疗或放疗中,特别是用于肿瘤学领域的放射治疗或放疗中。特别是针对癌症的治疗,使用了电离辐射对带有诸如肿瘤的病理结构或组织的身体部位进行治疗。由此,肿瘤作为治疗患体部位的示例。
优选地,对治疗束进行控制使其能穿过治疗患体部位。
治疗患体部位的移动具体是由于下文中称为“关键移动”的动作所导致。关于这方面还可参考上述申请EP08169422.6和EP09160153.4,在这些申请中对这些关键移动进行了详细讨论。为了确定治疗患体部位的位置,使用诸如X射线装置、CT装置或MRT装置的分析装置来产生身体的分析图像(例如X射线图像或MRT图像)。具体说,分析装置采用成像方法,具体说分析装置是例如通过使用波和/或辐射和/或能量束(具体为电磁波和/或辐射、超声波和/或粒子束)来分析患者身体的装置。分析装置具体说是通过分析患者身体来产生患者身体(特别是患者身体的内部结构和/或解剖部位)的图像(例如,二维或三维图像)的装置。分析装置特别地用于医学诊断中,特别是放射学中。然而,可能难以在分析图像中识别出治疗患体部位。具体说,可能更容易识别出与治疗患体部位的位置变化、特别是治疗患体部位的移动变化相关联的指示性身体部位。因此,在已知指示性身体部位(也称为“替代物(surrogate)”)与治疗患体部位的位置变化(特别是移动变化)之间的相关性的基础上,跟踪指示性身体部位允许对治疗患体部位的移动进行跟踪。作为指示性身体部位的跟踪的备选方式或附加方式,标记装置(其可用作指示器,并称为“标记指示器”)可通过使用标记检测装置来进行跟踪。标记指示器(例如代理标记(fudiciary marker))的位置具体说与会由于关键移动而改变位置的指示器结构(诸如胸壁(如真肋或假肋)、隔膜或肠壁等)的位置具有已知的(预定的)相关性,具体说是,相对于该指示器结构的位置而言具有固定的相对位置。
在医学领域,成像方法用于产生人体解剖结构(例如,软组织、骨骼、内脏等)的图像数据(例如,二维或三维图像数据)。医学成像方法可理解为是指有利地基于仪器(apparatus-based)的成像方法(即所谓的医学成像模式和/或辐射成像方法),诸如计算机层析成像(CT)和锥束计算机层析成像(CBCT,特别是测定体积式CBCT)、X射线层析成像、磁共振成像(MRT或MRI)、常规X线、超声波和/或超声检查以及正电子发射层析成像。分析装置特别地用于以基于仪器的成像方法来产生图像数据。成像方法特别用于医学诊断,以分析解剖体从而产生可由图像数据来描述的图像。成像方法特别用于检测人体中的病理变化。肿瘤例如代表了解剖结构变化的示例。如果肿瘤生长,那么可以说成是代表了扩张的解剖结构。例如,当使用造影剂渗入肿瘤时,早期/晚期脑部肿瘤通常在MRI扫描中可见。MRI扫描代表了成像方法的示例。在脑肿瘤的MRI扫描情况下,认为MRI图像中的信号增强(由于造影剂渗入肿瘤而导致)代表了固体肿瘤块。由此,肿瘤在由成像方法产生的图像中是可检测的,特别是可辨别的。
待由标记检测装置(例如相机、超声接收器或诸如CT或MRI的分析装置)检测的标记物的功能是使它的空间位置(即,其空间位置和/或排列)能够被确定。检测装置具体地根据导航系统(IGS,图像引导手术)而部分已知。标记物可以是主动式标记物。主动式标记物可例如发射电磁辐射和/或波,其中所述辐射可位于红外、可见光和/或紫外光谱范围内。然而,标记物也可以是被动式的,即例如能反射位于红外线、可见光和/或紫外光谱范围内的电磁辐射、或者能阻断X射线辐射。为此,标记物可设有具有相应反射特性的表面,也可以由金属制成以阻断X射线辐射。标记物也可以在无线电频率范围内或以超声波的波长来反射和/或发射电磁辐射和/或波。标记物优选具有球形和/或球体形状,并因此可称为标记球体;然而,标记物还可展现为例如立方体的带角形状。
标记装置可以是例如参考星标或指针,或者是一个标记物,或者优选成预定的空间关系一个以上的(单个)标记物。标记装置包括一个、两个、三个或在两个或多个标记物成预定的空间关系的情况下的多个标记物。这种预定的空间关系具体通过数据处理方法获得,并且例如存储在治疗束系统的计算机中。
本发明的目的在于确定目标的位置变化。
该目的通过独立权利要求的保护主题来解决。从属权利要求涉及优选的实施方式。本发明的优点、有利特征、优选实施例以及有利方面在下文中公开,并包含在从属权利要求的保护主题中。只要技术上合理并可行,则不同的有利特征可根据本发明进行组合。特别地,与一个实施例的某特征具有相同或相似的功能的另一实施例的另一特征之间可以互相交换。特别地,能为一个实施例补充又一功能的另一实施例的某特性可添加于那一实施例中。
本发明的发明人已经认识到,目标的关键移动是一个以上的执行关键移动的身体结构(本文中称为指示器结构)所导致的结果,并且这些指示器结构的移动与目标之间存在复杂的关系。本发明人得到的结论是,为了分别建立目标与各自的指示器结构之间的独立相关性,单独处理指示器的做法是不恰当的。例如,肺部的位置和几何构造(形状和/或尺寸)不仅受到一个指示器结构(例如隔膜)的影响,而且还受到其他指示器结构(例如真肋和假肋)的影响。隔膜执行颅尾移动,真肋执行前后移动,而假肋则执行横向移动。然而,真肋的移动例如也具有颅尾方向上的分量,由此该移动会干扰隔膜的颅尾移动。而且,由于指示器结构的不同移动所导致的肺部变形又与不同方向上的移动耦合。因此,根据本发明,优选考虑若干不同指示器结构的位置变化,以确定目标的位置(和/或几何构造)的变化。目标的位置(和/或几何构造)变化可具体地基于对嵌有目标的结构(称为“目标结构”)的位置变化和/或几何构造变化的确定来确定。嵌有肿瘤的目标结构以及任何类型的身体结构的示例可以例如是像肺、肝、乳腺、胰腺、肠、子宫、卵巢和前列腺这样的器官。指示器结构的移动还可经由中间身体结构(称为“耦接结构”)传递至目标结构。例如,尽管存在中间身体结构,但前列腺中的目标位置(即,目标结构)还是受到隔膜(作为第一指示器结构的示例)和肠壁(作为第二指示器结构的示例)的影响。同时,这些耦接结构导致了一组不同指示器结构的位置变化与目标的位置变化之间的复杂关系。
指示器结构具体来说是那些通过由脑干进行特定的无意识控制的肌肉(如隔膜)来表示其位置的结构。在本发明的意义上,指示器结构也可以是直接联结(具体说附着)于肌肉的身体结构。具体地,指示器结构的移动直接联结于肌肉(具体说是被无意识控制的肌肉)的移动。例如,外部肋间肌抬升胸廓的顶边界。还例如,胸锁乳突肌和胸小肌可影响胸腔和胸骨的位置。
只要下文中涉及到对象的位置变化的确定,那么这就同时包括在执行位置变化之后的该对象的新位置的确定。新位置的确定可基于先前的位置(在位置变化之前)和所确定的位置变化来执行。具体地,本发明的数据处理方法获取不同对象(例如目标、指示器、指示器结构等)的起始位置,并使用本发明的方法从已知的起始位置开始跟踪对象的变化。就确定目标的起始位置而言,例如可使用上述专利申请EP08169422.6、EP09160153.4或者WO2011/107145中所描述的任何方法。具体地,可以使用4D成像方法(具体说是四维CBCT,即4D CBCT)来确定目标和/或指示器结构的起始位置(具体在放射治疗之前确定)。本发明的方法获取起始位置,以便在放射治疗过程中跟踪目标。具体地,指示器可被跟踪,并且可基于执行对患者的身体的4D成像方法(例如,4D CT或4D MR)所得到的测量而获取所要求的确定规则数据,以确定目标位置,具体说用以确定指示性身体部位的位置。优选地,针对获取所述判断规则数据,作为对指示性患体部分位置的确定的备选或附加方式,还通过标记检测器装置来获取标记物指示器的位置,从而建立指示器(即指示性身体部位和/或指示性标记物)与目标之间的已知相对位置。通过这种方法,优选地确定出描述了指示器(标记物指示器和/或指示性身体部位)与目标之间的相对位置的位置集合,该位置集合可用作根据本发明的用于跟踪的起始点。
更详细地,本发明涉及一种用于确定目标的位置变化以便为患者身体执行放射治疗的方法(具体为跟踪方法)。该方法具体在治疗期间执行,但也可以在治疗(特别是启动治疗)之前执行,以特别地用于训练目的或者用于将本发明的跟踪方法锁定在目标上(并且随后再在治疗期间继续跟踪)。锁定在目标上的做法具体可通过四维成像方法(如4D CT或4DMR)来协助进行。
该方法优选包括获取确定规则数据的步骤。所述确定规则数据描述了指示器的位置变化与目标的位置变化之间的关系,具体说,指示器位置与目标位置之间的关系。所述确定规则数据具体说通过将(多于一个的)不同位置变化映射于目标的(一个)位置变化上的方式来描述一组不同指示器(多于一个的不同指示器)的位置变化与(单个)目标的位置变化之间的相关性。确定规则数据可基于使用分析装置和/或标记检测器而产生的测量数据来确定。特别地,使用四维成像方法,从而依据时间来确定目标的位置,并优选地通过分析装置和/或标记检测器来同时确定对应的指示器位置。优选地使用分析装置和/或标记检测器,以便将多个可能的目标结构位置与多个(多于一个的)不同的指示器位置(即,一组位置)(其特别是针对指示器位置的多个不同集合而言)相互关联。如下面将更详细描述的那样,所述确定规则数据可基于从一个患者或多个患者所产生的的测量数据来确定。所述确定规则数据具体描述了用于根据不同指示器(其与不同的指示器结构联结)的位置变化而导出目标的位置变化的规则。指示器的位置变化表示指示器变化集合,并通过称为集合数据的数据来描述。
优选地是,获取描述了称为“指示器变化集合”的集合的集合数据。指示器变化集合包括多于一个的元素,并且每个元素表示关于指示器(例如,标记物指示器或指示性身体部位)的位置变化的信息。也就是,根据本发明所使用的不仅仅是一个指示器的位置变化,而是一组不同的信息(物理信息),具体来说是可由具体代表值或向量的n元组合来描述的参数,特别是能具体描述指示器位置或指示器的位置变化的物理参数。优选地,所述指示器变化集合包括能描述关于不同指示器结构的位置信息的元素。“不同”在本文中意指:指示器结构由不同的肌肉来移动或者表示不同的肌肉(其具体由于身体的重要功能而移动)。特别地,“不同”是指指示器结构以不同的方向移动,具体是说当经受关键移动时特别地遵循不同的轨迹。
映射规则将指示器变化集合映射到一个特定目标变化数据上,所述特定目标变化数据描述了该(一个)目标的位置的一个变化以及特别是变化之后的目标位置。也就是,将能表示不同的指示器结构的所述变化集合的元素映射到相同的目标变化数据上,具体说是映射到相同的目标位置变化上。指示器结构的位置变化在本文中也称为“指示器结构变化”。映射规则取决于由指示器变化集合所描述的位置变化。针对变化的规则还取决于变化之前、变化期间、变化之后的标识位置。
上述确定规则数据优选地基于上述对患者和/或患者池进行的测量来建立。所获取的确定规则数据描述了用于将所获取的指示器变化集合映射到目标变化上的规则。也就是,例如,指示器的位置变化集合可用来通过应用由确定规则数据所描述的规则来确定目标的变化。所描述的规则优选为用于将指示器变化集合映射到目标变化上的规则。具体地,所述规则描述了用于将多个指示器变化集合中的各个集合映射到多个目标变化中的各个变化上的规则。用于映射的规则导致了:将该规则应用于所获取的指示器变化集合的结果是确定出一个目标变化。换言之,虽然指示器变化集合中的两个或多个元素表示不同的指示器结构的变化,但却将应用映射规则的指示器变化集合的元素都映射到相同的目标变化上。具体地说,映射规则可由满射函数来表示,具体说由将多个指示器变化集合中的一项转换成多个目标变化中的一项的单射函数来表示。“多个”可以是多个不同的指示器变化集合和/或目标变化,或者也可以是指示器变化集合和/或目标变化的(连续的)多项。所述规则具体说构造成适于所获取的某范围的指示器变化集合,和/或适于相应范围的目标变化。所述范围具体限定为例如对目标位置变化的覆盖小于5cm。
根据该方法的另一步骤,通过将由确定规则数据所描述的规则应用于所获取的指示器变化集合来基于所述确定规则数据和集合数据来确定目标变化。
不同的指示器结构可引起目标的不同移动。例如,如果肺部为目标结构,那么真肋的关键移动可引起目标的前后移动,假肋的关键移动可引起目标的横向移动,而隔膜可引起目标的颅尾移动。因此,目标可在特定方向上进行移动,其中可以这样理解特定方向:如同可通过将指向不同的方向的若干向量相加来描述向量那样,该特定方向可由不同的方向分量组成。根据本发明,这些相加向量中的一个向量特别地与表示不同的指示器结构的变化集合中的多于一个的元素相关联。根据本发明的一方面,映射规则构造成将目标的变化方向与指示器变化集合中多于一个的元素相关联。也就是,变化集合中的两个或多个元素可影响所述目标在目标变化的方向上的变化。假设变化集合的元素由例如表示位置变化向量或沿轨迹的往复移动的幅值的值和/或向量来表示的话,那么根据映射规则,允许并具体描述了(至少针对目标变化的子范围而言):使用表示关于两个或多个不同指示器结构的位置变化的信息(具体说表示位置变化)的两个或多个不同元素的两个或多个值和/或向量来确定沿主要由不同的指示器结构中的一个指示器结构的移动所导致(具体说是主要与不同的指示器结构中的一个指示器结构的移动相联结)的目标变化的那一方向上的目标变化。
具体地,联结两个或多个元素的数学关系(由值和/或向量所述描述的)用于描述所述映射规则。例如,可使用沿着颅尾方向的真肋移动分量与沿着颅尾方向的隔膜移动分量之间的数学关系来确定目标变化。数学关系具体通过如加、减、乘、除的数学运算符来描述。所述数学关系可通过数学运算符的组合来进行说明。例如,n元组(具体为向量)由指示器变化集合来描述并表示不同指示器结构的位置变化,而是n元组,具体为描述了目标的位置变化的向量,由此,映射规则可表示为,例如:
其中,是矩阵,该矩阵的元素可以是目标位置变化之前和/或之后的目标位置的函数,和/或是指示器Ii的位置(在指示器位置变化之前和/或之后)的函数。具体地,所述矩阵的元素可以是本文所提及的关于指示器位置的信息的任何参数的函数,这些信息具体说是一个或多个指示器的往复移动的频率和/或振幅以及下文提及的关于补充信息的参数。n元组描述了目标的位置变化的信息,并可包括仅描述了目标的位置变化的元素,或者可包括例如既描述了目标的位置变化又描述了目标位置的元素。n元组可包括仅描述了指示器的位置变化的元素,或者是既描述了指示器的位置变化又描述了指示器位置的元素。n元组描述了关于指示器的位置变化的信息。优选地,如上所述,映射规则考虑了复杂影响,具体说是指示器结构的移动对目标移动的复杂干扰影响。因此,矩阵优选构造为使得不存在其中矩阵是对角矩阵的参考系统。作为数学关系的又一示例,目标的位置变化量或者的一个或多个分量例如Δx可以是标识变化集合的一个或多个元素Δs的函数,例如
或者
特别地其中k是常数,并且Δsi和Δsj是描述了指示器Ii和另一不同指示器Ij(其遵循不同的轨迹并且表示不同指示器结构的移动)的位置变化量的指示器变化集合的元素。
关于指示器位置变化的信息(可由n元组描述)具体包括如下信息中的至少一项:指示器的位置变化、指示器所遵循的轨迹内的位置变化、指示器移动的方向变化以及指示器往复移动的振幅变化。上述信息也称为位置信息。备选地,关于指示器位置变化的信息(由n元组描述)还可包括与时间有关的信息,诸如在通过指示器变化集合来表示的指示器位置变化的期间所经过的时间。
具体地,关于指示器位置变化的信息(由n元组描述)可包括关于位置变化之前或之后的指示器位置的信息,具体是关于变化之前或之后所给出的指示器往复移动的振幅的信息。具体地,关于变化的信息可包括往复移动的当前频率,具体说是往复移动在变化时刻的当前最大振幅。特别地,目标变化根据具体在各个轨迹中的给定的指示器位置而不同,或者根据当指示器发生位置变化时给定的目标位置而不同。具体地,根据不同的位置指示器变化集合中元素的相同量的位置变化对目标变化具有不同影响。具体地,取决于元素所表示的指示器是否接近最大振幅或位于两个最大往复移动中间。还取决于指示器的最大振幅。例如,表示胸部呼吸的一个指示器的相同振幅变化可根据胸部呼吸是否是由腹部呼吸来补充而导致不同的目标变化。
根据本发明的又一实施例,能被确定规则数据所考量的、还能通过关于指示器位置变化的信息来描述的附加信息(称为“补偿信息”)是:指示器沿轨迹的往复移动的最大振幅、往复移动的频率、目标及指示器的医学信息(诸如目标类型(例如肿瘤类型))、目标结构内的目标位置、由指示器来表示的指示器结构的类型(具体说是指示器结构的主要移动方向),以及指示器的医学分类(诸如指示器结构(诸如假肋或隔膜)的医学名称,或者肌肉结构、假肋、真肋或涉及关键移动的器官(诸如心脏或肠)的名称)。
根据本发明的实施例,映射规则基于描述了指示器变化集合与相应的目标变化之间的对应关系的对应数据。对应数据可实现为一个或多个查找表。具体地,如果所获取的指示器变化集合与查找表中条目中的一项不同,则在查找表的条目之间进行插补。所述对应数据可具体地通过测量针对指示器往复移动的不同的最大振幅和/或频率的指示器的位置变化(称为“指示器变化”)以及通过测量针对沿往复轨迹的移动中的多个位置的指示器位置来产生。具体地,映射规则可包括针对指示器往复移动的不同最大振幅集合和/或不同频率的不同查找表。选择对应于由指示器变化集合所描述的最大振幅集合或频率的那个查找表。如果没有准确对应关系,那么在查找表之间进行插补。
指示器变化的测量优选地结合目标位置测量来执行,从而将每个目标位置(并且优选为变化之前的目标位置)分配给一组指示器位置,并由此使得每个目标位置变化能分配给指示器变化集合(所述指示器变化集合具体除了包括指示器位置变化的信息以外,还包括关于变化之前的指示器位置的信息)。对应数据可根据个体患者、具体通过使用上述四维成像方法和使用上述用于检测目标位置并任选地用于检测指示器位置的分析装置来产生。针对检测标记物的位置和位置变化,可以使用指示器和标记检测器。上述能描述指示器变化集合与目标变化之间的对应关系的对应数据称为“指示器目标对应数据”。根据该数据是否是基于对携带待治疗目标的个体患者所执行的测量,将该数据称为“个体指示器目标对应数据”,而如果数据是基于具有类似目标的多个不同患者(患者池),那么该数据则称为“通用目标对应数据”。“类似目标”具体是指,目标属于相同的医学类型(例如,相同类型的肿瘤),具体是说,目标在相同的目标结构中的具有相同的相对位置(例如,当肺部是目标结构时,在肺部中的相同位置),特别是患者具有大约相同的医学条件(例如,相同的年龄、相同的呼吸模式等)的情况以及特别是目标是否从一个(第一)身体结构向下一个(第二)身体结构延伸的情况。在第一种情况下,目标结构由第一身体结构构成。在后者的情况下,目标结构由所述第一和第二身体结构构成,并且当第一和第二身体结构由于关键移动而进行不同的移动时,所述延伸会对目标的位置变化造成影响。例如,肺癌会从肺部(第一身体结构)向胸壁(第二身体结构)生长。可采用类似措施和/或类似阈值来确定目标是否相似。
根据又一实施例,目标结构的位置的变化和/或几何构造的变化是在中间步骤中基于指示器变化集合来确定的,然后,在下一步骤中,基于目标结构的位置的变化和/或几何构造的变化来确定目标的位置变化(目标变化)。目标结构的几何构造变化可对目标的位置造成影响,因为几何构造的变化可导致目标结构的某些部分收缩而目标结构的其他部分扩张,这会导致目标的移动。这种几何构造对位置的影响例如可通过使用目标结构的模型(例如,目标结构的粘弹性模型或目标结构的弹性模型)来确定。所述模型具体说模仿了目标结构的内部结构对目标结构的边界变形和/或位移的机械响应。所述模型具体描述了目标结构(其内部结构可例如具有粘弹性性质)对目标结构边界上的变化的机械响应,其中所述变化与指示器变化集合相关联。所述模型具体描述了目标结构的内部结构,并且特别地可由目标结构的应力-应变曲线、目标结构的应力松弛特性和蠕变特性来描述。关于目标结构的几何构造(具体说是边界)的数据可通过由如上所述的分析装置所执行的成像方法来获取。关于如何模仿器官的方法例如记载在James C Ross等人发表于2008的《MedicalImage Computing and Computer-Assisted Intervention》(MICCAI2008)中的“Real-TimeLiver Motion Compensation for MrgFUS”;2008年9月6日Springer Berlin于海德堡出版的《计算机科学讲义》,ISBN:978-3-540-85989-5,XP19105238,总页数806-813,Bd.5242;以及Pierre-Frederik Villard,Wesley Bourne和Fernando Bello发表的“Modelling OrganDeformation using Mass-springs and Tensional Integrity”。这些论文的公开内容通过引用并入本文。
上述模型可具体包括关于弹性模块的假设。所述模型可具体基于个体患者或具有相似目标的多个患者的位置变化和目标结构的几何构造变化的多个测量来生成。
所述模型具体取决于目标结构的类型。目标结构的类型具体还取决于身体结构的功能类型、器官的类型、身体结构的密度、身体内身体结构的位置、身体结构的变形、身体结构的弹性(刚硬或不刚硬)等。
根据上述实施例,优选的是,目标结构中的目标位置通过获取描述了该位置的、称为“目标位置数据”的数据而获取。当本发明的跟踪方法开始时,所述目标位置数据优选地基于上述起始位置,并能被获取以使得所述数据与在考虑当前的指示器位置变化之前由该跟踪方法所确定的目标位置相对应。而且,优选地,体内的目标结构的几何构造和/或位置由目标几何数据来描述。优选地,映射规则考虑了目标位置数据和/或目标几何数据。
根据本实施例,映射规则进一步基于描述了指示器变化集合与目标结构的相应变化(即,位置和/或几何构造的变化)之间的对应关系的对应数据。目标结构的变化可例如通过目标结构的边界的位置和/或几何构造(形状和/或尺寸)的变化来描述。所述对应数据称为目标结构对应数据。如果所述目标结构对应数据基于对包括有待治疗目标的个体患者的测量,那么该对应数据称为“个体目标结构对应数据”。个体目标结构对应数据可具体通过测量针对不同的指示器变化集合的目标结构的位置和/或几何构造的方式产生,其中所述不同的指示器变化集合可通过引起不同振幅和/或频率的呼吸和/或心跳来引起,以及通过引起诸如腹式呼吸或胸部呼吸的不同类型的呼吸来引起。
作为个体目标结构对应数据的备选或附加,映射规则也可以基于通用目标结构对应数据,其描述了指示器变化集合与通过测量多个不同患者的对应关系所产生的相应的目标结构变化之间的对应关系。如上所述,这些患者优选具有相同类型的目标(具体说是目标结构内的相同位置、相同目标结构以及相同目标类型(例如相同的肿瘤类型、具体说是相同的肿瘤分类))。
通过使用目标结构模型,可基于所获取的指示器变化集合和(个体和/或通用的)目标结构对应数据来确定目标变化。具体地说,基于目标结构对应数据和指示器变化集合来确定目标结构变化。然后,基于所确定的目标结构变化、目标结构模型和任选的目标位置数据来确定目标变化。
通常,根据目标结构的位置和/或几何构造的变化(称为“目标结构变化”)的目标变化确定可通过将目标结构变化映射到目标变化上的、称为“子规则”的规则来规定。目标结构变化不仅可描述目标结构的位置变化和/或几何构造变化,而且还可描述目标结构的位置和/或几何构造。
所述子规则可基于称为目标结构模型的上述模型。目标结构模型优选描述了目标结构的物理参数,该参数规定了由于目标结构变化而导致的目标结构的内部结构的变化。例如,目标结构可通过包含相联结的节点的网格来进行建模。这种联结可例如表现为弹性或粘弹性形式。所述目标结构模型具体模仿了由于内部结构变化所导致的目标结构内的目标位置(具体说是几何构造)变化。优选根据目标位置数据进行目标变化的建模,例如根据上述网格内的目标位置或者也可根据目标几何数据来进行建模。
根据另一实施例,子规则基于弹性融合。具体地,可将表示了目标结构变化之前的目标结构的目标结构模型作为起始点。该模型例如可以是由上述一种成像方法产生的二维或三维的目标结构图像。因此,弹性融合具体为图像融合。可基于目标结构对应数据和指示器变化集合来确定目标结构的边界的位置变化和/或几何构造变化。在确定改变的边界之后,目标结构模型可进行弹性融合,以适应改变的边界。因此,目标结构变化之前的目标结构模型可融合到目标结构变化之后的目标结构中。目标结构模型优选地包括目标。因此,通过融合目标结构的模型以适应改变的边界,使得目标结构中的目标位置同样改变。因此,可通过进行弹性融合来确定目标变化。为此,优选地是,使用最好是描述了目标结构变化之前的目标结构模型中的目标位置的目标位置数据。
根据另一实施例,子规则基于描述了多个不同的目标结构变化与相应的目标变化之间的对应关系的对应数据。如上所述,这种数据可通过对个体患者执行测量以得到目标结构的多个不同的位置和/或几何构造来产生,具体说是通过使用由分析装置执行的成像方法来产生。这些测量可针对携带有待治疗目标的个体患者来执行,和/或针对满足了上述目标类型的相似性条件的患者池中的多个患者来执行。
当然,可以结合上述子规则来确定目标变化。上述目标变化表示描述了关于目标的位置变化的信息的数据。此信息具体描述目标从容纳了变化前的目标的先前位置向容纳了变化后的目标的随后位置的位置变化。称为目标变化数据的数据还可描述先前和/或随后位置。具体地,根据一实施例,所述目标变化数据描述了在指示器位置变化之后所容纳的目标的位置。也就是,位置的变化通过描述先前位置和随后位置来说明。而且,在本发明跟踪方法的下一步骤中,目标变化数据可用作目标位置数据。相应地,指示器变化集合可描述在变化之前和变化之后分别由指示器所占据的先前位置和随后位置。映射规则具体可将指示器的先前位置映射到目标的先前位置上,并将指示器的随后位置映射到目标的随后位置上。换言之,本发明不仅仅覆盖了位置的相对变化,而且还构造为能基于指示器变化集合、映射规则以及前一步骤(在考虑当前变化之前)中所确定的目标变化数据来确定目标的绝对位置。
指示器和/或目标的位置可以在参考系统中描述。该参考系统例如是用于进行放射治疗的治疗装置所处的参考系统,或者是患者的患体部位(例如目标、目标结构、指示器结构或像脊椎这样的其他身体部位)所处的参考系统,或者例如是用于执行放射治疗的房间所处的参考系统。在这一方面,还可参考上面已提及的申请EP 2189940、EP 2189943和WO2011/107145。
本发明还涉及一种程序,所述程序当运行于计算机上或加载到计算机上时使计算机执行本文所述的一个、多个或所有的方法步骤;和/或涉及一种存储有(具体是以非暂时性形式)所述程序的程序存储介质;和/或涉及一种运行所述程序的计算机或其存储器加载有所述程序的计算机;和/或涉及一种携带有可表示程序(特别是上述程序)的信息的信号波,具体说是数字信号波形,所述程序具体说包含适于执行本文所述的任何或全部方法步骤的代码装置。
在本发明的框架内,计算机程序单元可通过硬件和/或软件(这包括固件、驻留软件、微代码等)来实施。在本发明的框架内,计算机程序单元可采取通过计算机可用、特别是计算机可读的数据存储介质来实施的计算机程序产品的形式,具体说可采取所述数据存储介质中包含的指令执行系统上采用的或与之结合使用的计算机可读的程序指令、“代码”或“计算机程序”的形式。这种系统可以是计算机;计算机可以是包括用于执行根据本发明的计算机程序单元和/或程序的装置的数据处理设备,特别是包括执行计算机程序单元的数字处理器(中央处理单元,即CPU)以及可选地包括用于存储执行计算机程序单元所使用的和/或产生的数据的易失性存储器(具体说是随机存取存储器,即RAM)的数据处理设备。在本发明的框架内,计算机可用、特别是计算机可读的数据存储介质可以是任何数据存储介质,其可以包含、存储、传送、传播或传输那些指令执行系统、装置或设备上使用或与之结合使用的程序。计算机可用、特别是计算机可读的数据存储介质例如可以是,但不限于电子、磁、光、电磁、红外或半导体系统、仪器或设备,或者是诸如互联网的传播介质。计算机可用或计算机可读的数据存储介质甚至可以是例如可打印所述程序的纸或其他合适介质,因为程序可以通过电的方式而获取,例如通过光学扫描该纸或其他合适介质,然后再编译、解码或以适当方式另行处理。优选地,数据存储介质是非易失数据存储介质。本文所述的计算机程序产品以及任何软件和/或硬件形成用于执行本发明的示例性实施例中的功能的各种装置。计算机和/或数据处理设备可特别地包括具有可输出引导信息的装置的引导信息设备。引导信息可例如通过视觉指示装置(如监视器和/或灯)可视地输出给用户,和/或通过听觉指示装置(如扬声器和/或数字语音输出装置)可听地输出给用户,和/或通过触觉指示装置(如振动元件或并入到仪器中的振动元件)可触摸地输出给用户。
根据本发明的方法具体而言是一种数据处理方法。所述数据处理方法优选使用具体说是计算机的技术手段来执行。特别地,所述数据处理方法由计算机执行或在计算机上执行。为了处理数据,特别是为了以电的方式和/或光的方式处理数据,该计算机具体说还包括处理器和存储器。所描述的计算步骤具体说通过计算机来执行。确定步骤或计算步骤具体说是在该技术性的数据处理方法框架内(特别是在程序框架内)对数据进行确定的步骤。计算机可以具体为任何种类的数据处理装置,尤其是电子数据处理装置。计算机可以是那些通常被认为是例如台式个人计算机PC、笔记本、上网本等的装置,但计算机也可以是任何可编程装置,例如移动电话或嵌入式处理器。计算机可特别地包括“子计算机”系统(网络),其中每个子计算机代表其本身的计算机。术语“计算机”包含云计算机,特别是云服务器。术语“云计算机”包含云计算机系统,特别地包括具有至少一个云计算机的系统,特别是具有多个可操作式互连的云计算机(例如服务器场)的系统。优选地,云计算机与诸如万维网(WWW)的广域网连接。这种云计算机位于连接于万维网的计算机的所谓“云”中。这种基础设施用于云计算,所述云计算描述了不要求终端用户获知那些传送特定服务的计算机的物理位置及配置的计算、软件、数据访问和存储服务。特别地,术语“云”用来比喻互联网(万维网)。特别地,云提供作为服务(LaaS)的计算基础设施。云计算机可以用作用于执行本发明方法的操作系统和/或数据处理应用程序的虚拟主机。优选地,云计算机是由Amazon WebServicesTM提供的弹性计算云(EC2)。为了接收或输出数据和/或执行模数转换,计算机特别地还包括接口。该数据特别地是可表示物理特性的数据和/或可从技术信号生成的数据。该技术信号特别地通过(技术性)检测装置(例如,用于检测标记装置的装置)和/或(技术性)分析装置(例如,用于执行成像方法的装置)产生,其中所述技术信号特别地是电或光信号。技术信号特别地表示由计算机接收或输出的数据。
优选地是,本发明的方法至少部分地通过计算机来执行。也就是说,本发明方法的所有步骤或仅某些步骤(即少于总步骤数量的步骤)可以由计算机来执行。
表述“获取数据”具体说(在数据处理方法的框架内)包含数据是由数据处理方法或程序确定的情况。确定数据具体说包括:测量物理量,以及将所测量的值转换成具体说是数字数据和/或通过计算机来计算该数据,特别是在本发明的方法范围内计算该数据。“获取数据”的含义还特别地包含这样的情况,即通过数据处理方法或程序从例如另一程序、前一方法步骤或数据存储介质中接收或检索数据,进而特别地用于通过数据处理方法或程序进行进一步处理。因此,“获取数据”还可例如表示等待接收数据和/或正在接收数据。所接收的数据可例如经由接口来输入。“获取数据”还可指代数据处理方法或程序为了从数据源(主动)或经由接口(例如,从另一计算机或网络)接收或检索数据所执行的步骤,所述数据源例如为数据存储介质(例如,ROM、RAM、数据库、硬盘等)。数据可通过在获取步骤之前执行附加步骤的方式来实现“就绪”状态。根据附加步骤来生成数据,以用于获取。对数据进行特定的探测或捕获(例如,通过分析装置)。作为备选或附加,数据是根据所述附加步骤例如通过接口来输入的。可在特定况下输入所生成的数据(例如,输入到计算机中)。根据附加步骤(其在获取步骤之前进行),数据也可以通过执行将数据存储于数据存储介质(例如,ROM、RAM、CD和/或硬盘驱动器)的附加步骤的方式来提供,从而使数据在根据本发明的方法或程序的框架内就绪。因此,“获取数据”还可涉及操控装置以获得和/或提供待获取的数据。特别地,获取步骤不涉及侵入性步骤,侵入性步骤表示对身体进行实质上的生理干扰,其需要对身体采取专业医疗措施,并且即使采取了所需的专业医疗措施或护理,该身体还是会承受极大的健康风险。特别地,获取数据具体说是确定数据不涉及外科手术步骤,并且尤其不涉及使用外科手术或外科疗法来治疗人体或动物体的步骤。这也特别地适用于涉及确定数据的任何步骤。为了区分本发明的方法所使用的不同数据,数据被标记为(即简称为)“XY数据”等,并由其所描述的、称为“XY信息”的信息进行定义。
本发明还涉及一种基于所确定的目标变化数据来控制治疗束的方法。
本发明还针对一种治疗束系统。所述治疗束系统包括治疗装置。所述治疗装置构造成可发射治疗束。而且,所述治疗束系统包括运行上述程序或其中加载有上述程序的计算机。由此,所述计算机执行本发明的方法。具体地,所述计算机基于所获取的指示器集合数据来确定目标变化数据。而且,所述治疗束系统包括构造成可发射治疗束、并进一步构造成可根据所述计算机确定的目标变化数据来控制治疗束位置的治疗装置。例如,所述计算机将目标变化数据传输至治疗装置。
本发明的附加特征和优点将在下面关于本发明实施例的详细描述中公开。
附图说明
图1示出了治疗束系统;
图2示出了指示器I1和12的轨迹、位置和位置变化;
图3示出了用于依据指示器位置来确定目标位置的查找表;
图4示出了由于目标结构的变形而导致的目标模型化移动。
具体实施方式
图1示出了根据本发明的治疗束系统。治疗装置100发射治疗束110。装置100可按照箭头所指示的那样移动,从而引导治疗束110。计算机200与处理装置100电连接,以便通过控制信号来控制对治疗束110的引导。屏幕201、键盘202和鼠标203连接到计算机200上。计算机200还控制分析装置310和320。分析装置310和320具体来说为拍摄患者400的x射线图像的x射线装置。计算机200与分析设备310和320连接,以接收来自分析装置310和320的预先图像数据(在治疗之前)和处理图象数据(在治疗过程中)。具体地,计算机200还配置为控制分析装置310和320所执行的图像采集。在分析装置310和320与计算机200之间优选具有电连接,以使计算机200向分析装置310和320发送控制信号和/或从分析装置310和320接收上述数据。患者400躺在手术台500上。分析装置310和320以及治疗装置100由支承结构600来支承。根据又一实施例,分析装置310和320与治疗装置分离,具体说是分开安装。为了控制治疗束110与患者400之间的相对位置,也可移动手术台500,从而相对于治疗束110对患者400重新定位。总之,治疗束系统优选配置为控制治疗束110的位置,具体是相对于患者400的位置。这具体是由计算机200来执行。此外,辅助控制装置可例如集成到治疗装置100和/或手术台500中,并最好与计算机200配合操作。
计算机200优选还具有关于分析装置310和320相对于治疗束110的位置的成像几何构造的信息。
关于成像几何构造的信息优选包括如下信息:即,如果待分析的对象(解剖体部分)是已知的,在给出成像几何构造分析装置与待由X射线辐射进行分析的分析对象(即解剖患体部位)之间的已知相对位置的情况下,允许对分析图像(X射线图像)进行计算的信息,其中“已知”是指该分析对象的空间几何构造(尺寸及形状)是已知的。这具体意味着关于分析对象(解剖患体部为)与分析辐射(X射线辐射)之间的交互的三维“空间解析”的信息是已知的,其中“交互”意指例如分析辐射被阻断或者部分或完全被允许穿过分析对象。成像几何构造的位置具体说是取向,具体由X射线装置的位置来定义,特别是由X射线源和X射线检测器的位置来定义,和/或特别是由穿过分析对象并由X射线检测器检测到的多种(多重)X射线束的位置来定义。成像几何构造具体描述了所述多种(多重)射线束的位置(具体说是取向)和形状(例如,呈现特定倾斜角度的圆锥形状)。所述位置可具体地由穿过所述多重射线束的中心的X射线束的位置来表示,或者由能代表X射线束多重体(多支体)的几何对象(诸如截头圆锥)的位置来表示。关于上述交互的信息优选地例如通过三维CT而在三维上已知,并且以空间解析的方式描述了针对分析对象的点和/或区域的交互,具体来说是分析对象的所有点和/或区域中的力。具体地,成像几何构造方面的知识允许相对于图像平面(例如,X射线检测器的平面)来计算辐射源(例如,X射线源)的位置。就三维分析对象与通过成像几何构造定义的二维分析图像之间的关联而言,具体参考以下出版物:
1.Roger Y.Tsai于1986年佛罗里达州迈阿密海滩举办的题为Computer Visionand Pattern Recognition的IEEE会议记录第364-374页刊载的“An Efficient andAccurate Camera Calibration Technique for 3D Machine Vision”。
2.Roger Y.Tsai于1987年8月的IEEE期刊《Robotics and Automation》第RA-3卷,第4号,第323-344页刊载的“A Versatile Camera Calibration Technique for High-Accuracy 3D Machine Vision Metrology Using Off-the-Shelf TV Cameras andLenses”。也可参见http://www.cs.cmu.edu/~rgw/TsaiDesc.html
3.Ziv Yaniv发表的“Fluoroscopic X-ray Image Processing andRegistration for Computer-Aided Orthopedic Surgery”
4.EP 08156293.6
5.US 61/054,187
优选地,基于成像几何构造信息,计算机计算出治疗束110相对于目标(治疗患体部位)的位置。通过比较针对治疗束所确定的位置与针对目标所确定的位置,能够对治疗束的位置进行控制,从而击中目标。
根据本文所述的方法所确定的目标变化数据优选用于正确地控制治疗束110相对于患者400的位置。表示所述方法的实施方案的程序优选存储和/或在运行在计算机200上。指示器I1可以是附接在身体某部位(例如胸部以跟踪真肋的移动)上的标记物。另一个指示器I2例如是隔膜,其位置通过分析装置310和320来跟踪。目标T被治疗束110击中。标记物I1的位置例如可通过(立体)相机来检测,如根据导航系统得知的。优选地是,使用分别表示指示器结构的不同关键移动的至少两个指示器。为了确定肿瘤移动的起始位置以及为了将肿瘤运动的起始位置与指示器(I1,I2)的位置相关联,例如可以产生CBCT图像,从而确定具体说指示器轨迹和/或目标所处的参考系统中肿瘤的位置。优选地,在通过分析装置(例如CBCT)确定肿瘤位置的同时,确定指示器(如I1和I2)的位置。根据本发明肿瘤的位置变化随后可通过确定目标变化数据的方式来确定,并且可基于目标的起始位置以及所确定的目标位置变化(即目标变化数据)来确定目标的位置。
图2示出了左侧的第一指示器I1的循环轨迹10和右侧的第二指示器I2的循环轨迹20。第一指示器I1的轨迹路径上的位置由向量和描述。n元向量(和)是针对的示例,并且描述了关于指示器I1和I2的位置变化的信息。指示器I1在上述不同位置之间的位置变化由和来表示。指示器沿着循环轨迹的不同位置可描述在其中设有轨迹的参考系统中。类似地,图2的右侧示出了第二指示器I2的循环轨迹。沿着轨迹20的不同位置由和来表示。这些位置之间的位置变化是由如下向量来表示: 再次,描述位置的向量可以在其中设有轨迹20的参考系统中来描述。优选地,指示器(如I1和I2)的轨迹在其中设有轨迹的参考系统中进行描述。
根据一实施例,基于对沿轨迹的指示器位置的测量和对目标位置的相应测量(具体说是同时测量),产生描述了指示器I1和I2以及备选的其它指示器的位置与目标位置之间的对应关系的对应数据。所述测量具体说在治疗之前例如通过使用四维成像法来执行。图3中的简化查找表示出了指示器-目标对应数据的简单示例,在上面的行中,向量表示指示器I1的位置。在表中,目标I1的四个位置是:和这些位置记录在查找表的上面的行中。图3的查找表的左列示出了指示器I2的位置和图3的查找表显示了目标的哪个位置与两个指示器指标I1和I2的哪些位置相对应。例如,位置对应于指示器位置和如果,例如目标的先前位置为而当前位置为那么描述了目标变化的信息的n元向量例如是()。
指示器I1的位置变化标记为指示器I1的位置变化例如标记为而指示器I2的位置变化则标记为也就是,指示器变化集合例如可描述为( )。变化集合还可包括关于指示器位置的信息,针对指示器变化集合的示例也可以是:()。针对指示器变化集合的特定示例例如( )或()。位置变化还可以根据后面的示例导出,例如,
为了生成查找表,例如也可使用具有相同类型的待治疗目标(具体说是,在相同目标结构内的相同位置处的待治疗目标)的多个患者(患者池)。优选地,相同的指示器(具体说是,相对于身体而定位在相同位置处的指示器)用于确定指示器的轨迹。优选地,对这些指示器的轨迹求平均,从而针对不同类型的指示器中的每一种指示器生成通用的轨迹。基于所生成的通用轨迹,可产生通用的对应数据。图3所示的示例表示了指示器目标对应数据的示例。如果那些指示器目标对应数据是基于对待治疗患者进行的测量而产生的,那么它们则称为个体指示器目标对应数据。如果它们是通过采用患者池产生的,那么它们则称为通用指示器目标对应数据。
图4a)示出了目标结构的模型。该模型包括表示目标结构的内部结构的网格。当边界被压缩或扩张时,假设网格线具有弹性性能(例如,类似弹簧)。目标T放置在模型内。优选地,对网格线进行选择,从而使目标放置在由网格线建立的节点上。模型优选为三维的。为简单起见,图4中仅示出了二维简化模型。网格线在称为边界点位置或仅称为边界位置的位置至处与模型边界接触。根据一实施例,例如通过使用查找表来提供目标结构对应数据。查找表描述了指示器位置与边界位置至之间的关系。在图4a)给出的示例中,其示出了变化前的目标结构的模型边界,而图4b)则示出了变化后的目标结构的模型边界。在图4给出的示例中,模型在垂直方向上压缩并在水平方向上扩张。这导致了目标T向下且向右的位置变化。目标结构的变形可例如由位于目标结构上方且向下移动的相邻身体结构引起。由此可见,这种向下移动将不仅仅导致目标的向下移动,而且还将导致目标向右移动。这是由于目标结构的机械响应造成的,并且解释了目标结构的此类移动与指示器结构及指示器的移动之间的复杂关系。
在图4b)中,边界位置至标记为至边界位置至表示变化之后的边界点的位置。附图标记T’表示变化之后的目标位置。
假设图4a)所示的网格表示目标结构的图像,那么图4b)所示的网格则作为将图4a)的图像融合到图4b)所示的目标结构的边界之内的示例。由于这种弹性融合,还使得目标位置从T变化到T’。
Claims (10)
1.一种包括计算机的装置,所述计算机执行用于确定称为目标变化数据的方法,所述目标变化数据能用于执行放射治疗,所述目标变化数据描述了包含在患者体内的目标的位置变化的信息,其中所述方法通过所述计算机来执行并包括以下步骤:
·获取描述映射规则的确定规则数据,其中所述映射规则用于将指示器变化集合映射到所述目标变化数据上,所述指示器变化集合是包括多于一个元素的集合,所述指示器变化集合的元素分别表示指示器位置的各个变化的信息,至少两个指示器分别指示称为指示器结构的不同身体结构的位置变化,
·获取描述所获取的指示器变化集合的集合数据;以及
·通过将所述映射规则应用于所获取的指示器变化集合上,基于所述确定规则数据和所述集合数据来确定所述目标变化数据,
·其中,嵌有目标的身体结构称为目标结构,并且为了确定所述目标变化数据,通过使用所述目标结构的模型,根据指示器变化集合来确定所述目标结构的几何构造变化和/或位置变化,
其中所述模型描述目标结构对目标结构边界上的变化的机械响应,所述变化与指示器变化集合相关联。
2.根据权利要求1所述的装置,其中,所确定的目标变化数据描述了以下至少一项:
a)所述目标的位置变化;
b)所述目标的位置变化和变化前的目标位置;以及
c)变化后的目标位置。
3.根据权利要求1或2所述的装置,其中,关于指示器位置的各个变化的信息描述了变化发生时的指示器位置变化和指示器位置,其中所述映射规则取决于所述指示器位置,和/或
还包括获取目标位置数据的步骤,所述目标位置数据描述了目标位置变化之前的目标位置,并且所述映射规则取决于所述目标位置数据;和/或
其中,关于指示器位置的各个变化的信息描述指示器的位置变化以及在所述指示器经历往复移动的情况下的位置变化的最大振幅和/或位置变化的频率,其中所述映射规则取决于指示器的位置变化和所述位置变化的最大振幅和/或所述位置变化的频率。
4.根据权利要求1或2所述的装置,其中,所述映射规则包括所获取的、待映射的指示器变化集合的元素之间的数学关系。
5.根据权利要求1或2所述的装置,
其中所述目标结构的位置变化和/或几何构造变化称为目标结构变化,并且
所述映射规则还基于:
个体目标结构对应数据,其描述了多个不同的指示器变化集合与相应的目标结构变化之间的对应关系,所述个体目标结构对应数据已通过分析待治疗的患者身体而产生;并且
其中,所述目标结构变化基于所获取的指示器变化集合和所述个体目标结构对应数据来确定,以及
其中,所述映射规则还基于子规则,所述子规则是基于所述目标结构的模型并且描述了将所确定的目标结构变化映射到所述目标变化数据上的规则;以及
其中,所述目标变化数据基于所确定的目标结构变化和所述子规则来确定。
6.根据权利要求1所述的装置,
其中,所述目标结构的位置变化和/或几何构造变化称为目标结构变化,并且
所述映射规则还基于:
通用目标结构对应数据,其描述了多个不同的指示器变化集合与相应的目标结构变化之间的对应关系,所述通用目标结构对应数据已通过分析多个不同患者的身体而产生;并且
其中,
所述目标结构变化基于所获取的指示器变化集合和所述通用目标结构对应数据来确定,以及
其中,所述映射规则还基于子规则,所述子规则是基于所述目标结构的模型并且描述了将所确定的目标结构变化映射到所述目标变化数据上的规则;以及
其中,所述目标变化数据基于所确定的目标结构变化和所述子规则来确定。
7.根据权利要求6所述的装置,所述装置所执行的方法还包括以下步骤:
获取目标位置数据,所述目标位置数据描述了变化之前的所述目标结构中的目标位置;
其中所述子规则包括描述了所述目标结构的物理参数的所述目标结构模型,所述物理参数规定了由于所述目标结构的边界变化所导致的所述目标结构的内部结构的变化,所述模型基于目标位置数据模仿了由于所述内部结构的变化所导致的所述目标结构中的目标位置的变化;并且
其中,所述目标变化数据的确定还基于所述目标位置数据。
8.根据权利要求6或7所述的装置,其中,其中所述模型是所述目标结构的图像并且所述映射规则包括子规则,所述子规则基于将表示目标结构变化前含有所述目标的目标结构的目标结构图像弹性融合到根据所确定的目标结构变化而产生的目标结构图像中。
9.根据权利要求1或2所述的装置,其中,由所获取的变化集合的元素来表示其移动的指示器包括以下指示器中的至少两个:根据隔膜运动而移动的第一指示器、根据假肋而移动的第二指示器、根据真肋而移动的第三指示器、根据心脏跳动而移动的第四指示器以及根据肠运动而移动的第五指示器。
10.一种治疗束系统,包括:
如权利要求1到9中任一项所述的装置,其构造成基于所获取的指示器变化集合来确定所述目标变化数据;和
治疗装置(100),其构造成根据所确定的目标变化数据来发射治疗束并控制所述治疗束的位置。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/EP2012/053718 WO2013131547A1 (en) | 2012-03-05 | 2012-03-05 | Using different indicators for determining positional changes of a radiotherapy target |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104145295A CN104145295A (zh) | 2014-11-12 |
CN104145295B true CN104145295B (zh) | 2017-08-25 |
Family
ID=45787214
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201280071164.2A Active CN104145295B (zh) | 2012-03-05 | 2012-03-05 | 使用不同指示器来确定放射治疗目标的位置变化 |
Country Status (5)
Country | Link |
---|---|
US (1) | US9911199B2 (zh) |
EP (1) | EP2823466B1 (zh) |
JP (1) | JP6099061B2 (zh) |
CN (1) | CN104145295B (zh) |
WO (1) | WO2013131547A1 (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015087217A1 (en) | 2013-12-12 | 2015-06-18 | Koninklijke Philips N.V. | Real-time fusion of anatomical ultrasound information and radiation delivery information for radiation therapies |
US10750980B2 (en) * | 2014-12-02 | 2020-08-25 | Brainlab Ag | Human body measurement using thermographic images |
WO2016096038A1 (en) * | 2014-12-19 | 2016-06-23 | Brainlab Ag | Method for optimising the position of a patient's body part relative to an irradiation source |
EP3411118B1 (en) | 2016-02-02 | 2020-09-16 | Elekta Ltd. | Three dimensional respiratory motion management in image-guided radiotherapy |
CN105816196A (zh) * | 2016-05-13 | 2016-08-03 | 上海联影医疗科技有限公司 | 用于4dct成像的标记带及4dct成像方法 |
JP6727644B2 (ja) * | 2016-05-23 | 2020-07-22 | 株式会社日立製作所 | 動体追跡装置および放射線照射システム |
US11446520B2 (en) * | 2017-03-14 | 2022-09-20 | National Univ. Corporation Hokkaido Univ. | Radiation therapy apparatus configured to track a tracking object moving in an irradiation object |
CN111428709B (zh) * | 2020-03-13 | 2023-10-24 | 平安科技(深圳)有限公司 | 图像处理方法、装置、计算机设备和存储介质 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5418708A (en) | 1994-06-10 | 1995-05-23 | The United States Of America As Represented By The Secretary Of The Air Force | Constant power load bank |
US6144875A (en) * | 1999-03-16 | 2000-11-07 | Accuray Incorporated | Apparatus and method for compensating for respiratory and patient motion during treatment |
US8989349B2 (en) * | 2004-09-30 | 2015-03-24 | Accuray, Inc. | Dynamic tracking of moving targets |
US8747382B2 (en) * | 2005-04-13 | 2014-06-10 | University Of Maryland, Baltimore | Techniques for compensating movement of a treatment target in a patient |
US10279196B2 (en) * | 2006-09-28 | 2019-05-07 | Accuray Incorporated | Radiation treatment planning using four-dimensional imaging data |
CN101861185B (zh) * | 2007-12-07 | 2013-06-12 | 三菱重工业株式会社 | 放射线治疗计划装置以及放射线治疗计划方法 |
EP2119397B1 (de) | 2008-05-15 | 2013-12-18 | Brainlab AG | Bestimmung einer Kalibrier-Information für ein Röntgengerät |
EP2189943B1 (de) | 2008-11-19 | 2014-09-24 | Brainlab AG | Bestimmung vital-bewegter Regionen eines Analysebildes |
WO2011107145A1 (en) | 2010-03-02 | 2011-09-09 | Brainlab Ag | Tracking representations of indicator body parts |
-
2012
- 2012-03-05 WO PCT/EP2012/053718 patent/WO2013131547A1/en active Application Filing
- 2012-03-05 JP JP2014560255A patent/JP6099061B2/ja active Active
- 2012-03-05 US US14/383,121 patent/US9911199B2/en active Active
- 2012-03-05 CN CN201280071164.2A patent/CN104145295B/zh active Active
- 2012-03-05 EP EP12706856.7A patent/EP2823466B1/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN104145295A (zh) | 2014-11-12 |
US9911199B2 (en) | 2018-03-06 |
WO2013131547A1 (en) | 2013-09-12 |
US20150038765A1 (en) | 2015-02-05 |
JP2015509409A (ja) | 2015-03-30 |
JP6099061B2 (ja) | 2017-03-22 |
EP2823466B1 (en) | 2017-05-03 |
EP2823466A1 (en) | 2015-01-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104145295B (zh) | 使用不同指示器来确定放射治疗目标的位置变化 | |
US11663755B2 (en) | Determination of dynamic DRRs | |
EP3461296B1 (en) | Multi image fusion based positioning verification | |
EP3468668B1 (en) | Soft tissue tracking using physiologic volume rendering | |
JP2019511268A (ja) | 脳深部刺激療法の電極の三次元画像における回転配向の決定 | |
US11715212B2 (en) | Heatmap and atlas | |
KR20240013724A (ko) | 다중 펄스 x-선 소스 이동 단층합성 영상화 시스템을 사용한 인공 지능 훈련 | |
CN108885781A (zh) | 用于根据低剂量或低kV计算机断层扫描图像合成虚拟高剂量或高kV计算机断层扫描图像的方法和系统 | |
EP3790626B1 (en) | Computation of a breathing curve for medical applications | |
US9367926B2 (en) | Determining a four-dimensional CT image based on three-dimensional CT data and four-dimensional model data | |
JP2021532903A (ja) | 医療装置を撮像するためのコンセンサス面の決定 | |
Özbek et al. | Optimizing external surface sensor locations for respiratory tumor motion prediction | |
Wen | Registration of 3D models to 2D X-ray images using fast X-ray simulation and global optimisation algorithms |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information |
Address after: Munich, Germany Applicant after: BRAINLAB AG Address before: Feldkirchen Applicant before: BRAINLAB AG |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |