CN107530552B - 用于自适应放射治疗的运动靶的三维定位 - Google Patents

用于自适应放射治疗的运动靶的三维定位 Download PDF

Info

Publication number
CN107530552B
CN107530552B CN201680013998.6A CN201680013998A CN107530552B CN 107530552 B CN107530552 B CN 107530552B CN 201680013998 A CN201680013998 A CN 201680013998A CN 107530552 B CN107530552 B CN 107530552B
Authority
CN
China
Prior art keywords
target
slice
processor
adaptive filter
patch
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
Application number
CN201680013998.6A
Other languages
English (en)
Other versions
CN107530552A (zh
Inventor
韩骁
周燕
Original Assignee
医科达有限公司
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 医科达有限公司 filed Critical 医科达有限公司
Publication of CN107530552A publication Critical patent/CN107530552A/zh
Application granted granted Critical
Publication of CN107530552B publication Critical patent/CN107530552B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/103Treatment planning systems
    • A61N5/1037Treatment planning systems taking into account the movement of the target, e.g. 4D-image based planning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/20Analysis of motion
    • G06T7/277Analysis of motion involving stochastic approaches, e.g. using Kalman filters
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/70Determining position or orientation of objects or cameras
    • G06T7/73Determining position or orientation of objects or cameras using feature-based methods
    • G06T7/74Determining position or orientation of objects or cameras using feature-based methods involving reference images or patches
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/103Treatment planning systems
    • A61N5/1039Treatment planning systems using functional images, e.g. PET or MRI
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00Image generation
    • G06T2211/40Computed tomography
    • G06T2211/428Real-time

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Biomedical Technology (AREA)
  • Radiology & Medical Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Public Health (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Multimedia (AREA)
  • Medical Informatics (AREA)
  • Quality & Reliability (AREA)
  • Radiation-Therapy Devices (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本公开涉及用于分割医学图像的系统、方法以及计算机可读存储介质。本公开的实施例可以在三维(3D)体积中对靶进行定位。例如,图像采集装置可以提供包含靶的关注区域的3D医学图像。然后,处理器可以从3D图像中提取多个二维(2D)切片。处理器还可以为每个2D切片确定2D图像块,其中所述2D图像块对应于与所述靶相关联的2D切片的区域。处理器还可以将2D图像块转换到自适应滤波器模型以确定所述关注区域的位置。

Description

用于自适应放射治疗的运动靶的三维定位
相关专利申请的交叉引用
本专利申请与2015年1月28日提交的题目为“Three Dimensional Localizationand Tracking for Adaptive Radiation Therapy”(“用于自适应放射治疗的三维定位及追踪”)且代理人案件号为12475.0059-00000的案件相关,该案件的全部内容通过引用并入本文。
技术领域
本公开总体上涉及放射治疗或放射疗法。更具体地,本公开涉及在制定在放射疗法期间要使用的放射疗法治疗计划时使用的用于肿瘤的三维追踪的系统和方法。
背景技术
放射疗法用于治疗哺乳动物(例如人和动物)组织中的癌症或其它疾病。一种这样的放射疗法技术是“伽马刀”,通过“伽马刀”,患者被以较高强度和高精度汇聚到靶区(例如肿瘤)处的大量较低强度的γ射线照射。在另一个实施例中,使用直线加速器提供放射疗法,由此靶区被高能粒子(例如电子、质子、离子等)照射。必须准确地控制放射束的放置和剂量以确保肿瘤接受处方的放射,并且放射束的放置应最小化对周围健康组织(通常被称为“危及器官”(OAR))的损害。
例如,放射束可以通过使用多叶准直器而被成形为匹配肿瘤的形状(例如,多叶准直器可以包括可以彼此独立地运动以生成定制的放射束形状的多个钨叶)。(放射被称为“处方的”,因为医师为肿瘤和周围的器官开了类似于药物处方的预定义的放射量)。
传统上,对于每个患者,可以使用基于临床和剂量测定目标和约束(例如,对肿瘤和关键器官的最大、最小和平均放射剂量)的优化技术来创建放射疗法治疗计划(“治疗计划”)。治疗计划过程可以包括使用患者的三维图像来识别靶区(例如,肿瘤),并且识别肿瘤附近的关键器官。治疗计划的创建可能是一个耗时的过程,其中,计划者试图遵守各种治疗目的或约束(例如,剂量体积直方图(DVH)目的),考虑它们各自的重要性(例如,权重),以生成临床可接受的治疗计划。这个任务可能是由于各种危及器官(OAR)而复杂化的耗时的试错过程,因为随着OAR数量增大(例如,对于头颈部治疗达到十三处),该过程的复杂性也增大。远离肿瘤的OAR可以更容易地免于放射,但靠近靶肿瘤或与靶肿瘤重叠的OAR可能更难以免于放射。
计算机断层摄影(CT)成像传统上用作用于确定放射治疗的治疗计划的图像数据的主要来源。CT图像提供患者几何形状的精确表示,并且CT值可以直接转换为电子密度(例如,亨斯菲尔德单位)用于放射剂量计算。然而,使用CT使患者暴露于额外的放射剂量。除CT图像之外,核磁共振成像(MRI)扫描可以被用于放射治疗,由于其相比于CT图像的更好的软组织对比度,MRI没有电离放射,并且可以被用于捕捉诸如组织代谢和功能的人体的功能信息。
可以使用诸如计算机断层摄影(CT)、荧光透视、以及核磁共振成像(“MRI”)的成像系统来确定靶的定位和追踪靶(例如,器官、肿瘤等)。由于MRI可以在不使用如CT所使用的电离放射的情况下提供极好的软组织对比度,因而通常使用MRI。与成像系统集成的放射疗法治疗系统的示例可以包括MRI-Linac,其可以使用靶(例如肿瘤)的三维(3D)图像。MRI-LINAC(直线加速器)的MRI装置可以提供对应于患者的组织中氢核的局部图的多个图像。该患者的图像可以在二维(2D)平面或三维体积中获取。由于器官和肿瘤在患者的身体内移动,因此对靶进行快速和准确的3D定位是很重要的。例如,靶器官或肿瘤可能由于多种类型的运动(例如,呼吸、心脏的、蠕动或其它类型的患者运动)而移动。然而,2D MR切片通常是在患者的身体的特定位置获取的,并且由于靶器官或肿瘤的运动,2D MR切片可能不包括肿瘤或靶。因此,需要一种获取靶肿瘤并且追踪该肿瘤的结构和运动的系统和方法,使得所获取的2D MR切片就位置、取向和/或厚度而言包含靶肿瘤并且是在临床环境中在2D MR切片上是可见的。
本概述是本申请的一些教导的概述,并不旨在作为本主题的排他或穷尽的治疗。关于本主题的进一步细节在具体实施方式和所附权利要求中找到。在阅读并理解以下具体实施方式并参见形成具体实施方式的一部分的附图之后,本发明的其它方案对于本领域技术人员将是显而易见的,其中每一个不被视为限制性的。本发明的范围由所附权利要求及其依照法律的等同物限定。
发明内容
本公开的某些实施例涉及一种用于医疗图像处理的方法。所述方法可以由执行多个计算机可执行指令的处理器装置来实现。所述可以由处理器来实现的方法用于在三维(3D)体积中定位靶,并且所述方法包括:从图像采集装置接收包含所述靶的关注区域的3D医学图像;由所述处理器从所述3D图像中提取多个二维(2D)切片;由所述处理器为每个2D切片确定2D图像块,其中所述2D图像块对应于与所述靶相关联的所述2D切片的区域;以及由所述处理器将所述2D图像块转换到自适应滤波器模型,用于确定所述关注区域的位置。
本公开的某些实施例涉及一种医学图像处理系统。所述系统可以使用体积的二维成像切片追踪所述体积中的三维(3D)靶,并且所述系统包括:处理器;以及存储器,其可操作地联接到所述处理器并且存储计算机可执行指令,当所述计算机可执行指令由所述处理器执行时,使所述处理器执行方法,所述方法包括:从图像采集装置接收包含所述靶的关注区域的所述3D体积的3D医学图像,其中所述3D图像被存储在所述存储器中;由所述处理器从所述3D图像中提取多个二维(2D)切片;由所述处理器为每个2D切片确定2D图像块,其中所述2D图像块对应于与所述靶相关联的所述2D切片的区域;以及由所述处理器将所述2D图像块转换到自适应滤波器模型,用于确定所述关注区域的位置。
本公开的其它目的和优点将在下面的具体实施例中进行部分地阐述,并且将部分地从本说明书中显而易见,或者可以通过本公开的实践来获知。本公开的目的和优点可以通过所附权利要求中特别指出的元件和组合来实现和获得。
应当理解的是,以上的概述以及下面的具体实施方式仅仅是示例性和说明性的,而不是对所要求保护的本发明的限制。
附图说明
在不一定按比例绘制的附图中,相同的附图标记可以在不同的视图中描述相似的构件。具有不同字母后缀的相同附图标记可以表示相似构件的不同实例。附图大体上通过举例而不是限制的方式示出各种实施例,并且与说明书以及权利要求书一起用于对所公开的实施例进行说明。在适当的时候,在所有附图中使用相同的附图标记指代同一或相似的部分。这样的实施例是示例性的,并且不旨在作为本装置、系统或方法的穷尽性的或排他性的实施例。
图1示出了对用于放射治疗的三维靶进行定位和追踪的示例性放射治疗系统;
图2示出了图1的放射治疗系统中使用的放射治疗装置伽玛刀;
图3示出了图1的放射治疗系统中使用的放射治疗装置直线加速器;
图4示出了用于训练模块以生成自适应滤波器模型的示例性流程图;
图5示出了在患者的放射治疗期间使用一个或更多训练好的自适应滤波器模型来定位和追踪肿瘤的示例性流程图;
图6为自适应滤波器模型、响应图、以及追踪患者中的靶的图像的图示。
具体实施方式
在本公开中,介绍了用于对经受放射疗法治疗的患者的解剖结构中的肿瘤进行定位和追踪的放射治疗系统和相应的方法。具体地,在放射疗法治疗计划期间,捕获多个3D核磁共振图像(MRI)或4D MRI图像。可以通过放射线治疗系统使用本领域中公知的分割技术来确定靶(例如患者的解剖结构中的肿瘤)的确切位置。在肿瘤的位置被确定之后,系统可以生成多个图像块,并且可以确定从图像块的中心到肿瘤的偏移。放射治疗系统可以使用诸如相关性或使用外观模型的技术的多种方法来增强多个图像块的选择特征,其导致多个自适应滤波器模型。接着,在患者的实际“在线”治疗期间放射治疗系统可以使用这些自适应滤波器模型。
这种方法的一个优点是在放射疗法治疗期间可以对靶肿瘤进行定位和追踪。由于诸如呼吸、心脏运动、蠕动引起的运动、患者的不自主运动(例如,咳嗽、打喷嚏等)、或由放射治疗期间治疗台上的患者的自主运动的多种类型的运动,肿瘤通常在患者的解剖结构内移动。当使用MRI引导的放射治疗时,只有患者的解剖结构的图像的2D切片是可用的。这种方法允许确定哪些(一个或数个)2D切片包括肿瘤,并且进一步有利地提供了通过估计可能的未来位置来追踪肿瘤的能力。
图1示出用于在放射疗法治疗期间执行靶定位和追踪的示例性放射治疗系统100。放射治疗系统100可以包括连接到网络120的放射治疗装置110,网络120连接到因特网130。网络120可以将放射治疗装置110与数据库140、医院数据库142、肿瘤信息系统(OIS)150(例如,其可提供患者信息)、治疗计划系统(TPS)160(例如,用于生成待由所述放射治疗装置110使用的放射疗法治疗计划)、图像采集装置170、显示装置180和用户接口190连接。
放射治疗装置110可以包括处理器112、存储装置116、以及通信接口114。存储装置116可以存储用于操作系统118、治疗计划软件120、生成自适应滤波器126的训练模块124、以及靶定位模块120的计算机可执行指令,以及由处理器240执行的任何其他计算机可执行指令。
处理器112可以通信地联接到存储装置116,并且处理器112可以被配置为执行存储于其上的计算机可执行指令,例如,处理器112可以在放射治疗的执行期间执行训练模块124以实现训练模块124的功能以及靶定位模块128的功能,以便确定患者内的靶的定位。此外,处理器112可以执行治疗计划软件120(例如,诸如由Elekta公司制造的
Figure BDA0001399441100000051
软件),其可与训练模块124以及靶定位模块128交互。
处理器112可以是处理装置,包括例如微处理器、中央处理单元(CPU)、图形处理单元(GPU)、加速处理单元(APU)等的一个或更多通用处理装置。更具体地,处理器112可以是复杂指令集计算(CISC)微处理器、精简指令集计算(RISC)微处理器、超长指令字(VLIW)微处理器、实现其它指令集的处理器、或实现指令集的组合的处理器。处理器112还可以是一个或更多专用处理装置,例如专用集成电路(ASIC)、现场可编程门阵列(FPGA)、数字信号处理器(DSP)、片上系统(SOC)等。如本领域技术人员将领会的,在一些实施例中,处理器112可以是专用处理器而不是通用处理器。处理器112可以包括一个或更多已知的处理装置,例如来自由IntelTM制造的PentiumTM、CoreTM、XeonTM、或
Figure BDA0001399441100000052
系列的微处理器,由AMDTM制造的TurionTM、AthlonTM、SempronTM、OpteronTM、FXTM、PhenomTM系列的微处理器,或由太阳微系统公司制造的任何各种处理器。处理器112还可以包括图形处理单元,例如来自由NvidiaTM制造的
Figure BDA0001399441100000053
系列的图形处理单元,由IntelTM制造的GMA、IrisTM系列的图形处理单元,或由AMDTM制造的RadeonTM系列的图形处理单元。处理器112还可以包括加速处理单元,例如由AMDTM制造的台式A-4(6,8)系列TM,IntelTM制造的Xeon PhiTM系列。公开的实施例不限于任何类型的(一个或数个)处理器,其以其他方式被配置为满足识别、分析、保持、生成和/或提供大量成像数据或处理这样的成像数据以对靶进行追踪或定位的计算需求,或以其他方式被配置为处理与所公开的实施例一致的任何其它类型的数据。此外,术语“处理器”可包括多于一个处理器,例如,多核设计或每一个都具有多核设计的多个处理器。处理器112可以执行存储在存储器116的计算机程序指令序列以执行将在下面更详细地说明的各种操作、处理、方法。
存储装置116可存储从图像采集装置179接收的图像数据122(例如,3D MRI、4DMRI、2D切片等)、或放射治疗装置110可以使用以执行与所公开的实施例一致的操作的任何格式的任何其它类型的数据/信息。存储装置210可包括只读存储器(ROM)、闪存、随机存取存储器(RAM)、诸如同步DRAM(SDRAM)或Rambus DRAM的动态随机存取存储器(DRAM)、静态存储器(例如,闪速存储器、静态随机存取存储器)等,在其上存储有任何格式的计算机可执行指令。该计算机程序指令可以由处理器112访问,从ROM或任何其它合适的存储位置读取,加载到RAM以供处理器112执行。例如,存储器116可以存储一个或更多软件应用程序。例如,存储在存储器116的软件应用程序可以包括用于常见计算机系统以及用于软件控制的装置的操作系统118。此外,存储器116可以存储整个软件应用程序或可由处理器112执行的仅一部分软件应用程序。例如,存储装置116可以存储由治疗计划系统160生成的作为治疗计划软件120的一个或更多放射疗法治疗计划。
在一些实施例中,存储装置116可以包括机器可读的存储介质,而在实施例中,机器可读存储介质可以是单一介质,但是术语“机器可读存储介质”应当被理解为包括的单一介质或多介质(例如,集中式或分布式数据库,和/或相关的缓存和服务器),其存储一个或更多的计算机可执行指令集或数据集。术语“机器可读存储介质”还应当被理解为包括能够对由机器执行并且使得机器执行本公开的方法中的任何一个或更多的指令集进行存储或编码的任何介质。因此,术语“机器可读存储介质”应当被理解为包括但不限于固态存储器、光学介质和磁性介质。例如,存储器/存储装置210可以是一个或更多的易失性的、非暂时性的、或非易失性有形计算机可读介质。
放射治疗装置110可以经由通信接口114与网络130通信,所述通信接口114通信地联接到处理器112和存储装置116。例如,通信接口114可以包括网络适配器、电缆连接器、串行连接器、USB连接器、并行连接器、高速数据传输适配器(例如,光纤适配器、USB 3.0、雷电等)、无线网络适配器(例如,WiFi适配器)、电信适配器(例如,3G,4G/LTE等)等。通信接口114可以包括一个或更多允许放射治疗装置110经由网络130与其他机器和装置(如远程组件)通信的数字和/或模拟通信装置。
网络130可以提供局域网(LAN)、无线网络、云计算环境(例如,服务软件、服务平台、服务基础设施等)、客户端-服务器、广域网(WAN)等的功能。因此,网络130可以允许放射治疗装置110与若干各种其它系统和装置之间的数据传输,例如:治疗计划系统160,肿瘤信息系统150和图像采集装置170。此外,由治疗计划系统160、OIS 150和图像采集装置170生成的数据可以被存储在存储器116、数据库140、或医院数据库142中。该数据可以经由网络130被发送/接收,根据需要,通过通信接口114由处理器112进行访问。
另外,网络130可以被连接到互联网132以与位于远程并连接到互联网的服务器或客户端进行通信。如所描述的,网络130可以包括其它系统S1(134)、S2(136)、S3(138)。系统S1、S2和/或S3可以与系统100相同,或者可以是不同的系统。在一些实施例中,连接到网络130的一个或更多系统可形成协作执行图像采集、靶定位、和靶追踪以及向患者提供放射治疗的其他方面的分布式计算/仿真环境。
另外,放射治疗系统100可以与数据库140或医院数据库142进行通信,以便执行一个或更多远程存储的程序。举例来说,数据库140、医院数据库142、或二者可以包括诸如OracleTM数据库、SybaseTM数据库、或其他的关系数据库,并且可以包括诸如Hadoop序列文件、HBase、Cassandra或其他的非关系数据库。例如,这样的远程程序可以包括肿瘤信息系统(OIS)软件或治疗计划软件。例如,该OIS软件可以被存储在医院数据库142、数据库140、或OIS 150中。例如,治疗计划软件可以被存储在数据库140、医院数据库142、治疗计划系统160或OIS 150上。因此,例如,放射治疗装置110可以与医院数据库142进行通信以实现肿瘤信息系统150的功能。
然而,公开的实施例的系统和方法并不限于单独的数据库。在一个方面中,放射治疗系统100可以包括数据库220或医院数据库230。可替代地,数据库220和/或医院数据库230可以位于距放射治疗系统100远程处。数据库140和医院数据库142可以包括计算部件(例如,数据库管理系统,数据库服务器等),其配置为接收和处理对存储在数据库140或医院数据库142的存储装置中的数据的请求,并从数据库220或(一个或数个)医院数据库230提供数据。本领域技术人员将理解数据库140、142可以包括多个以集中式或分布式的方式定位的装置。
此外,放射治疗装置110可以通过网络130与数据库140进行通信以发送/接收存储在数据库140中多个各种类型的数据。例如,在一些实施例中,数据库140可以被配置为存储来自图像采集装置140的多个图像(例如,3D MRI、4D MRI、2D MRI切片图像、CT图像、2D荧光透视图像、X射线图像、来自MR扫描或CT扫描的原始数据、数字成像以及医学通信(DIMCOM)数据等)。数据库140可以存储待由靶定位模块128、训练模块124和治疗计划软件120使用的数据。如下所述,放射治疗装置110可以从数据库140接收成像数据(例如,3D MRI图像、4DMRI图像)以生成多个自适应滤波器模型。
进一步地,放射治疗系统100可以包括图像采集装置170,其可以采集患者的医学图像(例如,核磁共振成像(MRI)图像、3D MRI、2D流化MRI、4D容积MRI、计算机断层摄影(CT)图像、锥形束CT、正电子发射断层摄影(PET)图像、功能MRI图像(例如,fMRI、DCE-MRI和扩散MRI)、X射线图像、荧光图像、超声图像、放射治疗射野图像、单光子发射计算机断层摄影(SPECT)图像等)。例如,图像采集装置170可以是MRI成像装置、CT成像装置、PET成像装置、超声装置、荧光镜装置、SPECT成像装置、或用于获得患者的一个或更多医学图像的其他医学成像装置。通过图像采集装置170采集的图像可以作为成像数据和/或测试数据而被存储的数据库140内。作为示例,通过图像采集装置170采集的图像也可以由放射治疗装置110存储在存储器116中。
在一个实施例中,例如,图像采集装置170可以与放射治疗装置110集成作为单一的装置(例如,与直线加速器组合的MRI装置,也被称为“MRI-Linac”,或与伽玛刀相组合的集成的MRI装置)。例如,可以使用这样的MRI-Linac来确定患者中的靶器官或靶肿瘤的位置,从而根据放射疗法治疗计划将放射治疗导引到预定靶。
图像采集装置170可以被配置为获取对于关注区域(例如,靶器官,靶肿瘤或两者)的患者的解剖结构的一个或更多图像。每个2D切片可包括一个或更多参数(例如,2D切片厚度、取向和位置等等)。通过使用处理器112可以对该一个或更多参数进行调整以包含靶。例如,通过调节梯度或RF波形特性可以操纵2D切片的选择的特征。例如,通过改变RF脉冲的基本频率并保持相同的梯度强度可以改变切片的位置。另外,例如,通过使用身体上不同的梯度轴(例如,所选择的切片可以正交于应用的梯度)可以改变切片的取向。在一个实例中,图像采集装置170(例如,在MRI或MRI-Linac)可以获取任何取向的2D切片。例如,2D切片的取向可以包括矢状取向、冠状取向或轴状取向。这些取向可以对应于与MRI或MRI-Linac相关联的磁场梯度(例如,分别地Gx、Gy或Gz)。处理器112可以调整诸如2D切片的厚度的参数以包含靶器官或靶肿瘤。在一个实施例中,例如可以通过使用等式1来计算和调节切片的厚度:
Figure BDA0001399441100000091
在等式1中,T可以表示的2D切片厚度,如可以以距离为单位(例如,毫米)进行测量。Δωs可以代表以指定的频率(例如,赫兹,“Hz”)施加的射频脉冲所对应的激励带宽。字母γ可以代表称为旋磁比的常数(例如,对于质子,γ具有267.54MHz/特斯拉的值)。Gslice可以表示磁场梯度(例如,以特斯拉/米来测量)。在一个示例中,2D切片可从诸如3D MRI体积的信息来确定。当使用放射治疗装置110时,例如,在患者接受放射疗法治疗的同时,这种2D切片可以由图像采集装置170“实时”采集。
治疗计划系统160可以生成并存储用于待治疗的特定患者的放射疗法治疗计划、对于其他患者的放射疗法治疗计划、以及其他放射治疗信息(例如,射束角度、剂量直方图体积信息、治疗过程中使用的放射束的数量、射束角度、每束剂量等)。例如,治疗计划系统160可以提供关于要被施加到患者的特定的放射剂量的信息和其他放射治疗相关的信息(例如,诸如图像引导放射治疗(IGRT)、调强放射疗法(IMRT)、立体定向放射治疗等的治疗的类型)。
生成治疗计划可包括与所述图像采集装置170进行通信(例如,CT装置、MRI装置、PET装置、X射线装置、超声装置等),以便访问患者的图像并勾画诸如肿瘤的靶。在一些实施例中,可能需要对例如肿瘤周围或紧密接近肿瘤的健康组织的一个或更多危及器官进行勾画。因此,当OAR接近靶肿瘤时可以执行OAR的分割。另外,如果靶肿瘤靠近OAR(例如,靠近接近膀胱和直肠的前列腺),OAR的分割,所述治疗计划系统160可允许研究不仅靶中的剂量分布,还允许研究OAR中的剂量分布。
为了相对于OAR勾画出靶器官或靶肿瘤,可以通过图像采集装置170非侵入性地获得正接受放射治疗的患者的诸如MRI图像、CT图像、PET图像、fMRI图像、X射线图像、超声图像、放射治疗射野图像、SPECT图像等医学图像,以揭露身体部分的内部结构。基于来自医学图像中的信息,可以获取3D结构。此外,在治疗计划过程中,可以考虑许多参数以实现对靶肿瘤的有效治疗(例如,使得靶肿瘤接收对于有效疗法足够的放射剂量)与对(一个或数个)OAR的低照射(例如(一个或数个)OAR接收尽可能低的放射剂量)之间的平衡,该许多参数诸如靶器官和靶肿瘤的位置、OAR的位置、以及靶关于OAR的运动。例如,通过绘制MRI或CT图像的每个2D层或切片内的靶或OAR的轮廓并对每个2D层或切片的轮廓进行组合可以获取该3D结构。轮廓可以手动地生成(例如,由医师、剂量师、或医护人员)或自动地生成(例如,使用诸如由瑞典斯德哥尔摩的Elekta AB公司制造的基于图谱的自动分割软件,
Figure BDA0001399441100000101
的程序)。在某些实施方案中,可以由治疗计划系统160自动地生成靶肿瘤或OAR的3D结构。
在靶肿瘤和(一个或数个)OAR已被定位并勾画之后,剂量师、医生或医护人员可以确定要被施加到靶肿瘤和邻近肿瘤的任何OAR(例如,左和右腮腺、视神经、眼睛、晶状体、内耳、脊髓、脑干等)的放射剂量。在为每个解剖结构(例如,靶肿瘤、OAR)确定放射剂量后,可以执行称为逆向计划的过程以确定一个或更多治疗计划的参数,诸如体积勾画(例如,限定靶区、绘制敏感结构轮廓)、靶肿瘤和OAR周围的边界、剂量约束(例如,对肿瘤靶全剂量和对任何OAR零剂量;对PTV 95%的剂量,而脊髓≤45Gy、脑干≤55Gy、以及视神经结构<54Gy分别接收所述剂量)、射束角选择、准直器设置、以及射束开启时间。逆向计划的结果可以构成可以被存储在治疗计划系统160或数据库140中的放射治疗计划。这些处理参数中的一些可能相关的。例如,在改变治疗计划的尝试中调谐一个参数(例如,用于不同的目的的权重,如对靶肿瘤增大剂量)可能影响至少一个其它参数,这进而可能导致开发不同的治疗计划。因此,治疗计划系统160可以生成具有这些参数的定制的放射疗法治疗计划,以便使放射治疗装置110向患者提供放射治疗。
此外,该放射治疗系统100包括显示装置180和用户接口190。显示装置180可以包括将医学图像、接口信息、治疗计划参数(例如,轮廓、剂量、射束角度等)、治疗计划、靶、定位靶和/或追踪靶、或者任何相关的信息显示给用户的一个或更多显示屏幕。用户接口190可以是键盘、小型键盘、触摸屏或可以用户将信息输入到放射治疗系统100的任何类型的装置。
为了使放射治疗装置110对关注的解剖区域中的靶肿瘤进行定位和追踪,例如,处理器112可以执行靶定位模块128。如本公开中稍后描述的,靶获取模块120可利用由训练模块124生成自适应滤波器模型126。此外,自适应滤波器模型可以包括特定的域信息(例如,空间域或频域),位置信息(例如,体积中的3D位置、3D偏移信息、或2D位移信息)。如下所述,靶定位模块128使用由训练模块124生成的至少一个自适应滤波器模型126。
此外,在一个实施例中,放射治疗系统100的任何和所有部件可以被实现为虚拟机(例如,VM Ware,Hyper-V等)。例如,虚拟机可以是用作硬件的软件。因此,虚拟机可以包括共同起到硬件作用的至少一个或更多虚拟处理器、一个或更多虚拟存储器、以及一个或更多虚拟通信接口。例如,OIS 150、TPS 160、图像采集装置170可以被实现为虚拟机。考虑到可用的处理能力、存储器以及计算能力时,整个放射治疗系统可以被实现为虚拟机。
图2示出一种类型的放射疗法治疗装置200的示例,例如,由瑞典斯德哥尔摩的Elekta AB公司制造的Leksell伽玛刀。伽玛刀可以被配置为利用靶定位模块128(在图l中示出),其可以远程访问MRI图像(例如,从图像采集装置170)以在脑中定位靶肿瘤。在实施例中,作为图像采集装置170的MRI装置可以与伽玛刀集成。如图2所示,放射疗法分次治疗期间,患者210可以佩戴坐标架220,以使进行手术或放射疗法的患者的身体部分(例如,头部)保持稳定。坐标架220和患者定位系统230可以建立一个空间坐标系,其可以在对患者进行成像时或在放射外科手术期间使用。放射治疗装置200可以包括一个保护壳体240以包围多个放射源250。放射源250可以通过射束通道260生成多个放射束(例如,子束)。多个放射束可以被配置为从不同的方向聚集于等中心270。虽然每个单独的放射束可以具有相对低的强度,但当来自不同放射束的多个剂量积聚在等中心270时,等中心270可以接收相对高水平的放射。在某些实施方式中,等中心270可以对应于手术或治疗中的靶,例如肿瘤。
图3示出一种类型的放射治疗装置300的另一个示例(例如,直线加速器,被称为LINAC,由瑞典斯德哥尔摩的Elekta AB公司生产)。使用直线加速器300,患者302可以被定位在患者台304上以接收由放射由治疗计划系统160(图1中所示)生成的放射疗法治疗计划所确定的放射剂量。放射治疗计划可被用于定位和追踪体积中的3D靶,例如位于患者302的解剖结构内的靶器官或靶肿瘤。
直线加速器300可以包括连接到围绕患者302旋转的台架308的放射头306。放射头306生成指向靶器官或靶肿瘤的放射束310。当台架308旋转时,放射头306可以围绕患者302旋转。在旋转的同时,放射头306可以根据由治疗计划系统160(图1中所示)生成的治疗计划基于肿瘤的角度、形状和尺寸向患者302提供多种不同剂量的放射。由于器官和肿瘤在患者的身体内移动,靶的快速和准确的3D定位是很重要的。例如,靶器官或肿瘤可能由于各种类型的运动(例如,呼吸,心脏,蠕动或其它类型的患者运动)而移动。因此,直线加速器300可以被配置为通过使用靶定位模块128在放射疗法治疗期间定位所述靶(例如,器官或肿瘤)并追踪所述靶。
此外,患者台304的下面,可以设置平板闪烁体检测器312,其可以围绕定位在患者32体内的靶器官或靶肿瘤上的等中心314与放射头306同步旋转。平板闪烁器可以获取具有最高可达到信噪比的图像,并且可以用于验证在任何特定的放射疗法治疗分次(例如,放射疗法治疗可能需要多个分次的放射治疗,其中每个分次通常被称为“分次”)期间由患者302接收的放射量。此外,这样的图像被用于确定相对于放射头306的患者布置的几何精度。
轴线316与放射头306生成的射束310的中心的交点通常被称为“等中心”。患者台304可以被机动化,使得患者302可以被定位为使得肿瘤部位位于或接近等中心314。例如,患者台304可以相对于直线加速器300的一个或更多其他部件而改变位置,诸如相对于位于该放射头306的治疗放射源升高、改变患者302的纵向位置或横向位置。
在一个实施例中,直线加速器300可以与诸如磁共振成像装置的图像采集装置170(图1所示)集成为单个装置(例如,MRI-Linac)。在这样的情况下,当放射疗法治疗期间由于患者台304通过其插入的孔大小的有限尺寸导致患者台304被配置为不运动时,MRI-LINAC可以包括“虚拟检查台”,其可以相对于所述放射源“虚拟地”调整患者302的对准。在这样的实施例中,MRI-LINAC可以用于使用靶定位模块128确定靶的定位并追踪患者302中的靶,以便将放射治疗指向到预定的靶。
图4示出了用于训练模块124生成自适应滤波器模型126的过程400。靶定位模块128可以利用自适应滤波器模型126(图1所示)来确定靶的定位,然后对靶进行追踪。
在402,训练模块124(示于图1)从图像采集装置170接收特定患者的关注区域的多个图像(例如,3D MRI,4D MRI等)。例如,该多个图像可以是包含靶(例如,靶器官、靶肿瘤等)的关注区域的3D MRI图像或4D MRI图像。
在404,训练模块110利用处理器112继续从所述多个MRI图像提取包含关注区域中的靶(例如,靶器官、靶肿瘤等)的多个切片(例如,2D切片)。例如,可以在开始放射疗法治疗之前,预先确定2D切片的厚度(例如,由基于所述3D MRI体积或4D MRI体积的信息来确定)。该信息可以包括2D切片是应该在频域还是空间域。可替换地,训练模块110可以当正在用放射疗法对患者402进行治疗时(例如实时地)接收2D切片信息。例如,MRI-Linac可以用于利用放射疗法治疗患者402。放射疗法的这个过程期间,MRI-Linac可以在处理过程中获得多个3D MRI图像或4D MRI图像。提取的切片代表包围靶的患者的解剖结构的2D切片。此外,所提取的2D切片可以平行或正交于靶的运动。
在406,处理器112为每个提取的2D切片确定一个或更多2D图像块。在一个示例中,2D图像块可对应于可以包含靶的2D切片图像的区域。2D图像块可以被构造为任何形状(例如,正方形、矩形、圆形、多边形形状等),并且可以改变尺寸(例如,32×32像素的图像块、64×64像素的图像块等)。为了说明,例如,对应于512×512像素的2D切片中的靶的正方形图像块可以是是32×32像素。在示例中,2D图像块可以包括关于靶的信息以及对应于不包括靶的区域的信息(例如,背景信息)。
在408,处理器112确定的偏移值。例如,该偏移值是从2D图像块的中心到待追踪的靶的中心的距离。2D图像块可以包括域信息(例如,关于2D图像块的空间域或频域的信息)。此外,偏移信息可以包括例如从基准点(例如,基准点对应于体积中的靶)起2D图像块的位移、以及关于靶的形状变化(例如,变形)的信息。在示例中,基准点可以包含靶的计算的中心(例如,肿瘤的质心)。对每个包含靶的2D图像块确定偏移值。因此,每个2D图像块具有其自己的关联的偏移值。多个偏移值及其相关联的图像块被存储在存储器116中。
例如假设2D图像块的中心的坐标为(例如,[xp,yp,zp]),靶肿瘤的中心的坐标为(例如,[xt,yt,zt]),则与靶肿瘤的3D位置相关的2D图像块中心的偏移可以是向量(例如,[ox,oy,oz]=[xt-xp,yt-yp,zt-zp])。因此,处理器112可以追踪靶肿瘤。在追踪期间,处理器112可定位2D图像块的中心用于新患者摆位(例如,“在线/实时”诸如为[xpnew,ypnew,zpnew])。然后,处理器112可确定肿瘤的位置(例如,[xtnew,ytnew,ztnew]=[xpnew+ox,ypnew+oy,zpnew+oz])。
在410,处理器112可以使用一个或更多类型的自适应滤波器的设计技术将每个图像块转换到自适应滤波器模型126中。例如,可使用以下类型的自适应滤波器设计:匹配滤波器、最大间隔相关滤波器、综合判别函数滤波器、最小均方滤波器等。此外,自适应滤波器模型126可以包括域信息和从对应于该体积中的靶的2D图像块确定的偏移信息。
一般地,例如,自适应滤波器模型126“模拟”待定位的对象(例如,肿瘤)。可以生成多个自适应滤波器模型126。例如,可以有用于肿瘤的顶部的自适应滤波器模型126、用于肿瘤的中部的自适应滤波器模型126以及用于肿瘤的底部的自适应滤波器模型126。此外,可以为人体解剖结构的不同部分(例如,一个或更多关注的靶器官,诸如前列腺、乳腺、肺、心脏,脑等)创建自适应滤波器模型126。此外,可以为每个关注的靶(例如,在特定器官中的一种或更多肿瘤--例如在肺中多肿瘤;或在癌症已经转移且一个或更多器官中存在一种或更多肿瘤的情况等)生成自适应滤波器模型126。
此外,自适应滤波器模型126既可以设计在空间域中也可以设计在频域中的。在一个示例中,在频域中将自适应滤波器模型126应用到2D切片中可以是计算效率更高的。在患者的放射疗法治疗期间,自适应滤波器模型126可以存储在存储器116中(在图1中示出)并由处理器112取得以被应用到随后采集的2D切片(例如,“实时”/“在线”地采集的2D切片)来预测靶的定位。
图5示出了在患者的“实时”/“在线”放射疗法治疗期间使用一个或更多自适应滤波器模型126对靶进行追踪的工作流程过程500的示例性流程图。
在502,该过程500开始,放射治疗装置10(在图1中示出)使用处理器116访问靶定位模块128,所述靶定位模块128取得为包含特定患者的靶肿瘤的关注区域(例如,前列腺、乳腺、肺等)所设计的多个自适应滤波器模型126。
在504,处理器116接收采集协议以提供给图像采集装置170,以生成一组初始的2D切片。例如,该协议可以包括要取得的切片的定位(例如,肺、前列腺、脑、肾等)、切片的取向(例如,基于预定的可能的器官运动–诸如平行于或正交于要取得的切片)、以及切片厚度(例如,5mm切片、10mm的切片、15mm切片等)。用户可以提供对靶的定位的初始估计,该初始估计可以通过用户接口190来提供。处理器116可以利用该初始估计作为何处开始获取2D切片的初始位置。可替换地,可以由处理器116自动地确定这样的位置信息。
在506,放射治疗装置根据协议进一步接收来自图像采集装置170(例如,MRI装置或MRI-Linac装置)的多个2D MRI切片。该2D图像可以对应于具有肿瘤的关注区域。该2D切片可以平行于或正交于肿瘤。2D切片可以围绕肿瘤。另外,2D切片可以是肿瘤周围的区域。典型地,MRI装置提供频域中的2D切片。
如下所述,在508至514,处理器116确定肿瘤在2D切片上的位置及其在3D体积中的位置。
在508,处理器116可以将2D切片转换到频域或空间域。该2D切片被转换以匹配自适应滤波器模型126的域。例如,如果自适应滤波器模型12在空间域中被创建,则2D切片,例如,处理器116可以将2D切片转换到空间域。可选地,如果自适应滤波器模型126在频域中被创建,处理器116可以将2D切片转换到频域。
在510,处理器116可以应用自适应滤波器模型126。如上文所讨论,自适应滤波器模型126是由多个2D图像块生成的已经由训练模块126训练好的多个模型。处理器112将多个自适应滤波器模型126应用于2D切片。例如,2D切片和自适应滤波器模型可以具有相同的取向(例如,正交于运动的方向、平行于运动方向或两者)。在一个示例中,在频域中将自适应滤波器模型126应用到2D切片中可以是计算效率更高的。如下是自适应滤波器模型126可以如何应用于2D切片的一个示例:
在空间域中,自适应滤波器模型126可以由2D模板T(x,y)表示,并且2D切片(例如,2D图像)可由I(x,y)表示。将自适应滤波器模型126应用到2D切片的特定位置可提供相关性相关的“置信度得分”。例如,该置信度得分提供特定的自适应滤波器模型与2D切片的特定位置匹配得如何的指示。
在一个实例中,自适应滤波器模型126相对于2D切片的特定位置具有越好的匹配,则置信度得分越高。置信度得分可以被用来预测靶肿瘤是否位于2D切片内以及位于2D切片内何处。许多不同类型的相关可以被用于根据情况确定置信度得分。例如,可以使用以下类型的相关:互相关、归一化互相关、或相关系数,其可以被定义如下:
互相关:R(x,y)=∑i,j∈PT(i,j)·I(x-i,y-j);
归一化互相关:
Figure BDA0001399441100000161
相关系数:
Figure BDA0001399441100000162
(例如,相关系数等于平均校正的模板和图像的归一化互相关)。
在上述等式中:P表示模板的空间域;
Figure BDA0001399441100000163
Figure BDA0001399441100000164
分别表示模板T和切片I的平均值。此外,N(x,y)和C(x,y)二者都是有界的,例如,1≤N(x,y),C(x,y)≤+1,使得它们的值可以被解释为相对于完美匹配结果+1的置信度得分。例如,该置信度得分可以包括诸如上述定义的R,N或C的信息。
也可以通过应用卷积定理在频域中进行这三个相关计算,例如:
Figure BDA0001399441100000165
Figure BDA0001399441100000166
其中
Figure BDA0001399441100000167
Figure BDA0001399441100000168
分别表示正向傅里叶变换和逆向傅里叶变换。
在512,处理器112可确定响应图,作为将自适应滤波器模型126应用到2D切片的结果。对于单一2D切片可以创建多个响应图,每个响应图可以对应于特定自适应滤波器模型。自适应滤波器模型126和2D切片二者可以具有相同的取向。该响应图可以在空间域或频域中生成。2D切片和自适应滤波器模型126二者必须处于空间域中以生成在空间域中的响应图。类似地,2D切片和自适应滤波器模型126二者必须都处于频域中以生成频域中的响应图。可以为每个2D切片来创建响应图。所生成的响应图可以是指示特定的自适应滤波器模型126与给定的2D切片内不同位置的匹配程度的图像。
例如,自适应滤波器模型A(未示出)可以由处理器112应用到2D切片的每个位置上。对于每个位置,例如,则处理器112可以确定指示自适应滤波器A与2D切片的特定位置匹配得如何的置信度得分。自适应滤波器模型A相对于2D切片上的位置之间匹配得越好,则该特定位置的置信度得分越高(例如,越接近+1)。响应图图像上更亮的区域可以指示更好的匹配,并且因此指示比暗区更高的置信度得分(例如在自适应滤波器模型A与2D切片的特定位置之间)。在响应图中,较亮的区域可以具有比较暗的区域更高的置信度得分。之后,自适应滤波器模型A被应用到2D切片,可以生成响应图A。处理器112可以取得另一个自适应滤波器模型B(未示出),例如,并重复将自适应滤波器模型B应用到2D切片上的每一个位置以生成响应图B的过程。因此,处理器使用靶定位模块可以将多个自适应滤波器模型应用于2D切片以生成多个响应图,其可以被存储在存储器116中。在2D切片的每个位置处,生成对应于所使用的自适应滤波器模型的置信度得分。高置信度得分可以对应于2D切片上靶肿瘤可能处于的位置。
在514,处理器112可以基于置信度得分预测2D切片上靶肿瘤的位置。例如,对于响应图上的特定位置,可以有多个置信度得分值,其中每个置信度得分可以对应于某一特定自适应滤波器模型。为了预测肿瘤的位置,则处理器112可以选择对于2D切片上特定位置的最大置信度得分值。可替换地,处理器112可以使用对于2D切片上特定位置的所有置信度得分的加权平均值。在一个示例中,处理器112可以计算最大的N个置信度得分的加权平均值,其中N是一个预定数量。利用通过使用多个自适应滤波器模型生成的置信度得分可以提高预测的精确度。因此,对于特定位置的高置信度得分可以指示2D切片上靶的位置,而低置信度得分可以指示肿瘤不在该位置。置信度得分最接近1可以指示一部分肿瘤或整个肿瘤位于该2D切片上。
如果处理器112确定该2D切片不包含肿瘤,过程500追随501返回到506以获得另一2D切片。过程500可以使处理器112调整采集协议以采集下一2D切片。可替代地,根据切片厚度、肿瘤的尺寸、肿瘤的取向和其他因素,肿瘤可能被包含在一个或更多2D切片内。例如,如果肿瘤的底部被包含在当前2D切片中,处理器112可以使用偏移信息来确定随后应采集哪些2D切片(例如,包含肿瘤的质心的2D切片,或者可以是包含肿瘤的顶部的2D切片)。
例如,使用来自特定自适应滤波器模型及其相应的响应图的信息,处理器112可以确定先前2D切片上肿瘤的位置。另外,使用与自适应滤波器模型相关联的偏移信息,处理器112还可以估计3D体积中肿瘤的位置(例如,中心)。例如,处理器112可确定该肿瘤在例如第(5,10)像素位置与2D切片相交,则肿瘤的中心可能沿切片法线方向距2D切片的中心5mm。为了估计肿瘤的下一个位置(例如,随着肿瘤运动而对其进行追踪)或者为了获取关于肿瘤的形象的更完整的信息(例如,肿瘤的形状是否已经改变),处理器112可以调整各种参数(例如,随后的2D切片将被采集的位置、采集平行于先前的切片但处于不同位置的随后的2D切片、采集正交于先前的切片的2D切片、调整随后的切片的厚度等)。另外,使用来自先前的自适应滤波器模型的偏移信息,处理器112还可以估计肿瘤的下一位置。因此,当患者经历放射治疗时,处理器112可以实时地追踪肿瘤的运动。
在516,当患者“实时”经历放射疗法治疗时,处理器112可以对肿瘤进行追踪。此外,当多个2D切片被利用时,处理器112可以估计肿瘤的形状。
图6图示地示出应用于2D MRI切片以确定靶(例如,肝脏)的位置的示例性自适应滤波器模型。在602处,形象地表示在空间域中的自适应滤波器模型的一个示例。处理器112(示于图1)可以使用靶定位模块129以将602中所示的自适应滤波器模型应用于也处于空间域的特定的2D切片。通过将示例性自适应滤波器模型602应用到2D切片,处理器112可以为2D切片的每个位置确定,示例性自适应滤波器模型602与2D切片的特定位置匹配得如何。有利地,处理器112可以在约300毫秒内计算每个2D切片,靶定位模块129生成示例性响应图,其图示地示出在604中。示例性响应图604可以表示自适应滤波器模型602与2D切片的每个位置匹配得如何。在示例性响应图604示出的较亮的区域表示在2D切片的特定位置处与自适应滤波器模型602之间相较较暗的区域更好的匹配,并且因此,更高的置信度得分。在响应图604中,较亮的区域可以具有比较暗的区域更高的置信度得分。如图所示,响应图604指示如处理器112所确定的可以对应于该2D切片的最高置信度得分的亮区域608。因此,响应图604上的亮区域608可以对应于这一特定2D切片上的靶的位置。因此,处理器112已经确定位于该特定的2D切片上的靶608的位置。此外,由于通过图像采集装置170从其生成的2D切片的患者的解剖结构的位置,则处理器112可以确定606中示出的患者的解剖结构中靶608的位置。例如,606中的红色区域对应于覆盖患者的解剖结构的自适应滤波器模型602。例如,在另一个实施方案中,靶可以表示肿瘤。
可替代地,例如,根据2D切片和自适应滤波器模型是否处于频域中,可以在频域中表示响应图。为了在频域中生成响应图,处理器112可对2D切片进行快速傅立叶变换(FFT)。处理器112通过在频域中执行逐点相乘而应用自适应滤波器模型。为了生成响应图,处理器112可以进行快速傅立叶逆变换。
其他注意事项
本文件中提及的所有出版物、专利和专利文献通过引用整体并入本文,如同通过引用单独并入。如果本文献与通过引用并入的那些文献之间的用法不一致,则引用的参考文献中的用法应被视为对本文档的补充;对于不可调和的不一致,则以本文件中的用法为准。
示例性公开的实施例描述了当患者经历放射疗法治疗时用于对患者的解剖结构中的靶(例如,肿瘤)进行追踪的系统和方法。前面的描述已经出于说明的目的而示出。这并非穷尽的且不限于所公开的准确的形式或实施例。考虑到说明书及所公开的实施例的实践,变形例和实施例的修改将是显而易见的。
在本文献中,使用了术语“a”或“an”,如专利文件中常用的,其包括一个或多于一个,独立于“至少一个”或“一个或更多”的任何其它实例或用途。在本文献中,术语“或”用于表示非排他性的“或”,除非另有说明,否则“A或B”包括“有A没有B”、“有B没有A”、以及“A和B”。在本文献中,术语“包括”和“其中(in which)”用作各个术语“包括”和“其中(wherein)”的通俗易懂的等效。此外,在下列权利要求中,术语“包括”是开放式的,即,包括除了在权利要求中的这样的术语之后列出的那些元件之外的元件的设备、系统、装置、物品、构成、构想或过程仍被视为落入该权利要求的范围。此外,在下面的权利要求中,术语“第一”、“第二”和“第三”等仅用作标记,并不意图对其对象施加数值要求。
本文描述的方法示例至少部分地可以是机器或计算机实现的。一些示例可以包括编码有指令的计算机可读介质或机器可读介质,该指令可经操作以配置电子装置来执行如上述示例中所描述的方法。这些方法的实现可以包括例如微代码、汇编语言代码、更高级语言代码等的软件代码。可以使用各种软件编程技术来创建各种程序或程序模块。例如,可以用或通过Java、Python、C、C++、汇编语言、或任何已知的编程语言来设计程序段或程序模块。一个或更多这样的软件段或模块可以被集成到计算机系统和/或计算机可读介质中。这样的软件代码可以包括用于执行各种方法的计算机可读指令。这种软件代码可以形成计算机程序产品或计算机程序模块的部分。此外,在一个示例中,例如在运行期间或在其它时间,代码可以有形地存储在一个或更多易失性、非暂时性或非易失性有形计算机可读介质上。这些有形的计算机可读介质的示例可以包括但不限于硬盘、可移动磁盘、可移动光盘(例如,压缩光盘和数字视频盘)、磁带、存储卡或存储棒、随机存取存储器(RAM)、只读存储器(ROM)等。
此外,尽管已经在本文中描述了示例性实施例,其范围包括任何和所有基于本公开的具有等同元件、修改、省略、组合(例如,各种实施例的方案的)、改编或改变的实施例。权利要求书中的元件将被基于权利要求中采用语言广泛地解释,并不限于在本说明书中或本申请的进行期间所描述的示例,其示例将被解释为非排他性的。此外,可以以任何方式对所公开的方法的步骤进行修改,包括通过对步骤重新排序或插入或删除步骤。因此,本说明书和示例旨在仅被认为是示例,真正的范围和精神由以下权利要求以及其等同物的全部范围所指示。
以上描述旨在是说明性的而不是限制性的。例如,上述示例(或其一个或更多方案)可以彼此组合使用。例如本领域普通技术人员在阅读上述描述时可以使用其它实施例。提供摘要以符合37 C.F.R.§1.72(b),从而允许读者快速确定技术公开的实质。提交时应该理解是:它不会用于解释或限制权利要求的范围或含义。另外,在上述具体实施方式中,各种特征可以被分组在一起以简化本公开。这不应解释为一种使不要求保护的公开的特征对于任一权利要求都是必要的意图。而是,本发明的主题可以少于公开的特定实施例的全部特征。因此,以下权利要求书作为示例或实施例并入具体实施方式中,其中每个权利要求自身独立地作为单独的实施例,并且认为这些实施例可以以各种组合或排列彼此组合。本发明的范围应参照所附权利要求以及这些权利要求赋权的等同形式的全部范围来确定。

Claims (26)

1.一种由处理器实现的用于在三维(3D)体积中定位靶的方法,所述方法包括:
从图像采集装置接收包含所述靶的关注区域的3D医学图像;
由所述处理器从所述3D医学图像中提取多个二维(2D)切片;
由所述处理器为每个2D切片确定2D图像块,其中所述2D图像块对应于与所述靶相关联的所述2D切片的区域;
由所述处理器确定反映所述2D图像块的位置和所述靶的位置之间的距离的偏移值;以及
由所述处理器将所述2D图像块转换到自适应滤波器模型,用于确定所述关注区域中所述靶的位置,其中,所述自适应滤波器模型包括所述偏移值。
2.根据权利要求1所述的方法,其中转换所述2D图像块包括使用所述自适应滤波器模型,并且其中所述自适应滤波器模型包括匹配滤波器、最大间隔相关滤波器、综合判别函数滤波器和最小均方滤波器中的至少一种。
3.根据权利要求1所述的方法,其中所述偏移值包括在所述2D图像块上的点和所述3D体积内的所述靶内的基准点之间的距离向量。
4.根据权利要求1所述的方法,其中所述2D切片大致平行于所述3D体积内的靶运动。
5.根据权利要求1所述的方法,其中所述2D切片大致正交于所述3D体积内的靶运动。
6.根据权利要求1所述的方法,其中所述自适应滤波器模型和所述2D切片二者具有相同的取向。
7.根据权利要求1所述的方法,其中所述2D切片处于频域中或空间域中。
8.根据权利要求1所述的方法,其中所述自适应滤波器模型处于频域中或处于空间域中。
9.根据权利要求1所述的方法,其中所述自适应滤波器模型和所述2D切片二者处于相同的域中,其中,所述域是频域或空间域。
10.根据权利要求1所述的方法,进一步包括:由所述处理器将所述2D切片转换到所述自适应滤波器模型的域,其中所述域是频域或空间域中的至少一个;并且其中所述自适应滤波器模型处于与所述2D切片相似的取向。
11.根据权利要求1所述的方法,进一步包括:在从3D医学图像中提取所述多个二维(2D)切片之前,由所述处理器调节参数。
12.根据权利要求11所述的方法,其中所述参数包含靶位置参数。
13.根据权利要求11所述的方法,其中所述参数包括初始靶位置参数。
14.根据权利要求11所述的方法,其中所述参数包括切片取向参数。
15.权利要求14所述的方法,其中所述切片取向参数包括矢状取向、冠状取向或轴状取向中的至少一种。
16.根据权利要求11所述的方法,其中所述参数包括切片厚度参数。
17.根据权利要求11所述的方法,其中所述参数包括3D体积中靶运动的方向参数。
18.根据权利要求1所述的方法,其中所述自适应滤波器模型包括置信度得分,其中所述置信度得分指示所述自适应滤波器模型与所述2D切片的特定位置之间的相对匹配程度。
19.根据权利要求1所述的方法,其中所述处理器生成多个自适应滤波器模型。
20.根据权利要求19所述的方法,其中所述多个自适应滤波器模型中的每一个被针对包括所述靶的特定的关注区域而设计。
21.根据权利要求19所述的方法,其中所述处理器将两个或更多自适应滤波器模型进行组合,并且每个自适应滤波器模型对应于不同的切片取向。
22.根据权利要求1所述的方法,其中所述2D图像块对应于包括所述靶的所述2D切片的区域。
23.根据权利要求1所述的方法,其中所述2D图像块被配置为包括对应于围绕所述靶但不包括所述靶的区域的信息。
24.根据权利要求1所述的方法,其中所述2D图像块处于空间域中或处于频域中。
25.根据权利要求1所述的方法,其中每个2D图像块包括反映所述2D图像块的位置与肿瘤的位置之间的距离的相应的偏移值。
26.一种使用体积的二维成像切片追踪所述体积中的三维(3D)靶的系统,所述系统包括:
处理器;以及
存储器,其可操作地联接到所述处理器并且存储计算机可执行指令,所述计算机可执行指令在由所述处理器执行时使所述处理器执行方法,所述方法包括:
从图像采集装置接收包含所述靶的关注区域的所述体积的3D医学图像,其中所述3D医学图像被存储在所述存储器中;
由所述处理器从所述3D医学图像中提取多个二维(2D)切片;
由所述处理器为每个2D切片确定2D图像块,其中所述2D图像块对应于与所述靶相关联的所述2D切片的区域;
由所述处理器确定反映所述2D图像块的位置和所述靶的位置之间的距离的偏移值;以及
由所述处理器将所述2D图像块转换到自适应滤波器模型,用于确定所述关注区域中所述靶的位置,其中,所述自适应滤波器模型包括所述偏移值。
CN201680013998.6A 2015-01-28 2016-01-21 用于自适应放射治疗的运动靶的三维定位 Active CN107530552B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US14/607,626 2015-01-28
US14/607,626 US9652871B2 (en) 2015-01-28 2015-01-28 Three dimensional localization of a moving target for adaptive radiation therapy
PCT/US2016/014302 WO2016122961A1 (en) 2015-01-28 2016-01-21 Three dimensional localization of a moving target for adaptive radiation therapy

Publications (2)

Publication Number Publication Date
CN107530552A CN107530552A (zh) 2018-01-02
CN107530552B true CN107530552B (zh) 2020-06-19

Family

ID=55485282

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201680013998.6A Active CN107530552B (zh) 2015-01-28 2016-01-21 用于自适应放射治疗的运动靶的三维定位

Country Status (7)

Country Link
US (2) US9652871B2 (zh)
EP (1) EP3251089A1 (zh)
JP (1) JP6624695B2 (zh)
CN (1) CN107530552B (zh)
AU (1) AU2016211820B2 (zh)
RU (1) RU2699052C2 (zh)
WO (1) WO2016122961A1 (zh)

Families Citing this family (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10092251B2 (en) * 2013-03-15 2018-10-09 Varian Medical Systems, Inc. Prospective evaluation of tumor visibility for IGRT using templates generated from planning CT and contours
US9652871B2 (en) 2015-01-28 2017-05-16 Impac Medical Systems, Inc. Three dimensional localization of a moving target for adaptive radiation therapy
US11182907B2 (en) * 2016-02-02 2021-11-23 Elekta Ltd. Motion management in image-guided radiotherapy using contour elements
EP3436148B1 (en) * 2016-03-30 2021-11-17 Koninklijke Philips N.V. Adaptive radiation therapy planning
US11517768B2 (en) * 2017-07-25 2022-12-06 Elekta, Inc. Systems and methods for determining radiation therapy machine parameter settings
EP3684463A4 (en) 2017-09-19 2021-06-23 Neuroenhancement Lab, LLC NEURO-ACTIVATION PROCESS AND APPARATUS
US11717686B2 (en) 2017-12-04 2023-08-08 Neuroenhancement Lab, LLC Method and apparatus for neuroenhancement to facilitate learning and performance
EP3731749A4 (en) 2017-12-31 2022-07-27 Neuroenhancement Lab, LLC NEURO-ACTIVATION SYSTEM AND METHOD FOR ENHANCING EMOTIONAL RESPONSE
EP3543723A1 (en) * 2018-03-20 2019-09-25 Siemens Healthcare GmbH Method for acquiring magnetic resonance image data for image-guided radiotherapy
US11364361B2 (en) 2018-04-20 2022-06-21 Neuroenhancement Lab, LLC System and method for inducing sleep by transplanting mental states
CN108875080B (zh) * 2018-07-12 2022-12-13 百度在线网络技术(北京)有限公司 一种图像搜索方法、装置、服务器及存储介质
CA3112564A1 (en) 2018-09-14 2020-03-19 Neuroenhancement Lab, LLC System and method of improving sleep
US10918885B2 (en) * 2018-09-27 2021-02-16 Varian Medical Systems International Ag Systems, methods and devices for automated target volume generation
US20210316156A1 (en) * 2018-12-26 2021-10-14 Our United Corporation Positioning method realized by computer, and radiotherapy system
CN110490851B (zh) * 2019-02-15 2021-05-11 腾讯科技(深圳)有限公司 基于人工智能的乳腺图像分割方法、装置及系统
US11141609B2 (en) 2019-05-15 2021-10-12 Elekta Ab (Publ) Dose guided real-time adaptive radiotherapy
CA3156832A1 (en) 2019-10-03 2021-04-08 Cornell University Optimizing bra sizing according to the 3d shape of breasts
US11077320B1 (en) 2020-02-07 2021-08-03 Elekta, Inc. Adversarial prediction of radiotherapy treatment plans
US20210268313A1 (en) * 2020-02-27 2021-09-02 Elekta ltd Image Registration in Treatment Planning
US11690581B2 (en) * 2020-09-07 2023-07-04 C-Rad Positioning Ab Tumor position determination
US11904186B2 (en) * 2022-01-12 2024-02-20 Capital Medical University Radiotherapeutical or radiosurgical system comprising two or more rotatable high-intensity radiation sources and a ring-shaped imager, and methods thereof
CN114332090B (zh) * 2022-03-16 2022-05-10 中南大学 多源域自适应脑网络分类方法、系统、设备及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101553281A (zh) * 2006-11-02 2009-10-07 艾可瑞公司 使用直接目标配准的目标跟踪
WO2014144469A1 (en) * 2013-03-15 2014-09-18 Varian Medical Systems, Inc. Prospective evaluation of tumor visibility for igrt using templates generated from planning ct and contours
CN104718563A (zh) * 2012-10-01 2015-06-17 三菱电机株式会社 追踪三维物体的方法

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002165894A (ja) * 2000-12-04 2002-06-11 Mitsubishi Electric Corp 放射線治療システム
US7577282B2 (en) * 2002-11-27 2009-08-18 Hologic, Inc. Image handling and display in X-ray mammography and tomosynthesis
US20050053267A1 (en) * 2003-09-05 2005-03-10 Varian Medical Systems Technologies, Inc. Systems and methods for tracking moving targets and monitoring object positions
JP4549093B2 (ja) * 2004-04-12 2010-09-22 キヤノン株式会社 画像処理装置及びその方法、プログラム
US8989349B2 (en) * 2004-09-30 2015-03-24 Accuray, Inc. Dynamic tracking of moving targets
US8042209B2 (en) * 2005-04-13 2011-10-25 University Of Maryland Techniques for compensating movement of a treatment target in a patient
WO2013116694A1 (en) * 2012-02-03 2013-08-08 The Trustees Of Dartmouth College Method and apparatus for determining tumor shift during surgery using a stereo-optical three-dimensional surface-mapping system
JP5028528B2 (ja) * 2008-09-30 2012-09-19 株式会社日立メディコ X線ct装置
US9433389B2 (en) * 2011-07-12 2016-09-06 University Of Maryland, Baltimore Method for monitoring the accuracy of tissue motion prediction from surrogates
US8837778B1 (en) * 2012-06-01 2014-09-16 Rawles Llc Pose tracking
RU2504329C1 (ru) * 2012-10-02 2014-01-20 Федеральное государственное бюджетное учреждение "Научный центр неврологии" Российской академии медицинских наук Способ выявления в коре головного мозга сенсомоторных зон, ответственных за локомоцию
JP6351017B2 (ja) * 2014-02-24 2018-07-04 国立研究開発法人量子科学技術研究開発機構 放射線治療用動体追跡装置、放射線治療用照射領域決定装置および放射線治療装置
US10043284B2 (en) * 2014-05-07 2018-08-07 Varian Medical Systems, Inc. Systems and methods for real-time tumor tracking
US20160213948A1 (en) * 2015-01-22 2016-07-28 Impac Medical Systems, Inc. System and method to create test objects for radiation therapy
US9652871B2 (en) 2015-01-28 2017-05-16 Impac Medical Systems, Inc. Three dimensional localization of a moving target for adaptive radiation therapy

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101553281A (zh) * 2006-11-02 2009-10-07 艾可瑞公司 使用直接目标配准的目标跟踪
CN104718563A (zh) * 2012-10-01 2015-06-17 三菱电机株式会社 追踪三维物体的方法
WO2014144469A1 (en) * 2013-03-15 2014-09-18 Varian Medical Systems, Inc. Prospective evaluation of tumor visibility for igrt using templates generated from planning ct and contours

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
4D tumor centroid tracking using orthogonal 2D dynamic MRI:implications for radiotherapy planning;Erik Tryggestad,Aaron Flammang;《medical physics》;20130820;第40卷(第9期);摘要,第2.C部分,第3-4页,附图1 *
Evaluation of template matching for tumor motion management with cine-MR image in lung cancer patients;Xiutao Shi,Tejan Diwanji;《medical physics 》;20140430;第41卷(第5期);摘要,第2.C部分,第3-4页,附图1 *
Laura I cerviño,Jiang Du,Steve B Jiang.MRI-guided tumor tracking in lung cancer radiotherapy.《physics in medicine and biology》.2011,第56卷(第13期), *
MRI-guided tumor tracking in lung cancer radiotherapy;Laura I cerviño,Jiang Du,Steve B Jiang;《physics in medicine and biology》;20110531;第56卷(第13期);摘要,第3776页最后一段,第2.5部分,第3778-3779页,第2.6部分,附图3 *
Xiutao Shi,Tejan Diwanji.Evaluation of template matching for tumor motion management with cine-MR image in lung cancer patients.《medical physics 》.2014,第41卷(第5期), *

Also Published As

Publication number Publication date
US10152790B2 (en) 2018-12-11
US20160217595A1 (en) 2016-07-28
JP2018506349A (ja) 2018-03-08
CN107530552A (zh) 2018-01-02
US20170221206A1 (en) 2017-08-03
RU2699052C2 (ru) 2019-09-03
EP3251089A1 (en) 2017-12-06
RU2017130025A3 (zh) 2019-03-22
WO2016122961A1 (en) 2016-08-04
US9652871B2 (en) 2017-05-16
RU2017130025A (ru) 2019-02-28
AU2016211820A1 (en) 2017-08-17
JP6624695B2 (ja) 2019-12-25
AU2016211820B2 (en) 2019-06-06

Similar Documents

Publication Publication Date Title
US10987522B2 (en) Three dimensional localization and tracking for adaptive radiation therapy
CN107530552B (zh) 用于自适应放射治疗的运动靶的三维定位
CN107708808B (zh) 用于在放射治疗期间使用质量指标进行靶追踪的系统和方法
JP6725773B2 (ja) 放射線治療中の形状変化を考慮するためのシステムおよび方法

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