获取高质量、标准化的图像是整个流程的基础, 解决影像特征数据采集与重建参数的标准化是关键^[9]。由于医学影像的种类繁多, 需要集中对不同模态影像数据进行整合挖掘, 为临床指标提供精确预测, 这对于影像数据的质量提出了严苛要求, 因为劣质影像会影响数据分析的准确性。数据包含的样本应尽量多, 为避免由于特征过多造成过拟合问题, 集合应达到一定规模^{[10, 11]}。

采用兴趣区(region of interest, ROI)把病变轮廓逐层勾画出来, 通过三维容积重组生成三维兴趣容积(volume of interest, OI)^[12]。图像的分割有手动、半自动、全自动三种方式。前两者是目前的主流方式, 但随着计算机技术的进步及分割算法的完善, 全自动将成为将来的主流手段^{[13, 14]}。肿瘤分割算法主要可以分为基于区域的方法, 如区域生长算法^[15]; 基于边缘约束的方法, 如主动轮廓模型和水平集分割算法; 分类或聚类方法, 包括支持向量机、深度神经网络、随机森林^[16]等; 当前也有一些混合方法, 尝试利用各种附加信息来改进分割结果。当前已经有多款较成熟的影像分割工具可以在科研工作中使用, 常用的工具包包括3D Slicer、ITK-SNAP、MITK以及中国科学院研发的3D Med和MITK等。

肿瘤ROI一旦被确定, 影像特征便可以通过计算机进行提取。特征的选择和提取也是图像定量分析的重要一步, 目的是降低图像特征描述子(descriptor)的数量进而降低计算量并提高计算效率^[17]。目前常用的影像组学特征包括一阶统计特征、空间几何特征、纹理特征和小波特征4大类^[18]。常用的纹理测度包括以下3类:①灰度共生矩阵, 其不仅包含灰度统计信息, 而且反映灰度分布的空间信息。包括角二阶矩(能量)、惯量、熵。②分形维, 图像的几何形状(粗糙程度)也可用分形维描述。分形维作为分形的重要特征和度量, 可以作为描述物体的一个稳定的特征量。③行程统计。由于纹理特征描述相似(不相似)体素间对比值的内在统计关系, 因此通过纹理分析可以测量肿瘤内的异质性。

基于上述分析提取的影像特征非常庞大。因此, 高维特征数据的降维是特征提取过程的重要一步。其中有效的降维方式有主成分分析、互信息、Fisher判别、带惩罚项回归分析和线性判别等机器学习方法, 获得的参数具有高重复性、信息丰富和冗余度低的特点^[19]。

针对研究目标完成对模型的建立和训练。常用以下六种数学方法:①多层前馈型神经网络(back-propagation, BP)法, 是目前临床科研应用最为广泛、亦最为成熟的神经网络模型。②偏最小二乘回归法(partial least squares regression, PLSR), 一种新型的多元统计数据分析方法, 它主要研究的是多因变量对多自变量的回归建模, 特别当各变量内部高度线性相关时, 用PLSR更有效。③随机森林法(Random Forest, RF), 该算法下包含多个决策树, 可以有数以千计的输入变量, 同时可以按变量对分类结果的重要性进行排序。④支持向量机(support vector machine, SVM), 其基础是统计学习理论, 采用结构风险最小化准则, 根据有限的样本信息在模型的复杂性和学习能为之间寻求最佳折衷, 获得较好的推广能力。⑤Boost算法, 其主要原理是, 通过改变数据分布的模型, 将用不同的训练集获得的弱分类器组合成一个强分类器。 其基本思路是, 根据往的训练结果(包括总体分类的准确率), 赋予每个变量相应的权值。该方法的优点是能够将重点放在那些较难分开的样本上, 且不容易出现过拟合的现象。⑥Logistic回归判别分析(logistic discriminant analysis), 相比较于一般的线性判别分析, 是一种更为灵活的判别方法, 可以处理连续型、 离散型与混合型等多种类型的数据, logistic比一般的线性判别有更广的应用范围。近年来一些更精密更复杂的分类方法, 诸如支持向量机、 随机森林等不断迅速发展^[20]。

影像组学对于鉴别良恶性病变与组织分化程度具有重要价值。Sidhu等^[21]认为前列腺移行区的纹理分析能够鉴别出前列腺移行区显著癌。Wibmer等^[22]通过对147例经活检证实的前列腺癌患者的MRI图像的研究发现, T₂加权成像(T₂-weighted imaging, T₂WI)和表观扩散系数(apparent diffusion coefficient, ADC)图像的Haralick纹理特征有利于区分癌变组织与非癌变组织:ADC 图像中, 癌变组织的熵和惯量值明显高于非癌变组织, 而能量、相关性和均质性则明显低于非癌变组织; T₂WI 图像中, 癌变组织的惯量值明显高于非癌变组织, 而相关性则明显低于非癌变组织, 而两者能量、熵和均质性均的差异无统计学意义。Liang等^[23]学者研究了基于CT图像的影像组学特征对结直肠分期的价值, 该研究结果提示影像组学分析可作为结直肠术前分期的一个重要的辅助手段。将前列腺癌的诊断带入新的高度, 提高诊断的敏感性与特异性, 以及稳定性, 大大降低依赖个人诊断水平的因素。

影像组学分析可被用于评价前列腺癌的疗效与预后分析。Rosenkrantz 等^[24]认为影像组学能够区分Gleason3+4与Gleason4+3的前列腺癌。Nie等^[25]利用治疗前多参数核磁图像(DWI、T₁/T₂、DCE)的影像组学研究评价直肠癌新辅助放化疗效果, 结果显示通过纹理参数建立的模型相较于传统的检查策略可以提高预测的准确性。Aerts等^[26]分析了440个纹理特征, 这些特征主要评价肿瘤图像的灰度、形状和纹理类参数, 这些图像来自于1 019个肺部和头颈部肿瘤; 分析显示纹理特征与潜在的基因表达相关, 并与预后生存紧密相关。Gnep等^[27]认为T₂WI Haralick 影像组学特征与前列腺癌放射治疗后的生化复发密切相关。所以随着影像组学的发展, 以及对它的认识越来越深刻, 影像组学将越来越多地应用于高级别前列腺癌的确定和前列腺癌侵袭性的分析, 以及前列腺癌内分泌治疗及放疗后的疗效评价。

近年较多研究揭示了影像基因组学在肿瘤影像学取得一些鼓舞人心的成果。Mazurowski等^[28]研究认为分子亚型为luminal B型的乳腺癌与肿瘤动态增强曲线特征之间有显著的相关性。Lee等^[29]通过研究214例Ⅰ 期非小细胞肺癌手术患者EGFR蛋白表达情况及其影像学特征之间的关系, 认为EGFR过表达阳性组肿瘤直径显著大于阴性组; 磨玻璃密度所占比例> 50%者在 EGFR 蛋白表达阴性组较阳性组更多见(P=0.007)。Yamamoto等^[30]研究172例非小细胞肺癌的基因突变与影像特征之间的关系, 认为影像特征(肿瘤的中心性生长、胸膜尾征、胸腔积液)与ALK基因突变之间有明显相关性。而Rizzo等^[31]研究认为伴胸膜凹陷征的周围型非小细胞癌的EGFR基因突变率为27.5%, 高于不伴胸膜凹陷征。前列腺癌发生是一个多基因参与、多因素作用及经过多阶段演变才最终形成复杂的生物现象。自2005年Tomlins等^[32]首次报导前列腺癌的跨膜丝氨酸蛋白酶编码基因TMPRSS2(transmembrane protease serine 2)与ETS转录因子家族成员ERG、ETVl等之间发生融合, 认为融合基因与前列腺癌的发生、淋巴结转移有关。所以可以利用影像信息与基因组关联研究肿瘤分子基因及其亚型的临床发病机制、表型的相关性, 用于指导临床治疗方案的制定。利用影像组学特征诠释疾病的分子基因表型特征, 将影像学与遗传学的不同优势有机结合, 通过定量化的技术, 将更多的疾病特征更加有效地描述出来, 弥补基因诊断技术在临床应用的局限性。