一种用于预测增殖型肝细胞癌的模型和诺模图的构建方法

1.本发明属于生物医学工程领域，尤其涉及将三维多频磁共振弹性成像用于预测增殖型肝细胞癌的模型和诺模图的构建方法。


背景技术：

2.肝细胞癌(hepatocellular carcinoma，hcc)占原发性肝癌的90％，是全球第三大癌症相关死亡原因。
3.hcc在分子和组织学水平上表现出高度异质的表型。通过整合形态学和分子改变，hcc可分为增殖型(～50％)或非增殖型(约50％)。每种增殖类型的特征在于激活与细胞增殖和存活相关的各种基因组途径(例如akt/mtor、met、tgf-b和胰岛素样生长因子途径)、染色体不稳定率高和异常表观遗传变化。与非增殖性肝癌相比，增殖性肝癌表现出侵袭性表型，具有中度到低度的细胞分化、频繁的血管侵袭、高肿瘤复发和不良预后。
4.因此，在治疗前检查中识别侵袭性hcc亚型可能具有很强的预后和治疗意义。
5.申请公布日为2022.06.28，申请公布号为cn 114678062 a的发明专利申请，公开了一种“基于多组学特征的肝细胞癌预后预测系统及其预测方法”，其预测系统包括数据输入模块、预测模块和显示模块；所述数据输入模块用于输入肝细胞癌患者数据并上传至所述预测模块，肝细胞癌患者数据包括ges评分水平、is评分水平和tnm分期；所述预测模块为列线图模型；所述显示模块用于显示肝细胞癌患者24和/或36个月总体生存概率；通过同时使用癌和癌旁组织，并整合使用两种组学数据以及肝细胞癌常用临床指标，通过列线图模型对肝细胞癌预后进行预测，可广泛用于肝细胞癌患者的临床预后评估。该预测系统通过同时利用来自癌组织和癌旁组织的多组学数据，包括基因表达量数据，免疫细胞相对浸润水平数据，提供了一个分析肝细胞癌患者预后生存的评估系统，相对而言，准确性更高，相较于单一组学数据构建的模型而言，该技术方案整合了多组学数据构建的模型，能更好地预测病人预后。但是由于其患者数据包括ges评分水平、is评分水平和tnm分期，且需要使用两种组学数据以及肝细胞癌常用临床指标，通过列线图模型对肝细胞癌预后进行预测，使得其计算过程较繁琐，建模过程涉及到的影响因素较多(模型中包含的预测变量越多，列线图的解释难度就越大)，故其实施临床预后评估的过程较为繁琐，在准确性评估中，未经外部评估(未经证实具有临床益处的“新”模型不一定能为日常决策提供帮助)，缺乏前瞻性，不利于在临床推广。
6.授权公告日为2022.05.13，授权公告号为cn 113571194 b的发明专利，公开了一种“肝细胞癌远期预后预测的建模方法及装置”，包括：获取肝细胞癌病例资料；对所述病例资料进行标准化处理，得到标准化的病例资料；对所述标准化的病例资料进行处理，得到肝细胞癌影像组学特征；对所述影像组学特征进行筛选，得到用于肝细胞癌远期预后预测的模型建构的特征；基于所述用于肝细胞癌远期预后预测的模型建构的特征建构肝细胞癌远期预后的预测模型。该技术方案根据肝细胞癌远期预后对术后影像特征进行筛选，将筛选出的特征建立预后模型，从而实现肝细胞癌远期预后的精准预测和提前进行治疗干预，对
肝细胞癌远期预后预测提供全流程影像数据的支撑。但是由于其所获得的肝细胞癌影像组学特征需要放射科医生进行人工勾画，非常耗时，且受不同放射科医生的个人主观因素影响大，勾画结果不稳定；其数据标准化程度不足(现有的数据集有各种标准，缺乏对疾病的统一理解)，模型泛化能力不足(由于数据样本量和样本多样性不足，数据标记的一致性较差，特征筛选方法较差，无法找到反映肿瘤异质性的稳定和通用特征)；此外，由于缺乏对数据和标记数据的统一和明确的标准化描述，导致机器学习和数据之间的交互存在障碍，以及机器对数据的真实含义的误解，使得其临床易用性不够理想(操作复杂、学习成本高、运算速度慢，而临床决策需要客观、严谨、高效)。
7.授权公告日为2021.08.10，授权公告号为cn 110197701 b的发明专利，公开了一种“新型多发性骨髓瘤诺模图构建方法”，包括：从geo数据库中获取多发性骨髓瘤患者的基因表达谱gse24080，并对所述基因表达谱gse24080中的基因进行预处理，获得表达值方差最大的前25％的5413个基因；对所述5413个基因进行wgcna基因共表达网络分析，以识别共表达的功能模块；通过pearson相关检验将所述功能模块与临床信息之间的相关性进行评估，以确定最显著的模块；所述最显著的模块包括：聚集于黑色模块中的240个基因，其分别与efs时间和os时间性状的相关性最强；使用cox比例危险度模型对最显著的模块中的基因进行单变量生存分析，并通过lasso回归筛选出由10个最佳基因组成的评分模型。该技术方案基于所述评分模型、血清β2m和高ldh的新型诺模图，来预测新诊断多发性骨髓瘤患者的3年总生存率和5年总生存率，有助于定量预测患者的生存率。但是其采用基因测试的模式，所涉及到的测试装置或测试工具的种类及数量较多，对于操作者或实施人员的技术要求较高，且整个基因预处理、测试过程中的工作量以及计算量较大，不能满足相关病房、科室实际诊断的工作需要。
8.授权公告日为2021.06.01，授权公告号为cn 112233112 b的发明专利，公开了“一种基于深度学习分析的肝细胞癌的ck19的预测方法”，包括将待检测图像分割为若干子图像；基于若干子图像以及经过训练的分割网络模型，确定所述待检测图像对应的特征图；基于所述特征图以及经过训练的分类模型，确定所述待检测图像对应的ck19预测类别。该技术方案通过采用分割网络模型确定特征图作为分类模型的输入项，并基于该特征图确定所述待检测图像对应的ck19预测类别，这样通过深度学习方式确定待检测图像对应的ck19预测类别，可以提高ck19预测类别的准确性。但是在实施过程中，其需要依赖放射科医生将待检测图像分割为若干子图像，进行肿瘤分割，分割过程缺乏统一的标准和识别标准，无法避免各个医生之间的个人识别差异，且需要使用复杂的机器学习过程，以及不同类型的医院数据之间的建模结果缺乏统一性和鲁棒性，无法得到简单稳定的定量参数，难以进行泛化和推广、普及；此外，基于机器学习的影像组学特征具有“黑盒子”的特点，没有准确、完整的公式与定义，其深度学习的过程或结果通常缺乏可解释性(或称缺乏生物学可解释性)，使得其学习结果不具有生物学可解释性和不易于被理解，这必然会阻碍机器学习在临床实践中的实现。
9.最近的研究表明，影像学检查结果可能与hcc的特定分子特征相关。
10.组织学上，增殖性肝细胞癌包括祖细胞、大小梁、硬化型、肉瘤样和中性粒细胞丰富的类型。祖细胞型hcc通常表现为环形强化(lr-m)，肝胆期低信号更明显，表观扩散系数较低，肿瘤边缘不光滑。有研究表明，实质性坏死、高血清afp水平和巴塞罗那临床癌症
(bclc)b或c期是大小梁(mtm)-hcc亚型的独立预测因子；实质性坏死有助于以65％的敏感性识别mtm-hcc。
11.近期radiology杂志上发表的一篇文章，其研究发现，边缘动脉期高强化(rim arterial-phase hyperenhancement，aphe)和高甲胎蛋白(α-fetoprotein，afp)(》100ng/ml)可以独立预测增殖型hcc。然而，26.2％的敏感性可能不足以进行临床诊断。此外，该样本没有考虑ck7阳性的hcc(另一种重要的增殖类型)，而且用的是传统的mri影像特征。
12.总的来说，对增殖性hcc的影像学特征的全面了解仍然有限。
13.磁共振弹性成像(magnetic resonance elastography，mre)允许无创定量体内组织的力学特性(或称为生物力学特征)，并为肿瘤生物学提供了新的见解。
14.mre已被开发用于检测和表征癌症，评估治疗反应，并研究与恶性转化相关的潜在生物物理机制。三维mre可用于非侵入性鉴别良性和恶性肝脏病变，并检测胰腺癌和前列腺癌。现在，三维mre对预测增殖性肝癌患者的价值尚不清楚。
15.目前，尚无基于肿瘤生物力学特性预测增殖型hcc的表型的报道。


技术实现要素：

16.本发明所要解决的技术问题是提供一种用于预测增殖型肝细胞癌的模型和诺模图的构建方法。其采用多因素逻辑回归分析，筛选与增殖型肝细胞癌hcc相关的变量；利用三维多频磁共振弹性成像(亦称3d多频mre)参数联合传统影像和临床特征，建立增殖型肝细胞癌hcc的预测模型和预测诺模图，基于3d多频mre的生物力学特征，来预测增殖型肝细胞癌hcc。其通过加入3d多频mre的生物力学特征参数，构建预测模型和预测诺模图，确定肿瘤生物力学特性预测增殖型肝细胞癌hcc的价值，通过整合3d多频mre的生物力学特征参数，提高接受者操作特性曲线(receiver operating characteristic curve，简称roc曲线)下面积area under curve(auc)，可以显著提高增殖型肝细胞癌hcc的诊断准确性。
17.本发明的技术方案是：提供了一种用于预测增殖型肝细胞癌hcc的模型和诺模图的构建方法，包括以下步骤：
18.步骤1：进行三维多频磁共振弹性成像扫描，即3d多频mre扫描，以及图像处理，获得增殖型肝细胞癌的生物力学特征；
19.步骤2：采用多因素逻辑回归分析筛选与增殖型肝细胞癌相关的变量，包括：
20.对传统影像特征、3d多频mre参数、临床特征，采用多因素逻辑回归分析的方法，筛选与增殖型肝细胞癌相关的变量，确定增殖型肝细胞癌的预测因素；
21.步骤3，用筛选出的传统影像和临床特征、传统影像和临床特征基础上联合3d多频mre参数，分别建立增殖型肝细胞癌的预测模型和诺模图，并对模型进行外部验证。
22.具体的，所述的生物力学特征包括：肿瘤剪切波速c和肿瘤复剪切模量损失角以及肝脏剪切波速c和肝脏复剪切模量损失角
23.所述的多频包括30hz、40hz、50hz和60hz。
24.具体的，步骤1中所述的图像处理包括：
25.1)对3d多频mre获得的图像，使用bioqic软件进行处理，生成剪切波速c和复剪切模量损失角的视场图；
26.2)在剪切波速c和复剪切模量损失角的视场图上，绘制感兴趣区域roi；
27.3)测量剪切波速c和复剪切模量损失角值。
28.进一步的，在绘制增殖型肝细胞癌病灶的感兴趣区域roi时，选择一个显示最大横断面延伸的原发病变部位的主层面和两个相邻的层面来确定感兴趣区域roi，以获得所有肿瘤的生物力学特征；
29.对于肝脏，图像分析人员在绘制感兴趣区域roi时，只包括肝实质，同时避免肝脏边缘、大血管和伪影。
30.具体的，在步骤2中进行所述传统影像特征的评估时，由两位图像分析人员对于每个病变部位，独立评估每个增殖型肝细胞癌的以下传统影像特征：
31.(a)环形动脉期高强化；
32.(b)非周边廓清；
33.(c)包膜；
34.(d)瘤周增强；
35.(e)肿瘤边缘。
36.进一步的的，步骤2中所述的临床特征包括：性别、年龄、身体质量指数、肝硬化、病毒性肝炎，以及术前实验室结果；
37.其所述术前实验室结果至少包括：甲胎蛋白afp、癌胚抗原、糖类抗原125、糖类抗原199。
38.具体的，步骤3中所述的建立增殖型肝细胞癌的预测模型和诺模图，包括以下步骤：
39.3.1)用筛选出来的传统影像特征和临床特征建立第一预测模型：model 1；用筛选出来的3d多频mre生物力学特征—mre参数，联合筛选出来的传统影像特征和临床特征建立第二预测模型：model 2；
40.3.2)评估模型的预测性能，进行接受者操作特性曲线roc分析；使用roc曲线下面积auc、敏感性、特异性和精确性指标来评估预测增殖型hcc的诊断性能；
41.3.3)诺模图建立：基于增殖型hcc的预测模型建立诺模图：用筛选出来的传统影像特征和临床特征建立model 1的诺模图；在筛选出来的传统影像特征和临床特征的基础上加入筛选出来的3d多频mre的生物力学特征，建立model 2的诺模图，以用于展示model2比model 1的诊断性能提高；
42.其中筛选出来的的传统影像特征包括：环形动脉期高强化、瘤周增强和肿瘤边缘；筛选出来的临床特征包括：肝硬化、病毒性肝炎；
43.筛选出来的3d多频mre的生物力学特征—mre参数：包括肿瘤c和肿瘤
44.3.4)采用外部验证对模型进行验证，以验证训练模型的稳定性：验证模型结果和训练模型一样，显示整合肿瘤c和肿瘤参数指标可以显著提高术前对增殖型hcc的诊断性能；
45.3.5)采用delong检验，比较第一预测模型model 1和第二预测模型model 2预测增殖型肝细胞癌的auc。
46.进一步的，筛选出的传统影像特征包括：环形动脉期高强化、瘤周增强以及肿瘤边缘不光滑；
47.筛选出的临床特征包括：肝硬化、病毒性肝炎；
48.筛选出的3d多频mre生物力学特征—mre参数包括：肿瘤剪切波速c和肿瘤复剪切模量损失角
49.进一步的，在训练组和验证组中，第二预测模型model 2的auc均大于第一预测模型model 1的auc。
50.本发明的技术方案，还提供了一种计算机可读存储介质，其所述的计算机可读存储介质存储有一个或者多个程序，所述的一个或者多个程序可被一个或者多个处理器执行，以实现如上述任一项所述用于预测增殖型肝细胞癌的模型和诺模图的构建方法的步骤。
51.本发明的技术方案，又提供了一种终端设备，其所述的终端设备至少包括：处理器、存储器及通信总线；
52.所述存储器上存储有可被所述处理器执行的计算机可读程序；
53.所述通信总线实现处理器和存储器之间的连接通信；
54.所述处理器执行所述计算机可读程序，来实现上述任一项所述用于预测增殖型肝细胞癌的模型和诺模图的构建方法的步骤。
55.与现有技术比较，本发明的优点是：
56.1.本发明的技术方案，采用多因素逻辑回归分析，筛选与增殖型hcc相关的变量；利用3d多频mre参数联合传统影像和临床特征，建立增殖型hcc的预测模型和预测诺模图，基于3d多频mre的生物力学特征，来预测增殖型hcc，其所建立模型的稳定性好；
57.2.本发明的技术方案，在传统影像和临床特征的技术上，加入3d多频mre参数的方法构建增殖型hcc的预测模型和诺模图，确定肿瘤生物力学特性预测增殖型hcc的价值，其通过整合3d多频mre参数，可以显著提高增殖型hcc的诊断准确性；
58.3.本发明的技术方案，首先参考传统影像和临床特征，然后加入3d多频mre参数，建立增殖型hcc的预测模型和诺模图，提高了auc，进而达到提高诊断性能的效果。
59.4.对于3d多频mre的感兴趣区域(roi)，由两个放射科医生进行人工标识，人工标识比机器学习更加可靠；两名放射科医生的分别单独测量，使得测量结果更加客观和稳定，避免了个人主观因素的影响。
60.5.本发明的技术方案，在传统影像和临床特征的技术上，整合了3d多频mre参数，构成了一种成熟的、高诊断性能的定量方法；
61.6.本发明的技术方案，通过3d多频mre，提供了肿瘤内部力学特性的信息，而且所获得的肿瘤生物力学特性指标稳定、客观；
62.7.本发明的技术方案中所使用的是3d多频mre，其多频包括30、40、50和60hz，有助于更加全面、真实地获得肿瘤的生物力学特征。
附图说明
63.图1a是本发明用于增殖型hcc的预测模型和诺模图的构建方法流程示意图；
64.图1b本发明3d多频mre图像处理步骤的方法流程示意图；
65.图2是本发明用于增殖型hcc的预测模型和诺模图的技术路线示意图；
66.图3a和图3b是训练组和验证组病人的入选流程图；
67.图4是筛选出来的变量的森林图；
68.图5a是本发明预测增殖型hcc的训练组的两种模型roc曲线的示意图；
69.图5b是本发明预测增殖型hcc的外部验证组的两种模型roc曲线的示意图；
70.图6a和图6b是本发明预测增殖型hcc的两种模型的诺模图。
具体实施方式
71.下面结合附图和实施例对本发明做进一步说明。
72.基于三维多频磁共振弹性成像(亦称3d多频mre)技术的断层弹性成像，使用多频率数据采集和基于波数的反演算法。它分别产生剪切波速(c，m/s)和复剪切模量损失角(rad)的高分辨率定量图，作为组织硬度和流动性的替代物(亦称为生物力学特征)。
73.肝脏中的力学变化(即生物力学特征的变化)，可能是导致肿瘤发展的遗传不稳定性的来源，这可以用于从生物力学角度探讨肝脏肿瘤的生物学机制。
74.由于增殖型肝细胞癌具有潜在的组织学特征，如丰富的纤维间质、丰富的坏死成分和少量的肿瘤假性包膜，假设肝脏肿瘤组织的生物力学特征对结构组成和排列敏感，并且增殖型和非增殖型肝细胞癌之间可能存在差异。
75.本发明的技术方案，旨在利用3d多频mre参数来识别增殖型hcc的生物力学特征，并开发了基于3d多频mre的预测诺模图，以此来显著提高术前对增殖型hcc的诊断性能。
76.如图1a、图1b以及图2中所示，本发明提供的一种用于预测增殖型hcc的模型和诺模图的构建方法，具体包括以下步骤：
77.步骤1，进行3d多频mre扫描及图像处理，获得肿瘤的生物力学特征；
78.其图像处理包括对3d多频mre获得的图像，使用bioqic软件进行处理，生成了剪切波速(shear wave speed，c)和复剪切模量损失角(loss angle，)的视场图(field-of-view maps)：剪切波速c和复剪切模量损失角的视场图，在剪切波速c和复剪切模量损失角的视场图上，绘制感兴趣区域(region of interest，roi)以及测量剪切波速c和复剪切模量损失角值：
79.其所述的生物力学特征包括：肿瘤剪切波速c和肿瘤复剪切模量损失角以及肝脏剪切波速c和肝脏复剪切模量损失角
80.在实施例中，从2020年7月至2021年11月，本技术方案连续入选了273例局灶性肝脏病变患者。其中152例被排除的原因如下：非hcc病变(n＝86)、缺乏病理结果(n＝45)，次优图像质量(n＝15)、肝外原发性恶性肿瘤史(n＝3)和局部或全身治疗史(n＝3)。最后，纳入121例患者，共纳入124例hcc病变，作为训练组。符合标准的2021年12月至2022年9月的33例病人被纳入外部验证组。
81.图3a和图3b中所示，给出了训练组和外部验证组病人的入选流程示意图。
82.步骤2，采用多因素逻辑回归分析筛选与增殖型hcc相关的变量；
83.对传统影像特征、3d多频mre参数、临床特征，采用多因素逻辑回归分析的方法，筛选与增殖型肝细胞癌相关的变量，确定增殖型肝细胞癌的预测因素；
84.其中的传统影像特征：由两位图像分析人员对于每个病变部位，独立评估每个hcc的以下传统影像特征：
85.(a)环形动脉期高强化；
86.(b)非周边廓清；
87.(c)包膜；
88.(d)瘤周增强；
89.(e)肿瘤边缘。
90.其中的临床特征包括：性别、年龄、身体质量指数、肝硬化、病毒性肝炎和术前实验室结果(afp(甲胎蛋白)、cea(癌胚抗原，一种肿瘤标志物)、ca125(糖类抗原125，一种糖蛋白性肿瘤相关抗原)、ca199(糖类抗原199)等)。
91.步骤3，用筛选出的传统影像和临床特征，传统影像和临床特征基础上联合3d多频mre参数，分别建立增殖型hcc的预测模型和诺模图，并对模型进行外部验证。
92.其筛选出的传统影像特征包括：环形动脉期高强化、瘤周增强和肿瘤边缘。
93.其筛选出的临床特征包括：病人信息、包括性别、年龄、身体质量指数、肝硬化、病毒性肝炎。
94.筛选出的3d多频mre的生物力学特征—mre参数包括肿瘤c和肿瘤
95.进一步，所述步骤1具体包括：
96.3d多频mre图像扫描：患者接受3d多频mre检查，用以获得患者的断层弹性成像图像。
97.具体的，3d多频mre图像扫描包括：
98.3d多频mre扫描在1.5t mri扫描仪(magnetom aera,siemens,erlangen,germany)上进行。
99.使用四个由压缩空气驱动的基于表面的压力垫，产生了30、40、50和60hz的机械振动，并将其传递到肝脏。
100.两个后部气动致动器(0.6bar)和一个前部气动致动器(0.4bar)放置在肝脏附近。
101.采用流量补偿运动编码梯度(flow-compensated motion-coding gradients，megs)的单次自旋回波平面成像(epi)序列，获得三维波场。
102.在被测试者自由呼吸时获得了15个连续的横向切片，fov(field of view)为312
×
384mm2，矩阵大小为104
×
128，分辨率为3
×3×
5mm3。
103.进一步的成像参数为重复时间(repetition time，tr)＝2050ms；回波时间(echo time，te)＝59ms；grappa因子为2的并行成像；meg振幅为30mt/m；30、40和50hz振动频率为43.48hz，60hz振动频率为44.88hz。
104.一组完整的多频mre数据的总采集时间约为3.5分钟。
105.具体的，本发明技术方案中对3d多频mre的图像处理包括：
106.对3d多频mre获得的图像，使用https://bioqic-apps.com提供的专用软件bioqic进行处理，生成了剪切波速(shear wave speed，c)和复剪切模量损失角(loss angle，)的视场图(field-of-view maps)：c和的视场图。
107.绘制感兴趣区及测量c和值：在c和的视场图上(使用常规t2w mr图像进行肿瘤的解剖定位)，用image j软件手动绘制感兴趣区域(region of interest，roi)。
108.由于c与storage modulus(复剪切模量的实部)的平方根成正比，而从0(纯固体性质)连续变化为π/2(纯流体性质)，这两个参数也分别被认为是硬度(stiffness)和组织流动性(tissue fluidity)的替代品。
109.对于hcc病灶，在绘制感兴趣区域roi时，选择一个显示最大横断面延伸的原发病
变的主层面和两个相邻层面来确定roi。
110.对于肝脏，图像分析人员绘制感兴趣区域(roi)，只包括肝实质，同时避免肝脏边缘、大血管和伪影。
111.mre的结果由两名放射科医生(代号为图像分析人员1和图像分析人员2)独立测量。最后获得所有肿瘤的生物力学特征(肿瘤c和肿瘤以及肝脏c和肝脏)。
112.传统影像扫描：所有患者均接受常规多参数肝脏mri检查(siemens healthcare；philips；united imaging)，配备了专用的18通道系统。成像方案包括t1加权图像(t1wi)、t2加权图像(t2wi)、使用三个b值(0、50和800s/mm2)的弥散加权成像(dwi)和多期动态对比增强(dce)成像(gd-dtpa)。
113.其中传统影像特征的评估包括：两位经验丰富的放射科医生(图像分析人员1和图像分析人员2)对所有术前mri特征进行了一致的审查。
114.对于每个病变部位，每个图像分析人员独立评估每个hcc的以下影像学特征：
115.(a)环形动脉期高强化：定义为边缘样强化，在动脉期具有相对低的中心区域；
116.(b)非周边廓清：根据2018年li-rads版本，定义为相对于复合肝实质，非外周视觉评估的整体或部分增强减少；
117.(c)包膜：根据2018年li-rads版本，定义为门静脉或延迟期平滑过度增强的线性、薄且增强的外周边缘；
118.(d)瘤周增强：定义为在肿瘤边缘外可检测到的多边形或新月形强化，在动脉期与肿瘤边界广泛接触，在平衡期或门静脉期与背景肝实质具有相同的强度；
119.(e)肿瘤边缘：在平衡期或门静脉期定义的肿瘤边缘可分为：平滑肿瘤边缘，表现为具有平滑轮廓的结节状肿瘤；非平滑肿瘤边缘，表现为非平滑结节性肿瘤，伴有局灶性结节外生长。
120.如果无法达成共识，则咨询第三位图像分析人员(具有多年经验的图像分析人员)进行最终决定。
121.其临床特征包括：(1)病人信息：通过查阅医院病史系统，搜集病人基本信息：性别、年龄、身体质量指数、肝硬化、病毒性肝炎；(2)术前实验室结果：afp(甲胎蛋白)、cea(癌胚抗原，一种肿瘤标志物、ca125(糖类抗原125，一种糖蛋白性肿瘤相关抗原)、ca199(糖类抗原199)等。
122.进一步的，所述步骤2具体包括：
123.对传统影像、临床特征、3d多频mre参数，采用统计方法—多因素逻辑回归分析的方法，筛选与增殖型hcc相关的变量，来确定增殖型hcc的预测因素。
124.多因素逻辑回归分析：逻辑回归分析属于概率型非线性回归，它是研究二分类(甚至多分类)观察结果与一些影响因素之间关系的一种多变量分析方法。多因素逻辑回归分析是分析多个因变量与一个自变量之间的逻辑关系。
125.在临床研究中，各种指标或因素不是相互独立的，而是相互作用和影响的。如果需要考虑多个因素对同一结果的影响，则需要进行多因素逻辑回归分析。
126.多因素逻辑回归分析的用途：
127.1)多因素逻辑回归可以识别和控制混杂因素；
128.2)多因素逻辑回归可用于分析影响因素，探讨自变量是否为独立影响因素；
129.3)多因素逻辑回归可用于预测分析，以探讨自变量对结果的预测作用。
130.多因素分析中包含的变量，是在单因素逻辑回归分析中与增殖型hcc显著相关的变量，以及被认为影响临床结果的相关变量，具体的相关变量参见表1中所示。
131.表1与增殖型hcc相关的变量的单变量和多变量分析
132.[0133][0134]
最后筛选出与增殖型hcc相关的独立预测因素：3d多频mre的生物力学特征—mre参数：肿瘤c、肿瘤临床特征：肝硬化、病毒性肝炎；传统影像特征：环形动脉期高强化、瘤周增强、肿瘤边缘不光滑。
[0135]
图4是采用本发明技术方案筛选出来的与增殖型hcc相关的变量的森林图。
[0136]
进一步，所述步骤3中，建立预测模型和诺模图，包括以下步骤：
[0137]
3.1)用筛选出来的传统影像特征和临床特征建立第一预测模型：model 1；用3d多频mre生物力学特征—mre参数，联合传统影像和临床特征建立第二预测模型：model 2；
[0138]
其中筛选出的传统影像特征包括：环形动脉期高强化、瘤周增强以及肿瘤边缘不光滑；
[0139]
筛选出的临床特征包括：肝硬化和病毒性肝炎；
[0140]
筛选出的3d多频mre参数包括：肿瘤c和肿瘤
[0141]
3.2)评估模型的预测性能，进行接受者操作特性曲线roc分析；使用roc曲线下面积auc、敏感性、特异性和精确性指标来评估预测增殖型hcc的诊断性能；
[0142]
3.3)诺模图建立：基于预测模型建立诺模图：用筛选出来的传统影像特征(环形动脉期高强化、瘤周增强、肿瘤边缘不光滑)和筛选出来的临床特征(肝硬化、病毒性肝炎)建立model 1的诺模图；加入筛选出来的3d多频mre的生物力学特征(肿瘤c和肿瘤)，建立model 2的诺模图，以用于展示model 2比model 1的诊断性能提高；
[0143]
3.4)采用外部验证对模型进行验证，以验证训练模型的稳定性：验证模型结果和训练模型一样，显示整合mre参数(肿瘤c和肿瘤)可以显著提高术前对增殖型hcc的诊断性
能；
[0144]
3.5)采用delong检验，比较第一预测模型model 1和第二预测模型model 2预测增殖型hcc的auc。
[0145]
在训练组和验证组中，model 2的auc均大于model 1(均为p《0.05)。
[0146]
图5a为第一预测模型model 1的roc：结合肝硬化、病毒性肝炎、环形动脉期高强化、瘤周增强、肿瘤边缘不光滑；
[0147]
图5b为第二预测模型model 2的roc：结合肿瘤c、肿瘤肝硬化、病毒性肝炎、环形动脉期高强化、瘤周增强、肿瘤边缘不光滑。
[0148]
本发明的技术方案，通过整合3d多频mre参数，能够显著提高增殖型hcc的诊断性能(auc从0.72增加到0.81，p＝0.01)，提高诊断准确性。
[0149]
本发明的技术方案，还进行了外部验证，通过外部验证也有同样的发现(auc从0.62提高到0.77，p＝0.021)，证明所建立的模型的稳定性。
[0150]
使用roc曲线下面积auc、敏感性、特异性、精确性来评估预测增殖型hcc的诊断性能，其比较结果参见表2中所示。
[0151]
表2预测增殖型hcc各模型的诊断性能的比较
[0152][0153]
从上表所示数据可知，在model 1中加入3d多频mre参数后，获得model 2，auc显著提高了。
[0154]
3.3)诺模图建立：基于预测模型建立诺模图。
[0155]
图6a所示为第一预测模型model 1的诺模图：结合肝硬化、病毒性肝炎、环形动脉期高强化、瘤周增强和肿瘤边缘不光滑得到的诺模图；
[0156]
图6b所示为第二预测模型model 2的诺模图：结合肿瘤c、肿瘤肝硬化、病毒性肝炎、环形动脉期高强化、瘤周增强和肿瘤边缘不光滑得到的诺模图。
[0157]
图中对于肝硬化、病毒性肝炎、环形动脉期高强化、瘤周增强和肿瘤边缘不光滑，“1”是指存在。
[0158]
由图6a和图6b所示可知，加入3d多频mre参数后，诺模图的c-index由0.72提高到0.81，即显著提高了模型的诊断性能。
[0159]
本发明的技术方案，还提供了一种计算机可读存储介质，其所述的计算机可读存储介质存储有一个或者多个程序，所述的一个或者多个程序可被一个或者多个处理器执行，以实现如上述任一项所述用于预测增殖型肝细胞癌的模型和诺模图的构建方法的步骤。
[0160]
本发明的技术方案，又提供了一种终端设备，其所述的终端设备至少包括：处理器、存储器及通信总线；
[0161]
所述存储器上存储有可被所述处理器执行的计算机可读程序；
[0162]
所述通信总线实现处理器和存储器之间的连接通信；
[0163]
所述处理器执行所述计算机可读程序，来实现上述任一项所述用于预测增殖型肝细胞癌的模型和诺模图的构建方法的步骤。
[0164]
本发明的技术方案，采用多因素逻辑回归分析，筛选与增殖型hcc相关的变量，利用3d多频mre参数联合传统影像和临床特征，建立增殖型hcc的预测模型和预测诺模图，基于3d多频mre的生物力学特征，来预测增殖型hcc；其通过加入3d多频mre参数的方法构建预测模型和预测诺模图，确定肿瘤生物力学特性预测增殖型hcc的价值，通过整合3d多频mre参数，可以显著提高增殖型hcc的诊断准确性；其所建立模型的稳定性好，提高了auc，进而达到提高诊断性能的效果。
[0165]
本发明可广泛用于增殖型hcc的临床诊断领域。

技术特征：
1.一种用于预测增殖型肝细胞癌的模型和诺模图的构建方法，其特征是包括以下步骤：步骤1：进行三维多频磁共振弹性成像扫描，即3d多频mre扫描，以及图像处理，获得增殖型肝细胞癌的生物力学特征；步骤2：采用多因素逻辑回归分析筛选与增殖型肝细胞癌相关的变量，包括：对传统影像特征、3d多频mre参数、临床特征，采用多因素逻辑回归分析的方法，筛选与增殖型肝细胞癌相关的变量，确定增殖型肝细胞癌的预测因素；步骤3，用筛选出的传统影像和临床特征、传统影像和临床特征基础上联合3d多频mre参数，分别建立增殖型肝细胞癌的预测模型和诺模图，并对模型进行外部验证。2.按照权利要求1所述的用于预测增殖型肝细胞癌的模型和诺模图的构建方法，其特征是步骤1中所述的生物力学特征包括：肿瘤剪切波速c和肿瘤复剪切模量损失角φ以及肝脏剪切波速c和肝脏复剪切模量损失角φ；所述的多频包括30hz、40hz、50hz和60hz。3.按照权利要求1所述的用于预测增殖型肝细胞癌的模型和诺模图的构建方法，其特征是步骤1中所述的图像处理包括：1)对3d多频mre获得的图像，使用bioqic软件进行处理，生成剪切波速c和复剪切模量损失角φ的视场图；2)在剪切波速c和复剪切模量损失角φ的视场图上，绘制感兴趣区域roi；3)测量剪切波速c和复剪切模量损失角φ值。4.按照权利要求3所述的用于预测增殖型肝细胞癌的模型和诺模图的构建方法，其特征是在绘制增殖型肝细胞癌病灶的感兴趣区域roi时，选择一个显示最大横断面延伸的原发病变部位的主层面和两个相邻的层面来确定感兴趣区域roi，以获得所有肿瘤的生物力学特征；对于肝脏，图像分析人员在绘制感兴趣区域roi时，只包括肝实质，同时避免肝脏边缘、大血管和伪影。5.按照权利要求1所述的用于预测增殖型肝细胞癌的模型和诺模图的构建方法，其特征是在步骤2中进行所述传统影像特征的评估时，由两位图像分析人员对于每个病变部位，独立评估每个增殖型肝细胞癌的以下传统影像特征：(a)环形动脉期高强化；(b)非周边廓清；(c)包膜；(d)瘤周增强；(e)肿瘤边缘。6.按照权利要求1所述的用于预测增殖型肝细胞癌的模型和诺模图的构建方法，其特征是步骤2中所述的临床特征包括：性别、年龄、身体质量指数、肝硬化、病毒性肝炎，以及术前实验室结果；其所述术前实验室结果至少包括：甲胎蛋白afp、癌胚抗原、糖类抗原125、糖类抗原199。7.按照权利要求1所述的用于预测增殖型肝细胞癌的模型和诺模图的构建方法，其特
征是步骤3中所述的建立增殖型肝细胞癌的预测模型和诺模图，包括以下步骤：3.1)用筛选出来的传统影像特征和临床特征建立第一预测模型：model 1；用筛选出来的3d多频mre生物力学特征—mre参数，联合筛选出来的传统影像和临床特征建立第二预测模型：model 2；3.2)评估模型的预测性能，进行接受者操作特性曲线roc分析；使用roc曲线下面积auc、敏感性、特异性和精确性指标来评估预测增殖型hcc的诊断性能；3.3)诺模图建立：基于增殖型肝细胞癌的预测模型建立诺模图：用筛选出来的传统影像特征和临床特征建立model 1的诺模图；在筛选出来的传统影像特征和临床特征的基础上加入筛选出来的3d多频mre的生物力学特征，建立model 2的诺模图，以用于展示model 2比model 1的诊断性能提高；3.5)采用delong检验，比较第一预测模型model 1和第二预测模型model 2预测增殖型肝细胞癌的auc。8.按照权利要求7所述的用于预测增殖型肝细胞癌的模型和诺模图的构建方法，其特征是所述筛选出来的传统影像特征包括：环形动脉期高强化、瘤周增强以及肿瘤边缘不光滑；筛选出的临床特征包括：肝硬化、病毒性肝炎；筛选出的3d多频mre生物力学特征—mre参数包括：肿瘤剪切波速c和肿瘤复剪切模量损失角φ。9.一种计算机可读存储介质，其特征是：所述的计算机可读存储介质存储有一个或者多个程序，所述的一个或者多个程序可被一个或者多个处理器执行，以实现如权利要求1至8之任一项所述用于预测增殖型肝细胞癌的模型和诺模图的构建方法的步骤。10.一种终端设备，其特征是：所述的终端设备至少包括：处理器、存储器及通信总线；所述存储器上存储有可被所述处理器执行的计算机可读程序；所述通信总线实现处理器和存储器之间的连接通信；所述处理器执行所述计算机可读程序，来实现如权利要求1至8之任一项所述用于预测增殖型肝细胞癌的模型和诺模图的构建方法的步骤。

技术总结
一种用于预测增殖型肝细胞癌的模型和诺模图的构建方法，属生物医学工程领域。包括进行三维多频磁共振弹性成像扫描及图像处理，获得增殖型肝细胞癌肿瘤的生物力学特征；采用多因素逻辑回归分析的方法，筛选与增殖型肝细胞癌相关的变量，确定增殖型肝细胞癌的预测因素；用筛选出的传统影像和临床特征、传统影像和临床特征基础上联合3D多频MRE参数，分别建立增殖型肝细胞癌的预测模型和诺模图，并对模型进行外部验证。其在传统影像和临床特征的技术上加入3D多频MRE参数的方法，来构建增殖型HCC的预测模型和诺模图；通过整合3D多频MRE参数，提高了AUC，构成了一种成熟、高诊断性能的定量方法，可以显著提高增殖型HCC的诊断准确性。性。性。


技术研发人员：严福华 李若坤 刘桂雪 周家豪
受保护的技术使用者：上海交通大学医学院附属瑞金医院
技术研发日：2023.04.25
技术公布日：2023/8/4