[发明专利]一种针对CT影像的肺结节恶性度分级系统

权利要求书

1.一种针对CT影像的肺结节恶性度分级系统，系统采用的方法：

步骤1：数据准备，具体为：

步骤1.1：肺内小结节影像提取，肺结节影像集合可表示为

G＝{G₁,G₂,...,G_N},i＝1,2,...,N，

其中，N为肺结节集合的总数，G_i为根据肺部CT影像病例集中的医生对结节中心点的标注信息，上下左右分别加g个像素得到包含结节的矩形区域，可表示为G_i＝{x+g,x-g,y+g,y-g}；

步骤1.2：划分肺结节影像集合，根据医生对每个结节标注的放射学语义征象和恶性度信息将肺结节影像集合G进行划分，其结构可表示为

I＝{I¹,I²,...,I^rsn},n＝1,2,...,rsn，

其中，I代表rsn种放射学语义征象集合，Iⁿ为第n个放射学语义征象的肺结节图像数据与征象评分的集合，Y_iⁿ为Iⁿ中肺结节影像G_i的评分，E为肺结节影像数据与恶性度评分的集合，y_i为肺结节影像G_i恶性度对应的评分；

步骤2：放射学语义征象的特征提取，将I中的每个征象集合从影像组学和CNN两个方面提取特征，具体为：

步骤2.1：影像组学特征提取，影像组学提取的特征集可表示为

其中，为第n个征象的影像组学特征集，的提取方式具体为：

步骤2.1.1：图像预处理，对Iⁿ中每个肺结节影像G_i通过自适应阈值以及形态组学膨胀操作得到分割后的影像G_i′；

步骤2.1.2：灰度特征提取，将每个肺结节影像G_i′通过灰度直方图提取肺结节影像的最大值、最小值、均值、方差、偏度、峰度、直方图熵、能量、标准差、中值、均方根11维统计特征，最终形成灰度特征集

步骤2.1.3：纹理特征提取，将每个肺结节影像G_i′采用灰度共生矩阵在(0，45，90，135)上提取纹理能量、最大概率、纹理熵、对比度、相关性、差异性、同质性、协方差、纹理均值和逆差分矩阵的均值40维特征，最终形成纹理特征集

步骤2.1.4：几何形状特征提取，将每个肺结节影像G_i′提取了粗糙度、直径、椭圆的长短轴比例、离心率、周长、面积、紧凑度，矩形度7维特征，最终形成几何特征集

步骤2.1.5：影像组学特征合并，将灰度、纹理、几何特征集进行组合形成58维影像组学特征集，可表示为

步骤2.2：CNN高阶特征提取，CNN提取的特征集可表示为

其中，为第n个征象的CNN特征集，的提取方式具体为：

步骤2.2.1：构建CNN模型，网络结构为两个卷积层，卷积核大小分别为32×3×3，64×3×3，两个池化层，池化层大小分别为32×2×2，64×2×2，一个全连接层其大小为1×200，一个dropout层，设置丢弃率为0.5，一个softmax层；

步骤2.2.2：训练CNN，对Iⁿ的影像集采用70％的数据进行训练，采用30％的数据进行测试，每次批量学习处理的影像个数为63，训练迭代次数为1000，使用adam优化器将学习率设置为0.001，模型的准确率达到85％以上；

步骤2.2.3：提取CNN特征，将droput层和softmax层去掉，将Iⁿ的影像集中的每个肺结节图像G_i作为输入通过该CNN的全连接层提取出200维特征，形成CNN特征集

步骤2.3：影像组学和CNN特征融合，融合后的特征集可表示为

其中，为第n个征象的融合特征集，其融合特征集的维度为258维，融合后的特征集可提高对肺结节信息的表达能力；

步骤2.4：特征选择，根据特征重要性排名通过阈值α对特征总集FI_total中的每个征象特征集进行选择，以减少特征集合中无用冗余的特征，选择后的特征集可表示为

其中，表示第n个征象经过特征选择后的特征，表示第n个征象中排名小于阈值α的特征；

步骤3：放射学语义征象的分类模型，rsn种放射学语义征象的分类模型可表示为

MI＝{MI¹,MI²,...,MI^rsn},n＝1,2,...,rsn，

其中，MIⁿ为第n个征象的分类模型，MIⁿ的训练方式具体为：

步骤3.1：数据准备，将特征集与Iⁿ中的标签Yⁿ相对应，生成第n个征象的数据集，可表示为

步骤3.2:初始化种群，种群由NP个随机生成的森林构成，可表示为

pop＝{RF₁,RF₂,....,RF_NP},j＝1,2,...,NP，

其中，NP为种群规模，RF_j为tn个CART分类树的组合，tn为可变参数；为了增加种群的多样性，每次有放回的随机采样40％的DIⁿ作为个体RF_j初始化的数据集，置当前迭代数值gen＝0；

步骤3.3：搜索运算，种群进化主要依靠交叉变异算子实现，按照交叉概率P_c和变异概率P_m选择本次搜索为交叉或者变异操作，并转到相应的搜索步骤，具体为：

步骤3.3.1：变异算子，变异运算定义为RF′＝⊙RF，其中RF为从种群pop^gen中随机挑选的个体，RF′为变异后的个体，⊙为变异运算符；由于个体的性能与每颗决策树的性能以及决策树之间组合的结构有很大关系，故采取2个层面的变异策略对个体进行变异；变异策略1为个体中树层面的变异，首先从个体RF中随机挑选两颗树，分别随机选取两棵树的交换点，再交换以交换点为根节点的子树，生成新的个体RF′；变异策略2为个体结构层面的变异，将从个体中增加树和删除树两个方面考虑，如果随机数小于0.5，则向RF中添加一颗随机生成的树，否则将从RF中删除一颗随机被选择的树，然后生成新的个体RF′；按照概率P_ms将2种不同的变异策略进行有效的结合，不仅对个体的组合结构进行搜索，还对个体中的树进行搜索；

步骤3.3.2：交叉算子，交叉运算可定义为其中RF₁和RF₂为从种群pop^gen中随机的挑选两个个体，RF₁′和RF₂′为交叉后的个体，为交叉运算符；由于个体中树的个数为可变大小，故交叉操作具体为：首先比较RF₁和RF₂中树的个数tn₁和tn₂，从树个数较多的RF₁中无放回的随机挑选出tn₂个树组成RF₃，按照概率P_cg将RF₃和RF₂中对应位置的树进行互换，生成新的个体RF₃′和RF₂′，将RF₃′中的树放回RF₁中得到RF₁′；

步骤3.4：选择阶段，首先对步骤3.3中产生的子代个体RF′与父代个体RF进个体进行评估，然后选择性能较好的个体遗传到下一代种群pop^gen+1中，具体为：

步骤3.4.1：个体评估，采用准确率作为个体评估方式，计算表达式如下

其中，TP表示实际为正被预测为正的样本数量，FP表示实际为负但被预测为正的样本数量，TN表示实际为负被预测为负的样本数量，FN表示实际为正但被预测为负的样本数量；

步骤3.4.2：个体选择，对子代个体RF′与父代个体RF进行竞争，准确率高的个体加入到pop^gen+1中；

步骤3.5：若下一代种群pop^gen+1中个体数量小于种群规模NP，则转步骤3.3，否则转步骤3.6；

步骤3.6：更新全局最优个体，若全局最优个体B_best为空时，则否则在B_best和中根据步骤3.4选择获胜的个体更新全局最优个体，为种群pop^gen+1中准确率最高的个体，置gen＝gen+1；

步骤3.7：判断当前迭代次数gen是否达到最大迭代数值NG，若满足输出最优的个体B_best作为第n个征象的分类器MIⁿ，否则转步骤3.3；

步骤4：恶性度分级识别模型，具体为：

步骤4.1：数据准备，将数据集E中的G_i图像集通过步骤2和步骤3得到放射学语义征象的评分集X_i，将X_i与E中的标签y_i相对应，生成恶性度分级的数据集，可表示为

X_i＝{x¹,x²,...,x^rsn}，

其中，X_i代表第i个肺结节图片rsn种不同放射学语义征象的评分，y_i为第i个肺结节图片评分集X_i对应的恶性度级数；

步骤4.2：恶性度分级模型，采用softmax模型进行分级，模型可表示为

其中,W＝{w₁,w₂,...,w_K}是一个[rsn,K]型矩阵，rsn为数据集ED特征维度，X_i为数据集中第i个样本，K为恶性度级别；建立softmax模型后，采用差分进化算法对模型的参数W进行优化，差分进化算法优化的目标是最大化样本恶性度分级准确的数量，目标函数可表示为

其中，F_i为样本X_i经过softmax函数分级的结果是否正确，为ED中trn个样本中分级结果准确的数量；在优化过程中将参数W转成rsn×K维的行向量，在评价阶段中将rsn×K维的行向量转成[rsn,K]型矩阵。