基于多参数磁共振影像的前列腺癌计算机辅助识别系统的制作方法
【专利摘要】本发明公开了一种基于多参数磁共振影像的前列腺癌计算机辅助识别系统及方法，包括图像预处理模块、参数处理模块、预测与评价模块三个部分。该系统综合利用了T2加权图像、扩散加权图像和动态增强图像的特征，利用人工神经网络的结构达到识别前列腺癌病灶的目的。经测试该系统在前列腺外周带的ROC(受试者工作特性)曲线下面积为0.931，识别准确率为0.887，在中央腺体中的ROC曲线下面积为0.909，识别准确率为0.915。采用本发明的技术方案，能够很好地将常规磁共振扫描序列所得的影像信息整合起来，最大限度地利用了这些图像中的定量参数，较为客观地给出了前列腺癌的识别结果。且操作简便，能较为直观地给医生一个参考，为后续制定诊断方案提供重要依据。
【专利说明】基于多参数磁共振影像的前列腺癌计算机辅助识别系统
【技术领域】
[0001]本发明属于医学图像的计算机辅助分析【技术领域】，具体是一种基于磁共振图像(MRI)中感兴趣区域(ROI)的特征参数来辅助识别前列腺癌病灶的软件系统。
【背景技术】
[0002]前列腺癌是男性中常见的一种恶性肿瘤，是男性癌症患者中第二多发的疾病，在整个人群中发病率排在第五位。近些年来，前列腺癌受到了潜在患者和泌尿科医师越来越多的关注，临床病例数也在快速积累。磁共振技术应用于前列腺疾病有十多年时间，目前已成为前列腺癌识别、诊断、指导治疗、随访最常用的影像学方法。近几年，由于磁共振技术的进步和诊断经验的积累，同时包含解剖信息和功能信息的多参数磁共振成像大大提高了磁共振在前列腺癌各方面的应用价值。
[0003]目前的前列腺磁共振检查主要包括常规检查，如Tl加权成像(TlWI)和T2加权成像(T2WI)，以及能够提供功能信息的多种磁共振功能成像，如扩散加权成像(DWI)、动态增强图像(DCE)、磁共振波谱成像(MRS)等。各个序列图像之间的信息能够相互补充，但同时临床医生需要综合利用所有的影像结果，对病人的情况做出判断，而判断正确与否依赖于临床医生的个人经验。计算机辅助识别系统可有效避免医生基于单一指标判断，或过度依赖于个人经验的主观判断，这种计算机辅助技术可以综合多个参数并提供多变量分析，不仅使医生的诊断过程更加客观，而且提高了诊断效率。
[0004]人工神经网络(ANN)以其能够以计算方法为基础进行多因子分析的特性而被广泛应用于前列腺癌的识别。与传统的统计方法或是单一指标预测相比，ANN可含多个输入参数，所有指标可以同时作为ANN的输入，它可以有效地利用这些具有非线性关系的变量，进行训练和预测，从而实现辅助诊断，大大提高了诊断效率与诊断准确率，在不同程度上减少了不必要的穿刺。
[0005]由于技术的复杂性，多参数磁共振影像的信息不能很好地传递给临床医生，不易被临床医生所解读和利用。本计算机辅助识别系统正可作为一项辅助技术，解决临床医师难以有效利用多参数磁共振成像信息的难题。将多种磁共振序列所得影像的解剖及功能信息用于前列腺癌病灶的识别，可辅助医师更敏感地检出前列腺癌。这样的计算机辅助识别系统将使得多参数磁共振信息得到充分利用，在计算机辅助识别系统的帮助下临床医生的诊断过程更加客观准确，操作起来也非常方便，使前列腺癌的各项检查更有效率。
【发明内容】
[0006]本发明提出了一种通过计算机辅助的方法来识别前列腺癌的系统。通过提取多参数磁共振图像中的定量信息，选取合适的特征参数，应用人工神经网络的方法，对前列腺癌病灶进行了识别。
[0007]本发明能够通过计算机的定量分析，完成以下任务:
[0008]1.从磁共振图像中挖掘有用信息，识别前列腺癌病灶，给出患癌概率，具有较高的准确率、敏感性、特异性、阳性预测率和阴性预测率；
[0009]2.对前列腺区分外周带和中央腺体分别进行分析，比较了二者在磁共振影像上特征的差异；
[0010]3.通过特征筛选，选取对识别前列腺癌有重要意义的特征参数，为临床诊断提供依据；
[0011]4.对比不同磁共振序列的组合效果，找出最佳的组合形式，尝试为临床扫描提供更合理的方案和参考。
[0012]为实现上述目的，本发明采取以下技术方案:基于多参数磁共振影像的前列腺癌计算机辅助识别系统，包括以下三个步骤:
[0013]1.将多参数磁共振影像的数据导入该计算机辅助识别系统，从图像中选取可疑区域；
[0014]2.提取所选图像区域中的特征参数，选取对辅助识别前列腺癌贡献较大的特征作为后续步骤的输入；
[0015]3.将步骤2中得到的特征参数作为人工神经网络的输入参数，通过调节神经网络的结构，使得输出的识别结果达到最佳。
[0016]本发明通过采取以上技术方案，很好地将常规磁共振扫描序列所得的图像信息整合起来，利用这些图像中的定量参数，客观地给出了前列腺癌的识别结果。通过比较发现，综合利用磁共振影像中三个常规序列的扫描结果能更准确地识别前列腺癌，并为临床决策提供更加丰富的信息。该技术方案操作简单，只需临床医生从扫描图像中选取可疑的前列腺癌区域，即可给出该区域的患癌概率，能给医生一个直观的参考，为后续制定诊断方案提供重要依据。
[0017]本发明已在71例病人的影像数据中做过测试，所有前列腺癌患者为穿刺验证，非前列腺癌患者同时得到穿刺排除和三年以上随访。病人在接受穿刺之前均接受过多参数磁共振扫描，包括T2WI，DffI和DCE，并经过常规临床检查。
[0018]通过对已有的71个病例分析，目前该系统的识别准确率在外周带中为0.887，在中央腺体中为0.915，均能满足临床的需求，并且通过比较发现，该计算机辅助系统的识别结果与有经验的放射学专家(分别有7年和9年前列腺癌影像诊断经验)的识别结果相当。在外周带中，通过该系统，三个扫描序列单独的识别结果与有经验放射学专家的识别结果的相关系数分别为T2W1:0.706，Dff1:0.750，DCE:0.554，在中央腺体内三个序列的识别结果与有经验放射学专家的识别结果的相关系数分别为T2W1:0.492，Dff1:0.650，DCE:0.348。
【专利附图】
【附图说明】
[0019]图1是本计算机辅助识别系统的数据处理流程图。
[0020]图2是两个病人(一个前列腺癌患者，另一个为非前列腺癌病人)的磁共振T2加权图像、扩散加权图像和动态增强图像及曲线，图中箭头所指为前列腺癌区域。
[0021]图3是人工神经网络的结构示意图。
[0022]图4是本计算机辅助识别系统对测试结果统计后得到的ROC(受试者工作特性)曲线。
【具体实施方式】
[0023]下面结合附图对本发明作进一步详细说明，本发明共包含三个主要模块，如图1所示，具体如下:
[0024]1.图像预处理模块:将来自同一病人的轴位T2加权图像、扩散加权图像和动态增强图像分别导入到系统内，手动选取感兴趣区域(可固定大小也可任意大小)，计算扩散加权图像的表观扩散系数图(ADC)和动态增强图像的动态增强曲线，为后续提取特征做准备。整个选取过程由两位有经验的放射学专家(分别有7年和9年的前列腺磁共振诊断经验)来评定，共计从外周带中选取了 238个R0I，从中央腺体中选取了 188个R0I。
[0025]2.参数处理模块:分别对三个序列扫描得到的图像提取特征参数，如图2。其中T2加权图像中提取特征参数12个(包括感兴趣区的信号均值、标准差、变异系数、最小值、10%分位数、偏态系数、峰态系数和基于灰度共生矩阵提取的能量值、对比度、相关系数、一致性、熵等)，扩散加权图像特征12个(与T2加权图像中的类似)，动态增强图像特征5个(曲线流入率、流出率、峰值、达峰时间和曲线下面积等)。对所提取的特征使用顺序选择方法进行优化，得到贡献较大的特征参数作为人工神经网络的输入参数。
[0026]3.预测与评价模块:
[0027]I)前列腺癌识别:通过人工神经网络模型达到识别前列腺癌病灶的目的。该神经网络结构如图3所示，共包括三层，分别是输入层、隐层和输出层，其中输入层各节点代表从步骤2中选择的特征参数，将这些参数归一化后作为人工神经网络的输入。该网络隐层共有节点14个(通过经验获取)，输出层为识别结果(O到I之间，O代表不可能是前列腺癌，I代表确定是前列腺癌)；
[0028]2)识别性能评价:通过交叉检验的方法轮流训练和测试数据，利用ROC曲线的曲线下面积(AUC)作为性能评价的指标，同时也计算该系统识别的准确率、敏感性、特异性、阳性预测率和阴性预测率等。
[0029]在测试阶段，按照不同的序列组合方式对数据进行了分组，分组情况如下:
[0030]第一组:T2WI
[0031]第二组:T2WI+DWI
[0032]第三组:T2WI+DCE
[0033]第四组:T2WI+DWI+DCE
[0034]按照图1所示的流程，对每组数据进行测试后发现，当结合三个序列(即第四组)时，ROC曲线下的面积达到最大，在外周带和中央腺体中分别是0.931和0.909，这个结果说明了利用人工神经网络结构来综合多参数磁共振影像中不同序列的结果是可行且非常有必要的，三个序列的影像数据能够相互补充，为临床决策提供更加丰富的信息。
【权利要求】
1.一种基于磁共振图像的计算机辅助识别系统，具体是一种利用多参数磁共振影像中感兴趣区域的特征参数来辅助识别前列腺癌病灶的软件系统，其特征是包括以下部分: 1)图像预处理:将来自同一病人的多参数磁共振影像数据导入到系统内，手动提取感兴趣区域(可固定大小也可任意大小)，计算扩散加权图像的表观扩散系数和动态增强图像的动态增强曲线； 2)参数处理:分别对三个序列扫描得到的图像提取特征参数。对所提取的特征使用顺序选择方法进行优化，得到贡献较大的特征参数作为人工神经网络的输入参数； 3)预测与评价:通过人工神经网络模型达到识别前列腺癌病灶的目的。使用交叉检验的方法轮流训练和测试数据，利用ROC曲线下面积作为性能评价的指标，同时也计算该系统识别的准确率、敏感性、特异性、阳性预测率和阴性预测率等。
2.如权利要求1所述的方法，其中还包括在对图像处理之前，使用临床常规的六分区穿刺方法获得各分区内的病理信息，为本系统的性能评价作为金标准。
3.如权利要求1所述的方法，其中多参数磁共振影像包括轴位T2加权图像、扩散加权图像和动态增强图像，以及对这些影像数据的任意组合。
4.如权利要求1所述的方法，其中对参数提取的特征包括一般灰度特征、灰度直方图特征、纹理特征和功能参数，以及对这些特征的任意组合。
5.如权利要求1所述的方法，其中利用人工神经网络模型识别前列腺癌病灶的步骤包括将参数处理完之后的信息进行归一化后输入到该网络模型中，利用机器学习的方法不断训练和测试识别结果，最后比较得出最佳的识别前列腺癌的网络参数等。
6.如权利要求1所述的方法，其中评价方法中使用的交叉检验可以是多折交叉检验方法或者是leave-one-out (从数据组中去掉一个，以所有其它的数据建立模型来预测去掉的这个)检验方法。
7.如权利要求5所述的方法，该人工神经网络模型的输出结果为O到I之间的预测概率，O代表不可能是前列腺癌，I代表确定是前列腺癌。
8.如权利要求5所述的方法，其中包括连续或定期地利用新增病例数据对所述识别系统进行在训练，补充该系统的样本量，完善神经网络的结构参数等。
【文档编号】A61B5/055GK104424386SQ201310370531
【公开日】2015年3月18日 申请日期:2013年8月23日 优先权日:2013年8月23日
【发明者】王成彦, 王鹤, 胡娟, 陈双娟, 张珏, 王霄英, 方竞 申请人:北京大学