CN112927799A - 融合多示例学习和多任务深度影像组学的生存期分析系统 - Google Patents
融合多示例学习和多任务深度影像组学的生存期分析系统 Download PDFInfo
- Publication number
- CN112927799A CN112927799A CN202110393908.3A CN202110393908A CN112927799A CN 112927799 A CN112927799 A CN 112927799A CN 202110393908 A CN202110393908 A CN 202110393908A CN 112927799 A CN112927799 A CN 112927799A
- Authority
- CN
- China
- Prior art keywords
- clinical
- risk
- life cycle
- mri image
- image
- 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.)
- Granted
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 78
- 238000003384 imaging method Methods 0.000 title claims abstract description 29
- 238000007781 pre-processing Methods 0.000 claims abstract description 21
- 230000004083 survival effect Effects 0.000 claims description 56
- 238000000034 method Methods 0.000 claims description 50
- 206010028980 Neoplasm Diseases 0.000 claims description 17
- 238000004422 calculation algorithm Methods 0.000 claims description 13
- 230000003902 lesion Effects 0.000 claims description 13
- 206010027476 Metastases Diseases 0.000 claims description 12
- 238000013527 convolutional neural network Methods 0.000 claims description 12
- 230000009401 metastasis Effects 0.000 claims description 12
- 230000036961 partial effect Effects 0.000 claims description 12
- 238000012549 training Methods 0.000 claims description 12
- 238000012937 correction Methods 0.000 claims description 10
- 238000004393 prognosis Methods 0.000 claims description 10
- 238000009534 blood test Methods 0.000 claims description 8
- 230000035622 drinking Effects 0.000 claims description 8
- 230000000391 smoking effect Effects 0.000 claims description 8
- 239000000126 substance Substances 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 6
- 239000003550 marker Substances 0.000 claims description 4
- 238000010606 normalization Methods 0.000 claims description 4
- 238000000137 annealing Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 230000002829 reductive effect Effects 0.000 claims description 3
- 239000004576 sand Substances 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 2
- 238000000605 extraction Methods 0.000 abstract description 3
- 238000002595 magnetic resonance imaging Methods 0.000 description 53
- 230000008569 process Effects 0.000 description 11
- 238000010586 diagram Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 238000004891 communication Methods 0.000 description 4
- 208000002454 Nasopharyngeal Carcinoma Diseases 0.000 description 3
- 206010061306 Nasopharyngeal cancer Diseases 0.000 description 3
- 201000011216 nasopharynx carcinoma Diseases 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000004927 fusion Effects 0.000 description 2
- 230000000670 limiting effect Effects 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 210000001165 lymph node Anatomy 0.000 description 1
- 230000000873 masking effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000007637 random forest analysis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 210000001154 skull base Anatomy 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/20—ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/50—ICT 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Public Health (AREA)
- Medical Informatics (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Epidemiology (AREA)
- Data Mining & Analysis (AREA)
- Pathology (AREA)
- Databases & Information Systems (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computational Linguistics (AREA)
- Evolutionary Computation (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Molecular Biology (AREA)
- General Physics & Mathematics (AREA)
- Biophysics (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Image Analysis (AREA)
Abstract
本发明属于生物医学工程领域，具体涉及一种融合多示例学习和多任务深度影像组学的生存期分析系统，旨在解决传统的手工影像特征要求高精准、高重现性的病灶勾画，只适合提取单张影像图像的病灶特征，且建模时并未考虑多个临床结果的相关性，造成提取的影像特征鲁棒性较差，进而导致生存期分析结果准确性较低的问题。本发明系统包括：获取模块，配置为获取待分析的MR影像及对应的临床风险因子；预处理模块，配置为勾画病灶掩膜，并选取以病灶掩膜为中心的固定尺寸区域作为感兴趣区域；临床终点预测模块，配置为得到各个临床终点的风险预测值；分析模块，配置为得到生存期分析结果。本发明提高了影像特征提取的鲁棒性以及生存期分析结果的准确性。
Description
技术领域
本发明属于生物医学工程领域，具体涉及一种融合多示例学习和多任务深度影像组学的生存期分析系统、方法、设备。
背景技术
磁共振影像(Magnetic Resonance Imaging，MRI)和CT(Computed Tomograph)图像作为一种常规的诊疗工具和筛查工具，在临床上发挥着越来越重要的作用。比如，MRI图像作为鼻咽癌诊疗中一种常规、无创的影像分期工具，在软组织造影方面具有优势。研究表明，影像组学对鼻咽癌的生存期预测具有潜在价值。但由于有些影像图片的体素在不同维度之间的尺寸存在大的差异(比如鼻咽癌的MRI图片的像素横截面上尺寸为0.5mm，但层厚为5mm)，目前常见的方法是提取单张影像图像的病灶特征进行分析，例如，“一种基于影像组学的生存期预测方法及装置”，该专利采用的是手工的影像组学特征来预测患者的生存期，要求精准且重现性高的病灶勾画，适合提取单张影像图像的病灶特征，但没有考虑多个临床结果的相关性，导致了影像组学特征只能片面地表征肿瘤。
综上，由于传统的手工影像特征要求高精准、高重现性的病灶勾画，只适合提取单张影像图像的病灶特征，且建模时并未考虑多个临床结果的相关性，造成提取的影像特征鲁棒性较差，进而导致生存期分析结果准确性较低。基于此，本发明提出了一种融合多示例学习和多任务深度影像组学的生存期分析系统。
发明内容
为了解决现有技术中的上述问题，即为了解决传统的手工影像特征要求高精准、高重现性的病灶勾画，提取单张影像图像的病灶特征的片面性，建模时并未考虑多个临床结果的相关性，造成提取的影像特征鲁棒性较差、代表性弱，进而导致生存期分析结果准确性较低的问题，本发明提出了一种融合多示例学习和多任务深度影像组学的生存期分析系统，该系统包括：获取模块、预处理模块、临床终点预测模块、分析模块；
所述获取模块，配置为获取待分析的MR影像、对应的临床风险因子；所述临床风险因子包括患者的性别、年龄、TNM肿瘤分期、肿瘤家族史、吸烟史、饮酒史及血检指标；
所述预处理模块，配置为在MR影像的横断位MRI图像中勾画病灶掩膜；勾画后，对横断位MRI图像进行预处理，并在预处理后的横断位MRI图像中选取以勾画的病灶掩膜为中心的固定尺寸区域作为感兴趣区域；
所述临床终点预测模块，配置为将感兴趣区域以及相应的病灶掩膜输入预构建的生存数据预测模型，得到各个临床终点的风险预测值；所述临床终点包括死亡、转移、复发以及临床进展；
所述分析模块，配置为将各个临床终点的风险预测值和临床风险因子输入预构建的生存期分析模型，得到生存期分析结果；
其中，
所述生存数据预测模型是基于卷积神经网络和三个并行的全连接层构建；
所述生存期分析模型由基于LASSO的Cox比例风险模型构建。
在一些优选的实施方式中，所述预处理模块中的“对横断位MRI图像进行预处理”，其方法为：
通过双线性插值算法对横断位MRI图像进行像素尺寸调整；
利用直方图匹配方法对调整后的横断位MRI图像进行灰度校正；
通过Z-score方法对灰度校正后的横断位MRI图像进行灰度归一化处理。
在一些优选的实施方式中，“利用直方图匹配方法对调整后的横断位MRI图像进行灰度校正”，其方法为：
其中，x表示原横断位MRI图像的灰度值，xnew表示校正后的横断位MRI图像的灰度值，μ1i，μi和μ2i分别表示原横断位MRI图像灰度范围内的三个标志点，μ1s，μs和μ2s为目标域中的三个标志点。
在一些优选的实施方式中，所述临床终点预测模块“将感兴趣区域以及相应的病灶掩膜输入预构建的生存数据预测模型，得到各个临床终点的风险预测值”，其方法为：
将感兴趣区域以及相应的病灶掩膜输入生存数据预测模型的卷积神经网络，提取与生存期相关的高维影像特征；
基于提取的高维影像特征，通过生存数据预测模型的三个并行的全连接层，得到复发、转移、死亡三个临床终点的风险预测值；并将复发、转移、死亡三个临床终点的风险预测值的最大值作为临床进展对应的临床终点的风险预测值。
在一些优选的实施方式中，所述分析模块中的“将各个临床终点的风险预测值和临床风险因子输入预构建的生存期分析模型，得到生存期分析结果”，其方法为：
通过LASSO回归模型从各个临床终点的风险预测值和临床风险因子中选取独立的预后变量，作为风险变量；
将各风险变量输入Cox比例风险模型获取生存期分析结果；
其中，所述独立的预后变量的选取方法为：使用偏似然值作为LASSO回归模型的优化指标，当偏似然值达到最大值时停止筛选，余下的变量作为独立的预后变量。
在一些优选的实施方式中，所述生存数据预测模型，其训练方法为：
A10，采集多位患者治疗前的MR影像及其对应的临床风险因子、临床终点的随访时间和状态，构建训练数据集；
A20，在MR影像的横断位MRI图像中勾画病灶掩膜；勾画后，对横断位MRI图像进行预处理，并在预处理后的横断位MRI图像中选取以勾画的病灶掩膜为中心的固定尺寸区域作为感兴趣区域；
A30，分批的将感兴趣区域以及相应的病灶掩膜输入预构建的生存数据预测模型，得到各个临床终点的风险预测值；
A40，基于各个临床终点的风险预测值、所述临床终点的随访时间和状态，计算负的Cox偏似然损失值，并使用误差反传算法更新生存数据预测模型的模型参数；
A50，循环步骤A30-A40，直至负的Cox偏似然损失值不再变小，得到训练好的生存数据预测模型。
在一些优选的实施方式中，所述生存数据预测模型其在训练时，学习率策略为余弦退火下降策略，误差反传算法为随机梯度下降算法，batch size设置为128。
本发明的第二方面，提出了一种融合多示例学习和多任务深度影像组学的生存期分析方法，该方法包括：
S10，获取待分析的MR影像、对应的临床风险因子；所述临床风险因子包括患者的性别、年龄、TNM肿瘤分期、肿瘤家族史、吸烟史、饮酒史及血检指标；
S20，在MR影像的横断位MRI图像中勾画病灶掩膜；勾画后，对横断位MRI图像进行预处理，并在预处理后的横断位MRI图像中选取以勾画的病灶掩膜为中心的固定尺寸区域作为感兴趣区域；
S30，将感兴趣区域以及相应的病灶掩膜输入预构建的生存数据预测模型，得到各个临床终点的风险预测值；所述临床终点包括死亡、转移、复发以及临床进展；
S40，将各个临床终点的风险预测值和临床风险因子输入预构建的生存期分析模型，得到生存期分析结果；
其中，
所述生存数据预测模型是基于卷积神经网络和三个并行的全连接层构建；
所述生存期分析模型由基于LASSO的Cox比例风险模型构建。
本发明的第三方面，提出了一种包括：至少一个处理器；以及与至少一个所述处理器通信连接的存储器；其中，所述存储器存储有可被所述处理器执行的指令，所述指令用于被所述处理器执行以实现上述的融合多示例学习和多任务深度影像组学的生存期分析方法。
本发明的第四方面，提出了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于被所述计算机执行以实现上述的融合多示例学习和多任务深度影像组学的生存期分析方法。
本发明的有益效果：
本发明提高了影像特征提取的鲁棒性以及生存期分析结果的准确性。
1)本发明提出的生存数据预测模型以病灶为中心的图像块为输入，不要求精准的病灶勾画，同时输入病灶掩模，告知模型应重点关注的区域，且模型能自动地、端到端地输出预测不同临床终点的生存数据的预测标签。这样在避免了影像体素尺寸各向异性的同时，充分利用病灶的三维信息，提高了影像特征提取的鲁棒性。
2)本发明利用了多示例学习来挖掘与多个临床终点的生存数据的预测标签，结合独立的临床风险因子可同时建立针对不同临床终点的生存期预测模型，可无创地得到患者的针对多个临床终点的生存期评估，提升了生存期分析结果的准确性。
附图说明
通过阅读参照以下附图所做的对非限制性实施例所做的详细描述，本申请的其他特征、目的和优点将会变得更明显。
图1是本发明一种实施例的融合多示例学习和多任务深度影像组学的生存期分析系统的框架示意图；
图2为本发明一种实施例的生存数据预测模型的训练流程示意图；
图3是本发明一种实施例的融合多示例学习和多任务深度影像组学的生存期分析方法的流程示意图；
图4是本发明一种实施例的生存数据预测模型的结构示意图；
图5是本发明一种实施例的适于用来实现本申请实施例的电子设备的计算机系统的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
下面结合附图和实施例对本申请作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释相关发明，而非对该发明的限定。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。
本发明第一实施例的一种融合多示例学习和多任务深度影像组学的生存期分析系统，如图1所示，该系统包括：获取模块、预处理模块、临床终点预测模块、分析模块；
所述获取模块，配置为获取待分析的MR影像、对应的临床风险因子；所述临床风险因子包括患者的性别、年龄、TNM肿瘤分期、肿瘤家族史、吸烟史、饮酒史及血检指标；
所述预处理模块，配置为在MR影像的横断位MRI图像中勾画病灶掩膜；勾画后，对横断位MRI图像进行预处理，并在预处理后的横断位MRI图像中选取以勾画的病灶掩膜为中心的固定尺寸区域作为感兴趣区域；
所述临床终点预测模块，配置为将感兴趣区域以及相应的病灶掩膜输入预构建的生存数据预测模型，得到各个临床终点的风险预测值；所述临床终点包括死亡、转移、复发以及临床进展；
所述分析模块，配置为将各个临床终点的风险预测值和临床风险因子输入预构建的生存期分析模型，得到生存期分析结果；
其中，
所述生存数据预测模型是基于卷积神经网络和三个并行的全连接层构建；
所述生存期分析模型由基于LASSO的Cox比例风险模型构建。
为了更清晰地对本发明融合多示例学习和多任务深度影像组学的生存期分析系统进行说明，下面对本发明系统一种实施例中各模块进行展开详述。
在下述实施例中，先对生存数据预测模型的训练过程进行详述，再对基于融合多示例学习和多任务深度影像组学的生存期分析系统获取生存期分析结果的过程进行详述。
1、生存数据预测模型的训练过程，如图2所示
A10，采集多位患者治疗前的MR影像及其对应的临床风险因子、临床终点的随访时间和状态，构建训练数据集；
在本实施例中，先采集多位患者(例如鼻咽癌患者)治疗前的MR影像(本发明中以每名患者的横断位的MR影像为包)、临床风险因子(如性别、年龄、TNM肿瘤分期、肿瘤家族史、吸烟史、饮酒史及血检指标等)和多个临床终点(包括死亡、转移、复发、临床进展(即病情进展))的随访时间和状态，作为训练数据，构建训练数据集。
其中，各个临床终点的随访时间为从某个时间点开始到发生该临床事件或最后的随访时间点为止的时间长度，起始时间点由使用者定义，本发明优选以开始接受治疗时间点为起始时间点；临床终点的状态是一个二值变量，代表在随访时间内是否发生该临床事件。
A20，在MR影像的横断位MRI图像中勾画病灶掩膜；勾画后，对横断位MRI图像进行预处理，并在预处理后的横断位MRI图像中选取以勾画的病灶掩膜为中心的固定尺寸区域作为感兴趣区域；
在本实施例中，让有经验的医生在横断位的MRI图像中病灶掩膜，比如鼻咽层面和颅底层面的原发灶，若有必要也可勾画咽后淋巴结进行分析。在勾画的过程中，若有冠状位和矢状位的MRI图像，可以用作勾画的参考。
勾画完病灶掩膜后，我们还需要对横断位的MRI图像进行预处理，包括像素尺寸调整、灰度校正和灰度归一化处理。具体如下：
使用双线性插值算法对MRI横截面图像的进行像素尺寸调整；即将像素尺寸换到相同的大小尺度；
使用直方图匹配方法对调整后的横断位MRI图像进行灰度校正，将各患者的影像变换到同一目标域，以削弱机型和机型参数对影像灰度的影响。校正后的横断位MRI图像灰度值和原来的横断位MRI图像灰度值的关系为：xnew＝f(x)。
在本发明中优选为Sun等人提出的直方图匹配方法，如公式(1)所示：
其中，x表示原横断位MRI图像的灰度值，xnew表示校正后的横断位MRI图像的灰度值，μ1i，μi和μ2i分别表示原横断位MRI图像灰度范围内的三个标志点，μ1s，μs和μ2s为目标域中的三个标志点。
为了加快深度学习网络的收敛，我们使用Z-score方法来归一化输入的感兴趣区域的灰度值，Z-score的定义为：
f(x)＝(x-x_mean)/x_std (2)
其中，x_mean为用于训练的感兴趣区域的平均灰度值，x_std为用于训练的感兴趣区域的灰度值的方差。
预处理结束后，取以勾画的病灶掩膜为中心的固定尺寸的图像方块为感兴趣区域。
A30，分批的将感兴趣区域以及相应的病灶掩膜输入预构建的生存数据预测模型，得到各个临床终点的风险预测值；
在本实施例中，训练生存数据预测模型(基于卷积神经网络、三个并行的全连接层构建，卷积神经网络可以使用但不局限于VGG，ResNet和DenseNet框架，网络的深度由训练集的数据规模决定，一般不超过50层)时，输入样本有两个通道，一个通道为感兴趣区域，另一个通道为感兴趣区域的病灶掩模，如图4所示，输入的标签为四个临床终点的随访时间和状态。
卷积神经网络的终端使用三个并列的全连接层分别同时预测死亡、转移和复发的风险值，然后使用三个风险值中的最大值作为临床进展的风险值；在本发明中，每名患者的所有的横断位图像及其病灶掩模为包，每张横断位图像及其病灶掩模为一个示例(即上文的一个输入样本)，根据多实例学习假设，如果一个包的中有一个以上的示例被标定为正的，那么这个包就被标定为正的，因此本发明中若某名患者有一个输入样本被判断为高风险值，则该患者属于高危病例。
生存数据预测模型的主干网络，即卷积神经网络在本发明中优选用ResNet-18，使用验证集方法来选择最优的网络参数，得到各个临床终点的风险预测值，作为预测标签。
A40，基于各个临床终点的风险预测值、所述临床终点的随访时间和状态，计算负的Cox偏似然损失值，并使用误差反传算法更新生存数据预测模型的模型参数；
在本实施例中，基于临床进展的风险预测值，我们使用多示例学习选择每位患者的风险预测值最大的示例，然后用这些示例来更新网络权重。生存数据预测模型的batchsize可以设置为128，学习率策略为余弦退火下降策略，误差反传算法为随机梯度下降算法。
A50，循环步骤A30-A40，直至负的Cox偏似然损失值不再变小，得到训练好的生存数据预测模型。
在本实施例中，将感兴趣区域以及相应的病灶掩膜分批输入生存数据预测模型，对生存数据预测模型进行训练，直至得到训练好的生存数据预测模型。
2、融合多示例学习和多任务深度影像组学的生存期分析系统
所述获取模块，配置为获取待分析的MR影像、对应的临床风险因子；所述临床风险因子包括患者的性别、年龄、TNM肿瘤分期、肿瘤家族史、吸烟史、饮酒史及血检指标；
在本实施例中，获取待分析的MR影像及其对应的临床风险因子。
所述预处理模块，配置为在MR影像的横断位MRI图像中勾画病灶掩膜；勾画后，对横断位MRI图像进行预处理，并在预处理后的横断位MRI图像中选取以勾画的病灶掩膜为中心的固定尺寸区域作为感兴趣区域；
在本实施例中，对横断位MRI图像进行预处理，具体为：
通过双线性插值算法对横断位MRI图像进行像素尺寸调整；
利用直方图匹配方法对调整后的横断位MRI图像进行灰度校正；
通过Z-score方法对灰度校正后的横断位MRI图像进行灰度归一化处理。
所述临床终点预测模块，配置为将感兴趣区域以及相应的病灶掩膜输入预构建的生存数据预测模型，得到各个临床终点的风险预测值；所述临床终点包括死亡、转移、复发以及临床进展；
在本实施例中，基于感兴趣区域以及相应的病灶掩膜，通过上述训练好的生存数据预测模型，得到各个临床终点的风险预测值。即使用三个并列的全连接层同时预测死亡、转移和复发的风险，然后使用三个风险中的最大值作为临床进展的风险预测值。
所述分析模块，配置为将各个临床终点的风险预测值和临床风险因子输入预构建的生存期分析模型，得到生存期分析结果。
在本实施例中，为融合不同层次(如影像，TNM分期，年龄和治疗前的EBV-DNA水平等)的预后风险信息，我们使用LASSO(Least absolute shrinkage and selectionoperator)回归模型，也可以替换为随机森林等特征选择方法，从各个临床终点的风险预测值和临床风险因子中选择独立的预后变量(即从各个临床终点的风险预测值和临床风险因子包含的变量中选择出一部分变量，作为预后变量)，然后使用Cox比例风险模型建立完善的生存期预测模型，即生存期预测模型基于LASSO回归模型、Cox比例风险模型构建。使用预测不同临床终点的风险预测值，再结合独立的临床风险因子，通过构建的生存期分析模型，可以预测相应的生存期分析结果。
另外，本发明中，预后变量的选取方法为：以偏似然值作为LASSO回归模型的优化指标，偏似然值达到最大值时保留的特征变量为独立的预后变量。将选取后的独立的预后变量送入Cox比例风险模型确定各变量的权重系数以得到生存期分析模型。
需要说明的是，上述实施例提供的融合多示例学习和多任务深度影像组学的生存期分析系统，仅以上述各功能模块的划分进行举例说明，在实际应用中，可以根据需要而将上述功能分配由不同的功能模块来完成，即将本发明实施例中的模块或者步骤再分解或者组合，例如，上述实施例的模块可以合并为一个模块，也可以进一步拆分成多个子模块，以完成以上描述的全部或者部分功能。对于本发明实施例中涉及的模块、步骤的名称，仅仅是为了区分各个模块或者步骤，不视为对本发明的不当限定。
本发明第二实施例的一种融合多示例学习和多任务深度影像组学的生存期分析方法，如图3所示，包括：
S10，获取待分析的MR影像、对应的临床风险因子；所述临床风险因子包括患者的性别、年龄、TNM肿瘤分期、肿瘤家族史、吸烟史、饮酒史及血检指标；
S20，在MR影像的横断位MRI图像中勾画病灶区域以及病灶掩膜；勾画后，对横断位MRI图像进行预处理，并在预处理后的横断位MRI图像中选取以勾画的病灶区域为中心的固定尺寸区域作为感兴趣区域；
S30，将感兴趣区域以及感兴趣区域的病灶掩膜输入预构建的生存数据预测模型，得到各个临床终点的风险预测值；所述临床终点包括死亡、转移、复发以及临床进展；
S40，将各个临床终点的风险预测值和临床风险因子输入预构建的生存期分析模型，得到生存期分析结果；
其中，
所述生存数据预测模型是基于卷积神经网络和三个并行的全连接层构建；
所述生存期分析模型由基于LASSO的Cox比例风险模型构建。
所述技术领域的技术人员可以清楚的了解到，为描述的方便和简洁，上述描述的方法的具体的工作过程及有关说明，可以参考前述系统实施例中的对应过程，在此不再赘述。
本发明第三实施例的一种设备，包括：至少一个处理器；以及与至少一个所述处理器通信连接的存储器；其中，所述存储器存储有可被所述处理器执行的指令，所述指令用于被所述处理器执行以实现上述的融合多示例学习和多任务深度影像组学的生存期分析方法。
本发明第四实施例的一种计算机可读存储介质，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于被所述计算机执行以实现上述的融合多示例学习和多任务深度影像组学的生存期分析方法。
所述技术领域的技术人员可以清楚的了解到，为描述的方便和简洁，上述描述的存储装置、处理装置的具体工作过程及有关说明，可以参考前述方法实例中的对应过程，在此不再赘述。
下面参考图5，其示出了适于用来实现本申请方法、系统、设备实施例的服务器的计算机系统的结构示意图。图5示出的服务器仅仅是一个示例，不应对本申请实施例的功能和使用范围带来任何限制。
如图5所示，计算机系统包括中央处理单元(CPU，Central Processing Unit)501，其可以根据存储在只读存储器(ROM，Read Only Memory)502中的程序或者从存储部分508加载到随机访问存储器(RAM，Random Access Memory)503中的程序而执行各种适当的动作和处理。在RAM503中，还存储有系统操作所需的各种程序和数据。CPU501、ROM502以及RAM503通过总线504彼此相连。输入/输出(I/O，Input/Output)接口505也连接至总线504。
以下部件连接至I/O接口505：包括键盘、鼠标等的输入部分506；包括诸如阴极射线管(CRT，Cathode Ray Tube)、液晶显示器(LCD，Liquid Crystal Display)等以及扬声器等的输出部分507；包括硬盘等的存储部分508；以及包括诸如LAN(局域网，Local AreaNetwork)卡、调制解调器等的网络接口卡的通讯部分509。通讯部分509经由诸如因特网的网络执行通信处理。驱动器510也根据需要连接至I/O接口505。可拆卸介质511，诸如磁盘、光盘、磁光盘、半导体存储器等等，根据需要安装在驱动器510上，以便于从其上读出的计算机程序根据需要被安装入存储部分508。
特别地，根据本公开的实施例，上文参考流程图描述的过程可以被实现为计算机软件程序。例如，本公开的实施例包括一种计算机程序产品，其包括承载在计算机可读介质上的计算机程序，该计算机程序包含用于执行流程图所示的方法的程序代码。在这样的实施例中，该计算机程序可以通过通讯部分509从网络上被下载和安装，和/或从可拆卸介质511被安装。在该计算机程序被中央处理单元(CPU501执行时，执行本申请的方法中限定的上述功能。需要说明的是，本申请上述的计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质或者是上述两者的任意组合。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件，或者任意以上的组合。计算机可读存储介质的更具体的例子可以包括但不限于：具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机访问存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑磁盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本申请中，计算机可读存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。而在本申请中，计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号，其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式，包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质，该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。计算机可读介质上包含的程序代码可以用任何适当的介质传输，包括但不限于：无线、电线、光缆、RF等等，或者上述的任意合适的组合。
可以以一种或多种程序设计语言或其组合来编写用于执行本申请的操作的计算机程序代码，上述程序设计语言包括面向对象的程序设计语言-诸如Java、Smalltalk、C++，还包括常规的过程式程序设计语言-诸如“C”语言或类似的程序设计语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络——包括局域网(LAN)或广域网(WAN)-连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
附图中的流程图和框图，图示了按照本申请各种实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段、或代码的一部分，该模块、程序段、或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个接连地表示的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或操作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。
术语“第一”、“第二”等是用于区别类似的对象，而不是用于描述或表示特定的顺序或先后次序。
术语“包括”或者任何其它类似用语旨在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备/装置不仅包括那些要素，而且还包括没有明确列出的其它要素，或者还包括这些过程、方法、物品或者设备/装置所固有的要素。
至此，已经结合附图所示的优选实施方式描述了本发明的技术方案，但是，本领域技术人员容易理解的是，本发明的保护范围显然不局限于这些具体实施方式。在不偏离本发明的原理的前提下，本领域技术人员可以对相关技术特征作出等同的更改或替换，这些更改或替换之后的技术方案都将落入本发明的保护范围之内。
Claims (10)
1.一种融合多示例学习和多任务深度影像组学的生存期分析系统，其特征在于，该系统包括：获取模块、预处理模块、临床终点预测模块、分析模块；
所述获取模块，配置为获取待分析的MR影像、对应的临床风险因子；所述临床风险因子包括患者的性别、年龄、TNM肿瘤分期、肿瘤家族史、吸烟史、饮酒史及血检指标；
所述预处理模块，配置为在MR影像的横断位MRI图像中勾画病灶掩膜；勾画后，对横断位MRI图像进行预处理，并在预处理后的横断位MRI图像中选取以勾画的病灶掩膜为中心的固定尺寸区域作为感兴趣区域；
所述临床终点预测模块，配置为将感兴趣区域以及相应的病灶掩膜输入预构建的生存数据预测模型，得到各个临床终点的风险预测值；所述临床终点包括死亡、转移、复发以及临床进展；
所述分析模块，配置为将各个临床终点的风险预测值和临床风险因子输入预构建的生存期分析模型，得到生存期分析结果；
其中，
所述生存数据预测模型是基于卷积神经网络和三个并行的全连接层构建；
所述生存期分析模型由基于LASSO的Cox比例风险模型构建。
2.根据权利要求1所述的融合多示例学习和多任务深度影像组学的生存期分析系统，其特征在于，所述预处理模块中的“对横断位MRI图像进行预处理”，其方法为：
通过双线性插值算法对横断位MRI图像进行像素尺寸调整；
利用直方图匹配方法对调整后的横断位MRI图像进行灰度校正；
通过Z-score方法对灰度校正后的横断位MRI图像进行灰度归一化处理。
3.根据权利要求2所述的融合多示例学习和多任务深度影像组学的生存期分析系统，其特征在于，“利用直方图匹配方法对调整后的横断位MRI图像进行灰度校正”，其方法为：
其中，x表示原横断位MRI图像的灰度值，xnew表示校正后的横断位MRI图像的灰度值，μ1i，μi和μ2i分别表示原横断位MRI图像灰度范围内的三个标志点，μ1s，μs和μ2s为目标域中的三个标志点。
4.根据权利要求1所述的融合多示例学习和多任务深度影像组学的生存期分析系统，所述临床终点预测模块“将感兴趣区域以及感兴趣区域的病灶掩膜输入预构建的生存数据预测模型，得到各个临床终点的风险预测值”，其方法为：
将感兴趣区域以及相应的病灶掩膜输入生存数据预测模型的卷积神经网络，提取与生存期相关的高维影像特征；
基于提取的高维影像特征，通过生存数据预测模型的三个并行的全连接层，得到复发、转移、死亡三个临床终点的风险预测值；并将复发、转移、死亡三个临床终点的风险预测值的最大值作为临床进展对应的临床终点的风险预测值。
5.根据权利要求1所述的融合多示例学习和多任务深度影像组学的生存期分析系统，其特征在于，所述分析模块中的“将各个临床终点的风险预测值和临床风险因子输入预构建的生存期分析模型，得到生存期分析结果”，其方法为：
通过LASSO回归模型从各个临床终点的风险预测值和临床风险因子中选取独立的预后变量，作为风险变量；
将各风险变量输入Cox比例风险模型获取生存期分析结果；
其中，所述独立的预后变量的选取方法为：使用偏似然值作为LASSO回归模型的优化指标，当偏似然值达到最大值时停止筛选，余下的变量作为独立的预后变量。
6.根据权利要求1所述的融合多示例学习和多任务深度影像组学的生存期分析系统，其特征在于，所述生存数据预测模型，其训练方法为：
A10，采集多位患者治疗前的MR影像及其对应的临床风险因子、临床终点的随访时间和状态，构建训练数据集；
A20，在MR影像的横断位MRI图像中勾画病灶掩膜；勾画后，对横断位MRI图像进行预处理，并在预处理后的横断位MRI图像中选取以勾画的病灶掩膜为中心的固定尺寸区域作为感兴趣区域；
A30，分批的将感兴趣区域以及相应的病灶掩膜输入预构建的生存数据预测模型，得到各个临床终点的风险预测值；
A40，基于各个临床终点的风险预测值、所述临床终点的随访时间和状态，计算负的Cox偏似然损失值，并使用误差反传算法更新生存数据预测模型的模型参数；
A50，循环步骤A30-A40，直至负的Cox偏似然损失值不再变小，得到训练好的生存数据预测模型。
7.根据权利要求6所述的融合多示例学习和多任务深度影像组学的生存期分析系统，其特征在于，所述生存数据预测模型其在训练时，学习率策略为余弦退火下降策略，误差反传算法为随机梯度下降算法，batch size设置为128。
8.一种融合多示例学习和多任务深度影像组学的生存期分析方法，其特征在于，该方法包括：
S10，获取待分析的MR影像、对应的临床风险因子；所述临床风险因子包括患者的性别、年龄、TNM肿瘤分期、肿瘤家族史、吸烟史、饮酒史及血检指标；
S20，在MR影像的横断位MRI图像中勾画病灶掩膜；勾画后，对横断位MRI图像进行预处理，并在预处理后的横断位MRI图像中选取以勾画的病灶掩膜为中心的固定尺寸区域作为感兴趣区域；
S30，将感兴趣区域以及相应的病灶掩膜输入预构建的生存数据预测模型，得到各个临床终点的风险预测值；所述临床终点包括死亡、转移、复发以及临床进展；
S40，将各个临床终点的风险预测值输入预构建的生存期分析模型，得到生存期分析结果；
其中，
所述生存数据预测模型是基于卷积神经网络和三个并行的全连接层构建；
所述生存期分析模型基于LASSO回归模型、Cox比例风险模型构建。
9.一种设备，其特征在于，包括：
至少一个处理器；以及
与至少一个所述处理器通信连接的存储器；其中，
所述存储器存储有可被所述处理器执行的指令，所述指令用于被所述处理器执行以实现权利要求8所述的融合多示例学习和多任务深度影像组学的生存期分析方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于被所述计算机执行以实现权利要求8所述的融合多示例学习和多任务深度影像组学的生存期分析方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110393908.3A CN112927799B (zh) | 2021-04-13 | 2021-04-13 | 融合多示例学习和多任务深度影像组学的生存期分析系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110393908.3A CN112927799B (zh) | 2021-04-13 | 2021-04-13 | 融合多示例学习和多任务深度影像组学的生存期分析系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112927799A true CN112927799A (zh) | 2021-06-08 |
| CN112927799B CN112927799B (zh) | 2023-06-27 |
Family
ID=76174231
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110393908.3A Active CN112927799B (zh) | 2021-04-13 | 2021-04-13 | 融合多示例学习和多任务深度影像组学的生存期分析系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112927799B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113689382A (zh) * | 2021-07-26 | 2021-11-23 | 北京知见生命科技有限公司 | 基于医学影像和病理图像的肿瘤术后生存期预测方法及系统 |
| CN113723415A (zh) * | 2021-08-26 | 2021-11-30 | 泰康保险集团股份有限公司 | 一种生存时长预测方法、装置、设备及介质 |
| CN115188470A (zh) * | 2022-06-29 | 2022-10-14 | 山东大学 | 一种基于多任务Cox学习模型的多慢性病预测系统 |
| CN115439473A (zh) * | 2022-11-04 | 2022-12-06 | 北京精诊医疗科技有限公司 | 一种基于交互分组注意机制的多期相占位分类方法 |
| CN115830001A (zh) * | 2022-12-22 | 2023-03-21 | 抖音视界有限公司 | 肠道图像处理方法、装置、存储介质及电子设备 |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108257135A (zh) * | 2018-02-01 | 2018-07-06 | 浙江德尚韵兴图像科技有限公司 | 基于深度学习方法解读医学图像特征的辅助诊断系统 |
| CN109949288A (zh) * | 2019-03-15 | 2019-06-28 | 上海联影智能医疗科技有限公司 | 肿瘤类型确定系统、方法及存储介质 |
| EP3576100A1 (en) * | 2018-05-30 | 2019-12-04 | Siemens Healthcare GmbH | Decision support system for medical therapy planning |
| CN110705555A (zh) * | 2019-09-17 | 2020-01-17 | 中山大学 | 基于fcn的腹部多器官核磁共振图像分割方法、系统及介质 |
| CN110974288A (zh) * | 2019-12-26 | 2020-04-10 | 北京大学口腔医学院 | 一种牙周病cbct纵向数据记录及分析方法 |
| CN111128328A (zh) * | 2019-10-25 | 2020-05-08 | 中山大学肿瘤防治中心（中山大学附属肿瘤医院、中山大学肿瘤研究所） | 鼻咽癌结构化影像报告及数据处理系统和方法 |
| CN111178316A (zh) * | 2020-01-06 | 2020-05-19 | 武汉大学 | 一种基于深度架构自动搜索的高分辨率遥感影像土地覆盖分类方法 |
| US20210027890A1 (en) * | 2019-07-24 | 2021-01-28 | ConnSante BioTech, Inc. | Detecting, evaluating and predicting system for cancer risk |
| CN112309576A (zh) * | 2020-09-22 | 2021-02-02 | 江南大学 | 基于深度学习ct影像组学的结直肠癌生存期预测方法 |
-
2021
- 2021-04-13 CN CN202110393908.3A patent/CN112927799B/zh active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108257135A (zh) * | 2018-02-01 | 2018-07-06 | 浙江德尚韵兴图像科技有限公司 | 基于深度学习方法解读医学图像特征的辅助诊断系统 |
| EP3576100A1 (en) * | 2018-05-30 | 2019-12-04 | Siemens Healthcare GmbH | Decision support system for medical therapy planning |
| CN109949288A (zh) * | 2019-03-15 | 2019-06-28 | 上海联影智能医疗科技有限公司 | 肿瘤类型确定系统、方法及存储介质 |
| US20210027890A1 (en) * | 2019-07-24 | 2021-01-28 | ConnSante BioTech, Inc. | Detecting, evaluating and predicting system for cancer risk |
| CN110705555A (zh) * | 2019-09-17 | 2020-01-17 | 中山大学 | 基于fcn的腹部多器官核磁共振图像分割方法、系统及介质 |
| CN111128328A (zh) * | 2019-10-25 | 2020-05-08 | 中山大学肿瘤防治中心（中山大学附属肿瘤医院、中山大学肿瘤研究所） | 鼻咽癌结构化影像报告及数据处理系统和方法 |
| CN110974288A (zh) * | 2019-12-26 | 2020-04-10 | 北京大学口腔医学院 | 一种牙周病cbct纵向数据记录及分析方法 |
| CN111178316A (zh) * | 2020-01-06 | 2020-05-19 | 武汉大学 | 一种基于深度架构自动搜索的高分辨率遥感影像土地覆盖分类方法 |
| CN112309576A (zh) * | 2020-09-22 | 2021-02-02 | 江南大学 | 基于深度学习ct影像组学的结直肠癌生存期预测方法 |
Non-Patent Citations (4)
| Title |
|---|
| DONG NIE等: "3D deep learning for multi-modal imaging-guided survival time prediction of brain tumor patients", 《MEDICAL IMAGE COMPUTING AND COMPUTER-ASSISTED INTERVENTION – MICCAI 2016 》 * |
| KAOUTAR B. AHMED 等: "Fine-tuning convolutional deep features for MRI based brain tumor classification", 《MEDICAL IMAGING 2017: COMPUTER-AIDED DIAGNOSIS》 * |
| LIWEN ZHANG 等: "Multi-Focus Network to Decode Imaging Phenotype for Overall Survival Prediction of Gastric Cancer Patients", 《IEEE JOURNAL OF BIOMEDICAL AND HEALTH INFORMATICS》 * |
| 崔达华;赵莹;刘爱连;武敬君;郭妍;李昕;吴艇帆;崔景景;左盼莉;: "基于动态增强MRI的影像组学列线图预测肝细胞癌切除术后3年复发的价值", 中国临床医学影像杂志, no. 12 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113689382A (zh) * | 2021-07-26 | 2021-11-23 | 北京知见生命科技有限公司 | 基于医学影像和病理图像的肿瘤术后生存期预测方法及系统 |
| CN113689382B (zh) * | 2021-07-26 | 2023-12-01 | 北京知见生命科技有限公司 | 基于医学影像和病理图像的肿瘤术后生存期预测方法及系统 |
| CN113723415A (zh) * | 2021-08-26 | 2021-11-30 | 泰康保险集团股份有限公司 | 一种生存时长预测方法、装置、设备及介质 |
| CN113723415B (zh) * | 2021-08-26 | 2023-07-18 | 泰康保险集团股份有限公司 | 一种生存时长预测方法、装置、设备及介质 |
| CN115188470A (zh) * | 2022-06-29 | 2022-10-14 | 山东大学 | 一种基于多任务Cox学习模型的多慢性病预测系统 |
| CN115439473A (zh) * | 2022-11-04 | 2022-12-06 | 北京精诊医疗科技有限公司 | 一种基于交互分组注意机制的多期相占位分类方法 |
| CN115830001A (zh) * | 2022-12-22 | 2023-03-21 | 抖音视界有限公司 | 肠道图像处理方法、装置、存储介质及电子设备 |
| CN115830001B (zh) * | 2022-12-22 | 2023-09-08 | 抖音视界有限公司 | 肠道图像处理方法、装置、存储介质及电子设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112927799B (zh) | 2023-06-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112927799B (zh) | 融合多示例学习和多任务深度影像组学的生存期分析系统 | |
| Chen et al. | Source-free domain adaptive fundus image segmentation with denoised pseudo-labeling | |
| CN113571203B (zh) | 多中心基于联邦学习的脑肿瘤预后生存期预测方法及系统 | |
| CN106815481B (zh) | 一种基于影像组学的生存期预测方法及装置 | |
| US20180144244A1 (en) | Distributed clinical workflow training of deep learning neural networks | |
| BR112021015324A2 (pt) | Mascaramento de sombra e nuvem para aplicações agrícolas com o uso de redes neurais convolucionais | |
| CN107481254A (zh) | 医学图像的处理方法、装置、介质和电子设备 | |
| Dorent et al. | Inter extreme points geodesics for end-to-end weakly supervised image segmentation | |
| CN110490927A (zh) | 用于为图像中的对象生成中心线的方法、装置和系统 | |
| CN105608687A (zh) | 医用图像处理方法和医用图像处理装置 | |
| CN110717518A (zh) | 基于3d卷积神经网络的持续性肺结节识别方法及装置 | |
| CN112396606B (zh) | 一种基于用户交互的医学图像分割方法、系统和装置 | |
| TW202125345A (zh) | 用於組織分割之機器學習模型的交互式訓練 | |
| CN114944229A (zh) | 基于深度学习和磁共振结构脑影像的脑龄预测方法 | |
| CN111863247B (zh) | 基于结构磁共振影像的大脑年龄级联精炼预测方法及系统 | |
| CN111968130A (zh) | 脑部造影图像处理方法、装置、介质及电子设备 | |
| CN115330669A (zh) | 预测解剖结构的疾病量化参数的计算机实现的方法、系统及存储介质 | |
| CN116309571A (zh) | 一种基于半监督学习的三维脑血管分割方法和装置 | |
| CN114332132A (zh) | 图像分割方法、装置和计算机设备 | |
| CN110599444B (zh) | 预测血管树的血流储备分数的设备、系统以及非暂时性可读存储介质 | |
| CN109949300B (zh) | 解剖树结构分析的方法、系统和计算机可读介质 | |
| CN112801999A (zh) | 确定心脏冠脉优势型的方法及装置 | |
| CN114972220B (zh) | 一种图像处理方法、装置、电子设备及可读存储介质 | |
| CN108198213B (zh) | 一种标记管状结构的方法及装置、管状结构的处理终端 | |
| CN115861255A (zh) | 用于图像处理的模型训练方法、装置、设备、介质及产品 |
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 |