[发明专利]尘肺病大体成像中病灶区域的提取与分析方法

说明书

技术领域

本发明涉及图像中的区域提取与分析方法，具体而言是尘肺病大体成像中病灶区域的提取与分析方法。

背景技术

尘肺病(pneumoconiosis)是由于在职业活动中长期吸入生产性粉尘并在肺内潴留而引起的以肺组织弥漫性纤维化为主的全身性疾病，我国目前尘肺病人数目之多、增长速度之快位居世界之首。

目前尘肺病诊断主要有X射线透射成像（胸片）诊断、病理学诊断两种形式。胸片诊断虽然是一种无损诊断方式，但由于成像方式的原因，其判读结果偏差较大，只适用于尘肺病的筛查，尘肺病的最终确诊则需要依靠病理学诊断，而尘肺病病理学诊断的关键在于病理大体图像分析。目前，病理大体图像分析完全依靠人工操作来完成，不仅操作繁琐、效率低下,而且计量误差大。因此提供一种效率高、误差小且无需人工交互的尘肺病大体成像中病灶区域的提取与分析方法十分必要。

发明内容

本发明的目的是提供一种效率高、误差小且无需人工交互的尘肺病大体成像中病灶区域的提取与分析方法。

为实现上述目的，本发明采用如下技术方案：尘肺病大体成像中病灶区域的提取与分析方法，其步骤如下：

S1.提取特征：

S11.将尘肺病大体成像数据库与计算机连接，

S12.分别计算出尘肺病大体成像数据库中已经标注的尘肺正常组织区域和病灶区域的中心坐标(c_x,c_y)，以区域中心为起点，每隔10°均匀画出36条射线，分别交于区域边界上的对应点，得到36条线段，将每条线段均匀分成k等分，将各条线段上相应的等分点分别相连，将区域分为k个环形部分，算出每个环形部分的灰度平均值，形成k维特征向量，对k维特征向量进行差分运算，得到k-1维特征向量，将k维特征向量和k-1维特征向量合并为一个2k-1维的特征向量F₁，

S13.利用尘肺病大体成像数据库中尘肺大体成像对卷积神经网络进行预训练和学习，利用尘肺病大体成像数据库中已经标注的病灶区域对经过预训练的卷积神经网络进行微调，算出特征向量F₂；

S2.训练分类器：

S21.通过步骤S12得到的F₁，训练一个用于判别一个区域是否为病灶区域的支持向量机分类器h₀，

S22.将步骤S13得到的F₂随机分割成k个大小相同的组，分别利用每一组特征训练一个用于判别一个病灶区域是否为尘斑、结节或纤维化区域的多分类支持向量机分类器h₁～h_k，h₁～h_k组成随机多分类支持向量机队列H；

S3.病灶区域提取：

S31.输入待提取分析的尘肺病大体成像，

S32.利用图论分割算法对步骤S31输入的尘肺病大体成像进行分割，得到区域D₁～D_n，

S33.用步骤S12的方法分别计算出区域D₁～D_n的特征E₁～E_n，利用分类器h₀判别是否为病灶区域，对非病灶区域标记为正常组织区域，对病灶区域进入下一步骤，

S34.用步骤S13的方法计算出特征f₂，利用分类器H中h₁～h_k通过加权投票进行判定，若分类器H的置信度大于设定阈值，则标记该区域病灶类型，否则进入下一步骤进行精细处理，

S35.对需要精细处理的区域的面积进行开平方运算，得到等效正方形边长d，根据等效正方形边长d扩大待精细处理区域面积，重复步骤S32～S34；

S4.病灶区域统计特征计算：

S41.根据提取的尘斑病灶，计算出尘斑病灶的面积，找出尘斑病灶面积最大值M₁，将M₁等分为s₁个区间，计算出不同区间的尘斑在图像中所占的比例，建立尺度比例特征向量代入公式进行积分求和，得到尺度比例积分向量

S42.根据提取的结节病灶，计算出结节病灶的面积，找出结节病灶面积最大值M₂，将M₂等分为s₂个区间，计算出不同区间的尘斑在图像中所占的比例，建立尺度比例特征向量代入公式进行积分求和，得到尺度比例积分向量