CN113892967B - 一种超声图像的处理方法和超声仪器 - Google Patents
一种超声图像的处理方法和超声仪器 Download PDFInfo
- Publication number
- CN113892967B CN113892967B CN202111041389.0A CN202111041389A CN113892967B CN 113892967 B CN113892967 B CN 113892967B CN 202111041389 A CN202111041389 A CN 202111041389A CN 113892967 B CN113892967 B CN 113892967B
- Authority
- CN
- China
- Prior art keywords
- tissue
- image
- target
- preset
- gray
- 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
- 238000003672 processing method Methods 0.000 title abstract description 6
- 238000003384 imaging method Methods 0.000 claims abstract description 137
- 230000011218 segmentation Effects 0.000 claims abstract description 98
- 238000009499 grossing Methods 0.000 claims abstract description 53
- 238000002604 ultrasonography Methods 0.000 claims description 65
- 238000012545 processing Methods 0.000 claims description 60
- 238000000034 method Methods 0.000 claims description 59
- 230000006835 compression Effects 0.000 claims description 29
- 238000007906 compression Methods 0.000 claims description 29
- 239000011159 matrix material Substances 0.000 claims description 20
- 238000012937 correction Methods 0.000 claims description 10
- 210000000056 organ Anatomy 0.000 claims description 10
- 230000008569 process Effects 0.000 description 27
- 238000010586 diagram Methods 0.000 description 10
- 238000004590 computer program Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 7
- 238000003860 storage Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 4
- 238000000265 homogenisation Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 238000003708 edge detection Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 101100313373 Aliivibrio fischeri (strain ATCC 700601 / ES114) tfoX1 gene Proteins 0.000 description 1
- 101100313375 Aliivibrio fischeri (strain ATCC 700601 / ES114) tfoX2 gene Proteins 0.000 description 1
- 101100043727 Caenorhabditis elegans syx-2 gene Proteins 0.000 description 1
- 230000003187 abdominal effect Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000001364 causal effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000249 desinfective effect Effects 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000002203 pretreatment Methods 0.000 description 1
- 230000001954 sterilising effect Effects 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
Classifications
-
- 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
-
- 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)
- Radiology & Medical Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Biomedical Technology (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开了一种超声图像的处理方法和超声仪器,该设备的数据获取单元被配置为:获取超声图像;处理器被配置为:对获取到的超声图像中各个像素点在预设方向上的梯度进行平滑处理,得到超声图像对应的目标梯度图像;基于目标梯度图像中各个像素点的灰度值、预设灰度阈值和超声图像中包含的组织匹配的组织占比常数确定组织分割阈值;应用组织分割阈值对超声图像进行组织分割得到组织图像;根据预先设置的每组成像参数和目标组织亮度值,分别确定每组成像参数对应的组织图像的亮度误差,根据亮度误差选择一组成像参数为目标成像参数。将目标成像参数配置到超声仪器中,实现了超声图像亮度的均匀化,操作简单,准确性高。
Description
技术领域
本发明涉及图像处理技术领域,尤其涉及一种超声图像的处理方法和超声仪器。
背景技术
超声波是一种波长极短的机械波,易于获得各向异性的声能,可用于清洗、碎石、杀菌消毒等,在医学、工业上有很多的应用。其应用原理通常是,使用超声探头发射超声波给物体,记录物体内部结构的回波,将回波进行处理而形成灰度图像,以反映物体的内部结构。
在医学领域,通常应用超声波来反应人体内部的结构,比如组织或者器官的状态。而增益调整是超声成像过程中的重要环节,通常情况下,超声***中为每个部位或组织设置了最佳的增益以及动态范围,而最佳的增益和动态范围适用于绝大部分患者,但是在实际的检测过程中,由于个体差异较大(比如声波能量衰减不同),会存在一部分人在默认设置下的图像均匀性较差,表现为图像亮暗不均,此时,需要医生手动调节增益和动态范围。
鉴于医生手动调节过程中复杂且效率低,相关技术中,大多数的超声***可以自动根据患者的超声图像进行自动增益调整,以达到图像亮度均匀化的目的。但是,相关技术中的做法主要有如下问题:采用求取结构张量的方式进行组织区域判断,通过对结构张量和构建函数值判定噪声和非噪声区域、组织和非组织区域,该种方法计算较为复杂,单纯通过特征值并不一定能完全判断出组织区域;简单采用固定阈值方式对图像分块,导致对组织的区域判定存在误差,不能最大程度的考虑所有的组织区域的亮度,与理想的组织亮度会存在误差。
发明内容
本发明示例性的实施方式中提供一种超声图像的处理方法和超声仪器,将处理后的成像参数配置到超声仪器,可以实现超声图像亮度的均匀化,操作简单,准确性高。
根据示例性的实施方式中的第一方面,提供一种超声图像的处理方法,包括:
对获取到的超声图像中各个像素点在预设方向上的梯度进行平滑处理,得到所述超声图像对应的目标梯度图像;其中,所述超声图像为包含人体器官或组织的图像;
基于所述目标梯度图像中各个像素点的灰度值、预设灰度阈值和所述超声图像中包含的组织匹配的组织占比常数确定组织分割阈值;
应用所述组织分割阈值对所述超声图像进行组织分割,得到组织图像;
根据预先设置的每组成像参数和目标组织亮度值,分别确定每组成像参数对应的所述组织图像的亮度误差,根据所述亮度误差选择一组成像参数为目标成像参数;其中,每组成像参数包括动态范围与增益。
本申请实施例,通过平滑处理可以去除获取到的超声图像各个预设方向上的噪声,进而增加组织分割的准确性;另外,在确定组织分割阈值时,考虑了不同部位的组织特性,应用与当前超声图像中包含的组织匹配的组织占比常数,这样得到的组织分割阈值更准确,进而也增加了组织分割的准确性;另外,通过引入目标组织亮度值作为期望的亮度值以及比较不同成像参数下的组织图像的亮度误差来确定目标成像参数,进而再用目标成像参数作为超声仪器的参数输出超声图像,实现了超声图像亮度的均匀化,操作简单,准确性高。
在一些示例性的实施方式中,所述对获取到的超声图像中各个所述像素点在预设方向上的梯度进行平滑处理,得到所述超声图像对应的目标梯度图像,包括:
针对每个预设方向,采用邻域均值平滑的方式,对所述预设方向上的各个所述像素点的梯度在水平方向和竖直方向依次进行平滑处理,得到所述预设方向上的平滑后的梯度图像;
对于所述超声图像的像素矩阵中每一个位置处的像素点,将各个所述预设方向上的平滑后的梯度图像中所述位置处的最大灰度值更新为所述像素点的灰度值;
确定各个更新后的所述像素点构成的图像为目标梯度图像。
上述实施例,在平滑处理的过程中,考虑多个预设方向,针对每个预设方向上的像素点的梯度均依次做水平和竖直两个方向上的平滑处理,与考虑一个方向相比,更适用于复杂的组织;另外,用每个预设方向上的最大灰度值作为目标梯度图像的相应位置处的像素点的灰度值,这样,应用目标梯度图像算出来的组织分割阈值更准确,组织分割更准确。
在一些示例性的实施方式中,所述基于所述目标梯度图像中各个像素点的灰度值、预设灰度阈值和所述超声图像中包含的组织匹配的组织占比常数确定组织分割阈值,包括:
确定目标像素点的个数与所述超声图像中包含的组织匹配的组织占比常数的乘积为组织占比阈值;其中,所述目标像素点为所述目标梯度图像中灰度值大于预设灰度阈值的像素点;
从大于所述预设灰度阈值的灰度值范围中选择一个第一目标灰度值作为组织分割阈值;其中所述目标梯度图像中,大于所述第一目标灰度值的其他灰度值对应的像素点的个数大于等于所述组织占比阈值,且大于第二目标灰度值的其他灰度值对应的像素点的个数小于所述组织占比阈值,所述第二目标灰度值与所述第一目标灰度值相邻,且大于所述第一目标灰度值。
上述实施例,在确定组织分割阈值时,不是简单的采用固定阈值分割的方式,而是考虑了不同部位的组织不同的特性,应用了超声图像中包含的组织匹配的组织占比常数,组织不同,该常数也不同;再根据大于等于组织占比阈值的像素点的个数以及小于组织占比阈值的像素点的个数来确定组织分割阈值,这样确定的组织分割阈值更准确。
在一些示例性的实施方式中,所述根据预先设置的每组成像参数和目标组织亮度值,分别确定每组成像参数对应的所述组织图像的亮度误差,包括:
针对任意一组成像参数,根据所述成像参数对所述组织图像中的每一个灰度值大于第一预设灰度值的像素点进行参数压缩,并将压缩后的值与目标组织亮度值进行做差并平方处理;
将各个像素点进行做差并平方处理后的结果加和,得到所述成像参数对应的所述组织图像的亮度误差。
上述实施例,通过计算不同成像参数下的组织图像的亮度误差,可以确定出较优的成像参数,应用该较优的成像参数得到的超声图像的亮度更均匀。而通过上述处理方式得到的组织图像的亮度误差更准确,进而确定的组织分割阈值更准确。
在一些示例性的实施方式中,所述应用所述组织分割阈值对所述超声图像进行组织分割,包括:
保持所述超声图像中的各个像素点中大于所述组织分割阈值的像素点的灰度值不变,并将小于所述组织分割阈值的像素点的灰度值设置为第二预设灰度值。
上述实施例,为了使组织区域更明显,在对超声图像进行组织分割时,将小于组织分割阈值的像素点的灰度值统一设置为第二预设灰度值,这样,使组织区域的边界更明显,进而组织分割更准确。
在一些示例性的实施方式中,通过如下方式获取超声图像:
检测用户的图像调节操作;
通过所述图像调节操作出发开启的第二线程获取在第一线程中输出的包络数据;其中,所述第一线程为按照预设模式输出预设超声图像时应用到的线程;
通过所述第二线程应用所述预设模式中的默认设置参数对所述包络数据进行参数压缩;
通过所述第二线程对压缩处理后的所述包络数据进行噪声矫正,得到超声图像。
上述实施例,开启两个线程,应用第二线程从第一线程去获取包括数据,并进行压缩与去噪等预处理等过程,而不是直接应用第一线程去实现成像参数的处理过程,避免了用户在参数处理前进行其他操作而影响成像效果。
在一些示例性的实施方式中,所述根据所述亮度误差选择一组成像参数为目标成像参数,包括:
在各组成像参数中,确定亮度误差最小时应用的成像参数为目标成像参数。
上述实施例,亮度误差越小,表明应用的成像参数的优化效果越好,因此,选用亮度误差最小时应用的成像参数作为目标成像参数,将该目标成像参数配置到超声仪器中,输出的超声图像的亮度更均匀。
在一些示例性的实施方式中,通过如下方式确定各个像素点在预设方向上的梯度:
根据与预设方向匹配的预设矩阵对超声图像中的各个像素点进行卷积处理,得到各个像素点在所述预设方向上的梯度。
上述实施例,应用卷积的方式确定各个像素点在各个预设方向上的梯度,操作简单,计算准确。
根据示例性的实施方式中的第二方面,提供一种超声仪器,该超声仪器包括处理器和数据获取单元,其中:
所述数据获取单元被配置为:
获取超声图像;其中,所述超声图像为包含人体器官或组织的图像;
所述处理器被配置为:
对获取到的超声图像中各个像素点在预设方向上的梯度进行平滑处理,得到所述超声图像对应的目标梯度图像;基于所述目标梯度图像中各个像素点的灰度值、预设灰度阈值和所述超声图像中包含的组织匹配的组织占比常数确定组织分割阈值;
应用所述组织分割阈值对所述超声图像进行组织分割,得到组织图像;
根据预先设置的每组成像参数和目标组织亮度值,分别确定每组成像参数对应的所述组织图像的亮度误差,根据所述亮度误差选择一组成像参数为目标成像参数;其中,每组成像参数包括动态范围与增益。
在一些示例性的实施方式中,所述处理器被配置为:
针对每个预设方向,采用邻域均值平滑的方式,对所述预设方向上的各个所述像素点的梯度在水平方向和竖直方向依次进行平滑处理,得到所述预设方向上的平滑后的梯度图像;
对于所述超声图像的像素矩阵中每一个位置处的像素点,将各个所述预设方向上的平滑后的梯度图像中所述位置处的最大灰度值更新为所述像素点的灰度值;
确定各个更新后的所述像素点构成的图像为目标梯度图像。
在一些示例性的实施方式中,所述处理器被配置为:
确定目标像素点的个数与所述超声图像中包含的组织匹配的组织占比常数的乘积为组织占比阈值;其中,所述目标像素点为所述目标梯度图像中灰度值大于预设灰度阈值的像素点;
从大于所述预设灰度阈值的灰度值范围中选择一个第一目标灰度值作为组织分割阈值;其中所述目标梯度图像中,大于所述第一目标灰度值的其他灰度值对应的像素点的个数大于等于所述组织占比阈值,且大于第二目标灰度值的其他灰度值对应的像素点的个数小于所述组织占比阈值,所述第二目标灰度值与所述第一目标灰度值相邻,且大于所述第一目标灰度值。
在一些示例性的实施方式中,所述处理器被配置为:
针对任意一组成像参数,根据所述成像参数对所述组织图像中的每一个灰度值大于第一预设灰度值的像素点进行参数压缩,并将压缩后的值与目标组织亮度值进行做差并平方处理;
将各个像素点进行做差并平方处理后的结果加和,得到所述成像参数对应的所述组织图像的亮度误差。
在一些示例性的实施方式中,所述处理器被配置为:
保持所述超声图像中的各个像素点中大于所述组织分割阈值的像素点的灰度值不变,并将小于所述组织分割阈值的像素点的灰度值设置为第二预设灰度值。
在一些示例性的实施方式中,所述处理器被配置为通过如下方式获取超声图像:
检测用户的图像调节操作;
通过所述图像调节操作出发开启的第二线程获取在第一线程中输出的包络数据;其中,所述第一线程为按照预设模式输出预设超声图像时应用到的线程;
通过所述第二线程应用所述预设模式中的默认设置参数对所述包络数据进行参数压缩;
通过所述第二线程对压缩处理后的所述包络数据进行噪声矫正,得到超声图像。
在一些示例性的实施方式中,所述处理器被配置为:
在各组成像参数中,确定亮度误差最小时应用的成像参数为目标成像参数。
在一些示例性的实施方式中,所述处理器被配置为通过如下方式确定各个像素点在预设方向上的梯度:
根据与预设方向匹配的预设矩阵对超声图像中的各个像素点进行卷积处理,得到各个像素点在所述预设方向上的梯度。
根据示例性的实施方式中的第三方面,提供一种超声图像的处理装置,该装置包括:
平滑处理模块,用于对获取到的超声图像中各个像素点在预设方向上的梯度进行平滑处理,得到所述超声图像对应的目标梯度图像;其中,所述超声图像为包含人体器官或组织的图像;
分割阈值确定模块,用于基于所述目标梯度图像中各个像素点的灰度值、预设灰度阈值和所述超声图像中包含的组织匹配的组织占比常数确定组织分割阈值;
组织分割模块,用于应用所述组织分割阈值对所述超声图像进行组织分割,得到组织图像;
参数确定模块,用于根据预先设置的每组成像参数和目标组织亮度值,分别确定每组成像参数对应的所述组织图像的亮度误差,根据所述亮度误差选择一组成像参数为目标成像参数;其中,每组成像参数包括动态范围与增益。
在一些示例性的实施方式中,所述平滑处理模块具体用于:
针对每个预设方向,采用邻域均值平滑的方式,对所述预设方向上的各个所述像素点的梯度在水平方向和竖直方向依次进行平滑处理,得到所述预设方向上的平滑后的梯度图像;
对于所述超声图像的像素矩阵中每一个位置处的像素点,将各个所述预设方向上的平滑后的梯度图像中所述位置处的最大灰度值更新为所述像素点的灰度值;
确定各个更新后的所述像素点构成的图像为目标梯度图像。
在一些示例性的实施方式中,所述分割阈值确定模块具体用于:
确定目标像素点的个数与所述超声图像中包含的组织匹配的组织占比常数的乘积为组织占比阈值;其中,所述目标像素点为所述目标梯度图像中灰度值大于预设灰度阈值的像素点;
从大于所述预设灰度阈值的灰度值范围中选择一个第一目标灰度值作为组织分割阈值;其中所述目标梯度图像中,大于所述第一目标灰度值的其他灰度值对应的像素点的个数大于等于所述组织占比阈值,且大于第二目标灰度值的其他灰度值对应的像素点的个数小于所述组织占比阈值,所述第二目标灰度值与所述第一目标灰度值相邻,且大于所述第一目标灰度值。
在一些示例性的实施方式中,所述参数确定模块具体用于:
针对任意一组成像参数,根据所述成像参数对所述组织图像中的每一个灰度值大于第一预设灰度值的像素点进行参数压缩,并将压缩后的值与目标组织亮度值进行做差并平方处理;
将各个像素点进行做差并平方处理后的结果加和,得到所述成像参数对应的所述组织图像的亮度误差。
在一些示例性的实施方式中,所述组织分割模块具体用于:
保持所述超声图像中的各个像素点中大于所述组织分割阈值的像素点的灰度值不变,并将小于所述组织分割阈值的像素点的灰度值设置为第二预设灰度值。
在一些示例性的实施方式中,还包括图像获取模块,用于通过如下方式获取超声图像:
检测用户的图像调节操作;
通过所述图像调节操作出发开启的第二线程获取在第一线程中输出的包络数据;其中,所述第一线程为按照预设模式输出预设超声图像时应用到的线程;
通过所述第二线程应用所述预设模式中的默认设置参数对所述包络数据进行参数压缩;
通过所述第二线程对压缩处理后的所述包络数据进行噪声矫正,得到超声图像。
在一些示例性的实施方式中,所述参数确定模块具体用于:
在各组成像参数中,确定亮度误差最小时应用的成像参数为目标成像参数。
在一些示例性的实施方式中,还包括梯度确定模块,用于通过如下方式确定各个像素点在预设方向上的梯度:
根据与预设方向匹配的预设矩阵对超声图像中的各个像素点进行卷积处理,得到各个像素点在所述预设方向上的梯度。
根据示例性的实施方式中的第四方面,提供一种计算机存储介质,所述计算机存储介质中存储有计算机程序指令,当所述指令在计算机上运行时,使得所述计算机执行如第一方面所述的超声图像的处理方法。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简要介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域的普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1示例性示出了本发明实施例提供的一种超声仪器的示意图;
图2示例性示出了本发明实施例提供的一种超声图像的处理方法的流程图;
图3示例性示出了本发明实施例提供的一种应用处理前的成像参数得到的超声图像的示意图;
图4示例性示出了本发明实施例提供的一种应用处理后的成像参数得到的超声图像的示意图;
图5示例性示出了本发明实施例提供的一种超声图像的处理方法的流程图;
图6示例性示出了本发明实施例提供的一种超声图像的处理装置的结构示意图;
图7示例性示出了本发明实施例提供的一种超声仪器的结构示意图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述。
附图中的任何元素数量均用于示例而非限制,以及任何命名都仅用于区分,而不具有任何限制含义。
超声波类型通常包括A型、B型、D型和M型,超声波类型不同,其应用场景通常有差别。本申请实施例中,对超声波的类型不具体限定,在实际的应用过程中的超声波类型可根据实际情况确定。
在具体的实践过程中,比如在医学领域,通常应用超声波来反应人体内部的结构,比如组织或者器官的状态。而增益调整是超声成像过程中的重要环节,通常情况下,超声***中为每个部位设置了最佳的增益(体现在图像亮度上)以及动态范围(体现在对比度上),而最佳的增益和动态范围适用于绝大部分患者,也即,图像有着良好的亮度和对比度,体现在图像中近场、中场、远场保持均匀一致。但是在实际的检测过程中(医生打图时),由于个体差异较大(比如声波能量衰减不同),会存在一部分人在默认设置下的图像均匀性较差,表现为图像亮暗不均,有的较亮、有的较暗,此时,为了图像的均匀性,需要医生手动调节增益和动态范围。比如在操作界面上显示一键操作按钮,用户可以根据观察喜好,调节不同档位来达到不同的对比度和亮度。
鉴于医生手动调节过程中复杂且效率低,相关技术中,大多数的超声***可以自动根据患者的超声图像进行自动增益调整,以达到图像亮度均匀化的目的。但是,相关技术中的做法主要有如下几种:采用求取结构张量的方式进行组织区域判断,通过对结构张量和构建函数值判定噪声和非噪声区域、组织和非组织区域,该种方法计算较为复杂,单纯通过特征值并不一定能完全判断出组织区域;简单采用固定阈值方式对图像分块,导致对组织的区域判定存在误差,不能最大程度的考虑所有的组织区域的亮度,与理想的组织亮度会存在误差。
为此,本申请提供了一种超声图像的处理方法,该方法中:对获取到的超声图像中各个像素点在预设方向上的梯度进行平滑处理,得到所述超声图像对应的目标梯度图像;基于所述目标梯度图像中各个像素点的灰度值、预设灰度阈值和所述超声图像中包含的组织匹配的组织占比常数确定组织分割阈值;应用所述组织分割阈值对所述超声图像进行组织分割,得到组织图像;根据预先设置的每组成像参数和目标组织亮度值,分别确定每组成像参数对应的所述组织图像的亮度误差,根据所述亮度误差选择一组成像参数为目标成像参数;其中,每组成像参数包括动态范围与增益。充分考虑不同人体部位的组织特性得到处理后的动态范围和增益,操作简单;应用处理后的动态范围和增益生成超声图像,实现了亮度的均匀化,操作简单、准确性高。
在介绍完本申请实施例的设计思想之后,下面对本申请实施例的技术方案能够适用的应用场景做一些简单介绍,需要说明的是,以下介绍的应用场景仅用于说明本申请实施例而非限定。在具体实施时,可以根据实际需要灵活地应用本申请实施例提供的技术方案。
参考图1,其为本申请实施例提供的一种超声仪器的示意图,该超声仪器为医用超声仪器(比如用户是医生),其中,包括成像显示区域11和操作台12,操作台中包括多个操作区域,每个操作区域可以包括多个操作按钮,图1中用三个操作按钮来示例。本申请实施例中,设置有一键操作按钮121,该一键操作按钮121可以触发本申请实施例中的超声图像的处理方法的执行,该过程中,实现成像参数的优化,进而应用处理后的动态范围和增益调节超声图像,以便将亮度均匀的图像展示给医生。
当然,本申请实施例提供的方法并不限用于图1所示的应用场景中,还可以用于其它可能的应用场景,本申请实施例并不进行限制。对于图1所示的应用场景的各个设备所能实现的功能将在后续的方法实施例中一并进行描述,在此先不过多赘述。
为进一步说明本申请实施例提供的技术方案,下面结合附图以及具体实施方式对此进行详细的说明。虽然本申请实施例提供了如下述实施例或附图所示的方法操作步骤,但基于常规或者无需创造性的劳动在方法中可以包括更多或者更少的操作步骤。在逻辑上不存在必要因果关系的步骤中,这些步骤的执行顺序不限于本申请实施例提供的执行顺序。
下面结合图1所示的应用场景,以医学领域中,对人体的器官或者组织的成像为例,对本申请实施例提供的技术方案进行说明。
参考图2,本申请实施例提供一种超声图像的处理方法,包括以下步骤:
S201、对获取到的超声图像中各个像素点在预设方向上的梯度进行平滑处理,得到所述超声图像对应的目标梯度图像;其中,超声图像为包含人体器官或组织的图像。
S202、基于所述目标梯度图像中各个像素点的灰度值、预设灰度阈值和所述超声图像中包含的组织匹配的组织占比常数确定组织分割阈值。
S203、应用所述组织分割阈值对所述超声图像进行组织分割,得到组织图像。
S204、根据预先设置的每组成像参数和目标组织亮度值,分别确定每组成像参数对应的所述组织图像的亮度误差,根据所述亮度误差选择一组成像参数为目标成像参数;其中,每组成像参数包括动态范围与增益。
本申请实施例,通过平滑处理可以去除获取到的超声图像各个预设方向上的噪声,进而增加组织分割的准确性;另外,在确定组织分割阈值时,考虑了不同部位的组织特性,应用与当前超声图像中包含的组织匹配的组织占比常数,这样得到的组织分割阈值更准确,进而也增加了组织分割的准确性;另外,通过引入目标组织亮度值作为期望的亮度值以及比较不同成像参数下的组织图像的亮度误差来确定目标成像参数,进而再用目标成像参数作为超声仪器的参数输出超声图像,实现了超声图像亮度的均匀化,操作简单,准确性高。
由于超声仪器中有预先配置的动态范围和增益,正常情况下,应用预先配置的动态范围和增益输出的超声图像的亮度和对比度都比较均匀,这样,可以直接应用该超声图像,该过程是在第一线程中完成的。而比如当输出的超声图像的亮度和对比度较差时,可以通过操作超声仪器上的一键操作按钮来实现对该超声图像的亮度和对比度的调节过程,而该调节过程即为本申请实施例的超声图像的处理方法的执行过程,这样,应用处理后的动态范围和增益,即可得到亮度和对比度比较均匀的超声图像。而在实际的应用过程中,该处理过程可以在第二线程中完成的,第二线程是在检测到用户的触发操作(比如医生操作了一键操作按钮)后开启的。
为了进行参数优化,首先获取到超声图像,再基于获取到的该超声图像来确定目标成像参数。具体的,检测用户的图像调节操作;通过图像调节操作出发开启的第二线程获取在第一线程中输出的包络数据;其中,第一线程为按照预设模式输出预设超声图像时应用到的线程;通过第二线程应用预设模式中的默认设置参数对包络数据进行参数压缩;通过第二线程对压缩处理后的包络数据进行噪声矫正,得到超声图像。
具体的,第一线程中输出的包络数据,是在第一线程中求包络后得到的,该数据是默认全局增益和TGC(Time Gain Compensation,时间增益补偿)下的包络数据,是不受用户增益调节后影响的。由于第二线程是在用户一键处理后出发的,因此,从第一线程去获取包络数据,可以避免用于在按下一键处理之前对全局增益和TGC进行调节造成的干扰成像的问题。
而后,第二线程会对包络数据进行参数压缩,比如是log压缩,在压缩时选用的是预设模式中的默认设置参数,这样避免一键处理之前对log压缩参数进行其他操作造成的干扰成像的问题。
另外,在对压缩后的数据进行噪声校正,得到超声图像。示例性的,噪声矫正过程中,比如,根据不同的探头类型以及获取的***位来生成噪声曲线,将噪声曲线预存于文件中,可以设置为固定的点数,不需要再次计算,这样可以节约时间。如果当深度不一样时可以通过插值获取,或者预存不同***位下不同深度方向的噪声曲线。
如上,噪声处理后得到的结果为超声图像,也即,本申请实施例超声图像的处理方法的对象。上述处理过程可以避免在一键处理前进行其他操作或者调节造成的干扰问题。
涉及到S201,当图像中存在边缘时,通常有较大的梯度值,相反,当图像中有比较平滑的部分时,灰度值变化较小,则相应的梯度也较小,由图像中的各个像素点的梯度构成的图像成为梯度图像。由于平滑操作是对梯度进行处理的,接下来首先对梯度的计算过程进行说明。
为了提高准确性,可以对超声图像的四个方向进行边缘检测,这四个方向比如可以是水平方向、竖直方向、左对角线方向和右对角线方向,这四个方向可以统称为预设方向,是基于超声图像的像素矩阵中的阵列式排列的像素点来说的。具体的,可以采用不同方向的sobel模板来进行边缘检测,实际应用的公式如下:
其中,I为输入图像,也即去噪后得到的超声图像,其为一个像素矩阵,矩阵中的各个元素为各个像素点的灰度值;*表示卷积;x表示水平方向,y表示竖直方向,xy表示左对角线方向,yx表示右对角线方向;Sx为水平方向的梯度输出结果,Sy为竖直方向的梯度输出结果,Sxy为左对角线方向的梯度输出结果,Syx为右对角线方向的梯度输出结果。其中,每个梯度输出结果均为像素矩阵的形式,元素为各个像素点的梯度。在得到各个像素点在预设方向上的梯度后,对各个梯度进行平滑处理,可以得到超声图像对应的目标梯度图像。具体的,目标梯度图像可以通过如下三个步骤得到:
第一步、针对每个预设方向,采用邻域均值平滑的方式,对预设方向上的各个像素点的梯度在水平方向和竖直方向依次进行平滑处理,得到预设方向上的平滑后的梯度图像。
分别对Sx、Sy、Sxy、Syx进行水平方向(x方向)平滑,以此来去除x方向噪声的影响,增加组织分割的准确性,采用邻域均值平滑的方式,通过如下公式实现:
其中,radius为滤波半径,Sx1(x,y)为x方向上对Sx的平滑结果,Sy1(x,y)为x方向上对Sy的平滑结果;Sxy1(x,y)为x方向上对Sxy的平滑结果;Syx1(x,y)为x方向上对Syx的平滑结果;i为x方向上像素点的序号。
基于在x方向上的平滑结果,分别进行y方向平滑,以此来去除y方向噪声的影响,增加组织分割的准确性,采用邻域均值平滑的方式,通过如下公式实现:
其中,radius为滤波半径,Sx2(x,y)为y方向上对Sx1的平滑结果,Sy2(x,y)为y方向上对Sy1的平滑结果;Sxy2(x,y)为y方向上对Sxy1的平滑结果;Syx2(x,y)为y方向上对Syx1的平滑结果;j为y方向上像素点的序号。
这样,每个像素点均有四个梯度,分别在4个预设方向。
第二步、对于超声图像的像素矩阵中每一个位置处的像素点,将各个预设方向上的平滑后的梯度图像中位置处的最大灰度值更新为像素点的灰度值。
在一个具体的例子中,比如对于第1行第1列的像素点,在平滑后的梯度结果中,取水平方向、竖直方向、左对角线方向和右对角线方向中的最大的灰度值为该点的灰度值;同理,对于其他像素点,也这样处理,这样,各个像素点的灰度值得到了更新。
第三步、确定各个更新后的像素点构成的图像为目标梯度图像。
由于各个像素点的灰度值用预设方向的梯度中的最大值进行更新后,各个更新后的像素点构成的图像为目标梯度图像。
涉及到S202,在得到目标梯度图像后,基于目标梯度图像中各个像素点的灰度值、预设灰度阈值和超声图像中包含的组织匹配的组织占比常数确定组织分割阈值。
其中,通常情况下,预设灰度阈值是根据当前组织类型确定的;超声图像中包含的组织匹配的组织占比常数是和组织相关的一个常数,比如组织是腹部组织,该常数为0.1,预设灰度阈值是250。
接下来说明,如何针对目标梯度图像中的各个像素点的灰度值,通过预设灰度阈值和超声图像中包含的组织匹配的组织占比常数来确定组织分割阈值。
具体的,将目标梯度图像中灰度值大于预设灰度阈值的像素点称为目标像素点,统计目标像素点的个数;并确定目标像素点的个数与超声图像中包含的组织匹配的组织占比常数的乘积为组织占比阈值。
在一个具体的例子中,目标梯度图像中,灰度值大于250的像素点一共有280个,超声图像中包含的组织匹配的组织占比常数为0.1,则确定的组织占比阈值为28。
在这个例子中,比如,灰度值为255的像素点的个数为9、灰度值为254的像素点的个数为10、灰度值为253的像素点的个数为12、灰度值为252的像素点的个数为7、灰度值为251的像素点的个数为9。这样,灰度值大于254的像素点的个数为9、灰度值大于253的像素点的个数为19、灰度值大于252的像素点的个数为26、灰度值大于251的像素点的个数为35、灰度值大于250的像素点的个数为48。
从大于预设灰度阈值的灰度值范围中选择一个第一目标灰度值作为组织分割阈值;其中目标梯度图像中,大于第一目标灰度值的其他灰度值对应的像素点的个数大于等于组织占比阈值,且大于第二目标灰度值的其他灰度值对应的像素点的个数小于组织占比阈值,第二目标灰度值与第一目标灰度值相邻,且大于第一目标灰度值。
在上述示例中,大于预设灰度值的灰度值范围为251-255,大于第一目标灰度值的其他灰度值对应的像素点的个数大于等于28,此时,第一目标灰度值可以取251和252;大于第二目标灰度值的其他灰度值对应的像素点的个数小于28,则第二目标灰度阈值可以取253、254和255;第二目标灰度值与第一目标灰度值相邻,且大于第一目标灰度值,则确定第一目标灰度值为251,也即,组织分割阈值为251。
在实际的应用过程中,组织分割阈值可以通过如下方式确定:
按照像素值,对目标梯度图像进行直方图统计,得到的直方图用Hist表示,Hist统计了0-255的每个灰度值的像素点的个数,对Hist进行累计统计,统计每个灰度值大于预设灰度阈值thr的个数总和,根据每个部位的组织占比阈值segthresh确定组织分割阈值Thresh。
segthresh=StructRatio*H(250);
Thresh=min(H>segthresh);
其中,n表示灰度值,k为全部像素点的编号;预设灰度阈值是用来剔除灰度值较小的部分,StructRatio为超声图像中包含的组织匹配的组织占比常数。
涉及到S203,在确定组织分割阈值后,应用该组织分割阈值对超声图像进行组织分割。具体的,保持超声图像中的各个像素点中大于组织分割阈值的像素点的灰度值不变,并将小于组织分割阈值的像素点的灰度值设置为第二预设灰度值。
在一个具体的例子中,为了提高组织图像分割的准确度,第二预设灰度值可以取0,;分割的图像是原始输入图像,也即,超声图像,通过如下方式实现:
其中,G(x,y)即为分割出的组织,称为组织图像;I(x,y)为输入的超声图像。以上完成了对组织区域的分割,该分割方法能从四个方向确定梯度最大值。另外,这里的分割为分块的意思,是将超声图像中的组织区域分离出来,而不是将超声图像分割成多个图像。
涉及到S204,在得到组织图像后,应用每组成像参数(动态范围和增益),均可以确定组织图像的亮度误差。具体的,通过如下方式计算组织图像的亮度误差:
假设一共有10组成像参数,针对任意一组成像参数(动态范围D1和增益G1),根据成像参数对组织图像中的每一个灰度值大于第一预设灰度值(例如可以取0)的像素点进行参数压缩,并将压缩后的值与目标组织亮度值进行做差并平方处理;将各个像素点进行做差并平方处理后的结果加和,得到成像参数对应的组织图像的亮度误差。
在一个具体的例子中,应用D*log10(data)+G公式,在该公式中,针对组织图像中的每一个灰度值大于0的像素点,将其灰度值表示为data作为上述公式的输入,D和G分别取当前组的成像参数的取值,得到的out即为当前像素点的亮度值;则组织图像的亮度误差为其中,ap为像素点p的亮度值,b为目标亮度值,P为像素点的总个数。同理,针对每组成像参数,均可以得到对应的组织亮度误差。再根据各个亮度误差的大小,从多组成像参数中选择一组成像参数作为目标成像参数。在实际的应用过程中,通常选用亮度误差最小时应用的成像参数为目标成像参数,其中,亮度误差越小,表明选用的成像参数得到亮度值越逼近期望的组织亮度值。
为了使本申请的技术方案更容易理解,下面用一个具体的示例说明组织分割的示意图,其中,图3示出了一种应用处理前的成像参数得到的超声图像的示意图;图4示出了一种应用处理后的成像参数得到的超声图像的示意图。通过图3和图4可以看出,应用处理后的成像参数得到的超声图像的亮度更均匀。
下面用一个完整的流程图对本申请实施例的超声图像的处理方法进行说明,参见图5。
S501、检测用户的图像调节操作。
S502、通过图像调节操作触发的第二线程获取在第一线程中输出的包络数据。
其中,第一线程为按照预设模式输出预设超声图像时应用到的线程。
S503、应用预设模式中的默认设置参数对包络数据进行参数压缩。
S504、对压缩处理后的包络数据进行噪声矫正,得到超声图像。
S505、根据与预设方向匹配的预设矩阵对超声图像中的各个像素点进行卷积处理,得到各个像素点在预设方向上的梯度。
S506、针对每个预设方向,采用邻域均值平滑的方式,对预设方向上的各个像素点的梯度在水平方向和竖直方向依次进行平滑处理,得到预设方向上的平滑后的梯度图像。
S507、对于超声图像的像素矩阵中每一个位置处的像素点,将各个预设方向上的平滑后的梯度图像中位置处的最大灰度值更新为像素点的灰度值。
S508、确定各个更新后的像素点构成的图像为目标梯度图像。
S509、确定目标像素点的个数与超声图像中包含的组织匹配的组织占比常数的乘积为组织占比阈值;其中,目标像素点为目标梯度图像中灰度值大于预设灰度阈值的像素点。
S510、从大于预设灰度阈值的灰度值范围中选择一个第一目标灰度值作为组织分割阈值。
其中目标梯度图像中,大于第一目标灰度值的其他灰度值对应的像素点的个数大于等于组织占比阈值,且大于第二目标灰度值的其他灰度值对应的像素点的个数小于组织占比阈值,第二目标灰度值与第一目标灰度值相邻,且大于第一目标灰度值;
S511、保持超声图像中的各个像素点中大于组织分割阈值的像素点的灰度值不变,并将小于组织分割阈值的像素点的灰度值设置为第二预设灰度值,进行组织分割,得到组织图像。
S512、针对任意一组成像参数,根据成像参数对组织图像中的每一个灰度值大于第一预设灰度值的像素点进行参数压缩,并将压缩后的值与目标组织亮度值进行做差并平方处理。
其中,每组成像参数包括动态范围与增益。
S513、将各个像素点进行做差并平方处理后的结果加和,得到成像参数对应的组织图像的亮度误差。
S514、在各组成像参数中,确定亮度值误差最小时应用的成像参数为目标成像参数。
该方法可以实现对噪声图像中组织区域的自动识别,并根据部位的不同进行组织的准确判断,通过将组织图像的亮度与目标亮度值的不断逼近实现将增益与期望增益的不断逼近达到理想动态范围及增益的目的,该方法计算简单,准确性高,能充分考虑到不同人体部位的组织特性。
如图6所示,基于相同的发明构思,本发明实施例提供一种超声图像的处理装置,包括平滑处理模块61、分割阈值确定模块62、组织分割模块63和参数确定模块64。
平滑处理模块61,用于对获取到的超声图像中各个像素点在预设方向上的梯度进行平滑处理,得到超声图像对应的目标梯度图像;其中,超声图像为包含人体器官或组织的图像;
分割阈值确定模块62,用于基于目标梯度图像中各个像素点的灰度值、预设灰度阈值和超声图像中包含的组织匹配的组织占比常数确定组织分割阈值;
组织分割模块63,用于应用组织分割阈值对超声图像进行组织分割,得到组织图像;
参数确定模块64,用于根据预先设置的每组成像参数和目标组织亮度值,分别确定每组成像参数对应的组织图像的亮度误差,根据亮度误差选择一组成像参数为目标成像参数;其中,每组成像参数包括动态范围与增益。
在一些示例性的实施方式中,平滑处理模块61具体用于:
针对每个预设方向,采用邻域均值平滑的方式,对预设方向上的各个像素点的梯度在水平方向和竖直方向依次进行平滑处理,得到预设方向上的平滑后的梯度图像;
对于超声图像的像素矩阵中每一个位置处的像素点,将各个预设方向上的平滑后的梯度图像中位置处的最大灰度值更新为像素点的灰度值;
确定各个更新后的像素点构成的图像为目标梯度图像。
在一些示例性的实施方式中,分割阈值确定模块62具体用于:
确定目标像素点的个数与超声图像中包含的组织匹配的组织占比常数的乘积为组织占比阈值;其中,目标像素点为目标梯度图像中灰度值大于预设灰度阈值的像素点;
从大于预设灰度阈值的灰度值范围中选择一个第一目标灰度值作为组织分割阈值;其中目标梯度图像中,大于第一目标灰度值的其他灰度值对应的像素点的个数大于等于组织占比阈值,且大于第二目标灰度值的其他灰度值对应的像素点的个数小于组织占比阈值,第二目标灰度值与第一目标灰度值相邻,且大于第一目标灰度值。
在一些示例性的实施方式中,参数确定模块64具体用于:
针对任意一组成像参数,根据成像参数对组织图像中的每一个灰度值大于第一预设灰度值的像素点进行参数压缩,并将压缩后的值与目标组织亮度值进行做差并平方处理;
将各个像素点进行做差并平方处理后的结果加和,得到成像参数对应的组织图像的亮度误差。
在一些示例性的实施方式中,组织分割模块63具体用于:
保持超声图像中的各个像素点中大于组织分割阈值的像素点的灰度值不变,并将小于组织分割阈值的像素点的灰度值设置为第二预设灰度值。
在一些示例性的实施方式中,还包括图像获取模块,用于通过如下方式获取超声图像:
检测用户的图像调节操作;
通过图像调节操作出发开启的第二线程获取在第一线程中输出的包络数据;其中,第一线程为按照预设模式输出预设超声图像时应用到的线程;
通过第二线程应用预设模式中的默认设置参数对包络数据进行参数压缩;
通过第二线程对压缩处理后的包络数据进行噪声矫正,得到超声图像。
在一些示例性的实施方式中,参数确定模块64具体用于:
在各组成像参数中,确定亮度误差最小时应用的成像参数为目标成像参数。
在一些示例性的实施方式中,还包括梯度确定模块,用于通过如下方式确定各个像素点在预设方向上的梯度:
根据与预设方向匹配的预设矩阵对超声图像中的各个像素点进行卷积处理,得到各个像素点在预设方向上的梯度。
由于该装置即是本发明实施例中的方法中的装置,并且该装置解决问题的原理与该方法相似,因此该装置的实施可以参见方法的实施,重复之处不再赘述。
如图7所示,基于相同的发明构思,本发明实施例提供一种超声仪器,该超声仪器包括:处理器701以及数据获取单元702。
数据获取单元702被配置为:
获取超声图像;
处理器701被配置为:
对获取到的超声图像中各个像素点在预设方向上的梯度进行平滑处理,得到超声图像对应的目标梯度图像;
基于目标梯度图像中各个像素点的灰度值、预设灰度阈值和超声图像中包含的组织匹配的组织占比常数确定组织分割阈值;
应用组织分割阈值对超声图像进行组织分割,得到组织图像;
根据预先设置的每组成像参数和目标组织亮度值,分别确定每组成像参数对应的组织图像的亮度误差,根据亮度误差选择一组成像参数为目标成像参数;其中,每组成像参数包括动态范围与增益。
一种可选的实施方式为,处理器701被配置为:
针对每个预设方向,采用邻域均值平滑的方式,对预设方向上的各个像素点的梯度在水平方向和竖直方向依次进行平滑处理,得到预设方向上的平滑后的梯度图像;
对于超声图像的像素矩阵中每一个位置处的像素点,将各个预设方向上的平滑后的梯度图像中位置处的最大灰度值更新为像素点的灰度值;
确定各个更新后的像素点构成的图像为目标梯度图像。
一种可选的实施方式为,处理器701被配置为:
确定目标像素点的个数与超声图像中包含的组织匹配的组织占比常数的乘积为组织占比阈值;其中,目标像素点为目标梯度图像中灰度值大于预设灰度阈值的像素点;
从大于预设灰度阈值的灰度值范围中选择一个第一目标灰度值作为组织分割阈值;其中目标梯度图像中,大于第一目标灰度值的其他灰度值对应的像素点的个数大于等于组织占比阈值,且大于第二目标灰度值的其他灰度值对应的像素点的个数小于组织占比阈值,第二目标灰度值与第一目标灰度值相邻,且大于第一目标灰度值。
一种可选的实施方式为,处理器701被配置为:
针对任意一组成像参数,根据成像参数对组织图像中的每一个灰度值大于第一预设灰度值的像素点进行参数压缩,并将压缩后的值与目标组织亮度值进行做差并平方处理;
将各个像素点进行做差并平方处理后的结果加和,得到成像参数对应的组织图像的亮度误差。
一种可选的实施方式为,处理器701被配置为:
保持超声图像中的各个像素点中大于组织分割阈值的像素点的灰度值不变,并将小于组织分割阈值的像素点的灰度值设置为第二预设灰度值。
一种可选的实施方式为,处理器701被配置为通过如下方式获取超声图像:
检测用户的图像调节操作;
通过图像调节操作出发开启的第二线程获取在第一线程中输出的包络数据;其中,第一线程为按照预设模式输出预设超声图像时应用到的线程;
通过第二线程应用预设模式中的默认设置参数对包络数据进行参数压缩;
通过第二线程对压缩处理后的包络数据进行噪声矫正,得到超声图像。
一种可选的实施方式为,处理器701被配置为:
在各组成像参数中,确定亮度误差最小时应用的成像参数为目标成像参数。
一种可选的实施方式为,处理器701被配置为通过如下方式确定各个像素点在预设方向上的梯度:
根据与预设方向匹配的预设矩阵对超声图像中的各个像素点进行卷积处理,得到各个像素点在预设方向上的梯度。
本发明实施例还提供一种计算机存储介质,计算机存储介质中存储有计算机程序指令,当指令在计算机上运行时,使得计算机执行上述电子家居设备的配网方法的步骤。
本领域内的技术人员应明白,本申请的实施例可提供为方法、***、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
显然,本领域的技术人员可以对本申请进行各种改动和变型而不脱离本申请的精神和范围。这样,倘若本申请的这些修改和变型属于本申请权利要求及其等同技术的范围之内,则本申请也意图包含这些改动和变型在内。
Claims (9)
1.一种超声仪器,其特征在于,包括处理器和数据获取单元:
所述数据获取单元被配置为:
获取超声图像;其中,所述超声图像为包含人体器官或组织的图像;
所述处理器被配置为:
对获取到的超声图像中各个像素点在预设方向上的梯度进行平滑处理,得到所述超声图像对应的目标梯度图像;
基于所述目标梯度图像中各个像素点的灰度值、预设灰度阈值和所述超声图像中包含的组织匹配的组织占比常数确定组织分割阈值;
应用所述组织分割阈值对所述超声图像进行组织分割,得到组织图像;
根据预先设置的每组成像参数和目标组织亮度值,分别确定每组成像参数对应的所述组织图像的亮度误差,根据所述亮度误差选择一组成像参数为目标成像参数;其中,每组成像参数包括动态范围与增益;
所述处理器具体被配置为:
确定目标像素点的个数与所述超声图像中包含的组织匹配的组织占比常数的乘积为组织占比阈值;其中,所述目标像素点为所述目标梯度图像中灰度值大于预设灰度阈值的像素点;
从大于所述预设灰度阈值的灰度值范围中选择一个第一目标灰度值作为组织分割阈值;其中所述目标梯度图像中,大于所述第一目标灰度值的其他灰度值对应的像素点的个数大于等于所述组织占比阈值,且大于第二目标灰度值的其他灰度值对应的像素点的个数小于所述组织占比阈值,所述第二目标灰度值与所述第一目标灰度值相邻,且大于所述第一目标灰度值。
2.根据权利要求1所述的超声仪器,其特征在于,所述处理器被配置为:
针对每个预设方向,采用邻域均值平滑的方式,对所述预设方向上的各个所述像素点的梯度在水平方向和竖直方向依次进行平滑处理,得到所述预设方向上的平滑后的梯度图像;
对于所述超声图像的像素矩阵中每一个位置处的像素点,将各个所述预设方向上的平滑后的梯度图像中所述位置处的最大灰度值更新为所述像素点的灰度值;
确定各个更新后的所述像素点构成的图像为目标梯度图像。
3.根据权利要求1所述的超声仪器,其特征在于,所述处理器被配置为:
针对任意一组成像参数,根据所述成像参数对所述组织图像中的每一个灰度值大于第一预设灰度值的像素点进行参数压缩,并将压缩后的值与目标组织亮度值进行做差并平方处理;
将各个像素点进行做差并平方处理后的结果加和,得到所述成像参数对应的所述组织图像的亮度误差。
4.根据权利要求1所述的超声仪器,其特征在于,所述处理器被配置为:
保持所述超声图像中的各个像素点中大于所述组织分割阈值的像素点的灰度值不变,并将小于所述组织分割阈值的像素点的灰度值设置为第二预设灰度值。
5.根据权利要求1所述的超声仪器,其特征在于,所述处理器被配置为通过如下方式获取超声图像:
检测用户的图像调节操作;
通过所述图像调节操作出发开启的第二线程获取在第一线程中输出的包络数据;其中,所述第一线程为按照预设模式输出预设超声图像时应用到的线程;
通过所述第二线程应用所述预设模式中的默认设置参数对所述包络数据进行参数压缩;
通过所述第二线程对压缩处理后的所述包络数据进行噪声矫正,得到超声图像。
6.根据权利要求1所述的超声仪器,其特征在于,所述处理器被配置为:
在各组成像参数中,确定亮度误差最小时应用的成像参数为目标成像参数。
7.根据权利要求1~6任一项所述的超声仪器,其特征在于,所述处理器被配置为通过如下方式确定各个像素点在预设方向上的梯度:
根据与预设方向匹配的预设矩阵对超声图像中的各个像素点进行卷积处理,得到各个像素点在所述预设方向上的梯度。
8.一种超声图像的处理方法,其特征在于,包括:
对获取到的超声图像中各个像素点在预设方向上的梯度进行平滑处理,得到所述超声图像对应的目标梯度图像;其中,所述超声图像为包含人体器官或组织的图像;
基于所述目标梯度图像中各个像素点的灰度值、预设灰度阈值和所述超声图像中包含的组织匹配的组织占比常数确定组织分割阈值;
应用所述组织分割阈值对所述超声图像进行组织分割,得到组织图像;
根据预先设置的每组成像参数和目标组织亮度值,分别确定每组成像参数对应的所述组织图像的亮度误差,根据所述亮度误差选择一组成像参数为目标成像参数;其中,每组成像参数包括动态范围与增益;
所述基于所述目标梯度图像中各个像素点的灰度值、预设灰度阈值和所述超声图像中包含的组织匹配的组织占比常数确定组织分割阈值,包括:
确定目标像素点的个数与所述超声图像中包含的组织匹配的组织占比常数的乘积为组织占比阈值;其中,所述目标像素点为所述目标梯度图像中灰度值大于预设灰度阈值的像素点;
从大于所述预设灰度阈值的灰度值范围中选择一个第一目标灰度值作为组织分割阈值;其中所述目标梯度图像中,大于所述第一目标灰度值的其他灰度值对应的像素点的个数大于等于所述组织占比阈值,且大于第二目标灰度值的其他灰度值对应的像素点的个数小于所述组织占比阈值,所述第二目标灰度值与所述第一目标灰度值相邻,且大于所述第一目标灰度值。
9.根据权利要求8所述的方法,其特征在于,所述对获取到的超声图像中各个所述像素点在预设方向上的梯度进行平滑处理,得到所述超声图像对应的目标梯度图像,包括:
针对每个预设方向,采用邻域均值平滑的方式,对所述预设方向上的各个所述像素点的梯度在水平方向和竖直方向依次进行平滑处理,得到所述预设方向上的平滑后的梯度图像;
对于所述超声图像的像素矩阵中每一个位置处的像素点,将各个所述预设方向上的平滑后的梯度图像中所述位置处的最大灰度值更新为所述像素点的灰度值;
确定各个更新后的所述像素点构成的图像为目标梯度图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111041389.0A CN113892967B (zh) | 2021-09-07 | 2021-09-07 | 一种超声图像的处理方法和超声仪器 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111041389.0A CN113892967B (zh) | 2021-09-07 | 2021-09-07 | 一种超声图像的处理方法和超声仪器 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113892967A CN113892967A (zh) | 2022-01-07 |
CN113892967B true CN113892967B (zh) | 2023-10-27 |
Family
ID=79188540
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111041389.0A Active CN113892967B (zh) | 2021-09-07 | 2021-09-07 | 一种超声图像的处理方法和超声仪器 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113892967B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114947951B (zh) * | 2022-08-01 | 2022-10-25 | 深圳华声医疗技术股份有限公司 | 超声成像控制方法、装置、超声设备及存储介质 |
CN115393230B (zh) * | 2022-10-28 | 2023-02-03 | 武汉楚精灵医疗科技有限公司 | 超声内镜图像标准化方法、装置及其相关装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1919144A (zh) * | 2005-08-26 | 2007-02-28 | 深圳迈瑞生物医疗电子股份有限公司 | 超声图像增强与斑点抑制方法 |
CN101849240A (zh) * | 2007-05-21 | 2010-09-29 | 康奈尔大学 | 图像物体分割方法 |
CN103845077A (zh) * | 2012-12-05 | 2014-06-11 | 深圳迈瑞生物医疗电子股份有限公司 | 超声图像增益优化方法及超声成像增益自动优化装置 |
CN104367331A (zh) * | 2013-08-15 | 2015-02-25 | 深圳市蓝韵实业有限公司 | 用于数字乳腺x射线机的全数字自动曝光方法 |
CN104574361A (zh) * | 2014-11-27 | 2015-04-29 | 沈阳东软医疗***有限公司 | 一种乳腺外周组织均衡的图像处理方法和装置 |
CN111415308A (zh) * | 2020-03-13 | 2020-07-14 | 青岛海信医疗设备股份有限公司 | 一种超声图像处理方法和通信终端 |
CN113040878A (zh) * | 2021-03-25 | 2021-06-29 | 青岛海信医疗设备股份有限公司 | 超声穿刺针的位置信息处理方法、超声设备及存储介质 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030048931A1 (en) * | 2001-03-23 | 2003-03-13 | Peter Johnson | Quantification and differentiation of tissue based upon quantitative image analysis |
JP5632680B2 (ja) * | 2010-08-25 | 2014-11-26 | 日立アロカメディカル株式会社 | 超音波画像処理装置 |
WO2013183432A1 (ja) * | 2012-06-07 | 2013-12-12 | 日立アロカメディカル株式会社 | 関心領域設定方法及び超音波診断装置 |
-
2021
- 2021-09-07 CN CN202111041389.0A patent/CN113892967B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1919144A (zh) * | 2005-08-26 | 2007-02-28 | 深圳迈瑞生物医疗电子股份有限公司 | 超声图像增强与斑点抑制方法 |
CN101849240A (zh) * | 2007-05-21 | 2010-09-29 | 康奈尔大学 | 图像物体分割方法 |
CN103845077A (zh) * | 2012-12-05 | 2014-06-11 | 深圳迈瑞生物医疗电子股份有限公司 | 超声图像增益优化方法及超声成像增益自动优化装置 |
CN104367331A (zh) * | 2013-08-15 | 2015-02-25 | 深圳市蓝韵实业有限公司 | 用于数字乳腺x射线机的全数字自动曝光方法 |
CN104574361A (zh) * | 2014-11-27 | 2015-04-29 | 沈阳东软医疗***有限公司 | 一种乳腺外周组织均衡的图像处理方法和装置 |
CN111415308A (zh) * | 2020-03-13 | 2020-07-14 | 青岛海信医疗设备股份有限公司 | 一种超声图像处理方法和通信终端 |
CN113040878A (zh) * | 2021-03-25 | 2021-06-29 | 青岛海信医疗设备股份有限公司 | 超声穿刺针的位置信息处理方法、超声设备及存储介质 |
Non-Patent Citations (3)
Title |
---|
SWI、DWI联合评价颈动脉斑块易损性的应用研究;周东晓等;《中国医学计算机成像杂志》;第23卷(第6期);第508-512页 * |
基于人工神经网络的肝部CT图像识别方法;肖海慧等;《信息技术与信息化》(第11期);第28-31页 * |
基于超声弹性成像的女性盆底肌肉定量评估方法;朱超超等;《中国生物医学工程学报》;第36卷(第4期);第401-409页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113892967A (zh) | 2022-01-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110505459B (zh) | 适用于内窥镜的图像颜色校正方法、装置和存储介质 | |
CN113892967B (zh) | 一种超声图像的处理方法和超声仪器 | |
US10813595B2 (en) | Fully automated image optimization based on automated organ recognition | |
JP6416582B2 (ja) | 医用画像における金属アーチファクト除去のための方法および装置 | |
EP2267655A2 (en) | Image processing method and image processing apparatus | |
US10217201B2 (en) | Image processing method, image processing system, and imaging system | |
JPH11155834A (ja) | 医用イメージング・システム及び医用画像を形成する方法 | |
JP4938427B2 (ja) | 脳出血部位体積計算装置 | |
CN109949274B (zh) | 一种图像处理方法、装置及*** | |
JP7034120B2 (ja) | X線装置を用い取得した画像データセットの自動修正方法、コンピュータプログラム、データメモリ及びx線装置 | |
US20110305382A1 (en) | Image processing method and radiographic apparatus using the same | |
CN109431584A (zh) | 超声成像的方法及*** | |
CN109330594B (zh) | 获取饱和带的参数的方法、波谱成像方法、***和介质 | |
WO2008056305A1 (en) | Noise reduction of an image signal | |
Ee et al. | Estimation of window width setting for CT scan brain images using mean of greyscale level to standard deviation ratio | |
US11253219B2 (en) | Image contrast enhancement of an x-ray image | |
JP6358244B2 (ja) | 医用画像処理装置、医用画像処理装置の制御方法、およびプログラム | |
JP6041957B1 (ja) | 超音波診断装置 | |
JP5937286B1 (ja) | 画像処理装置、画像処理方法、及び画像処理プログラム | |
KR101341576B1 (ko) | 등고선 기반 관심영역 결정방법 및 장치 | |
US20060062485A1 (en) | Method and system for filtering in medical imaging | |
CN116883270B (zh) | 一种碎石手术软镜清晰化成像*** | |
JP2006000449A (ja) | 超音波画像処理装置及び超音波画像処理方法 | |
US20240203565A1 (en) | Method for digital image processing of medical images recorded with a medical imaging device and complete medical system | |
CN117522722A (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 |