CN108186051B - 一种从超声图像中自动测量胎儿双顶径长度的图像处理方法及处理系统 - Google Patents
一种从超声图像中自动测量胎儿双顶径长度的图像处理方法及处理系统 Download PDFInfo
- Publication number
- CN108186051B CN108186051B CN201711435550.6A CN201711435550A CN108186051B CN 108186051 B CN108186051 B CN 108186051B CN 201711435550 A CN201711435550 A CN 201711435550A CN 108186051 B CN108186051 B CN 108186051B
- Authority
- CN
- China
- Prior art keywords
- image
- edge
- area
- skull
- value
- 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
- 210000003754 fetus Anatomy 0.000 title claims abstract description 22
- 238000003672 processing method Methods 0.000 title claims abstract description 11
- 230000001605 fetal effect Effects 0.000 claims abstract description 57
- 210000003625 skull Anatomy 0.000 claims abstract description 53
- 238000012545 processing Methods 0.000 claims abstract description 27
- 238000005259 measurement Methods 0.000 claims abstract description 26
- 238000000034 method Methods 0.000 claims abstract description 19
- 210000004556 brain Anatomy 0.000 claims abstract description 13
- 230000002708 enhancing effect Effects 0.000 claims abstract description 9
- 238000012549 training Methods 0.000 claims description 29
- 210000001103 thalamus Anatomy 0.000 claims description 16
- 238000001514 detection method Methods 0.000 claims description 15
- 238000004422 calculation algorithm Methods 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 11
- 238000005516 engineering process Methods 0.000 claims description 10
- 238000012937 correction Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 6
- 238000012015 optical character recognition Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 5
- 210000004720 cerebrum Anatomy 0.000 claims description 4
- 210000001638 cerebellum Anatomy 0.000 claims description 3
- 238000003708 edge detection Methods 0.000 claims description 2
- 230000000542 thalamic effect Effects 0.000 claims description 2
- 238000002604 ultrasonography Methods 0.000 description 18
- 238000010586 diagram Methods 0.000 description 7
- 238000003745 diagnosis Methods 0.000 description 4
- 230000000877 morphologic effect Effects 0.000 description 4
- 230000035935 pregnancy Effects 0.000 description 3
- 238000010845 search algorithm Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 230000003187 abdominal effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 241001270131 Agaricus moelleri Species 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002490 cerebral effect Effects 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000003623 enhancer Substances 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000009546 growth abnormality Effects 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0866—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving foetal diagnosis; pre-natal or peri-natal diagnosis of the baby
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4411—Device being modular
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Gynecology & Obstetrics (AREA)
- Pregnancy & Childbirth (AREA)
- Physiology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种从超声图像中自动测量胎儿双顶径(BPD)长度的图像处理方法及处理系统，实现过程为：输入待处理的超声图像，获取超声图像常规拍摄参数，图像边缘增强使颅骨边缘及脑中线边缘清晰化，以及胎儿颅骨边缘拟合，得到双顶径的长度。本发明实现了胎儿双顶径的自动化处理，相较于主观测量，其精确性和稳定性有了明显的提高。
Description
技术领域
本发明涉及医疗超声诊断技术领域，涉及产科超声图像的自动化分析和测量技术，尤其涉及一种从超声图像中自动测量胎儿双顶径长度的图像处理方法及处理系统。
背景技术
超声诊断是目前最重要的医学影像诊断方式之一，具有实时、无损、非侵入等特点，在产科诊断中有广泛的应用。随着图像处理和模式识别技术的发展，将超声图像的处理与识别技术应用于临床诊断成为当前研究的热点。
中孕期常规胎儿超声检查的主要目的是为医务人员提供准确的诊断信息，尽可能给母亲和胎儿带来最佳的产前处理和最好的妊娠结局。通过检查可以确定孕周以及进行胎儿大小的测量，以便在后期的妊娠中及时发现生长异常。通过测量胎儿双顶径(biparietaldiameter，BPD)、头围(Head circumference，HC)、腹围(Abdominal circumference，AC)或腹径(Abdominal diameter，AD)，以及股骨骨干长度(Femoral diaphysis length， FDL)可以估计孕龄及胎儿大小。
由于超声图像信噪比低，目前的检测和测量以肉眼观察、手动定位为主。人工测量主要依靠超声医师的经验进行判断，主观性较强。随着计算机技术的发展，医学图像的处理技术逐渐得到了广泛的应用，如何通过上述技术实现胎儿双顶径的自动检测，提高检测准确性，是需要解决的技术问题。
发明内容
为了克服上述问题，本发明人进行了锐意研究，提出一种从超声图像中自动测量胎儿双顶径长度的图像处理方法及处理系统。通过HOG及Adaboost算法建立识别分类器对输入图像进行胎儿头部区域检测，筛选得到满足要求的头部丘脑水平横切面图像，降低后续图像处理无效率；通过对图像区域进行增强和去噪，可以获得颅骨和脑中线清晰的图像，提高颅骨椭圆拟合准确性；采用最小二乘法对颅骨边缘进行椭圆拟合，结合修正公式，生成BPD；通过科学严谨模型建立、图像处理，以及有效的图像评估方法，得到的BPD的准确度极高，从而完成本发明。
本发明的目的在于提供以下技术方案：
(1)一种从超声图像中自动测量胎儿双顶径长度的图像处理方法，该方法包括以下步骤：
步骤1)，输入待处理的超声图像；
步骤2)，获取超声图像常规拍摄参数；
步骤3)，图像边缘增强：对图像区域进行增强和去噪，使图像边缘即颅骨边缘及脑中线边缘清晰化；
步骤4)，生成双顶径值：对胎儿颅骨外缘进行拟合，求出拟合曲线，得到双顶径长度值。
(2)一种用于实施上述从超声图像中自动测量胎儿双顶径长度的图像处理方法的系统，该系统包括：
模型训练模块，建立训练模型，用于对输入图像进行胎儿头部区域检测；
图像处理模块，提取输入图像的常规拍摄参数，实施图像边缘增强处理；
数据生成模块，利用椭圆模型对颅骨外缘进行拟合，求出胎儿颅骨边缘拟合曲线，结合图像像素距离与实际距离的转换关系，得到双顶径长度。
根据本发明提供的一种从超声图像中自动测量胎儿双顶径长度的图像处理方法及处理系统，具有以下有益效果：
(1)本发明中，获取刻度尺信息时通过对刻度尺区域图像进行二值化操作降低了图像信息获取难度；而对该区域的图像进行投影，将刻度点转化为具有最大值的曲线，可方便准确地获得刻度尺上刻度。
(2)本发明中，通过卷积运算获得增强图像，然后依次采用高斯滤波、二值化操作和形态学开运算进行噪声去除，以及使用连通域搜索算法并将连通域面积过小区域内的像素灰度值置为背景灰度值，可以获得优质的用于椭圆拟合的图像，该图像颅骨边缘和脑中线极其清晰，降低了椭圆拟合难度，极大提高了BPD值准确度。
(3)本发明中，通过HOG及Adaboost算法建立识别分类器，获得训练模型，对输入图像进行胎儿头部区域检测；从基础上降低误检率，避免了经验不足的医师将不具备测量条件的超声图像选为待测的超声图像的情况。
(4)本发明建立的基于超声图像视觉特征自动测量胎儿双顶径大小的方法和系统，实现了双顶径的自动化处理，由于不具有主观性，其精确性和稳定性均较高。
附图说明
图1示出本发明一种实施方式中从超声图像中自动测量胎儿双顶径长度的图像处理方法的流程示意图；
图2示出本发明一种实施方式中超声图像区域划分示意图；
图3示出本发明一种优选实施方式中刻度尺区域竖直方向投影示意图；
图4a示出胎儿头部超声图像取反操作后示意图；
图4b示出胎儿头部超声图像经增强处理后示意图；
图4c示出胎儿头部超声图像经增强、去噪(二值化后)处理后示意图；
图5示出8方向增强算子的算法；
图6示出本发明中胎儿头部超声图像进一步经连通域搜索算法处理后示意图；
图7示出本发明中胎儿颅骨边缘拟合曲线示意图；
图8示出本发明一种实施方式中从超声图像中自动测量胎儿双顶径长度的图像处理系统结构示意图；
图9示出本发明一种优选实施方式中ROI区域。
具体实施方式
下面通过附图对本发明进一步详细说明。通过这些说明，本发明的特点和优点将变得更为清楚明确。
如图1所示，本发明的目的是提供了一种从超声图像中自动测量胎儿双顶径长度的图像处理方法，所述方法包括以下步骤：
步骤1)，输入待处理的超声图像；
步骤2)，获取超声图像常规拍摄参数；
步骤3)，图像边缘增强：对图像区域进行增强和去噪，使图像边缘(颅骨边缘及脑中线边缘)清晰化；
步骤4)，生成双顶径(BPD)值：对胎儿颅骨外缘进行拟合，求出拟合曲线，进而得到双顶径长度值。
步骤1)，输入待处理的超声图像。
图像来源为医学超声仪器采集，图像格式为常用光栅图像，如JPG、BMP格式。
在一种优选的实施方式中，可以将超声图像调整为某一设定尺寸(如64*128像素)，便于后续步骤2)中拍摄参数的提取，或图像处理。
步骤2)，获取超声图像常规拍摄参数。
在一种优选的实施方式中，步骤2)包括以下子步骤：
步骤2.1)，对输入的超声图像进行划分，获得刻度尺区域、图像区域和图像放大率区域；
步骤2.2)，获取图像刻度尺和图像放大率，结合刻度尺及放大率求得图像像素距离与实际距离的转换关系。其中，刻度尺用于表示图像中像素距离对应的实际物理尺寸。
对于步骤2.1)中区域的划分，由于不同厂商的超声仪器所显示的界面略有不同，实际使用中划分为几个区域可根据实际情况灵活划分，但前提是可以清楚得到刻度尺区域、图像区域和图像放大率区域，为图像处理以及参数提取提供基础。如图2 所示，该超声图像划分有五个区域，分别为顶部文字区域、底部文字区域、中部图像区域、左侧刻度尺区域，以及右侧文字 (图像放大率)区域，满足图像处理以及参数提取的要求。
在一种优选的实施方式中，步骤2.2)包括以下子步骤：
步骤2.2.1)，对刻度尺区域进行图像截取，对截取的图像进行二值化操作，并将该区域的图像的每一行所有像素点的灰度值进行累加，得到累加值曲线，如图3中曲线所示，累加值曲线最大值对应的位置即为原刻度点对应位置，将相邻累加值曲线波峰位置的纵坐标相减，可以获得任意相邻两刻度点之间的像素距离，为了提高计算准确性，可以取上述距离的平均值δ。该获得累加值曲线的过程定义为投影。
本发明人了解到，要直接探测刻度尺上的刻度是比较困难的，因为刻度点本身就比较小，加上有图像噪声影响，有时会误将噪声判为刻度点。在该步骤中，通过二值化操作将原始的处理起来比较复杂的彩色图像，转化为仅有黑白两种灰度的图像，降低图像处理难度；而经刻度尺投影后将刻度点转化为具有最大值的曲线，可方便准确地获得刻度尺上刻度。
步骤2.2.2)，截取图像放大率区域，利用光学字符识别技术(OCR)由图像放大率区域获得图像放大率f。如图2中所示，当前图像放大率为66％。
步骤2.2.3)，结合刻度尺及放大率信息，可求得图像像素距离d与实际距离D的转换关系，
其中C为刻度尺单位刻度代表的实际物理尺寸，C可以由用户自由配置。
在一种优选的实施方式中，在获取图像刻度尺和图像放大率前，对刻度尺区域和图像放大率区域进行锐化处理。通过锐化处理补偿文字/图像的轮廓，增强文字/图像的边缘及灰度跳变的部分，使图像变得清晰，利于后续刻度尺投影以及OCR技术提取文字的操作。
步骤3)，图像边缘增强：对图像区域进行增强和去噪，使图像边缘(颅骨边缘及脑中线边缘)清晰化。
因超声图像中固有的边缘缺失、斑点噪声和伪影等影响，在超声图像中进行胎儿颅骨边缘检测存在一定困难。将胎儿头部区域超声图像取反操作(使图像中黑白色彩相互转换)后可以看到，图像上较暗区域大部分是头骨特征，而背景区域则较亮。因此检测头骨区域就等价于检测图像中环形暗区，如图4a 所示。
为准确检测到颅骨边缘及脑中线边缘，步骤3)通过以下子步骤进行图像处理操作：
步骤3.1)截取图像区域，使用8方向增强算子对头围图像进行增强，利用8个方向上9×9大小的模板同该像素点9×9邻域内的像素点进行卷积运算，将8个方向上卷积运算后的最大值作为该像素点的灰度值。R0°至R157.5°分别为0°至157.5°之间的卷积算子，如图5中所示。
经过8方向卷积运算后，颅骨所在位置的图像得到了增强，如图4b所示，但仍然可以看到增强后头围图像中依然存在大量噪声。
步骤3.2)采用一次高斯滤波，结合二值化操作和形态学开运算进行噪声的去除，二值化操作后结果如图4c所示。
高斯滤波是一种线性平滑滤波，用于消除高斯噪声；二值化操作将图像上的像素点的灰度值设置为0或255，使图像呈现出明显的只有黑白的视觉效果；形态学开运算，就是先腐蚀后膨胀的过程，该方式可以消除图像中的小物体，在纤细点处分离物体，平滑较大物体的边界的同时并不明显的改变其面积。
本发明人考虑到，颅骨边缘及脑中线边缘检测(步骤3) 涉及到较多的图像处理操作，若在该步骤发现图像质量不合格，则相当于无偿消耗了图像处理操作。因而，有必要在图像边缘增强步骤前，实施胎儿头部区域检测。
即，步骤3’)，对输入图像进行胎儿头部区域检测，获得合格的胎儿头部丘脑水平的横切面图像。
合格的胎儿头部丘脑水平的横切面图像标准为：
(i)胎儿头部区域图像应全部显示，并且头部区域占超声图像显示区域的60％以上；
(ii)显示为头部丘脑水平横切面图像，颅骨边缘图像清晰；
(iii)脑中线图像清晰且连接，仅在中间被透明隔腔和丘脑分隔；
(iv)两侧大脑半球对称；
(v)不应看到小脑。
在一种优选的实施方式中，为实施输入图像中胎儿头部区域的检测，可以通过建立训练模型(下述步骤0))，通过调用训练模型对输入图像进行胎儿头部区域检测。若未能检测到合格的胎儿头部丘脑水平的横切面图像，则终止计算，或者读入下一张图像，继续进行检测；若检测到合格的胎儿头部丘脑水平的横切面图像，则进行颅骨边缘及脑中线边缘检测步骤。
步骤0)，训练模型的建立包括以下步骤：
子步骤1)，在系统运行之初，建立胎儿头部区域标准图像库。标准图像格式为常用光栅图像，如JPG、BMP格式，标准图像来源为医学超声仪器采集，标准图像的选择标准与上述“合格的胎儿头部丘脑水平的横切面图像标准”一致。
子步骤2)，采用方向梯度直方图特征(HOG)及Adaboost 分类器算法，建立识别分类器，获得训练模型。分类器的训练在测量系统搭建之前完成。系统运行时，调用训练模型对输入图像进行胎儿头部区域检测。
HOG特征是一种局部区域描述子，能够很好地描述物体的边缘，并且对亮度变化和小量偏移不敏感。HOG特征的提取步骤如下：
①将标准图像划分为若干个单元，每个单元为8*8个像素。考虑到胎儿头围区域的超声图像近似为高亮椭圆形；以及在标注ROI区域(感兴趣区域)时，不可以避免地会选中超声图像左右两侧的背景区域。背景区域与超声图像的边界处会形成稳定地梯度方向，为消除这些影响，背景区域的将以临近区域的超声图像的灰度值进行填充。
②在每个单元内进行梯度统计，形成一维的加权梯度方向直方图。其中，直方图划分为9级，划分区间为0°～360°；
③将多个临近的单元组合成一个block块，求其梯度方向直方图向量。
④采用L2-Norm with Hysteresis threshold方法进行归一化，即将直方图向量中的最大值限制在0.2以下，然后再重新归一化一次。
Adaboost算法是一种分类器算法，其基本思想是利用大量的分类能力一般的简单分类器通过一定方法叠加起来，构成一个分类能力很强的强分类器。
在本发明中，由标准图像库中选取训练样本，首先从超声图像中截取头围ROI图像作为正样本，从非ROI区域中随机截取多副子图作为负样本；训练得到分类器后，即可应用分类器进行胎儿头围区域的定位。
步骤4)，生成双顶径(BPD)值：对胎儿颅骨外缘进行拟合，求出拟合曲线，进而得到双顶径长度值。
在一种优选的实施方式中，在生成BPD值前，对图像进行简洁化处理。经过二值化和形态学开运算得到颅骨边缘图像，可能存在一些孤立的点噪声，以及一些细小的空洞。为了提高后续计算的准确性，需要将其去除。此处使用连通域搜索算法，并将连通域面积过小区域(小于图像面积的百分之一)内的像素灰度值置为背景灰度值，如图6所示。
由于原始输入的图像为头部丘脑水平横切面图像，经后续处理的图像中必然包括完整颅骨边缘。此步骤对颅骨边缘进行椭圆拟合，椭圆的短轴的2倍则对应胎儿头部双顶径BPD的长度。通过图像像素距离与实际距离的转换关系，最终取得BPD 的实际值。拟合的原理采用最小二乘法，简述如下。
在平面坐标系中，椭圆方程的一般式为：
ax2+bxy+cy2+dx+ey+f＝0 (1)
其中，a、b、c、d、e、f分别表示椭圆方程的系数；x表示椭圆上点的横坐标；y表示椭圆上点的纵坐标。
约束条件为：a+c＝1 (2)
根据最小二乘原理，通过确定式(3)的最小值可以确定方程(1)中各个系数的大小。
其中(xi，yi)为图6颅骨区域上的边缘点(亮暗交界处点) 坐标，n为抽取的边缘点个数。
根据极值原理，当函数g取得极小值时，
从而可以得到一个线性方程组，结合约束条件，就可确定方程(1)中各个系数的大小，得到胎儿颅骨边缘拟合曲线，如图7所示。
根据胎儿颅骨边缘拟合曲线，通过计算即可得到BPD的长度lBPD。
本发明了解到，双顶径的测量有多个测量标准，主要有“外缘到内缘”，以及“外缘到外缘”等。其中，“外缘到内缘”是指颅骨的外缘到相对侧颅骨的内缘的测量标准；“外缘到外缘”是指颅骨的外缘到相对侧颅骨的外缘的测量标准。
本发明利用上述边缘点拟合出来的椭圆其实是头骨的中间层。在采用“外缘到外缘”的测量标准时，最终的测量结果需要进行修正，修正公式如下
lBPD＝2a’+t (5)
其中，a’为上述拟合椭圆的短轴，t为头骨边缘的平均厚度。
如图8所示，本发明的另一目的是，提供了一种从超声图像中自动测量胎儿双顶径长度的图像处理系统，该系统包括：
模型训练模块，建立训练模型，用于对输入图像进行胎儿头部区域检测；
图像处理模块，提取输入图像的常规拍摄参数，实施图像边缘增强处理；
数据生成模块，对胎儿颅骨外缘进行拟合，求出拟合曲线，结合图像像素距离与实际距离的转换关系，进而得到双顶径长度值。
在本发明中，模型训练模块包括标准图库子模块和模型训练子模块，其中，
标准图库子模块，用于依据标准图像选取标准，建立胎儿头部区域的标准图像库；
模型训练子模块，采用方向梯度直方图特征(HOG)及Adaboost分类器算法，利用标准图像进行训练，建立识别分类器，得到训练模型。
在本发明中，图像处理模块包括图像输入子模块、图像划分子模块、常规拍摄参数获取子模块、图像检测子模块和图像边缘增强子模块，其中，
图像输入子模块，用于输入待处理的超声图像；
图像划分子模块，对输入图像进行区域划分，获得刻度尺区域、图像区域和图像放大率区域；
常规拍摄参数获取子模块，用于获取图像刻度尺信息和放大率信息；
图像检测子模块，调用训练模型对输入图像进行胎儿头部区域检测；若未能检测到合格的胎儿头部丘脑水平的横切面图像，则终止计算；
图像边缘增强子模块，对图像区域进行增强和去噪，使图像边缘清晰化。
在一种优选的实施方式中，常规拍摄参数获取子模块包括刻度尺获取亚子模块和放大率获取亚子模块，其中，
刻度尺获取亚子模块，对刻度尺区域进行图像截取，对截取的图像进行二值化操作，并将该区域的图像每一行所有像素点的灰度值进行累加，得到累加值曲线，累加值曲线最大值对应的位置即为刻度点对应位置，从而获得任意相邻两刻度点之间的像素距离δ；
放大率获取亚子模块，截取图像放大率区域，利用光学字符识别技术(OCR)由图像放大率区域获得图像放大率值f；
转换关系亚子模块，结合刻度尺及放大率信息，可求得图像像素距离d与实际距离D的转换关系，
其中C为刻度尺单位刻度代表的实际物理尺寸。
在一种优选的实施方式中，图像边缘增强子模块包括图像增强亚子模块和噪声去除亚子模块，其中
图像增强亚子模块，截取图像区域，使用8方向增强算子对图像进行增强；所述8方向增强算子是指利用8个方向上9×9 大小的模板同该像素点9×9邻域内的像素点进行卷积运算，将8 个方向上卷积运算后的最大值作为该像素点的灰度值。
噪声去除亚子模块，采用一次高斯滤波，结合二值化操作和形态学开运算进行噪声的去除。
在本发明中，数据生成模块包括曲线拟合子模块和数据显示子模块，其中，
曲线拟合子模块，通过最小二乘法对颅骨边缘进行椭圆拟合；根据极值原理结合椭圆约束条件确定椭圆方程一般式中各系数值，得到胎儿颅骨边缘拟合曲线。根据胎儿颅骨边缘拟合曲线，通过计算(现有的公式)即可得到BPD的长度lBPD。
数据显示子模块，输出胎儿BPD的测量值。
数据生成模块，还包括数据修正子模块。本发明利用上述边缘点拟合出来的椭圆其实是头骨的中间层。当采用“外缘到外缘”的测量标准时，数据修正子模块，通过修正公式进行双顶径长度换算。修正公式为式(5)。
数据生成模块还包括图像前处理子模块，用于在对颅骨边缘进行椭圆拟合前，对图像中存在的孤立的点噪声和细小的空洞进行去除。具体地，使用连通域搜索算法，并将连通域面积过小区域(小于图像面积的百分之一)内的像素灰度值置为背景灰度值。
实施例
实施例1
模型建立：收集2013年1月～2015年12月间中山大学孙逸仙纪念医院妇产科超声工作站存储的孕中期(18～24+6)超声图像 3000幅。筛选出符合条件的图像，纳入研究的图像标准：①胎头丘脑水平的横切面；②理想的超声入射角与脑中线夹角为 90°；③两侧大脑半球对称；④脑中线回声(大脑镰)连接，仅在中间被透明隔腔和丘脑分隔；⑤不应看到小脑。
经过人工筛选出1321幅符合测量条件的图像，随机挑选出 800幅人工标注双顶径的测量位置及参数值，并作为训练用样本，剩余521幅作为测试样本。具体地，首先从训练用超声图像中截取头围ROI图像作为正样本，从非ROI区域中随机截取多幅子图作为负样本。使用HOG及Adaboost算法训练得到分类器后，即可应用分类器进行胎儿头围区域的定位。
如图9所示，用矩形框所标注用矩行所标注位置为所要研究的ROI区域，以矩形框的左上顶点坐标(x，y)和矩形框的长和宽共同描述。
采用521幅测试样本、以及50幅临床样本，对本发明中模型建立方法进行评估；并通过本发明中“模型建立”—“合格胎儿头部区域图像检测，以及获取超声图像常规拍摄参数”—“图像边缘增强”—“对颅骨外缘进行拟合，生成BPD”的完整步骤，上述样本进行BPD的测量。以人工测量为基准，本发明系统相较于人工测量的测量结果见表1。
作为对比例，通过传统的人工方法，对521幅测试样本、以及50幅临床样本进行BPD测量。测量结果见表1。
表1头围区域定位及双顶径测量结果
以上结合了优选的实施方式对本发明进行了说明，不过这些实施方式仅是范例性的，仅起到说明性的作用。在此基础上，可以对本发明进行多种替换和改进，这些均落入本发明的保护范围内。
Claims (2)
1.一种从超声图像中自动测量胎儿双顶径长度的图像处理方法，其特征在于，该方法包括以下步骤：
通过调用训练模型对输入图像进行胎儿头部区域检测；若未能检测到合格的胎儿头部丘脑水平的横切面图像，则终止计算，或者读入下一张图像，继续进行检测；若检测到合格的胎儿头部丘脑水平的横切面图像，则进行颅骨边缘及脑中线边缘检测步骤；
步骤0)，训练模型的建立包括以下步骤：
子步骤1)，在系统运行之初，建立胎儿头部区域标准图像库；
子步骤2)，采用方向梯度直方图特征及Adaboost分类器算法，建立识别分类器，获得训练模型；分类器的训练在测量系统搭建之前完成；系统运行时，调用训练模型对输入图像进行胎儿头部区域检测；由标准图像库中选取训练样本，首先从超声图像中截取头围ROI图像作为正样本，从非ROI区域中随机截取多副子图作为负样本；训练得到分类器后，应用分类器进行胎儿头围区域的定位；ROI为感兴趣区域；
方向梯度直方图特征的提取步骤如下：
①将标准图像划分为若干个单元，每个单元为8*8个像素；为消除背景区域与超声图像的边界处会形成稳定地梯度方向的影响，背景区域将以临近区域的超声图像的灰度值进行填充；
②在每个单元内进行梯度统计，形成一维的加权梯度方向直方图；其中，直方图划分为9级，划分区间为0°～360°；
③将多个临近的单元组合成一个block块，求其梯度方向直方图向量；
④采用L2-Norm with Hysteresis threshold方法进行归一化，将直方图向量中的最大值限制在0.2以下，然后再重新归一化一次；
步骤1)，输入待处理的超声图像；
步骤2)，获取超声图像常规拍摄参数；
步骤2)包括以下子步骤：
步骤2.1)，对输入的超声图像进行划分，获得刻度尺区域、图像区域和图像放大率区域；
步骤2.2)，获取图像刻度尺和图像放大率，结合刻度尺及放大率求得图像像素距离与实际距离的转换关系；
步骤2.2)包括以下子步骤：
步骤2.2.1)，对刻度尺区域进行图像截取，对截取的图像进行二值化操作，并将该区域的图像的每一行所有像素点的灰度值进行累加，得到累加值曲线，累加值曲线最大值对应的位置即为原刻度点对应位置，将相邻累加值曲线波峰位置的纵坐标相减，获得任意相邻两刻度点之间的像素距离或者上述距离的平均值δ；
步骤2.2.2)，截取图像放大率区域，利用光学字符识别技术由图像放大率区域获得图像放大率f；
步骤2.2.3)，结合刻度尺及放大率信息，求得图像像素距离d与实际距离D的转换关系，
其中C为刻度尺单位刻度代表的实际物理尺寸；
步骤2.2)中还包括，在获取图像刻度尺和图像放大率前，对刻度尺区域和图像放大率区域进行锐化处理；
步骤3)，图像边缘增强：对图像区域进行增强和去噪，使图像边缘即颅骨边缘及脑中线边缘清晰化；
步骤3)包括以下子步骤：
步骤3.1)，截取图像区域，使用8方向增强算子对头围图像进行增强：利用8个方向上9×9大小的模板与对应的头围图像中每一个像素点(i,j)进行卷积运算，将8个方向上卷积运算后的最大值作为该像素点的灰度值；
步骤3.2)，采用一次高斯滤波，结合二值化操作和形态学开运算进行噪声的去除；
步骤4)，生成双顶径值：对胎儿颅骨外缘进行拟合，求出拟合曲线，得到双顶径长度值；
步骤4)中，生成BPD值前，对图像进行简洁化处理，使用连通域搜索算法，并将连通域面积小于图像面积的百分之一的区域内的像素灰度值置为背景灰度值，
对颅骨边缘进行椭圆拟合，椭圆的短轴的2倍则对应胎儿头部双顶径BPD的长度，通过图像像素距离与实际距离的转换关系，最终取得BPD的实际值，
步骤4)中，通过最小二乘法对颅骨边缘进行椭圆拟合；根据极值原理结合椭圆约束条件确定椭圆方程一般式中各系数值，得到胎儿颅骨边缘拟合曲线；根据胎儿颅骨边缘拟合曲线即可得到双顶径的长度lBPD；具体如下：
在平面坐标系中，椭圆方程的一般式为：
ax2+bxy+cy2+dx+ey+f＝0 (1)
其中，a、b、c、d、e、f分别表示椭圆方程的系数；x表示椭圆上点的横坐标；y表示椭圆上点的纵坐标；
约束条件为：a+c＝1 (2)
通过确定式(3)的最小值确定方程(1)中各个系数的大小：
其中(xi，yi)为颅骨区域上的边缘点坐标，n为抽取的边缘点个数；
根据极值原理，当函数g取得极小值时，
结合约束条件，确定方程(1)中各个系数的大小，得到胎儿颅骨边缘拟合曲线；
采用“外缘到外缘”的测量标准时，最终的测量结果需要进行修正，修正公式如下：
lBPD＝2a’+t
其中，a’为上述拟合椭圆的短轴，t为头骨边缘的平均厚度。
2.根据权利要求1所述的方法，其特征在于，在步骤3)图像边缘增强处理之前，对输入图像进行胎儿头部区域检测，获得合格的胎儿头部丘脑水平的横切面图像；
合格的胎儿头部丘脑水平的横切面图像标准为：
(i)胎儿头部区域图像应全部显示，并且头部区域占超声图像显示区域的60％以上；
(ii)显示为头部丘脑水平横切面图像，颅骨边缘图像清晰；
(iii)脑中线图像清晰且连接，仅在中间被透明隔腔和丘脑分隔；
(iv)两侧大脑半球对称；
(v)不应看到小脑。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711435550.6A CN108186051B (zh) | 2017-12-26 | 2017-12-26 | 一种从超声图像中自动测量胎儿双顶径长度的图像处理方法及处理系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711435550.6A CN108186051B (zh) | 2017-12-26 | 2017-12-26 | 一种从超声图像中自动测量胎儿双顶径长度的图像处理方法及处理系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108186051A CN108186051A (zh) | 2018-06-22 |
| CN108186051B true CN108186051B (zh) | 2021-11-30 |
Family
ID=62584306
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711435550.6A Active CN108186051B (zh) | 2017-12-26 | 2017-12-26 | 一种从超声图像中自动测量胎儿双顶径长度的图像处理方法及处理系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108186051B (zh) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108876745B (zh) * | 2018-06-27 | 2020-09-01 | 厦门美图之家科技有限公司 | 图像处理方法及装置 |
| CN109886938B (zh) * | 2019-01-29 | 2023-07-18 | 深圳市科曼医疗设备有限公司 | 一种超声图像血管直径自动测量方法 |
| CN111820948B (zh) * | 2019-04-20 | 2022-03-18 | 无锡祥生医疗科技股份有限公司 | 胎儿生长参数测量方法、系统及超声设备 |
| WO2020215484A1 (zh) * | 2019-04-20 | 2020-10-29 | 无锡祥生医疗科技股份有限公司 | 胎儿颈部透明层厚度测量方法、设备及存储介质 |
| CN111899237B (zh) * | 2020-07-27 | 2024-06-04 | 深圳蓝湘智影科技有限公司 | 标尺精度测量方法、装置、计算机设备和存储介质 |
| CN112155604B (zh) * | 2020-09-24 | 2023-03-31 | 广州爱孕记信息科技有限公司 | 基于胎儿超声图像的胎儿严重畸形检测方法及装置 |
| CN116171131A (zh) * | 2020-11-16 | 2023-05-26 | 深圳迈瑞生物医疗电子股份有限公司 | 早孕期胎儿的超声成像方法和超声成像系统 |
| CN115049661B (zh) * | 2022-08-15 | 2022-11-25 | 深圳华声医疗技术股份有限公司 | 目标结构围度测量方法、装置、超声设备及存储介质 |
Family Cites Families (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5588435A (en) * | 1995-11-22 | 1996-12-31 | Siemens Medical Systems, Inc. | System and method for automatic measurement of body structures |
| CN1846233B (zh) * | 2003-08-06 | 2012-01-04 | 特里格医疗有限公司 | 监测分娩参数的方法和装置 |
| US7783095B2 (en) * | 2005-08-11 | 2010-08-24 | Siemens Medical Solutions Usa, Inc. | System and method for fetal biometric measurements from ultrasound data and fusion of same for estimation of fetal gestational age |
| US20080167581A1 (en) * | 2007-01-10 | 2008-07-10 | Yoav Paltieli | Determining parameters associated with a female pelvis and cervix |
| CN102113897B (zh) * | 2009-12-31 | 2014-10-15 | 深圳迈瑞生物医疗电子股份有限公司 | 一种在图像中提取及测量感兴趣目标的方法及其装置 |
| CN102151149B (zh) * | 2010-12-24 | 2013-01-23 | 深圳市理邦精密仪器股份有限公司 | 一种胎儿超声图像自动测量方法及系统 |
| CN103017672A (zh) * | 2012-12-06 | 2013-04-03 | 中铁大桥局集团武汉桥梁科学研究院有限公司 | 用于桥梁结构的非接触无损检测方法 |
| CN103239249B (zh) * | 2013-04-19 | 2015-04-22 | 深圳大学 | 一种胎儿超声图像的测量方法 |
| CN104463876B (zh) * | 2014-12-12 | 2017-05-03 | 湖南科技大学 | 一种基于自适应滤波的复杂背景下图像多圆快速检测方法 |
| CN105662474B (zh) * | 2016-01-22 | 2018-08-17 | 飞依诺科技(苏州)有限公司 | 胎儿头围超声图像的自动检测方法及检测系统 |
| CN106530237B (zh) * | 2016-09-19 | 2019-06-07 | 中山大学 | 一种图像增强方法 |
| CN106963378A (zh) * | 2017-03-03 | 2017-07-21 | 华南理工大学 | 一种基于电磁定位与超声成像的胎头方位测量方法 |
-
2017
- 2017-12-26 CN CN201711435550.6A patent/CN108186051B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN108186051A (zh) | 2018-06-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108186051B (zh) | 一种从超声图像中自动测量胎儿双顶径长度的图像处理方法及处理系统 | |
| CN108378869B (zh) | 一种从超声图像中自动测量胎儿头围长度的图像处理方法及处理系统 | |
| CN107180421B (zh) | 一种眼底图像病变检测方法及装置 | |
| WO2021139258A1 (zh) | 基于图像识别的细胞识别计数方法、装置和计算机设备 | |
| CN107798679B (zh) | 乳腺钼靶图像乳房区域分割与钙化点检测方法 | |
| CN110772286B (zh) | 一种基于超声造影识别肝脏局灶性病变的系统 | |
| JP6422198B2 (ja) | 画像処理装置、画像処理方法、及び画像処理プログラム | |
| CN106228161B (zh) | 一种指针式表盘自动读数方法 | |
| CN104075965B (zh) | 一种基于分水岭分割的显微图像颗粒粒度测量方法 | |
| US8831311B2 (en) | Methods and systems for automated soft tissue segmentation, circumference estimation and plane guidance in fetal abdominal ultrasound images | |
| CN105825516A (zh) | 一种dicom影像血流分析系统 | |
| CN102860814B (zh) | 一种oct合成眼底图视盘中心定位方法及设备 | |
| CN109035227A (zh) | 对ct图像进行肺部肿瘤检测与诊断的系统 | |
| CN104732520A (zh) | 一种胸部数字影像的心胸比测量算法及系统 | |
| CN115546232A (zh) | 一种肝脏超声图像工作区域提取方法、系统及电子设备 | |
| JP2019164450A (ja) | 画像処理方法、コンピュータプログラムおよび記録媒体 | |
| CN112734743B (zh) | 一种胎儿头围自动测量方法及装置 | |
| WO2020140380A1 (zh) | 一种光学相干断层扫描图像的快速分割方法及装置 | |
| CN110580697B (zh) | 一种从超声视频图像中测量胎儿颈背透明物厚度的视频图像处理方法及系统 | |
| CN110051384B (zh) | 结合医学统计信息的胎儿颈背透明物位置检测方法及系统 | |
| CN110543802A (zh) | 眼底图像中左右眼识别方法与装置 | |
| CN113940704A (zh) | 一种基于甲状腺的肌肉和筋膜检测装置 | |
| KR101334029B1 (ko) | 두경부 근육 추출 방법 및 ｒｏｉ분석 방법 | |
| CN113940702A (zh) | 一种甲状腺结节回声分析装置 | |
| CN111612776A (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 | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: An image processing method and system for automatically measuring fetal biparietal diameter length from ultrasound images Granted publication date: 20211130 Pledgee: Zhuhai Branch of China Construction Bank Co.,Ltd. Pledgor: ZHUHAI APPLETREE BIOTECHNOLOGY CO.,LTD. Registration number: Y2024980042018 |