CN109452954B - 超声成像方法及装置 - Google Patents
超声成像方法及装置 Download PDFInfo
- Publication number
- CN109452954B CN109452954B CN201710796923.6A CN201710796923A CN109452954B CN 109452954 B CN109452954 B CN 109452954B CN 201710796923 A CN201710796923 A CN 201710796923A CN 109452954 B CN109452954 B CN 109452954B
- Authority
- CN
- China
- Prior art keywords
- matched
- blood flow
- image information
- image
- target
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
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/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
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)
- Hematology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
公开了一种超声成像方法及装置,该方法包括:在当前帧所对应的超声图像的血流区域设定取样门;向所述取样门所在位置发射多普勒脉冲并接收回波信号,以得到第一频谱图像;检测所述血流区域的偏移量;根据所述血流区域的偏移量调整所述取样门所在位置;向调整后的取样门所在位置发射多普勒脉冲并接收回波信号,以得到第二频谱图像;从而采集到连续且更为准确的频谱图像。
Description
技术领域
本申请涉及医疗器械领域,特别涉及一种超声成像方法及装置。
背景技术
超声仪器一般用于医生观察人体的内部组织结构,医生将超声探头放在人体部位对应的皮肤表面,可以得到该部位的超声图像。超声由于其安全、方便、无损、廉价等特点,已经成为医生诊断的主要辅助手段之一。在超声诊断中,应用最广泛的技术包括灰阶成像技术、彩色血流成像技术、脉冲多普勒(Pulse Wave Doppler,PW)技术。
其中,PW技术在血流的定量分析和诊断中起到至关重要的作用。在超声仪器中,PW的操作方法为将PW的取样门置于待观测血流上(B模式或Color模式),然后启动PW扫查,系统将实时获取该取样门所在位置处血流的频谱信息。
图1为一种PW的工作模式:用户根据B图(B模式)或者C图(Color模式)将PW的取样门移到图示血流区域的待测血流处,启动PW扫描,系统将在该取样门所在位置处进行PW实时扫描,得到该位置PW频谱图。通常需要连续扫描若干秒才能获得满足临床要求的PW频谱图。但是,在扫查过程中,往往会因为探头移动、病人呼吸等因素导致血流位置发生改变,此时容易出现血流位置改变后采集不到连续PW图像的现象。实际临床中,往往需要病人的配合并采集多次才能获得符合诊断要求的连续PW图像。
发明内容
本申请提供了一种超声成像方法及装置,使得在超声扫描中,当血流位置发生改变时,能够采集到连续且更为准确的频谱图像。
本申请第一方面提供了一种超声成像方法,包括:
在当前帧所对应的超声图像的血流区域设定取样门;
向所述取样门所在位置发射多普勒脉冲并接收回波信号,以得到第一频谱图像;
检测所述血流区域的偏移量;
根据所述血流区域的偏移量调整所述取样门所在位置;
向调整后的取样门所在位置发射多普勒脉冲并接收回波信号,以得到第二频谱图像。
结合本申请的第一方面,在第一方面的第一种实现方式中,所述检测所述血流区域的偏移量包括:
获取待匹配帧所对应的超声图像,所述待匹配帧包括所述当前帧之前的至少一帧;
提取目标图像信息,所述目标图像信息为所述当前帧所对应的超声图像中的图像信息;
从所述待匹配帧中确定与所述目标图像信息的相似度位于第一预设阈值内的待匹配图像信息;
确定所述目标图像信息相对于所述待匹配图像信息的偏移量。
结合本申请的第一方面的第一种实现方式,在第一方面的第二种实现方式中,所述超声图像包括超声灰阶图像或者超声彩色血流图像,所述目标图像信息包括目标图像块,所述目标图像块包括所述取样门所在位置的图像块或者所述取样门所在位置的预设距离范围内的图像块。
结合本申请的第一方面的第二种实现方式,在第一方面的第三种实现方式中,所述从所述待匹配帧中确定出与所述目标图像信息的相似度位于第一预设阈值内的待匹配图像信息包括:
以所述目标图像块的中心为参考点,确定M×N个像素的搜索区域,其中,M、N均为自然数;
按照所述搜索区域从所述待匹配帧中提取待匹配图像块;
当所述待匹配图像块与所述目标图像块的相似度位于第一预设阈值内时,确定所述待匹配图像块为所述待匹配图像信息。
结合本申请的第一方面的第一种实现方式,在第一方面的第四种实现方式中,所述超声图像包括超声灰阶图像或者超声彩色血流图像,所述目标图像信息包括目标特征点,所述目标特征点为所述当前帧所对应的超声图像中的特征点。
结合本申请的第一方面的第四种实现方式,在第一方面的第五种实现方式中,所述当前帧所对应的超声图像中的特征点包括用于指示所述当前帧所对应的超声图像的梯度或者边界的特征点。
结合本申请的第一方面的第四种或第五种实现方式,在第一方面的第六种实现方式中,所述从所述待匹配帧中确定与所述目标图像信息的相似度位于第一预设阈值内的待匹配图像信息包括:
从所述待匹配帧所对应的超声图像中提取待匹配特征点;
以所述目标特征点为中心,确定M×N个像素的目标图像块,以及以所述待匹配特征点为中心,确定M×N个像素的待匹配图像块,其中,M、N均为自然数;
当所述待匹配图像块与所述目标图像块的相似度位于第一预设阈值内时,确定所述待匹配特征点为所述待匹配图像信息。
结合本申请的第一方面的第一种实现方式,在第一方面的第七种实现方式中,所述超声图像包括超声血流图像,所述目标图像信息包括目标血流边界轮廓,所述目标血流边界轮廓为所述当前帧所对应的超声血流图像的血流边界轮廓。
结合本申请的第一方面的第七种实现方式,在第一方面的第八种实现方式中,所述提取目标图像信息包括:
从当前帧所对应的超声血流图像中提取血流速度的绝对值超过第二预设阈值的血流点;
根据所述血流点提取所述当前帧所对应的超声血流图像的血流边界轮廓。
结合本申请的第一方面的第七种或第八种实现方式,在第一方面的第九种实现方式中,所述从所述待匹配帧中确定与所述目标图像信息的相似度位于第一预设阈值内的待匹配图像信息包括:
从所述待匹配帧所对应的超声图像中提取待匹配血流边界轮廓;
当所述待匹配血流边界轮廓与所述目标血流边界轮廓的相似度位于第一预设阈值内时,确定所述待匹配血流边界轮廓为所述待匹配图像信息。
本申请第二方面提供了一种超声成像装置,包括:
人机交互模块,用于在当前帧所对应的超声图像的血流区域设定取样门;
收发模块,用于向所述取样门所在位置发射多普勒脉冲并接收回波信号,以得到第一频谱图像;
处理模块,用于检测所述血流区域的偏移量,以及根据所述血流区域的偏移量调整所述取样门所在位置;
所述收发模块,还用于向调整后的取样门所在位置发射多普勒脉冲并接收回波信号,以得到第二频谱图像。
结合本申请的第二方面,在第二方面的第一种实现方式中,所述处理模块用于检测所述血流区域的偏移量包括:
所述处理模块,用于获取待匹配帧所对应的超声图像,所述待匹配帧包括所述当前帧之前的至少一帧;提取目标图像信息,所述目标图像信息为所述当前帧所对应的超声图像中的图像信息;从所述待匹配帧中确定与所述目标图像信息的相似度位于第一预设阈值内的待匹配图像信息;确定所述目标图像信息相对于所述待匹配图像信息的偏移量。
结合本申请的第二方面的第一种实现方式,在第二方面的第二种实现方式中,所述超声图像包括超声灰阶图像或者超声彩色血流图像,所述目标图像信息包括目标图像块,所述目标图像块包括所述取样门所在位置的图像块或者所述取样门所在位置的预设距离范围内的图像块。
结合本申请的第二方面的第二种实现方式,在第二方面的第三种实现方式中,所述处理模块用于从所述待匹配帧中确定出与所述目标图像信息的相似度位于第一预设阈值内的待匹配图像信息包括:
所述处理模块,用于以所述目标图像块的中心为参考点,确定M×N个像素的搜索区域;按照所述搜索区域从所述待匹配帧中提取待匹配图像块;以及当所述待匹配图像块与所述目标图像块的相似度位于第一预设阈值内时,确定所述待匹配图像块为所述待匹配图像信息;其中,M、N均为自然数。
结合本申请的第二方面的第一种实现方式,在第二方面的第四种实现方式中,所述超声图像包括超声灰阶图像或者超声彩色血流图像,所述目标图像信息包括目标特征点,所述目标特征点为所述当前帧所对应的超声图像中的特征点。
结合本申请的第二方面的第四种实现方式,在第二方面的第五种实现方式中,所述当前帧所对应的超声图像中的特征点包括用于指示所述当前帧所对应的超声图像的梯度或者边界的特征点。
结合本申请的第二方面的第四种或第五种实现方式,在第二方面的第六种实现方式中,所述处理模块用于从所述待匹配帧中确定与所述目标图像信息的相似度位于第一预设阈值内的待匹配图像信息包括:
所述处理模块,用于从所述待匹配帧所对应的超声图像中提取待匹配特征点;以所述目标特征点为中心,确定M×N个像素的目标图像块,以所述待匹配特征点为中心,确定M×N个像素的待匹配图像块;以及当所述待匹配图像块与所述目标图像块的相似度位于第一预设阈值内时,确定所述待匹配特征点为所述待匹配图像信息;其中,M、N均为自然数。
结合本申请的第二方面的第一种实现方式,在第二方面的第七种实现方式中,所述超声图像包括超声血流图像,所述目标图像信息包括目标血流边界轮廓,所述目标血流边界轮廓为所述当前帧所对应的超声血流图像的血流边界轮廓。
结合本申请的第二方面的第七种实现方式,在第二方面的第八种实现方式中,所述处理模块用于提取目标图像信息包括:
所述处理模块,用于从当前帧所对应的超声血流图像中提取血流速度的绝对值超过第二预设阈值的血流点;以及根据所述血流点提取所述当前帧所对应的超声血流图像的血流边界轮廓。
结合本申请的第二方面的第七种或第八种实现方式,在第二方面的第九种实现方式中,所述处理模块用于从所述待匹配帧中确定与所述目标图像信息的相似度位于第一预设阈值内的待匹配图像信息包括:
所述处理模块,用于从所述待匹配帧所对应的超声图像中提取待匹配血流边界轮廓;以及当所述待匹配血流边界轮廓与所述目标血流边界轮廓的相似度位于第一预设阈值内时,确定所述待匹配血流边界轮廓为所述待匹配图像信息。
本申请提供的技术方案,通过检测当前帧所对应的超声图像的血流区域的偏移量,并根据该偏移量调整取样门所在位置;使得在超声扫描中,当血流位置发生改变时,能够采集到连续且更为准确的频谱图像。
附图说明
图1为本申请提供的一个基于PW模式下的频谱图;
图2为本申请提供的一个超声成像系统结构图;
图3为本申请提供的一个超声成像方法的流程示意图;
图4为本申请提供的另一个基于PW模式下的频谱图;
图5为本申请提供的另一个基于PW模式下的频谱图;
图6为本申请提供的一个超声成像装置的结构图。
具体实施方式
如图2所示,为申请的提供的超声成像系统的结构示意图。该超声成像系统包括探头2、发射/接收选择开关3、发射电路4、接收电路5、波束合成模块6、处理模块7和显示器8。发射电路4将一组经过延迟聚焦的脉冲发送到探头2,探头2向受测机体组织(图2中未示出)发射超声波,经一定延时后接收从受测机体组织反射回来的带有组织信息的超声回波,并将此超声回波重新转换为电信号。接收电路5接收这些电信号,并将这些超声回波信号送入波束合成模块6。超声回波信号在波束合成模块6完成聚焦延时、加权和通道求和,再经过处理模块7进行处理,然后送入显示器8进行显示。
如3图所示,为本申请提供的超声成像方法的一个实施例流程示意图。该方法包括:
301、在当前帧所对应的超声图像的血流区域设定取样门。
需要说明的是,通过将超声探头置于被测对象的体表,可得到相应部位的超声图像。因此,该步骤中,在超声成像界面显示的当前帧所对应的超声图像中找出血流区域。为了定量计算血流速度,可在该血流区域内的某一位置处设定取样门(例如PW取样门)。该超声图像包括超声灰阶图像(即B图像)和/或超声彩色血流图像(C图像)。
302、向该取样门所在位置发射多普勒脉冲并接收回波信号,以得到第一频谱图像。
需要说明的是,通过步骤301将取样门置于待测血流的位置,然后启动PW扫查模式,该模式下,通过向被测对象发射多普勒脉冲并接收回波信号,经信号处理模块的处理,可以得到该待测血流的频谱图像。
303、检测该血流区域的偏移量。
该步骤通过实时检测血流区域所在位置的偏移量,从而确定取样门的偏移量。
可选的,检测该血流区域的偏移量包括:
获取待匹配帧所对应的超声图像,该待匹配帧包括当前帧之前的至少一帧;
提取目标图像信息,其中,该目标图像信息为当前帧所对应的超声图像中的图像信息;
从待匹配帧中确定与该目标图像信息的相似度位于第一预设阈值内的待匹配图像信息;
确定该目标图像信息相对于该待匹配图像信息的偏移量。
具体计算偏移量的方式可参考如下说明:
一方面,可以基于块匹配的方法利用当前帧和待匹配帧中B信号或血流信号的相关性计算取样门的偏移量。待匹配帧可以根据需要进行设定,可以是当前帧的上一帧或上若干帧的超声图像,也可以是起始帧的超声图像。
基于块匹配的方法,该目标图像信息包括目标图像块,该目标图像块包括该取样门所在位置的图像块或者该取样门所在位置的预设距离范围内的图像块。例如,取当前帧取样门或取样门附近1块或若干块血流或B图像块数据,块的大小可以根据需要设定,例如m*n,其中m、n均为自然数。
然后,以该目标图像块的中心为参考点,确定M×N个像素的搜索区域,其中,M、N均为自然数;该目标图像块为多个时,分别取每个图像块的中心位置为参考点,设定大小为M*N的区域为搜索区域。
最后,按照该搜索区域从该待匹配帧中提取与该目标图像块的相似度位于第一预设阈值内的待匹配图像块,计算该目标图像块相对于待匹配图像块的偏移量。
例如,对搜索区域内的每个图像块所选定的参考点,以该参考点为中心,取大小为m*n的图像块与目标图像块计算相似度,选取搜索区域内相似度最优的位置为目标图像块的最佳匹配块。目标图像块与其对应的最佳匹配块的中心位置之差(dx,dy)即为当前帧相对于待匹配帧的偏移量。
一种度量两个邻域块相似度的方法为计算两个邻域块对应像素差的绝对值之和,即
其中,IL,IR为从当前帧取出的图像块(即目标图像块)和从待匹配帧中取出的图像块各点的灰度值,Ω为图像块。从定义中可以看出,E值越小,说明形似度越好。
也可以采用其它方法进行相似度的计算,例如相关系数法,即:
此时,E值越大说明相似度越好。
相似度的计算包括但不局限于上述两种方法,其它类似的定义如欧式距离(Euclidean Distance),Cosine相似度(Cosine Similarity)等也可以达到相似的效果,此处不做详细描述。
采用取样门附近多个图像块、跟踪时也可采用B信号或血流信号(如速度、方差、能量等)中的一种或几种。因而,在计算取样门的偏移量时,会获得一组或多组偏移量。需要将多组偏移量合成最终的取样门的偏移量。合成的方法包括但不限于对多组偏移量取平均值、取中值、最小二乘估计进行拟合等等。
得到了取样门的偏移量,即也得到了取样门在当前帧的位置。将取样门自动移动到该位置并在该位置进行PW扫描。重复检测、调整、扫描的步骤,直至停止扫描,将可获得完整PW频谱图像。
值得注意的是,以上流程中,跟踪时也可采用B信号或血流信号(如速度、方差、能量等)中的一种或几种,其均属于本申请的保护范围。
另一方面,可基于特征点的匹配在当前帧和待匹配帧中B图或血流图像各自提取若干特征点,然后对两个图中的特征点进行匹配,找出其中一幅图中的特征点在另一个图中的对应点。
基于特征点的匹配,该目标图像信息包括目标特征点,该目标特征点为当前帧所对应的超声图像中的特征点。其中,当前帧所对应的超声图像中的特征点包括用于指示该当前帧所对应的超声图像的梯度或者边界的特征点。例如,对当前帧和待匹配帧进行特征点提取,特征点可以根据需要定义,例如是图像梯度较大的点、边界点等,或则对图像采用一些算子提取特征,例如Sobel算子、SIFT算子等,在此不再详述。
然后,以该目标特征点为中心,确定M×N个像素的目标图像块,以及以该待匹配特征点为中心,确定M×N个像素的待匹配图像块,其中,M、N均为自然数。
最后,从该待匹配帧中提取与该目标图像块的相似度位于第一预设阈值内的待匹配特征点,计算该目标特征点相对于待匹配特征点的偏移量。
需要说明的是,特征点匹配的目的是在当前帧和待匹配帧中找出相同位置的特征点。一种匹配方法为,对当前帧中的每个特征点(即目标特征点),取以该特征点为中心,大小为m*n的图像块,和待匹配帧中每个特征点,同样取特征点为中心,大小为m*n的图像块进行相似度匹配,取相似度最优(例如位于第一预设阈值内)的特征点为匹配点,相似度的定义如块匹配法中描述的类似。同时,还可根据相似度设定阈值,作为没有找到特征点匹配点的判断依据,例如,采用上述公式(1)作为相似度度量时,当计算得到的相似度大于某个阈值,则认为该特征点在待匹配图像没有对应点。
通过一系列特征点及其在另一个图中匹配点的位置即可得到一系列的偏移。
可将上述得到的所有偏移量通过取平均值、取中值、最小二乘估计进行拟合等方式得到最终的偏移量。也可取取样门附近轮廓的偏移量通过取平均值、取中值、最小二乘估计进行拟合等方式得到最终的偏移量。
得到了取样门的偏移量,即也得到了取样门在当前帧的位置。将取样门自动移动到该位置并在该位置进行PW扫描。重复检测、调整、扫描的步骤,直至停止扫描,将可获得完整PW频谱图像。
又一方面,可基于轮廓匹配的方法利用当前帧和待匹配帧中血流图像(可采用速度、方差或能量图像中的一个或多个)中血流的轮廓的相似性进行匹配。
基于轮廓匹配的方法,该目标图像信息包括目标血流边界轮廓,该目标血流边界轮廓为当前帧所对应的超声血流图像的血流边界轮廓。从当前帧所对应的超声血流图像中提取血流速度的绝对值超过第二预设阈值的血流点;根据该血流点提取当前帧所对应的超声血流图像的血流边界轮廓。例如,血流图像中,满足一定规则的点认为是血流,否则该点为非血流点。血流判断规则可根据实际需要设定,例如,速度的绝对值大于某个阈值的点认为是血流点,否则为非血流点。对当前帧和待匹配帧中的血流图像进行边界轮廓提取,获得每一块血流的边界轮廓。
然后,从该待匹配帧中提取与该目标血流边界轮廓的相似度位于第一预设阈值内的待匹配血流边界轮廓,计算该目标血流边界轮廓相对于待匹配血流边界轮廓的偏移量。
例如,将当前帧中的每条轮廓(即目标血流边界轮廓)和待匹配帧中的各条血流边界轮廓相比较,接受那些形状特征点差异小于某个阈值的轮廓,得到当前帧和待匹配帧一系列相匹配的血流边界轮廓,再以相匹配的血流边界轮廓的中心作为匹配点,得到一系列相对于当前帧的目标血流边界轮廓的偏移量。
形状特征点差异可以通过轮廓的一阶距或二阶距等来衡量,例如,一阶距为:
其中,(xi,yi)为轮廓坐标,n为轮廓点数,(xc,yc)为轮廓中心,从表达式可以看出,两个轮廓越接近,则一阶距的差异也会越小。
可将上述中得到的所有偏移量通过取平均值、取中值、最小二乘估计进行拟合等方式得到最终的偏移量。也可取取样门附近轮廓的偏移量通过取平均值、取中值、最小二乘估计进行拟合等方式得到最终的偏移量。
得到了取样门的偏移量,即也得到了取样门在当前帧的位置。将取样门自动移动到该位置并在该位置进行PW扫描。重复检测、调整、扫描的步骤,直至停止扫描,将可获得完整PW频谱图像。
304、根据所述血流区域的偏移量调整所述取样门所在位置。
需要说明的是,基于步骤303中血流区域的位置偏移量的检测结果,可实时调整取样门的位置,例如按照步骤303中血流区域的位置偏移量,相应的调整取样门(例如PW取样门)的位置,从而实现取样门对血流区域位置的实时跟踪,即随着血流区域的位置的变化,取样门也相应的改变位置,以得到连续且较为准确的频谱图像(如PW图像)。
305、向调整后的取样门所在位置发射多普勒脉冲并接收回波信号,以得到第二频谱图像。
需要说明的是,取样门的位置随血流区域的位置的变化做出了相应调整后,基于步骤302得到的频谱图像,继续向该取样门所在的新位置发射多普勒脉冲并接收回波信号,从而得到连续的且较为准确的频谱图像(如PW图像)。例如,如图1所示,为某一时刻血流区域及取样门的位置,由于探头移动、病人呼吸等因素的影响,血流区域发生了位移(如图4所示),传统PW无法获得连续的PW图像,基于本申请提供的超声成像方法,PW取样门随着血流位置的改变而改变,得到连续的PW图像(如图5所示)。
图6为本申请提供的一个超声成像装置的结构示意图。该超声成像装置包括:收发模块602、处理模块604。该收发模块602集成了图2所示发射/接收选择开关3、发射电路4、接收电路5的功能,该处理模块604集成了图2所示波束合成模块6、处理模块7的功能。
其中,收发模块602,用于在当前帧所对应的超声图像的血流区域设定取样门后,向该取样门所在位置发射多普勒脉冲并接收回波信号,以得到第一频谱图像;
处理模块604,用于检测该血流区域的偏移量,以及根据该血流区域的偏移量调整该取样门所在位置;
该收发模块602,还用于向调整后的取样门所在位置发射多普勒脉冲并接收回波信号,以得到第二频谱图像。
可选的,该处理模块604用于检测该血流区域的偏移量包括:
该处理模块604,用于获取待匹配帧所对应的超声图像,该待匹配帧包括该当前帧之前的至少一帧;提取目标图像信息,该目标图像信息为该当前帧所对应的超声图像中的图像信息;从该待匹配帧中确定与该目标图像信息的相似度位于第一预设阈值内的待匹配图像信息;确定该目标图像信息相对于该待匹配图像信息的偏移量。
可选的,该超声图像包括超声灰阶图像或者超声彩色血流图像,该目标图像信息包括目标图像块,该目标图像块包括该取样门所在位置的图像块或者该取样门所在位置的预设距离范围内的图像块。
可选的,该处理模块604用于从该待匹配帧中确定出与该目标图像信息的相似度位于第一预设阈值内的待匹配图像信息包括:
该处理模块604,用于以该目标图像块的中心为参考点,确定M×N个像素的搜索区域;按照该搜索区域从该待匹配帧中提取待匹配图像块;以及当该待匹配图像块与该目标图像块的相似度位于第一预设阈值内时,确定该待匹配图像块为该待匹配图像信息;其中,M、N均为自然数。
可选的,该超声图像包括超声灰阶图像或者超声彩色血流图像,该目标图像信息包括目标特征点,该目标特征点为该当前帧所对应的超声图像中的特征点。
可选的,该当前帧所对应的超声图像中的特征点包括用于指示该当前帧所对应的超声图像的梯度或者边界的特征点。
可选的,该处理模块604用于从该待匹配帧中确定与该目标图像信息的相似度位于第一预设阈值内的待匹配图像信息包括:
该处理模块604,用于从该待匹配帧所对应的超声图像中提取待匹配特征点;以该目标特征点为中心,确定M×N个像素的目标图像块,以该待匹配特征点为中心,确定M×N个像素的待匹配图像块;以及当该待匹配图像块与该目标图像块的相似度位于第一预设阈值内时,确定该待匹配特征点为该待匹配图像信息;其中,M、N均为自然数。
可选的,该超声图像包括超声血流图像,该目标图像信息包括目标血流边界轮廓,该目标血流边界轮廓为该当前帧所对应的超声血流图像的血流边界轮廓。
可选的,该处理模块604用于提取目标图像信息包括:
该处理模块604,用于从当前帧所对应的超声血流图像中提取血流速度的绝对值超过第二预设阈值的血流点;以及根据该血流点提取该当前帧所对应的超声血流图像的血流边界轮廓。
可选的,该处理模块604用于从该待匹配帧中确定与该目标图像信息的相似度位于第一预设阈值内的待匹配图像信息包括:
该处理模块604,用于从该待匹配帧所对应的超声图像中提取待匹配血流边界轮廓;以及当该待匹配血流边界轮廓与该目标血流边界轮廓的相似度位于第一预设阈值内时,确定该待匹配血流边界轮廓为该待匹配图像信息。
上述超声成像装置的相关描述及效果可对应参考方法进行理解,在此不再赘述。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,如果可以通过计算机程序来指令相关的硬件来完成,那么相应的计算机程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory,ROM)或随机存储记忆体(Random Access Memory,RAM)等。
以上内容是结合具体的优选实施方式对本申请所作的进一步详细说明,不能认定本申请的具体实施只局限于这些说明。对于本申请所属技术领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本申请的保护范围。
Claims (12)
1.一种超声成像方法,其特征在于,包括:
在当前帧所对应的超声图像的血流区域设定取样门;
向所述取样门所在位置发射多普勒脉冲并接收回波信号,以得到第一频谱图像;
检测所述血流区域的偏移量;
根据所述血流区域的偏移量调整所述取样门所在位置;
向调整后的取样门所在位置发射多普勒脉冲并接收回波信号,以得到第二频谱图像;
所述检测所述血流区域的偏移量包括:
获取待匹配帧所对应的超声图像,所述待匹配帧包括所述当前帧之前的多帧;
提取目标图像信息,所述目标图像信息为所述当前帧所对应的超声图像中的图像信息,所述目标图像信息包括目标血流边界轮廓,所述提取目标图像信息包括:从当前帧所对应的超声血流图像中提取血流速度的绝对值超过第二预设阈值的血流点;根据所述血流点提取所述当前帧所对应的超声血流图像的血流边界轮廓;
从所述待匹配帧中确定与所述目标图像信息的相似度位于第一预设阈值内的待匹配图像信息,所述从所述待匹配帧中确定与所述目标图像信息的相似度位于第一预设阈值内的待匹配图像信息包括:从所述待匹配帧所对应的多张超声图像中分别提取待匹配血流边界轮廓;当所述待匹配血流边界轮廓与所述目标血流边界轮廓的相似度位于第一预设阈值内时,确定所述待匹配血流边界轮廓为所述待匹配图像信息,所述待匹配图像信息包括与所述目标血流边界轮廓的相似度位于第一预设阈值内的一系列相匹配的血流边界轮廓,所述相似度根据所述待匹配血流边界轮廓与所述目标血流边界轮廓的一阶距或者二阶距确定;
确定所述目标图像信息相对于所述待匹配图像信息的偏移量,所述确定所述目标图像信息相对于所述待匹配图像信息的偏移量包括:分别确定所述一系列相匹配的血流边界轮廓相对于当前帧的目标血流边界轮廓的偏移量,对得到的多个偏移量进行取平均值、取中值或者最小二乘估计进行拟合确定所述目标图像信息相对于所述待匹配图像信息的偏移量。
2.根据权利要求1所述的方法,其特征在于,所述超声图像包括超声灰阶图像或者超声彩色血流图像,所述目标图像信息包括目标图像块,所述目标图像块包括所述取样门所在位置的图像块或者所述取样门所在位置的预设距离范围内的图像块。
3.根据权利要求2所述的方法,其特征在于,所述从所述待匹配帧中确定出与所述目标图像信息的相似度位于第一预设阈值内的待匹配图像信息包括:
以所述目标图像块的中心为参考点,确定M×N个像素的搜索区域,其中,M、N均为自然数;
按照所述搜索区域从所述待匹配帧中提取待匹配图像块;
当所述待匹配图像块与所述目标图像块的相似度位于第一预设阈值内时,确定所述待匹配图像块为所述待匹配图像信息。
4.根据权利要求1所述的方法,其特征在于,所述超声图像包括超声灰阶图像或者超声彩色血流图像,所述目标图像信息包括目标特征点,所述目标特征点为所述当前帧所对应的超声图像中的特征点。
5.根据权利要求4所述的方法,其特征在于,所述当前帧所对应的超声图像中的特征点包括用于指示所述当前帧所对应的超声图像的梯度或者边界的特征点。
6.根据权利要求4或5所述的方法,其特征在于,所述从所述待匹配帧中确定与所述目标图像信息的相似度位于第一预设阈值内的待匹配图像信息包括:
从所述待匹配帧所对应的超声图像中提取待匹配特征点;
以所述目标特征点为中心,确定M×N个像素的目标图像块,以及以所述待匹配特征点为中心,确定M×N个像素的待匹配图像块,其中,M、N均为自然数;
当所述待匹配图像块与所述目标图像块的相似度位于第一预设阈值内时,确定所述待匹配特征点为所述待匹配图像信息。
7.一种超声成像装置,其特征在于,包括:
收发模块,用于在当前帧所对应的超声图像的血流区域设定取样门后,向所述取样门所在位置发射多普勒脉冲并接收回波信号,以得到第一频谱图像;
处理模块,用于检测所述血流区域的偏移量,以及根据所述血流区域的偏移量调整所述取样门所在位置;
所述收发模块,还用于向调整后的取样门所在位置发射多普勒脉冲并接收回波信号,以得到第二频谱图像;
所述处理模块用于检测所述血流区域的偏移量包括:
所述处理模块,用于获取待匹配帧所对应的超声图像,所述待匹配帧包括所述当前帧之前的多帧;提取目标图像信息,所述目标图像信息为所述当前帧所对应的超声图像中的图像信息,所述目标图像信息包括目标血流边界轮廓,所述提取目标图像信息包括从当前帧所对应的超声血流图像中提取血流速度的绝对值超过第二预设阈值的血流点;以及根据所述血流点提取所述当前帧所对应的超声血流图像的血流边界轮廓;从所述待匹配帧中确定与所述目标图像信息的相似度位于第一预设阈值内的待匹配图像信息,所述从所述待匹配帧中确定与所述目标图像信息的相似度位于第一预设阈值内的待匹配图像信息包括:从所述待匹配帧所对应的超声图像中提取待匹配血流边界轮廓;以及当所述待匹配血流边界轮廓与所述目标血流边界轮廓的相似度位于第一预设阈值内时,确定所述待匹配血流边界轮廓为所述待匹配图像信息,所述待匹配图像信息包括与所述目标血流边界轮廓的相似度位于第一预设阈值内的一系列相匹配的血流边界轮廓,所述相似度根据所述待匹配血流边界轮廓与所述目标血流边界轮廓的一阶距或者二阶距确定;确定所述目标图像信息相对于所述待匹配图像信息的偏移量,所述确定所述目标图像信息相对于所述待匹配图像信息的偏移量包括:分别确定所述一系列相匹配的血流边界轮廓相对于当前帧的目标血流边界轮廓的偏移量,对得到的多个偏移量进行取平均值、取中值或者最小二乘估计进行拟合确定所述目标图像信息相对于所述待匹配图像信息的偏移量。
8.根据权利要求7所述的装置,其特征在于,所述超声图像包括超声灰阶图像或者超声彩色血流图像,所述目标图像信息包括目标图像块,所述目标图像块包括所述取样门所在位置的图像块或者所述取样门所在位置的预设距离范围内的图像块。
9.根据权利要求8所述的装置,其特征在于,所述处理模块用于从所述待匹配帧中确定出与所述目标图像信息的相似度位于第一预设阈值内的待匹配图像信息包括:
所述处理模块,用于以所述目标图像块的中心为参考点,确定M×N个像素的搜索区域;按照所述搜索区域从所述待匹配帧中提取待匹配图像块;以及当所述待匹配图像块与所述目标图像块的相似度位于第一预设阈值内时,确定所述待匹配图像块为所述待匹配图像信息;其中,M、N均为自然数。
10.根据权利要求7所述的装置,其特征在于,所述超声图像包括超声灰阶图像或者超声彩色血流图像,所述目标图像信息包括目标特征点,所述目标特征点为所述当前帧所对应的超声图像中的特征点。
11.根据权利要求10所述的装置,其特征在于,所述当前帧所对应的超声图像中的特征点包括用于指示所述当前帧所对应的超声图像的梯度或者边界的特征点。
12.根据权利要求10或11所述的装置,其特征在于,所述处理模块用于从所述待匹配帧中确定与所述目标图像信息的相似度位于第一预设阈值内的待匹配图像信息包括:
所述处理模块,用于从所述待匹配帧所对应的超声图像中提取待匹配特征点;以所述目标特征点为中心,确定M×N个像素的目标图像块,以所述待匹配特征点为中心,确定M×N个像素的待匹配图像块;以及当所述待匹配图像块与所述目标图像块的相似度位于第一预设阈值内时,确定所述待匹配特征点为所述待匹配图像信息;其中,M、N均为自然数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710796923.6A CN109452954B (zh) | 2017-09-06 | 2017-09-06 | 超声成像方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710796923.6A CN109452954B (zh) | 2017-09-06 | 2017-09-06 | 超声成像方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109452954A CN109452954A (zh) | 2019-03-12 |
CN109452954B true CN109452954B (zh) | 2022-05-03 |
Family
ID=65606039
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710796923.6A Active CN109452954B (zh) | 2017-09-06 | 2017-09-06 | 超声成像方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109452954B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114025671A (zh) * | 2019-09-04 | 2022-02-08 | 深圳迈瑞生物医疗电子股份有限公司 | 一种vti测量装置及方法 |
CN113925528B (zh) * | 2020-06-29 | 2023-11-03 | 青岛海信医疗设备股份有限公司 | 多普勒成像方法和超声设备 |
WO2022099698A1 (zh) * | 2020-11-16 | 2022-05-19 | 深圳迈瑞生物医疗电子股份有限公司 | 超声成像系统、方法和计算机存储介质 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6390984B1 (en) * | 2000-09-14 | 2002-05-21 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus for locking sample volume onto moving vessel in pulsed doppler ultrasound imaging |
JP3610347B2 (ja) * | 2002-05-15 | 2005-01-12 | ジーイー横河メディカルシステム株式会社 | サンプリングボリューム位置の補正方法および超音波診断装置 |
EP2255730A4 (en) * | 2008-03-03 | 2014-12-10 | Konica Minolta Inc | ULTRASOUND |
CN102631224B (zh) * | 2012-04-24 | 2014-07-02 | 无锡祥生医学影像有限责任公司 | 一种用于触摸屏超声诊断仪的多普勒自动优化方法及装置 |
-
2017
- 2017-09-06 CN CN201710796923.6A patent/CN109452954B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109452954A (zh) | 2019-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8094893B2 (en) | Segmentation tool for identifying flow regions in an image system | |
KR101121396B1 (ko) | 2차원 초음파 영상에 대응하는 2차원 ct 영상을 제공하는 시스템 및 방법 | |
US7833159B2 (en) | Image processing system and method of enhancing the quality of an ultrasound image | |
WO2016067072A1 (en) | Imaging methods and apparatuses for performing shear wave elastography imaging | |
CN109452954B (zh) | 超声成像方法及装置 | |
KR101182891B1 (ko) | 2차원 초음파 영상과 3차원 초음파 영상의 합성영상을 제공하는 초음파 시스템 및 방법 | |
CN105246415A (zh) | 超声波观测装置、超声波观测装置的动作方法以及超声波观测装置的动作程序 | |
JP2001286474A (ja) | 対象物パラメータの動的測定 | |
JP2005193017A (ja) | 乳房患部分類の方法及びシステム | |
JP5467783B2 (ja) | 超音波システム及びクラッタ信号フィルタリング方法 | |
EP3820374B1 (en) | Methods and systems for performing fetal weight estimations | |
JP2011115580A (ja) | 変化推移映像を提供する超音波システムおよび方法 | |
JP2000300557A (ja) | 超音波診断装置 | |
US11324487B2 (en) | Ultrasound diagnostic apparatus and control method of ultrasound diagnostic apparatus | |
EP0934724A1 (en) | Ultrasonic diagnostic apparatus | |
JP2004267464A (ja) | 超音波診断装置 | |
CN114025672B (zh) | 一种超声成像设备和子宫内膜蠕动的检测方法 | |
KR101024857B1 (ko) | 3차원 초음파 영상에 컬러 모델링 처리를 수행하는 초음파 시스템 및 방법 | |
KR101117900B1 (ko) | 고유벡터를 설정하는 초음파 시스템 및 방법 | |
JP4634872B2 (ja) | 超音波診断装置 | |
US20060122505A1 (en) | M-Mode presentation of an ultrasound scan | |
CN112672695B (zh) | 用于执行胎儿重量估计的方法和系统 | |
JP2013183982A (ja) | 超音波診断装置及び弾性画像生成方法 | |
EP4285836A1 (en) | Ultrasound diagnostic apparatus and control method of ultrasound diagnostic apparatus | |
EP4223227A1 (en) | A method and system for performing fetal weight estimations |
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 |