CN113040821B - 颅内脑血流三维成像去噪方法、装置、终端设备及存储介质 - Google Patents
颅内脑血流三维成像去噪方法、装置、终端设备及存储介质 Download PDFInfo
- Publication number
- CN113040821B CN113040821B CN202110229105.4A CN202110229105A CN113040821B CN 113040821 B CN113040821 B CN 113040821B CN 202110229105 A CN202110229105 A CN 202110229105A CN 113040821 B CN113040821 B CN 113040821B
- Authority
- CN
- China
- Prior art keywords
- dimensional imaging
- ultrasonic echo
- blood flow
- convex hull
- echo signals
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 145
- 238000007917 intracranial administration Methods 0.000 title claims abstract description 65
- 230000003727 cerebral blood flow Effects 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 35
- 230000017531 blood circulation Effects 0.000 claims abstract description 79
- 210000004204 blood vessel Anatomy 0.000 claims abstract description 48
- 239000000523 sample Substances 0.000 claims abstract description 25
- 238000001914 filtration Methods 0.000 claims description 15
- 238000002604 ultrasonography Methods 0.000 claims description 13
- 239000013589 supplement Substances 0.000 claims description 9
- 210000003625 skull Anatomy 0.000 claims description 7
- 230000002123 temporal effect Effects 0.000 claims description 6
- 230000001502 supplementing effect Effects 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims 1
- 230000002792 vascular Effects 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 9
- 238000002595 magnetic resonance imaging Methods 0.000 description 6
- 210000000709 aorta Anatomy 0.000 description 5
- 238000004590 computer program Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000008321 arterial blood flow Effects 0.000 description 2
- 210000004556 brain Anatomy 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 230000002980 postoperative effect Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 210000003743 erythrocyte Anatomy 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009469 supplementation Effects 0.000 description 1
- 210000001519 tissue Anatomy 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/06—Measuring blood flow
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/483—Diagnostic techniques involving the acquisition of a 3D volume of data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/488—Diagnostic techniques involving Doppler signals
-
- 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)
- Medical Informatics (AREA)
- Surgery (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Hematology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开了一种颅内脑血流三维成像去噪方法、装置、存储介质及终端设备，所述方法包括：利用经颅多普勒探头进行单波束或多波束超声波的颅内血管扫描，并接收超声回波信号，依据所接收到的超声回波信号计算得到血流信息；对所述血流信息进行数据处理，并构建具有所述血流信息的曲面，形成具有所述血流信息的三维成像数据；对所述三维成像数据进去噪处理，去除所述三维成像数据中的凸包噪声与细小血管噪声，输出去噪后的图像信息。本发明可对三维成像颅内血流三维成像数据进行去噪处理，提高数据精确度，以便得到更为准确的图像信息。
Description
技术领域
本发明涉及医学超声成像技术领域，尤其涉及一种颅内脑血流三维成像去噪方法、装置、终端设备及存储介质。
背景技术
目前针对患者脑部术后手术效果的评估，现有三维成像技术存在无法区分颅内血管主动脉与细小血管以及噪声过多的问题。虽然可以通过核磁共振(MRI)，然而MRI存在所需扫描时间较长、体内留有金属物品不宜接受MRI、危重病人不宜做以及成本较高(患者费用成本与医院采购设备成本)等缺点。现有技术中，在进行三维成像时，基本都是通过采集到的颅内血流信息来构建三维成像数据，但并未对其进行去噪处理，导致输出的三维成像数据存在干扰，影响数据分析的准确性。
因此，现有技术还有待改进和提高。
发明内容
本发明要解决的技术问题在于，针对现有技术的上述缺陷，提供一种颅内脑血流三维成像去噪方法、装置、终端设备及存储介质，旨在解决现有技术中未对三维成像数据进行去噪处理，导致输出的三维成像数据存在干扰，影响数据分析准确性的问题。
为了解决上述技术问题，本发明所采用的技术方案如下：
第一方面，本发明提供一种颅内脑血流三维成像去噪方法,其中，所述方法包括：
利用经颅多普勒探头进行颅内血管扫描，并接收超声回波信号，根据所接收到的超声回波信号得到血流信息；
对所述血流信息进行数据处理，并构建具有所述血流信息的三维成像数据；
根据所述三维成像数据，去除所述三维成像数据中的凸包噪声与细小血管噪声，并输出去噪后的图像信息。
在一种实现方式中，所述经颅多普勒探头设置在颅骨两侧颞窗上，以完成经颅骨两侧颞窗的颅内超声信号扫描。
在一种实现方式中，所述对所述血流信息进行数据处理，并构建具有所述血流信息的三维成像数据，包括：
获取所述血流信息中预设深度范围内的超声回波信号，并从所述预设范围内的超声回波信号中获取每一组处于同一深度信息的超声回波信号；
分别对每一组所述处于同一深度信息的超声回波信号进行插值处理，补充超声回波信号；
插值处理完成后，根据补充超声回波信号的血流信息，构建所述三维成像数据。
在一种实现方式中，所述获取所述血流信息中预设深度范围内的超声回波信号，包括：
获取所述血流信息中的深度信息，以及所述深度信息对应的超声回波信号；
根据所述深度信息，获取所述深度信息中预设深度范围的超声回波信号。
在一种实现方式中，所述分别对每一组所述处于同一深度信息的超声回波信号进行插值处理，补充超声回波信号，包括：
获取处于同一深度信息中相邻两个的超声回波信号的强度数据和方向数据；
计算相邻两个超声回波信号的强度数据与方向数据的均值，并将所述均值插入所述相邻两个超声回波信号之间。
在一种实现方式中，所述根据所述三维成像数据，去除所述三维成像数据中的凸包噪声与细小血管噪声，并输出去噪后的图像信息，包括：
根据所述三维成像数据，获取所述三维成像数据中的凸包以及所述凸包的尺寸；
若所述凸包尺寸小于预设值时，则确定所述凸包为所述凸包噪声，并将所述凸包噪声进行过滤。
在一种实现方式中，所述根据所述三维成像数据，去除所述三维成像数据中的凸包噪声与细小血管噪声，并输出去噪后的图像信息，包括：
在去除凸包噪声后，获取所述三维成像数据中的剩余凸包；
获取所述剩余凸包的平均能量值与剩余凸包的直径；
根据所述平均能量值与所述剩余凸包直径，确定所述剩余凸包中的细小血管，并将所述细小血管过滤。
第二方面，本实施例提供一种颅内脑血流三维成像去噪装置，其中，所述方法包括：
血流信息获取模块，用于利用经颅多普勒探头进行颅内血管扫描，并接收超声回波信号，根据所接收到的超声回波信号得到血流信息；
三维成像模块，用于对所述血流信息进行数据处理，并构建具有所述血流信息的三维成像数据；
数据去噪模块，用于根据所述三维成像数据，去除所述三维成像数据中的凸包噪声与细小血管噪声，并输出去噪后的图像信息，输出去噪后的图像信息。
第三方面，本发明实施例还提供一种终端设备，其中，所述终端设备包括存储器、处理器及存储在所述存储器中并可在所述处理器上运行的颅内脑血流三维成像去噪程序，所述处理器执行所述颅内脑血流三维成像去噪程序时，实现上述方案中任一项所述的颅内脑血流三维成像去噪方法的步骤。
第四方面，本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质上存储有颅内脑血流三维成像去噪程序，所述颅内脑血流三维成像去噪程序被处理器执行时，实现上述方案中任一项所述的颅内脑血流三维成像去噪方法的步骤。
有益效果：与现有技术相比，本发明提供了一种颅内脑血流三维成像去噪方法，本发明首先利用经颅多普勒探头进行颅内血管扫描，并接收超声回波信号，根据所接收到的超声回波信号得到血流信息；对所述血流信息进行数据处理，并构建具有所述血流信息的三维成像数据；根据所述三维成像数据，去除所述三维成像数据中的凸包噪声与细小血管噪声，并输出去噪后的图像信息。本发明可对三维成像颅内血流三维成像数据进行去噪处理，提高数据精确度，以便得到更为准确的图像信息。
附图说明
图1为本发明实施例提供的颅内脑血流三维成像去噪方法的具体实施方式的流程图。
图2为本发明实施例提供的颅内脑血流三维成像去噪方法中的超声探头的设置形式示意图。
图3为本发明实施例提供的颅内脑血流三维成像去噪方法中超声扫描形成的连续多个曲面示意图。
图4为本发明实施例提供的颅内脑血流三维成像去噪方法中进行插值处理的示意图。
图5为本发明实施例提供的颅内脑血流三维成像去噪方法中凸包噪声示意图。
图6为本发明实施例提供的颅内脑血流三维成像去噪方法中细小血管噪声示意图。
图7是本发明实施例提供的颅内脑血流三维成像去噪装置的原理框图。
图8是本发明实施例提供的终端设备的内部结构原理框图。
具体实施方式
为使本发明的目的、技术方案及效果更加清楚、明确，以下参照附图并举实施例对本发明进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
目前针对患者脑部术后手术效果的评估，现有三维成像技术存在无法区分颅内血管主动脉与细小血管以及噪声过多的问题。虽然可以通过核磁共振(MRI)，然而MRI存在所需扫描时间较长、体内留有金属物品不宜接受MRI、危重病人不宜做以及成本较高(患者费用成本与医院采购设备成本)等缺点。现有技术中，在进行三维成像时，基本都是通过采集到的颅内具有血流信息来构建三维成像数据，但并未对其进行去噪处理，导致医生在根据三维成像数据进行分析时容易出现误差。
为此，本实施例提供一种颅内脑血流三维成像去噪方法，本发明首先利用经颅多普勒探头进行颅内血管扫描，并接收超声回波信号，根据所接收到的超声回波信号得到血流信息；对所述血流信息进行数据处理，并构建具有所述血流信息的三维成像数据；根据所述三维成像数据，去除所述三维成像数据中的凸包噪声与细小血管噪声，并输出去噪后的图像信息。本实施例可对三维成像颅内血流三维成像数据进行去噪处理，提高数据精确度，以便医生可以根据去噪后的三维成像数据进行更为准确地判断。
具体实施时，如图1中所示，本实施例的颅内脑血流三维成像去噪方法包括如下步骤：
步骤S100、利用经颅多普勒探头进行颅内血管扫描，并接收超声回波信号，根据所接收到的超声回波信号得到血流信息。
具体实施时，超声在人体中传播会发生偏转，单振元或多振元超声探头在发射超声时，超声束在人体内的传播随着深度增加逐渐扩散。本实施例利用经颅多普勒探头进行颅内血管扫描，当单振元或多振元探头以极小的角度进行偏转扫描时，探头在连续轨迹上的移动中采集的回波信号近似一个扇形空间，具体如图2中所示。根据探头扫描的特性，决定了在靠近探头面部分，每束超声之间存在重叠，采集的血流信号存在冗余，而远离探头面部分，每个超声束之间存在间隙。当超声探头进行一次完整轨迹扫描时获得一个具有血流信息的三维锥体空间。本实施例所述经颅多普勒探头设置在颅骨两侧颞窗上，以完成经颅骨两侧颞窗的颅内超声信号扫描。
步骤S200、对所述血流信息进行数据处理，并构建具有所述血流信息的三维成像数据。
由于超声探头在进行扫描时，扫描的轨迹是一个扇形空间，并且该扇形空间中设置有不同的深度信息，每一个深度信息均进行了扫描。因此在每一个深度信息中均存在超声回波信号。具体地，本实施例为了实现三维成像数据后，本实施例将所述血流信息中的超声回波信号进行处理。由于在进行超声扫描时所呈现的扫描轨迹为弧形，并且形成一个扇形的扫描区域，在该扇形的扫描区域内随着扫描路径的延长以及改变，会出现相邻超声回波信号之间的距离过大，导致相邻两个超声回波信号之间存在未覆盖超声回波信号的区域。因此，为了使得构建的三维成像数据更为精确，本实施例对未覆盖超声回波信号的区域进行补充。具体实施时，本实施例并不是对所有的深度信息都需要进行超声回波信号补充，对于一些深度信息中，其中的超声回波信号是比较密集的，是不需要进行补充的。为此，本实施例首先需要获取预设深度范围内的超声回波信号，然后对该预设深度范围内的超声回波信号进行超声回波信号的补充。
在本实施例中，本实施例利用插值法对预设深度范围内的超声回波信号进行超声回波信号的补充。本实施例首先获取所述血流信息中的深度信息，以及所述深度信息对应的超声回波信号；根据所述深度信息，获取所述深度信息中预设深度范围的超声回波信号。具体地，本实施例首先获取超声回波信号的深度信息，并将所述深度信息与预设的第一阈值和第二阈值进行比较；将深度信息大于第一阈值以及深度信息小于第二阈值的超声回波信号舍弃。例如，当深度信息超过90mm后，超声回波信号通过插值得到的数据不能真实的表现颅内血流的情况，故超过90mm的超声回波信号会被舍弃。而超声回波信号在0～5mm间，由于存在皮肤、颅骨等组织，这些超声回波信号也将被舍弃；实际颅内脑血流三维空间分布建模使用的数据为5～90mm(即第一深度范围)之间的数据。而在5～90mm的深度信息中，有一些深度范围内的超声回波信号是密集的，无需进行补充。本实施例限定需要补充超声回波信号的深度范围为75～90mm的深度范围。具体地，在对血流信息进行分析时，实施例按照扫描轨迹点中从预设的深度范围75～90mm内获取处于同一深度的超声回波信号，并从所述预设范围内的超声回波信号中获取每一组处于同一深度信息的超声回波信号。由于本实施例在进行超声回波信号的插值时是对两两相邻的超声回波信号进行插值，因此本实施例分别对每一组所述处于同一深度信息的超声回波信号进行插值处理，补充超声回波信号，具体如图4中所示。在插值处理完成后，根据补充超声回波信号的血流信息，构建所述三维成像数据。
具体地，本实施例获取处于同一深度信息中相邻两个的超声回波信号的强度数据和方向数据；然后计算相邻两个超声回波信号的强度数据与方向数据的均值，并将所述均值插入所述相邻两个超声回波信号之间。进行数据插值时，在一个弧面上通过计算超声束所在深度范围(Gate)上超声回波信号的强度和方向量化此范围的超声表现，在两两相邻的超声束之间进行插值，所插入的数据计算方法为，Gate’＝(Gate1的强度和方向+Gate 2的强度和方向)/2；需要进行插值处理的深度限定于深度范围75～90mm之间，即5～74mm之间的超声回波数据无需进行插值处理。具体如图4中所示，图4中的
为实际测量得到的每个Gate深度范围)的强度值，/>
为两个实际Gate强度值之间插值得到的强度值。每个Gate的数据实际可认为是圆柱体内信号的血流信号的强度和方向的表示/>
从图中可以看出，在进行插值处理时，是计算相邻两个超声回波信号的强度数据和方向数据模的平均值，然后再将计算得到的平均值插入至相邻的两个超声回波信号之间。当插值结束后，本实施例根据插值完后的超声回波数据生成三维成像数据。
步骤S300、根据所述三维成像数据，去除所述三维成像数据中的凸包噪声与细小血管噪声，并输出去噪后的图像信息。
通过获取的具有血流信息的三维成像数据，构造空间中三维脑血流模型即可将这些信息表现出来。但如果直接将获取的三维图形数据进行三维重建显示出来的图像会发现其中包含了很多噪声及细小血管，较大程度的影响医生和研究人员对主动脉血管情况的判断。为此，本实施例需要对三维成像数据进行去噪处理，输出去噪后的图像信息，本实施例主要需要去除的噪声干扰为去除噪声干扰和去除细小血流信号的干扰。具体地，在完成三维脑血流数据中的空间位置量化后，得到三维空间中具有能量和方向信息的点。而这些具有空间信息的能量和点，聚集在一起即可认为一个凸包，如图5中所示。较强的血流信号在空间分布中具有连续性所形成的凸包通常体积较大，跨越多个深度区域，不同于包括固定强度阈值的白噪声和具有不连续特征的由细小血管形成的不连续强信号体积很小，在空间中离散分布。通过计算凸包尺寸的大小，当凸包尺寸小于预设值时，则可认为这个凸包为一个凸包噪声。将这些凸包噪声点的空间坐标记录后，在进行三维二次重建时进行有效的过滤，以便输出去噪后的图像信息。比如，当获取到所述三维成像数据中的凸包以及所述凸包的尺寸为2mm时，则将该凸包尺寸与预设值(比如2.5mm)进行比较，此时凸包尺寸小于所述预设值，因此就可以确定所述凸包为所述凸包噪声，并将所述凸包噪声进行过滤。
本实施例在去除凸包噪声后，获取所述三维成像数据中的剩余凸包，接着获取所述剩余凸包的平均能量值与剩余凸包的直径。然后根据所述平均能量值与所述剩余凸包直径，确定所述剩余凸包中的细小血管，并将所述细小血管过滤。在本实施例中，当将凸包噪声去除后，剩余凸包均为主要动脉血流信息在空间分布形成的凸包。细小血管中血流信息形成的凸包与主动脉之间的主要差别在两个方面：主动脉血管中红细胞的数量远大于细小血管，导致其超声信号能量更大，且主动脉血管直径大于细小血管直径，其间血流信息形成的空间分布具有不同深度的较大体积的连续分布。因此通过对比剩余凸包的平均能量和凸包直径，对能量小、体积小、空间分布离散的剩余凸包进行过滤。比如，当剩余凸包平均能量值以及剩余凸包的直径与预设的值进行比较，就可以确定出该剩余凸包是否是细小血管，如果是细小血管，则就对该细小血管进行过滤。通过以上三个维度的过滤后，则可获得需关注的主要动脉血流在空间的分布形态，具体如图6中所示。这样以便医生可以根据去噪后的三维成像数据进行更为准确地判断。
综上，本实施例首先利用经颅多普勒探头进行单波束或多波束超声波的颅内血管扫描，并接收超声回波信号，依据所接收到的超声回波信号计算得到血流信息；对所述血流信息进行数据处理，并构建具有所述血流信息的曲面，形成具有所述血流信息的三维成像数据；对所述三维成像数据进去噪处理，去除所述三维成像数据中的凸包噪声与细小血管噪声，输出去噪后的图像信息。本实施例可对三维成像颅内血流三维成像数据进行去噪处理，提高血流空间成像精确度，以便医生可以根据去噪后的三维成像数据进行更为准确地判断。
如图7中所示，本实施例还提供一种颅内脑血流三维成像去噪装置，该装置包括：血流信息获取模块10、三维成像模块20以及数据去噪模块30。具体地，所述血流信息获取模块10，用于利用经颅多普勒探头进行颅内血管扫描，并接收超声回波信号，根据所接收到的超声回波信号得到血流信息。所述三维成像模块20，用于对所述血流信息进行数据处理，并构建具有所述血流信息的三维成像数据。所述数据去噪模块30，用于根据所述三维成像数据，去除所述三维成像数据中的凸包噪声与细小血管噪声，并输出去噪后的图像信息。
在一种实现方式中，所述三维成像模块20包括：
超声回波信号获取单元，用于获取所述血流信息中预设深度范围内的超声回波信号，并从所述预设范围内的超声回波信号中获取每一组处于同一深度信息的超声回波信号；
超声回波信号补充单元，用于分别对每一组所述处于同一深度信息的超声回波信号进行插值处理，补充超声回波信号；
三维成像单元，用于插值处理完成后，根据补充超声回波信号的血流信息，构建所述三维成像数据。
在一种实现方式中，所述超声回波信号补充单元包括：
强度和方向获取子单元，用于获取处于同一深度信息中相邻两个的超声回波信号的强度数据和方向数据，
插值处理子单元，用于计算相邻两个超声回波信号的强度数据与方向数据的均值，并将所述均值插入所述相邻两个超声回波信号之间。
在一种实现方式中，所述数据去噪模块30包括：
凸包尺寸获取单元，用于根据所述三维成像数据，获取所述三维成像数据中的凸包以及所述凸包的尺寸；
凸包噪声过滤单元，用于若所述凸包尺寸小于预设值时，则确定所述凸包为所述凸包噪声，并将所述凸包噪声进行过滤。
在一种实现方式中，所述数据去噪模块30还包括：
剩余凸包获取单元，用于在去除凸包噪声后，获取所述三维成像数据中的剩余凸包；
能量值和直径获取单元，用于获取所述剩余凸包的平均能量值与剩余凸包的直径；
细小血管过滤单元，用于根据所述平均能量值与所述剩余凸包直径，确定所述剩余凸包中的细小血管，并将所述细小血管过滤。
基于上述实施例，本发明还提供了一种终端设备，其原理框图可以如图8所示。该终端设备包括通过系统总线连接的处理器、存储器、网络接口、显示屏、温度传感器。其中，该终端设备的处理器用于提供计算和控制能力。该终端设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该终端设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种颅内脑血流三维成像去噪方法。该终端设备的显示屏可以是液晶显示屏或者电子墨水显示屏，该终端设备的温度传感器是预先在终端设备内部设置，用于检测内部设备的运行温度。
本领域技术人员可以理解，图8中示出的原理框图，仅仅是与本发明方案相关的部分结构的框图，并不构成对本发明方案所应用于其上的终端设备的限定，具体的终端设备以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
在一个实施例中，提供了一种终端设备，终端设备包括存储器、处理器及存储在存储器中并可在处理器上运行的颅内脑血流三维成像去噪程序，处理器执行颅内脑血流三维成像去噪程序时，实现如下操作指令：
利用经颅多普勒探头进行颅内血管扫描，并接收超声回波信号，根据所接收到的超声回波信号得到血流信息；
对所述血流信息进行数据处理，并构建具有所述血流信息的三维成像数据；
根据所述三维成像数据，去除所述三维成像数据中的凸包噪声与细小血管噪声，并输出去噪后的图像信息。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本发明所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限，RAM以多种形式可得，诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
综上，本发明公开了一种颅内脑血流三维成像去噪方法、装置、存储介质及终端设备，所述方法包括：利用经颅多普勒探头进行颅内血管扫描，并接收超声回波信号，根据所接收到的超声回波信号得到血流信息；对所述血流信息进行数据处理，并构建具有所述血流信息的三维成像数据；根据所述三维成像数据，去除所述三维成像数据中的凸包噪声与细小血管噪声，并输出去噪后的图像信息。本发明可对三维成像颅内血流三维成像数据进行去噪处理，提高数据精确度，以便医生可以根据去噪后的三维成像数据进行更为准确地判断。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (5)
1.一种颅内脑血流三维成像去噪方法，其特征在于，所述方法包括：
利用经颅多普勒探头进行颅内血管扫描，并接收超声回波信号，根据所接收到的超声回波信号得到血流信息；
对所述血流信息进行数据处理，并构建具有所述血流信息的三维成像数据；
根据所述三维成像数据，去除所述三维成像数据中的凸包噪声与细小血管噪声，并输出去噪后的图像信息；
所述对所述血流信息进行数据处理，并构建具有所述血流信息的三维成像数据，包括：
获取所述血流信息中预设深度范围内的超声回波信号，并从所述预设深度范围内的超声回波信号中获取每一组处于同一深度信息的超声回波信号；
分别对每一组所述处于同一深度信息的超声回波信号进行插值处理，补充超声回波信号；
插值处理完成后，根据补充超声回波信号的血流信息，构建所述三维成像数据；
所述获取所述血流信息中预设深度范围内的超声回波信号，包括：
获取所述血流信息中的深度信息，以及所述深度信息对应的超声回波信号；
根据所述深度信息，获取所述深度信息中预设深度范围的超声回波信号；
所述根据所述深度信息，获取所述深度信息中预设深度范围的超声回波信号，包括：
获取超声回波信号的深度信息，并将所述深度信息与预设的第一阈值和第二阈值进行比较；
将深度信息大于第一阈值以及深度信息小于第二阈值的超声回波信号舍弃，得到所述预设深度范围的超声回波信号；
所述分别对每一组所述处于同一深度信息的超声回波信号进行插值处理，补充超声回波信号，包括：
获取处于同一深度信息中相邻两个的超声回波信号的强度数据和方向数据；
计算相邻两个超声回波信号的强度数据与方向数据的均值，并将所述均值插入所述相邻两个超声回波信号之间；
所述根据所述三维成像数据，去除所述三维成像数据中的凸包噪声与细小血管噪声，并输出去噪后的图像信息，包括：
根据所述三维成像数据，获取所述三维成像数据中的凸包以及所述凸包的尺寸；
若所述凸包尺寸小于预设值时，则确定所述凸包为所述凸包噪声，并将所述凸包噪声进行过滤；
所述根据所述三维成像数据，去除所述三维成像数据中的凸包噪声与细小血管噪声，并输出去噪后的图像信息，包括：
在去除凸包噪声后，获取所述三维成像数据中的剩余凸包；
获取所述剩余凸包的平均能量值与剩余凸包的直径；
根据所述平均能量值与所述剩余凸包直径，确定所述剩余凸包中的细小血管，并将所述细小血管过滤。
2.根据权利要求1所述的颅内脑血流三维成像去噪方法，其特征在于，所述经颅多普勒探头设置在颅骨两侧颞窗上，以完成经颅骨两侧颞窗的颅内超声信号扫描。
3.一种颅内脑血流三维成像去噪装置，其特征在于，所述装置包括：
血流信息获取模块，用于利用经颅多普勒探头进行颅内血管扫描，并接收超声回波信号，根据所接收到的超声回波信号得到血流信息；
三维成像模块，用于对所述血流信息进行数据处理，并构建具有所述血流信息的三维成像数据；
数据去噪模块，用于根据所述三维成像数据，去除所述三维成像数据中的凸包噪声与细小血管噪声，并输出去噪后的图像信息；
所述三维成像模块包括：
超声回波信号获取单元，用于获取所述血流信息中预设深度范围内的超声回波信号，并从所述预设深度范围内的超声回波信号中获取每一组处于同一深度信息的超声回波信号；
超声回波信号补充单元，用于分别对每一组所述处于同一深度信息的超声回波信号进行插值处理，补充超声回波信号；
三维成像单元，用于插值处理完成后，根据补充超声回波信号的血流信息，构建所述三维成像数据；
所述超声回波信号获取单元，包括：
获取所述血流信息中的深度信息，以及所述深度信息对应的超声回波信号；
根据所述深度信息，获取所述深度信息中预设深度范围的超声回波信号；
所述根据所述深度信息，获取所述深度信息中预设深度范围的超声回波信号，包括：
获取超声回波信号的深度信息，并将所述深度信息与预设的第一阈值和第二阈值进行比较；
将深度信息大于第一阈值以及深度信息小于第二阈值的超声回波信号舍弃，得到所述预设深度范围的超声回波信号；
所述超声回波信号补充单元包括：
强度和方向获取子单元，用于获取处于同一深度信息中相邻两个的超声回波信号的强度数据和方向数据，
插值处理子单元，用于计算相邻两个超声回波信号的强度数据与方向数据的均值，并将所述均值插入所述相邻两个超声回波信号之间；
所述数据去噪模包括：
凸包尺寸获取单元，用于根据所述三维成像数据，获取所述三维成像数据中的凸包以及所述凸包的尺寸；
凸包噪声过滤单元，用于若所述凸包尺寸小于预设值时，则确定所述凸包为所述凸包噪声，并将所述凸包噪声进行过滤；
所述数据去噪模块还包括：
剩余凸包获取单元，用于在去除凸包噪声后，获取所述三维成像数据中的剩余凸包；
能量值和直径获取单元，用于获取所述剩余凸包的平均能量值与剩余凸包的直径；
细小血管过滤单元，用于根据所述平均能量值与所述剩余凸包直径，确定所述剩余凸包中的细小血管，并将所述细小血管过滤。
4.一种终端设备，其特征在于，所述终端设备包括存储器、处理器及存储在所述存储器中并可在所述处理器上运行的颅内脑血流三维成像去噪程序，所述处理器执行所述颅内脑血流三维成像去噪程序时，实现如权利要求1-2任一项所述的颅内脑血流三维成像去噪方法的步骤。
5.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有颅内脑血流三维成像去噪程序，所述颅内脑血流三维成像去噪程序被处理器执行时，实现如权利要求1-2任一项所述的颅内脑血流三维成像去噪方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110229105.4A CN113040821B (zh) | 2021-03-02 | 2021-03-02 | 颅内脑血流三维成像去噪方法、装置、终端设备及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110229105.4A CN113040821B (zh) | 2021-03-02 | 2021-03-02 | 颅内脑血流三维成像去噪方法、装置、终端设备及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113040821A CN113040821A (zh) | 2021-06-29 |
| CN113040821B true CN113040821B (zh) | 2023-03-31 |
Family
ID=76509779
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110229105.4A Active CN113040821B (zh) | 2021-03-02 | 2021-03-02 | 颅内脑血流三维成像去噪方法、装置、终端设备及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113040821B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113616242A (zh) * | 2021-07-05 | 2021-11-09 | 青岛大学附属医院 | 一种超声经颅多普勒采集装置及系统 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105232086A (zh) * | 2015-10-29 | 2016-01-13 | 深圳市德力凯医疗设备股份有限公司 | 一种基于经颅多普勒的颅内血流三维信息显示方法及系统 |
| CN105816202B (zh) * | 2016-05-20 | 2018-08-21 | 深圳市德力凯医疗设备股份有限公司 | 一种颅内血流三维成像采集装置及方法 |
| CN108852415A (zh) * | 2018-05-07 | 2018-11-23 | 深圳市德力凯医疗设备股份有限公司 | 一种经颅三维脑血管复合成像方法及系统 |
| CN117731322A (zh) * | 2018-12-06 | 2024-03-22 | 深圳迈瑞生物医疗电子股份有限公司 | 超声成像方法、设备及可读存储介质 |
| CN112043308B (zh) * | 2020-08-31 | 2022-09-16 | 深圳市德力凯医疗设备股份有限公司 | 一种颅内三维脑血流重建方法、存储介质及超声设备 |
-
2021
- 2021-03-02 CN CN202110229105.4A patent/CN113040821B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN113040821A (zh) | 2021-06-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5605155A (en) | Ultrasound system for automatically measuring fetal head size | |
| US5795296A (en) | Pipeline process for automatically measuring object boundary from ultrasound image samples | |
| US9955951B2 (en) | Sensor coordinate calibration in an ultrasound system | |
| US10368833B2 (en) | Method and system for fetal visualization by computing and displaying an ultrasound measurement and graphical model | |
| CN110325119B (zh) | 卵巢卵泡计数和大小确定 | |
| EP1875867B1 (en) | Ultrasound diagnosis apparatus and method of processing ultrasound data | |
| EP2296011B1 (en) | Ultrasound system and method for providing multiple plane images for a plurality of views | |
| CN106659473B (zh) | 超声成像装置 | |
| JP7203823B2 (ja) | 画像とのタッチインタラクションを使用してボリュームデータから画像平面を抽出する超音波システム | |
| JPWO2010092919A1 (ja) | 医用画像表示方法、医用画像診断装置、及び医用画像表示装置 | |
| US9445783B2 (en) | Ultrasound diagnostic adapter, ultrasound diagnostic apparatus, and ultrasound diagnostic method | |
| CN113040821B (zh) | 颅内脑血流三维成像去噪方法、装置、终端设备及存储介质 | |
| EP3422953B1 (en) | Automated ultrasonic measurement of nuchal fold translucency | |
| KR20190022185A (ko) | 태아의 신체 계측 방법 및 이를 이용한 태아의 신체 계측 디바이스 | |
| US20200305837A1 (en) | System and method for guided ultrasound imaging | |
| CN115330681A (zh) | 一种超声图像测量方法、系统、智能终端及存储介质 | |
| CN112826535B (zh) | 一种超声成像中自动定位血管的方法和装置及设备 | |
| US20220280139A1 (en) | Ultrasound system acoustic output control using image data | |
| KR20130096051A (ko) | 탄성 데이터 생성 방법 및 장치, 및 이를 이용한 탄성 영상 생성 방법 및 장치 | |
| CN112638267B (zh) | 超声成像方法及系统、存储介质、处理器和计算机设备 | |
| CN114190972B (zh) | 超声波影像物件的体积计算方法及应用其的超声波系统 | |
| CN116616819B (zh) | 一种超声经颅多普勒血流分析系统、方法及存储介质 | |
| EP4226863A1 (en) | Fetal heart rate monitoring | |
| Yin et al. | 3D Single Vessel Fractional Moving Blood Volume (3D-svFMBV): Fully Automated Tissue Perfusion Estimation Using Ultrasound | |
| Xue et al. | Automated multivolume placental reconstruction using three‐dimensional power Doppler ultrasound and infrared camera tracking |
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 | ||
| CP03 | Change of name, title or address |
Address after: 518000, 18th Floor, Building B, High tech Innovation Center, Guangqiao Road, Tianliao Community, Yutang Street, Guangming District, Shenzhen, Guangdong Province Patentee after: Shenzhen Delikai Medical Electronics Co.,Ltd. Address before: 6 / F, building 10, Guanlong Village Second Industrial Zone, Xili street, Nanshan District, Shenzhen City, Guangdong Province Patentee before: SHENZHEN DELICA MEDICAL EQUIPMENT Co.,Ltd. |
|
| CP03 | Change of name, title or address |