CN108403148A - 一种基于mv自适应波束形成的超声ct成像方法 - Google Patents
一种基于mv自适应波束形成的超声ct成像方法 Download PDFInfo
- Publication number
- CN108403148A CN108403148A CN201810339093.9A CN201810339093A CN108403148A CN 108403148 A CN108403148 A CN 108403148A CN 201810339093 A CN201810339093 A CN 201810339093A CN 108403148 A CN108403148 A CN 108403148A
- Authority
- CN
- China
- Prior art keywords
- data
- signal
- ultrasonic
- covariance matrix
- array
- 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.)
- Pending
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/13—Tomography
- A61B8/14—Echo-tomography
-
- 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/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
-
- 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)
- Veterinary Medicine (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)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明涉及一种基于MV自适应波束形成的超声CT成像方法,包括数据采集、数据预处理、区域剖分、聚焦延时、数据处理和显示步骤;本发明根据接收的回波数据,计算不同通道数据的权重,将数据加权后再相加,得到最终成像点的值,该方法可有效的减少噪声的影响,提高图像分辨率和对比度。
Description
技术领域
本发明涉及超声成像技术领域,尤其涉及一种基于MV自适应波束形成的超声CT成像方法。
背景技术
超声CT是指在不损伤研究物体内部结构的基础上,通过超声设备测量物体在超声波的照射下投影数据,利用这些数据可以重建出二维或者三维的超声图像。超声检测具有指向性好、价格低廉、对人体无害、设备便于携带等优点,因此,用超声波作为发射源的检测技术取代射线来照射对象,已逐渐成为超声应用领域的研究者们追求的新目标之一。
超声成像的主要缺点是超声图像的分辨率低、对比度差。超声CT系统采用环阵探头,与线阵探头相比,它能够提供360°回波数据,对这些数据进行处理,可以进一步提高图像的分辨率,改善成像质量,有助于医生临床上对疾病的诊断和治疗。超声CT成像有两种方式,反射式成像和透射式成像,根据不同的成像方式选择对应的回波数据进行重建。
延时叠加波束形成是医学超声成像的标准技术,与传统的DAS相比,自适应波束形成的主要是计算每个成像点的值时增加了孔径的权重,该方法对于提高图像的空间分辨率有独特的优势。
目前国内外主要医学影像企业的的高端超声产品均没有将MV自适应波束的方法应用到超声CT的成像中。主要原因是该技术有几个难点尚未被攻克:环阵探头的阵元数大,通道数多,数据量大;MV自适应波束形成的重建算法虽然能够提高分辨率,但计算量巨大,重建需要花费大量的时间。
基于此,需要一种能提高超声CT图像的对比度和分辨率的基于MV自适应波束形成的超声CT成像方法被设计出来,利用超声CT成像系统采集的数据,应用MV自适应波束形成的算法,进行高分辨率图像重建。
发明内容
本发明的目的是为了克服现有技术的不足,得到一种能提高超声CT图像的对比度和分辨率的基于MV自适应波束形成的超声CT成像方法,利用超声CT成像系统采集的数据,应用MV自适应波束形成的算法,进行高分辨率图像重建。
本发明是通过以下技术方案实现:
一种基于MV自适应波束形成的超声CT成像方法,包括数据采集、数据预处理、区域剖分、聚焦延时、数据处理和显示步骤,
所述数据采集步骤包括:发射阵元按顺序依次发射信号,全部的阵元接收信号;
所述数据预处理步骤包括:保留发射阵元左右各45°接收阵元的回波数据,截掉发射阵元其余方向接收阵元的回波数据;
所述区域剖分步骤包括:根据成像区域的大小,选择合适的尺寸,将成像区域剖分为网格,每个格子中的区域作为成像点;
所述聚焦延时步骤包括:对经过数据预处理步骤得到的每个通道的回波数据进行聚焦延时;
所述数据处理步骤包括:
a1、计算各通道权重:采用MV自适应波束成形算法,对聚焦延时后的各通道回波数据计算权重;
a2、加权叠加:将计算的各通道权重进行叠加,波束形成器的输出为空间数据的加权和,即
其中,z(n)为加权和,wm(n)为阵元m的孔径权重,w(n)=[w0(n) w1(n) ... wM-1(n)]H为阵列的权重向量;
a3、加权平均:通过波束形成器寻找加权和的最小方差,采用空间平滑的方式,将阵列划分为多个重叠的子阵列,将得到的子阵列的协方差矩阵进行平均;
a4、包络检测:检出信号的上包络,提取出回波信号中携带的低频分量,得到被测物体信息。
a5、对数压缩:将包络检测后的信号取以10为底的对数,并乘以20的倍数,单位是dB,取完对数之后可以调整回波的动态范围,得到实时成像效果;
a6、灰度映射:采用线性映射,成比例地将最弱的信号映射到0dB,将最强的信号映射到255dB;
所述显示步骤包括:对经过灰度映射之后的数据进行图像显示。
作为优选,所述加权平均步骤中的具体算法为:
b1、通过波束形成器寻找z(n)的最小方差,同时保持聚焦点具有单位增益,它输出的方差为E[|z(n)|2]=w[n]HR[n]w[n],其中,E[|z(n)|2]表示期望算子,R[n]=E[X[n]X[n]H]为协方差矩阵;
b2、波束形成器根据接收的数据计算孔径的权重,寻找z(n)的最小方差来决定权重大小,当保证聚焦点具有单位增益时,最优化问题的数学表达式为subject to w[n]Ha=1,
其中,a为操纵向量,
b3、步骤b2中,协方差矩阵R[n]为样本协方差矩阵采用空间平滑的方式,将阵列划分为多个重叠的子阵列,将得到的子阵列的协方差矩阵进行平均,在某个采样点的协方差矩阵为
其中,L为子阵列长度,M为全部阵元个数,为第l个子阵列聚焦延时后的信号向量;
b4、利用样本协方差矩阵计算优化的孔径权重后,波束形成器最终的输出为
作为优选,所述加权平均步骤还包括采用对角加载的方法对协方差矩阵进行处理,将R[n]用R[n]+εI代替,ε为加载因子,可表示为
ε=Δ·tr{R[n]},
其中,tr{R[n]}为求矩阵的迹,Δ为常数,常设置为1/100L。
作为优选,所述包络检测方法为希尔伯特变换法,具体为,加权平均后的信号经过希尔伯特变换将得到加权平均后的信号的正交信号,以加权平均后的信号为实部,以希尔伯特变换得到的信号为虚部构造复信号,这个复信号的模就是所要求的实信号的包络。
作为优选,所述显示步骤通过调用MATLAB中的成像函数得到图像。
与现有的技术相比,本发明的有益效果是:本发明根据接收的回波数据,计算不同通道数据的权重,将数据加权后再相加,得到最终成像点的值,该方法可有效的减少噪声的影响,提高图像分辨率和对比度。
附图说明
图1为本发明采用的算法得到的重建图像;
图2为采用DAS算法得到的重建图像。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
一种基于MV自适应波束形成的超声CT成像方法,包括数据采集、数据预处理、区域剖分、聚焦延时、数据处理和显示步骤,
所述数据采集步骤包括:发射阵元按顺序依次发射信号,全部的阵元接收信号,本步骤中,采用单发全接的方式;
所述数据预处理步骤包括:保留发射阵元左右各45°接收阵元的回波数据,截掉发射阵元其余方向接收阵元的回波数据,本步骤中,由于数据量巨大,为了提高成像速度,所以要对数据先进行预处理;
所述区域剖分步骤包括:根据成像区域的大小,选择合适的尺寸,将成像区域剖分为网格,每个格子中的区域作为成像点,本步骤中,利用采集的数据成像,首先要对成像区域进行剖分,根据成像区域的大小,选择合适的尺寸剖分网格,再利用采集的原始回波数据,计算每个网格中发射阵元对他产生的影响,理论上剖分的网格数越多,分辨率越高,但由于衍射极限的影响,如果剖分的网格数过多,计算量加大,但成像效果并没有有效提高,因此要选择合适的网格大小,每个网格作为一个成像点,采用MV自适应波束成形算法计算成像点的值;
所述聚焦延时步骤包括:对经过数据预处理步骤得到的每个通道的回波数据进行聚焦延时;
所述数据处理步骤包括:
a1、计算各通道权重:采用MV自适应波束成形算法,对聚焦延时后的各通道回波数据计算权重,本步骤中,自适应波束形成的目标是在计算出最优的孔径权重对通道数据进行加权;
a2、加权叠加:将计算的各通道权重进行叠加,波束形成器的输出为空间数据的加权和,即
其中,z(n)为加权和,wm(n)为阵元m的孔径权重,w(n)=[w0(n) w1(n) ... wM-1(n)]H为阵列的权重向量;
a3、加权平均:通过波束形成器寻找加权和的最小方差,采用空间平滑的方式,将阵列划分为多个重叠的子阵列,将得到的子阵列的协方差矩阵进行平均;
a4、包络检测:检出信号的上包络,提取出回波信号中携带的低频分量,得到被测物体信息。
a5、对数压缩:将包络检测后的信号取以10为底的对数,并乘以20的倍数,单位是dB,取完对数之后可以调整回波的动态范围,得到实时成像效果,本步骤中,一般调整为40dB或60dB,其值越小对比度越高,具体方法,以60dB为例,是将信号中的最大值映射到60dB,将比最大信号小60dB的信号以及更小的信号映射到0dB;
a6、灰度映射:采用线性映射,成比例地将最弱的信号映射到0dB,将最强的信号映射到255dB;
所述显示步骤包括:对经过灰度映射之后的数据进行图像显示。
作为本发明一个较佳的实施例,所述加权平均步骤中的具体算法为:
b1、通过波束形成器寻找z(n)的最小方差,同时保持聚焦点具有单位增益,它输出的方差为E[|z(n)|2]=w[n]HR[n]w[n],其中,E[|z(n)|2]表示期望算子,R[n]=E[X[n]X[n]H]为协方差矩阵;
b2、波束形成器根据接收的数据计算孔径的权重,寻找z(n)的最小方差来决定权重大小,当保证聚焦点具有单位增益时,最优化问题的数学表达式为subject to w[n]Ha=1,
其中,a为操纵向量,本步骤中,由于数据进行了聚焦延时,这里a为简单的值全为1的向量,a包含了不同阵元的敏感程度,使所有的值不全为1,与预先确定的孔径变迹形成的固定波束样式相比,MV波束形成器可以自适应的记录数据,允许能量小的方向上有大的旁瓣,有干扰的方向上设置为0。对于复数数据,波束样式在操纵角度上是不对称的,由接收波束聚焦点有单位增益的约束条件,可以得到一个最优解;
b3、步骤b2中,协方差矩阵R[n]为样本协方差矩阵采用空间平滑的方式,将阵列划分为多个重叠的子阵列,将得到的子阵列的协方差矩阵进行平均,在某个采样点的协方差矩阵为
其中,L为子阵列长度,M为全部阵元个数,为第l个子阵列聚焦延时后的信号向量,在本步骤中,由于各路超声回波信号之间有较大的相关性,传统的MV波束形成不适合用于相干信号,为了获得一个好的评估值,可以采用空间平滑的方式,将阵列划分为多个重叠的子阵列,将得到的子阵列的协方差矩阵进行平均,应用样本协方差矩阵可以使权重系数w[n]的维数降低,限制了抑制相干信号和噪声的自由度,在精确估计协方差矩阵和选择空间滤波的长度上要有一个权衡,减小子阵列长度会导致分辨率的不稳定,分辨率会降低;
b4、利用样本协方差矩阵计算优化的孔径权重后,波束形成器最终的输出为
作为本发明一个较佳的实施例,所述加权平均步骤还包括采用对角加载的方法对协方差矩阵进行处理,将R[n]用R[n]+εI代替,ε为加载因子,可表示为
ε=Δ·tr{R[n]},
其中,tr{R[n]}为求矩阵的迹,Δ为常数,常设置为1/100L,本步骤中,MV波束形成器具有高分辨率,高信噪比的特点,但是它对扰动十分敏感,因此需要密集的角度采样增强它的鲁棒性,可以采用对角加载的方法对协方差矩阵进行处理。
根据上述的MV自适应波束形成的算法对各通道数据加权计算后,可以得到最终成像点的值,如果成像区域划分为N×N个格子,则用这种方式可以计算N×N个成像点的值。
作为本发明一个较佳的实施例,所述包络检测方法为希尔伯特变换法,具体为,加权平均后的信号经过希尔伯特变换将得到加权平均后的信号的正交信号,以加权平均后的信号为实部,以希尔伯特变换得到的信号为虚部构造复信号,这个复信号的模就是所要求的实信号的包络。
作为本发明一个较佳的实施例,所述显示步骤通过调用MATLAB中的成像函数得到图像。
本发明采用的方法的工作原理为:
UCT系统采用的环阵探头,采集一次即可获得大量原始数据,延时叠加波束形成的方法是进行反射重建的常用方法,该方法是将采集的数据通过延时后直接相加,结果为成像点的值。
而MV波束形成自适应算法可以根据接收的回波数据,计算不同通道数据的权重,将数据加权后再相加,得到最终成像点的值,该方法可有效的减少噪声的影响,提高图像分辨率。
采用超声反射接收的数据进行图像重建,假设超声在理想介质中传播,声速变化不大,反射型超声是根据回波信息来重建图象的。由于环阵探头阵元数较多,数据量较大,而实际的反射成像只需要用到发射阵元附近接收阵元的数据,因此,为了提高图像重建的速度,可先对采集到的数据进行初步的预处理。截掉远离发射阵元的数据,只保留发射阵元左右各45度接收阵元的回波数据。
利用预处理后的数据进行图像重建,采用的是MV(Minimum-variance最小方差)自适应波束形成的方法,该方法相当于对接收到的回波信号进行聚焦延时后再加权。首先将成像区域划分为N×N个网格,每个格子中的区域作为成像点,计算每个发射阵元对该成像点产生的影响,即对于每个发射阵元发出的信号来说,接收阵元在该点处的回波数据。根据MV自适应波束形成的算法计算每个通道上回波数据的权重,将所有的回波数据进行加权后再叠加,最后可得到该点数据值。在对每个格子中的区域进行计算完成后,可以得到N×N的矩阵数据,对数据进行对数压缩,灰度映射,得到重建的图像。
最小方差波束形成算法的核心是寻找一个最佳的加权矢量,在保持期望方向增益不变的条件下,阵列的输出能量最小。分辨率的提高是因为根据输出的最小方差,进而推导出最优权矢量,通道数据进行加权后再叠加,抵消了部分噪声的影响,提高了信噪比和对比度,如图1和图2所示,通过对比发现,获得了传统的合成孔径成像方法无法达到的高分辨率。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种基于MV自适应波束形成的超声CT成像方法,包括数据采集、数据预处理、区域剖分、聚焦延时、数据处理和显示步骤,其特征在于,
所述数据采集步骤包括:发射阵元按顺序依次发射信号,全部的阵元接收信号;
所述数据预处理步骤包括:保留发射阵元左右各45°接收阵元的回波数据,截掉发射阵元其余方向接收阵元的回波数据;
所述区域剖分步骤包括:根据成像区域的大小,选择合适的尺寸,将成像区域剖分为网格,每个格子中的区域作为成像点;
所述聚焦延时步骤包括:对经过数据预处理步骤得到的每个通道的回波数据进行聚焦延时;
所述数据处理步骤包括:
a1、计算各通道权重:采用MV自适应波束成形算法,对聚焦延时后的各通道回波数据计算权重;
a2、加权叠加:将计算的各通道权重进行叠加,波束形成器的输出为空间数据的加权和,即
其中,z(n)为加权和,wm(n)为阵元m的孔径权重,w(n)=[w0(n)w1(n)...wM-1(n)]H为阵列的权重向量;
a3、加权平均:通过波束形成器寻找加权和的最小方差,采用空间平滑的方式,将阵列划分为多个重叠的子阵列,将得到的子阵列的协方差矩阵进行平均;
a4、包络检测:检出信号的上包络,提取出回波信号中携带的低频分量,得到被测物体信息。
a5、对数压缩:将包络检测后的信号取以10为底的对数,并乘以20的倍数,单位是dB,取完对数之后可以调整回波的动态范围,得到实时成像效果;
a6、灰度映射:采用线性映射,成比例地将最弱的信号映射到0dB,将最强的信号映射到255dB;
所述显示步骤包括:对经过灰度映射之后的数据进行图像显示。
2.根据权利要求1所述的一种基于MV自适应波束形成的超声CT成像方法,其特征在于,所述加权平均步骤中的具体算法为:
b1、通过波束形成器寻找z(n)的最小方差,同时保持聚焦点具有单位增益,它输出的方差为E[|z(n)|2]=w[n]HR[n]w[n],其中,E[|z(n)|2]表示期望算子,R[n]=E[X[n]X[n]H]为协方差矩阵;
b2、波束形成器根据接收的数据计算孔径的权重,寻找z(n)的最小方差来决定权重大小,当保证聚焦点具有单位增益时,最优化问题的数学表达式为subject to w[n]Ha=1,
其中,a为操纵向量,
b3、步骤b2中,协方差矩阵R[n]为样本协方差矩阵采用空间平滑的方式,将阵列划分为多个重叠的子阵列,将得到的子阵列的协方差矩阵进行平均,在某个采样点的协方差矩阵为
其中,L为子阵列长度,M为全部阵元个数,为第l个子阵列聚焦延时后的信号向量;
b4、利用样本协方差矩阵计算优化的孔径权重后,波束形成器最终的输出为
3.根据权利要求2所述的一种基于MV自适应波束形成的超声CT成像方法,其特征在于,所述加权平均步骤还包括采用对角加载的方法对协方差矩阵进行处理,将R[n]用R[n]+εI代替,ε为加载因子,可表示为
ε=Δ·tr{R[n]},
其中,tr{R[n]}为求矩阵的迹,Δ为常数,常设置为1/100L。
4.根据权利要求1所述的一种基于MV自适应波束形成的超声CT成像方法,其特征在于:所述包络检测方法为希尔伯特变换法,具体为,加权平均后的信号经过希尔伯特变换将得到加权平均后的信号的正交信号,以加权平均后的信号为实部,以希尔伯特变换得到的信号为虚部构造复信号,这个复信号的模就是所要求的实信号的包络。
5.根据权利要求1所述的一种基于MV自适应波束形成的超声CT成像方法,其特征在于:所述显示步骤通过调用MATLAB中的成像函数得到图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810339093.9A CN108403148A (zh) | 2018-04-16 | 2018-04-16 | 一种基于mv自适应波束形成的超声ct成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810339093.9A CN108403148A (zh) | 2018-04-16 | 2018-04-16 | 一种基于mv自适应波束形成的超声ct成像方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108403148A true CN108403148A (zh) | 2018-08-17 |
Family
ID=63135526
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810339093.9A Pending CN108403148A (zh) | 2018-04-16 | 2018-04-16 | 一种基于mv自适应波束形成的超声ct成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108403148A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108618799A (zh) * | 2018-04-24 | 2018-10-09 | 华中科技大学 | 一种基于空间相干性的超声ct成像方法 |
CN110236593A (zh) * | 2019-07-11 | 2019-09-17 | 安阳工学院 | 模板匹配的医学超声波束形成方法 |
CN111466949A (zh) * | 2020-04-13 | 2020-07-31 | 剑桥大学南京科技创新中心有限公司 | 一种mmse波束形成器、mmse波束形成方法、计算机可读存储介质 |
CN113640809A (zh) * | 2021-08-23 | 2021-11-12 | 飞依诺科技(苏州)有限公司 | 超声成像降噪方法及超声成像设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102835975A (zh) * | 2012-09-19 | 2012-12-26 | 重庆博恩克医疗设备有限公司 | 最小方差波束形成与基于最小方差相干系数融合的方法 |
WO2014086695A1 (fr) * | 2012-12-07 | 2014-06-12 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Dispositif et procédé d'imagerie par ultrasons avec filtrage des artefacts dus aux interférences entre modes de reconstruction |
CN104970831A (zh) * | 2015-07-07 | 2015-10-14 | 重庆大学 | 一种基于特征结构的广义旁瓣相消超声成像波束合成方法 |
CN105411626A (zh) * | 2015-12-24 | 2016-03-23 | 华中科技大学 | 一种基于超声ct的合成孔径成像方法及系统 |
CN105760892A (zh) * | 2016-03-10 | 2016-07-13 | 重庆大学 | 一种改进的最小方差超声成像方法 |
-
2018
- 2018-04-16 CN CN201810339093.9A patent/CN108403148A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102835975A (zh) * | 2012-09-19 | 2012-12-26 | 重庆博恩克医疗设备有限公司 | 最小方差波束形成与基于最小方差相干系数融合的方法 |
WO2014086695A1 (fr) * | 2012-12-07 | 2014-06-12 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Dispositif et procédé d'imagerie par ultrasons avec filtrage des artefacts dus aux interférences entre modes de reconstruction |
CN104970831A (zh) * | 2015-07-07 | 2015-10-14 | 重庆大学 | 一种基于特征结构的广义旁瓣相消超声成像波束合成方法 |
CN105411626A (zh) * | 2015-12-24 | 2016-03-23 | 华中科技大学 | 一种基于超声ct的合成孔径成像方法及系统 |
CN105760892A (zh) * | 2016-03-10 | 2016-07-13 | 重庆大学 | 一种改进的最小方差超声成像方法 |
Non-Patent Citations (1)
Title |
---|
许梦玲: "医学超声成像中自适应波束形成方法研究", 《中国博士学位论文全文数据库 医药卫生科技辑》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108618799A (zh) * | 2018-04-24 | 2018-10-09 | 华中科技大学 | 一种基于空间相干性的超声ct成像方法 |
CN108618799B (zh) * | 2018-04-24 | 2020-06-02 | 华中科技大学 | 一种基于空间相干性的超声ct成像方法 |
CN110236593A (zh) * | 2019-07-11 | 2019-09-17 | 安阳工学院 | 模板匹配的医学超声波束形成方法 |
CN110236593B (zh) * | 2019-07-11 | 2021-12-28 | 安阳工学院 | 一种模板匹配的医学超声波束形成方法 |
CN111466949A (zh) * | 2020-04-13 | 2020-07-31 | 剑桥大学南京科技创新中心有限公司 | 一种mmse波束形成器、mmse波束形成方法、计算机可读存储介质 |
CN113640809A (zh) * | 2021-08-23 | 2021-11-12 | 飞依诺科技(苏州)有限公司 | 超声成像降噪方法及超声成像设备 |
CN113640809B (zh) * | 2021-08-23 | 2023-08-25 | 飞依诺科技股份有限公司 | 超声成像降噪方法及超声成像设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108403148A (zh) | 一种基于mv自适应波束形成的超声ct成像方法 | |
Montaldo et al. | Coherent plane-wave compounding for very high frame rate ultrasonography and transient elastography | |
US9451932B2 (en) | Clutter suppression in ultrasonic imaging systems | |
CN105411626B (zh) | 一种基于超声ct的合成孔径成像方法及系统 | |
US11030780B2 (en) | Ultrasound speckle reduction and image reconstruction using deep learning techniques | |
US6475150B2 (en) | System and method for ultrasonic tomography | |
CN102640014B (zh) | 图像生成装置和图像生成方法 | |
CN108836389B (zh) | 平面波相关点相干自适应波束合成成像方法 | |
US20130343627A1 (en) | Suppression of reverberations and/or clutter in ultrasonic imaging systems | |
EP2820445A2 (en) | Clutter suppression in ultrasonic imaging systems | |
CN103536316A (zh) | 一种空时平滑相干因子类自适应超声成像方法 | |
Agarwal et al. | Improving spatial resolution using incoherent subtraction of receive beams having different apodizations | |
Viola et al. | Time-domain optimized near-field estimator for ultrasound imaging: Initial development and results | |
CN103969651A (zh) | 自适应声学成像方法 | |
JP2022016414A (ja) | 超音波減衰パラメータを推定するための方法及びシステム | |
CN102764139A (zh) | 基于特征空间分析和区域判别的医学超声波束形成方法 | |
CN108309352A (zh) | 一种余弦变换域超声成像方法 | |
Hooi et al. | Analysis of tissue changes, measurement system effects, and motion artifacts in echo decorrelation imaging | |
Thon et al. | Detection of point scatterers in medical ultrasound | |
Andresen et al. | Synthetic aperture focusing for a single-element transducer undergoing helical motion | |
CN108700651A (zh) | 成像方法、实施所述方法的装置、计算机程序以及计算机可读存储介质 | |
Walker | C-and D-weighted ultrasonic imaging using the translating apertures algorithm | |
Peretz et al. | Deep learning applied to beamforming in synthetic aperture ultrasound | |
Yamaryo et al. | Range-point migration-based image expansion method exploiting fully polarimetric data for UWB short-range radar | |
Guenther et al. | Generalized cystic resolution: a metric for assessing the fundamental limits on beamformer performance |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180817 |