一种超声甲状腺结节自动定位方法技术

本发明专利技术公开了一种超声甲状腺结节自动定位方法，属于医学图像处理领域，包括：根据B超图像的水平方向平均灰度值的分布，得到甲状腺腺体层；在甲状腺腺体层中，利用甲状腺结节分类器，得到初始目标检测图；从B超图像获取相位不对称图，将初始目标检测图在相位不对称图中检测到的目标位置全部作为非目标，得到目标检测图；在目标检测图中，对各目标进行位置和尺寸约束，进而得到甲状腺结节定位结果。本发明专利技术能将与甲状腺结节内部具有相似纹理和强度信息的周围组织和结节很好地分离开来，从而准确定位到甲状腺结节区域。

【技术实现步骤摘要】
一种超声甲状腺结节自动定位方法
本专利技术属于医学图像处理领域，更具体地，涉及一种超声甲状腺结节自动定位方法。
技术介绍
研究表明，全世界范围内有越来越多的成年人出现甲状腺结节损伤病况，在过去30年，甲状腺癌症的发生率增加了2.4倍，已经对人类的健康造成了极大地威胁。目前，甲状腺结节的检查通常是医生通过B超进行的。B超灵活实时，但严重依赖于医生的经验与手法，比CT、MRI等其它影像更难以处理，所以三甲医院甲状腺结节的平均正确率也只有60％多，如果不进行组织切片检查，不同的医生常有不同的判断，因此，设计一套计算机辅助诊断系统来协助医生进行超声甲状腺结节的良恶性诊断是一件意义深远的事情。常见的超声甲状腺诊断系统中，往往需要有经验的医生手动给定一个包含完整甲状腺结节的感兴趣区域(ROI)，然后在该感兴趣区域内对结节进行良恶性识别。这种方法有两个缺点：(1)由于个人认知和经验不同，两个不同的医生给出的ROI可能不尽相同，这会对结节的识别造成不同程度的影响。(2)所有样本(包括训练和测试样本)都需要医生给定ROI，这无疑增加了医生的工作负担。由此可见，现有技术存在无法准确的定位到甲状腺结节区域的技术问题。
技术实现思路
针对现有技术的以上缺陷或改进需求，本专利技术提供了一种超声甲状腺结节自动定位方法，由此解决现有技术存在无法准确的定位到甲状腺结节区域的技术问题。为实现上述目的，按照本专利技术的一个方面，提供了一种超声甲状腺结节自动定位方法，包括：(1)根据B超图像的水平方向平均灰度值的分布，得到甲状腺腺体层；(2)在甲状腺腺体层中，利用甲状腺结节分类器，得到初始目标检测图；(3)从B超图像获取相位不对称图，将初始目标检测图在相位不对称图中检测到的目标位置全部作为非目标，得到目标检测图；(4)在目标检测图中，获取各目标的轮廓，计算各目标的中心位置和尺寸，对各目标按照中心位置离目标检测图的中心位置由近到远的顺序进行结节目标检测，当目标的尺寸大于阈值时，该目标为结节目标区域，其他目标为伪目标，进而得到甲状腺结节定位结果。进一步的，步骤(1)中的B超图像为对原始B超图像进行对比度增强处理后的图像。进一步的，步骤(1)的具体实现方式为：根据B超图像的水平方向平均灰度值的分布，得到水平方向最大平均灰度值的行索引值Rmax和从下到上遍历第一个水平方向最小平均灰度值的行索引值Rmin，进而得到Rmax与Rmin之差*90％的行索引值Rupper和Rmax与Rmin之差*25％的行索引值Rlower，取Rmax与Rupper之间的大者作为甲状腺腺体层的上边界，取Rmin与Rlower之间的小者作为腺体层的下边界，得到甲状腺腺体层。进一步的，甲状腺结节分类器的训练方法为：获取样本B超图像，根据样本B超图像的水平方向平均灰度值的分布，得到样本甲状腺腺体层，将样本B超图像中的甲状腺结节区域作为正样本，将样本甲状腺腺体层中除了甲状腺结节区域外的非结节区域作为负样本，利用正样本和负样本训练支持向量机，得到甲状腺结节分类器。进一步的，步骤(4)的具体实现方式为：在目标检测图中，获取各目标的轮廓，计算各目标的中心位置和尺寸，得到各目标中的最大尺寸值S，阈值为S/3，对各目标按照中心位置离目标检测图的中心位置由近到远的顺序进行排列，分别记录为O1，O2，O3，……On，先判断O1处的目标，如果其尺寸大于S/3，那么该目标为结节目标区域，其他目标为伪目标，否则，该目标作为伪目标，然后继续判断O2处的目标，如果其尺寸大于S/3，那么该目标为结节目标区域，其他目标为伪目标，否则，该目标为伪目标，然后继续判断O3处的目标，以此类推，进而得到甲状腺结节定位结果。进一步的，步骤(4)还包括：在目标检测图中，首先去除与甲状腺腺体层边界相连接的区域以及与B超图像左右边界相连接的区域，然后获取各目标的轮廓。总体而言，通过本专利技术所构思的以上技术方案与现有技术相比，能够取得下列有益效果：(1)由于在超声图像中，甲状腺的正常组织可能存在和结节内部相似的纹理和强度特征，因此，现有技术并没有将这种情况下的结节和周围组织完全分离开来，这样，就无法准确的定位到甲状腺结节的感兴趣区域。强度特征即为灰度值，本专利技术根据B超图像的水平方向平均灰度值的分布，得到甲状腺腺体层；在甲状腺腺体层中，得到初始目标检测图；从B超图像获取相位不对称图，将初始目标检测图与相位不对称图结合，得到目标检测图；在目标检测图中，对各目标进行位置和尺寸约束，进而得到甲状腺结节定位结果。本专利技术能将与甲状腺结节内部具有相似纹理和强度信息的周围组织和结节很好地分离开来，从而准确定位到甲状腺结节区域。(2)本专利技术首先根据B超图像的水平方向平均灰度值的分布，得到甲状腺腺体层；然后对甲状腺腺体层进行后续检测，有效减少了工作量，提高了定位甲状腺结节的准确率和效率。(3)由于甲状腺结节和一些非结节部分具有相似的纹理特征，所以初始目标检测图得到的目标并不仅仅是结节部分，还有很多其他的非结节部分。除此，有些B超图像中，结节和周围的非结节部分连接到了一块。因此，本专利技术从B超图像获取相位不对称图，将初始目标检测图在相位不对称图中检测到的目标位置全部作为非目标，由于相位不对称图其实是从频域位置分析图像灰度的突变位置，能够摆脱对图像亮度变化的依赖，所以能在边缘较模糊的B超图像中依然较好的检测到结节边界，因此通过将初始目标检测图与相位不对称图相结合，可以很好地分离结节与周围纹理，强度信息相似的组织。提高了定位甲状腺结节的准确率和效率。(4)本专利技术定位的甲状腺腺体层上下边界均比医生实际标注的甲状腺腺体层上下边界更靠外，这样是为了保证真实的腺体层被包含在内，所以在目标检测图中，首先去除与甲状腺腺体层边界相连接的区域以及与B超图像左右边界相连接的区域，然后获取各目标的轮廓。这样基本可以去除绝大多数的伪目标，部分可以直接去掉所有的伪目标。提高了定位甲状腺结节的准确率和效率。附图说明图1是本专利技术实施例提供的一种超声甲状腺结节自动定位方法的流程图；图2是本专利技术实施例提供的B超图像分层示意图；图3(a)是本专利技术实施例提供的水平方向平均灰度值分布图；图3(b)是本专利技术实施例提供的甲状腺腺体层定位结果图；图4(a)是本专利技术实施例提供的去除斑点噪声后的B超图像；图4(b)是本专利技术实施例提供的初始目标检测结果图；图5是本专利技术实施例提供的B超图像的相位不对称图；图6是本专利技术实施例提供的目标分离结果示意图；图7是本专利技术实施例提供的位置和尺寸约束进行目标选择的流程图；图8是本专利技术实施例提供的甲状腺结节定位结果示意图。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本专利技术，并不用于限定本专利技术。此外，下面所描述的本专利技术各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。图1是本专利技术实施例提供的一种超声甲状腺结节自动定位方法的流程图；包括：(1)在现有技术中，将图像水平方向最大平均强度的行索引值和从上到下遍历第一个水平方向最小水平方向平均强度的行索引值之间的部分作为甲状腺腺体层，经过实验，我们发现这样定位的甲状腺腺体层会出现过大或过小的情况，过大不能造本文档来自技高网...

【技术保护点】
一种超声甲状腺结节自动定位方法，其特征在于，包括：(1)根据B超图像的水平方向平均灰度值的分布，得到甲状腺腺体层；(2)在甲状腺腺体层中，利用甲状腺结节分类器，得到初始目标检测图；(3)从B超图像获取相位不对称图，将初始目标检测图在相位不对称图中检测到的目标位置全部作为非目标，得到目标检测图；(4)在目标检测图中，获取各目标的轮廓，计算各目标的中心位置和尺寸，对各目标按照中心位置离目标检测图的中心位置由近到远的顺序进行结节目标检测，当目标的尺寸大于阈值时，该目标为结节目标区域，其他目标为伪目标，进而得到甲状腺结节定位结果。

【技术特征摘要】
1.一种超声甲状腺结节自动定位方法，其特征在于，包括：(1)根据B超图像的水平方向平均灰度值的分布，得到甲状腺腺体层；(2)在甲状腺腺体层中，利用甲状腺结节分类器，得到初始目标检测图；(3)从B超图像获取相位不对称图，将初始目标检测图在相位不对称图中检测到的目标位置全部作为非目标，得到目标检测图；(4)在目标检测图中，获取各目标的轮廓，计算各目标的中心位置和尺寸，对各目标按照中心位置离目标检测图的中心位置由近到远的顺序进行结节目标检测，当目标的尺寸大于阈值时，该目标为结节目标区域，其他目标为伪目标，进而得到甲状腺结节定位结果。2.如权利要求1所述的一种超声甲状腺结节自动定位方法，其特征在于，所述步骤(1)中的B超图像为对原始B超图像进行对比度增强处理后的图像。3.如权利要求1或2所述的一种超声甲状腺结节自动定位方法，其特征在于，所述步骤(1)的具体实现方式为：根据B超图像的水平方向平均灰度值的分布，得到水平方向最大平均灰度值的行索引值Rmax和从下到上遍历第一个水平方向最小平均灰度值的行索引值Rmin，进而得到Rmax与Rmin之差*90％的行索引值Rupper和Rmax与Rmin之差*25％的行索引值Rlower，取Rmax与Rupper之间的大者作为甲状腺腺体层的上边界，取Rmin与Rlower之间的小者作为腺体层的下边界，得到甲状...

【专利技术属性】
技术研发人员：桑农，都文鹏，高常鑫，张士伟，邵远杰，皮智雄，韩楚楚，
申请(专利权)人：华中科技大学，
类型：发明
国别省市：湖北,42