CN115212475A - 用于评价直线加速器等中心质量的x射线透射图像分析 - Google Patents
用于评价直线加速器等中心质量的x射线透射图像分析 Download PDFInfo
- Publication number
- CN115212475A CN115212475A CN202210412279.9A CN202210412279A CN115212475A CN 115212475 A CN115212475 A CN 115212475A CN 202210412279 A CN202210412279 A CN 202210412279A CN 115212475 A CN115212475 A CN 115212475A
- Authority
- CN
- China
- Prior art keywords
- radiation
- radiation beam
- angle
- gantry
- beam axis
- 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
- 230000005540 biological transmission Effects 0.000 title claims abstract description 66
- 238000010191 image analysis Methods 0.000 title description 2
- 230000005855 radiation Effects 0.000 claims abstract description 530
- 239000003550 marker Substances 0.000 claims abstract description 145
- 238000000034 method Methods 0.000 claims abstract description 107
- 239000013598 vector Substances 0.000 claims description 50
- 230000033001 locomotion Effects 0.000 claims description 47
- 238000003384 imaging method Methods 0.000 claims description 34
- 206010028980 Neoplasm Diseases 0.000 claims description 18
- 238000004590 computer program Methods 0.000 description 14
- 238000001959 radiotherapy Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000013500 data storage Methods 0.000 description 2
- 210000003484 anatomy Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010894 electron beam technology Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000011275 oncology therapy Methods 0.000 description 1
- 238000004513 sizing Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- 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/1075—Monitoring, verifying, controlling systems and methods for testing, calibrating, or quality assurance of the radiation treatment apparatus
-
- 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
-
- 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
-
- 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
-
- 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/1077—Beam delivery systems
- A61N5/1081—Rotating beam systems with a specific mechanical construction, e.g. gantries
-
- G—PHYSICS
- G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
- G21K—TECHNIQUES FOR HANDLING PARTICLES OR IONISING RADIATION NOT OTHERWISE PROVIDED FOR; IRRADIATION DEVICES; GAMMA RAY OR X-RAY MICROSCOPES
- G21K5/00—Irradiation devices
- G21K5/08—Holders for targets or for other objects to be irradiated
-
- 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
- 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
- A61N2005/1054—Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam using a portal imaging system
-
- 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/1075—Monitoring, verifying, controlling systems and methods for testing, calibrating, or quality assurance of the radiation treatment apparatus
- A61N2005/1076—Monitoring, verifying, controlling systems and methods for testing, calibrating, or quality assurance of the radiation treatment apparatus using a dummy object placed in the radiation field, e.g. phantom
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Veterinary Medicine (AREA)
- Pathology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Radiology & Medical Imaging (AREA)
- Surgery (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- High Energy & Nuclear Physics (AREA)
- Radiation-Therapy Devices (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
用于确定直线加速器(LINAC)的辐射等中心的***和方法。确定辐射等中心可以包括根据二维(2D)辐射透射图像确定LINAC的一组三维(3D)辐射束轴线。可以至少基于一组3D辐射束轴线来确定辐射等中心。确定一组3D辐射束轴线可以包括基于确定的在LINAC的机架相对于参考机架角度以一定角度的情况下生成的辐射束的束轴线的位置、确定的辐射束的辐射场中的辐射不透过标记的阴影的中心以及机架角度来构造3D辐射束轴线。
Description
相关申请的交叉引用
本申请要求2021年4月19日提交的美国临时专利申请序列号为No.63/176,494的优先权,其通过引用整体并入本文。
技术领域
本发明大体上涉及辐射疗法,其中医用直线加速器(LINAC)将辐射束传递到患者体内的精确点。特别地,本发明涉及分析和显示用于评估LINAC的辐射等中心的质量的二维x射线透射图像。
背景技术
1.1.LINAC概述
辐射疗法是一种癌症疗法,其使用强高能辐射束的束来杀死癌细胞。附图说明1-3示出了包括机架(gantry)102和诊察台(couch)106的医用直线加速器(LINAC)100。如图4所示,机架102包括准直器头104,其限定LINAC 100生成的辐射束412。在治疗期间,LINAC 100的机架102将辐射束传送到支撑在诊察台106上的患者(未示出)体内的精确位置。为了定制辐射治疗计划,机架102可以围绕患者旋转,并且诊察台106可以相对于机架102的旋转平面改变患者的方位。
如图2所示,机架102能够围绕患者旋转整个360°,以便优化辐射束进入患者体内。如图2所示,机架102可以围绕机架旋转轴线208旋转。如图3所示,诊察台106可围绕诊察台旋转轴线310旋转。这允许相对于机架102的旋转平面调整患者的方位,以便进一步优化辐射束104进入患者体内。
在理想的LINAC 100中,当机架102旋转时,辐射束将保持聚焦在空间中的固定点上(即,在患者的肿瘤上)。类似地,当患者由于诊察台106的旋转而移动时,理想的LINAC100将保持肿瘤固定在辐射束内。
如图4所示,LINAC 100可以具有成像装置416(例如,电子入口成像装置(EPID)),该成像装置416用于通过测量在治疗过程期间从辐射端口透射通过患者的X射线强度来在治疗之前验证患者的位置。
1.2辐射等中心
辐射疗法的目标是向肿瘤递送最高可能的辐射剂量,同时使对周围健康组织的剂量最小化。为此,辐射疗法包括以下目标:(1)当机架102旋转时保持辐射束精确地聚焦在患者的肿瘤上,以及(2)当诊察台106旋转时保持患者的肿瘤固定在辐射束内。然而,这两个目标很少被完美地实现。机架102通常在其自身的重量下弯曲,从而有效地使辐射束模糊。诊察台106的旋转通常也不是完全圆形的和/或诊察台106的旋转轴线310通常不与机架102的旋转轴线208相交。
为了使这些常见的LINAC缺陷的影响最小化,肿瘤应该被放置在空间中使最大的束对靶误差最小化的点处。该位置被称为等中心点。
1.3等中心尺寸
对于机架和诊察台旋转的所有组合的最大束轴线错位距离(beam axis missdistance)是等中心尺寸。较小的等中心点将导致更精确的治疗。
1.4现有方法
用于分析等中心点的现有方法包括两个主要步骤:(1)利用不同的机架和诊察台旋转测量位于(或接近)等中心点处的辐射不透过标记(radiation-opaque marker)的EPIDX射线透射图像;以及(2)分析图像以确定等中心点的尺寸。
图4所示的设置通常用于完成第一步骤。在图4中,LINAC 100通过放置在等中心处的辐射不透过标记414传送辐射束412。成像装置416在辐射束412通过标记414时获取辐射束412的二维透射图像。
在图5中示出了包括不透过辐射标记414的标记组件500的详细视图。在图5中,高密度球形标记414连接到低密度支撑杆502,该低密度支撑杆502然后连接到基部504。基部504设置在诊察台106的顶部,使得标记414可位于(或接近)等中心。
图6示出了通过标记414的典型图像600(例如,EPID图像)。在图6中,暗方形区域602由曝光成像装置410的辐射束412的方形辐射场产生,而较亮的内部圆形604由位于辐射场内的辐射不透过标记414的阴影产生。
用于测量等中心点的现有方法包括将机架102旋转通过一组固定机架角度,同时诊察台106保持在0°,并获取每个机架位置的EPID图像。当机架102保持固定在0°并且诊察台106旋转通过一组固定诊察台角度时,获得另外的EPID图像。
现有方法的第二步骤包括分析上述图像以确定等中心点的尺寸。该步骤可以进一步分解为两个子过程:(2a)使用图像处理技术来找到每个图像(具有确定的场和标记位置的图像700的示例在图7中示出)中的场和标记414的中心;以及(2b)为所有获取的图像计算最大标记中心与场中心的距离,并将该最大值分配给等中心尺寸。
发明内容
现有的用于分析等中心而获取的所有EPID图像的场-标记距离的使用利用二维坐标系(CS)对复杂三维过程进行建模。现有的方法维护二维图像CS中的数据,该数据从EPID图像获取。通过这样做,当问题出现在LINAC等中心时,常常难以直观地理解问题的性质。本发明的各方面可以克服现有方法的一个或多个问题(例如,通过替换和改进现有方法的子过程步骤2b)。
本发明的方面可以包括(1)使用直线加速器(LINAC)来获取不同机架和/或诊察台旋转的等中心处或附近定位的辐射不透过标记的辐射透射图像,以及(2)分析图像以确定等中心的尺寸、标记放置误差(即,相对于当前标记位置的等中心位置)和/或诊察台走出量(couch walkout)。分析图像可以包括在图像中找到辐射场和标记的中心。
关于确定等中心尺寸,不是比较场和标记中心以找到最大差异并将最大距离指定为等中心尺寸,而是LINAC可以将场和标记中心从二维(2D)图像坐标系(CS)变换为三维(3D)(例如,真实世界)CS。从2D图像CS移动到3D CS可以提供允许在3D中更灵活的分析和直观的显示的优点,其更准确地反映了底层过程。
在一些方面,在分析图像以表征LINAC的等中心时,LINAC可以不同地处理机架和诊察台旋转。当机架旋转时,场移动,但标记不移动。类似地,当诊察台旋转时,标记移动,但场保持静止。在一些方面,LINAC因此可以将图像分组为两个类别:(1)用于确定等中心和/或确定标记放置误差的图像,其可以在诊察台保持在固定诊察台角度(例如,0°)时利用机架旋转来获取,以及(2)用于确定诊察台走出量的图像,其可以在机架保持在固定机架角度(例如,0°)时利用诊察台旋转来获取。
本发明的一个方面可以提供一种方法,包括根据二维(2D)辐射投射图像来确定直线加速器(LINAC)的一组三维(3D)辐射束轴线,其中该组3D辐射束轴线至少包括第一和第二3D辐射束轴线。该方法可以包括至少基于该组3D辐射束轴线来确定LINAC的辐射等中心。
在一些方面,确定第一3D辐射束轴线可以包括相对于参考机架角度以第一机架角度定位LINAC的机架,并且定位机架可以包括使机架围绕机架旋转轴线旋转。在一些方面,确定第一3D辐射束轴线可以包括,在以第一机架角度定位机架的情况下,使用LINAC生成第一辐射束。在一些方面,确定第一3D辐射束轴线可以包括,在以第一机架角度定位机架的情况下,使用LINAC的成像装置来获取第一二维(2D)辐射透射图像,第一二维(2D)辐射透射图像指示第一辐射束通过辐射不透过标记之后的辐射场。在一些方面,确定第一3D辐射束轴线可以包括基于第一2D辐射透射图像,确定第一辐射束的束轴线的位置和第一辐射束的辐射场中的辐射不透过标记的阴影的中心。在一些方面,确定第一3D辐射束轴线可以包括基于确定的第一辐射束的束轴线的位置、确定的第一辐射束的辐射场中的辐射不透过标记的阴影的中心以及第一机架角度,构造第一3D辐射束轴线。
在一些方面,构造第一3D辐射束轴线可以包括,使第一2D辐射透射图像围绕第一图像旋转轴线旋转,第一图像旋转轴线与确定的第一辐射束的辐射场中的辐射不透过标记的阴影的中心一致并且与机架旋转轴线平行。在一些方面,第一2D辐射透射图像可以旋转与第一机架角度相等的量。在一些方面,第一3D辐射束轴线可以与确定的旋转的第一2D辐射透射图像中的第一辐射束的束轴线的位置一致并且与旋转的第一2D辐射透射图像的平面垂直。
在一些替代方面,构造所述第一3D辐射束轴线可以包括:生成初始3D辐射束轴线,初始3D辐射束轴线与确定的第一2D辐射透射图像中的第一辐射束的束轴线的位置一致并且与第一2D辐射透射图像的平面垂直;和使初始3D辐射束轴线围绕第一图像旋转轴线旋转,第一图像旋转轴线与确定的第一辐射束的辐射场中的辐射不透过标记的阴影的中心一致并且与机架旋转轴线平行;并且初始3D辐射束轴线可以旋转与第一机架角度相等的量。
在一些方面,确定第二3D辐射束轴线可以包括相对于参考机架角度以第二机架角度定位LINAC的机架。在一些方面,确定第二3D辐射束轴线可以包括,在以第二机架角度定位机架的情况下,使用LINAC生成第二辐射束。在一些方面,确定第二3D辐射束轴线可以包括,在以第二机架角度定位机架的情况下,使用LINAC的成像装置获取第二2D辐射透射图像,第二2D辐射透射图像指示第二辐射束经过辐射不透过标记之后的辐射场。在一些方面,确定第二3D辐射束轴线可以包括基于第二2D辐射透射图像,确定第二辐射束的束轴线的位置和第二辐射束的辐射场中的辐射不透过标记的阴影的中心。在一些方面,确定第二3D辐射束轴线可以包括基于确定的第二辐射束的束轴线的位置、确定的第二辐射束的辐射场中的辐射不透过标记的阴影的中心以及第二机架角度,构造第二3D辐射束轴线。
在一些方面中,构造第二3D辐射束轴线可以包括使所述第二2D辐射透射图像围绕第二图像旋转轴线旋转,第二图像旋转轴线与确定的第二辐射束的辐射场中的辐射不透过标记的阴影的中心一致并且与机架旋转轴线平行。在一些方面,第二2D辐射透射图像可以旋转与第二机架角度相等的量。在一些方面,第二3D辐射束轴线可以与确定的旋转的第二2D辐射透射图像中的第二辐射束的束轴线的位置一致并且与旋转的第二2D辐射透射图像的平面垂直。
在一些替代方面,构造第二3D辐射束轴线可以包括:生成初始第二3D辐射束轴线,初始第二3D辐射束轴线与确定的第二2D辐射透射图像中的第二辐射束的束轴线的位置一致并且与第二2D辐射透射图像的平面垂直;和使初始第二3D辐射束轴线围绕第二图像旋转轴线旋转,第二图像旋转轴线与确定的第二辐射束的辐射场中的辐射不透过标记的阴影的中心一致并且与机架旋转轴线平行;并且初始第二3D辐射束轴线可以旋转与第二机架角度相等的量。
在一些方面,在LINAC用于生成第一辐射束和第二辐射束并且LINAC的成像装置用于获取第一辐射透射图像和第二2D辐射透射图像的同时,可以以固定诊察台角度定位诊察台。
在一些方面,确定辐射等中心可以包括,对于一组3D辐射束轴线中的每个3D辐射束轴线,确定3D空间中的位置和3D辐射束轴线之间的束轴线错位距离。在一些方面,确定辐射等中心可以包括确定对于一组3D辐射束轴线中的3D辐射束轴线确定的束轴线错位距离中哪一个最大,并且对于一组3D辐射束轴线中的3D辐射束轴线确定的束轴线错位距离中最大的可以是位置的最大束轴线错位距离。在一些方面,确定辐射等中心可以包括在3D空间中找到具有最小的最大束轴线错位距离的位置。在一些方面,3D空间中的位置和3D辐射束轴线之间的束轴线错位距离可以是该位置和3D辐射束轴线之间的最短距离。在一些方面,辐射等中心的尺寸可以是该位置的最大射束轴线错位距离。
在一些方面,该方法还可以包括将肿瘤定位在确定的辐射等中心。
在一些方面,该方法可以进一步包括对于一组诊察台角度中的每个诊察台角度确定标记移动矢量(marker movement vector),并且一组诊察台角度至少可以包括相对于参考诊察台角度的第一诊察台角度和相对于参考诊察台角度的第二诊察台角度。在一些方面,该方法可以包括基于确定的标记移动矢量,对于一组诊察台角度中的诊察台角度与一组3D辐射束轴线中的3D辐射束轴线的每个组合确定标记与3D辐射束轴线误差距离,并且标记与3D辐射束轴线误差距离表示辐射不透过标记的位置和3D辐射束轴线误差之间的最短距离。
在一些方面中,确定用于第一诊察台角度的标记移动矢量可以包括,在以参考诊察台角度定位LINAC的诊察台并且以参考机架角度定位LINAC的机架的情况下,使用LINAC生成第一参考辐射束。在一些方面,在以参考诊察台角度定位诊察台并且以参考机架角度定位机架的情况下,使用LINAC的成像装置获取参考诊察台角度二维(2D)辐射透射图像,参考诊察台角度二维(2D)辐射透射图像指示第一参考辐射束经过辐射不透过标记之后的辐射场。在一些方面,确定用于第一诊察台角度的标记移动矢量可以包括基于参考诊察台角度2D辐射透射图像,确定第一参考辐射束的辐射场中的辐射不透过标记的阴影的中心。在一些方面,确定用于第一诊察台角度的标记移动矢量可以包括以第一诊察台角度定位诊察台。在一些方面,确定用于第一诊察台角度的标记移动矢量可以包括,在以第一诊察台角度定位诊察台并且以参考机架角度定位机架的情况下,使用LINAC生成第二参考辐射束。在一些方面,确定用于第一诊察台角度的标记移动矢量可以包括,在以第一诊察台角度定位诊察台并且以参考机架角度定位机架的情况下,使用LINAC的成像装置获取第一诊察台角度2D辐射透射图像,第一诊察台角度2D辐射透射图像指示第二参考辐射束经过辐射不透过标记之后的辐射场。在一些方面,确定用于第一诊察台角度的标记移动矢量可以包括基于第一诊察台角度2D辐射透射图像,确定第二参考辐射束的辐射场中的辐射不透过标记的阴影的中心。在一些方面,确定用于第一诊察台角度的标记移动矢量可以包括基于第一参考辐射束的辐射场中的辐射不透过标记的阴影的中心与第二参考辐射束的辐射场中的辐射不透过标记的阴影的中心的比较,确定第一诊察台角度的标记移动矢量。
在一些方面,确定第二诊察台角度的标记移动矢量可以包括以第二诊察台角度定位诊察台。在一些方面,确定第二诊察台角度的标记移动矢量可以包括,在以第二诊察台角度定位诊察台并且以参考机架角度定位机架的情况下,使用LINAC生成第三参考辐射束。在一些方面,确定第二诊察台角度的标记移动矢量可以包括,在以第二诊察台角度定位诊察台并且以参考机架角度定位机架的情况下,使用LINAC的成像装置获取第二诊察台角度2D辐射透射图像,第二诊察台角度2D辐射透射图像指示第三参考辐射束经过辐射不透过标记之后的辐射场。在一些方面,确定第二诊察台角度的标记移动矢量可以包括基于第二诊察台角度2D辐射透射图像,确定第三参考辐射束的辐射场中的辐射不透过标记的阴影的中心。在一些方面,确定第二诊察台角度的标记移动矢量可以包括基于第一参考辐射束的辐射场中的辐射不透过标记的阴影的中心与第三参考辐射束的辐射场中的辐射不透过标记的阴影的中心的比较,确定第二诊察台角度的标记移动矢量。
在一些方面,该方法可以进一步包括确定临床等中心,其中,确定临床等中心包括确定使最大的标记与3D辐射束轴线误差距离最小化的空间位置。在一些方面,确定使最大的标记与3D辐射束轴线误差距离最小化的空间位置可以包括基于不同诊察台角度下的辐射不透过标记的阴影,预测诊察台旋转轴线。在一些方面,确定使最大的标记与3D辐射束轴线误差距离最小化的空间位置可以包括预测辐射不透过标记在新的参考位置处的一组位置,该一组位置通过围绕预测的诊察台旋转轴线的诊察台旋转而移位。在一些方面,确定使最大的标记与3D辐射束轴线误差距离最小化的空间位置可以包括基于确定的标记移动矢量,对于一组诊察台角度中的诊察台角度和一组3D辐射束轴线中的3D辐射束轴线的每个组合确定预测的标记与3D辐射束轴线误差距离,并且预测的标记与3D辐射束轴线误差距离可以表示辐射不透过标记的位置和3D辐射束轴线误差之间的最短距离。在一些方面,确定使最大的标记与3D辐射束轴线误差距离最小化的空间位置可以包括确定使最大的标记与束轴线误差距离最小化的参考标记位置。在一些方面,该方法可以进一步包括将肿瘤置于确定的临床等中心。
本发明的另一方面可以提供一种设备,该设备被配置为:根据二维(2D)辐射透射图像,确定直线加速器(LINAC)的一组三维(3D)辐射束轴线,并且一组3D辐射束轴线包括至少第一3D辐射束轴线和第二3D辐射束轴线。该设备可以被配置为至少基于一组3D辐射束轴线,确定LINAC的辐射等中心。
本发明的又一方面可以提供一种方法,该方法包括对于直线加速器(LINAC)的一组诊察台角度中的每个诊察台角度确定标记移动矢量,并且一组诊察台角度至少包括相对于参考诊察台角度的第一诊察台角度和相对于参考诊察台角度的第二诊察台角度。该方法卡伊包括基于确定的标记移动矢量,对于一组诊察台角度中的诊察台角度与一组3D辐射束轴线中的3D辐射束轴线的每个组合确定标记与3D辐射束轴线误差距离,并且标记与3D辐射束轴线误差距离表示辐射不透过标记的位置和3D辐射束轴线误差之间的最短距离。
本发明的又一方面可以提供一种设备,该设备被配置为对于直线加速器(LINAC)的一组诊察台角度中的每个诊察台角度确定标记移动矢量,并且一组诊察台角度可以至少包括相对于参考诊察台角度的第一诊察台角度和相对于参考诊察台角度的第二诊察台角度。该设备可以被配置为基于确定的标记移动矢量,对于一组诊察台角度中的诊察台角度与一组3D辐射束轴线中的3D辐射束轴线的每个组合确定标记与3D辐射束轴线误差距离,并且标记与3D辐射束轴线误差距离可以表示辐射不透过标记的位置和3D辐射束轴线误差之间的最短距离。
本发明的又一方面可以提供一种计算机程序,其包括用于使设备适于执行上述任何方法的指令。本发明的又一方面可以提供包含计算机程序的载体,并且载体可以是电信号、光信号、无线电信号或计算机可读存储介质中的一种。
本发明的又一方面可以提供一种包括处理电路和存储器的设备。该存储器可以包含可以由处理电路执行的指令,由此该设备可操作以执行上述任何方法。
本发明的又一方面可以提供一种适于上述任何方法的设备。
本发明的又一方面可以提供上述方面的任何组合。
在下面的本发明的详细描述中描述了包含在***和方法内的进一步变化。
附图说明
结合于此并形成说明书的一部分的附图示出了本发明的各种非限制性实施例。在附图中,相同的附图标记表示相同或功能相似的元件。
图1示出了医用直线加速器(LINAC)。
图2示出LINAC的机架的旋转。
图3示出了LINAC的诊察台的旋转。
图4示出了LINAC的辐射束和成像面板。
图5示出了根据一些方面的标记组件。
图6示出了根据一些方面的未处理的辐射透射图像。
图7示出了根据一些方面的经处理的辐射透射图像。
图8示出了根据一些方面的辐射不透过标记和辐射场中心的位置的差异。
图9示出了根据一些方面的经处理的二维图像,其示出了找到的辐射场、辐射场的中心、标记的找到的阴影以及辐射场中心轴线。
图10示出了根据一些方面的原始二维图像围绕标记的阴影的中心旋转到新位置的量和方向,该量和方向对应于图像采集时机架旋转的量和方向。
图11示出了根据一些方面的已经从二维图像坐标变换成三维真实世界坐标并限定等中心的束轴线的叠加。
图12示出了根据一些方面的在不同诊察台旋转角度处的诊察台标记走出位置以及用于机架和诊察台角度的组合的肿瘤束轴线错位距离的计算。
图13示出了根据一些方面的过程。
图14示出了根据一些方面的用于确定第一3D辐射束轴线的过程。
图15示出了根据一些方面的用于确定辐射等中心点的过程。
图16示出了根据一些方面的用于确定诊察台走出量的过程。
图17示出了根据一些方面的用于确定诊察台角度的标记移动矢量的过程。
图18示出了根据一些方面的LINAC的控制器。
图19示出了根据一些方面的设备。
具体实施方式
在本申请中,术语“束轴线”可以是表示辐射束中心的三维矢量。
在本申请中,术语“束轴线错位距离”可以是目标(例如,肿瘤中心)与沿着束轴线的任何点之间的最短距离。
在本申请中,术语“CS”或“坐标系”可以是例如二维(2D)CS,诸如图像(例如EPID图像)的CS,或三维(3D)CS,诸如直线加速器(LINAC)(例如LINAC 100)的3D CS。
在本申请中,术语“诊察台走出量”可以是当诊察台106旋转通过其整个运动范围时,患者的肿瘤(如果位于(或接近)等中心)移动的量。
在本申请中,术语“诊察台”可以是支撑患者的LINAC(例如LINAC 100)的部件(例如诊察台106)。
在本申请中,术语“成像装置”(例如,电子射野成像装置(EPID))可以是LINAC(例如,LINAC 100)的部件(例如,成像装置416),其测量从辐射端口(例如,在治疗期间)透射通过患者的x射线强度。例如,成像装置可以将辐射信号电子地转换为二维(2D)数字射线照相图像,以验证与患者的解剖结构相关的正确的束放置。
在本申请中,术语“机架”可以是LINAC(例如LINAC 100)的部件(例如机架102),其在传递辐射束的同时围绕患者旋转。
在本申请中,术语“等中心”可以是使所有机架角度的辐射束轴线错位距离最小化的空间位置。
在本申请中,术语“标记放置误差”可以是在获取用于分析等中心的EPID图像时真实LINAC等中心位置和标记的放置位置之间的三维矢量。
在一些方面,LINAC 100可以(1)利用不同的机架和/或诊察台旋转获取位于(或接近)等中心处的辐射不透过标记414的辐射透射图像,以及(2)分析图像以确定等中心的尺寸、标记放置误差(即,相对于当前标记位置的等中心位置)、和/或诊察台走出量。分析图像可以包括在图像中找到辐射场和标记的中心。
2.1等中心
在一些方面,确定等中心可以涉及分析图像(例如,EPID图像),其中诊察台106保持在固定诊察台角度(例如,0°),同时改变机架角度(例如,通过其整个运动范围)。在一些方面,确定等中心可以涉及仅分析图像,其中诊察台106保持在固定诊察台角度而机架角度变化。一些方面可以包括确定辐射束的束轴线在图像中的位置,以及在图像中找到标记414的中心(参见,例如,图6-8)。如图9和10中所示,一些方面可以包括将束轴线位置和标记中心从二维(2D)图像坐标系(CS)变换为三维(3D)CS。在一些方面,如图9和10所示,因为机架102旋转,同时标记414保持固定在空间中,从2D图像CS到3D(例如,真实世界)CS的变换可以使用标记414的阴影906作为固定参考点(例如,所有的束轴线确定将相对于标记位置,这是空间中的固定点)。
在一些方面,如图9所示,表征等中心可以包括,对于由LINAC 100的成像装置416(例如,EPID)在机架102旋转并且诊察台角度保持固定时获取的每个图像,使用图像处理技术来(i)确定使用LINAC 100的机架102生成的辐射束412的束轴线908的位置,以及(ii)在辐射束412的辐射场902中找到标记414的阴影906的中心。在一些方面,确定束轴线908的位置可以包括找到使用LINAC 100的机架102生成的辐射束412的辐射场902的中心904。在一些替代方面,确定束轴线908的位置可包括找到使用LINAC 100的机架102在(i)第一准直器角度和(ii)与第一准直器角度相差180度的第二准直器角度处生成的辐射场412的中心904的平均值(例如,以在辐射场未被良好校准的情况下考虑误差)。
在一些方面,如图10所示,表征等中心可以包括,对于由成像装置416在机架102旋转并且诊察台角度保持固定时获取的每个图像,构造与确定的图像中的束轴线908的位置一致并且垂直于成像平面的束轴线908。如图10所示,一些方面可以包括,对于由成像装置416在机架102旋转并且诊察台角度保持固定时获取的每个图像(除了以0°的机架角度获取的图像,对于0°的机架角度将不需要旋转),围绕标记414的阴影906的中心1010将构造的束轴线908旋转到新的位置,旋转的量和方向等于在获取图像时使用的机架旋转。如图10所示,一些替代方面可以包括,对于由成像装置416在机架102旋转并且诊察台角度保持固定时获取的每个图像(除了以0°的机架角度获取的图像,对于0°的机架角度将不需要旋转),围绕第一图像旋转轴线旋转第一2D辐射透射图像,该第一图像旋转轴线与确定的第一辐射束的辐射场中的辐射不透过标记的阴影的中心一致并且与机架旋转轴线线平行,第一2D辐射透射图像可以以等于机架旋转角度的量旋转,并且第一3D辐射束轴线可以与确定的旋转的第一2D辐射透射图像中的第一辐射束的束轴线的位置一致,并且垂直于旋转的第一2D辐射透射图像的平面。图10示出了根据一些方面的旋转辐射场1002、旋转辐射场1002的中心1004、旋转辐射场1002中的标记414的阴影1006、以及旋转束轴线1008。
在一些方面,如图11所示,以等于机架旋转角度的量旋转构造的束轴线(或者以等于机架旋转角度的量旋转图像并且使用旋转的图像构造束轴线)可以在辐射束412的辐射场1102中创建标记414的公共阴影1106的一组三维束轴线(例如,包括旋转的束轴线1008、1110和1112,以及以0°的机架旋转角度构造的束轴线1108),其中机架102处于0°的机架旋转角度(以及诊察台106处于固定的诊察台角度)作为参考点(例如,用于计算等中心1114)。一些方面可以包括使用三维束轴线(例如,包括3D束轴线1008、1108、1110和1112)来计算等中心1114。一些方面可以包括通过在空间中找到最小化该点与所有三维束轴线之间的最大束轴线错位距离的点来计算等中心1114。一些方面可以包括将等中心的尺寸确定为等中心1114的中心与所考虑的所有三维束轴线之间的最大错位距离。
2.2标记放置误差
在一些方面,在上面2.1节中描述的三维等中心分析可以具有将任何标记放置误差与所得到的等中心尺寸去耦的额外益处。例如,如图11所示,标记定位误差可以被确定为矢量1116,该矢量1116从(a)在辐射束412的辐射场1102中的标记414的阴影1106的中心,其中机架102处于0°的机架旋转角度(以及诊察台106处于固定的诊察台角度)到(b)找到的等中心1114的位置。
2.3诊察台走出量
一些方面可以包括执行诊察台走出量确定,以确定当诊察台106旋转时放置在等中心1114处的肿瘤将移动多少。在一些方面中,与其中标记414保持固定且辐射场随着机架102旋转而变化的等中心确定相反,LINAC 100可利用固定辐射场(例如,机架102以固定机架角度定位)和(潜在地)随着诊察台106旋转而移动的标记414进行诊察台走出量测量。
在一些方面中,测量诊察台走出量可包括,对于由LINAC 100的成像装置416(例如,EPID)在诊察台106旋转且机架角度保持固定时获取的每个图像,使用图像处理技术来在使用机架102生成的辐射束的辐射场中找到标记414的阴影的中心。例如,如图12所示,其中机架102处于固定机架角度(例如,0°),标记414的阴影1106的中心与诊察台106处于0°的诊察台角度,并且标记414的阴影(例如,阴影1202,1204,1206和1208)的中心可以在诊察台106处于不同的非零诊察台角度的情况下找到。
一些方面可以包括使用标记414的阴影1106在辐射场中的位置,其中诊察台106处于0°诊察台角度作为基线,并且利用非零诊察台旋转从每个图像计算标记414的阴影的移动矢量。图12用虚线示出了从诊察台106处于0°诊察台角度的标记414的阴影1106的中心到诊察台106处于非零诊察台角度的标记414的每个阴影1202,1204,1206和1208的中心的移动矢量。
对于机架旋转和诊察台旋转的不同组合,一些方面可以包括计算在诊察台角度处的标记414和在机架角度处的三维束轴线之间的误差距离。例如,如图12所示,可以计算在特定非零诊察台角度的标记414的阴影1208和在特定非零机架角度的三维束轴线1112之间的标记与束轴线距离1210。
在一些方面,计算对于机架旋转和诊察台旋转的所有组合的标记与束轴线距离可以允许理解对于机架和诊察台角度的所有组合,辐射束将错位肿瘤多少(即使对于诊察台旋转的图像可能已经全部获得,其中机架106固定在一个角度)。
2.4流程图
图13示出了根据一些方面的过程1300。在一些方面,过程1300的一个或多个步骤可以由LINAC 100(例如,LINAC 100的控制器)执行。在一些方面中,过程1300的一个或多个步骤可附加地或替代地由设备(例如,图19的设备1900)执行。在一些方面中,过程1300可包括确定一组三维(3D)辐射束轴线(例如,辐射束轴线1008,1108,1112和1114)的步骤1302。在一些方面,该一组3D辐射束轴线可以包括至少第一和第二3D辐射束轴线1008和1110。在一些方面,过程1300可以包括步骤1304,其基于至少一组3D辐射束轴线来确定辐射等中心1114。
图14示出了根据一些方面的用于确定第一3D辐射束轴线1008的过程1400。在一些方面,过程1400的一个或多个步骤可以由LINAC 100(例如,LINAC 100的控制器)执行。在一些方面中,过程1400的一个或多个步骤可额外地或替代地由设备(例如,图19的设备1900)执行。在一些方面中,过程1400可以在图13所示的过程1300的步骤1302中执行。
在一些方面中,过程1400可包括步骤1402,其将LINAC 100的机架102相对于参考机架角度(例如,0°)以第一机架角度(例如,90°)定位。在一些方面,定位机架102可包括围绕机架旋转轴线208旋转机架102。
在一些方面中,过程1400可包括步骤1404,在机架102以第一机架角度定位的情况下,使用LINAC 100来产生第一辐射束。在一些方面中,过程1400可包括步骤1406,在机架102以第一机架角度定位的情况下,使用LINAC 100的成像装置416来获取第一二维(2D)辐射透射图像,所述第一二维(2D)辐射透射图像指示第一辐射束在通过辐射不透过标记414之后的辐射场902。在一些方面,如图9所示,过程1400可以包括步骤1408,其基于第一2D辐射透射图像来确定第一辐射束的束轴线908的位置和第一辐射束的辐射场902中的辐射不透过标记414的阴影906的中心。
在一些方面,过程1400可以包括步骤1410,其基于确定的第一辐射束的束轴线908的位置、确定的第一辐射束的辐射场902中的辐射不透过标记414的阴影906的中心、以及第一机架角度,来构造第一3D辐射束轴线1008。在一些方面,如图10所示,构造第一3D辐射束轴线1008的步骤1410可以包括围绕第一图像旋转轴线旋转第一2D辐射透射图像,该第一图像旋转轴线与确定的第一辐射束的辐射场902中的辐射不透过标记414的阴影906的中心一致,并且与机架旋转轴线208平行。在一些方面,第一2D辐射透射图像可以以等于第一机架角度的量旋转。在一些方面,如图10所示,第一3D辐射束轴线1008可以与确定的旋转的第一2D辐射透射图像中的第一辐射束的束轴线的位置一致,并且垂直于旋转的2D第一辐射透射图像的平面。
在一些替代方面,如图10所示,构造第一3D辐射束轴线1008的步骤1410可以包括:(i)生成初始3D辐射束轴线908,其与与确定的第一2D辐射透射图像中的第一辐射束的束轴线908的位置一致,并且垂直于第一2D辐射透射图像的平面;以及(ii)围绕第一图像旋转轴线旋转初始3D辐射束轴线,第一图像旋转轴线与确定的第一辐射束的辐射场中的辐射不透过标记414的阴影906的中心一致并且平行于机架旋转轴线。在一些方面,初始3D辐射束轴线可以以等于第一机架角度的量旋转。
在一些方面中,确定第二3D辐射束轴线1110(例如,在过程1300的步骤1302中)可包括将LINAC 100的机架102相对于参考机架角度(例如,0°)以第二机架角度(例如,180°)定位。在一些方面中,确定第二3D辐射束轴线1110可包括,在机架102以第二机架角度定位的情况下,使用LINAC 100来产生第二辐射束。在一些方面中,确定第二3D辐射束轴线1110可包括,在机架102以第二机架角度定位的情况下,使用LINAC 100的成像装置416来获取第二2D辐射透射图像,所述第二2D辐射透射图像指示第二辐射束在通过辐射不透过标记414之后的辐射场。在一些方面,确定第二3D辐射束轴线1110可以包括基于第二2D辐射透射图像来确定第二辐射束的束轴线的位置和第二辐射束的辐射场中的辐射不透过标记414的阴影的中心。在一些方面中,确定第二3D辐射束轴线1110可包括基于确定的第二辐射束的束轴线的位置、确定的辐射不透过标记414的阴影在第二辐射束的辐射场中的中心、以及第二机架角度,来构造第二3D辐射束轴线1110。
在一些方面中,构造第二3D辐射束轴线1110(例如,在过程1300的步骤1302中)可包括围绕第二图像旋转轴线旋转第二2D辐射透射图像,所述第二图像旋转轴线与确定的第二辐射束的辐射场中的辐射不透过标记414的阴影的中心一致并平行于机架旋转轴线208。在一些方面,第二2D辐射透射图像可以以等于第二机架角度的量旋转。在一些方面,第二3D辐射束轴线可以与与确定的旋转的第二2D辐射透射图像中的第二辐射束的束轴线的位置一致,并且垂直于旋转的第二2D辐射透射图像的平面。
在一些方面中,构建第二3D辐射束轴线1110(例如,在过程1300的步骤1302中)可包括(i)产生初始第二3D辐射束轴线,其与与确定的第二2D辐射透射图像中的第二辐射束的束轴线的位置一致,且垂直于第二2D辐射透射图像的平面,和(ii)围绕第二图像旋转轴线旋转初始第二3D辐射束轴线,第二图像旋转轴线与确定的第二辐射束的辐射场中的辐射不透过标记的阴影的中心一致并且与机架旋转轴线平行。在一些方面,初始第二3D辐射束轴线可以以等于第二机架角度的量旋转。
在一些方面,诊察台106可以以固定的诊察台角度(例如,0°)定位,同时LINAC 100用于产生第一和第二辐射束,并且LINAC的成像装置用于获取第一和第二2D辐射透射图像(例如,在过程1300的步骤1302中和/或至少在过程1400的步骤1402,1404和1406中)。
图15示出了根据一些方面的用于确定辐射等中心1114的过程1500。在一些方面,过程1500的一个或多个步骤可以由LINAC 100(例如,LINAC 100的控制器)执行。在一些方面中,过程1500的一个或多个步骤可附加地或替代地由设备(例如,图19的设备1900)执行。在一些方面,过程1500可在图13所示的过程1300的步骤1304中执行。
在一些方面中,过程1500可包括步骤1502,其针对一组3D辐射束轴线(例如,包括3D辐射束轴线1008,1108,1112和1114)中的每个3D辐射束轴线,确定3D空间中的位置与3D辐射束轴线之间的束轴线错位距离。在一些方面,3D空间中的位置与3D辐射束轴线之间的束轴线错位距离可以是该位置与3D辐射束轴线之间的最短距离。在一些方面中,过程1500可包括步骤1504,其确定对于一组3D辐射束轴线的3D辐射束轴线确定的束轴线错位距离中哪一个最大。在一些方面,对于一组3D辐射束轴线中的3D辐射束轴线确定的束轴线错位距离中最大的是该位置的最大束轴线错位距离。在一些方面中,过程1500可包括步骤1506,其在3D空间中找到具有最小的最大束轴线错位距离的位置。在一些方面,辐射等中心1114的尺寸可以是该位置的最大束轴线错位距离。在一些替代方面中,代替找到使最大束轴线距离最小化的位置,可使用其它度量(例如,使束轴线错位距离的平均值最小化,使束轴线错位距离的均方根误差最小化,或其组合)。
在一些方面,如图13所示,过程1300可以包括确定标记放置误差的可选步骤1306。一些方面可以包括将标记放置误差确定为矢量,该矢量从辐射场中的辐射不透过标记414的阴影中心到所确定的辐射等中心1114。在一些方面中,过程1300可包括使用所确定的标记放置误差来将辐射不透过标记414移动到所确定的辐射等中心1114。
在一些方面,如图13所示,过程1300可以包括确定诊察台走出量的可选步骤1308。
在一些方面,如图13所示,过程1300可以包括将肿瘤定位在所确定的辐射等中心1114处的可选步骤1310。
图16示出了根据一些方面的用于确定诊察台走出量的过程1600。在一些方面,过程1600的一个或多个步骤可以由LINAC 100(例如,LINAC 100的控制器)执行。在一些方面中,过程1600的一个或多个步骤可额外地或替代地由设备(例如,图19的设备1900)执行。在一些方面,过程1600可在图13所示的过程1300的步骤1308中执行。
在一些方面中,过程1600可包括步骤1602,其确定一组诊察台角度中的每一诊察台角度的标记移动矢量。在一些方面,该一组诊察台角度可以包括相对于基准诊察台角度(例如,0°)的至少第一诊察台角度和相对于参考诊察台角度的第二诊察台角度。在一些方面,过程1600可包括步骤1604,其基于所确定的标记移动矢量,对于一组诊察台角度的诊察台角度与一组3D辐射束轴线的3D辐射束轴线的每个组合的来确定标记与3D辐射束轴线误差距离(例如,标记与3D辐射束轴线误差距离1210)。在一些方面,标记与3D辐射束轴线误差距离可以表示辐射不透过标记的位置和3D辐射束轴线误差之间的最短距离。
图17示出了用于针对一些方面确定诊察台角度(例如,第一诊察台角度)的标记移动矢量的过程1700。在一些方面,过程1700的一个或多个步骤可以由LINAC 100(例如,LINAC100的控制器)执行。在一些方面中,过程1700的一个或多个步骤可附加地或替代地由设备(例如,图19的设备置1900)执行。在一些方面,过程1700可以在过程1600的步骤1602和/或在图13所示的过程1300的步骤1308中执行。
在一些方面中,过程1700可包括步骤1702,其中LINAC 100的诊察台106以参考诊察台角度定位且LINAC 100的机架102以参考机架角度(例如,0°)定位,使用LINAC 100产生第一参考辐射束。在一些方面中,过程1700可包括步骤1704,在诊察台106以参考诊察台角度定位且机架102以参考机架角度定位的情况下,使用LINAC 100的成像装置416来获取参考诊察台角度二维(2D)辐射透射图像,所述参考诊察台角度二维(2D)辐射透射图像指示第一参考辐射束在通过辐射不透过标记414之后的辐射场1102。在一些方面,过程1700可以包括步骤1706,其基于参考诊察台角度2D辐射透射图像,来确定第一参考辐射束的辐射场1102中的辐射不透过标记414的阴影1106的中心。
在一些方面,过程1700可包括以第一诊察台角度定位诊察台的步骤1708。在一些方面中,过程1700可包括步骤1710,在诊察台106以第一诊察台角度定位且机架102以参考机架角度定位的情况下,使用LINAC 100来产生第二参考辐射束。在一些方面中,过程1700可包括步骤1712:在诊察台106以第一诊察台角度定位且机架102以参考机架角度定位的情况下,使用LINAC 100的成像装置416来获取第一诊察台角度2D辐射透射图像,所述第一诊察台角度2D辐射透射图像指示第二参考辐射束在通过辐射不透过标记之后的辐射场。在一些方面,过程1700可以包括步骤1714,其基于第一诊察台角度2D辐射透射图像来确定第二参考辐射束的辐射场中的辐射不透过标记414的阴影1208的中心。在一些方面,过程1700可以包括步骤1716,其基于第一参考辐射束的辐射场1102中的辐射不透过标记414的阴影1106的中心与第二参考辐射束的辐射场中的辐射不透过标记的阴影1208的中心的比较,来确定第一诊察台角度的标记移动矢量。
在一些方面,在步骤1602中确定第二诊察台角度的标记移动矢量可以包括将诊察台106定位在第二诊察台角度。在一些方面中,在步骤1602中确定第二诊察台角度的标记移动矢量可包括,在诊察台106位于第二诊察台角度且机架102位于参考机架角度的情况下,使用LINAC 100来产生第三参考辐射束。在一些方面,在步骤1602中确定第二诊察台角度的标记移动矢量可以包括,在诊察台106位于第二诊察台角度并且机架102位于参考机架角度的情况下,使用LINAC 100的成像装置416来获取第二诊察台角度2D辐射透射图像,该图像指示在通过辐射不透过标记414之后第三参考辐射束的辐射场。在一些方面,在步骤1602中确定第二诊察台角度的标记移动矢量可以包括基于第二诊察台角度2D辐射透射图像,确定第三参考辐射束的辐射场中的辐射不透过标记414的阴影1206的中心。在一些方面,在步骤1602中确定第二诊察台角度的标记移动矢量可以包括基于第一参考辐射束的辐射场1102中的辐射不透过标记414的阴影1106的中心与第三参考辐射束的辐射场中的辐射不透过标记的阴影1206的中心的比较,来确定第二诊察台角度的标记移动矢量。
在一些方面中,过程1600可进一步包括确定临诊察台等中心。在一些方面,确定临诊察台等中心可以包括确定使最大标记与射束轴线误差距离最小化的空间位置。在一些方面中,确定使最大标记与3D辐射束轴线误差距离最小化的空间位置可包括基于辐射不透过标记在不同诊察台角度下的阴影,来预测诊察台旋转轴线。在一些方面中,确定使最大标记与3D辐射束轴线误差距离最小化的空间位置可包括预测辐射不透过标记414在新的参考位置处的一组位置,该一组位置通过围绕预测的诊察台旋转轴线的诊察台旋转而移位。在一些方面中,确定使最大标记与束轴线误差距离最小化的空间位置可包含基于所确定的标记移动矢量,对一组诊察台角度的诊察台角度与一组3D辐射束轴线的3D辐射束轴线的每个组合来确定预测的标记与3D辐射束轴线误差距离,且预测的标记与3D辐射束轴线误差距离可表示辐射不透过标记的位置与3D辐射束轴线误差之间的最短距离。在一些方面,确定使最大标记与束轴线误差距离最小化的空间位置可以包括确定使最大标记与束轴线误差距离最小化的参考标记位置。在一些方面中,过程1600可进一步包括将肿瘤置于所确定的临诊察台等中心的步骤。
2.5直线加速器控制器
图18是根据一些方面的LINAC 100的控制器1800的框图。如图18所示,控制器1800可以包括:处理电路(PC)1802,其可以包括一个或多个处理器(P)1855(例如,一个或多个通用微处理器和/或一个或多个其他处理器,诸如专用集成电路(ASIC)、现场可编程门阵列(FPGA)等),这些处理器可以共同位于单个外壳中或单个数据中心中,或者可以是地理分布的(即,***可以是分布式计算设备);网络接口1868,包括发射机(Tx)1865和接收机(Rx)1867,用于使控制器1800能够向连接到网络1810(例如,因特网协议(IP)网络)的其他节点发送数据并从其接收数据,网络接口1868连接到该网络1810;机架旋转器1850,其被配置为围绕机架旋转轴线208旋转机架102;诊察台旋转器1852,其被构造成使诊察台106绕诊察台旋转轴线310旋转;辐射束发生器1854,其被配置为在LINAC 100的波导中产生电子束;以及本地存储单元(又名“数据存储***”)1808,其可包括一个或多个非易失性存储装置和/或一个或多个易失性存储装置。在PC 1802包括可编程处理器的方面,可以提供计算机程序产品(CPP)1841。在一些方面,CPP 1841可以包括存储计算机程序(CP)1843的计算机可读介质(CRM)1842,该计算机程序(CP)1843包括计算机可读指令(CRI)1844。在一些方面,CRM 1842可以是非暂时性计算机可读介质,例如磁介质(例如硬盘)、光介质、存储器装置(例如随机存取存储器,闪存)等。在一些方面中,计算机程序1843的CRI 1844可经配置以使得当由PC1802执行时,CRI致使LINAC 100执行本文所描述的步骤(例如,本文参考本文的流程图所描述的一个或多个步骤)。在其它方面中,控制器1800可经配置以在不需要代码的情况下执行本文所述的步骤。即,例如,PC 1802可以仅由一个或多个ASIC组成。因此,这里描述的方面的特征可以用硬件和/或软件来实现。
2.6设备
图19是根据一些方面的设备1900的框图。如图19所示,设备1900可以包括:处理电路(PC)1902,其可以包括一个或多个处理器(P)1955(例如,一个或多个通用微处理器和/或一个或多个其它处理器,诸如专用集成电路(ASIC),现场可编程门阵列(FPGA)等),所述处理器可以共同位于单个外壳中或单个数据中心中,或者可以是地理分布的(即,***可以是分布式计算设备);网络接口1968,包括发射机(Tx)1965和接收机(Rx)1967,用于使设备1900能够向连接到网络1910(例如,因特网协议(IP)网络)的其他节点发送数据并从其接收数据,网络接口1968连接到该网络1910;和本地存储单元(又名“数据存储***”)1908,其可包括一个或多个非易失性存储装置和/或一个或多个易失性存储装置。在PC 1902包括可编程处理器的方面,可以提供计算机程序产品(CPP)1941。在一些方面,CPP 1941可包括存储计算机程序(CP)1943的计算机可读介质(CRM)1942,该计算机程序(CP)1943包括计算机可读指令(CRI)1944。在一些方面,CRM 1942可以是非暂时性计算机可读介质,例如磁介质(例如硬盘)、光介质、存储器装置(例如随机存取存储器,闪存)等。在一些方面中,计算机程序1943的CRI 1944可经配置以使得当由PC 1902执行时,CRI致使设备1900执行本文描述的步骤(例如,本文参考本文的流程图描述的一个或多个步骤)。在其它方面中,设备1900可经配置以在不需要代码的情况下执行本文所述的步骤。即,例如,PC 1902可以仅由一个或多个ASIC组成。因此,这里描述的方面的特征可以用硬件和/或软件来实现。
虽然在这里描述了各种实施例,但是应当理解,它们仅以示例的方式给出,而不是限制。因此,本公开的广度和范围不应受上述示例性实施例中的任何一个的限制。此外,除非本文另有说明或与上下文明显矛盾,否则上述要素在其所有可能变化中的任何组合都被本公开涵盖。
另外,虽然上述和附图中所示的过程被示为一系列步骤,但是这仅仅是为了说明的目的而完成的。因此,可以设想可以添加一些步骤,可以省略一些步骤,可以重新安排步骤的顺序,并且可以并行地执行一些步骤。
Claims (19)
1.一种方法,其特征在于,包括:
根据二维(2D)辐射透射图像,确定直线加速器(LINAC)的一组三维(3D)辐射束轴线,其中,所述一组3D辐射束轴线包括至少第一3D辐射束轴线和第二3D辐射束轴线;和
至少基于所述一组3D辐射束轴线,确定所述LINAC的辐射等中心。
2.根据权利要求1所述的方法,其特征在于,其中确定所述第一3D辐射束轴线包括:
相对于参考机架角度以第一机架角度定位所述LINAC的机架,其中定位所述机架包括使所述机架围绕机架旋转轴线旋转;
在以所述第一机架角度定位所述机架的情况下,使用所述LINAC生成第一辐射束;
在以所述第一机架角度定位所述机架的情况下,使用所述LINAC的成像装置来获取第一二维(2D)辐射透射图像,所述第一2D辐射透射图像指示所述第一辐射束经过辐射不透过标记之后的辐射场;
基于所述第一2D辐射透射图像,确定所述第一辐射束的束轴线和所述辐射不透过标记在所述第一辐射束的所述辐射场中的阴影的中心的位置;和
基于确定的所述第一辐射束的所述束轴线的位置、确定的所述第一辐射束的所述辐射场中的所述辐射不透过标记的所述阴影的所述中心以及所述第一机架角度,构造所述第一3D辐射束轴线。
3.根据权利要求2所述的方法,其特征在于,其中:
构造所述第一3D辐射束轴线包括,使所述第一2D辐射透射图像围绕第一图像旋转轴线旋转,所述第一图像旋转轴线与确定的所述第一辐射束的所述辐射场中的所述辐射不透过标记的所述阴影的所述中心一致并且与所述机架旋转轴线平行;
所述第一2D辐射透射图像旋转与所述第一机架角度相等的量;并且
所述第一3D辐射束轴线与确定的旋转的所述第一2D辐射透射图像中的所述第一辐射束的所述束轴线的所述位置一致并且与旋转的所述第一2D辐射透射图像的平面垂直。
4.根据权利要求2所述的方法,其特征在于,其中,构造所述第一3D辐射束轴线包括:
生成初始3D辐射束轴线,所述初始3D辐射束轴线与确定的所述第一2D辐射透射图像中的所述第一辐射束的所述束轴线的所述位置一致并且与所述第一2D辐射透射图像的平面垂直;和
使所述初始3D辐射束轴线围绕第一图像旋转轴线旋转,所述第一图像旋转轴线与确定的所述第一辐射束的所述辐射场中的所述辐射不透过标记的所述阴影的所述中心一致并且与所述机架旋转轴线平行;
其中,所述初始3D辐射束轴线旋转与所述第一机架角度相等的量。
5.根据权利要求2-4中任一项所述的方法,其特征在于,其中,确定所述第二3D辐射束轴线包括:
相对于所述参考机架角度以第二机架角度定位所述LINAC的所述机架;
在以所述第二机架角度定位所述机架的情况下,使用LINAC生成第二辐射束;
在以所述第二机架角度定位所述机架的情况下,使用LINAC的所述成像装置获取第二2D辐射透射图像,所述第二2D辐射透射图像指示所述第二辐射束经过所述辐射不透过标记之后的辐射场;
基于所述第二2D辐射透射图像,确定所述第二辐射束的束轴线的位置和所述第二辐射束的所述辐射场中的所述辐射不透过标记的阴影的中心;和
基于确定的所述第二辐射束的所述束轴线的位置、确定的所述第二辐射束的所述辐射场中的所述辐射不透过标记的所述阴影的所述中心以及所述第二机架角度,构造所述第二3D辐射束轴线。
6.根据权利要求5的方法,其特征在于,其中:
构造所述第二3D辐射束轴线包括,使所述第二2D辐射透射图像围绕第二图像旋转轴线旋转,所述第二图像旋转轴线与确定的所述第二辐射束的所述辐射场中的所述辐射不透过标记的所述阴影的所述中心一致并且与所述机架旋转轴线平行;
所述第二2D辐射透射图像旋转与所述第二机架角度相等的量;并且
所述第二3D辐射束轴线与确定的旋转的所述第二2D辐射透射图像中的所述第二辐射束的所述束轴线的所述位置一致并且与旋转的所述第二2D辐射透射图像的平面垂直。
7.根据权利要求5所述的方法,其特征在于,其中,构造所述第二3D辐射束轴线包括:
生成初始第二3D辐射束轴线,所述初始第二3D辐射束轴线与确定的所述第二2D辐射透射图像中的所述第二辐射束的所述束轴线的所述位置一致并且与所述第二2D辐射透射图像的平面垂直;和
使所述初始第二3D辐射束轴线围绕第二图像旋转轴线旋转,所述第二图像旋转轴线与确定的所述第二辐射束的所述辐射场中的所述辐射不透过标记的所述阴影的所述中心一致并且与所述机架旋转轴线平行;
其中,所述初始第二3D辐射束轴线旋转与所述第二机架角度相等的量。
8.根据权利要求5-7中任一项所述的方法,其特征在于,其中,在所述LINAC用于生成所述第一辐射束和所述第二辐射束并且所述LINAC的所述成像装置用于获取所述第一2D辐射透射图像和所述第二2D辐射透射图像的同时,以固定诊察台角度定位诊察台。
9.根据权利要求1-8中任一项所述的方法,其特征在于,其中,确定所述辐射等中心包括:
对于所述一组3D辐射束轴线中的每个3D辐射束轴线,确定3D空间中的位置和所述3D辐射束轴线之间的束轴线错位距离;
确定对于所述一组3D辐射束轴线中的3D辐射束轴线确定的所述束轴线错位距离中哪一个最大,其中对于所述一组3D辐射束轴线中的3D辐射束轴线确定的所述束轴线错位距离中最大的是所述位置的最大束轴线错位距离;和
在3D空间中找到具有最小的最大束轴线错位距离的位置。
10.根据权利要求9所述的方法,其特征在于,其中,所述辐射等中心的尺寸是所述位置的最大束轴线错位距离。
11.根据权利要求1-10中任一项所述的方法,其特征在于,进一步包括将肿瘤定位在确定的辐射等中心处。
12.根据权利要求1-11中任一项所述的方法,其特征在于,进一步包括:
对于一组诊察台角度中的每个诊察台角度确定标记移动矢量,其中,所述一组诊察台角度至少包括相对于参考诊察台角度的第一诊察台角度和相对于所述参考诊察台角度的第二诊察台角度;和
基于确定的所述标记移动矢量,对于所述一组诊察台角度中的诊察台角度与所述一组3D辐射束轴线中的3D辐射束轴线的每个组合确定标记与3D辐射束轴线误差距离,其中所述标记与3D辐射束轴线误差距离表示辐射不透过标记的位置和所述3D辐射束轴线误差之间的最短距离。
13.根据权利要求12所述的方法,其特征在于,其中,确定用于所述第一诊察台角度的所述标记移动矢量包括:
在以所述参考诊察台角度定位所述LINAC的诊察台并且以所述参考机架角度定位所述LINAC的机架的情况下,使用所述LINAC生成第一参考辐射束;
在以所述参考诊察台角度定位所述诊察台并且以所述参考机架角度定位所述机架的情况下,使用所述LINAC的成像装置获取参考诊察台角度二维(2D)辐射透射图像,所述参考诊察台角度二维(2D)辐射透射图像指示所述第一参考辐射束经过所述辐射不透过标记之后的辐射场;
基于所述参考诊察台角度2D辐射透射图像,确定所述第一参考辐射束的所述辐射场中的所述辐射不透过标记的阴影的中心;
以所述第一诊察台角度定位所述诊察台;
在以所述第一诊察台角度定位所述诊察台并且以所述参考机架角度定位所述机架的情况下,使用所述LINAC生成第二参考辐射束;
在以所述第一诊察台角度定位所述诊察台并且以所述参考机架角度定位所述机架的情况下,使用所述LINAC的所述成像装置获取第一诊察台角度2D辐射透射图像,所述第一诊察台角度2D辐射透射图像指示所述第二参考辐射束经过所述辐射不透过标记之后的辐射场;
基于所述第一诊察台角度2D辐射透射图像,确定所述第二参考辐射束的所述辐射场中的所述辐射不透过标记的阴影的中心;和
基于所述第一参考辐射束的所述辐射场中的所述辐射不透过标记的所述阴影的所述中心与所述第二参考辐射束的所述辐射场中的所述辐射不透过标记的所述阴影的所述中心的比较,确定所述第一诊察台角度的所述标记移动矢量。
14.根据权利要求13所述的方法,其特征在于,其中确定用于所述第二诊察台角度的所述标记移动矢量包括:
以所述第二诊察台角度定位所述诊察台;
在以所述第二诊察台角度定位所述诊察台并且以所述参考机架角度定位所述机架的情况下,使用所述LINAC生成第三参考辐射束;
在以所述第二诊察台角度定位所述诊察台并且以所述参考机架角度定位所述机架的情况下,使用所述LINAC的所述成像装置获取第二诊察台角度2D辐射透射图像,所述第二诊察台角度2D辐射透射图像指示所述第三参考辐射束经过所述辐射不透过标记之后的辐射场;
基于所述第二诊察台角度2D辐射透射图像,确定所述第三参考辐射束的所述辐射场中的所述辐射不透过标记的阴影的中心;和
基于所述第一参考辐射束的所述辐射场中的所述辐射不透过标记的所述阴影的所述中心与所述第三参考辐射束的所述辐射场中的所述辐射不透过标记的所述阴影的所述中心的比较,确定所述第二诊察台角度的所述标记移动矢量。
15.根据权利要求12-14中任一项所述的方法,其特征在于,进一步包括确定临床等中心,其中,确定所述临床等中心包括确定使最大的标记与3D辐射束轴线误差距离最小化的空间位置。
16.根据权利要求15所述的方法,其特征在于,其中,确定使所述最大的标记与3D辐射束轴线误差距离最小化的所述空间位置包括:
基于不同诊察台角度下的所述辐射不透过标记的阴影,预测诊察台旋转轴线;
预测所述辐射不透过标记在新的参考位置处的一组位置,所述一组位置通过围绕预测的所述诊察台旋转轴线的诊察台旋转而移位;
基于确定的所述标记移动矢量,对于所述一组诊察台角度中的诊察台角度和所述一组3D辐射束轴线中的3D辐射束轴线的每个组合确定预测的标记与3D辐射束轴线误差距离,其中预测的所述标记与3D辐射束轴线误差距离表示所述辐射不透过标记的位置和所述3D辐射束轴线误差之间的最短距离;和
确定使所述最大的标记与束轴线误差距离最小化的参考标记位置。
17.根据权利要求15或16所述的方法,其特征在于,进一步包括将肿瘤置于确定的所述临床等中心处。
18.一种设备,其特征在于,所述设备被配置为:
根据二维(2D)辐射透射图像,确定直线加速器(LINAC)的一组三维(3D)辐射束轴线,其中所述一组3D辐射束轴线包括至少第一3D辐射束轴线和第二3D辐射束轴线;和
至少基于所述一组3D辐射束轴线,确定所述LINAC的辐射等中心。
19.一种方法,包括:
对于直线加速器(LINAC)的一组诊察台角度中的每个诊察台角度确定标记移动矢量,其中,所述一组诊察台角度至少包括相对于参考诊察台角度的第一诊察台角度和相对于所述参考诊察台角度的第二诊察台角度;和
基于确定的所述标记移动矢量,对于所述一组诊察台角度中的诊察台角度与所述一组3D辐射束轴线中的3D辐射束轴线的每个组合确定标记与3D辐射束轴线误差距离,其中所述标记与3D辐射束轴线误差距离表示辐射不透过标记的位置和所述3D辐射束轴线误差之间的最短距离。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US202163176494P | 2021-04-19 | 2021-04-19 | |
US63/176,494 | 2021-04-19 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115212475A true CN115212475A (zh) | 2022-10-21 |
Family
ID=81328008
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210412279.9A Pending CN115212475A (zh) | 2021-04-19 | 2022-04-19 | 用于评价直线加速器等中心质量的x射线透射图像分析 |
Country Status (6)
Country | Link |
---|---|
US (1) | US20220331608A1 (zh) |
EP (2) | EP4079372A3 (zh) |
JP (1) | JP2022165419A (zh) |
CN (1) | CN115212475A (zh) |
AU (1) | AU2022202565A1 (zh) |
GB (1) | GB2607702A (zh) |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7142634B2 (en) * | 2003-01-29 | 2006-11-28 | New England Medical Center Hospitals, Inc. | Radiation field detection |
EP2186542A1 (en) * | 2008-11-14 | 2010-05-19 | Ion Beam Applications S.A. | Method and device for determining a radiation isocenter in a particle therapy room |
CN102414759B (zh) * | 2009-06-09 | 2014-02-19 | 三菱电机株式会社 | 粒子射线照射装置 |
EP2450083B1 (fr) * | 2010-11-03 | 2013-08-07 | Qualiformed | Objet-test pour le contrôle qualité d'un appareil de traitement par radiothérapie et procédés de fabrication et d'utilisation dudit objet-test |
US9142016B2 (en) * | 2013-09-27 | 2015-09-22 | Varian Medical Systems, Inc. | Methods for image processing |
US9192784B1 (en) * | 2015-06-16 | 2015-11-24 | Radiological Imaging Technology, Inc. | Radiation delivery isocenter analysis tool |
US10682528B2 (en) * | 2017-03-03 | 2020-06-16 | Varian Medical Systems International Ag | Systems, methods, and devices for radiation beam asymmetry measurements using electronic portal imaging devices |
GB2589929B (en) * | 2019-12-13 | 2021-12-22 | Elekta ltd | Calibration of radiotherapy apparatus |
US11311747B2 (en) * | 2020-07-16 | 2022-04-26 | Uih America, Inc. | Systems and methods for isocenter calibration |
US11617902B2 (en) * | 2020-08-26 | 2023-04-04 | Aktina Corp. | Modulated radiation beam alignment for medical linear accelerator |
-
2022
- 2022-04-19 EP EP22168855.9A patent/EP4079372A3/en active Pending
- 2022-04-19 US US17/724,027 patent/US20220331608A1/en active Pending
- 2022-04-19 GB GB2205695.6A patent/GB2607702A/en active Pending
- 2022-04-19 EP EP23212978.3A patent/EP4316582A2/en active Pending
- 2022-04-19 AU AU2022202565A patent/AU2022202565A1/en active Pending
- 2022-04-19 JP JP2022068870A patent/JP2022165419A/ja active Pending
- 2022-04-19 CN CN202210412279.9A patent/CN115212475A/zh active Pending
Also Published As
Publication number | Publication date |
---|---|
EP4079372A3 (en) | 2023-01-11 |
GB202205695D0 (en) | 2022-06-01 |
GB2607702A (en) | 2022-12-14 |
AU2022202565A1 (en) | 2022-11-03 |
JP2022165419A (ja) | 2022-10-31 |
US20220331608A1 (en) | 2022-10-20 |
EP4079372A2 (en) | 2022-10-26 |
EP4316582A2 (en) | 2024-02-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11633629B2 (en) | Method of calibration of a stereoscopic camera system for use with a radio therapy treatment apparatus | |
EP3698848A1 (en) | Magnetic resonance projection imaging | |
US8471222B2 (en) | Radiotherapy apparatus control method and radiotherapy apparatus control apparatus | |
CN111132730B (zh) | 与放射治疗设备一起使用的患者监测***的校准方法 | |
JP6565080B2 (ja) | 放射線治療装置、その作動方法及びプログラム | |
JP6009705B1 (ja) | 放射線照射アイソセンタ分析器具 | |
CN109414235B (zh) | 使用辐射成像的体模设置和源到表面距离验证 | |
CN110075428A (zh) | 一种射束检验、测量方法及装置 | |
CN111432878B (zh) | 放射治疗装置的等中心的基于epid的验证、修正和最小化的方法 | |
CN103028195A (zh) | 用于辐射治疗计划的组合成像模式 | |
US20210052919A1 (en) | Techniques for skin illumination of radiation therapy treatment region using non-coaxial digital projector | |
US20210137485A1 (en) | Treatment system, calibration method, and storage medium | |
JP2017169627A (ja) | X線撮影機器のアライメント調整支援装置、方法及びプログラム | |
CN115212475A (zh) | 用于评价直线加速器等中心质量的x射线透射图像分析 | |
Yuan et al. | A Monte Carlo model and its commissioning for the Leksell Gamma Knife Perfexion radiosurgery system | |
EP4201333A1 (en) | Universal phantom for calibration and verification of optical and radiation systems | |
US20240173573A1 (en) | Method of identifying an isocenter | |
CN115400358A (zh) | 医用直线加速器的辐射束对准 | |
CN115965741A (zh) | 透视成像设备的控制方法、装置、电子设备和存储介质 | |
Torfeh et al. | 3D numerical test objects for the evaluation of a software used for an automatic analysis of a linear accelerator mechanical stability |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: DE Ref document number: 40081629 Country of ref document: HK |
|
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |