针对DR片的肺叶轮廓提取方法与流程

本发明涉及医疗信息化领域，特别涉及大规模筛查肺结核传染病的诊断辅助中肺叶轮廓提取方法。

背景技术：

DR片即数字化照影图像，是Digital Radiography的缩写。胸部DR片图像是传统的胸部透视成像的代替，其成像清晰度高而辐射低。X光成像是体检筛查肺部疾病的主要手段，这些肺部疾病包括肺部炎症、肿块、结核、肺癌等。胸透DR片利用人体组织的不同密度在X射线下成像差异，以观察厚度和密度差别较小部位的病变。但人体组织结构复杂，胸腔和腹腔包含了人体的关键器官，包含了高密度和低密度的各种脏器。因此其影像互相交叠，对观察和判断的影响尤其大。阅读和判别DR片对医生的经验要求很高，很难在基层医院和体检点展开大规模DR片检查。

在居民体检中，拍摄DR片的目的通常是做肺结核筛查。结核病是由结核分枝杆菌引起的，容易经空气中的飞沫甚至气溶胶传播。结核病患者主体是青壮年，这会导致家庭和社会劳动力损失。世界卫生组织指出结核病是全球重要的公共卫生问题。全球每年结核病夺取140万人生命。我国目前有500万活动性肺结核患者，每年有5万人死于结核病。根据世界卫生组织的估算，我国每年新发病人100万，年递降幅度为3％。在全球22个结核病高负担国家中，我国仅次于印度排在第二位。结核病是我国重点防控的主要传染病之一。中国疾病预防控制中心指出，我国是全球22个结核病高负担国家之一和全球27个耐多药肺结核高负担国家之一，耐多药肺结核患者人数位居全球首位；普通肺结核患者人数位居全球第二位，仅次于印度。

尽管我国肺结核疫情呈现感染人数多、发病人数多且现患人数多的严峻形势，但结核病本身是可防可治的。鉴于结核病的危害严重，防治工作难度大，地方各级结核病防治队伍规模仍较小，力量和经费尚不能适应防治需求，需要加强技术和资金投入，建立医防结合机制，形成切实有效的防治体系。当前，在结核病防治的实施上，存在着早期发现和治疗管理两个薄弱、困难环节。在此形势下，国内逐步开展利用基本公共卫生服务开展重点人群肺结核病人筛查项目，并对确诊病人全部予以免费抗结核治疗。

和肺结核相比，肺癌对病患身体健康的破坏严重程度胜于前者。众所周知，肺癌的死亡率远远高于其他癌症，而且近些年来还呈逐年增加的趋势。影像学检查是癌症诊断、检测和防治方面的重要技术手段之一。通常意义上的胸部影像学的主要研究对象就是肺部癌症，它一般通过建立肺窗来观察相应肺部影像。肺癌是与吸烟、大气污染和免疫功能低下有关的肿瘤性疾病。诸如近年来雾霾等因素对肺部反复的炎症刺激会带来慢性损害，影响支气管粘上皮的正常功能和机体的免疫抗病毒状态，对肺癌的发生有促进作用。未来可能导致肺癌高发。此外，肺癌的致癌潜伏期很长，一般需要10～30年。因此有些人年青时得过肺结核，治愈后肺部留下钙化灶，而在以后体检时即使发现肺部有阴影，也可能会被误认为是肺结核留下的钙化灶，从而可能错过了早期发现肺癌的机会。

与肺结核合并肺癌相比，肺癌并发肺结核或许更值得关注。统计显示，临床上约有10％～15％的肺癌患者会并发肺结核，尸检病例更高达30％左右。这是因为部分肺癌和肺结核在影像学上难以区别，临床上极易发生误诊、误治或延误治疗。临床上有的病人明明只是肺结核，但却被个别基层医院单凭一个影像学检查就诊断成了肺癌，甚至盲目地做起了放化疗，给病人的身体造成了严重的危害。相反，也有些病人明明是肺癌合并了肺结核，但却被漏诊了肺癌，仅仅将其当成单纯的肺结核进行治疗。造成这种现象的原因主要是因为结核病与肺癌的临床表现很相似，就是咳嗽、咳血等，而且两者通过影像学检查均可表现为胸部X线片上有阴影。

肺部DR片智能化处理的一大难点是如何确定肺叶的区域。若能准确确定肺叶范围，减弱或消除肺部以外的干扰，就能更有利于发现细微的病变。此外，肺叶轮廓的形状本身也是判断体检者相关生理指标的重要因素；可靠的肺叶轮廓提取算法可以使体检者减少复检的次数，拍摄更少的胸片，在较低辐射剂量的代价下获取医生对病灶确诊。

技术实现要素：

本发明所要解决的技术问题是提供一种针对DR片的肺叶轮廓提取方法，此方法无需人工干预，全自动化的处理各种胸透DR片，提取肺叶轮廓。

为解决上述技术问题，本发明采用的技术方案是：

一种针对DR片的肺叶轮廓提取方法，包括以下步骤：

步骤1：通过离线训练获取肺叶轮廓的代表性模板，包括左肺上部分模板、左肺下部分模板、右肺上部分模板以及右肺下部分模板；

步骤2：胸部DR片肺叶区域提取系统初始化，从DR片数据库中读取一幅DICOM图像；

步骤3：根据DICOM图像的尺寸，对图像做三层金字塔分解；进行降采样过程，将DICOM图像在高斯平滑的基础上不断缩小为原来图像高度的1/2，得到一系列图像，即为高斯图像金字塔；进行升采样过程，将每个图放大为原来高度的2倍，并与高斯金字塔下原始图像求差，得到拉普拉斯图像金字塔；

步骤4：用Gabor滤波器组重建待处理图像，将Gabor滤波后重建图像的残差，用自适应的局部二值化算法，将其转换为黑白图像；

步骤5：用Zhan-Suen细化算法将黑白图像细化，得到DICOM图像中包含其他干扰的关键轮廓图像；

步骤6：调用离线训练的左肺上部分模板、左肺下部分模板、右肺上部分模板以及右肺下部分模板，对每一个模板，用其作为卷积核算子，对步骤5所述的轮廓图像进行卷积，记录每一处的卷积计算结果；

步骤7：从卷积结果中筛选最佳可能性的局部最优卷积值，并组合评价；

步骤8：根据最佳匹配的上下模板及其最佳匹配位置组合生成肺叶轮廓形状。

根据上述方案，所述步骤1具体包括如下子步骤：

步骤1.1：采集具有代表性的DICOM图像，作为样板数据集；

步骤1.2：通过标注软件，人工手动绘制与样板相对应的黑白图像，其中的白色区域对应肺叶区域，而黑色区域为其他区域；保存黑白图像到文件，作为数据集中的样板肺叶的轮廓标注；

步骤1.3：遍历数据集中所有标注图像，提取每张轮廓标注图像中的轮廓，根据轮廓的重心，区分左肺轮廓和右肺轮廓，分别将左肺轮廓和右肺轮廓的像素按逆时针顺序加入左肺轮廓点集和右肺轮廓点集；

步骤1.4：遍历左肺轮廓点集，对每个左肺的轮廓，读取其对应的像素点，判别x和y坐标的最大值x_max，y_max和最小值x_min，y_min，并根据y_max和y_min提取肺叶上部分的区域和下部分区域，其中，上部分由参数x_max，x_min，y_min共同决定，下部分由参数x_max，x_min，y_max共同决定；

步骤1.5：用K-means算法对左肺叶上部分、左肺叶下部分轮廓、右肺叶上部分的轮廓以及右肺叶下部分的轮廓进行聚类，获取代表性模板。

根据上述方案，所述步骤7中组合评价的评价指标为：上部分模板匹配值、下部分模板匹配值和上下模板连通性匹配值的加权求和。

根据上述方案，所述上下模板连通性匹配值是通过上下模板端点及其切线构建三阶贝塞尔曲线的拐点个数来评价。

与现有技术相比，本发明的有益效果是：

1)通过提取肺部轮廓，使后续的病灶阴影提取模块更具有针对性，提升了后续处理的准确性。将处理后的肺结核肺部DR片给医生检验，可对肺叶区域做伪彩色提示，从而能够降低医生的视觉负荷，提高整体的识别准确率和处理效率；能够降低由于医生经验不足和经验差异对病情判别的影响，使分析结论更加客观和稳定。

2)适应各种胸透DR片，适应不同身形和年龄被拍摄者，实现全自动处理。

3)有效利用网络资源，实现远程会诊的功能，提高疑难病症的会诊可靠性。

4)集成了现有医疗设备和信息化网络化资源，无需购置额外专用设备，工作模式完全兼容传统方式，使迁移工作容易接受。同时提高了降低了设备的使用率，避免了设备闲置和资源浪费。

附图说明

图1为本发明中设备连接示意图。

图2为本发明中DR片的原始内容。

图3为本发明中DR片在Gabor滤波重建后的残差内容。

图4为本发明中滤波重建后DR片的细化效果。

图5为本发明中胸透DR片的肺部轮廓自动化提取方法的示意图。

图6为本发明中胸透DR片的肺部轮廓自动化提取方法的处理流程示意图。

具体实施方式

下面结合附图和具体实施方式对本发明作进一步详细的说明。本发明的功能包括：读取数字化胸透影像DICOM数据(Digital Imaging and Communications in Medicine，医学数字成像和通信)，利用数字图像处理技术检测肺部轮廓，并确定其区域范围；使标注内容为易于医生判别肺部病灶阴影的肺部区域影像，处理结果由医生确认以保证可靠性，同时医生反馈信息可在线改进系统性能。本发明方法基本思想为：面向胸透DR片的肺部肋骨自动化抑制方法方法，应用于基础规模化体检点及重大传染病筛查，主要由计算机自动化处理功能组成。设备通过读取数字化胸片透视DICOM图像数据，再将数据传输给肺部轮廓提取模块，包括以下检测步骤：

1)体检点工作人员将装配有DR片管理模块的计算机连接至DR片数据库，并配置DICOM图像文件读取的参数。

2)将体检点计算机连接至远程的医疗图像数据服务器。

3)在完成体检后，体检点计算机在空余时间将当日新增DR片自动上传到远程服务器。

4)系统服务器接收新增DR片，并给予其标签，并将其标签加如待处理队列，处理的优先级由时间先后确定。

5)系统服务器在扫描到未处理队列中存在数据时，自动运行肺叶轮廓自动化提取模块，并将提取结果保存到文件。

6)系统服务器在批量化处理一定数量DR片之后，生成包含DR原始图像和轮廓提取结果的融合图像形的处理报告。

7)根据系统登记的医师记录，系统服务器将汇报发送给不同医生，以人工确认分割结果的可靠性。

8)医生使用智能终端打开汇报，对认可的肺叶轮廓提取接口，可点击确认按钮；对不认可的分割结果，可根据DR片的难易程度，选择人工处理或者延后处理。其中，人工处理一般是针对疑难病例，用手工方式以打开人工标注的程序界面标注出肺叶区域。而延后处理是将图像重新加入待处理队列，但降低其处理优先级，使系统延后处理。

9)系统同时支持交叉评估的方式，由多名医生共同判别分割可靠性较低的肺部轮廓。

10)在未对医生的用户终端设备服务的时段，系统会根据医生的反馈，自动调用自适应状态更新模块，以根据医生新手动标记的图像特征对现有处理模块的参数改进和深化训练。

11)系统服务在进行状态更新后，会对待处理列表中的剩余图像进行重新处理。由于系统参数的优化，原来不能处理的DR片在改进后的系统中可以正确提取肺野形状。

12)体检点DR片的某一批次处理结束之后，由系统服务器向体检点计算机发出处理结束的消息，并有体检点信息员接收处理结果。

在以上各步骤中，系统会以图形化的方式提示医生操作，通过计算机自动识别和动态学习，以减少医生需操作键盘和鼠标的频率，从而提高处理效率和改善用户体验，使枯燥的标注和验证工作变得容易让人接受。此外，系统采用BS架构，使医生只要有用户名和密码即可在任意连接互联网络的计算机上进行肺结核影像的标注和评估，使工作平台从局部化的专用网络扩展到广域通用网络。不仅利于医生的工作和协调，而且有利于地方卫士部门和疾控单位对基层工作的把握及数据分析和挖掘。

本发明是一种针对胸透影像并利用计算机视觉技术进行肺部轮廓提取的方法，以DICOM数据为处理对象，以现有医学影像设备和计算机及互联网络为基础，不涉及专门设计硬件。

本发明面向居民健康体检的胸透DR片图像处理，提供辅助医生开展肺结核等关键肺部疾病筛查工作的关键核心算法。采用计算机图像处理技术自动获取DICOM图像中可能的肺轮廓边界的后续区域，从高可能性候选区域中进一步筛选最佳匹配方案，以解决目前大规模居民体检产生数据量太大，而医生在有限时间逐个人工标注且难以保持高精度检测的难题。

本发明利用医疗信息化的优势，能适应医疗人员的主观因素导致偏差、体检点的设备变更、操作人员的计算机水平差异等问题。整个处理过程简单方便，从基础上提高了肺结核筛查的处理效率，同时降低了体检点医护人员的工作负担，因而可使计算机辅助筛查及时推广到缺乏肺结核胸透影像评估经验的基层医疗机构，更有利于针对重点传染病的居民大规模体检进一步常态化和规范化。