CN113499091B - 一种患者体表和体内肿瘤运动相关性及肿瘤内部动度的预测方法和系统 - Google Patents
一种患者体表和体内肿瘤运动相关性及肿瘤内部动度的预测方法和系统 Download PDFInfo
- Publication number
- CN113499091B CN113499091B CN202110956747.4A CN202110956747A CN113499091B CN 113499091 B CN113499091 B CN 113499091B CN 202110956747 A CN202110956747 A CN 202110956747A CN 113499091 B CN113499091 B CN 113499091B
- Authority
- CN
- China
- Prior art keywords
- tumor
- patient
- scanning
- roi
- feature
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 206010028980 Neoplasm Diseases 0.000 title claims abstract description 137
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000000605 extraction Methods 0.000 claims abstract description 11
- 230000000241 respiratory effect Effects 0.000 claims description 27
- 238000001727 in vivo Methods 0.000 claims description 23
- 230000029058 respiratory gaseous exchange Effects 0.000 claims description 23
- 238000010801 machine learning Methods 0.000 claims description 22
- 238000010187 selection method Methods 0.000 claims description 16
- 238000004422 calculation algorithm Methods 0.000 claims description 13
- 238000013528 artificial neural network Methods 0.000 claims description 12
- 238000012549 training Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 8
- 210000000056 organ Anatomy 0.000 claims description 7
- 210000000779 thoracic wall Anatomy 0.000 claims description 7
- 210000002569 neuron Anatomy 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- 230000000877 morphologic effect Effects 0.000 claims description 5
- 238000012216 screening Methods 0.000 claims description 5
- 210000000038 chest Anatomy 0.000 claims description 4
- 238000002591 computed tomography Methods 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims description 4
- 210000001015 abdomen Anatomy 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 claims description 3
- 201000011510 cancer Diseases 0.000 claims description 2
- 238000000546 chi-square test Methods 0.000 claims description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 claims description 2
- 230000008030 elimination Effects 0.000 claims description 2
- 238000003379 elimination reaction Methods 0.000 claims description 2
- 238000010606 normalization Methods 0.000 claims description 2
- 238000012847 principal component analysis method Methods 0.000 claims description 2
- 238000007637 random forest analysis Methods 0.000 claims description 2
- 230000000391 smoking effect Effects 0.000 claims description 2
- 239000000126 substance Substances 0.000 claims description 2
- 238000012706 support-vector machine Methods 0.000 claims description 2
- 230000037230 mobility Effects 0.000 description 20
- 238000001959 radiotherapy Methods 0.000 description 11
- 238000011156 evaluation Methods 0.000 description 6
- 238000007726 management method Methods 0.000 description 6
- 230000003287 optical effect Effects 0.000 description 5
- 230000035945 sensitivity Effects 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- 208000014018 liver neoplasm Diseases 0.000 description 4
- 208000020816 lung neoplasm Diseases 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 206010019695 Hepatic neoplasm Diseases 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 208000037841 lung tumor Diseases 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 230000004899 motility Effects 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 208000006678 Abdominal Neoplasms Diseases 0.000 description 2
- 206010058467 Lung neoplasm malignant Diseases 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 210000004185 liver Anatomy 0.000 description 2
- 210000004072 lung Anatomy 0.000 description 2
- 201000005202 lung cancer Diseases 0.000 description 2
- 230000005019 pattern of movement Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 206010011224 Cough Diseases 0.000 description 1
- 206010016654 Fibrosis Diseases 0.000 description 1
- 206010037394 Pulmonary haemorrhage Diseases 0.000 description 1
- 238000012952 Resampling Methods 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007882 cirrhosis Effects 0.000 description 1
- 208000019425 cirrhosis of liver Diseases 0.000 description 1
- 238000007408 cone-beam computed tomography Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000002594 fluoroscopy Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000012905 input function Methods 0.000 description 1
- 201000007270 liver cancer Diseases 0.000 description 1
- 230000004199 lung function Effects 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 230000001394 metastastic effect Effects 0.000 description 1
- 206010061289 metastatic neoplasm Diseases 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 208000002154 non-small cell lung carcinoma Diseases 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000004393 prognosis Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 208000029729 tumor suppressor gene on chromosome 11 Diseases 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/46—Arrangements for interfacing with the operator or the patient
- A61B6/461—Displaying means of special interest
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
-
- 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/1001—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy using radiation sources introduced into or applied onto the body; brachytherapy
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1048—Monitoring, verifying, controlling systems and methods
- A61N5/1049—Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/40—Analysis of texture
- G06T7/41—Analysis of texture based on statistical description of texture
- G06T7/42—Analysis of texture based on statistical description of texture using transform domain methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/40—Analysis of texture
- G06T7/41—Analysis of texture based on statistical description of texture
- G06T7/44—Analysis of texture based on statistical description of texture using image operators, e.g. filters, edge density metrics or local histograms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/40—Analysis of texture
- G06T7/41—Analysis of texture based on statistical description of texture
- G06T7/45—Analysis of texture based on statistical description of texture using co-occurrence matrix computation
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10076—4D tomography; Time-sequential 3D tomography
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20081—Training; Learning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20084—Artificial neural networks [ANN]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30061—Lung
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30096—Tumor; Lesion
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biomedical Technology (AREA)
- Medical Informatics (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Theoretical Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Pathology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Optics & Photonics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Heart & Thoracic Surgery (AREA)
- Probability & Statistics with Applications (AREA)
- High Energy & Nuclear Physics (AREA)
- Biophysics (AREA)
- Quality & Reliability (AREA)
- Human Computer Interaction (AREA)
- Pulmonology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Image Analysis (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明提供了一种患者体表和体内肿瘤运动相关性及肿瘤内部动度的预测方法和系统。通过病人临床信息和CT数据采集、4DCT图像重建、影响组学特征提取，经机器学习模型预测患者体表和体内肿瘤运动相关性系数及肿瘤AP、SI、LR方向的内部动度，可为临床决策提供参考，具有很好的应用前景。
Description
技术领域
本发明属于放疗影像学领域，具体涉及一种患者体表和体内肿瘤运动相关性及肿瘤内部动度的预测方法和系统。
背景技术
肺癌和肝癌均为世界范围内的高发癌症，放射治疗是治疗非小细胞肺癌占肺癌诊断85％、原发性和转移性肝肿瘤的有效手段，肺、肝等胸腹部肿瘤受呼吸运动影响位置不断变化，体积较小的肿瘤区甚至有“脱靶”风险，给放疗精度带来巨大影响，CBCT辅助摆位能减少分次间位置变化，但不能在治疗中提供实时成像。肿瘤位置监测和管理对于保证肿瘤治理效果、保护正常器官具有重要意义。
为了提高治疗精确性，不同的肿瘤位置实时监控方法被提出。呼吸运动门控通过在特定的呼吸时相进行照射来减少正常组织剂量，深吸气屏气(Deep inspirationbreath-hold)通过减少/消除患者呼吸运动提高放疗的精度，由于呼吸控制受到患者耐受性的影响，多数胸部肿瘤患者肺功能较差，伴有咳嗽、身体虚弱等状况无法配合；植入标记物的肿瘤位置实时监测允许加速器自动根据探测到的位置反馈持续调节照射束(通过MLC、治疗床位置)，精度最高，但将对患者造成侵入性伤害，并伴随肺出血、标记物移位等副作用风险。非侵入性实时成像技术也被广泛研究，基于荧光透视的实时成像效果理想，但病人将暴露于额外剂量。患者解剖结构、体表标记和人工神经网络也被用于实时肿瘤跟踪，在一定程度上降低了肿瘤跟踪误差的发生率。然而，当肿瘤边界模糊或被骨性结构遮盖时，上述方法通常会失败、降低治疗精确性。
为了开发理想的高精度无标记实时肿瘤跟踪方法，腹部肿瘤体表标记与体内肿瘤之间的呼吸运动相关性被研究，DAVID报告了体表标记与体内肿瘤之间运动频率的一致性(误差小于等于6mm)。美国医学物理学家协会报告指出肺肿瘤的运动模式与位置密切相关，位于下叶且不依附刚性结构(如胸壁、椎骨)的肿瘤运动幅度明显高于其他位置，对于肝肿瘤而言，肿瘤位置、是否肝硬化、肝脏手术史都在不同程度影响肿瘤的运动模式。基于以上原因，体表标记与体内肿瘤运动之间的相关系数无法一概而论，对于指定患者的相关性高低判定至关重要。
影像组学特征提取能够从CT、MRI等影像中提取感兴趣区的统计特征和组织异质性、表形等高通量纹理特征，避免主观性且省时，结合机器学习对这些数据信息进行深层次的分析，找到某些特征在不同群体上的分布规律，从而实现肿瘤的辅助诊断、放射损伤评估、预后预测等，在放射肿瘤学研究中广泛应用。肺、肝肿瘤患者的4DCT分为十个呼吸时相CT、平均密度和最大密度投影CT，包含肿瘤的运动信息，利用影像组学提取、筛选可得到包含肿瘤运动信息的相关特征，这些特征可以用来研究肿瘤运动。目前尚无利用影像组学特征进行肿瘤运动与体表变化相关性或体内肿瘤动度的研究。
发明内容
本发明的目的在于提供一种患者体表和体内肿瘤运动相关性及肿瘤内部动度的预测方法和系统。
本发明提供了一种患者肿瘤内部动度的预测方法，包括如下步骤：
(1)数据采集：对患者肿瘤部位进行CT扫描，扫描方式为电影扫描，得到10个扫描周期的扫描数据，每个扫描周期包含一个完整的呼吸周期；
(2)数据处理：将步骤(1)得到的扫描数据进行4DCT图像重建，得到与所述的10个扫描周期的扫描数据对应的10个呼吸时相的CT，对步骤(1)获得的扫描数据进行平均密度投影，得到AIP CT；
(3)勾画ROI：在步骤(2)的10个呼吸时相的CT勾画ROI信息，结合AIP CT得到包含运动信息的ROI；
(4)特征提取：从步骤(3)得到的ROI中提取影像组学特征，所述特影像组学特征包括描述肿瘤3D物理外观的形态学特征、描述ROI强度的一阶分布的强度特征、原始纹理特征、高阶特征；
(5)将影像组学特征和患者的临床信息特征输入机器学习模型处理，得到预测值；所述预测值包括：肿瘤在体内AP、SI、LR三个方向的最大运动值。
进一步地，步骤(4)所述的描述ROI强度的一阶分布的强度特征包括描述ROI区域内体素强度的一阶统计特征值；所述高阶特征包括小波转换特征和高斯滤波转换特征。
进一步地，步骤(5)所述的临床信息特征包括：性别、年龄、体重、吸烟史、癌症病种、TNM分期、ITV体积、目标患瘤脏器体积、肿瘤质心至前胸壁距离、肿瘤质心至后胸壁距离、患瘤脏器AP方向维度、胸腹部肿瘤SI方向维度、肿瘤质心至肿瘤所在脏器边缘距离、肿瘤质心至隔膜距离、肿瘤是否靠胸壁、肿瘤位置。
更进一步地，步骤(5)所述机器学习模型通过如下方法构建而成：
1)基本模型的构建：
(a)数据预处理：对不属于同一量纲的影像组学特征进行无量纲化，剔除冗余信息，将定性特征进行哑编码转化为定量特征，并进行归一化处理；
公式如下：
上式fnor表示归一化的特征,f表示未处理的特征,fmean表示每类特征的平均值，SD表示每类特征的标准差；
(b)特征筛选：通过特征选择方法筛选出与AP、SI、LR三方向的肿瘤内部动度相关系数最高的特征；
(c)基于步骤(b)筛选出的特征和患者的临床信息特征构建得到机器学习模型，并采用机器学习算法对机器学习模型进行训练。
更进一步地，步骤(b)所述特征选择方法是：单变量特征选择方法、方差选择法、卡方检验、递归特征消除法、基于惩罚项的特征选择法、基于树模型的特征选择法或主成分分析法，优选为单变量特征选择方法。
更进一步地，步骤(c)所述机器学习算法为神经网络算法、支持向量机、随机森林、XGBoost或广义线性模型，优选为神经网络算法。
更进一步地，上述神经网络算法包括分支网络和神经元，所述分支网络包括手工特征网络、影像组学特征网络、全特征网络；所述手工特征网络学习患者的临床信息特征，所述影像组学特征网络学习影像组学特征；所述全特征网络学习临床信息特征与影像组学特征之间的知识；所述神经元学习各分支网络的权重。
本发明还提供了一种患者体表和体内肿瘤运动相关性的预测方法，包括如下步骤：
1)体表数据采集：使用Varian实时位置管理系统通过红外摄像机跟踪安装在患者体表的红外反射标记，获取患者自由呼吸状态下AP方向的体表变化，得到10个呼吸时相对应的AP方向体表运动值；
2)肿瘤体内动度分析：采用权利要求1～6任一项所述的预测方法预测肿瘤内部动度，得到10个呼吸时相对应的肿瘤在体内AP、SI、LR三个方向的运动值；
3)将步骤1)得到的AP方向体表运动值与步骤2)得到的肿瘤在体内AP、SI、LR三个方向的最大运动值通过如下公式计算，得到预测的患者体表和体内肿瘤AP、SI、LR三个方向运动值的spearman相关性系数：
xi表示肿瘤在第i个呼吸时相对应的AP方向的体表运动值，表示AP方向所有呼吸时相对应的体表运动值的平均值；
yi表示第i个呼吸时相对应的肿瘤在体内AP、SI或LR方向的运动值，表示所有呼吸时相对应的肿瘤在体内AP、SI或LR方向的运动值的平均值；i为呼吸时相的序号，i为1到10的整数。
本发明还提供了一种患者体表和体内肿瘤运动相关性及肿瘤内部动度的预测系统，包括如下模块：
数据采集和输入模块：对患者肿瘤部位进行CT扫描，扫描方式为电影扫描，得到10个扫描周期的扫描数据，每个扫描周期包含一个完整的呼吸周期，输入患者的临床信息特征；
数据处理模块：将数据采集模块获得的扫描数据进行4DCT图像重建，得到与所述的10个扫描周期的扫描数据对应的10个呼吸时相的CT，并对数据采集模块获得的扫描数据进行平均密度投影，得到AIP CT；
ROI勾画模块：对数据处理模块获得的10个呼吸时相的CT勾画ROI信息，结合AIPCT得到包含运动信息的ROI；
特征提取模块：从ROI勾画模块获得的ROI中提取影像组学特征，所述特影像组学特征包括描述肿瘤3D物理外观的形态学特征、描述ROI强度的一阶分布的强度特征、原始纹理特征、高阶特征；
预测模块：将数据采集和输入模块获得的患者的临床信息特征，以及特征提取模块获得的影像组学特征输入机器学习模型，机器学习模型输出预测值；所述预测值包括：肿瘤在体内AP、SI、LR三个方向的最大运动值。
本发明还提供了一种介质，所述介质存储有计算机可执行指令，所述指令在被处理器执行时可用于实现上述预测肿瘤内部动度及体表运动和体内肿瘤运动之间的相关性的方法。
本发明特别针对胸腹部肿瘤的动度以及体表运动和体内肿瘤运动之间的相关性预测，尤其是肺部肿瘤的预测。
本方法可以实现准确预测肿瘤在体内AP、SI、LR三个方向的内部动度，准确预测体表运动和体内肿瘤运动之间的相关性，为放射肿瘤医生确定是否对肿瘤患者进行主动运动管理策略提供参考。
对预测得到的肿瘤内部动度和控制限值进行比较，如果超过控制限值，说明肿瘤体内动度较大，建议采用主动呼吸控制技术；如果小于等于控制限值，说明肿瘤体内动度较小，建议可不采用主动呼吸控制技术。
对预测得到的体表运动和体内肿瘤运动的相关性系数和控制限值进行比较，如果超过控制限值，说明患者体表运动和体内肿瘤运动具有较强相关性，建议采用光学体表引导放疗技术；如果该相关性系数小于等于控制限值，说明患者体表运动和体内肿瘤运动具有较小相关性，不建议采用光学体表引导放疗技术。
显然，根据本发明的上述内容，按照本领域的普通技术知识和惯用手段，在不脱离本发明上述基本技术思想前提下，还可以做出其它多种形式的修改、替换或变更。
以下通过实施例形式的具体实施方式，对本发明的上述内容再作进一步的详细说明。但不应将此理解为本发明上述主题的范围仅限于以下的实例。凡基于本发明上述内容所实现的技术均属于本发明的范围。
附图说明
图1为本发明神经网络的结构示意图。
图2为本发明的预测流程图。
具体实施方式
本发明实施例的CT模拟机是GE Healthcare(Waukesha,WI)，所用实时位置管理系统是瓦里安公司的RPM系统。放疗计划由Raystation计划系统(Raystation MedicalLaboratories AB,Stockholm,Sweden)创建、计算并导出，并在配备Agility多叶准直器的医科达Versa HDTM(Elekta,Crawley,UK)医用电子直线加速器上进行照射。
实施例1、本发明预测方法
1、数据采集：患者在自由呼吸状态下使用CT模拟机(GE Healthcare；Waukesha,WI)进行“电影”模式扫描，每个床位扫描持续时间为患者的一个呼吸周期加一秒，同时使用Varian实时位置管理(RPM)系统通过红外摄像机跟踪安装在患者体表的红外反射标记，获取患者自由呼吸状态下AP方向的体表变化。
2、数据预处理：扫描后对所得CT图像进行4DCT重建，重建厚度为3毫米，平面分辨率为0.9*0.9毫米，得到包含完整呼吸周期的十个时相CT。再对扫描数据进行平均密度投影，生成AIP CT。
3、ROI勾画：首先由一个放射肿瘤医生使用Raystation计划系统(TPS)V4.7.6(RaySearch Laboratories AB,Stockholm,Sweden)分别在十个呼吸时相CT上手动勾画感兴趣区(ROI)，然后在AIP CT上融合以上十个ROI轮廓生成包含运动信息的ROI，并审查，编辑，最后由具有10多年经验的放射肿瘤学家批准。
4、特征提取：提取高通量的特征是影像组学研究的基础，是将影像数据与临床指标联系起来的桥梁。使用pyradiomics v3.0.1(http://www.radiomics.io/pyradiomics.html)从勾画好的ROI中提取影像组学特征，YAML配置文件参数的设置参考pyradiomics开发示例(https://github.com/AIM-Harvard/pyradiomics/tree/master/examples/exampleSe ttings)，重采样网格为[1,1,1]，箱宽设置为25Hounsfield以提高模式灵敏度、降低图像噪声，voxelArrayShift为1000，提取特征参数共四类，包括描述肿瘤3D物理外观的形态学特征14个；描述ROI强度的一阶分布的强度特征18个，描述ROI区域内体素强度的一阶统计特征值如平均值、最大值和最小值等；原始纹理特征68个，用四种不同方法描述了体素强度之间的空间关系，包括灰度共生矩阵(GLCM)、灰度游程矩阵(GLRLM)、灰度区域大小矩阵(GLSZM)、灰度相关矩阵(GLDM)；高阶特征1118个(包括小波转换特征688个和高斯滤波转换特征430个)。
5、将提取出的影像组学特征和临床信息输入机器学习模型。机器学习模型构建包括如下步骤：
5.1数据预处理，包括对不属于同一量纲的数据进行无量纲化，剔除冗余信息，将定性特征进行哑编码转化为定量特征，对数据进行标准化，公式如下。
5.2特征选择，采用单一变量的相关性特征选择法，剔除冗余特征。具体步骤为计算各个特征对目标值的相关系数以及相关系数的p值，并且选择得到相关系数较高的特征。
5.3采用神经网络算法，对输入特征进行训练，进行超参数调整，以获得最优的预测模型。
a)模型结构：本方法采用了两种特征进行训练，包括手工提取的临床特征和影像组学特征，本方法采用了神经网络进行训练。
通过手工特征网络学习手工提取的临床特征，通过影像组学特征网络学习CT的影像组学特征，通过全特征网络学习两种特征之间的知识。之后，神经元被设计用来学习每个分支网络的权重来进行预测。
F:fNN(X)×WNN→P
fNN表示神经网络，WNN表示fNN的权重。手工临床特征和CT影像组学特征用作输入层。隐藏层用于将输入功能映射到高维空间以进行完整表达，然后压缩到低维空间进行预测。如图1所示，隐藏的网络层有3层。网络中的神经元数量分别为512、1024和128。
b)模型参数：使用均方误差(MSE)作为成本函数。优化结果为随机梯度下降(SGD)，学习率为0.001。学习率以每5个循环乘以0.98的速度衰减。批次大小设置为200。为了减少过拟合的副作用，在最后一个隐藏层中设置了dropout，概率为0.6。
6、将临床信息特征(见表1)和基于CT的影像组学特征输入步骤5所得到模型预测预测肿瘤在SI、LR、AP三个方向的内部动度(如图2所示)。
7、将步骤1得到的AP方向体表运动值与步骤6得到的肿瘤在体内AP、SI、LR三个方向的运动值通过如下公式计算，得到预测的患者体表和体内肿瘤AP、SI、LR三个方向运动值的spearman相关性系数：
xi表示肿瘤在第i个呼吸时相时AP方向的体表最大运动值，表示AP方向所有呼吸时相体表最大运动值的平均值；
yi表示第i个呼吸时相肿瘤在体内AP、SI或LR方向的最大运动值，表示所有呼吸时相肿瘤在体内AP、SI或LR方向的最大运动值的平均值；i为呼吸时相的序号，i为1到10的整数。
对预测得到的肿瘤内部动度和控制限值(控制限值由临床医师根据临床标准设定，一般标准为8mm)进行比较，如果超过控制限值，说明肿瘤体内动度较大，建议采用主动呼吸控制技术，请患者协助控制呼吸以提高放疗精度；如果小于等于控制限值，说明肿瘤体内动度较小，建议可不采用主动呼吸控制技术。
对预测得到的体表运动和体内肿瘤运动的相关性系数和控制限值(本发明将控制限值定为0.7)进行比较，如果超过控制限值(spearman相关性系数大于等于0.7)，说明患者体表运动和体内肿瘤运动具有较强相关性，建议采用光学体表引导放疗技术；如果该相关性系数小于等于控制限值，说明患者体表运动和体内肿瘤运动具有较小相关性，建议不采用光学体表引导放疗技术。
表1.临床特征汇总
以下通过实验例证明本发明的有益效果。
实验例1、本发明预测系统的预测准确性
1、实验方法
(1)本发明方法预测患者肿瘤内部动度的准确度评估
实验数据共有165例，随机分为126例的训练集和39例的测试集，训练集用于机器学习模型的训练，测试集用于实施例1实验结果的评估。
上述方法的平均绝对误差视为预测误差分别进行计算，公式如下：
其中，N表示测试集的数量，pi表示模型预测的肿瘤动度，yi表示使用4DCT获取的实际肿瘤动度。
(2)本发明方法预测患者体表和体内肿瘤运动相关性的准确性评估
实验数据共有181例，随机分为86例的训练集和95例的测试集，训练集用于机器学习模型的训练，测试集用于实施例1实验结果的评估。
上述方法的评价指标为敏感性(Sensitivity)、特异性(Specificity)，公式如下：
其中，N表示测试集数量，TP表示预测为体表体内弱相关性(预测得到的spearman相关系数小于0.7)且实际也为弱相关性(实际spearman相关系数也小于0.7)的类别，TN表示预测结果为强相关性(预测得到的spearman相关系数大于等于0.7)且实际也为强相关性(实际spearman相关系数也大于等于0.7)的类别，FP表示预测为强相关性而实际为弱相关性的类别。
2、实验结果
(1)本发明方法预测患者肿瘤内部动度的准确度评估
本发明在测试集中实现的平均预测误差仅为1.23mm，误差很小，说明本发明方法预测患者肿瘤内部动度准确度高。
(2)本发明方法预测患者体表和体内肿瘤运动相关性的准确性评估结果见表2。
表2预测患者体表和体内肿瘤运动相关性的准确性评估结果
| 评价指标 | 实施例1 |
| Sensitivity | 0.788 |
| Specificity | 0.930 |
从表2可以看出，本发明方法预测患者体表和体内肿瘤运动相关性的敏感性和特异性高，说明预测准确度高。
综上，本发明提供了一种患者体表和体内肿瘤运动相关性及肿瘤内部动度的预测方法和系统，可以准确预测肿瘤在体内AP、SI、LR三个方向的内部动度，为放射肿瘤医生确定是否对肿瘤患者进行主动运动管理策略提供参考，还可准确预测体表运动和体内肿瘤运动之间的相关性，为放射肿瘤医生决策是否对患者使用光学表面引导放疗提供参考。
Claims (12)
1.一种患者肿瘤内部动度的预测方法，其特征在于，包括如下步骤：
(1)数据采集：对患者肿瘤部位进行CT扫描，扫描方式为电影扫描，得到10个扫描周期的扫描数据，每个扫描周期包含一个完整的呼吸周期；
(2)数据处理：将步骤(1)得到的扫描数据进行4DCT图像重建，得到与所述的10个扫描周期的扫描数据对应的10个呼吸时相的CT，对步骤(1)获得的扫描数据进行平均密度投影，得到AIP CT；
(3)勾画ROI：在步骤(2)的10个呼吸时相的CT勾画ROI信息，结合AIP CT得到包含运动信息的ROI；
(4)特征提取：从步骤(3)得到的ROI中提取影像组学特征，所述影像组学特征包括描述肿瘤3D物理外观的形态学特征、描述ROI强度的一阶分布的强度特征、原始纹理特征、高阶特征；
(5)将影像组学特征和患者的临床信息特征输入机器学习模型处理，得到预测值；所述预测值包括：肿瘤在体内AP、SI、LR三个方向的最大运动值。
2.如权利要求1所述的预测方法，其特征在于，步骤(4)所述的描述ROI强度的一阶分布的强度特征包括描述ROI区域内体素强度的一阶统计特征值；所述高阶特征包括小波转换特征和高斯滤波转换特征。
3.如权利要求1所述的预测方法，其特征在于，步骤(5)所述的临床信息特征包括：性别、年龄、体重、吸烟史、癌症病种、TNM分期、ITV体积、目标患瘤脏器体积、肿瘤质心至前胸壁距离、肿瘤质心至后胸壁距离、患瘤脏器AP方向维度、胸腹部肿瘤SI方向维度、肿瘤质心至肿瘤所在脏器边缘距离、肿瘤质心至隔膜距离、肿瘤是否靠胸壁、肿瘤位置。
4.如权利要求1～3任一项所述的预测方法，其特征在于，步骤(5)所述机器学习模型通过如下方法构建而成：
1)基本模型的构建：
(a)数据预处理：对不属于同一量纲的影像组学特征进行无量纲化，剔除冗余信息，将定性特征进行哑编码转化为定量特征，并进行归一化处理；
公式如下：
上式fnor表示归一化的特征,f表示未处理的特征,fmean表示每类特征的平均值，SD表示每类特征的标准差；
(b)特征筛选：通过特征选择方法筛选出与AP、SI、LR三方向的肿瘤内部动度相关系数最高的特征；
(c)基于步骤(b)筛选出的特征和患者的临床信息特征构建得到机器学习模型，并采用机器学习算法对机器学习模型进行训练。
5.如权利要求4所述的预测方法，其特征在于，步骤(b)所述特征选择方法是：单变量特征选择方法、方差选择法、卡方检验、递归特征消除法、基于惩罚项的特征选择法、基于树模型的特征选择法或主成分分析法。
6.如权利要求4所述的预测方法，其特征在于，步骤(b)所述特征选择方法是单变量特征选择方法。
7.如权利要求4所述的预测方法，其特征在于，步骤(c)所述机器学习算法为神经网络算法、支持向量机、随机森林、XGBoost或广义线性模型，优选为神经网络算法。
8.如权利要求4所述的预测方法，其特征在于，步骤(c)所述机器学习算法为神经网络算法。
9.如权利要求7所述的预测方法，其特征在于，所述神经网络算法包括分支网络和神经元，所述分支网络包括手工特征网络、影像组学特征网络、全特征网络；所述手工特征网络学习患者的临床信息特征，所述影像组学特征网络学习影像组学特征；所述全特征网络学习临床信息特征与影像组学特征之间的知识；所述神经元学习各分支网络的权重。
10.一种患者体表和体内肿瘤运动相关性的预测方法，其特征在于，包括如下步骤：
1)体表数据采集：使用Varian实时位置管理系统通过红外摄像机跟踪安装在患者体表的红外反射标记，获取患者自由呼吸状态下AP方向的体表变化，得到10个呼吸时相对应的AP方向体表运动值；
2)肿瘤体内动度分析：采用权利要求1～7任一项所述的预测方法预测肿瘤内部动度，得到10个呼吸时相对应的肿瘤在体内AP、SI、LR三个方向的运动值；
3)将步骤1)得到的AP方向体表运动值与步骤2)得到的肿瘤在体内AP、SI、LR三个方向的最大运动值通过如下公式计算，得到预测的患者体表和体内肿瘤AP、SI、LR三个方向运动值的spearman相关性系数：
xi表示肿瘤在第i个呼吸时相对应的AP方向的体表运动值，表示AP方向所有呼吸时相对应的体表运动值的平均值；
yi表示第i个呼吸时相对应的肿瘤在体内AP、SI或LR方向的运动值，表示所有呼吸时相对应的肿瘤在体内AP、SI或LR方向的运动值的平均值；i为呼吸时相的序号，i为1到10的整数。
11.一种患者体表和体内肿瘤运动相关性及肿瘤内部动度的预测系统，其特征在于，包括如下模块：
数据采集和输入模块：对患者肿瘤部位进行CT扫描，扫描方式为电影扫描，得到10个扫描周期的扫描数据，每个扫描周期包含一个完整的呼吸周期，输入患者的临床信息特征；
数据处理模块：将数据采集模块获得的扫描数据进行4DCT图像重建，得到与所述的10个扫描周期的扫描数据对应的10个呼吸时相的CT，并对数据采集模块获得的扫描数据进行平均密度投影，得到AIP CT；
ROI勾画模块：对数据处理模块获得的10个呼吸时相的CT勾画ROI信息，结合AIP CT得到包含运动信息的ROI；
特征提取模块：从ROI勾画模块获得的ROI中提取影像组学特征，所述影像组学特征包括描述肿瘤3D物理外观的形态学特征、描述ROI强度的一阶分布的强度特征、原始纹理特征、高阶特征；
预测模块：将数据采集和输入模块获得的患者的临床信息特征，以及特征提取模块获得的影像组学特征输入机器学习模型，机器学习模型输出预测值；所述预测值包括：肿瘤在体内AP、SI、LR三个方向的最大运动值。
12.一种介质，其特征在于，所述介质存储有计算机可执行指令，所述指令在被处理器执行时可用于实现权利要求1～9任一项所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110956747.4A CN113499091B (zh) | 2021-08-19 | 2021-08-19 | 一种患者体表和体内肿瘤运动相关性及肿瘤内部动度的预测方法和系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110956747.4A CN113499091B (zh) | 2021-08-19 | 2021-08-19 | 一种患者体表和体内肿瘤运动相关性及肿瘤内部动度的预测方法和系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113499091A CN113499091A (zh) | 2021-10-15 |
| CN113499091B true CN113499091B (zh) | 2023-08-15 |
Family
ID=78015863
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110956747.4A Active CN113499091B (zh) | 2021-08-19 | 2021-08-19 | 一种患者体表和体内肿瘤运动相关性及肿瘤内部动度的预测方法和系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113499091B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114028738B (zh) * | 2021-11-16 | 2023-11-24 | 上海理工大学 | 一种三维有限元网格精准定位放射治疗系统 |
| CN117476110B (zh) * | 2023-10-10 | 2024-04-02 | 中山大学肿瘤防治中心（中山大学附属肿瘤医院、中山大学肿瘤研究所） | 基于人工智能的多尺度生物标记物发现系统 |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101623198A (zh) * | 2008-07-08 | 2010-01-13 | 深圳市海博科技有限公司 | 动态肿瘤实时跟踪方法 |
| WO2011133606A2 (en) * | 2010-04-19 | 2011-10-27 | The Regents Of The University Of California | Real-time volumetric image reconstruction and 3d tumor localization based on a single x-ray projection image for lung cancer radiotherapy |
| CN103530890A (zh) * | 2012-07-02 | 2014-01-22 | 株式会社东芝 | 运动校正装置以及方法 |
| WO2015193776A1 (en) * | 2014-06-18 | 2015-12-23 | Elekta Ab | System and method for automatic treatment planning |
| CN105828875A (zh) * | 2013-10-16 | 2016-08-03 | 皇家飞利浦有限公司 | 利用生物机械模型估计器官的位置 |
| CN106777976A (zh) * | 2016-12-15 | 2017-05-31 | 苏州大学 | 基于粒子滤波的放疗机器人肿瘤运动估计预测系统及方法 |
| CN109727672A (zh) * | 2018-12-28 | 2019-05-07 | 江苏瑞尔医疗科技有限公司 | 患者胸腹部肿瘤呼吸运动预测跟踪方法 |
| WO2019219387A1 (en) * | 2018-05-16 | 2019-11-21 | Koninklijke Philips N.V. | Automated tumor identification during surgery using machine-learning |
| CN110604550A (zh) * | 2019-09-24 | 2019-12-24 | 广州医科大学附属肿瘤医院 | 一种肿瘤放疗后正常组织器官并发症的预测方法 |
| CN110880195A (zh) * | 2019-10-23 | 2020-03-13 | 李夏东 | 基于肺部肿瘤4dct图像进行影像组学特征提取的专用图像的重建方法、介质和装置 |
| CN112581475A (zh) * | 2021-02-25 | 2021-03-30 | 四川大学华西医院 | 一种预测放疗计划的伽马通过率的方法及其应用 |
| CN113041515A (zh) * | 2021-03-25 | 2021-06-29 | 中国科学院近代物理研究所 | 三维图像引导运动器官定位方法、系统及存储介质 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140270448A1 (en) * | 2013-03-15 | 2014-09-18 | University Of Macau | System and method for attenuation correction in emission computed tomography |
| US10157462B2 (en) * | 2016-06-27 | 2018-12-18 | University Of Central Florida Research Foundation, Inc. | System and method for image-based quantification of white and brown adipose tissue at the whole-body, organ and body-region levels |
-
2021
- 2021-08-19 CN CN202110956747.4A patent/CN113499091B/zh active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101623198A (zh) * | 2008-07-08 | 2010-01-13 | 深圳市海博科技有限公司 | 动态肿瘤实时跟踪方法 |
| WO2011133606A2 (en) * | 2010-04-19 | 2011-10-27 | The Regents Of The University Of California | Real-time volumetric image reconstruction and 3d tumor localization based on a single x-ray projection image for lung cancer radiotherapy |
| CN103530890A (zh) * | 2012-07-02 | 2014-01-22 | 株式会社东芝 | 运动校正装置以及方法 |
| CN105828875A (zh) * | 2013-10-16 | 2016-08-03 | 皇家飞利浦有限公司 | 利用生物机械模型估计器官的位置 |
| WO2015193776A1 (en) * | 2014-06-18 | 2015-12-23 | Elekta Ab | System and method for automatic treatment planning |
| CN106777976A (zh) * | 2016-12-15 | 2017-05-31 | 苏州大学 | 基于粒子滤波的放疗机器人肿瘤运动估计预测系统及方法 |
| WO2019219387A1 (en) * | 2018-05-16 | 2019-11-21 | Koninklijke Philips N.V. | Automated tumor identification during surgery using machine-learning |
| CN109727672A (zh) * | 2018-12-28 | 2019-05-07 | 江苏瑞尔医疗科技有限公司 | 患者胸腹部肿瘤呼吸运动预测跟踪方法 |
| CN110604550A (zh) * | 2019-09-24 | 2019-12-24 | 广州医科大学附属肿瘤医院 | 一种肿瘤放疗后正常组织器官并发症的预测方法 |
| CN110880195A (zh) * | 2019-10-23 | 2020-03-13 | 李夏东 | 基于肺部肿瘤4dct图像进行影像组学特征提取的专用图像的重建方法、介质和装置 |
| CN112581475A (zh) * | 2021-02-25 | 2021-03-30 | 四川大学华西医院 | 一种预测放疗计划的伽马通过率的方法及其应用 |
| CN113041515A (zh) * | 2021-03-25 | 2021-06-29 | 中国科学院近代物理研究所 | 三维图像引导运动器官定位方法、系统及存储介质 |
Non-Patent Citations (1)
| Title |
|---|
| 黄志业.肿瘤放疗中基于GA-BP神经网络的呼吸运动预测技术.《中国优秀硕士论文全文数据库医药卫生科技辑》.2013,(第4期),I140-100. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113499091A (zh) | 2021-10-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11944463B2 (en) | Pseudo-CT generation from MR data using a feature regression model | |
| US10664723B2 (en) | Pseudo-CT generation from MR data using tissue parameter estimation | |
| CN108815721B (zh) | 一种照射剂量确定方法及系统 | |
| CN104408734B (zh) | 联合图像分割与形变配准技术的自适应目标区域转换方法 | |
| CN113499091B (zh) | 一种患者体表和体内肿瘤运动相关性及肿瘤内部动度的预测方法和系统 | |
| Hostettler et al. | A real-time predictive simulation of abdominal viscera positions during quiet free breathing | |
| CN107530552A (zh) | 用于自适应放射治疗的运动靶的三维定位 | |
| Rios et al. | Population model of bladder motion and deformation based on dominant eigenmodes and mixed-effects models in prostate cancer radiotherapy | |
| EP3468668A1 (en) | Soft tissue tracking using physiologic volume rendering | |
| Jafari et al. | In-vivo lung biomechanical modeling for effective tumor motion tracking in external beam radiation therapy | |
| Nasief et al. | Predicting necessity of daily online adaptive replanning based on wavelet image features for MRI guided adaptive radiation therapy | |
| Ono et al. | Estimation of lung tumor position from multiple anatomical features on 4D‐CT using multiple regression analysis | |
| Li et al. | Machine learning for predicting accuracy of lung and liver tumor motion tracking using radiomic features | |
| CN112997216B (zh) | 一种定位图像的转化系统 | |
| Ertan et al. | Assessment of the Reproducibility of Deep Inspiration Breath Hold Technique During Left-Sided Breast Cancer Radiotherapy | |
| Grova et al. | Validation of MRI/SPECT similarity-based registration methods using realistic simulations of normal and pathological SPECT data | |
| Dick | Fiducial-Less Real-Time Tracking of the Radiation Therapy of Liver Tumors Using Artificial Neural Networks | |
| Zhang | Application of online image guidance for moving tumour treatment using scanned proton therapy | |
| Geneser et al. | Incorporating patient breathing variability into a stochastic model of dose deposition for stereotactic body radiation therapy | |
| Ceberg et al. | Correlation between Surface and Tumour Motion in Lung Cancer-including Deep Learning Perspectives Caisa Kjellström | |
| Samavati | Hybrid Deformable Image Registration for Advanced Applications in Image-guided Interventions | |
| Banik et al. | Introduction to Medical Image Analysis | |
| Abeygunasekera | Framework for Respiratory-induced Motion Management and Modeling for the Reduction of the PTV Margin in EBRT | |
| Ramadaan | Validation of Deformable Image Registration for Head & Neck Cancer Adaptive Radiotherapy |
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 |