尘肺病大体成像中病灶区域的提取与分析方法
【专利摘要】本发明涉及尘肺病大体成像中病灶区域的提取与分析方法，其步骤如下：算出一个2k-1维的特征向量F1，进行预训练并算出特征向量F2，训练一个用于判别一个区域是否为病灶区域的支持向量机分类器h0，训练一个用于判别一个病灶区域类别的多分类支持向量机分类器h1～hk，输入尘肺病大体成像，提取病灶区域特征，对病灶区域进行分类，病灶区域统计特征计算。本发明通过对图像像素进行空间拓展形成面元小区域，利用训练好的分类器对面元区域进行分类，并对分类置信度低的面元区域进行精细处理，从而实现不同病灶区域的准确提取。本发明效率高、误差小且无需人工交互。
【专利说明】尘肺病大体成像中病灶区域的提取与分析方法
【技术领域】
[0001]本发明涉及图像中的区域提取与分析方法，具体而言是尘肺病大体成像中病灶区域的提取与分析方法。
【背景技术】
[0002]尘肺病(pneumoconiosis)是由于在职业活动中长期吸入生产性粉尘并在肺内潴留而引起的以肺组织弥漫性纤维化为主的全身性疾病，我国目前尘肺病人数目之多、增长速度之快位居世界之首。
[0003]目前尘肺病诊断主要有X射线透射成像(胸片)诊断、病理学诊断两种形式。胸片诊断虽然是一种无损诊断方式，但由于成像方式的原因，其判读结果偏差较大，只适用于尘肺病的筛查，尘肺病的最终确诊则需要依靠病理学诊断，而尘肺病病理学诊断的关键在于病理大体图像分析。目前，病理大体图像分析完全依靠人工操作来完成，不仅操作繁琐、效率低下，而且计量误差大。因此提供一种效率高、误差小且无需人工交互的尘肺病大体成像中病灶区域的提取与分析方法十分必要。
【发明内容】
[0004]本发明的目的是提供一种效率高、误差小且无需人工交互的尘肺病大体成像中病灶区域的提取与分析方法。
[0005]为实现上述目的，本发明采用如下技术方案:尘肺病大体成像中病灶区域的提取与分析方法，其步骤如下:
[0006]S1.提取特征:
[0007]Sll.将尘肺病大体成像数据库与计算机连接，
[0008]S12.分别计算出尘肺病大体成像数据库中已经标注的尘肺正常组织区域和病灶区域的中心坐标(cx，cy)，以区域中心为起点，每隔10°均匀画出36条射线，分别交于区域边界上的对应点，得到36条线段，将每条线段均匀分成k等分，将各条线段上相应的等分点分别相连，将区域分为k个环形部分，算出每个环形部分的灰度平均值，形成k维特征向量，对k维特征向量进行差分运算，得到k-Ι维特征向量，将k维特征向量和k-Ι维特征向量合并为一个2k-l维的特征向量F1,
[0009]S13.利用尘肺病大体成像数据库中尘肺大体成像对卷积神经网络进行预训练和学习，利用尘肺病大体成像数据库中已经标注的病灶区域对经过预训练的卷积神经网络进行微调，算出特征向量F2;
[0010]S2.训练分类器:
[0011]S21.通过步骤S12得到的F1，训练一个用于判别一个区域是否为病灶区域的支持向量机分类器Iv
[0012]S22.将步骤S13得到的F2随机分割成k个大小相同的组，分别利用每一组特征训练一个用于判别一个病灶区域是否为尘斑、结节或纤维化区域的多分类支持向量机分类器hi~hk, Ii1~hk组成随机多分类支持向量机队列H ；
[0013]S3.病灶区域提取:
[0014]S31.输入待提取分析的尘肺病大体成像，
[0015]S32.利用图论分割算法对步骤S31输入的尘肺病大体成像进行分割，得到区域D1 ~Dn,
[0016]S33.用步骤S12的方法分别计算出区域D1~Dn的特征E1~En，利用分类器&判别是否为病灶区域，对非病灶区域标记为正常组织区域，对病灶区域进入下一步骤，
[0017]S34.用步骤S13的方法计算出特征f2，利用分类器H中Ii1~hk通过加权投票进行判定，若分类器H的置信度大于设定阈值，则标记该区域病灶类型，否则进入下一步骤进行精细处理，
[0018]S35.对需要精细处理的区域的面积进行开平方运算，得到等效正方形边长d，根据等效正方形边长d扩大待精细处理区域面积，重复步骤S32~S34 ；
[0019]S4.病灶区域统计特征计算:
[0020]S41.根据提取的尘斑病灶，计算出尘斑病灶的面积，找出尘斑病灶面积最大值M1,将M1等分为S1个区间，计算出不同区间的尘斑在图像中所占的比例，建立尺度比例特征向
量?，代入公式化进行积分求和，得到尺度比例积分向量&~F1,，
[0021]S42.根据提取的结节病灶，计算出结节病灶的面积，找出结节病灶面积最大值M2,将M2等分为S2个区间，计算出不同区间的尘斑在图像中所占的比例，建立尺度比例特征向
量/!μ ?，代入公式
【权利要求】
1.尘肺病大体成像中病灶区域的提取与分析方法，其步骤如下: S1.提取特征: S11.将尘肺病大体成像数据库与计算机连接， S12.分别计算出尘肺病大体成像数据库中已经标注的尘肺正常组织区域和病灶区域的中心坐标(cx，cy)，以区域中心为起点，每隔10°均匀画出36条射线，分别交于区域边界上的对应点，得到36条线段，将每条线段均匀分成k等分，将各条线段上相应的等分点分别相连，将区域分为k个环形部分，算出每个环形部分的灰度平均值，形成k维特征向量，对k维特征向量进行差分运算，得到k-Ι维特征向量，将k维特征向量和k-Ι维特征向量合并为一个2k-l维的特征向量F1, S13.利用尘肺病大体成像数据库中尘肺大体成像对卷积神经网络进行预训练和学习，利用尘肺病大体成像数据库中已经标注的病灶区域对经过预训练的卷积神经网络进行微调，算出特征向量F2; S2.训练分类器: S21.通过步骤S12得到的F1，训练一个用于判别一个区域是否为病灶区域的支持向量机分类器&， S22.将步骤S13得到的F2随机分割成k个大小相同的组，分别利用每一组特征训练一个用于判别一个病灶区域是否为尘斑、结节或纤维化区域的多分类支持向量机分类器Ii1~hk, hi~hk组成随机多分类支持向量机队列H ； S3.病灶区域提取: S31.输入待提取分析的尘肺病大体成像， S32.利用图论分割算法对步骤S31输入的尘肺病大体成像进行分割，得到区域D1~Dn, S33.用步骤S12的方法分别计算出区域D1~Dn的特征E1~En，利用分类器&判别是否为病灶区域，对非病灶区域标记为正常组织区域，对病灶区域进入下一步骤， S34.用步骤S13的方法计算出特征f2，利用分类器H中Ill~hk通过加权投票进行判定，若分类器H的置信度大于设定阈值，则标记该区域病灶类型，否则进入下一步骤进行精细处理， S35.对需要精细处理的区域的面积进行开平方运算，得到等效正方形边长d，根据等效正方形边长d扩大待精细处理区域面积，重复步骤S32~S34 ； S4.病灶区域统计特征计算: S41.根据提取的尘斑病灶，计算出尘斑病灶的面积，找出尘斑病灶面积最大值M1，将M1等分为S1个区间，计算出不同区间的尘斑在图像中所占的比例，建立尺度比例特征向量
2.根据权利要求1所述的尘肺病大体成像中病灶区域的提取与分析方法，其特征在于:所述步骤S35中，根据等效正方形边长d扩大2-4倍待精细处理区域面积，重复步骤S32 ~S34。
3.根据权利要求2所述的尘肺病大体成像中病灶区域的提取与分析方法，其特征在于:所述步骤S35中，根据等效正方形边长d扩大3倍待精细处理区域面积，重复步骤S32~S34。
4.根据权利要求1~3所述的尘肺病大体成像中病灶区域的提取与分析方法，其特征在于:所述步骤S35中 ，步骤S32~S34只重复一次。
【文档编号】G06K9/62GK103886318SQ201410125612
【公开日】2014年6月25日 申请日期:2014年3月31日 优先权日:2014年3月31日
【发明者】明德烈 申请人:武汉天仁影像科技有限公司