[发明专利]一种基于核分数和组织病理全玻片图像的三阴性乳腺癌肿瘤突变负荷预测方法及系统

权利要求书

1.一种基于核分数和组织病理全玻片图像的三阴性乳腺癌肿瘤突变负荷预测方法，其特征在于，包括：

S1：从乳腺癌组织病理图像中筛选出三阴性乳腺癌组织病理全玻片图像WSI；

S2：根据每个三阴性乳腺癌患者的基因突变情况，计算每个患者的肿瘤突变负荷值，并根据设定的阈值划分为肿瘤突变负荷高和低两组，分别记为TMB-H和TMB-L，作为每个患者的WSI对应的标签；

S3：将WSI切分为设定大小的块，并进行预处理；

S4：根据核分数函数从预处理后的块中筛选出核分数值满足阈值的块；

S5：搭建卷积神经网络分类模型，随机初始化分类模型的参数；

S6：对核分数值满足阈值的块进行颜色标准化，将进行颜色标准化后的块与对应的标签输入分类模型中训练出TMB分类器，其中，每一个块属于对应的WSI，块对应的标签为与块对应的WSI的标签；

S7：利用训练好的TMB分类器对三阴性乳腺癌肿瘤突变负荷进行预测。

2.如权利要求1所述的三阴性乳腺癌肿瘤突变负荷预测方法，其特征在于，步骤S2中根据每个三阴性乳腺癌患者的基因突变情况，计算每个患者的肿瘤突变负荷值，包括：将患者的体细胞蛋白编码区的非同义突变的肿瘤除以蛋白编码区的总长度得到每个患者的肿瘤突变负荷值，单位为mutations/mb，用以表征蛋白编码区的非同义突变分布的密度。

3.如权利要求1所述的三阴性乳腺癌肿瘤突变负荷预测方法，其特征在于，步骤S2中划分肿瘤突变负荷高和低两组时，采用中值划分法，阈值记为M，当患者的TMB值大于M时，该患者为TMB-H组，否则为TMB-L组。

4.如权利要求1所述的三阴性乳腺癌肿瘤突变负荷预测方法，其特征在于，步骤S3包括：

首先选定WSI的层数，基于该层依次保存设定大小的图像，以将图像切分成块；

从切分后的块中去除空白和不规则的块，其中，去除空白块的方法为：对每个patch进行像素均值计算，当该patch的像素均值小于设定阈值时，则保留该patch，否则丢弃；去除不规则块的方法为：计算每个patch的长和宽是否等于设定的patch大小，如果相等，则保留，否则丢弃。

5.如权利要1所述的三阴性乳腺癌肿瘤突变负荷预测方法，其特征在于，步骤S4包括：

S4.1：将RGB图像转为HED空间，提取出H通道数值；

S4.2：用H通道数值分别生成初步的掩码与用于清洗的掩码，其中，初步的掩码通过在H通道上通过多级图像阈值划分得到，用于清洗的掩码通过在H通道上进行多级图像阈值划分和形态学变换操作得到；

S4.3：将初步的掩码与用于清洗的掩码相减得到细胞核的掩码；

S4.4：计算每张patch的细胞核比率N_t，该比率为细胞核的掩码的非零像素个数与该掩码的像素总个数的比值；

S4.5：生成组织区域的掩码；

S4.6：计算组织的比率T_t，该比率为组织区域的掩码的非零像素个数与整个掩码像数总个数的比值；

S4.7：根据每张patch的细胞核比率、组织的比率T_t，通过核分数函数计算每个patch的核分数值s_t；

S4.8：根据得到的核分数值进行排序，筛选出核分数满足阈值的块。

6.如权利要5所述的三阴性乳腺癌肿瘤突变负荷预测方法，其特征在于，步骤S4.7中核分数函数为：

s_t＝N_t·tanh(T_t),0≤s_t＜1

其中，s_t代表第t个patch的细胞核的分数，N_t表示在patcht上的细胞核的比率，T_t表示patcht上组织的比率，patcht表示第t个patch，即第t个块。

7.如权利要1所述的三阴性乳腺癌肿瘤突变负荷预测方法，其特征在于，步骤S5中的分类模型使用resnet18作为特征提取模块，并将最后一层全连接层的输出修改为2。