CN110163907A - 胎儿颈部透明层厚度测量方法、设备及存储介质 - Google Patents

胎儿颈部透明层厚度测量方法、设备及存储介质 Download PDF

Info

Publication number
CN110163907A
CN110163907A CN201910451627.1A CN201910451627A CN110163907A CN 110163907 A CN110163907 A CN 110163907A CN 201910451627 A CN201910451627 A CN 201910451627A CN 110163907 A CN110163907 A CN 110163907A
Authority
CN
China
Prior art keywords
distributed areas
neck
pixel
fetus
hyaline layer
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.)
Granted
Application number
CN201910451627.1A
Other languages
English (en)
Other versions
CN110163907B (zh
Inventor
殷晨
李璐
赵明昌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuxi Chison Medical Technologies Co Ltd
Original Assignee
Wuxi Chison Medical Technologies Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuxi Chison Medical Technologies Co Ltd filed Critical Wuxi Chison Medical Technologies Co Ltd
Priority to CN201910451627.1A priority Critical patent/CN110163907B/zh
Priority to PCT/CN2019/093710 priority patent/WO2020215484A1/zh
Publication of CN110163907A publication Critical patent/CN110163907A/zh
Application granted granted Critical
Publication of CN110163907B publication Critical patent/CN110163907B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/0858—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving measuring tissue layers, e.g. skin, interfaces
    • 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
    • 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
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • G06T7/66—Analysis of geometric attributes of image moments or centre of gravity
    • 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/10132—Ultrasound 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/30044—Fetus; Embryo

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Surgery (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Geometry (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Gynecology & Obstetrics (AREA)
  • Pregnancy & Childbirth (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明涉及超声图像处理技术领域，具体公开了一种胎儿颈部透明层厚度测量方法、设备及存储介质。测量方法包括：获取至少包含胎儿颈部透明层的超声图像；通过卷积神经网络模型识别超声图像中颈部透明层的第一分布区域；根据梯度方式从颈部透明层的第一分布区域中分割出包含颈部透明层的第二分布区域，第二分布区域的识别精度大于第一分布区域的识别精度；遍历第二分布区域中的像素单元，以获取第二分布区域沿宽度方向的最大宽度值。本发明了还提供了一种超声设备及存储介质。本发明实现了颈部透明层厚度的自动测量，准确度高，大大提高了医生的工作效率。

Description

胎儿颈部透明层厚度测量方法、设备及存储介质
技术领域
本发明涉及超声图像处理技术领域，尤其是一种胎儿颈部透明层厚度测量方法、设备及存储介质。
背景技术
随着产前超声诊断技术的不断发展，许多胎儿结构异常在产前得到发现。超声检查作为一项无创、无致畸、方便、快捷及安全的检查方法，在产前诊断与筛查胎儿畸形中发挥着重要的作用。超声检查作为一种监测胎儿宫内生长情况的重要手段之一，现已在临床中运用广泛，从二维超声直至目前所运用的四维超声，不论是在医师的操作技能上，还是在仪器的功能和分辨率上都有很大程度的提高。通过超声测量孕期胎儿的颈部透明物厚度，能协助医生在孕期诊断胎儿染色体异常及胎儿心脏畸形等疾病，是孕期一个关键的检查项目。
现有技术中主要通过超声设备获取动态超声图像，由医生通过肉眼和经验观察，选定一个最佳切面，然后手动标定并测量胎儿颈部透明层的厚度。一方面最佳切面选取很难规范化，医生和医生之间存在差异，即使同一个医生在不同的时间段也会存在差异，另外一方面，手动标定测量的准确度低，因此迫切需要一种更规范化并且更精准的测量方式。
发明内容
本发明旨在至少解决现有技术中存在的技术问题之一，提供一种胎儿颈部透明层厚度测量方法、设备及存储介质，以实现胎儿颈部透明层最佳切面的自动选取、颈部透明层的自动识别和厚度自动测量，提高了测量的规范化、准确度和效率。
本发明的第一方面提供了一种胎儿颈部透明层厚度测量方法，包括：
获取至少包含胎儿颈部透明层的超声图像；
通过卷积神经网络模型识别超声图像中颈部透明层的第一分布区域；
根据梯度方式从颈部透明层的第一分布区域中分割出包含颈部透明层的第二分布区域，第二分布区域的识别精度大于第一分布区域的识别精度；
遍历第二分布区域中的像素单元，以获取第二分布区域沿宽度方向的最大宽度值。
进一步地，遍历第二分布区域中的像素单元，包括：
突出显示第二分布区域的轮廓；
通过像素值方式遍历第二分布区域轮廓内的像素单元。
进一步地，根据梯度方式从颈部透明层的第一分布区域中分割出包含颈部透明层的第二分布区域，包括：
计算第一分布区域旋转至水平位置所需的旋转角度；
若旋转角度不为零，则根据旋转角度对第一分布区域进行仿射变换，将第一分布区域旋转至水平位置。
进一步地，计算第一分布区域旋转至水平位置所需的旋转角度，包括：
获取第一分布区域的像素坐标Pn；
根据像素坐标Pn标记第一分布区域的外接矩形Rn；
将外接矩形Rn沿长度方向平均分割为第一外接矩形Rl和第二外接矩形Rr；
计算第一外接矩形Rl的第一质心坐标Cl和第二外接矩形Rr的第二质心坐标Cr；
根据第一质心坐标Cl和第二质心坐标Cr计算第一分布区域所需旋转角度。
进一步地，根据梯度方式从颈部透明层的第一分布区域中分割出包含颈部透明层的第二分布区域，包括：
获取第一分布区域上下边缘像素坐标BN；
根据第一损失函数计算第一分布区域内每个像素点坐标对应的第一损失值，选取最小第一损失值；
根据第二损失函数计算第一分布区域内每个像素点坐标对应的第二损失值，选取最小第二损失值；
根据最小第一损失值与最小第二损失值对应的像素点坐标通过动态规划算法计算出第二分布区域的上下边缘的坐标。
进一步地，动态规划算法为通过迭代方式反推所有第二分布区域的上下边缘的最优像素点坐标。
进一步地，遍历第二分布区域中的像素单元之前，包括：
将第二分布区域内的颜色像素值填充为255，第二分布区域外的颜色像素值填充为0。
进一步地，通过像素值方式遍历第二分布区域轮廓内的像素单元，包括：
沿第二分布区域长度方向将第二分布区域分割为若干等宽度值的分割像素单元；
遍历若干等间距的分割像素单元，选取面积最大的分割像素单元；
根据选取的面积最大的分割像素单元获取对应的上边缘的最厚像素点坐标和下边缘的最厚像素点坐标；
根据上下边缘的最厚像素点坐标计算胎儿的颈部透明层厚度。
进一步地，根据上下边缘的最厚像素点坐标计算胎儿的颈部透明层厚度，包括：
将上下边缘的最厚像素点坐标还原到原始超声图像中，得到上下边缘最厚像素点的实际坐标；
根据上下边缘最厚像素点的实际坐标计算胎儿的颈部透明层厚度。
本发明的第二个方面还提供了另一种胎儿颈部透明层厚度测量方法，包括：
获取至少包含胎儿颈部透明层的动态超声视频；
通过训练后的第一卷积测量模型从动态超声视频中分割出颈部透明层的第一分布区域，第一卷积测量模型至少包括分类神经网络、检测神经网络和分割神经网络；
根据梯度方式从颈部透明层的第一分布区域中分割出包含颈部透明层的第二分布区域，第二分布区域的识别精度大于第一分布区域的识别精度；
遍历第二分布区域中的像素单元，以获取第二分布区域沿宽度方向的最大宽度值。
进一步地，通过训练后的第一卷积测量模型从动态超声视频中分割出颈部透明层的第一分布区域，包括：
通过分类神经网络从动态超声视频中识别出测量胎儿颈部透明层的最佳单帧超声图像；
通过检测神经网络检测最佳超声图像中像素点为透明层的像素点概率，根据像素点概率标记颈部透明层的位置区域；
通过分割神经网络从颈部透明层的位置区域中分割出颈部透明层的第一分布区域。
进一步地，分类神经网络设有交叉熵损失函数，以从动态超声视频中识别出测量胎儿颈部透明层的最佳单帧超声图像。
进一步地，分割神经网络设有像素级交叉熵损失函数。
本发明的第三个方面提供了一种超声设备，包括:
存储器，用于存储计算机程序；
处理器，用于执行计算机程序，使得处理器执行上述的胎儿颈部透明层厚度测量方法。
本发明的第四个方面提供了一种计算机可读存储介质，计算机可读存储介质中存储有计算机程序，计算机程序被处理器执行时用以实现上述的胎儿颈部透明层厚度测量方法的步骤。
本发明的胎儿颈部透明层厚度测量方法能够通过卷积神经网络模型从超声图像中自动识别胎儿颈部透明层的第一分布区域，再通过梯度方式从颈部透明层的第一分布区域中分割出最佳颈部透明层的第二分布区域，提高了颈部透明层的识别精度，最后遍历第二分布区域中的像素单元，以获取第二分布区域沿宽度方向的最大宽度值，实现了颈部透明层厚度的自动测量，准确度高，大大提高了医生的工作效率。
进一步地，本发明的胎儿颈部透明层厚度测量方法能够针对识别出的颈部透明层的第一分布区域不在水平位置进行旋转处理，不需要医生在采集超声图像时采集水平位置的超声图像，降低了对包含胎儿颈部透明层的超声图像的要求，提高了测量的速度和准确度。
本发明的超声设备能够自动识别出胎儿颈部透明层的第一分布区域，再通过梯度方式从颈部透明层的第一分布区域中分割出最佳颈部透明层的第二分布区域，提高了颈部透明层的识别精度，最后通过遍历第二分布区域中的像素单元，以获取第二分布区域沿宽度方向的最大宽度值，实现了颈部透明层厚度的自动测量，准确度高，大大提高了医生的工作效率。
附图说明
图1为本发明的胎儿颈部透明层厚度测量方法的流程示意图。
图2为本发明的旋转获得水平第一分布区域的流程示意图。
图3为本发明的第一分布区域旋转角度计算的流程示意图。
图4为本发明的计算第二分布区域上下边缘坐标的流程示意图。
图5为本发明的一实施例胎儿颈部透明层厚度测量方法的流程示意图。
图6a为本发明获取的超声图像。
图6b为本发明获取的超声图像已识别出颈部透明层第一分布区域的示意图。
图6c为本发明获取的超声图像已标出第一质心和第二质心的示意图。
图6d为已经旋转获得水平位置的颈部透明层第一分布区域的示意图。
图7为本发明获取的超声图像已识别出颈部透明层第二分布区域的示意图。
图8为本发明获取的超声图像的掩膜区。
图9为本发明的另一胎儿颈部透明层厚度测量方法的流程示意图
图10为本发明一实施例卷积神经网络结构示意图。
图11为本发明的卷积神经网络模型中分割神经网络的结构示意图。
具体实施方式
下面通过具体实施方式结合附图对本发明作进一步详细说明。其中不同实施方式中类似元件采用了相关联的类似的元件标号。在以下的实施方式中，很多细节描述是为了使得本申请能被更好的理解。然而，本领域技术人员可以毫不费力的认识到，其中部分特征在不同情况下是可以省略的，或者可以由其他元件、材料、方法所替代。在某些情况下，本申请相关的一些操作并没有在说明书中显示或者描述，这是为了避免本申请的核心部分被过多的描述所淹没，而对于本领域技术人员而言，详细描述这些相关操作并不是必要的，他们根据说明书中的描述以及本领域的一般技术知识即可完整了解相关操作。另外，说明书中所描述的特点、操作或者特征可以以任意适当的方式结合形成各种实施方式。同时，方法描述中的各步骤或者动作也可以按照本领域技术人员所能显而易见的方式进行顺序调换或调整。因此，说明书和附图中的各种顺序只是为了清楚描述某一个实施例，并不意味着是必须的顺序，除非另有说明其中某个顺序是必须遵循的。本文中为部件所编序号本身，例如“第一”、“第二”等，仅用于区分所描述的对象，不具有任何顺序或技术含义。
目前测量胎儿颈部透明层厚度主要通过超声设备获取超声图像后由医生通过肉眼和经验观察，手动标定后测量胎儿颈部透明层的厚度，准确度低，且依靠医生肉眼和经验工作的方式效率低。
本发明的第一方面提供了一种胎儿颈部透明层厚度测量方法，如图1所示，包括以下步骤：
S100，获取至少包含胎儿颈部透明层的超声图像；
S200，通过卷积神经网络模型识别超声图像中颈部透明层的第一分布区域600；
S300，根据梯度方式从颈部透明层的第一分布区域600中分割出包含颈部透明层的第二分布区域700，第二分布区域的识别精度大于第一分布区域的识别精度；
S400，遍历第二分布区域700中的像素单元，以获取第二分布区域沿宽度方向的最大宽度值。
本发明的胎儿颈部透明层厚度测量方法能够通过卷积神经网络模型从超声图像中自动识别胎儿颈部透明层的第一分布区域600，再通过梯度方式从颈部透明层的第一分布区域600中分割出最佳颈部透明层的第二分布区域700，提高了颈部透明层的识别精度，最后通过遍历第二分布区域700中的像素单元，以获取第二分布区域沿宽度方向的最大宽度值，实现了颈部透明层厚度的自动测量，准确度高，大大提高了医生的工作效率。可以理解的是第一分布区域为颈部透明层的大概分布区域的图像，第二分布区域为颈部透明层的精确分布区域的图像。
步骤S100获取至少包含胎儿颈部透明层的超声图像主要通过超声成像设备的换能器采集。超声成像设备至少包括换能器、超声主机、输入单元、控制单元、和存储器。超声成像设备可以包括显示屏，超声成像设备的显示屏可以为识别系统的显示器。换能器用于发射和接收超声波，换能器受发射脉冲的激励，向目标组织(例如，人体或者动物体内的器官、组织、血管等等)发射超声波，经一定延时后接收从目标区域反射回来的带有目标组织的信息的超声回波，并将此超声回波重新转换为电信号，以获得超声图像或者视频。换能器可以通过有线或无线的方式连接到超声主机。
输入单元用于输入操作人员的控制指令。输入单元可以为键盘、跟踪球、鼠标、触摸面板、手柄、拨盘、操纵杆以及脚踏开关中的至少一个。输入单元也可以输入非接触型信号，例如声音、手势、视线或脑波信号。
控制单元至少可以控制焦点信息、驱动频率信息、驱动电压信息以及成像模式等扫描信息。控制单元根据用户所需成像模式的不同，对信号进行不同的处理，获得不同模式的超声图像数据，然后经对数压缩、动态范围调整、数字扫描变换等处理形成不同模式的超声图像，如B图像，C图像，D图像，多普勒血流图像，包含组织弹性特性的弹性图像等等，或者其他类型的二维超声图像或三维超声图像。
需要理解的是，在一实施例中获取至少包含胎儿颈部透明层的超声图像也可以为存储在存储介质中的超声图像，例如，云服务器，U盘或者硬盘等。
医生通常会根据肉眼和工作经验对获取的超声图像进行观察判断对胎儿颈部透明层进行标定测量，准确度和效率都很低。
本发明步骤S200通过卷积神经网络模型识别超声图像中颈部透明层的第一分布区域600，需要理解的是卷积神经网络模型通过卷积神经网络对已标记的胎儿颈部透明层的若干超声图像进行训练确定。利用训练后的卷积神经网络模型可以自动识别出超声图像中颈部透明层的第一分布区域，第一分布区域600是卷积神经网络模型识别出的大概的分布区域。
本发明训练后的卷积神经网络模型包括输入层、隐含层和输出层；其中隐含层包括若干卷积层、下采样层、上采样层；输入的超声图像先经过若干卷积层和下采样层，分别进行卷积操作和下采样操作，再经过若干卷积层和上采样层，分别进行卷积操作和上采样操作；神经网络的输入层与隐含层、各隐含层之间、隐含层与输出层之间通过权重参数相连接；卷积神经网络模型中的卷积层用于自动提取超声图像中的特征。更优地，卷积神经网络模型的神经网络每次对超声图像进行卷积或采样处理时，从较浅卷积层复制特征到较深的卷积层，复制的特征和较深卷积层的特征对应像素相加后进入下一层卷积层。
针对步骤S200识别出的第一分布区域600并非都是水平的，为了提高计算精度方便计算，首先需要计算出第一分布区域600需要旋转的旋转角度，然后将第一分布区域600旋转至水平位置。如图2所示参见图6a-图6d，具体包括以下步骤：
S210，计算第一分布区域旋转至水平位置所需的旋转角度；
S220，若旋转角度不为零，则根据旋转角度对第一分布区域进行仿射变换，将第一分布区域旋转至水平位置。
在另一实施例中可以先判断颈部透明层的第一分布区域600是否处于水平位置；若第一分布区域不处于水平位置，则计算第一分布区域600旋转至水平位置所需的旋转角度。判断颈部透明层的第一分布区域是否处于水平位置可以通过训练好的神经网络进行判断，也可以通过医护人员对第一分布区域进行人工判断的方式进行判断。
如图3所示，计算旋转角度的方法如下:
S2201，获取第一分布区域600的像素坐标Pn；
S2202，根据像素坐标Pn标记第一分布区域600的外接矩形Rn；
S2203，将外接矩形Rn沿长度方向平均分割为第一外接矩形Rl和第二外接矩形Rr；
S2204，计算第一外接矩形Rl的第一质心坐标Cl和第二外接矩形Rr的第二质心坐标Cr；
S2205，根据第一质心坐标Cl和第二质心坐标Cr计算第一分布区域600所需旋转角度，计算方式如下：
其中，θ表示旋转角度，表示第二质心坐标Cr的y轴坐标值，表示第二质心坐标Cr的x轴坐标值，表示第一质心坐标Cl的y轴坐标值，表示第一质心坐标Cl的x轴坐标值。
S230，根据旋转角度对第一分布区域600进行仿射变换，将第一分布区域旋转至水平位置。需要理解的是，仿射变换，又称仿射映射，是指在几何中一个向量空间进行一次线性变换并接上一个平移，变换为另一个向量空间。根据计算出的旋转角度进行仿射变换得到水平的第一分布区域600。
为了进一步地提高颈部透明层分布区域识别的精度，本发明根据梯度方式从颈部透明层的第一分布区域600中分割出最佳颈部透明层的第二分布区域700。可以理解的是，第二分布区域700的精度大于第一分布区域600的精度。如图4，具体地包括
S310，获取第一分布区域600上下边缘像素坐标BN；
BN＝{p1,p2,p3,pN-1,pN}
N表示超声图像的宽度，p1,p2,p3,pN-1,pN表示超声图像上下边缘第一个像素点的坐标至第N个像素点的坐标，需要理解的时，像素点的坐标包括x轴和y轴。
S320，根据第一损失函数计算第一分布区域600内每个像素点坐标对应的第一损失值，选取最小第一损失值；
第一损失函数：
其中，pj表示第j个像素点的坐标，Zl(pj)表示利用图像拉普拉斯的双边效应，使寻找的像素坐标接近第一分布区域600的下边缘。fadj(pj,pj-1)函数限制了相邻两个像素间距必须要小，从而保证了边缘的连续性；q为p前一个像素点的坐标，t为p后一个像素点的坐标，θ这里为90度，表示只计算y方向的二阶导数。
S330，根据第二损失函数计算第一分布区域600内每个像素点坐标对应的第二损失值，选取最小第二损失值；
第二损失函数：
其中，
Zu(pj)表示利用图像拉普拉斯的双边效应，使寻找的像素坐标接近第一分布区域600的上边缘。q为p前一个像素点的坐标，θ这里为90度，表示计算y方向的二阶导数，为两个像素点坐标距离的sigmoid函数。
S340，根据最小第一损失值与最小第二损失值对应的像素点坐标通过动态规划算法计算出第二分布区域700的上下边缘的坐标。动态规划算法为通过迭代方式反推所有第二分布区域700的上下边缘的最优像素点坐标。
动态规划算法为：
对于下边缘采用如下迭代方式：
对于上边缘采用如下迭代方式：
和分别为上下边缘第一个像素点的损失函数，通过计算每一行N个点的损失函数，得到超声图像最后一列所有像素点的损失函数，可以反推所有第二分布区域700的上下边缘的最优像素点坐标。如图7所示，通过反推出的所有第二分布区域700的上下边缘的最优像素点坐标可以突出显示第二分布区域700的轮廓。突出显示可以采用线条或曲线勾勒出第二部分区域的轮廓，也可以高亮显示第二分布区域700的轮廓。
需要理解的是，胎儿颈部透明层的厚度为透明层最厚的位置的厚度，本发明遍历第二分布区域中的像素单元，以获取第二分布区域沿宽度方向的最大宽度值。在一实施中，先突出显示第二分布区域的轮廓，在通过像素值方式遍历第二分布区域轮廓内的像素单元，以获取第二分布区域沿宽度方向的最大宽度值。即将第二分布区域700的上下边缘之间的最大距离值确定为胎儿的颈部透明层厚度。如图8所示，首先将第二分布区域700内的颜色像素值填充为255，第二分布区域700外的颜色像素值填充为0，以获得第二分布区域700的掩膜区。
如图5所示，通过像素值方式遍历第二分布区域轮廓内的像素单元，包括：
S410，沿第二分布区域700长度方向将第二分布区域700分割为若干等宽度值的分割像素单元；
S420，遍历若干等间距的分割像素单元，选取面积最大的分割像素单元，等宽度值为1个像素；
S430，根据选取的最大分割像素单元获取对应的上边缘的最厚像素点坐标和下边缘的最厚像素点坐标；通过列公式计算：
Ex＝Sx×Wx
高度为 为掩膜中当前x位置对应的像素值为255的最小的纵坐标,则为最大的纵坐标,Wx代表权重,为当前Sx分割像素单元区域在原始图片上面像素的平均值,px,y是原始坐标的像素值，x的范围为第二分布区域700的宽度，Ex表示计算的数值结果。此做法避免了某些情况因为边缘轮廓计算错误导致最终结果在非实际颈部透明层内，有些面积最大的区域不一定是实际厚度最大的位置，该方法理论上也考虑了人眼观察规律，也就是颈部透明带最厚的区域肯定是个低回声区域。
S440，根据上下边缘的最厚像素点坐标计算胎儿的颈部透明层厚度。
上述计算出最大的Ex对应的x位置和其对应的上下两个纵坐标px和即为颈部透明带最厚的区域位置的坐标，由于该坐标旋转过后的，将该坐标的位置反变换到原始空间中得到最终的实际坐标，并根据两点距离公式计算出上下两点的距离，得到最终的厚度。
本发明的第二个方面还提供了另一种胎儿颈部透明层厚度测量方法，如图9所示，包括：
S510，获取至少包含胎儿颈部透明层的动态超声视频；
S520，通过训练后的第一卷积测量模型从动态超声视频中分割出颈部透明层的第一分布区域600，第一卷积测量模型至少包括分类神经网络、检测神经网络和分割神经网络；
S530，根据梯度方式从颈部透明层的第一分布区域中分割出包含颈部透明层的第二分布区域，第二分布区域的识别精度大于第一分布区域的识别精度；
S540，遍历第二分布区域中的像素单元，以获取第二分布区域沿宽度方向的最大宽度值。
本发明能够通过训练后的第一卷积测量模型从动态超声视频中分割出颈部透明层的第一分布区域600。再通过梯度方式从颈部透明层的第一分布区域600中分割出最佳颈部透明层的第二分布区域700，提高了颈部透明层的识别精度，最后遍历第二分布区域中的像素单元识别出透明层最厚的位置测量出胎儿颈部透明层的厚度，实现了颈部透明层厚度的自动测量，准确度高，大大提高了医生的工作效率。可以理解的是第一分布区域为颈部透明层的大概分布区域，第二分布区域为颈部透明层的精确分布区域。
在进行颈部透明层测量之前需要识别出动态超声视频中测量胎儿颈部透明层的最佳单帧超声图像。本发明一实施了通过分类神经网络用于从动态超声视频识别出测量胎儿颈部透明层的最佳单帧超声图像。分类神经网络的输入是动态超声视频中所有完整的超声图像，分类神经网络经过若干卷积层和下采样层后，根据超声图像中的特征输出超声图像的预测结果。分类神经网络的多个卷积层用来自动提取超声图像中的特征；分类神经网络的各卷积层之间、输入和卷积层之间、卷积层和输出之间通过权重参数相连接；设置输入层尺寸，以和输入神经网络的超声图像的尺寸相适配。
分类神经网络至少包括卷积层(conv)、最大池化层(max-pooling)、平均层(avg)、逻辑回归层(softmax)和滤波器。如图10所示，在一实施了中分类神经网络包括5个卷积层、4个最大池化层(max-pooling)、1个平均层(avg)和1个逻辑回归层(softmax)。可选的，设置分类神经网络的输入层尺寸为416*416*1，经过若干个3*3的卷积操作和最大池化操作，然后再经过对每组特征进行平均操作得到输入的超声图像是最佳测量图像与否的概率，最终进行softmax操作，计算方法如下：
其中，i表示图10中第10层输出的第一个数值或第二个数值，等号右边分母表示以e为底数，第10层输出的两个数值分别为指数计算得到的结果的合计；softmaxi表示逻辑回归层操作后输出的概率结果；根据概率结果识别输入的动态超声视频中测量胎儿颈部透明层的最佳超声图像。
通过分类神经网络从动态超声视频中识别出胎儿颈部透明层的最佳单帧超声图像后，通过检测神经网络检测最佳超声图像中像素点为透明层的像素点概率，根据像素点概率标记颈部透明层的位置区域。检查神经网络将像素点概率大于和/或等于预设概率的像素确定为透明层。检测神经网络包括输入层、隐含层、输出层，检测神经网络的输入层与隐含层、各隐含层之间、隐含层与输出层之间通过权重参数相连接；隐含层包括卷积层、最大池化层和结合层；首先是数个卷积层和数个最大池化层交替连接，然后再连接若干卷积层，随后再连接一个结合层，将结合层之前相连接的高级特征层与该高级特征层之前的一层或数层隐含层相结合；该高级特征层与相结合的隐含层的输出图像的长和宽一致；该高级特征层与之前的一层或数层隐含层相结合后一起输入到最后一个卷积层，最后一个卷积层作为输出层。
通过检测神经网络检测最佳超声图像中像素点为透明层的像素点概率，根据像素点概率标记颈部透明层的位置区域后，通过分割神经网络从颈部透明层的位置区域中分割出颈部透明层的第一分布区域600。突出显示可以采用线条或曲线勾勒出第一部分区域的轮廓，也可以高亮显示一分布区域的轮廓。分割神经网络通过突出显示从颈部透明层的位置区域中分割出颈部透明层的第一分布区域600。在一实施了分割神经网络包括输入层、隐含层和输出层；其中隐含层包括若干卷积层、下采样层、上采样层；输入的超声图像先经过若干卷积层和下采样层，分别进行卷积操作和下采样操作，再经过若干卷积层和上采样层，分别进行卷积操作和上采样操作；更优地，分割神经网络对超声图像进行卷积或采样处理时，从较浅卷积层复制特征到较深的卷积层，复制的特征和较深卷积层的特征对应像素相加后进入下一层卷积层；分割神经网络的输入层与隐含层、各隐含层之间、隐含层与输出层之间通过权重参数相连接；卷积层用于自动提取超声图像中的特征。分割出第一分布区域后根据梯度方式从颈部透明层的第一分布区域600中分割出最佳颈部透明层的第二分布区域700；通过像素值加权方式遍历第二分布区域700中的像素单元，将第二分布区域700的上下边缘之间的最大距离值确定为胎儿的颈部透明层厚度，即遍历第二分布区域中的像素单元，以获取第二分布区域沿宽度方向的最大宽度值。
第一卷积测量模型的训练过程如下：
步骤S1，收集胎儿的动态超声视频，并对动态超声视频进行标注，本实施例中优选的标注方法是标注出视频中用于测量胎儿颈部透明层厚度的最佳单帧超声图像；收集测量胎儿颈部透明层厚度的最佳图像，用连续折线形成的闭合曲线标注出胎儿的颈部透明层处；
步骤S2，基于收集到的动态超声视频、超声图像及标注信息建立分类神经网络、检测神经网络和分割神经网络，训练得到第一卷积测量模型。
神经网络的处理流程包括：
步骤S21，将收集到的超声图像划分为训练集、验证集和测试集；
收集的用于分类神经网络的超声图像视频中的超声图像中随机选取3/5的图像作为训练集，随机选取1/5的图像作为验证集；剩余1/5的超声图像作为测试集使用；收集的用于检测神经网络和分割神经网络训练的超声图像同样分为三个集合。训练集用于训练神经网络；验证集用于验证神经网络的效果并帮助选择最优的神经网络模型参数；测试集超用于测神经网络的使用效果。优选地，训练集、验证集和测试集选取的比例为3/5、1/5，1/5。
不同品牌的超声设备或同一品牌不同型号的超声设备采集的动态超声视频中拆分的单帧超声图像的大小是不一样的，需要对超声图像进行预处理。将训练集和验证集的超声图像固定到一定尺寸，并归一化同样尺寸的超声图像；如预处理后的超声图像为256*256*3；256*256表示预处理后超声图像的长和宽，即256像素长，256像素宽；可选地，将超声图像固定到一定尺寸时，保持原始图像的长宽比例，或者改变原始图像的长宽比例；对超声图像进行归一化操作的具体处理方法为将超声图像中每个像素值减去图像像素的均值后除以图像像素的方差；归一化后超声图像的均值是0，方差是1；由于超声图像预处理时超声图像的尺寸发生了变化，超声图像的模板也需要进行相应比例的改变；
步骤S22，建立神经网络结构；
本发明先建立了分类神经网络，用于预测收集到的动态超声视频中哪一帧用来测量胎儿颈部透明层厚度最佳。通过分类神经网络从动态超声视频识别出测量胎儿颈部透明层的最佳单帧超声图像。分类神经网络的输入是动态超声视频中所有完整的超声图像，分类神经网络经过若干卷积层和下采样层后，根据超声图像中的特征输出超声图像的预测结果。分类神经网络的多个卷积层用来自动提取超声图像中的特征；分类神经网络的各卷积层之间、输入和卷积层之间、卷积层和输出之间通过权重参数相连接；设置输入层尺寸，以和输入神经网络的超声图像的尺寸相适配。
在以上分类神经网络之后，本发明建立检测神经网络，用来检测出最佳颈部透明层厚度测量的超声图像中的颈部透明层的部分。检测神经网络的输入是收集到的最佳颈部透明层厚度测量的超声图像，超声图像的曲线标注转化为标注框后用于训练该神经网络；检测神经网络将卷积神经网络末端的高级特征层与前一层或者前几层的低级细粒度特征结合，增加卷积神经网络对偏小的目标对象的检测效果。可选的，检测神经网络的输入设置为416*416*1，输出设置为13*13*35；检测神经网络输出了超声图像中颈部透明层可能的标记框的坐标信息和概率信息；
在一实施中分割神经网络设置的输入层尺寸为256*256*1，经过两个3*3的卷积操作分别得到两个256*256*16的特征，再经过下采样操作得到128*128*16的特征；再继续进行若干3*3卷积和下采样得到64*64*64的特征；再继续进行若干上采样和3*3卷积得到256*256*16的特征；最终进行1*1卷积操作得到256*256*2的预测结果，预测结果中数值的取值范围是0到1，表示超声图像中对应的像素点是胎儿颈部透明层范围内的概率，需要理解的是预测结果为概率。如图11所示，的灰色长方形即表示图像经过每次卷积或者采样操作后提取到的特征，白色长方形表示复制得到的特征；优选的，分割神经网络的卷积选择适合尺度的扩张卷积，以提高网络的感受野，提高网络预测准确率。可选地，网络中的上采样层和下采样层可以去除，仍保证网络的输入层和输出层的长和宽一致；可选的，分割神经网络的输入超声图像可以在检测网络预测的检测框的基础上适当扩大，如上下左右分别扩张20个像素点。
步骤S23，初始化神经网络：将神经网络的权重参数设置为0到1之间的随机数；
步骤S24，计算神经网络的损失函数；
以上设计的分类神经网络的损失函数为交叉熵损失；以上涉及的检测网络的损失函数包括检测框位置和检测框的预测概率的损失；以上涉及的分割神经网络的损失函数选择像素级的交叉熵损失；
分类神经网络的交叉熵计算公式为：
其中，N表示超声图像的总数量；pi表示第i张超声图像是否为最佳测量图像；ti表示分类神经网络预测的第i张超声图像为最佳测量图像的概率；pi和ti的数值越接近，交叉熵损失越小。
检测网络的损失函数包含两大部分，分别为含有目标检测对象的检测框的概率预测的误差和中心坐标、高度、宽度预测的误差；不含有目标检测对象的检测框的概率预测的误差，计算公式为：
其中，λ1-λ3表示各项误差在总的损失函数中占的比重，各项误差都选用平方误差的形式。
Loss_function损失函数的第一项表示含有目标颈部透明层部分对象的检测框的概率预测的误差。其中，S2表示将超声图像划分成S×S个网格单元；B表示每个网格单元设置多少个检测框；表示第i个网格单元的第j个检测框是否含有目标检测对象，若检测框与标注框的交集较大，则认为检测框含有目标检测对象，否则，认为检测框不含有目标检测对象， 表示检测网络对该网格单元当前的第j个检测框的预测概率；第二项表示含有目标对象的检测框的位置和长宽的预测误差。其中xi，yi，hi，wi表示第i个网格单元的标注框的中心位置和长度、宽度信息，表示预测的边界框相应的信息；第三项是不含有目标对象的检测框的概率预测的误差，因为不含有目标对象的边界框占多数，所以λ3通常会设置得比λ1小，否则无法训练得到识别效果较好的网络。可选的，λ1＝5，λ2＝λ3＝1。
分割神经网络的像素级交叉熵损失函数为：
像素级交叉熵损失函数是对超声图像中每个像素上的预测误差的总和；其中，x,y为分割神经网络的输入图像的长度和宽度，pij是超声图像中第i行第j列像素被分割神经网络预测的该像素为预测部位的概率，tij是超声图像中第i行第j列像素在超声图像模板中对应的数值，若该像素为预测部位，取值为1，否则，取值为0；分割神经网络输出的预测的概率与超声图像模板越接近，交叉熵损失函数越小；
步骤S25，训练神经网络，得到第一卷积测量模型；
此步骤中，随机选择训练集中的超声图像并对其进行随机的变换，然后输入神经网络，选择合适的训练迭代次数和批处理大小对神经网络进行训练；可选的，变换操作有旋转、缩放、裁剪、弹性形变等；优选的，本发明中只进行了随机旋转操作；可选的，根据神经网络的损失函数更新网络参数的机制使用自适应矩估计优化方法；
步骤S26，选择最优的神经网络模型参数。
计算不同参数下的三种神经网络在验证集上的预测结果，计算分类准确率，准确率最高对应的参数作为分类神经网络的参数；计算检测网络的预测框和验证集标注转化得到的标注框之间的交并比，交并比的计算方法为：
交并比＝(预测结果∩图形模板)/(预测结果∪图形模板)
交并比的可能取值范围为[0,1]
选择交并比最大情况下的参数作为检测网络的最优参数；计算分割神经网络的预测结果和标注转化得到的验证图像模板之间的交并比，选择交并比最大情况下的参数作为分割神经网络的最优参数；两个对象间的交并比即两者的交集除以两者的并集；
通过经过训练得到的第一卷积测量模型从动态超声视频中分割出颈部透明层的第一分布区域600包括：
步骤S31，将超声图像视频中的所有单帧超声图像固定到与分类神经网络输入层相适配的同样尺寸，对超声图像进行归一化；
步骤S32，将归一化后的单帧超声图像输入训练好的分类神经网络模型，采用最优参数，得到分类神经网络输出的预测为最佳测量帧的超声图像；将最佳测量帧图像输入检测神经网络，检测出超声图像中胎儿颈部透明层的位置区域，再经过分割神经网络得到胎儿颈部透明层第一分布区域。
本发明的第三个方面提供了一种超声设备，包括:存储器，用于存储计算机程序；处理器，用于执行计算机程序，使得处理器执行使得处理器执行上述的胎儿颈部透明层厚度测量方法。
本发明的超声设备能够自动识别出胎儿颈部透明层的第一分布区域600，再通过梯度方式从颈部透明层的第一分布区域600中分割出最佳颈部透明层的第二分布区域700，提高了颈部透明层的识别精度，最后通过像素值加权方式识别出透明层最厚的位置测量出胎儿颈部透明层的厚度，实现了颈部透明层厚度的自动测量，准确度高，大大提高了医生的工作效率。
本发明的第四个方面提供了一种计算机可读存储介质，计算机可读存储介质中存储有计算机程序，计算机程序被处理器执行时用以实现上述的胎儿颈部透明层厚度测量方法的步骤。
最后所应说明的是，以上具体实施方式仅用以说明本发明的技术方案而非限制，尽管参照实例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。

Claims (15)

1.一种胎儿颈部透明层厚度测量方法，其特征在于，包括：
获取至少包含胎儿颈部透明层的超声图像；
通过卷积神经网络模型识别所述超声图像中颈部透明层的第一分布区域；
根据梯度方式从颈部透明层的第一分布区域中分割出包含颈部透明层的第二分布区域，所述第二分布区域的识别精度大于所述第一分布区域的识别精度；
遍历所述第二分布区域中的像素单元，以获取所述第二分布区域沿宽度方向的最大宽度值。
2.如权利要求1所述的胎儿颈部透明层厚度测量方法，其特征在于，所述遍历所述第二分布区域中的像素单元，包括：
突出显示所述第二分布区域的轮廓；
通过像素值方式遍历第二分布区域轮廓内的像素单元。
3.如权利要求1所述的胎儿颈部透明层厚度测量方法，其特征在于，所述根据梯度方式从颈部透明层的第一分布区域中分割出包含颈部透明层的第二分布区域，包括：
计算所述第一分布区域旋转至水平位置所需的旋转角度；
若所述旋转角度不为零，则根据所述旋转角度对所述第一分布区域进行仿射变换，将所述第一分布区域旋转至水平位置。
4.如权利要求3所述的胎儿颈部透明层厚度测量方法，其特征在于，所述计算所述第一分布区域旋转至水平位置所需的旋转角度，包括：
获取所述第一分布区域的像素坐标Pn；
根据所述像素坐标Pn标记所述第一分布区域的外接矩形Rn；
将所述外接矩形Rn沿长度方向平均分割为第一外接矩形Rl和第二外接矩形Rr；
计算所述第一外接矩形Rl的第一质心坐标Cl和所述第二外接矩形Rr的第二质心坐标Cr；
根据所述第一质心坐标Cl和第二质心坐标Cr计算所述第一分布区域所需旋转角度。
5.如权利要求1-4中任一项所述的胎儿颈部透明层厚度测量方法，其特征在于，所述根据梯度方式从颈部透明层的第一分布区域中分割出包含颈部透明层的第二分布区域，包括：
获取第一分布区域上下边缘像素坐标；
根据第一损失函数计算所述第一分布区域内每个像素点坐标对应的第一损失值，选取最小第一损失值；
根据第二损失函数计算所述第一分布区域内每个像素点坐标对应的第二损失值，选取最小第二损失值；
根据所述最小第一损失值与最小第二损失值对应的像素点坐标通过动态规划算法计算出所述第二分布区域的上下边缘的坐标。
6.如权利要求5所述的胎儿颈部透明层厚度测量方法，其特征在于，所述动态规划算法为通过迭代方式反推所有所述第二分布区域的上下边缘的最优像素点坐标。
7.如权利要求1-4中任一项所述的胎儿颈部透明层厚度测量方法，其特征在于，所述遍历所述第二分布区域中的像素单元之前，包括：
将所述第二分布区域内的颜色像素值填充为255，所述第二分布区域外的颜色像素值填充为0。
8.如权利要求2所述的胎儿颈部透明层厚度测量方法，其特征在于，所述通过像素值方式遍历第二分布区域轮廓内的像素单元，包括：
沿所述第二分布区域长度方向将所述第二分布区域分割为若干等宽度值的分割像素单元；
遍历所述若干等间距的分割像素单元，选取面积最大的分割像素单元；
根据选取的面积最大的分割像素单元获取对应的上边缘的最厚像素点坐标和下边缘的最厚像素点坐标；
根据上下边缘的最厚像素点坐标计算胎儿的颈部透明层厚度。
9.如权利要求8所述的胎儿颈部透明层厚度测量方法，其特征在于，所述根据上下边缘的最厚像素点坐标计算胎儿的颈部透明层厚度，包括：
将所述上下边缘的最厚像素点坐标还原到原始超声图像中，得到上下边缘最厚像素点的实际坐标；
根据上下边缘最厚像素点的实际坐标计算胎儿的颈部透明层厚度。
10.一种胎儿颈部透明层厚度测量方法，其特征在于，包括：
获取至少包含胎儿颈部透明层的动态超声视频；
通过训练后的第一卷积测量模型从所述动态超声视频中分割出颈部透明层的第一分布区域，所述第一卷积测量模型至少包括分类神经网络、检测神经网络和分割神经网络；
根据梯度方式从颈部透明层的第一分布区域中分割出包含颈部透明层的第二分布区域，所述第二分布区域的识别精度大于所述第一分布区域的识别精度；
遍历所述第二分布区域中的像素单元，以获取所述第二分布区域沿宽度方向的最大宽度值。
11.根据权利要求10所述的胎儿颈部透明层厚度测量方法，其特征在于，所述通过训练后的第一卷积测量模型从所述动态超声视频中分割出颈部透明层的第一分布区域，包括：
通过所述分类神经网络从动态超声视频中识别出测量胎儿颈部透明层的最佳单帧超声图像；
通过所述检测神经网络检测所述最佳超声图像中像素点为透明层的像素点概率，根据所述像素点概率标记颈部透明层的位置区域；
通过分割神经网络从所述颈部透明层的位置区域中分割出颈部透明层的第一分布区域。
12.根据权利要求11所述的胎儿颈部透明层厚度测量方法，其特征在于，所述分类神经网络设有交叉熵损失函数，以从动态超声视频中识别出测量胎儿颈部透明层的最佳单帧超声图像。
13.根据权利要求10或11所述的胎儿颈部透明层厚度测量方法，其特征在于，所述分割神经网络设有像素级交叉熵损失函数。
14.一种超声设备，其特征在于，包括:
存储器，用于存储计算机程序；
处理器，用于执行所述计算机程序，使得处理器执行如权利要求1至13中任一项所述的胎儿颈部透明层厚度测量方法。
15.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质中存储有计算机程序，所述计算机程序被处理器执行时用以实现如权利要求1至13中任一项所述的胎儿颈部透明层厚度测量方法的步骤。
CN201910451627.1A 2019-04-20 2019-05-28 胎儿颈部透明层厚度测量方法、设备及存储介质 Active CN110163907B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201910451627.1A CN110163907B (zh) 2019-05-28 2019-05-28 胎儿颈部透明层厚度测量方法、设备及存储介质
PCT/CN2019/093710 WO2020215484A1 (zh) 2019-04-20 2019-06-28 胎儿颈部透明层厚度测量方法、设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910451627.1A CN110163907B (zh) 2019-05-28 2019-05-28 胎儿颈部透明层厚度测量方法、设备及存储介质

Publications (2)

Publication Number Publication Date
CN110163907A true CN110163907A (zh) 2019-08-23
CN110163907B CN110163907B (zh) 2021-06-29

Family

ID=67629666

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910451627.1A Active CN110163907B (zh) 2019-04-20 2019-05-28 胎儿颈部透明层厚度测量方法、设备及存储介质

Country Status (1)

Country Link
CN (1) CN110163907B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111062948A (zh) * 2019-11-18 2020-04-24 北京航空航天大学合肥创新研究院 一种基于胎儿四腔心切面图像的多组织分割方法
CN111127305A (zh) * 2019-10-14 2020-05-08 李胜利 基于早孕期胎儿颅面部三维容积自动获取标准切面的方法
CN116012283A (zh) * 2022-09-28 2023-04-25 逸超医疗科技(北京)有限公司 一种全自动超声图像量测方法、设备及存储介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120101383A1 (en) * 2010-10-22 2012-04-26 Samsung Medison Co., Ltd. 3d ultrasound system for intuitive displaying to check abnormality of object and method for operating 3d ultrasound system
CN103263278A (zh) * 2013-01-23 2013-08-28 郑末晶 一种从超声图像上自动测量胎儿颈背透明物厚度的图像处理方法
CN104156967A (zh) * 2014-08-18 2014-11-19 深圳市开立科技有限公司 一种胎儿颈部透明层图像分割方法、装置及系统
CN104268538A (zh) * 2014-10-13 2015-01-07 江南大学 一种易拉罐点阵喷码字符在线视觉检测方法
CN105405119A (zh) * 2015-10-21 2016-03-16 复旦大学 基于深度置信网络和三维模型的胎儿正中矢状面自动检测方法
CN107103587A (zh) * 2017-06-05 2017-08-29 新疆大学 一种生物芯片图像倾斜的校正方法及装置
CN108888294A (zh) * 2018-03-30 2018-11-27 杭州依图医疗技术有限公司 一种检测颈部透明带宽度的方法及装置

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120101383A1 (en) * 2010-10-22 2012-04-26 Samsung Medison Co., Ltd. 3d ultrasound system for intuitive displaying to check abnormality of object and method for operating 3d ultrasound system
CN103263278A (zh) * 2013-01-23 2013-08-28 郑末晶 一种从超声图像上自动测量胎儿颈背透明物厚度的图像处理方法
CN104156967A (zh) * 2014-08-18 2014-11-19 深圳市开立科技有限公司 一种胎儿颈部透明层图像分割方法、装置及系统
CN104268538A (zh) * 2014-10-13 2015-01-07 江南大学 一种易拉罐点阵喷码字符在线视觉检测方法
CN105405119A (zh) * 2015-10-21 2016-03-16 复旦大学 基于深度置信网络和三维模型的胎儿正中矢状面自动检测方法
CN107103587A (zh) * 2017-06-05 2017-08-29 新疆大学 一种生物芯片图像倾斜的校正方法及装置
CN108888294A (zh) * 2018-03-30 2018-11-27 杭州依图医疗技术有限公司 一种检测颈部透明带宽度的方法及装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
A.ANZALONE ET AL.: "A system for the automatic measuremthe nuchal translucency thickness from ultrasound video stream of the foetus", 《CBMS 2013》 *
LAKSHMI. PS ET AL.: "Automated Screening for Trisomy 21 by measuring Nuchal Translucency anf Frontomaillary Facial Angle", 《2018 INTERNATIONAL CONFERENCE ON ADVANCES IN COMPUTING, COMMUNICATIONS》 *
田学东 等: "基于多特征模糊模式是被的公式符号关系判定", 《计算机工程与应用》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111127305A (zh) * 2019-10-14 2020-05-08 李胜利 基于早孕期胎儿颅面部三维容积自动获取标准切面的方法
CN111062948A (zh) * 2019-11-18 2020-04-24 北京航空航天大学合肥创新研究院 一种基于胎儿四腔心切面图像的多组织分割方法
CN111062948B (zh) * 2019-11-18 2022-09-13 北京航空航天大学合肥创新研究院 一种基于胎儿四腔心切面图像的多组织分割方法
CN116012283A (zh) * 2022-09-28 2023-04-25 逸超医疗科技(北京)有限公司 一种全自动超声图像量测方法、设备及存储介质
CN116012283B (zh) * 2022-09-28 2023-10-13 逸超医疗科技(北京)有限公司 一种全自动超声图像量测方法、设备及存储介质

Also Published As

Publication number Publication date
CN110163907B (zh) 2021-06-29

Similar Documents

Publication Publication Date Title
JP6467041B2 (ja) 超音波診断装置、及び画像処理方法
EP1593094B1 (en) Image analysis for the purpose of assessing cancer
CN110163907A (zh) 胎儿颈部透明层厚度测量方法、设备及存储介质
Han et al. A texture feature analysis for diagnosis of pulmonary nodules using LIDC-IDRI database
CN110189255A (zh) 基于两级检测的人脸检测方法
CN107016665A (zh) 一种基于深度卷积神经网络的ct肺结节检测方法
CN102920477A (zh) 医学影像的目标区域边界确定装置和方法
CN102247172A (zh) 自动化胎儿胎龄评估的系统和方法
CN111768366A (zh) 超声成像系统、bi-rads分级方法及模型训练方法
CN105678734B (zh) 一种图像匹配系统的异源测试图像标定方法
CN107832838A (zh) 评价细胞涂片标本满意度的方法和装置
CN102939616B (zh) 组织分类
CN107708550A (zh) 用于检测和测量解剖异常的分段回声结构的表面建模
Shi et al. Automatic estimation of dairy cow body condition score based on attention-guided 3D point cloud feature extraction
CN111820948B (zh) 胎儿生长参数测量方法、系统及超声设备
CN108230313A (zh) 基于部件自适应选择与判别模型的sar图像目标检测方法
US20220061810A1 (en) Systems and methods for placing a gate and/or a color box during ultrasound imaging
Chaudhari et al. Ultrasound image based fully-automated nuchal translucency segmentation and thickness measurement
CN106355578B (zh) 基于卷积神经网络的超声颈动脉远端识别装置和方法
CN116989694A (zh) 基于3d轮廓扫描的尺寸定位检测系统
CN105528791B (zh) 一种面向触摸屏手绘图像的质量评价装置及其评价方法
CN115619941A (zh) 一种超声成像方法及超声设备
CN110033042A (zh) 一种基于深度神经网络的直肠癌环周切缘mri图像自动识别方法及系统
Fang et al. Reliable thyroid carcinoma detection with real-time intelligent analysis of ultrasound images
Ke et al. An automatic instrument recognition approach based on deep convolutional neural network

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