CN106815481B - 一种基于影像组学的生存期预测方法及装置 - Google Patents
一种基于影像组学的生存期预测方法及装置 Download PDFInfo
- Publication number
- CN106815481B CN106815481B CN201710038465.XA CN201710038465A CN106815481B CN 106815481 B CN106815481 B CN 106815481B CN 201710038465 A CN201710038465 A CN 201710038465A CN 106815481 B CN106815481 B CN 106815481B
- Authority
- CN
- China
- Prior art keywords
- patient
- image
- region
- sub
- tumor
- 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
Images
Classifications
-
- 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
- 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/70—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for mining of medical data, e.g. analysing previous cases of other patients
-
- 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/30016—Brain
-
- 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)
- Medical Informatics (AREA)
- Public Health (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Epidemiology (AREA)
- Pathology (AREA)
- Databases & Information Systems (AREA)
- Biomedical Technology (AREA)
- Primary Health Care (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
Abstract
本申请提供了一种基于影像组学的生存期预测方法及装置,其中,方法包括获取多个患者的影像数据及各患者的生存期;从各患者的影像数据中分别分割出各患者的肿瘤区域及肿瘤区域的子区域;对各患者的肿瘤区域及肿瘤区域的子区域进行影像特征提取;将各患者影像特征中的冗余特征去除以得到各患者筛选后的影像特征;根据各患者的生存期及各患者筛选后的影像特征得到影像特征与生存期之间的关联关系。本申请通过划分出肿瘤区域及肿瘤区域的子区域,对肿瘤区域及肿瘤区域的子区域进行特征提取,从而获得大量的影像特征,为寻求影像特征和患者生存期的关系提供更多有力支持。
Description
技术领域
本申请属于生物医学工程领域,特别涉及一种基于影像组学的生存期预测方法及装置。
背景技术
近年来,随着模式识别工具的增加和肿瘤个体化治疗的发展,影像组学应运而生。影像组学从影像、病理、基因等海量数据中利用自动化高通量的数据特征提取算法将影像数据转化为具有高分辨率的可挖掘的影像特征数据。通过这些数据挖掘描述组织特性。有研究报道,影像组学数据可以判断组织特性对治疗的反应,并预测患者的预后。应用影像组学数据对患者或肿瘤表型进行精确区分,可以成为对临床参数的一种有效补充。
影像组学是一种可以根据患者医学影像中的肿瘤表型对患者进行分类的方法。它通过对肿瘤图像提取大量先进的影像特征来评估肿瘤表型,然后使用可靠,重现性好的方法论,将这些特征与临床结果进行综合分析作为潜在的预后指标,从而提供一种非侵入性的精准诊疗方法。影像组学生成一个独特的肿瘤数据集,它是一个肿瘤表型的量化,可以提供比目前临床应用的成像指标更高的预测能力。
现有技术采用影像组学方法预测患者生存期的方法中分割出了完整的肿瘤区域,根据该完整的肿瘤区域提取了影像特征,建立了完整肿瘤区域影像特征和生存期之间的关联性。
发明内容
现有技术中,基于影像组学的分析只分割出了完整的肿瘤区域,没有考虑肿瘤区域的子区域具有不同的病理特征,另只提取了肿瘤区域的影像特征,忽略了肿瘤区域的子区域具有不同的病理的影像特征,从而导致所提取的影像特征不能全面的代表感兴趣区域的特点,极大的限制了基于感兴趣区域所提取的影像特征的数量与质量的问题。
为了解决上述技术问题,本申请的一技术方案为提供一种基于影像组学的生存期预测方法,包括:
获取多个患者的影像数据及各患者的生存期;
从各患者的影像数据中分别分割出各患者的肿瘤区域及肿瘤区域的子区域;
对各患者的肿瘤区域及肿瘤区域的子区域进行影像特征提取;
将各患者影像特征中的冗余特征去除以得到各患者筛选后的影像特征;
根据各患者的生存期及各患者筛选后的影像特征得到影像特征与生存期之间的关联关系。
本申请另一技术方案为提供一种基于影像组学的生存期预测装置,包括:
获取模块,用于获取多个患者的影像数据及各患者的生存期;
分割模块,用于从各患者的影像数据中分别分割出各患者的肿瘤区域及肿瘤区域的子区域;
特征提取模块,用于对各患者的肿瘤区域及肿瘤区域的子区域进行影像特征提取;
筛选模块,用于将各患者影像特征中的冗余特征去除以得到各患者筛选后的影像特征;
关联模块,用于根据各患者的生存期及各患者筛选后的影像特征得到影像特征与生存期之间的关联关系。
本申请提供的基于影像组学的生存期预测方法及装置能够细化肿瘤区域得到肿瘤区域的子区域,对肿瘤区域及肿瘤区域的子区域进行特征提取,从而获得大量的影像特征,为寻求影像特征和(患者)生存期的关系提供更多有力支持。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本申请实施例的基于影像组学的患者生存期预测方法的流程图;
图2为本申请实施例的基于影像组学的患者生存期预测装置的结构图;
图3为本申请实施例的基于影像组学的患者生存期预测装置的结构图;
图4为本申请实施例的脑瘤区域的子区域划分的示意图。
具体实施方式
为了使本申请的技术特点及效果更加明显,下面结合附图对本申请的技术方案做进一步说明,本申请也可有其他不同的具体实例来加以说明或实施,任何本领域技术人员在权利要求范围内做的等同变换均属于本申请的保护范畴。
在本说明书的描述中,参考术语“一个实施例”、“一个具体实施例”、“一些实施例”、“例如”、“示例”、“具体示例”或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本申请的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。各实施例中涉及的步骤顺序用于示意性说明本申请的实施,其中的步骤顺序不作限定,可根据需要作适当调整。
如图1所示,图1为本申请实施例的基于影像组学的患者生存期预测方法的流程图。本实施例中,考虑了肿瘤组织空间-时间异质性,提取出的感兴趣区域不局限于肿瘤区域,还提取出了肿瘤区域的子区域,并对肿瘤区域的子区域进行了影像特征提取,使得获得的影像特征更全面,建立的影像特征与生存期之间的关联关系精度更高,能够提高患者生存期预测的准确度。具体的,包括:
步骤101:获取多个患者的影像数据及各患者的生存期。
多个患者的影像数据为患者的影像集合,例如表示为V={vi,i=1,…,N},N为患者的个数,每个体数据vi表示一个影像样本。
本步骤中所述的多个患者患有同一种肿瘤疾病,如脑瘤。影像数据为患者进行治疗前的影像图像,包括但不限于包括PET,CT及MRI等影像。生存期为影像数据产生至患者死亡之间的时间段。实施时,可从医院获取多个患者的影像数据及各患者的生存期。
本步骤中获得的患者的影像数据例如来源于TCIA(The Cancer ImagingArchive,癌症影像存档)中的57例患者影像数据集。本申请对患者的影像数据个数不做限定,一般情况下,个数越多,关联关系确定越准确。
步骤102:从各患者的影像数据中分别分割出各患者的肿瘤区域及肿瘤区域的子区域,肿瘤区域为整个肿瘤区域,肿瘤区域的子区域为整个肿瘤区域内部划分的子区域。以脑瘤为例,脑瘤区域的子区域包括但不限于坏死区、增强区及水肿区,划分结果如图4所示。
肿瘤区域及肿瘤区域的子区域构成感兴趣区域,肿瘤区域除了包括真实肿瘤区域外,还包括疑似肿瘤区域。
病理学已证实:肿瘤水肿区内有异常形态的毛细血管,间质水肿和散在的肿瘤细胞在新生血管或扩张血管的浸润生长;坏死区则是因肿瘤生长过快,养料供应不够,导致肿瘤内部坏死,它可以间接的反应肿瘤的增长速度。将肿瘤区域的子区域作为感兴趣区域可以详细的具有代表性的反应肿瘤所具有的性质特征,为下一步的高通量特征提取提供了更大的提取区域。
步骤103:对各患者的肿瘤区域及肿瘤区域的子区域进行影像特征提取。
肿瘤区域的子区域有其独特的病理特征,提取这些区域的特征可更全面的反应出肿瘤的性质,同时,进一步获得大量的影像学特征数据以极大的提高后期预测患者生存期的准确性。
肿瘤区域、肿瘤区域的子区域提取的影像特征可以相同,也可以不同,本申请对此不作具体限定。
步骤104:将各患者影像特征中的冗余特征去除以得到各患者筛选后的影像特征。冗余特征包括可以通过其他特征推演出来的特征及与预测无关的特征。
步骤105:根据各患者的生存期及各患者筛选后的影像特征得到影像特征与生存期之间的关联关系。
本实施例能够细化肿瘤区域得到肿瘤区域的子区域,对肿瘤区域及肿瘤区域的子区域进行影像特征提取,从而获得大量的影像特征,为寻求影像特征和(患者)生存期的关系提供更多有力支持。
一些实施方式中,上述步骤102之前还包括:对各患者的影像数据进行预处理。
因患者的影像数据获取参数不同,通过本步骤的处理能够标准化患者的影像数据。具体的,对各患者的影像数据进行预处理包括对各患者的影像数据进行图像配准、图像平滑及数据标准化处理。下面将对预处理过程进行详细的说明:
图像配准:将不同时间,不同成像设备或不同条件(如照度,摄像位置和角度)获取的两幅或多幅影像进行匹配,叠加,从而统一各患者影像数据的坐标系。
具体的配准方法有相对配准和绝对配准。相对配准是指选择多影像中的一张影像作为参考影像,将其它的相关影像与之配准,配准时可任意选择坐标系统。绝对配准是指先定义一个控制网格,所有的影像相对于该控制网络来进行配准,也就是分别完成各影像的几何校正来实现坐标系的统一。具体实施时,可根据需求选择图像配准方法,本申请对此不作具体限定。
图像平滑:在影像数据获取、形态学处理过程中,可能会产生不平滑的毛刺、锋利的边缘等情况。这通常不是物体原始特征,而是人为造成的,需要去除,图像平滑用于去除各患者影像数据中的毛刺。影像数据平滑的方法有很多种,例如均值滤波器,中值滤波器等,这些方法均可通过现有技术实现,此处不再详细叙述。
数据标准化:标准化即将获取的多个患者的影像数据变为统一标准的影像数据。实施时,影像数据标准化的方法有多种,例如最小-最大标准化,Z-score标准化,小数定标标准化,线性变换等等。这些方法均可通过现有技术实现,下面以线性变换为例进行说明。
设样本A为一患者的影像数据,maxA为样本A像素的最大值,A为样本A的任意一原像素值,A′为与A相对的标准化的像素值,M为所有患者影像数据像素最大值。
其公式为:A′=(A/maxA)×(M+100)。
一些实施方式中,上述步骤103中提取的影像特征包括但不限于直方图特征数据,形状特征数据及纹理特征数据。
直方图特征用于描述肿瘤区域及各个肿瘤区域的子区域的所有像素点的灰度值,包括均值,中值,最大值,最小值,极差,能量,熵,偏斜度,峰度,标准差,方差,平均绝对差,均方根等等。
形状特征用于描述肿瘤区域及各个肿瘤区域的子区域的三维特征,包括体积,最长径,表面积,硬度,密度,等周商,球形不均衡度,曲率,偏心率,表面积体积比等等。
纹理特征用于描述肿瘤区域及各个肿瘤区域的子区域的纹理特征来量化肿瘤内部的异质性,包括灰度共生矩阵类特征,灰度游程矩阵类特征,灰度尺寸区域矩阵类特征,邻域灰度调差矩阵类特征,小波变换,拉普拉斯变换,高斯变换等等。
直方图特征数据,形状特征数据及纹理特征数据具体包括的特征参见表一。需要说明的是,表一中的所有特征并非都得提取,可根据肿瘤类型进行选取,例如,脑瘤分区后,对脑瘤整个区域及内部子区域分别提取14个直方图特征,28个形状特征,52个纹理特征。
表一
上述步骤104可以减轻后期学习中的维数灾难问题,也可以降低后期学习的难度,实施时,可采用现有的方法去除各患者影像特征中的冗余特征,具体的方法包括:
1)广度优先搜索:枚举了所有的特征组合,广度优先遍历特征子空间,即从一个顶点开始,辐射状地优先遍历其周围较广的特征子空间;
2)分支限界搜索:在穷举搜索的基础上加入分支限界。例如:若断定某些分支不可能搜索出比当前找到的最优解更优的解,则可以剪掉这些分支。
3)定向搜索:首先选择N个得分最高的特征作为特征子集,将其加入一个限制最大长度的优先队列,每次从队列中取出得分最高的子集,然后穷举向该子集加入1个特征后产生的所有特征集,将这些特征集加入队列。特征的得分即特征与生存期的一致性指标(CI:Harrell’s Concordance Index)值的大小,CI值的高低由所选特征和生存期的相关性有关,相关性越高,特征的分数越高。
4)最优优先搜索:与定向搜索类似,唯一的不同点是不限制优先队列的长度。
5)序列前向选择:特征子集X从空集开始,每次选择一个特征x加入特征子集X,使得特征函数J(X)最优。简单说就是,每次都选择一个使得评价函数的取值达到最优的特征加入,其实就是一种简单的贪心算法。
6)序列后向选择:从特征全集O开始,每次从特征集O中剔除一个特征x,使得剔除特征x后评价函数值达到最优。
7)增L去R选择算法:该算法有两种形式:(1)算法从空集开始,每轮先加入L个特征,然后从中去除R个特征,使得评价函数值最优。(L>R);(2)算法从全集开始,每轮先去除R个特征,然后加入L个特征,使得评价函数值最优。(L<R)。
8)序列浮动选择:序列浮动选择由增L去R选择算法发展而来,该算法与增L去R选择算法的不同之处在于:序列浮动选择的L与R不是固定的,而是“浮动”的,也就是会变化的。序列浮动选择根据搜索方向的不同,有以下两种变种。(1)序列浮动前向选择:空集开始,每轮在未选择的特征中选择一个子集x,使加入子集x后评价函数达到最优,然后在已选择的特征中选择子集z,使剔除子集z后评价函数达到最优。(2)序列浮动后向选择:与SFFS类似,不同之处在于SFBS是从全集开始,每轮先剔除特征,然后加入特征。
9)决策树:在训练样本集上运行C4.5或其他决策树生成算法,待决策树充分生长后,再在树上运行剪枝算法。则最终决策树各分支处的特征就是选出来的特征子集了。决策树方法一般使用信息增益作为评价函数。
10)遗传算法:首先随机产生一批特征子集,并用评价函数给这些特征子集评分,然后通过交叉、突变等操作繁殖出下一代的特征子集,并且评分越高的特征子集被选中参加繁殖的概率越高。这样经过N代的繁殖和优胜劣汰后,种群中就可能产生了评价函数值最高的特征子集。
11)聚类算法:基于K均值聚类的特征选择方法,其基本思想就是对每一个特征子集利用K均值聚类算法确定其最佳分类数,从而删除掉相关性较大的特征之一。本发明采用的是一致性聚类方法:计算在多次运行的情况下两个子样本聚集在一起的频率,并根据一致率的结果做出视觉评价:比较稳定性和确定最优聚类数(K)。基本假设是,如果存在最优聚类数K,K所对应的子样本稳定性会达到最佳。一致性聚类是用(0-1)来描述稳定。我们用层次聚类算法与基于皮尔森相关的不同性进行测量,经过2000次重采样迭代,为了确定最佳聚类数,我们首先计算不同聚类数的累积分布函数。最佳聚类数即为累积分布函数中函数收敛时所对应的数值。
一些实施方式中,上述步骤102从各患者的影像数据中分别分割出各患者的肿瘤区域及肿瘤区域的子区域的过程包括:
从预先分割好的影像数据中提取影像特征集;将该提取出的影像特征集训练分类器,通过该分类器来分割(预处理后的)各患者的影像数据以得到肿瘤区域及肿瘤区域的子区域。
实施时,预先分割好的影像数据可由有经验的专家按照肿瘤区域及肿瘤区域的子区域的解剖结构将患者的影像数据划分成肿瘤区域和肿瘤区域的子区域。以脑瘤为例,分割好的脑部影像包括水肿区域、增强区域、坏死区域及非增强区域。
提取的影像特征集包括但不限于直方图特征,形状特征,纹理特征。直方图特征,形状特征及纹理特征具体包括的内容如表一。影像特征集具体包括的特征由预先分割好的影像数据各区域确定。对各分割出的子区域分别提取特征,例如提取的特征包括23个低阶特征:14个灰度值特征(每个像素点和与它相邻的6个像素点的灰度值及以这7个像素点为中心用3×3模块提取的均值灰度值);6个一阶纹理特征(均值,方差,偏斜度,峰值,能量,熵);3个位置特征(X,Y,Z三个方向)。
本实施方式中所述的学习方法包括但不限于支持向量机,随机森林,卷积神经网络等等,本申请对学习方法不做具体限定。
具体实施时,上述步骤102还可采用手动分割方式实现肿瘤区域及肿瘤区域的子区域的分割。例如,让有经验的专家按照肿瘤区域及肿瘤区域的子区域的解剖结构将患者的影像划分为肿瘤区域及肿瘤区域的子区域。亦或采样人工和计算机处理结合的方式进行分割处理。
一实施方式中,上述步骤102之后还包括:针对每一患者,对该患者的部分肿瘤区域的子区域进行合并处理以得到该患者的肿瘤区域的合并子区域。
实施时,肿瘤区域的子区域合并方式有很多种,例如随机合并、具有相似病理意义的肿瘤区域的子区域合并,本申请对具体的合并方法不作限定。以脑部胶质母细胞瘤为例,该胶质母细胞瘤包括坏死区、增强区、非增强区及水肿区,因坏死区、增强区和非增强区具有相似病理意义(致命性),所以可以将坏死区、增强区和非增强区进行合并。
上述步骤103对各患者的肿瘤区域及肿瘤区域的子区域进行影像特征提取进一步为:对各患者的肿瘤区域、肿瘤区域的子区域及肿瘤区域的合并子区域进行影像特征提取。
本实施方式更适用于实际需求,提取肿瘤区域的合并子区域的特征可以更全面的、多角度的反应出肿瘤的性质。同时这一步获得的大量影像学特征可以极大的提高后期预测患者生存期的准确性。
一些实施方式中,上述步骤105得到的影像特征与生存期之间的关联关系可以表示为:生存期为筛选后的多个影像特征或全部筛选后的影像特征进行线性组合的结果。假设每个影像特征代表一个因子,则生存期就是一个多因子方程,该方程的构建流程如下:
步骤201:随机确认一个生存期的方程,该方程可以是单因子或多因子,线性或非线性,本申请对此不作限定。例如,方程可表示为:
Y=W*X+W1*X1+...+Wi+Xi+...+Wn*Xn+z,
Y表示生存期;X为肿瘤区域筛选出的影像特征;m0为肿瘤区域筛选出的影像特征个数;W为X的系数,为未知量;X1、Xi及Xn为肿瘤区域的子区域筛选出的影像特征;m1、mi及mn为肿瘤区域的子区域筛选出的影像特征个数;W1、Wi及Wn分别为X1、Xi及Xn的系数,为未知量;z表示常数,为未知量。W、W1、Wi及Wn为常数系数。
具体构建方程时,可选用多参数COX逻辑回归模型构建一个线性的多因子方程,并用梯度下降的算法寻求上述方程的最优参数组合,即计算得到未知量的值。
步骤202:分别将各患者筛选后的影像特征值及生存期代入步骤201中的方程,建模得到影像特征的系数W={Wi,i=1,…,m}及常数z,将得到的影像特征的系数及常数z代入至步骤201中得到构建好的生存期的方程。
进一步的,构建好的生存期的方程还可由时间-生存概率曲线表示,该时间-生存概率曲线的生成过程包括:
1)对于步骤101中多个患者,分别通过构建好的方程计算方程值。例如,步骤101中包括59个病例样本,即包括59个患者,则需要将这59个病例筛选后的影像特征代入构建好的方程中,共得到59个方程值。
2)提取计算出的方程值中的中间值,将该中间值为界,计算出的方程值若大于或等于中间值则赋为1,计算出的方程值若小于中间值则赋为0。
同时,记录患者的生存状态,若最后随访时患者已死亡,则该患者生存状态赋值为1,若最后随访时患者活着或者失访,则患者的生存状态赋值为0。
继续步骤1)中的例子,假设中间值对应第3个患者,则步骤2)得到的结果如表二,
表二:
3)将各患者方程值的赋值及其对应患者的生存期、患者生存状态作为输入,通过SPSS软件绘制得到时间-生存概率的曲线。
基于同一发明构思,本申请实施例中还提供了一种基于影像组学的患者生存期预测装置,如下面的实施例所述。由于该装置解决问题的原理与基于影像组学的患者生存期预测方法相似,因此该装置的实施可以参见方法的实施,重复之处不再赘述。
如图2所示,图2为本申请实施例的基于影像组学的患者生存期预测装置的结构图。该装置可以通过逻辑电路实现运行于智能终端,例如手机、平板电脑、计算机等设备中,或者以功能模块的方式由软件实现各部件的功能,运行于所述智能终端上。具体的,该装置包括:
获取模块201,用于获取多个患者的影像数据及各患者的生存期。
分割模块202,用于从各患者的影像数据中分别分割出各患者的肿瘤区域及肿瘤区域的子区域。详细的说,分割模块203具体用于从预先分割好的影像数据中提取影像特征;将该提取出的影像特征作为训练集,通过学习方法来分割预处理后的各患者的影像数据以得到肿瘤区域及肿瘤区域的子区域。
特征提取模块203,用于对各患者的肿瘤区域及肿瘤区域的子区域进行影像特征提取。详细的说,影像特征包括直方图特征数据,形状特征数据及纹理特征数据。
筛选模块204,用于将各患者影像特征中的冗余特征去除以得到各患者筛选后的影像特征。
关联模块205,用于根据各患者的生存期及各患者筛选后的影像特征得到影像特征与生存期之间的关联关系。
本申请提供的基于影像组学的生存期预测装置能够细化肿瘤区域得到肿瘤区域的子区域,对肿瘤区域及肿瘤区域的子区域进行影像特征提取,从而获得大量的影像特征,为寻求影像特征和(患者)生存期的关系提供更多有力支持。
一些实施方式中,如图3所示,基于影像组学的患者生存期预测装置还包括预处理模块206,用于对各患者的影像数据进行预处理。预处理模块206具体用于对各患者的影像数据进行图像配准、图像平滑及数据标准化处理。
一些实施方式中,复请参阅图3,基于影像组学的患者生存期预测装置还包括合并模块207,用于针对每一患者,对该患者的部分肿瘤区域的子区域进行合并处理以得到该患者的肿瘤区域的合并子区域。
特征提取模块203进一步用于对各患者的肿瘤区域、肿瘤区域的子区域及肿瘤区域的合并子区域进行影像特征提取。
一些实施方式中,关联模块205得到的影像特征与生存期之间的关联关系的表达公式为:
Y=W*X+W1*X1+...+Wi+Xi+...+Wn*Xn+z,
Y表示生存期,X表示肿瘤区域筛选出的影像特征,m0表示肿瘤区域筛选出的影像特征个数,W表示X的系数,X1、Xi及Xn表示肿瘤区域的子区域筛选出的影像特征,m1、mi及mn表示肿瘤区域的子区域筛选出的影像特征个数,W1、Wi及Wn分别表示X1、Xi及Xn的系数,z表示常数。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述仅用于说明本申请的技术方案,任何本领域普通技术人员均可在不违背本申请的精神及范畴下,对上述实施例进行修饰与改变。因此,本申请的权利保护范围应视权利要求范围为准。
Claims (8)
1.一种基于影像组学的患者生存期预测方法,其特征在于,包括:
获取多个患者的影像数据及各患者的生存期;
从各患者的影像数据中分别分割出各患者的肿瘤区域及肿瘤区域的子区域;
对各患者的肿瘤区域及肿瘤区域的子区域进行影像特征提取;
将各患者影像特征中的冗余特征去除以得到各患者筛选后的影像特征;
根据各患者的生存期及各患者筛选后的影像特征得到影像特征与生存期之间的关联关系,表示为:
Y=W*X+W1*X1+...+Wi+Xi+...+Wn*Xn+z,
Y表示生存期,X表示肿瘤区域筛选出的影像特征,m0表示肿瘤区域筛选出的影像特征个数,W表示X的系数,X1、Xi及Xn表示肿瘤区域的子区域筛选出的影像特征,m1、mi及mn表示肿瘤区域的子区域筛选出的影像特征个数,W1、Wi及Wn分别表示X1、Xi及Xn的系数,z表示常数。
2.如权利要求1所述的基于影像组学的患者生存期预测方法,其特征在于,从各患者的影像数据中分别分割出各患者的肿瘤区域及肿瘤区域的子区域之前还包括:
对各患者的影像数据进行预处理,其中,预处理包括图像配准、图像平滑及数据标准化处理。
3.如权利要求1所述的基于影像组学的患者生存期预测方法,其特征在于,所述影像特征包括直方图特征数据,形状特征数据及纹理特征数据。
4.如权利要求1所述的基于影像组学的患者生存期预测方法,其特征在于,从各患者的影像数据中分别分割出各患者的肿瘤区域及肿瘤区域的子区域之后还包括:
针对每一患者,对该患者的部分肿瘤区域的子区域进行合并处理以得到该患者的肿瘤区域的合并子区域;
对各患者的肿瘤区域及肿瘤区域的子区域进行影像特征提取进一步为对各患者的肿瘤区域、肿瘤区域的子区域及肿瘤区域的合并子区域进行影像特征提取。
5.一种基于影像组学的患者生存期预测装置,其特征在于,包括:
获取模块,用于获取多个患者的影像数据及各患者的生存期;
分割模块,用于从各患者的影像数据中分别分割出各患者的肿瘤区域及肿瘤区域的子区域;
特征提取模块,用于对各患者的肿瘤区域及肿瘤区域的子区域进行影像特征提取;
筛选模块,用于将各患者影像特征中的冗余特征去除以得到各患者筛选后的影像特征;
关联模块,用于根据各患者的生存期及各患者筛选后的影像特征得到影像特征与生存期之间的关联关系,表示为:
Y=W*X+W1*X1+...+Wi+Xi+...+Wn*Xn+z,
Y表示生存期,X表示肿瘤区域筛选出的影像特征,m0表示肿瘤区域筛选出的影像特征个数,W表示X的系数,X1、Xi及Xn表示肿瘤区域的子区域筛选出的影像特征,m1、mi及mn表示肿瘤区域的子区域筛选出的影像特征个数,W1、Wi及Wn分别表示X1、Xi及Xn的系数,z表示常数。
6.如权利要求5所述的基于影像组学的患者生存期预测装置,其特征在于,还包括预处理模块,用于对各患者的影像数据进行预处理,其中,预处理包括图像配准、图像平滑及数据标准化处理。
7.如权利要求5所述的基于影像组学的患者生存期预测装置,其特征在于,所述影像特征包括直方图特征数据,形状特征数据及纹理特征数据。
8.如权利要求5所述的基于影像组学的患者生存期预测装置,其特征在于,还包括合并模块,用于针对每一患者,对该患者的部分肿瘤区域的子区域进行合并处理以得到该患者的肿瘤区域的合并子区域;
所述特征提取模块进一步用于对各患者的肿瘤区域、肿瘤区域的子区域及肿瘤区域的合并子区域进行影像特征提取。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710038465.XA CN106815481B (zh) | 2017-01-19 | 2017-01-19 | 一种基于影像组学的生存期预测方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710038465.XA CN106815481B (zh) | 2017-01-19 | 2017-01-19 | 一种基于影像组学的生存期预测方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106815481A CN106815481A (zh) | 2017-06-09 |
CN106815481B true CN106815481B (zh) | 2020-07-17 |
Family
ID=59112288
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710038465.XA Active CN106815481B (zh) | 2017-01-19 | 2017-01-19 | 一种基于影像组学的生存期预测方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106815481B (zh) |
Families Citing this family (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108090507A (zh) * | 2017-10-19 | 2018-05-29 | 电子科技大学 | 一种基于集成方法的医疗影像纹理特征处理方法 |
CN108319969B (zh) * | 2018-01-12 | 2021-06-22 | 复旦大学 | 基于稀疏表示框架的脑胶质瘤生存期预测方法和系统 |
CN108389109B (zh) * | 2018-02-11 | 2022-02-01 | 中国民航信息网络股份有限公司 | 一种基于混合特征选择算法的民航可疑订单特征提取方法 |
CN108897984A (zh) * | 2018-05-07 | 2018-11-27 | 上海理工大学 | 基于ct影像组学特征与肺癌基因表达间相关性分析方法 |
CN108766563A (zh) * | 2018-05-25 | 2018-11-06 | 戴建荣 | 基于剂量组学的放射治疗结果预测方法和系统 |
CN108921821A (zh) * | 2018-06-01 | 2018-11-30 | 中国人民解放军战略支援部队信息工程大学 | 基于lasso回归的乳腺癌腋窝淋巴结转移状态的判别方法 |
CN110197236A (zh) * | 2018-10-17 | 2019-09-03 | 中山大学附属第三医院 | 一种基于影像组学对糖皮质激素疗效的预测和验证方法 |
CN109686426A (zh) * | 2018-12-29 | 2019-04-26 | 上海商汤智能科技有限公司 | 医疗影像处理方法及装置、电子设备及存储介质 |
CN109872394B (zh) * | 2019-01-10 | 2021-04-27 | 重庆大学 | 基于最小二乘支持向量机的狭长三角形网格优化方法 |
CN109949268A (zh) * | 2019-01-24 | 2019-06-28 | 郑州大学第一附属医院 | 一种基于机器学习的肝细胞癌分化水平分级方法 |
CN110148115A (zh) * | 2019-04-04 | 2019-08-20 | 中国科学院深圳先进技术研究院 | 一种癌转移预测影像特征的筛选方法、装置和存储介质 |
CN110211089A (zh) * | 2019-04-23 | 2019-09-06 | 中国科学院深圳先进技术研究院 | 肾透明细胞癌转移判断特征的提取方法、装置和存储介质 |
CN110391015B (zh) * | 2019-06-14 | 2021-08-13 | 广东省人民医院(广东省医学科学院) | 一种基于影像组学量化肿瘤免疫状态的方法 |
CN112750528A (zh) * | 2019-10-30 | 2021-05-04 | 中国医药大学附设医院 | 预测肿瘤的特性参数的计算机辅助预测系统、方法及计算机程序产品 |
CN110993106A (zh) * | 2019-12-11 | 2020-04-10 | 深圳市华嘉生物智能科技有限公司 | 一种结合病理图像和临床信息的肝癌术后复发风险预测方法 |
CN111626980A (zh) * | 2020-03-27 | 2020-09-04 | 杭州依图医疗技术有限公司 | 一种医学影像的显示方法及装置、计算设备及存储介质 |
CN111798465A (zh) * | 2020-07-02 | 2020-10-20 | 中国人民解放军空军军医大学 | 基于医学影像的异质肿瘤高危区域检测方法及系统 |
CN111814868A (zh) * | 2020-07-03 | 2020-10-23 | 苏州动影信息科技有限公司 | 一种基于影像组学特征选择的模型、构建方法和应用 |
CN112329876A (zh) * | 2020-11-16 | 2021-02-05 | 中山大学附属第六医院 | 一种基于影像组学的结直肠癌预后预测方法及设备 |
CN112907555B (zh) * | 2021-03-11 | 2023-01-17 | 中国科学院深圳先进技术研究院 | 一种基于影像基因组学的生存预测方法和系统 |
CN113643805B (zh) * | 2021-08-10 | 2024-04-02 | 复旦大学附属华山医院 | 基于影像组学的脑膜瘤伽马刀治疗后水肿预测系统 |
CN116228753B (zh) * | 2023-05-06 | 2023-08-01 | 中山大学孙逸仙纪念医院 | 肿瘤预后评估方法、装置、计算机设备和存储介质 |
CN116385441B (zh) * | 2023-06-05 | 2023-09-05 | 中国科学院深圳先进技术研究院 | 一种基于mri对少突胶质细胞瘤进行风险分层的方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6182069B1 (en) * | 1992-11-09 | 2001-01-30 | International Business Machines Corporation | Video query system and method |
CN103973968A (zh) * | 2013-02-01 | 2014-08-06 | 宏达国际电子股份有限公司 | 电子装置及其影像合成方法 |
CN104698168A (zh) * | 2008-09-01 | 2015-06-10 | 国立大学法人名古屋大学 | 抗癌剂作用增强剂及其应用、以及癌患者预后推定用生物标记物及其应用 |
CN104878012A (zh) * | 2015-05-22 | 2015-09-02 | 中国人民解放军第二军医大学 | Hsa-miR-3200-5p在制备早期筛查或诊断Brachyury阳性肿瘤试剂或试剂盒中的应用 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1789923A1 (en) * | 2004-08-11 | 2007-05-30 | Aureon Laboratories, Inc. | Systems and methods for automated diagnosis and grading of tissue images |
CN101087358A (zh) * | 2006-06-07 | 2007-12-12 | 英业达股份有限公司 | 影像处理系统及方法 |
TW201007451A (en) * | 2008-08-12 | 2010-02-16 | Phison Electronics Corp | Block management method for flash memory and storage system and controller using the same |
CN102855483B (zh) * | 2011-06-30 | 2017-09-12 | 北京三星通信技术研究有限公司 | 处理超声图像的方法和装置以及乳腺癌诊断设备 |
US9974512B2 (en) * | 2014-03-13 | 2018-05-22 | Convergence Medical, Llc | Method, system, and computer program product for determining a patient radiation and diagnostic study score |
CN104143101A (zh) * | 2014-07-01 | 2014-11-12 | 华南理工大学 | 一种基于超声图像的乳腺肿瘤区域自动识别方法 |
CN104881568A (zh) * | 2015-04-27 | 2015-09-02 | 苏州敏宇医疗科技有限公司 | 一种基于云计算的肿瘤治疗疗效早期评价的系统和方法 |
CN105653858A (zh) * | 2015-12-31 | 2016-06-08 | 中国科学院自动化研究所 | 一种基于影像组学的病变组织辅助预后系统和方法 |
-
2017
- 2017-01-19 CN CN201710038465.XA patent/CN106815481B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6182069B1 (en) * | 1992-11-09 | 2001-01-30 | International Business Machines Corporation | Video query system and method |
CN104698168A (zh) * | 2008-09-01 | 2015-06-10 | 国立大学法人名古屋大学 | 抗癌剂作用增强剂及其应用、以及癌患者预后推定用生物标记物及其应用 |
CN103973968A (zh) * | 2013-02-01 | 2014-08-06 | 宏达国际电子股份有限公司 | 电子装置及其影像合成方法 |
CN104878012A (zh) * | 2015-05-22 | 2015-09-02 | 中国人民解放军第二军医大学 | Hsa-miR-3200-5p在制备早期筛查或诊断Brachyury阳性肿瘤试剂或试剂盒中的应用 |
Also Published As
Publication number | Publication date |
---|---|
CN106815481A (zh) | 2017-06-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106815481B (zh) | 一种基于影像组学的生存期预测方法及装置 | |
Kallenberg et al. | Unsupervised deep learning applied to breast density segmentation and mammographic risk scoring | |
CN110110617B (zh) | 医学影像分割方法、装置、电子设备和存储介质 | |
CN104751178B (zh) | 基于形状模板匹配结合分类器的肺结节检测装置及方法 | |
CN112768072B (zh) | 基于影像组学定性算法构建癌症临床指标评估系统 | |
CN110956626B (zh) | 一种基于图像的预后评估方法及装置 | |
Aranguren et al. | Improving the segmentation of magnetic resonance brain images using the LSHADE optimization algorithm | |
JP6168426B2 (ja) | 疾患分析装置、制御方法、及びプログラム | |
CN112927799B (zh) | 融合多示例学习和多任务深度影像组学的生存期分析系统 | |
US20230056839A1 (en) | Cancer prognosis | |
CN111798424A (zh) | 一种基于医学图像的结节检测方法、装置及电子设备 | |
EP3817648A1 (en) | Method for diagnosing, predicting, determining prognosis, monitoring, or staging disease based on vascularization patterns | |
CN114332132A (zh) | 图像分割方法、装置和计算机设备 | |
Messadi et al. | Extraction of specific parameters for skin tumour classification | |
Meswal et al. | A weighted ensemble transfer learning approach for melanoma classification from skin lesion images | |
Ma et al. | Automatic pulmonary ground‐glass opacity nodules detection and classification based on 3D neural network | |
CN116563651A (zh) | 鼻咽癌预后特征确定方法、系统、装置及存储介质 | |
CN116703880A (zh) | 一种免疫状态预测模型构建方法、预测方法及构建装置 | |
CN113763343B (zh) | 基于深度学习的阿兹海默症检测方法及计算机可读介质 | |
WO2018132997A1 (zh) | 一种基于影像组学的生存期预测方法及装置 | |
US20240087751A1 (en) | Systems and methods for organ shape analysis for disease diagnosis and risk assessment | |
Demidova et al. | Skin cancer classification computer system development with deep learning. | |
Sharma et al. | Improved Segmentation algorithm using PSO and K-means for Basal Cell Carcinoma Classification from Skin Lesions | |
Jafari et al. | Nuclei Segmentation in Hematoxylin & Eosin-stained Tissue Images Using U-Net-based Models | |
Jonaitytė | Survival analysis incorporating medical imaging data |
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 |