CN113674834A - 基于剂量分布预览系统的放疗靶区确立和修正方法 - Google Patents

基于剂量分布预览系统的放疗靶区确立和修正方法 Download PDF

Info

Publication number
CN113674834A
CN113674834A CN202110938598.9A CN202110938598A CN113674834A CN 113674834 A CN113674834 A CN 113674834A CN 202110938598 A CN202110938598 A CN 202110938598A CN 113674834 A CN113674834 A CN 113674834A
Authority
CN
China
Prior art keywords
dose
data
dose distribution
ptv
module
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
Application number
CN202110938598.9A
Other languages
English (en)
Inventor
于金明
张功森
王琳琳
朱健
穆向魁
李兆斌
刘新超
贺小龙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai United Imaging Healthcare Co Ltd
Original Assignee
Beijing Huaguang Putai Biotechnology Co ltd
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 Beijing Huaguang Putai Biotechnology Co ltd filed Critical Beijing Huaguang Putai Biotechnology Co ltd
Priority to CN202110938598.9A priority Critical patent/CN113674834A/zh
Publication of CN113674834A publication Critical patent/CN113674834A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H20/00ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • G16H20/40ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to mechanical, radiation or invasive therapies, e.g. surgery, laser therapy, dialysis or acupuncture
    • 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/1031Treatment planning systems using a specific method of dose optimization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Public Health (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Epidemiology (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Computational Linguistics (AREA)
  • Pathology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Primary Health Care (AREA)
  • Artificial Intelligence (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Urology & Nephrology (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Databases & Information Systems (AREA)
  • Radiation-Therapy Devices (AREA)

Abstract

本发明实施例公开了一种基于剂量分布预览系统的放疗靶区确立和修正方法,包括:读取CT影像数据、PTV轮廓和OARs轮廓;对CT影像数据和PTV轮廓进行处理,得到处理结果;所述处理结果包括3D剂量分布图、DVH图、所关注的剂量指数和HI/CI指标;基于处理结果,对PTV轮廓进行“雕刻式”的修正,使得修正后的PTV轮廓在保证病灶受量和降低危机器官风险之间达到平衡。本发明基于深度学习设计剂量分布预览系统,通过该系统简化了靶区修改和计划制定流程,提高了放疗靶区确立的合理性和放疗计划质量,缩短了肿瘤患者的治疗延迟时间,降低了工作人员负担的同时提高患者放疗疗效。

Description

基于剂量分布预览系统的放疗靶区确立和修正方法
技术领域
本发明涉及计算机技术领域,具体涉及一种基于剂量分布预览系统的放疗靶区确立和修正方法。
背景技术
随着放射生物学和计算机等技术的发展,放射治疗(Radiation Therapy,RT)在肿瘤治疗中日益发挥不可替代的作用。放疗计划(RT Plan)的合理制定是确保放疗精确性和良好疗效的基础,而在此过程中,放疗靶区的确定是计划制定的重要前提。靶区通常分为三个层次:GTV(Gross Tumor Volume,大致肿瘤靶区)、CTV(Clinical Target Volume,临床靶区)和PTV(Planning Target Volume,计划靶区)。其中,PTV为射野布局、剂量计算和计划评价的关键参考对象,除了GTV所代表的临床诊断证实肿瘤部位以及CTV所包含的亚临床病变范围之外,肿瘤医生在确定PTV区域时也会结合体积变化、器官运动以及影像误差等因素来对CTV进行扩展。
实际上,PTV的确定是一个反复修正的过程,且此过程难以被机器替代。首先,肿瘤医师需要考虑患者的年龄、病情等个体化差异,例如考虑到儿童脊柱生长情况,需要使其两侧剂量尽量对称;第二,患者靶区周围的解剖学差异可能导致危机器官(Organ at risk,OARs)受量缺乏规律性;第三,医师水平和工作经验不一,导致PTV确定过程的合理性有所差异;此外,在同步推量(Simultaneous integrated boost,SIB)等调强放疗技术中,需要考虑不同限定区域的多层级三维剂量分布情况,使得靶区的确定和修改工作更为复杂化。
目前,放射治疗中靶区确定和修正流程主要依靠肿瘤医师的临床经验和技巧,且在原始病灶较大或位置特殊以及SIB(同步推量)等病例计划中,需要结合物理师的放疗计划优化情况,对靶区进行反复调整,现阶段放疗计划从指定到实施的具体流程如下:
(1)肿瘤医师基于CT或多模态(MRI、PET-CT等)融合的医学影像数据,结合检查或检验结果以及自身临床经验确定GTV位置和区域,并在此基础上外扩得到CTV区域。
(2)在所确立CTV区域的基础上,外扩得到初始PTV轮廓,在所确立CTV区域的基础上,外扩得到初始PTV轮廓,同时勾画OARs轮廓。医师根据自身临床经验以及肿瘤患者实际情况对初始PTV轮廓进行修改。
(3)在富有经验或高级医师的指导下,对病灶及周围组织的治疗或预防剂量和危机器官剂量分布情况进行综合预判,寻求保证理想治疗剂量和保护危机器官之间的平衡点,并基于此对PTV进行修改。
(4)PTV、OARs及相关影像数据上传至TPS,由物理师进行设备选择、射野布局、约束条件设置、剂量计算和计划优化及预评估等工作,计划设计完成后反馈给医师评估。
(5)肿瘤医师根据物理师反馈的三维剂量分布情况和PTV与OARs的剂量指数(或DVH图)以及HI\CI等指标对计划进行评估,并根据临床规定,要求物理师进一步优化计划,计划需反复调整直至最优状态。
(6)对于反复优化但仍无法满足临床规定(靶区受量足够高、均匀性和适形度较好,危机器官风险足够低)的放疗计划,由肿瘤医师合理修改PTV轮廓,并反馈给物理师重新制定放疗计划。
(7)物理师重新制定放疗计划,并不断重复(5)-(6)步骤,直至放疗计划满足临床要求,并且保证病灶(及周围必要区域)接受足够治疗剂量的情况下,充分保护危机器官。
(8)放疗计划经由医师、高级物理师评估和审核后交付临床验证和实施。
上述方案主要存在以下缺点:
(1)肿瘤医师在修改放疗靶区过程中,由于缺乏当前PTV所对应的整体剂量分布预判、剂量指数(或DVH图)、HI\CI等预估手段,导致PTV区域的确立难以在保证处方治疗剂量(对于病灶及其周围必要区域)和降低危机器官(正常组织器官)风险之间找到平衡。一方面,为了保证病灶受量,同时考虑消除器官运动、形状变化、摆位误差及影像误差等因素所带来的照射剂量不确定性,可能会使PTV边界外扩过度;另一方面,对于原始病灶体积较大,紧邻、重叠危机器官或者SIB(同步推量)的病例,由于考虑最大程度降低危机器官受量,可能会导致PTV边界外扩不足。
(2)不同的临床工作者由于临床经验和业务水平差异,在PTV区域修正的速度和效果都存在差异,相应地,其将对计划制定的周期和质量造成影响。
(3)对于反复优化但仍未获得理想剂量分布、剂量指数(或DVH图)以及HI\CI等指标仍不满足临床要求的放疗计划,物理师难以明确其原因来自于PTV的不合理或是来自于计划设计流程失误。
(4)肿瘤医师在每次在PTV进行修改之后,只能根据临床经验预估剂量分布、剂量指数(或DVH图)以及HI\CI指标的相应改变,其需要物理师通过TPS重新进行计划设计以对其验证,若仍达不到临床要求,则需反馈至医师再次对PTV进行修改,而后重复此过程。因此,对于某些特殊(病灶体积或位置)病例,PTV的修正和计划重新制定过程将耗费大量时间和精力,患者治疗时间的延迟及在此期间导致的病灶体积变化降低放疗精度,造成实际病灶剂量不足或正常组织器官受损等。
发明内容
本发明实施例的目的在于提供一种基于剂量分布预览系统的放疗靶区确立和修正方法、剂量分布预览系统,以简化靶区修改和计划制定流程,提高放疗靶区确立的合理性和放疗计划质量,缩短肿瘤患者的治疗延迟时间,降低工作人员负担的同时提高患者放疗疗效。
为实现上述目的,第一方面,本发明实施例提供了一种基于剂量分布预览系统的放疗靶区确立和修正方法。其中,所述剂量分布系统包括数据读取与处理模块、病例设定与剂量指数选择模块、3D剂量分布图预览模块、计算模块和通信模块,所述方法包括:
通过所述数据读取与处理模块读取CT影像数据、PTV轮廓和OARs轮廓;所述轮廓数据由肿瘤医师根据所述CT影像数据所勾画;
接收肿瘤医师对所述病例设定与剂量指数选择模块的操作,实现病例设定和对应于所设定病例的剂量指数指标的选择,得到患者情况;
采用所述3D剂量分布图预览模块、计算模块对所述CT影像数据和PTV轮廓进行处理,得到处理结果;所述处理结果包括3D剂量分布图、DVH图、所关注的剂量指数和HI/CI指标;
基于所述处理结果,接收肿瘤医师对所述PTV轮廓的修正,使得修正后的PTV轮廓在保证病灶受量和降低危机器官风险之间达到平衡;
通过所述计算模块对修正后的处理结果进行计算,得到靶区轮廓的多目标帕累托最优解集,使得肿瘤医师根据靶区轮廓的多目标帕累托最优解集、结合患者情况和临床要求进行辅助决策,得到靶区勾画结果;
通过所述通信模块将CT影像数据、PTV轮廓和OARs轮廓上传至放疗计划系统,使得物理师进行设备选择、射野布局、约束条件设置、剂量计算、计划优化及预评估,得到放疗计划,由肿瘤医师与物理师对放疗计划进行多次反复调整后交付临床验证和实施。
第二方面,本发明实施例提供了一种剂量分布预览系统,包括:
数据读取与处理模块,用于读取CT影像数据、PTV轮廓和OARs轮廓,以及进行剂量预测模型训练;所述轮廓数据由肿瘤医师根据所述CT影像数据所勾画;
病例设定与剂量指数选择模块,用于接收肿瘤医师对所述病例设定与剂量指数选择模块的操作,实现病例设定和对应于所设定病例的剂量指数指标的选择;
3D剂量分布图预览模块,用于采用训练完成的剂量预测模型,对所述CT影像数据和PTV数据进行处理,得到3D剂量分布图,以供肿瘤医师预览;
计算模块,用于根据所述3D剂量分布图计算并绘制靶区和各危机器官的DVH曲线,并解析出所关注的剂量指数和HI/CI指标,以供肿瘤医师预判评估;
通信模块,用于实现所述剂量分布系统与外部放疗计划系统之间的通信。
作为本申请的一种具体实施方式,所述数据读取与处理模块用于训练剂量预测模型具体过程包括:
收集病例数据;所述病例数据包括放疗模拟定位CT医学影像数据、PTV勾画轮廓数据、OARs勾画轮廓数据、临床剂量分布数据;
对所述病例数据进行格式转化、2D/3D初始转化和图像预处理,并将处理后的病例数据按5:1的比例将所述病例数据随机划分为训练集和验证集;
搭建深度学习网络模型架构;所述深度学习网络模型架构包括生成网络和判别网络;
针对所述生成网络和判别网络分别设置不同的目标函数组合;
采用所述训练集和验证集对所述深度学习网络模型进行训练和测试,得到剂量预测模型。
作为本申请的一种具体实施方式,所述病例设定与剂量指数选择模块具体用于:
接收肿瘤医师设置的病情数据,所述病例数据包括患病部位及具体病种;
根据所述病情数据选取默认的剂量指数。
实施本发明实施例,简化了靶区修改和计划制定流程,提高了放疗靶区确立的合理性和放疗计划质量,缩短了肿瘤患者的治疗延迟时间,降低了工作人员负担的同时提高患者放疗疗效。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍。
图1是本发明实施例提供的剂量分布预览系统的结构图;
图2是剂量分布预览系统开发中的数据预处理流程图;
图3是剂量预测整体深度网络架构示意图;
图4是剂量预测判别网络结构示意图;
图5是剂量预测生成网络结构示意图;
图6是模型训练过程中数据增强方法示意图;
图7是本发明实施例提供的基于剂量分布预览系统指导进行放疗靶区确立和修正流程、放疗计划指定到实施的流程示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例所涉及的缩略语及关键术语定义如下表:
Figure BDA0003213863260000061
Figure BDA0003213863260000071
本发明的发明构思是:近年来深度学习技术在放射治疗领域的应用日益广泛,主要涉及自动诊断、自动分割和剂量预测等方面。本发明先依托大量临床病例数据收集、深度学习方法以及相关软件工程技术开发一套剂量分布预览系统。再基于该系统指导肿瘤医师进行放疗靶区确立和修正,流程大致如下:
首先,基于自主搭建的深度学习(Deep Learning,DL)网络结构开发一套剂量分布预览系统,然后,医师可实时预览并获取基于当前靶区轮廓的剂量分布、HI\CI以及剂量指数(或DVH图)等评估指标情况,以对放疗靶区进行“雕刻式”修改,系统自动存储每次修改后的靶区轮廓及其对应的各项评估指标,并基于多项指标为医师提供靶区轮廓的“帕累托最优(Pareto Optimality)解集”,最终在保证理想治疗剂量和保护危机器官之间寻求平衡。
基于前述发明构思,本发明实施例提供了一种剂量分布预览系统。该系统的开发阶段在整体上大致可分为病例数据收集、数据预处理,深度学习网络模型结构搭建,目标函数设置,模型训练、验证和测试评估,剂量预览系统流程设计等。
如图1所示,该剂量分布预览系统主要包括:
数据读取与处理模块,用于读取CT影像数据、PTV轮廓和OARs轮廓,以及进行剂量预测模型训练;所述轮廓数据由肿瘤医师根据所述CT影像数据所勾画;
病例设定与剂量指数选择模块,用于接收肿瘤医师对所述病例设定与剂量指数选择模块的操作,实现病例设定和对应于所设定病例的剂量指数指标的选择;肿瘤医师设置患病部位(头部、胸部、腹部、盆腔或四肢)以及具体病种(脑瘤、鼻咽癌、食管癌、肺癌、直肠癌、宫颈癌和前列腺癌),系统自动根据病例情况选取默认剂量指数;医师也可根据临床要求或个人偏好设置靶区或不同危机器官的剂量指数指标;
3D剂量分布图预览模块,用于采用训练完成的剂量预测模型,对所述CT影像数据和PTV数据进行处理,得到3D剂量分布图,以供肿瘤医师预览;
计算模块,用于根据所述3D剂量分布图计算并绘制靶区和各危机器官的DVH曲线,并解析出所关注的剂量指数和HI/CI指标,以供肿瘤医师预判评估;
通信模块,用于实现所述剂量分布系统与外部放疗计划系统之间的通信。
其中,训练剂量预测模型具体过程包括:
(1)病例数据收集
·本发明中收集的原始病例数据均为配对数据,即对于某患者病例,应完整收集其相应的DICOM(Digital Imaging and Communications in Medicine)通讯协议医学数据,包括:
·医学影像数据,主要为放疗模拟定位CT医学影像数据,格式为RT Image
·PTV,OARs勾画轮廓信息数据,格式为RT Structure
·临床剂量分布(Dose distribution map)数据,格式为RT Dose
其中,所有原始医学影像数据均来自于放疗科实际临床使用的模拟定位设备,并采取与临床情景一致参数扫描获取;所有勾画轮廓均由至少两位富有经验的肿瘤医师协力完成,并通过高级医师的评估与审核;基于上述勾画轮廓数据所得到的放疗计划均由至少两位富有经验的物理师协力完成,所有计划通过高级物理师和医师审核,并已交付临床使用。本发明中所涉及病种包括脑瘤、鼻咽癌、食管癌、肺癌、直肠癌、宫颈癌和前列腺癌,对于上述所有病种,原始病例收集数量均为100~120例,按照约5:1的比例被随机划分为训练-验证集(Train-Valid Dataset)和测试集(Test Dataset)。
(2)数据预处理
原始数据预处理是深度学习网络训练与验证的前提,所有数据通过TPS导出并进行格式(.dcm/.nii)及2D/3D初始转化,之后进行图像预处理,本发明中医学影像预处理步骤如图2所示,具体包括:
·医学影像重采样(Image Resampling):初始转化后的CT,Structures和Dosedistribution map等3D图像体素间距(Voxel spacing)不一致,设其分别为(SCT_X,SCT_Y,SCT_Z),(SStruc_X,SStruc_Y,SStruc_Z),(SDose_X,SDose_Y,SDose_Z);设其3D图像分辨率分别为(DCT_X,DCT_Y,DCT_Z),(DStruc_X,DStruc_Y,DStruc_Z),(DDose_X,DDose_Y,DDose_Z);设其3D图像大小分别为(FCT_X,FCT_Y,FCT_Z),(FStruc_X,FStruc_Y,FStruc_Z),(FDose_X,FDose_Y,FDose_Z);为保证不同类型3D图像单体素所映射的真实物理空间大小一致,需使其体素间距保持统一为S0,根据在单一维度中图像大小(mm),F=S*D且F为定值,可求得重采样后所有CT,Structures和Dose distribution map的分辨率。本发明中,对于Structures二值化图像(0/255)重采样采取最近邻域从插值(Nearest NeighborInterpolation),对于CT和Dose distribution map采取B-样条插值(B-SplineInterpolation)
·Structures图像设置标签值:针对PTV以及不同OARs和身体外轮廓分别设置不同的标签值L1,L2…Ln,空气部分均保持为0值,使其分别形成对应的二值化图像,其中,标签值上下界与CT灰度区间大致保持一致。
·Dose distribution map像素值-剂量值转化:基于Python-pydicom或ITK等相关工具读取初始Dose distribution map的DoseGridScaling参数,基于此将剂量图像像素值矩阵转化为剂量值矩阵。
·图像裁剪及填充(Image Cropping And Padding):将所有CT,Structures和Dose distribution map图像的3D分辨率裁剪并统一到相同大小,即所有图像矩阵的行、列数均保持一致,为便于网络训练和验证,一般取值128或256,对于实际分辨率小于此值的,以0值对其进行填充。
·图像归一化。将CT,Structures和Dose distribution map图像矩阵按照其对应的区间归一化至[-1,1].
·多通道(Multi-Channels)数据堆叠:本发明中神经网络输入数据采取多通道数据形式输入以提高网络参数量和预测效果。整体多通道数据Channel_1,Channel_2,Channel_3,Channel_4…Channel_n分别为CT影像数据、三维剂量数据、靶区数据和不同OARs二值化图像数据。
(3)深度学习网络模型结构搭建
本发明中所涉及的剂量预测主要基于深度学习图像生成方法,所采用的网络模型结构在整体架构上融合了两种深度学习思想,其一是条件生成对抗网络(ConditionalGenerative Adversarial Network,CGAN)思想,其二是多尺度/多通道特征融合(Multi-scale/channels feature fusion)的编码-解码网络思想。总体网络结构包含两部分,如图3所示,即生成网络和判别网络,两者具有相反的优化目标,因此在每个训练周期中都存在相互对抗、相互促进的过程。在图3中,X为由CT和轮廓数据组成的多通道输入数据,G(X)为生成网络的输出即预测剂量分布图,“Fake”、“Real”代表判别器的分类结果为“假”、“真”。
·判别网络(Discriminator),如图4所示,理想优化目标为可精准区分临床剂量分布图(Clinical dose distribution map)和预测剂量分布(Predicted dosedistribution map),即输入为上述两者时,网络输出分别为“真(Real)”和“假(Fake)”。网络结构采用由卷积/归一/激活层组成的PatchGAN判别网络,其中归一层采用批归一化(Batch Normalize),激活层以LeakyReLU(Leaky Rectified Linear Unit,带泄露线性整流函数)作为激活函数。基于对输入图像重叠块的二分类输出概率值矩阵,设Fn为第n层感受野大小,KS为卷积核尺寸(Kernel Size),S为步长,则第n层感受野大小:
Fn=KS*Fn+1-(KS-S)*(Fn+1-1) (1)
由此可以求得由后往前感受野分别为:4*4*4,7*7*7,16*16*16,34*34*34,70*70*70…,根据所需求感受野大小则可求得判别器卷积层数这一超参。
在图4中,X为由CT和轮廓数据组成的多通道输入数据,G(X)为生成网络的输出即预测剂量分布图,Cli为临床剂量分布图,MCli、MPre分别判别器针对临床剂量分布和预测剂量分布的概率输出矩阵。
·生成网络(Generator),如附图5所示,理想优化目标为输出与临床剂量分布图没有差异的预测剂量分布,使得自身预测结果输入至判别器时输出结果为“真(Real)”。编码部分采取多个卷积/归一/激活层组成的编码块和一系列下采样操作在多尺度提取来自于多通道的图像特征,其中归一层采用批归一化(Batch Normalize),激活层以ReLU(Rectified Linear Unit,线性整流函数)作为激活函数对于每个解码层,其输入都包含对不同尺度、不同通道的特征融合,并对特征图进行卷积/归一/激活处理,解码部分通过上采样实现将高维特征图逐步转变为3D预测剂量分布图。
在图5中,
Figure BDA0003213863260000121
Figure BDA0003213863260000122
分别表示不同等级的编码层和解码层,N@S3表示当前特征图数量为N,三维尺寸大小为S3
(4)目标函数设置
本发明中针对判别网络和生成网络分别设置不同的目标函数组合。
判别网络目标函数主要采用常规的二值交叉熵(Binary Cross Entropy,BCE)和Sigmoid(σ)函数相结合的方式,如下所示:
Figure BDA0003213863260000123
Figure BDA0003213863260000124
Figure BDA0003213863260000125
其中X表示多通道输入数据,Cli和Pre分别代表临床剂量分布和预测剂量分布,LD_Cli和LD_Pre分别表示判别网络对于临床剂量分布和预测剂量分布的损失项;LD表示判别器总的损失函数;MCli和MPre分别为判别网络针对临床剂量分布和预测剂量分布的分类输出概率矩阵。
生成网络目标函数由以下多类损失项按照不同权重组成:
A.对抗损失项(与判别器)(Adversarial loss)
Figure BDA0003213863260000126
其中G(X)表示生成器输出的预测剂量图,MG(X)表示由概率值构成的矩阵,LG_CGAN表示生成网络对抗损失项。
B.近邻体素差异损失项(Adjacent Voxels Difference Loss)
Figure BDA0003213863260000131
其中mT、mC、mS分别表示3D剂量分布图像的横断面、冠状面和矢状面层数,Clii、Clij、、Clik分别表示临床剂量图中某一层的剂量值矩阵,而G(X)i、G(X)j、G(X)k为生成器预测剂量分布图对应的剂量矩阵,LG_AVD表示通过最小化相邻体素剂量值差异来尽量保持图像边缘特征的损失项。
C.L1范数损失项
LG_L1=||Cli-G(X)||1 (7)
其中,LG_L1表示生成网络L1范数损失项。
生成网络总目标函数为:
LG=λCGAN·LG_CGANG_AVD·LG_AVDL1·LG_L1 (8)
其中,λCGAN,λG_AVD,λL1超参数分别表示对抗损失项、近邻体素差异损失项和L1范数损失项的权重。
(5)模型训练、验证和测试评估
所有病例按照总量约5:1比例随机划分为训练-验证集和测试集,采取五折交叉验证(5-Fold Valid)的方式进行模型的训练和验证,输入采取数据增强方法(如图6所示,图中D0表示原始图像尺寸,DX表示随机裁剪后的图像尺寸),每次迭代对多通道图像进行随机旋转和随机裁剪。训练期间不断调整模型超参数,并每隔五个训练周期(epoch)对模型进行存储,以便后期测试,测试后的评估指标主要包括以下方面:
·剂量分布差值图和差值统计直方图
·MAE(平均绝对误差)或MSE(均方误差)
·DVH(Dose-Volume Histograms,剂量体积直方图)
·剂量指数统计:PTV的D99%,D98%,D95%,D50%,D2%,不同OARs的Dmax和Dmean等
·PTV剂量分布的均匀性指数(Homogeneity index,HI):
Figure BDA0003213863260000141
·PTV剂量分布的适形度指数(Conformity index,CI):
Figure BDA0003213863260000142
其中VT表示PTV体积,VP表示处方剂量覆盖体积,VTP表示处方剂量覆盖PTV区域体积
·DICE系数:
Figure BDA0003213863260000143
其中Pre和Cli分别表示预测剂量分布图和临床剂量分布图,Prescription_Dose表示处方剂量
进一步地,上述通信模块具体用于:
将CT影像数据、PTV轮廓和OARs轮廓上传至放疗计划系统,使得物理师进行设备选择、射野布局、约束条件设置、剂量计算、计划优化及预评估,得到放疗计划,由肿瘤医师与物理师对放疗计划进行多次反复调整后交付临床验证和实施。
需要说明的是,本实施例中的剂量分布预览系统基于UDP/IP传输层协议并充分结合DICOM协议实现网络通信功能;基于.Net Core的ASP.Net Core应用开发框架实现客户端。该剂量分布预览系统主要为肿瘤医师提供剂量分布预览、剂量指数(或DVH图)及HI\CI指标计算等功能。
基于上述剂量分布预览系统,请参考图7,本发明实施例提供的放疗靶区确立和修正方法流程如下:
(1)肿瘤医师基于CT或多模态(MRI、PET-CT等)融合的医学影像数据,结合检查或检验结果以及自身临床经验确定GTV位置和区域,并在此基础上外扩得到CTV区域。
(2)在所确定CTV区域的基础上,外扩到初始PTV轮廓,同时勾画OARs轮廓。肿瘤医师可通过前述病例设定与剂量指数选择模块进行病例设定和剂量指数选择。
(3)CT和轮廓数据读取至剂量预览系统,系统处理原始数据,并为肿瘤医师提供与当前PTV所对应的3D剂量分布图、DVH图、所关注的剂量指数和HI\CI等指标计算结果。即,将所述CT影像数据和PTV数据输入3D剂量分布图预览模块进行处理,得到3D剂量分布图;通过所述计算模块根据所述3D剂量分布图计算并绘制靶区和各危机器官的DVH曲线,并解析出所关注的剂量指数和HI/CI指标。
(4)医师参考预览系统并基于当前轮廓数据所得结果,对PTV进行“雕刻式”修正,每次做出细小修正之后均可实时获取三维剂量分布图、DVH和各项剂量指数,尽量使所确定的PTV轮廓在保证病灶(及周围必要区域)受量和降低危机器官风险之间达到平衡。即,基于处理结果(包括3D剂量分布图、DVH图、所关注的剂量指数和HI/CI指标),接收肿瘤医师对所述PTV轮廓的修正,使得修正后的PTV轮廓在保证病灶受量和降低危机器官风险之间达到平衡。
(5)系统保存每次修改轮廓及其所对应的三维剂量分布图、DVH、PTV\OARs的各项剂量指数,并自动基于上述指标给出靶区轮廓的多目标“帕累托最优(Pareto Optimality)解集”,肿瘤医师结合患者情况和临床要求进行决策。即,通过所述计算模块对修正后的处理结果进行计算,得到靶区轮廓的多目标帕累托最优解集,使得肿瘤医师根据靶区轮廓的多目标帕累托最优解集、结合患者情况和临床要求进行辅助决策,得到靶区勾画结果。
(6)将靶区勾画结果交付给高级医师审核。
(7)PTV、OARs及相关影像数据上传至TPS,由物理师进行设备选择、射野布局、约束条件设置、剂量计算和计划优化及预评估等工作,计划设计完成后反馈给医师评估。
(8)肿瘤医师根据物理师反馈的三维剂量分布情况、PTV与OARs的剂量指数(或DVH图)、HI\CI等对计划进行评估,并根据临床规定,要求物理师进一步优化计划,计划需反复调整直至最优状态。
(9)放疗计划经由医师、高级物理师评估和审核后交付临床验证和实施。
从以上描述可以得知,与当前放疗靶区的确立和修正工作方式相比,本发明所提出的基于深度学习的剂量预览系统及其辅助指导下的放疗靶区确立和修正方法具备以下优点:
(1)在肿瘤医师对进行靶区确立和修正过程中,本发明所提出方法针对修改中的靶区轮廓,在无需制定放疗计划的情况下,可提供实时的三维剂量分布预览、相关剂量指数(或DVH图)和HI\CI等指标的提前计算评估功能,以指导肿瘤医师保证处方治疗剂量(对于病灶及其周围必要区域)和降低危机器官(正常组织器官)风险之间找到平衡,使其对靶区轮廓进行“雕刻式”的修改。
(2)本发明所提出方法可使肿瘤医师在靶区确立和修正过程中(尤其是针对病灶体积大、位置特殊或者SIB病例),以系统提供的剂量分布和剂量指数(或DVH图)为标准,从而弥补临床经验的不足和业务水平的差异。
(3)对于放疗物理师,一方面,由于医师基于本发明所提出方法可得到剂量分布、剂量指数、HI\CI等指标情况,因此在每次靶区修改之后,无需物理师进行计划制定以向医师反馈各项评估指标;另一方面,物理师制定放疗计划是基于医师所提供的理想靶区轮廓数据,当计划结果不理想时,可排除靶区勾画的不合理性,从而只关注计划制定过程的失误;因此,本发明所提供方法可在较大程度上减小物理师的工作负担。
(4)本发明所提供的方法在医师和物理师相配合的靶区确立和修正流程所带来的冗余时间,较大程度上减少了计划制定过程中的无用功,有利于缩短患者接受治疗的延迟时间,从而避免延迟治疗所带来的疗效降低和放疗毒副作用风险。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到各种等效的修改或替换,这些修改或替换都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求的保护范围为准。

Claims (10)

1.一种基于剂量分布预览系统的放疗靶区确立和修正方法,其特征在于,所述剂量分布系统包括数据读取与处理模块、病例设定与剂量指数选择模块、3D剂量分布图预览模块、计算模块和通信模块,所述方法包括:
通过所述数据读取与处理模块读取CT影像数据、PTV轮廓和OARs轮廓;所述轮廓数据由肿瘤医师根据所述CT影像数据所勾画;
接收肿瘤医师对所述病例设定与剂量指数选择模块的操作,实现病例设定和对应于所设定病例的剂量指数指标的选择,得到患者情况;
采用所述3D剂量分布图预览模块、计算模块对所述CT影像数据和PTV轮廓进行处理,得到处理结果;所述处理结果包括3D剂量分布图、DVH图、所关注的剂量指数和HI/CI指标;
基于所述处理结果,接收肿瘤医师对所述PTV轮廓的修正,使得修正后的PTV轮廓在保证病灶受量和降低危机器官风险之间达到平衡;
通过所述计算模块对修正后的处理结果进行计算,得到靶区轮廓的多目标帕累托最优解集,使得肿瘤医师根据靶区轮廓的多目标帕累托最优解集、结合患者情况和临床要求进行辅助决策,得到靶区勾画结果;
通过所述通信模块将CT影像数据、PTV轮廓和OARs轮廓上传至放疗计划系统,使得物理师进行设备选择、射野布局、约束条件设置、剂量计算、计划优化及预评估,得到放疗计划,由肿瘤医师与物理师对放疗计划进行多次反复调整后交付临床验证和实施。
2.如权利要求1所述的方法,其特征在于,得到处理结果具体包括:
将所述CT影像数据和PTV数据输入3D剂量分布图预览模块进行处理,得到3D剂量分布图;
通过所述计算模块根据所述3D剂量分布图计算并绘制靶区和各危机器官的DVH曲线,并解析出所关注的剂量指数和HI/CI指标。
3.一种剂量分布预览系统,其特征在于,包括:
数据读取与处理模块,用于读取CT影像数据、PTV轮廓和OARs轮廓,以及进行剂量预测模型训练;所述轮廓数据由肿瘤医师根据所述CT影像数据所勾画;
病例设定与剂量指数选择模块,用于接收肿瘤医师对所述病例设定与剂量指数选择模块的操作,实现病例设定和对应于所设定病例的剂量指数指标的选择;
3D剂量分布图预览模块,用于采用训练完成的剂量预测模型,对所述CT影像数据和PTV数据进行处理,得到3D剂量分布图,以供肿瘤医师预览;
计算模块,用于根据所述3D剂量分布图计算并绘制靶区和各危机器官的DVH曲线,并解析出所关注的剂量指数和HI/CI指标,以供肿瘤医师预判评估;
通信模块,用于实现所述剂量分布系统与外部放疗计划系统之间的通信。
4.如权利要求3所述的系统,其特征在于,所述数据读取与处理模块用于训练剂量预测模型具体过程包括:
收集病例数据;所述病例数据包括放疗模拟定位CT医学影像数据、PTV勾画轮廓数据、OARs勾画轮廓数据、临床剂量分布数据;
对所述病例数据进行格式转化、2D/3D初始转化和图像预处理,并将处理后的病例数据按5:1的比例将所述病例数据随机划分为训练集和验证集;
搭建深度学习网络模型架构;所述深度学习网络模型架构包括生成网络和判别网络;
针对所述生成网络和判别网络分别设置不同的目标函数组合;
采用所述训练集和验证集对所述深度学习网络模型进行训练和测试,得到剂量预测模型。
5.如权利要求4所述的系统,其特征在于,对所述病例数据进行图像预处理具体包括:医学影像重采样、二值化图像标签值设置、剂量分布图像素值-剂量值转化、图像裁剪及填充、图像归一化及多通道数据堆叠。
6.如权利要求4所述的系统,其特征在于,所述判别网络采用由卷积/归一/激活层组成的PatchGAN判别网络。
7.如权利要求4所述的系统,其特征在于,所述判别网络的目标函数采用值交叉熵和Sigmoid(σ)函数相结合的方式;所述生成网络的目标函数由多项损失项按照不同权重组成,多项损失项包括:对抗损失项、近邻体素差异损失项和L1范数损失项。
8.如权利要求4所述的系统,其特征在于,采用所述训练集和验证集对所述深度学习网络模型进行训练和测试,得到剂量预测模型,具体包括:
采取五折交叉验证的方式进行模型的训练和验证,输入采取数据增强方法,每次迭代对多通道图像进行随机旋转和随机裁剪,训练期间不断调整模型超参数,并每隔五个训练周期对模型进行存储,以便后期测试。
9.如权利要求4所述的系统,其特征在于,所述病例设定与剂量指数选择模块具体用于:
接收肿瘤医师设置的病情数据,所述病例数据包括患病部位及具体病种;
根据所述病情数据选取默认的剂量指数。
10.如权利要求9所述的系统,其特征在于,所述通信模块具体用于:
将CT影像数据、PTV轮廓和OARs轮廓上传至放疗计划系统,使得物理师进行设备选择、射野布局、约束条件设置、剂量计算、计划优化及预评估,得到放疗计划,由肿瘤医师与物理师对放疗计划进行多次反复调整后交付临床验证和实施。
CN202110938598.9A 2021-08-16 2021-08-16 基于剂量分布预览系统的放疗靶区确立和修正方法 Pending CN113674834A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110938598.9A CN113674834A (zh) 2021-08-16 2021-08-16 基于剂量分布预览系统的放疗靶区确立和修正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110938598.9A CN113674834A (zh) 2021-08-16 2021-08-16 基于剂量分布预览系统的放疗靶区确立和修正方法

Publications (1)

Publication Number Publication Date
CN113674834A true CN113674834A (zh) 2021-11-19

Family

ID=78543119

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110938598.9A Pending CN113674834A (zh) 2021-08-16 2021-08-16 基于剂量分布预览系统的放疗靶区确立和修正方法

Country Status (1)

Country Link
CN (1) CN113674834A (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114344737A (zh) * 2022-01-07 2022-04-15 北京大学第一医院 一种肿物放射治疗控制系统及存储介质
CN114796901A (zh) * 2022-05-30 2022-07-29 北京大学第一医院 一种腰骶神经根的自动勾画方法、设备及存储介质
CN114904164A (zh) * 2022-05-10 2022-08-16 苏州雷泰医疗科技有限公司 基于人工智能的放疗自动计划生成方法及装置
CN115300811A (zh) * 2022-08-08 2022-11-08 中山大学肿瘤防治中心(中山大学附属肿瘤医院、中山大学肿瘤研究所) 一种基于机器学习的剂量分布确定方法及装置
CN117577275A (zh) * 2024-01-17 2024-02-20 福建自贸试验区厦门片区Manteia数据科技有限公司 基于壳结构的剂量信息确定装置、电子设备及存储介质
CN117695534A (zh) * 2023-12-14 2024-03-15 陕西省人民医院(陕西省临床医学研究院) 一种宫颈癌放疗精准定位与剂量控制系统
CN118039062A (zh) * 2024-04-12 2024-05-14 四川省肿瘤医院 一种基于大数据分析的个体化化疗剂量远程控制方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105825073A (zh) * 2016-06-17 2016-08-03 张家港赛提菲克医疗器械有限公司 一种在线放疗计划质量控制软件
CN107441637A (zh) * 2017-08-30 2017-12-08 南方医科大学 调强放疗计划中三维剂量分布的预测方法及其应用
CN111462916A (zh) * 2020-04-03 2020-07-28 戴建荣 一种放疗计划危及器官剂量体积直方图的预测方法和装置
CN112635024A (zh) * 2021-03-10 2021-04-09 四川大学 一种放射治疗自动计划设计系统及其构建方法
CN113241151A (zh) * 2021-02-25 2021-08-10 中科超精(南京)科技有限公司 一种智慧放疗云系统及实现方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105825073A (zh) * 2016-06-17 2016-08-03 张家港赛提菲克医疗器械有限公司 一种在线放疗计划质量控制软件
CN107441637A (zh) * 2017-08-30 2017-12-08 南方医科大学 调强放疗计划中三维剂量分布的预测方法及其应用
CN111462916A (zh) * 2020-04-03 2020-07-28 戴建荣 一种放疗计划危及器官剂量体积直方图的预测方法和装置
CN113241151A (zh) * 2021-02-25 2021-08-10 中科超精(南京)科技有限公司 一种智慧放疗云系统及实现方法
CN112635024A (zh) * 2021-03-10 2021-04-09 四川大学 一种放射治疗自动计划设计系统及其构建方法

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114344737A (zh) * 2022-01-07 2022-04-15 北京大学第一医院 一种肿物放射治疗控制系统及存储介质
CN114904164A (zh) * 2022-05-10 2022-08-16 苏州雷泰医疗科技有限公司 基于人工智能的放疗自动计划生成方法及装置
CN114796901A (zh) * 2022-05-30 2022-07-29 北京大学第一医院 一种腰骶神经根的自动勾画方法、设备及存储介质
CN115300811A (zh) * 2022-08-08 2022-11-08 中山大学肿瘤防治中心(中山大学附属肿瘤医院、中山大学肿瘤研究所) 一种基于机器学习的剂量分布确定方法及装置
CN115300811B (zh) * 2022-08-08 2024-01-05 中山大学肿瘤防治中心(中山大学附属肿瘤医院、中山大学肿瘤研究所) 一种基于机器学习的剂量分布确定方法及装置
CN117695534A (zh) * 2023-12-14 2024-03-15 陕西省人民医院(陕西省临床医学研究院) 一种宫颈癌放疗精准定位与剂量控制系统
CN117695534B (zh) * 2023-12-14 2024-05-28 陕西省人民医院(陕西省临床医学研究院) 一种宫颈癌放疗精准定位与剂量控制系统
CN117577275A (zh) * 2024-01-17 2024-02-20 福建自贸试验区厦门片区Manteia数据科技有限公司 基于壳结构的剂量信息确定装置、电子设备及存储介质
CN117577275B (zh) * 2024-01-17 2024-04-19 福建自贸试验区厦门片区Manteia数据科技有限公司 基于壳结构的剂量信息确定装置、电子设备及存储介质
CN118039062A (zh) * 2024-04-12 2024-05-14 四川省肿瘤医院 一种基于大数据分析的个体化化疗剂量远程控制方法

Similar Documents

Publication Publication Date Title
CN113674834A (zh) 基于剂量分布预览系统的放疗靶区确立和修正方法
CN114681813B (zh) 放射治疗自动计划系统、自动计划方法及存储介质
US20210056687A1 (en) Systems and methods for segmentation of intra-patient medical images
CN110197709B (zh) 一种基于深度学习与先验计划的三维剂量预测方法
CN108717866B (zh) 一种预测放疗计划剂量分布的方法、装置、设备及存储介质
CN110211664B (zh) 一种基于机器学习自动设计放射治疗方案的系统
US11964170B2 (en) Standardized artificial intelligence automatic radiation therapy planning method and system
CN111028914B (zh) 人工智能引导的剂量预测方法与系统
CN108815721B (zh) 一种照射剂量确定方法及系统
US20230020911A1 (en) Methods and systems for adaptive radiotherapy treatment planning using deep learning engines
CN108778416A (zh) 使用组织参数估计根据mr数据的伪ct生成
CN108770373A (zh) 使用特征回归模型根据mr数据的伪ct生成
CN110085298A (zh) 基于深度网络学习的调强放疗计划三维剂量分布预测方法
CN112820377B (zh) 一种基于深度学习的放疗计划自动生成方法
CN115829972B (zh) 一种放疗计划三维剂量分布与通量的同步预测方法及装置
Wang et al. DeepEC: An error correction framework for dose prediction and organ segmentation using deep neural networks
CN113178242B (zh) 一种基于耦合式生成对抗网络的自动计划优化系统
CN110706779A (zh) 一种放射治疗计划精确目标函数的自动生成方法
CN111888665B (zh) 基于自适应对抗网络的三维剂量分布预测模型的构建方法
Kallis et al. Automated treatment planning framework for brachytherapy of cervical cancer using 3D dose predictions
US20230065196A1 (en) Patient-specific organ dose quantification and inverse optimization for ct
CN113941100A (zh) 根据三维空间剂量分布生成可交付放疗计划的方法及装置
US20240001149A1 (en) Methods and systems for dose distribution prediction
US20220241614A1 (en) Fluence map prediction and treatment plan generation for automatic radiation treatment planning
Ni et al. Fast VMAT planning for prostate radiotherapy: dosimetric validation of a deep learning-based initial segment generation method

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
TA01 Transfer of patent application right

Effective date of registration: 20230906

Address after: 250014 No. 1301, unit 1, building 9, No. 89, Jingshi Road, Lixia District, Jinan City, Shandong Province

Applicant after: Yu Jinming

Applicant after: Affiliated Tumor Hospital of Shandong First Medical University (Shandong cancer prevention and treatment institute Shandong Cancer Hospital)

Address before: 250014 No. 1301, unit 1, building 9, No. 89, Jingshi Road, Lixia District, Jinan City, Shandong Province

Applicant before: Yu Jinming

Applicant before: Beijing Huaguang Putai Biotechnology Co.,Ltd.

TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20231007

Address after: 201807 2258 Chengbei Road, Jiading District, Shanghai

Applicant after: Shanghai Lianying Medical Technology Co.,Ltd.

Address before: 250014 No. 1301, unit 1, building 9, No. 89, Jingshi Road, Lixia District, Jinan City, Shandong Province

Applicant before: Yu Jinming

Applicant before: Affiliated Tumor Hospital of Shandong First Medical University (Shandong cancer prevention and treatment institute Shandong Cancer Hospital)

TA01 Transfer of patent application right