CN110956632B - 钼靶图像中胸大肌区域自动检测方法及装置 - Google Patents
钼靶图像中胸大肌区域自动检测方法及装置 Download PDFInfo
- Publication number
- CN110956632B CN110956632B CN202010002800.2A CN202010002800A CN110956632B CN 110956632 B CN110956632 B CN 110956632B CN 202010002800 A CN202010002800 A CN 202010002800A CN 110956632 B CN110956632 B CN 110956632B
- Authority
- CN
- China
- Prior art keywords
- image
- feature map
- module
- pectoralis major
- molybdenum target
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 210000000062 pectoralis major Anatomy 0.000 title claims abstract description 171
- ZOKXTWBITQBERF-UHFFFAOYSA-N Molybdenum Chemical compound [Mo] ZOKXTWBITQBERF-UHFFFAOYSA-N 0.000 title claims abstract description 139
- 229910052750 molybdenum Inorganic materials 0.000 title claims abstract description 139
- 239000011733 molybdenum Substances 0.000 title claims abstract description 139
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000001514 detection method Methods 0.000 claims abstract description 82
- 210000003205 muscle Anatomy 0.000 claims abstract description 42
- 230000011218 segmentation Effects 0.000 claims abstract description 40
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 claims abstract description 21
- 238000012549 training Methods 0.000 claims abstract description 18
- 238000012805 post-processing Methods 0.000 claims abstract description 17
- 238000007781 pre-processing Methods 0.000 claims abstract description 17
- 210000005075 mammary gland Anatomy 0.000 claims abstract description 8
- 238000010586 diagram Methods 0.000 claims description 79
- 238000005070 sampling Methods 0.000 claims description 32
- 230000009466 transformation Effects 0.000 claims description 14
- 238000004422 calculation algorithm Methods 0.000 claims description 12
- 230000006870 function Effects 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 9
- 238000013527 convolutional neural network Methods 0.000 claims description 8
- 238000010606 normalization Methods 0.000 claims description 8
- 230000002194 synthesizing effect Effects 0.000 claims description 8
- 238000001914 filtration Methods 0.000 claims description 5
- 238000009499 grossing Methods 0.000 claims description 5
- 238000000605 extraction Methods 0.000 claims description 4
- 230000004927 fusion Effects 0.000 claims description 4
- 238000013519 translation Methods 0.000 claims description 4
- 238000003745 diagnosis Methods 0.000 abstract description 7
- 238000004458 analytical method Methods 0.000 abstract description 6
- 238000004364 calculation method Methods 0.000 abstract description 6
- 201000010099 disease Diseases 0.000 abstract description 4
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 abstract description 4
- 210000000481 breast Anatomy 0.000 description 13
- 238000004195 computer-aided diagnosis Methods 0.000 description 9
- 238000004891 communication Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000010339 dilation Effects 0.000 description 3
- 210000001519 tissue Anatomy 0.000 description 3
- 206010006187 Breast cancer Diseases 0.000 description 2
- 208000026310 Breast neoplasm Diseases 0.000 description 2
- 210000000038 chest Anatomy 0.000 description 2
- 238000002156 mixing Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 210000002976 pectoralis muscle Anatomy 0.000 description 2
- 206010028980 Neoplasm Diseases 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 238000010205 computational analysis Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30068—Mammography; Breast
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Biophysics (AREA)
- Software Systems (AREA)
- Evolutionary Computation (AREA)
- Computational Linguistics (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Analysis (AREA)
Abstract
本发明实施例提供一种钼靶图像中胸大肌区域自动检测方法及装置,所述方法包括:对待检测的钼靶图像进行预处理,获得经过预处理的钼靶图像;将所述经过预处理的钼靶图像输入至预先训练好的胸大肌区域检测模型中,获得胸大肌区域的概率图像;对所述胸大肌区域的概率图像进行后处理,获得胸大肌区域分割结果;其中,所述胸大肌区域检测模型是基于经过预处理的钼靶图像样本和所述钼靶图像样本对应的胸大肌分割金标准进行训练获得的。本发明实施例可为现有的乳腺CAD系统提供准确有效地钼靶图像胸大肌检测结果,从而提高计算分析结果的准确度,为放射科医生提供更好的辅助诊断,最终减少病例病情的漏报和误报。
Description
技术领域
本发明涉及图像处理技术领域,更具体地,涉及一种钼靶图像中胸大肌区域自动检测方法及装置。
背景技术
乳腺癌是女性癌症死亡的主要原因之一。目前,乳腺癌早期检测的最有效的方法是乳腺钼靶图像的筛查。然而,目前放射科医生每天分析数百张乳腺钼靶图像,任务繁重耗时耗力,这导致了部分病例病情的漏报和误报。计算机辅助诊断(CAD)系统作为“第二读者意见”有其一致性、可靠性且速度快的优点,因此越来越受欢迎。CAD系统可以为放射科医生提供有效的辅助诊断意见,这能提高了乳腺钼靶图像检测的准确性。
在乳腺的CAD系统中,胸大肌区域的自动检测是钼靶图像计算分析的关键的预处理步骤。准确的胸大肌区域分割对CAD系统诊断有许多帮助,例如,能减少钼靶图像乳腺密度估计的误差,为乳腺钼靶图像多视图分析的图像配准提供参数,提高CAD系统后续计算分析辅助诊断的效果等。然而,胸大肌区域的准确自动分割是一个极具挑战的任务,特别是在筛查的钼靶图像中伪影、乳房皮肤线的低对比度,胸部和乳房组织之间的均匀性等,都极大增加了胸大肌区域准确分割的难度。
目前,钼靶图像中胸大肌区域检测方法中,现有的技术都是基于经典计算机视觉算法实现的。这些胸大肌区域检测方法分为以下几类。
1)基于灰度的胸大肌检测。这一类检测方法应用乳房组织和胸大肌之间的灰度差异,其原理是假定胸大肌的灰度高于其周围组织,利用此信息来进行胸大肌的分割,如简单的使用阈值分割,使用梯度跟踪等等。
2)基于线检测技术的胸大肌检测。这一类检测方法原理是假定胸大肌边界是一条直线,然后使用检测直线的方法进行胸大肌的分割。常用的直线检测方法比如霍夫变换被用来找霍夫空间中正确梯度幅度和方向,投影长度和相应的胸部面积。
3)基于统计分析方法的胸大肌检测。乳房结构根据其纹理、外观和密度被分为不同的类别进行研究。利用此信息,可以推断得到几个统计参数,从而可以运用一些统计技术对乳房密度模式进行分类。比如使用纹理特性来识别胸部边界,方法中可用支持向量回归来识别胸肌的端点。
4)基于小波的胸大肌检测方法。这一类方法运用小波剧透从图形中获取空间频率信息的能力来对钼靶图像胸大肌进行检测。比如使用Dyadic小波分解进行胸肌检测。
5)其他检测方法。比如基于主动轮廓的分割方法用来检测胸大肌。
从上面的阐述可知,目前,钼靶图像中胸大肌区域检测方法中,现有的技术都是基于经典算法实现的,在大多数情况下,这些方法都是具有较强的使用前提假设的,它们更多的是针对一组特定类型的数据或者特定的问题,具有较强的使用条件,无法很好地实现胸大肌的自动检测。
发明内容
本发明的主要目的在于克服现有技术的缺点与不足,提供一种通用且准确有效的钼靶图像中胸大肌区域自动检测方法及装置,为现有的乳腺CAD系统提供准确的钼靶图像胸大肌检测结果,从而提高计算分析结果的准确度,为放射科医生提供更好的辅助诊断,最终减少病例病情的漏报和误报。
第一方面,本发明实施例提供一种钼靶图像中胸大肌区域自动检测方法,包括:
对待检测的钼靶图像进行预处理,获得经过预处理的钼靶图像;
将所述经过预处理的钼靶图像输入至预先训练好的胸大肌区域检测模型中,获得胸大肌区域的概率图像;
对所述胸大肌区域的概率图像进行后处理,获得胸大肌区域分割结果;
其中,所述胸大肌区域检测模型是基于经过预处理的钼靶图像样本和所述钼靶图像样本对应的胸大肌分割金标准进行训练获得的。
第二方面,本发明实施例提供一种钼靶图像中胸大肌区域自动检测装置,包括:
预处理模块,用于对待检测的钼靶图像进行预处理,获得经过预处理的钼靶图像;
特征提取模块,用于将所述经过预处理的钼靶图像输入至预先训练好的胸大肌区域检测模型中,获得胸大肌区域的概率图像;
后处理模块,用于对所述胸大肌区域的概率图像进行后处理,获得胸大肌区域分割结果;
其中,所述胸大肌区域检测模型是基于经过预处理的钼靶图像样本和所述钼靶图像样本对应的胸大肌分割金标准进行训练获得的。
第三方面,本发明实施例提供一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现如第一方面所提供的钼靶图像中胸大肌区域自动检测方法的步骤。
第四方面,本发明实施例提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现如第一方面所提供的钼靶图像中胸大肌区域自动检测方法的步骤。
本发明实施例提供的钼靶图像中胸大肌区域自动检测方法及装置,可为现有的乳腺CAD系统提供准确有效地钼靶图像胸大肌检测结果,从而提高计算分析结果的准确度,为放射科医生提供更好的辅助诊断,最终减少病例病情的漏报和误报。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的钼靶图像中胸大肌区域自动检测方法的流程示意图;
图2为本发明实施例提供的对待检测的钼靶图像进行预处理的流程示意图;
图3为本发明实施例提供的生成拉普拉斯金字塔残差图的流程示意图;
图4为本发明实施例提供的对所述胸大肌区域的概率图像进行后处理的流程示意图;
图5为本发明实施例提供的训练获得所述胸大肌区域检测模型的流程示意图
图6为本发明实施例所构建的基于深度卷积神经网络的胸大肌区域检测模型的结构示意图;
图7为本发明实施例提供的SE模块的结构示意图;
图8为本发明实施例提供的轻量化卷积模块的结构示意图;
图9为本发明实施例提供的注意力网络模块的结构示意图;
图10为本发明实施例提供的胸大肌区域自动检测方法最终实现的效果图;
图11为本发明实施例提供的钼靶图像中胸大肌区域自动检测装置的结构示意图;
图12为本发明实施例提供的电子设备的实体结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,为本发明实施例提供的钼靶图像中胸大肌区域自动检测方法的流程示意图,包括:
步骤100、对待检测的钼靶图像进行预处理,获得经过预处理的钼靶图像;
具体地,图2为本发明实施例提供的对待检测的钼靶图像进行预处理的流程示意图,具体包括:
步骤200、利用边界检测算法对所述待检测的钼靶图像进行乳腺边界区域提取和修正,获得修正后的最小包围矩形区域图像;
在一个实施例中,可以采用基于Canny算子的边界检测算法提取乳腺边界区域,获得该区域的最小包围矩形区域mask,对原始图像(即待检测的钼靶图像)进行裁剪,得到修正的mask区域,修正的目的是保留图像的有效尺寸以减少输入图像的尺寸。
步骤201、对所述修正后的最小包围矩形区域图像分别进行高斯噪声扰动、滤波平滑、基于不同比例的图像增强、基于指数的图像灰度变换和基于对数的灰度变换,得到第一至第五图像p1,p2,p3,p4,p5,合成所述第一图像至第五图像,生成第六图像p;
具体地,对修正后的最小包围矩形区域图像分别依概率阈值选择施加高斯噪声扰动、滤波平滑、基于不同比例的图像增强、基于指数的图像灰度变换和基于对数的灰度变换的操作,分别得到第一至第五图像p1,p2,p3,p4,p5,合成所述第一图像至第五图像,生成第六图像p:
p=ω0p0+ω1p1+ω2p2+ω3p3+ω4p4+ω5p5
其中,ω0,ω1,ω2,ω3,ω4,ω5权重生成方式为:随机生成5个(0,1)的数。
步骤202、对所述第六图像p依次进行平移、旋转变换、尺寸归一化和灰度域归一化,获得第七图像q;
在一个实施例中将图像尺寸归一化到512x512x1,将图像灰度域归一化到[0,1]。
步骤203、计算所述第七图像q分别在x方向、y方向和对角线方向的梯度图,获得第八至第十图像g1,g2,g3;
步骤204、计算所述第七图像q的拉普拉斯金字塔残差图,获得第十一至第十六图像L1~L6;
图3为本发明实施例提供的生成拉普拉斯金字塔残差图的流程示意图。参考图3即可计算获得第十一至第十六图像L1~L6。
步骤205、将所述第七图像q、第八至第十图像g1,g2,g3以及第十一图像L1按照通道进行合成,获得经过预处理的钼靶图像n。
在一个实施例中,钼靶图像n的尺寸为512x512x5。
步骤101、将所述经过预处理的钼靶图像输入至预先训练好的胸大肌区域检测模型中,获得胸大肌区域的概率图像;
具体地,本发明实施例中,所述胸大肌区域检测模型具体为深度卷积神经网络。
所述胸大肌区域检测模型是利用经过预处理的钼靶图像样本作为输入,利用所述钼靶图像样本对应的胸大肌分割金标准作为理想输出,进行训练获得的。其中,胸大肌分割金标准是具有丰富经验的医生手工分割获得的。训练结束后,胸大肌区域检测模型具备了对钼靶图像的胸大肌区域进行自动检测的能力,能够提供准确的钼靶图像胸大肌检测结果。将经过预处理的钼靶图像输入至预先训练好的胸大肌区域检测模型中,即可以获得胸大肌区域的概率图像。
步骤102、对所述胸大肌区域的概率图像进行后处理,获得胸大肌区域分割结果;
具体地,胸大肌区域检测模型输出的胸大肌区域的概率图像的边缘往往存在锯齿,不够平滑,因此,还需要对胸大肌区域检测模型输出的胸大肌区域的概率图像进行后处理。
图4为本发明实施例提供的对所述胸大肌区域的概率图像进行后处理的流程示意图,后处理包括如下步骤:
将所述胸大肌区域的概率图像转换为二值图像,并提取所述二值图像的边缘轮廓;
基于原图(即钼靶图像)增强后的图像,利用主动轮廓算法对胸大肌边缘进行微调。
可以利用canny算子提取边缘轮廓。
本发明实施例提供了一种通用且准确有效的钼靶图像中胸大肌区域自动检测方法,可为现有的乳腺CAD系统提供准确的钼靶图像胸大肌检测结果,从而提高计算分析结果的准确度,为放射科医生提供更好的辅助诊断,最终减少病例病情的漏报和误报。
基于上述实施例的内容,如图5所示,为本发明实施例提供的训练获得所述胸大肌区域检测模型的流程示意图,具体包括:
步骤500、构建基于深度卷积神经网络的胸大肌区域检测模型;
图6为本发明实施例所构建的基于深度卷积神经网络的胸大肌区域检测模型的结构示意图。
步骤501、获取钼靶图像样本和各所述钼靶图像样本对应的胸大肌分割金标准,对所述钼靶图像样本进行预处理,获得经过预处理的钼靶图像样本;
其中,采用与图2所示的预处理流程相同的方法对钼靶图像样本进行预处理,即包括:
利用边界检测算法对所述钼靶图像样本进行乳腺边界区域提取和修正,获得修正后的最小包围矩形区域图像;
对所述修正后的最小包围矩形区域图像分别进行高斯噪声扰动、滤波平滑、基于不同比例的图像增强、基于指数的图像灰度变换和基于对数的灰度变换,得到第一至第五图像p1,p2,p3,p4,p5,合成所述第一图像至第五图像,生成第六图像p;
对所述第六图像p依次进行平移、旋转变换、尺寸归一化和灰度域归一化,获得第七图像q;
计算所述第七图像q分别在x方向、y方向和对角线方向的梯度图,获得第八至第十图像g1,g2,g3;
计算所述第七图像q的拉普拉斯金字塔残差图,获得第十一至第十六图像L1~L6;
将所述第七图像q、第八至第十图像g1,g2,g3以及第十一图像L1按照通道进行合成,获得经过预处理的钼靶图像样本n。
步骤502、将所述经过预处理的钼靶图像样本输入至所述胸大肌区域检测模型中,获得输出图像;
具体地,对所构建的基于深度卷积神经网络的胸大肌区域检测模型进行训练。
将所述经过预处理的钼靶图像样本输入至所述胸大肌区域检测模型中,获得输出图像,具体包括:
步骤a1,将所述经过预处理的钼靶图像样本n输入至膨胀卷积模块,获得第一特征图e1,将所述第一特征图e1输入至SE模块,利用所述SE模块的输出结果更新所述第一特征图e1;
其中,膨胀卷积模块用于执行以下操作:kernel为3x3,stride为1,膨胀系数为2,padding为2的膨胀卷积,然后执行BN和ReLU。
图7为本发明实施例提供的SE模块的结构示意图。在一个实施例中,第一特征图e1为512x512x64的特征图。
步骤a2,将所述第一特征图e1输入至轻量化卷积模块,得到第二特征图e2,对所述经过预处理的钼靶图像样本n进行隔行采样得到第一采样图像n1,将所述第一采样图像n1、第十二图像L2以及第二特征图e2依通道连接后输入至第一卷积模块,利用所述第一卷积模块的输出结果更新所述第二特征图e2;
其中,轻量化卷积模块具体为depthwise残差模块,用于减小网络,降低网络运行时的显存需求和加快推断速度。图8为本发明实施例提供的轻量化卷积模块的结构示意图。在一个实施例中,第二特征图e2为256x256x128的特征图。第一卷积模块执行1x1x134x128的卷积操作,在本发明实施例中卷积操作均默认包括以下操作:对特征图依次执行卷积、BN和ReLU。
步骤a3,将所述第二特征图e2输入至所述轻量化卷积模块,得到第三特征图e3,对所述第一采样图像n1进行隔行采样得到第二采样图像n2,将所述第二采样图像n2、第十三图像L3以及第三特征图e3依通道连接后输入至第二卷积模块,利用所述第二卷积模块的输出结果更新第三特征图e3,将所述第三特征图e3输入至第一反卷积模块,将所述第一反卷积模块输出的特征图与所述第二特征图e2依通道连接,将连接得到的特征图输入至第三卷积模块,并将所述第三卷积模块的输出结果输入至所述SE模块,利用所述SE模块的输出结果更新所述第二特征图e2;
在一个实施例中,第三特征图e3为128x128x256的特征图。第二卷积模块执行1x1x262x256的卷积操作;第一反卷积模块执行kernel为2,步长为2,输出通道为128的反卷积操作;第三卷积模块执行1x1x256x128的卷积操作。
步骤a4,将所述第三特征图e3输入至所述轻量化卷积模块,得到第四特征图e4,对所述第二采样图像n2进行隔行采样得到第三采样图像n3,将所述第三采样图像n3,第十四图像L4以及第四特征图e4依通道连接后输入至第四卷积模块,利用所述第四卷积模块的输出结果更新所述第四特征图e4,将所述第四特征图e4输入至第二反卷积模块,将所述第二反卷积模块输出的特征图与所述第三特征图e3依通道连接,将连接得到的特征图输入至第五卷积模块,并将所述第五卷积模块的输出结果输入至所述SE模块,利用所述SE模块的输出结果更新所述第三特征图e3;
在一个实施例中,所述第四特征图e4为64x64x512的特征图,所述第四卷积模块执行1x1x518x512的卷积操作,得到一个64x64x512的特征图;所述第二反卷积模块执行kernel为2,步长为2,输出通道为256的反卷积操作,连接得到的特征图尺寸为128x128x512,第五卷积模块执行1x1x512x256的卷积操作。
步骤a5,将所述第四特征图e4输入至所述轻量化卷积模块,得到第五特征图e5,对所述第三采样图像n3进行隔行采样得到第四采样图像n4,将所述第四采样图像n4,第十五图像L5以及第五特征图e5依通道连接后输入至第六卷积模块,利用所述第六卷积模块的输出结果更新所述第五特征图e5,将所述第五特征图e5输入至第三反卷积模块,将所述第三反卷积模块输出的特征图与所述第四特征图e4依通道连接,将连接得到的特征图输入至第七卷积模块,并将所述第七卷积模块的输出结果输入至所述SE模块,利用所述SE模块的输出结果更新所述第四特征图e4;
在一个实施例中,所述第五特征图e5为32x32x1024的特征图,所述第六卷积模块执行1x1x1030x1024的卷积操作,所述第三反卷积模块执行kernel为2,步长为2,输出通道为1024的反卷积操作,所述第七卷积模块执行1x1x1024x512卷积模块操作。
步骤a6,将所述第五特征图e5输入至所述轻量化卷积模块,得到第六特征图e6,将所述第六特征图输入至第四反卷积模块,将所述第四反卷积模块输出的特征图与所述第五特征图e5依通道连接后输入至第八卷积模块,将所述第八卷积模块的输出结果输入至所述SE模块,利用所述SE模块的输出结果更新所述第五特征图;
在一个实施例中,所述第六特征图e6为16x16x2048的特征图,所述第四反卷积模块执行kernel为2,步长为2,输出通道为1024的反卷积模块操作,所述第八卷积模块执行1x1x2048x1024的卷积操作。
步骤a7,对所述第四采样图像n4进行隔行采样得到第五采样图像n5,将所述第五采样图像n5、第十六图像L6以及第六特征图e6依通道连接后输入至第九卷积模块,将所述第九卷积模块的输出结果输入至所述SE模块,获得第七特征图d6;
在一个实施例中,所述第九卷积模块执行1x1x2054x2048的卷积操作,得到16x16x2048的特征图。
步骤a8,将所述第七特征图d6输入至第五反卷积模块,得到第八特征图d5;
在一个实施例中,所述第五反卷积模块执行kernel为2,步长为2,输出通道为1024的反卷积模块,得到的第八特征图d5为32x32x1024的特征图。
步骤a9,将所述第七特征图d6和第五特征图e5输入至注意力网络模块,利用所述注意力网络模块的输出结果更新所述第五特征图e5,将所述第五特征图e5和第八特征图d5依通道连接后输入至第十卷积模块,利用所述第十卷积模块的输出结果更新所述第八特征图d5;
图9为本发明实施例提供的注意力网络模块的结构示意图。注意力网络模块结构包括:对输入e执行一个We的运算(首先执行一个3x3xcex1(ce为e的通道数)的卷积操作,而后经过一个BN操作)得到特征图e1,对输入d执行一个Wd的运算(首先执行一个3x3xcd(cd为d的通道数),然后执行一个BN的操作)得到特征图d1,将e1和d1依通道相加,得到新的特征图,以该特征图更新e1;将e1执行一个ReLU和1x1x1x1的卷积操作,然后执行BN,最后经过一个Sigmoid输出,以该概率图输出更新e1;将e1依通道和e中的每个通道的特征相乘,得到的特征图和e依通道相加,得到新的特征图,以该特征图更新e。
在一个实施例中,所述注意力网络模块的输出结果为一个32x32x1024的特征图,所述第十卷积模块执行3x3x2048x1024的卷积操作。
步骤a10,将所述第八特征图d5输入至第六反卷积模块,得到第九特征图d4;
其中,所述第六反卷积模块执行kernel为2,步长为2,输出通道为512的反卷积操作。
步骤a11,将所述第四特征图e4和第八特征图d5输入至所述注意力网络模块,利用所述注意力网络模块的输出结果更新所述第四特征图e4,将所述第四特征图e4和第九特征图d4依通道连接后输入至第十一卷积模块,利用所述第十一卷积模块的输出结果更新所述第九特征图d4;
在一个实施例中,所述第十一卷积模块执行3x3x1024x512的卷积操作。
步骤a12,将所述第九特征图d4输入至第七反卷积模块,获得第十特征图d3;
其中,所述第七反卷积模块执行kernel为2,步长为2,输出通道为256的反卷积操作。
步骤a13,将所述第三特征图e3和第九特征图d4输入至所述注意力网络模块,利用所述注意力网络模块的输出结果更新所述第三特征图e3,将所述第三特征图e3和第十特征图d3依通道连接后输入至第十二卷积模块,利用所述第十二卷积模块的输出结果更新所述第十特征图d3;
在一个实施例中,所述注意力网络模型的输出结果为一个64x64x256的特征图,第十二卷积模块执行3x3x512x256的卷积操作。
步骤a14,将所述第十特征图d3输入至第八反卷积模块,获得第十一特征图d2;
在一个实施例中,所述第八反卷积模块执行kernel为2,步长为2,输出通道为128的反卷积操作。
步骤a15,将所述第二特征图e2和第十特征图d3输入至所述注意力网络模块,利用所述注意力网络模块的输出结果更新所述第二特征图e2,将所述第二特征图e2和第十一特征图d2依通道连接后输入至第十三卷积模块,利用所述第十三卷积模块的输出结果更新所述第十一特征图d2;
在一个实施例中,所述注意力网络模型的输出结果为一个256x256x128的特征图,第十三卷积模块执行3x3x256x128的卷积操作。
步骤a16,将所述第十一特征图d2输入至第九反卷积模块,获得第十二特征图d1;
在一个实施例中,所述第九反卷积模块执行kernel为2,步长为2,输出通道为64的反卷积操作。
步骤a17,将所述第一特征图e1和第十一特征图d2输入至所述注意力网络模块,利用所述注意力网络模块的输出结果更新所述第一特征图e1,将所述第一特征图e1和第十二特征图d1依通道连接后输入至第十四卷积模块,利用所述第十四卷积模块的输出结果更新所述第十二特征图d1;
在一个实施例中,所述注意力网络模型的输出结果为一个512x512x64的特征图,第十四卷积模块执行3x3x128x64的卷积操作。
步骤a18,将所述第七特征图d6输入至第六卷积操作模块T6,获得第六输出图像f6,将所述第八特征图d5输入至第五卷积操作模块T5,获得第五输出图像f5,将所述第九特征图d4输入至第四卷积操作模块T4,获得第四输出图像f4,将所述第十特征图d3分别输入至第三卷积操作模块T3_1,T3_2,获得第三输出图像f3,将所述第十一特征图d2分别输入至第二卷积操作模块T2_1,T2_2,获得第二输出图像f2,将所述第十二特征图d1分别输入至第一卷积操作模块T1_1,T1_2,获得第一输出图像f1,
其中,第六卷积操作模块T6、第五卷积操作模块T5、第四卷积操作模块T4、第三卷积操作模块T3_1,T3_2、第二卷积操作模块T2_1,T2_2以及第一卷积操作模块T1_1,T1_2执行的具体操作如表1所示。
表1 T系列卷积操作模块功能
步骤a19,将f2,f3,f4,f5,f6缩放到和f1相同的尺寸,f1,f2,f3,f4,f5,f6依权重融合得到输出图像f0;
f0=ω1f1+ω2f2+ω3f3+ω4f4+ω5f5+ω6f6
在一个实施例中,ω1,ω2,ω3,ω4,ω5,ω6设定为0.629,0.524,0.419,0.314,0.21,0.105。
在一个实施例中,ω1,ω2,ω3设定为0.801,0.535,0.267。
值得说明的是,在训练结束后,利用所述胸大肌区域检测模型进行胸大肌区域提取时不需要执行步骤a19和a20,胸大肌区域检测模型的最后输出为f1。
步骤503、基于所述输出图像和所述钼靶图像样本对应的胸大肌分割金标准,计算损失函数值;
具体地,本发明实施例对损失函数进行设计。
首先,将所述钼靶图像样本对应的胸大肌分割金标准m0放缩到和输出图像f1,f2,f3,f4,f5,f6相同的尺寸,得到m1,m2,m3,m4,m5,m6,对m0,m1,m2,m3提取胸大肌的边界得到h0,h1,h2,h3;
其中,y′为每像素输出的概率值;y为对应的像素点在胸大肌分割金标准上的标签值;α,γ分别设定为0.25,2;m为m0,m1,m2,m3,m4,m5,m6,h为h0,h1,h2,h3;f为f0,f1,f2,f3,f4,f5,f6,fe为fe,当y属于m0时,ω(m,h)=0.8;y属于m1点时,ω(m,h)=1.0;当y属于m2点时,ω(m,h)=0.6;当y属于m3时,ω(m,h)=0.5;当y属于h0时,ω(m,h)=1.6;当y属于h1时,ω(m,h)=2.0;当y属于h2时,ω(m,h)=1.2;当y属于h3时,ω(m,h)=0.8。
上述取值是为了让网络更加关注于边界,同时由于作为最终实际使用需要的输出,因此给的权重大一些。
步骤504、通过反向传播算法,从所述胸大肌区域检测模型的输出层开始动态调整所述胸大肌区域检测模型的各个参数,以使所述损失函数值朝最小化方向移动;
步骤505、判断是否达到训练结束条件,若是,则保存当前迭代所述胸大肌区域检测模型的参数,获得训练好的胸大肌区域检测模型。
本发明实施例提供的钼靶图像中胸大肌区域自动检测方法,利用深度卷积神经网络来进行钼靶图像中胸大肌区域的自动检测,准确性非常高。本发明在756个钼靶图像上对比有经验医生的手工分割结果显示本方法和医生的手工分割结果间的平均距离方差为0.77mm。图10为本发明实施例提供的胸大肌区域自动检测方法最终实现的效果图。其中,图10中左边为原图,中间为深度学习网络提取图,右边为后处理后的效果图。
如图11所示,为本发明实施例提供的钼靶图像中胸大肌区域自动检测装置的结构示意图,包括:
预处理模块1110,用于对待检测的钼靶图像进行预处理,获得经过预处理的钼靶图像;
特征提取模块1120,用于将所述经过预处理的钼靶图像输入至预先训练好的胸大肌区域检测模型中,获得胸大肌区域的概率图像;
后处理模块1130,用于对所述胸大肌区域的概率图像进行后处理,获得胸大肌区域分割结果;
其中,所述胸大肌区域检测模型是基于经过预处理的钼靶图像样本和所述钼靶图像样本对应的胸大肌分割金标准进行训练获得的。
所述钼靶图像中胸大肌区域自动检测装置用于实现前述钼靶图像中胸大肌区域自动检测方法,因此,所述钼靶图像中胸大肌区域自动检测装置中各功能的理解可参照前述实施例的描述,在此不再赘述。
图12为本发明实施例提供的电子设备的实体结构示意图,如图12所示,该电子设备可以包括:处理器(processor)1210、通信接口(Communications Interface)1220、存储器(memory)1230和通信总线1240,其中,处理器1210,通信接口1220,存储器1230通过通信总线1240完成相互间的通信。处理器1210可以调用存储在存储器1230上并可在处理器1210上运行的计算机程序,以执行上述各方法实施例所提供的钼靶图像中胸大肌区域自动检测方法,例如包括:对待检测的钼靶图像进行预处理,获得经过预处理的钼靶图像;将所述经过预处理的钼靶图像输入至预先训练好的胸大肌区域检测模型中,获得胸大肌区域的概率图像;对所述胸大肌区域的概率图像进行后处理,获得胸大肌区域分割结果;其中,所述胸大肌区域检测模型是基于经过预处理的钼靶图像样本和所述钼靶图像样本对应的胸大肌分割金标准进行训练获得的。
此外,上述的存储器1230中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明实施例的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
本发明实施例还提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现上述各方法实施例提供的钼靶图像中胸大肌区域自动检测方法,例如包括:对待检测的钼靶图像进行预处理,获得经过预处理的钼靶图像;将所述经过预处理的钼靶图像输入至预先训练好的胸大肌区域检测模型中,获得胸大肌区域的概率图像;对所述胸大肌区域的概率图像进行后处理,获得胸大肌区域分割结果;其中,所述胸大肌区域检测模型是基于经过预处理的钼靶图像样本和所述钼靶图像样本对应的胸大肌分割金标准进行训练获得的。
以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (7)
1.一种钼靶图像中胸大肌区域自动检测方法,其特征在于,包括:
对待检测的钼靶图像进行预处理,获得经过预处理的钼靶图像;
将所述经过预处理的钼靶图像输入至预先训练好的胸大肌区域检测模型中,获得胸大肌区域的概率图像;
对所述胸大肌区域的概率图像进行后处理,获得胸大肌区域分割结果;
其中,所述胸大肌区域检测模型是基于经过预处理的钼靶图像样本和所述钼靶图像样本对应的胸大肌分割金标准进行训练获得的;
训练获得所述胸大肌区域检测模型的步骤,具体包括:
构建基于深度卷积神经网络的胸大肌区域检测模型;
获取钼靶图像样本和各所述钼靶图像样本对应的胸大肌分割金标准,对所述钼靶图像样本进行预处理,获得经过预处理的钼靶图像样本;对所述钼靶图像样本进行预处理,获得经过预处理的钼靶图像样本,具体包括:
利用边界检测算法对所述钼靶图像样本进行乳腺边界区域提取和修正,获得修正后的最小包围矩形区域图像;
对所述修正后的最小包围矩形区域图像分别进行高斯噪声扰动、滤波平滑、基于不同比例的图像增强、基于指数的图像灰度变换和基于对数的灰度变换,得到第一至第五图像p1,p2,p3,p4,p5,合成所述第一图像至第五图像,生成第六图像p;
对所述第六图像p依次进行平移、旋转变换、尺寸归一化和灰度域归一化,获得第七图像q;
计算所述第七图像q分别在x方向、y方向和对角线方向的梯度图,获得第八至第十图像g1,g2,g3;
计算所述第七图像q的拉普拉斯金字塔残差图,获得第十一至第十六图像L1~L6;
将所述第七图像q、第八至第十图像g1,g2,g3以及第十一图像L1按照通道进行合成,获得经过预处理的钼靶图像样本n;
将所述经过预处理的钼靶图像样本输入至所述胸大肌区域检测模型中,获得输出图像;
基于所述输出图像和所述钼靶图像样本对应的胸大肌分割金标准,计算损失函数值;
通过反向传播算法,从所述胸大肌区域检测模型的输出层开始动态调整所述胸大肌区域检测模型的各个参数,以使所述损失函数值朝最小化方向移动;
判断是否达到训练结束条件,若是,则保存当前迭代所述胸大肌区域检测模型的参数,获得训练好的胸大肌区域检测模型。
2.根据权利要求1所述的钼靶图像中胸大肌区域自动检测方法,其特征在于,将所述经过预处理的钼靶图像样本输入至所述胸大肌区域检测模型中,获得输出图像,具体包括:
将所述经过预处理的钼靶图像样本n输入至膨胀卷积模块,获得第一特征图e1,将所述第一特征图e1输入至SE模块,利用所述SE模块的输出结果更新所述第一特征图e1;
将所述第一特征图e1输入至轻量化卷积模块,得到第二特征图e2,对所述经过预处理的钼靶图像样本n进行隔行采样得到第一采样图像n1,将所述第一采样图像n1、第十二图像L2以及第二特征图e2依通道连接后输入至第一卷积模块,利用所述第一卷积模块的输出结果更新所述第二特征图e2;
将所述第二特征图e2输入至所述轻量化卷积模块,得到第三特征图e3,对所述第一采样图像n1进行隔行采样得到第二采样图像n2,将所述第二采样图像n2、第十三图像L3以及第三特征图e3依通道连接后输入至第二卷积模块,利用所述第二卷积模块的输出结果更新第三特征图e3,将所述第三特征图e3输入至第一反卷积模块,将所述第一反卷积模块输出的特征图与所述第二特征图e2依通道连接,将连接得到的特征图输入至第三卷积模块,并将所述第三卷积模块的输出结果输入至所述SE模块,利用所述SE模块的输出结果更新所述第二特征图e2;
将所述第三特征图e3输入至所述轻量化卷积模块,得到第四特征图e4,对所述第二采样图像n2进行隔行采样得到第三采样图像n3,将所述第三采样图像n3,第十四图像L4以及第四特征图e4依通道连接后输入至第四卷积模块,利用所述第四卷积模块的输出结果更新所述第四特征图e4,将所述第四特征图e4输入至第二反卷积模块,将所述第二反卷积模块输出的特征图与所述第三特征图e3依通道连接,将连接得到的特征图输入至第五卷积模块,并将所述第五卷积模块的输出结果输入至所述SE模块,利用所述SE模块的输出结果更新所述第三特征图e3;
将所述第四特征图e4输入至所述轻量化卷积模块,得到第五特征图e5,对所述第三采样图像n3进行隔行采样得到第四采样图像n4,将所述第四采样图像n4,第十五图像L5以及第五特征图e5依通道连接后输入至第六卷积模块,利用所述第六卷积模块的输出结果更新所述第五特征图e5,将所述第五特征图e5输入至第三反卷积模块,将所述第三反卷积模块输出的特征图与所述第四特征图e4依通道连接,将连接得到的特征图输入至第七卷积模块,并将所述第七卷积模块的输出结果输入至所述SE模块,利用所述SE模块的输出结果更新所述第四特征图e4;
将所述第五特征图e5输入至所述轻量化卷积模块,得到第六特征图e6,将所述第六特征图输入至第四反卷积模块,将所述第四反卷积模块输出的特征图与所述第五特征图e5依通道连接后输入至第八卷积模块,将所述第八卷积模块的输出结果输入至所述SE模块,利用所述SE模块的输出结果更新所述第五特征图;
对所述第四采样图像n4进行隔行采样得到第五采样图像n5,将所述第五采样图像n5、第十六图像L6以及第六特征图e6依通道连接后输入至第九卷积模块,将所述第九卷积模块的输出结果输入至所述SE模块,获得第七特征图d6;
将所述第七特征图d6输入至第五反卷积模块,得到第八特征图d5;
将所述第七特征图d6和第五特征图e5输入至注意力网络模块,利用所述注意力网络模块的输出结果更新所述第五特征图e5,将所述第五特征图e5和第八特征图d5依通道连接后输入至第十卷积模块,利用所述第十卷积模块的输出结果更新所述第八特征图d5;
将所述第八特征图d5输入至第六反卷积模块,得到第九特征图d4;
将所述第四特征图e4和第八特征图d5输入至所述注意力网络模块,利用所述注意力网络模块的输出结果更新所述第四特征图e4,将所述第四特征图e4和第九特征图d4依通道连接后输入至第十一卷积模块,利用所述第十一卷积模块的输出结果更新所述第九特征图d4;
将所述第九特征图d4输入至第七反卷积模块,获得第十特征图d3;
将所述第三特征图e3和第九特征图d4输入至所述注意力网络模块,利用所述注意力网络模块的输出结果更新所述第三特征图e3,将所述第三特征图e3和第十特征图d3依通道连接后输入至第十二卷积模块,利用所述第十二卷积模块的输出结果更新所述第十特征图d3;
将所述第十特征图d3输入至第八反卷积模块,获得第十一特征图d2;
将所述第二特征图e2和第十特征图d3输入至所述注意力网络模块,利用所述注意力网络模块的输出结果更新所述第二特征图e2,将所述第二特征图e2和第十一特征图d2依通道连接后输入至第十三卷积模块,利用所述第十三卷积模块的输出结果更新所述第十一特征图d2;
将所述第十一特征图d2输入至第九反卷积模块,获得第十二特征图d1;
将所述第一特征图e1和第十一特征图d2输入至所述注意力网络模块,利用所述注意力网络模块的输出结果更新所述第一特征图e1,将所述第一特征图e1和第十二特征图d1依通道连接后输入至第十四卷积模块,利用所述第十四卷积模块的输出结果更新所述第十二特征图d1;
将所述第七特征图d6输入至第六卷积操作模块T6,获得第六输出图像f6,将所述第八特征图d5输入至第五卷积操作模块T5,获得第五输出图像f5,将所述第九特征图d4输入至第四卷积操作模块T4,获得第四输出图像f4,将所述第十特征图d3分别输入至第三卷积操作模块T3_1,T3_2,获得第三输出图像f3,将所述第十一特征图d2分别输入至第二卷积操作模块T2_1,T2_2,获得第二输出图像f2,将所述第十二特征图d1分别输入至第一卷积操作模块T1_1,T1_2,获得第一输出图像f1,
将f2,f3,f4,f5,f6缩放到和f1相同的尺寸,f1,f2,f3,f4,f5,f6依权重融合得到输出图像f0;
3.根据权利要求2所述的钼靶图像中胸大肌区域自动检测方法,其特征在于,基于所述输出图像和所述钼靶图像样本对应的胸大肌分割金标准,计算损失函数值,具体为:
将所述钼靶图像样本对应的胸大肌分割金标准m0放缩到和输出图像f1,f2,f3,f4,f5,相同的尺寸,得到m1,m2,m3,m4,m5,m6,对m0,m1,m2,m提取胸大肌的边界得到h0,h1,h2,h3;
4.根据权利要求1所述的钼靶图像中胸大肌区域自动检测方法,其特征在于,对所述胸大肌区域的概率图像进行后处理,获得胸大肌区域分割结果,具体为:
将所述胸大肌区域的概率图像转换为二值图像,并提取所述二值图像的边缘轮廓;
基于原图增强后的图像,利用主动轮廓算法对胸大肌边缘进行微调。
5.一种钼靶图像中胸大肌区域自动检测装置,其特征在于,包括:
预处理模块,用于对待检测的钼靶图像进行预处理,获得经过预处理的钼靶图像;
特征提取模块,用于将所述经过预处理的钼靶图像输入至预先训练好的胸大肌区域检测模型中,获得胸大肌区域的概率图像;
后处理模块,用于对所述胸大肌区域的概率图像进行后处理,获得胸大肌区域分割结果;
其中,所述胸大肌区域检测模型是基于经过预处理的钼靶图像样本和所述钼靶图像样本对应的胸大肌分割金标准进行训练获得的;
训练获得所述胸大肌区域检测模型的步骤,具体包括:
构建基于深度卷积神经网络的胸大肌区域检测模型;
获取钼靶图像样本和各所述钼靶图像样本对应的胸大肌分割金标准,对所述钼靶图像样本进行预处理,获得经过预处理的钼靶图像样本;对所述钼靶图像样本进行预处理,获得经过预处理的钼靶图像样本,具体包括:
利用边界检测算法对所述钼靶图像样本进行乳腺边界区域提取和修正,获得修正后的最小包围矩形区域图像;
对所述修正后的最小包围矩形区域图像分别进行高斯噪声扰动、滤波平滑、基于不同比例的图像增强、基于指数的图像灰度变换和基于对数的灰度变换,得到第一至第五图像p1,p2,p3,p4,p5,合成所述第一图像至第五图像,生成第六图像p;
对所述第六图像p依次进行平移、旋转变换、尺寸归一化和灰度域归一化,获得第七图像q;
计算所述第七图像q分别在x方向、y方向和对角线方向的梯度图,获得第八至第十图像g1,g2,g3;
计算所述第七图像q的拉普拉斯金字塔残差图,获得第十一至第十六图像L1~L6;
将所述第七图像q、第八至第十图像g1,g2,g3以及第十一图像L1按照通道进行合成,获得经过预处理的钼靶图像样本n;
将所述经过预处理的钼靶图像样本输入至所述胸大肌区域检测模型中,获得输出图像;
基于所述输出图像和所述钼靶图像样本对应的胸大肌分割金标准,计算损失函数值;
通过反向传播算法,从所述胸大肌区域检测模型的输出层开始动态调整所述胸大肌区域检测模型的各个参数,以使所述损失函数值朝最小化方向移动;
判断是否达到训练结束条件,若是,则保存当前迭代所述胸大肌区域检测模型的参数,获得训练好的胸大肌区域检测模型。
6.一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现如权利要求1至4任一项所述钼靶图像中胸大肌区域自动检测方法的步骤。
7.一种非暂态计算机可读存储介质,其上存储有计算机程序,其特征在于,该计算机程序被处理器执行时实现如权利要求1至4任一项所述钼靶图像中胸大肌区域自动检测方法的步骤。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010002800.2A CN110956632B (zh) | 2020-01-02 | 2020-01-02 | 钼靶图像中胸大肌区域自动检测方法及装置 |
PCT/CN2020/141367 WO2021136368A1 (zh) | 2020-01-02 | 2020-12-30 | 钼靶图像中胸大肌区域自动检测方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010002800.2A CN110956632B (zh) | 2020-01-02 | 2020-01-02 | 钼靶图像中胸大肌区域自动检测方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110956632A CN110956632A (zh) | 2020-04-03 |
CN110956632B true CN110956632B (zh) | 2023-04-18 |
Family
ID=69985504
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010002800.2A Active CN110956632B (zh) | 2020-01-02 | 2020-01-02 | 钼靶图像中胸大肌区域自动检测方法及装置 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN110956632B (zh) |
WO (1) | WO2021136368A1 (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110956632B (zh) * | 2020-01-02 | 2023-04-18 | 广州柏视医疗科技有限公司 | 钼靶图像中胸大肌区域自动检测方法及装置 |
CN112132006B (zh) * | 2020-09-21 | 2022-08-26 | 西南交通大学 | 一种面向耕地保护的林地和建筑物智能提取方法 |
KR102549694B1 (ko) * | 2021-01-15 | 2023-06-29 | 고려대학교 산학협력단 | 딥러닝 기반 자동 대흉근 측정 장치 및 그 방법을 실행하기 위한 프로그램이 기록된 기록매체 |
CN116363155B (zh) * | 2023-05-25 | 2023-08-15 | 南方医科大学南方医院 | 一种胸大肌区域智能分割方法、装置和存储介质 |
CN117876848B (zh) * | 2024-03-13 | 2024-05-07 | 成都理工大学 | 基于改进yolov5的复杂环境落石检测方法 |
CN117893530B (zh) * | 2024-03-14 | 2024-06-04 | 中国医科大学附属第一医院 | 一种基于人工智能的喉部影像分析系统 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101484051B1 (ko) * | 2013-06-21 | 2015-01-20 | 중앙대학교 산학협력단 | 능동적 윤곽선 방법을 이용한 cad 시스템을 위한 전처리 장치 및 방법 |
CN105405105B (zh) * | 2015-10-07 | 2017-07-21 | 南京巨鲨显示科技有限公司 | 针对乳腺钼靶图像的显示器灰度曲线校正系统及方法 |
CN108564561A (zh) * | 2017-12-29 | 2018-09-21 | 广州柏视医疗科技有限公司 | 一种钼靶图像中胸大肌区域自动检测方法 |
CN110211086B (zh) * | 2018-02-28 | 2023-04-28 | 腾讯科技(深圳)有限公司 | 图像分割方法、装置及存储介质 |
CN110136117A (zh) * | 2019-05-15 | 2019-08-16 | 东软医疗系统股份有限公司 | 一种钼靶图像的处理方法、系统、图像处理机及存储介质 |
CN110956632B (zh) * | 2020-01-02 | 2023-04-18 | 广州柏视医疗科技有限公司 | 钼靶图像中胸大肌区域自动检测方法及装置 |
-
2020
- 2020-01-02 CN CN202010002800.2A patent/CN110956632B/zh active Active
- 2020-12-30 WO PCT/CN2020/141367 patent/WO2021136368A1/zh active Application Filing
Also Published As
Publication number | Publication date |
---|---|
CN110956632A (zh) | 2020-04-03 |
WO2021136368A1 (zh) | 2021-07-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110956632B (zh) | 钼靶图像中胸大肌区域自动检测方法及装置 | |
CN109461495B (zh) | 一种医学图像的识别方法、模型训练的方法及服务器 | |
CN110310287B (zh) | 基于神经网络的危及器官自动勾画方法、设备和存储介质 | |
CN110706246B (zh) | 一种血管图像分割方法、装置、电子设备和存储介质 | |
US9256941B2 (en) | Microcalcification detection and classification in radiographic images | |
CN108154519A (zh) | 眼底图像中血管的分割方法、装置及存储介质 | |
EP1789920A1 (en) | Feature weighted medical object contouring using distance coordinates | |
JP2023540910A (ja) | 病変検出のための共同訓練を伴う接続機械学習モデル | |
Padmavathy et al. | Performance analysis of pre-cancerous mammographic image enhancement feature using non-subsampled shearlet transform | |
Wang et al. | Retinal vessel segmentation approach based on corrected morphological transformation and fractal dimension | |
CN114926374B (zh) | 一种基于ai的图像处理方法、装置、设备及可读存储介质 | |
CN113012086A (zh) | 一种跨模态图像的合成方法 | |
CN116758336A (zh) | 一种基于人工智能的医学图像智能分析系统 | |
Tummala et al. | Liver tumor segmentation from computed tomography images using multiscale residual dilated encoder‐decoder network | |
Sivanandan et al. | A new CNN architecture for efficient classification of ultrasound breast tumor images with activation map clustering based prediction validation | |
CN109460768B (zh) | 一种用于组织病理学显微图像的文本检测与去除方法 | |
Luo et al. | A liver segmentation algorithm based on wavelets and machine learning | |
Hashmi et al. | Contrast Enhancement in Mammograms Using Convolution Neural Networks for Edge Computing Systems | |
Prakoso et al. | Enhancement methods of brain MRI images: A Review | |
Abdushkour et al. | Enhancing fine retinal vessel segmentation: Morphological reconstruction and double thresholds filtering strategy | |
Singh et al. | A novel approach for brain MRI segmentation and image restoration under intensity inhomogeneity and noisy conditions | |
Gunjan et al. | An effective user interface image processing model for classification of Brain MRI to provide prolific healthcare | |
Nithya et al. | An Efficient PM-Multisampling Image Filtering with Enhanced CNN Architecture for Pneumonia Classfication | |
Ghatwary et al. | Liver CT enhancement using fractional differentiation and integration | |
CN112785580A (zh) | 确定血管流速的方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |