CN108836392A - 基于超声rf信号的超声成像方法、装置、设备及存储介质 - Google Patents
基于超声rf信号的超声成像方法、装置、设备及存储介质 Download PDFInfo
- Publication number
- CN108836392A CN108836392A CN201810286743.8A CN201810286743A CN108836392A CN 108836392 A CN108836392 A CN 108836392A CN 201810286743 A CN201810286743 A CN 201810286743A CN 108836392 A CN108836392 A CN 108836392A
- Authority
- CN
- China
- Prior art keywords
- image
- layered
- interpolation
- resolution
- signal
- 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
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/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/06—Measuring blood flow
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Hematology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明适用计算机技术领域,提供了一种基于超声RF信号的超声成像方法、装置、设备及存储介质,该方法包括:接收采集的超声回波射频信号,沿超声回波射频信号对应的原始超声图像的深度方向,对原始超声图像进行分层,以得到多个分层图像,对每个分层图像进行去噪处理,以从每个分层图像中提取对应的分层血流图像,对每个分层血流图像进行插值,以得到插值后的高分辨插值图像,将高分辨插值图像叠加重组成一幅图像,以得到超声回波射频信号对应的高分辨率超声图像,最终得到超声回波射频信号对应的高分辨率超声图像,从而提高了血管图像的重建效果。
Description
技术领域
本发明属于医疗设备技术领域,尤其涉及一种基于超声RF信号的超声成像方法、装置、设备及存储介质。
背景技术
人体各脏器的血管分布情况是临床上诊断一些病灶的重要指标。现在主要血流成像手段有CT血管造影技术(CTA)、磁共振血管成像技术(MRA)、数字减影血管造影技术(DSA)、超声多普勒成像技术(CDFI)。由于CT和X光对人体存在一定的辐射危害,而超声具有无创性、安全、成像速度快、使用方便等优点,因此,超声成像在临床诊断中大量被采用,是临床诊断的重要工具之一。当前,超声成像系统是基于发射聚焦超声波提高图像质量。一幅完整的扫描图像是由若干次发射若干条扫描线组成,如果采用多个聚焦位置,则需要在不同聚焦位置重新扫描,该过程极其耗时,大幅降低了超声图像的显示帧频。彩色超声多普勒血流成像技术多采用2.0MHz-6.0MHz频率范围内的超声,可以分辨直径1mm以上且流速为1-2cm/s以上的血流。这种成像方法对血流的分辨率远不及CTA、MRA、DSA等技术。因为声波在人体中穿透能力不及射线,且当成像物体尺寸小于波长时会发生衍射现象,导致波长越长,理其理论分辨率越高,这种现象成为衍射极限。
超声造影剂(UCA)也称超声微泡(MB)最初被用于诊断显像诊断,其原理是利用声波对气体反射比液体大近1000倍的原理,通过声衰减、声速的改变和增强后向散射等方式改变声波与组织的吸收、反射和折射等相互作用,使超声回声增强,从而得到更高的对比分辨力,达到有利于诊断疾病的目的。目前临床常用的第三代造影剂微泡直径较小且趋于一致,理化性能更稳定,稳定性和抗压性更高,例如SonoVue等。造影剂的利用提高了超声血流成像的质量,提高了分辨率,但是图像质量仍然受到衍射极限的限制。美国及德国三位科学家Eric Betzig、Stefan W.Hell和William E.Moerner凭借“研制出超分辨率荧光显微镜”突破了光学成像衍射极限,克服了光学显微镜中的“阿贝分变率”,因此获得2014年诺贝尔化学奖。赖利、克里斯滕森等人受到启发,提出了基于点扩散函数调制的超分辨技术,将多角度平面波成像技术与超声造影剂微泡相结合,通过反卷积点扩散函数求解单帧图像中微泡的位置。
超声平面波成像方法是利用平面波扫描检测区域,当所有超声换能器阵元同时发射相同的声波信号,由于各阵元间的间隙极小,因此,在传播过程中,声场叠加,波阵面可近似为平面,这种扫描方式一次发射控制即可完成扫描,扫描区域是完整的全幅感兴趣区,一幅图像花费的时间远低于传统聚焦逐线扫描。虽然一次发射平面波可以极大提高超声的帧频,但是由于平面波是非聚焦波,得到的图像对比度和分辨率较低。多角度平面波相干复合成像方法,将多个不同角度的平面波辐射感兴趣区域,相干叠加后可以提高图像质量。但是利用多角度平面波对较深的感兴趣区域进行成像时,近场的目标血流信号会受到远场信号的影响,导致近场信号被抑制,对于复杂的血管区域,例如脑血流,当微泡信号出现在近场区域时,由于被抑制,在微泡定位算法的过程中会被忽略而无法被提取出来,也使得每帧图像中能提取到的微泡数量降低,导致近场的血管重建效果不好。
发明内容
本发明的目的在于提供一种基于超声RF信号的超声成像方法、装置、设备及存储介质,旨在解决由于现有超声成像方法的血管重建效果不佳的问题。
一方面,本发明提供了一种基于超声RF信号的超声成像方法,所述方法包括下述步骤:
接收采集的超声回波射频信号,沿所述超声回波射频信号对应的原始超声图像的深度方向,对所述原始超声图像进行分层,以得到多个分层图像;
对每个所述分层图像进行去噪处理,以从每个所述分层图像中提取对应的分层血流图像;
对所述每个所述分层血流图像进行插值,以得到插值后的高分辨插值图像;
将所述高分辨插值图像叠加重组成一幅图像,以得到所述超声回波射频信号对应的高分辨率超声图像。
另一方面,本发明提供了一种基于超声RF信号的超声成像装置,所述装置包括:
图像分层单元,用于接收采集的超声回波射频信号,沿所述超声回波射频信号对应的原始超声图像的深度方向,对所述原始超声图像进行分层,以得到多个分层图像;
图像去噪单元,用于对每个所述分层图像进行去噪处理,以从每个所述分层图像中提取对应的分层血流图像;
图像插值单元,用于对所述每个所述分层血流图像进行插值,以得到插值后的高分辨插值图像;以及
图像重建单元,用于将所述高分辨插值图像叠加重组成一幅图像,以得到所述超声回波射频信号对应的高分辨率超声图像。。
另一方面,本发明还提供了一种医疗设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现如上所述方法的步骤。
另一方面,本发明还提供了一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现如上所述方法的步骤。
本发明在接收采集的超声回波射频信号,沿超声回波射频信号对应的原始超声图像的深度方向,对原始超声图像进行分层,以得到多个分层图像,对每个分层图像进行去噪处理,以从每个分层图像中提取对应的分层血流图像,对每个分层血流图像进行插值,以得到插值后的高分辨插值图像,将高分辨插值图像叠加重组成一幅图像,以得到超声回波射频信号对应的高分辨率超声图像,最终得到超声回波射频信号对应的高分辨率超声图像,从而提高了血管图像的重建效果。
附图说明
图1是本发明实施例一提供的基于超声RF信号的超声成像方法的实现流程图;
图2是本发明实施例一提供的基于超声RF信号的超声成像方法对原始超声图像进行分层的示意图;
图3是本发明实施例二提供的基于超声RF信号的超声成像方法的实现流程图;
图4是本发明实施例三提供的基于超声RF信号的超声成像装置的结构示意图;
图5是本发明实施例四提供的基于超声RF信号的超声成像装置的结构示意图;以及
图6是本发明实施例四提供的医疗设备的结构示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
以下结合具体实施例对本发明的具体实现进行详细描述:
实施例一:
图1示出了本发明实施例一提供的基于超声RF信号的超声成像方法的实现流程,为了便于说明,仅示出了与本发明实施例相关的部分,详述如下:
在步骤S101中,接收采集的超声回波射频信号,沿超声回波射频信号对应的原始超声图像的深度方向,对原始超声图像进行分层,以得到多个分层图像。
本发明实施例适用于医疗设备,例如,超声成像设备或系统,以根据超声成像设备或系统接收到的超声回波射频(RF)信号进行成像。在本发明实施例中,在对仿体或活体动物注射完微泡造影剂之后,通过调整超声成像设备或系统探头的位置获得目标成像区域。优选地,微泡造影剂为声诺维的六氟化硫微泡,以得到较好的回波射频信息,便于后续成像。通过对超声回波射频信号进行模-数转换,可得到超声回波射频信号数据,即三维射频数据,该三维射频数据可表示为(阵元方向,成像深度,时间)三元组。
在本发明实施例中,如图2所示,在得到三维射频数据表示的原始超声图像S(Nz,Nx,Nt)后,沿该原始超声图像的深度方向(Z轴)对原始超声图像进行分层,以得到多个分层图像,从而提高后续去噪处理和微泡定位处理的效果。
在步骤S102中,对每个分层图像进行去噪处理,以从每个分层图像中提取对应的分层血流图像。
在本发明实施例中,在对每个分层图像进行去噪处理时,优选地,通过差分处理去除每个分层图像中的部分静态组织和血管壁等信号,接着通过预设的壁滤波器对差分处理后的每个分层图像进行滤波,以去除差分处理后的每个分层图像中的杂波信号,以从每个分层图像中提取对应的分层血流信号,获得分层血流图像,从而提高分层血流图像的分辨率和信噪比。
在具体实施中,壁滤波器可以是无限冲激响应高通滤波器、自适应高通滤波器。优选地,壁滤波器为SVD(Singular Value Decomposition,奇异值分解)时频滤波器,通过对差分处理后的每个分层图像进行分解,获得差分处理后的每个分层图像的奇异值之后,对奇异值从大到小进行排列,获取预设数量个中间部分的奇异值,以从每个分层图像中提取对应的分层血流信号,获得分层血流图像,从而在去除杂波信号的同时,提高提取的分层血流图像的分辨率和信噪比,获得更佳的滤波效果。
在步骤S103中,对每个分层血流图像进行插值,以得到插值后的高分辨插值图像。
在本发明实施例中,对每个分层血流图像的每帧图像数据在X/Z平面进行插值,每层血流图像的插值方式和插值数相同,方便后面将各分层重建成全幅图像。具体地,插值方式可以为线性插值、临近插值、高斯插值等。
在步骤S104中,将高分辨插值图像叠加重组成一幅图像,以得到超声回波射频信号对应的高分辨率超声图像。
在本发明实施例中,沿超声回波射频信号对应的原始超声图像的深度方向,对原始超声图像进行分层,以得到多个分层图像,对每个分层图像进行去噪处理,以从每个分层图像中提取对应的分层血流图像,对每个分层血流图像进行插值,以得到插值后的高分辨插值图像,根据每个高分辨插值图像中的每个像素对应的射频信号强度值、以及预设的强度阈值,从每个高分辨插值图像中提取对应的微泡定位图像,将提取得到的多层微泡定位图像叠加重组成一幅图像,最终得到超声回波射频信号对应的高分辨率超声图像,从而通过对原始超声图像进行分层处理,提高了血管图像的重建效果。
实施例二:
图3示出了本发明实施例二提供的基于超声RF信号的超声成像方法的实现流程,为了便于说明,仅示出了与本发明实施例相关的部分,详述如下:
在步骤S301中,接收采集的超声回波射频信号,沿超声回波射频信号对应的原始超声图像的深度方向,对原始超声图像进行分层,以得到多个分层图像。
本发明实施例适用于医疗设备,例如,超声成像设备或系统,以根据超声成像设备或系统接收到的超声回波射频信号进行成像。在本发明实施例中,在对仿体或活体动物注射完微泡造影剂之后,通过调整超声成像设备或系统探头的位置获得目标成像区域。优选地,微泡造影剂为声诺维的六氟化硫微泡,以得到较好的回波射频信息,便于后续成像。通过对超声回波射频信号进行模-数转换,可得到超声回波射频信号数据,即三维射频数据,该三维射频数据可表示为(阵元方向,成像深度,时间)三元组。
在步骤S302中,对每个分层图像进行去噪处理,以从每个分层图像中提取对应的分层血流图像。
在本发明实施例中,在对每个分层图像进行去噪处理时,优选地,通过差分处理去除每个分层图像中的部分静态组织和血管壁等信号,接着通过预设的壁滤波器对差分处理后的每个分层图像进行滤波,以去除差分处理后的每个分层图像中的杂波信号,以从每个分层图像中提取对应的分层血流信号,获得分层血流图像,从而提高分层血流图像的分辨率和信噪比。
在具体实施中,壁滤波器可以是无限冲激响应高通滤波器、自适应高通滤波器。优选地,壁滤波器为SVD时频滤波器,通过对差分处理后的每个分层图像进行分解,获得差分处理后的每个分层图像的奇异值之后,对奇异值从大到小进行排列,获取预设数量个中间部分的奇异值,以从每个分层图像中提取对应的分层血流信号,获得分层血流图像,从而在去除杂波信号的同时,提高提取的分层血流图像的分辨率和信噪比,获得更佳的滤波效果。
在步骤S303中,对每个分层血流图像进行插值,以得到插值后的高分辨插值图像。
在本发明实施例中,对每个分层血流图像的每帧图像数据在X/Z平面进行插值,每层血流图像的插值方式和插值数相同,方便后面将各分层重建成全幅图像。具体地,插值方式可以为线性插值、临近插值、高斯插值等。
在步骤S304中,根据每个高分辨插值图像中的每个像素对应的射频信号强度值、以及预设的强度阈值,从每个高分辨插值图像中提取对应的微泡定位图像。
在本发明实施例中,由于微泡对超声的后向反射回波要远大于血流的反射回波,因此,可预先设置合理的射频信号阈值,以将微泡信号从每个分层血流图像提取出来,以去除掉多余的血流信号。优选地,强度阈值设置为每个高分辨插值图像中最大强度值的0.7-0.9倍,从而在不同深度的区域内采用不同的强度阈值进行微泡提取,使得微泡提取更准确,并在准确度和提取微泡数量之间达到较好的平衡。
进一步优选地,从每个高分辨插值图像中提取对应的微泡定位图像时,将每个像素对应的射频信号强度值与强度阈值进行相减,将小于强度阈值的射频信号强度值对应的像素从每个分层血流图像中删除,以得到每个分层血流图像对应的微泡信号分布图,对每个微泡信号分布图中的微泡进行定位,用单个像素对定位得到的每个微泡进行表示,以得到每个分层血流图像对应的微泡定位图像。具体地,在得到每个分层血流图像对应的微泡信号分布图后,可通过强度加权、最大值点、50%最大值点等算法,对每个微泡信号分布图中的微泡进行定位,以得到每个微泡信号的中心点,用单个像素表示该中心,从而得到微泡定位图像。
在步骤S305中,将提取得到的多层微泡定位图像叠加重组成一幅图像,以得到所述超声回波射频信号对应的超声分辨率图像。
在本发明实施例中,沿超声回波射频信号对应的原始超声图像的深度方向,对原始超声图像进行分层,以得到多个分层图像,对每个分层图像进行去噪处理,以从每个分层图像中提取对应的分层血流图像,对每个分层血流图像进行插值,以得到插值后的高分辨插值图像,根据每个高分辨插值图像中的每个像素对应的射频信号强度值、以及预设的强度阈值,从每个高分辨插值图像中提取对应的微泡定位图像,将提取得到的多层微泡定位图像叠加重组成一幅图像,最终得到超声回波射频信号对应的超声分辨率图像,从而通过对原始超声图像进行分层处理,提高了血管图像的重建效果。
实施例三:
图4示出了本发明实施例三提供的基于超声RF信号的超声成像装置的结构,为了便于说明,仅示出了与本发明实施例相关的部分,其中包括:
图像分层单元41,用于接收采集的超声回波射频信号,沿超声回波射频信号对应的原始超声图像的深度方向,对原始超声图像进行分层,以得到多个分层图像;
图像去噪单元42,用于对每个分层图像进行去噪处理,以从每个分层图像中提取对应的分层血流图像;
图像插值单元43,用于对每个分层血流图像进行插值,以得到插值后的高分辨插值图像;以及
图像重建单元44,用于将高分辨插值图像叠加重组成一幅图像,以得到超声回波射频信号对应的高分辨率超声图像。
优选地,图像去噪单元42包括第一去噪子单元和第二去噪子单元,第一去噪子单元用于通过差分处理去除每个分层图像中的部分静态组织和血管壁信号,第二去噪子单元用于通过预设的壁滤波器对差分处理后的每个分层图像进行滤波,以去除差分处理后的每个分层图像中的杂波信号,以从每个分层图像中提取对应的分层血流信号,获得分层血流图像。
在本发明实施例中,超声成像装置的各单元可由相应的硬件或软件单元实现,各单元可以为独立的软、硬件单元,也可以集成为一个软、硬件单元,在此不用以限制本发明。各单元的具体实施方式可参考实施一的描述,在此不再赘述。
实施例四:
图5示出了本发明实施例四提供的基于超声RF信号的超声成像装置的结构,为了便于说明,仅示出了与本发明实施例相关的部分,其中包括:
图像分层单元51,用于接收采集的超声回波射频信号,沿超声回波射频信号对应的原始超声图像的深度方向,对原始超声图像进行分层,以得到多个分层图像;
图像去噪单元52,用于对每个分层图像进行去噪处理,以从每个分层图像中提取对应的分层血流图像;
图像插值单元53,用于对每个分层血流图像进行插值,以得到插值后的高分辨插值图像;
图像提取单元54,用于根据每个高分辨插值图像中的每个像素对应的射频信号强度值、以及预设的强度阈值,从每个高分辨插值图像中提取对应的微泡定位图像;以及
图像重建单元55,用于将高分辨插值图像叠加重组成一幅图像,以得到超声回波射频信号对应的高分辨率超声图像,
在本发明实施例中,图像重建单元55包括超声分辨率图像重建单元551,以用于将提取到的多层微泡定位图像叠加重组成一幅图像,以得到超声回波射频信号对应的超声分辨率图像。
优选地,图像提取单元54包括微泡获取单元和微泡定位单元。其中,微泡获取单元用于将每个像素对应的射频信号强度值与强度阈值进行相减,将小于强度阈值的射频信号强度值对应的像素从每个分层血流图像中删除,以得到每个分层血流图像对应的微泡信号分布图。微泡定位单元用于对每个微泡信号分布图中的微泡进行定位,用单个像素对定位得到的每个微泡进行表示,以得到每个分层血流图像对应的微泡定位图像。
在本发明实施例中,超声成像装置的各单元可由相应的硬件或软件单元实现,各单元可以为独立的软、硬件单元,也可以集成为一个软、硬件单元,在此不用以限制本发明。各单元的具体实施方式可参考实施二的描述,在此不再赘述。
实施例五:
图6示出了本发明实施例五提供的医疗设备的结构,为了便于说明,仅示出了与本发明实施例相关的部分。
本发明实施例的医疗设备6包括处理器60、存储器61以及存储在存储器61中并可在处理器60上运行的计算机程序62。该处理器60执行计算机程序62时实现上述方法实施例中的步骤,例如图1所示的步骤S101至S104。或者,处理器60执行计算机程序62时实现上述各装置实施例中各单元的功能,例如图4所示单元41至44的功能。
该处理器60执行计算机程序62时,接收采集的超声回波射频信号,沿超声回波射频信号对应的原始超声图像的深度方向,对原始超声图像进行分层,以得到多个分层图像,对每个分层图像进行去噪处理,以从每个分层图像中提取对应的分层血流图像,对每个分层血流图像进行插值,以得到插值后的高分辨插值图像,将高分辨插值图像叠加重组成一幅图像,以得到超声回波射频信号对应的高分辨率超声图像,最终得到超声回波射频信号对应的高分辨率超声图像,从而提高了血管图像的重建效果。
本发明实施例的医疗设备可以为彩超等超声诊断设备。该医疗设备6中处理器60执行计算机程序62时实现基于超声RF信号的超声成像方法时实现的步骤可参考前述方法实施例的描述,在此不再赘述。
实施例六:
在本发明实施例中,提供了一种计算机可读存储介质,该计算机可读存储介质存储有计算机程序,该计算机程序被处理器执行时实现上述方法实施例中的步骤,例如,图1所示的步骤S101至S104。或者,该计算机程序被处理器执行时实现上述各装置实施例中各单元的功能,例如图4所示单元41至44的功能。
该计算机程序被处理器执行时,接收采集的超声回波射频信号,沿超声回波射频信号对应的原始超声图像的深度方向,对原始超声图像进行分层,以得到多个分层图像,对每个分层图像进行去噪处理,以从每个分层图像中提取对应的分层血流图像,对每个分层血流图像进行插值,以得到插值后的高分辨插值图像,将高分辨插值图像叠加重组成一幅图像,以得到超声回波射频信号对应的高分辨率超声图像,最终得到超声回波射频信号对应的高分辨率超声图像,从而提高了血管图像的重建效果。
本发明实施例的计算机可读存储介质可以包括能够携带计算机程序代码的任何实体或装置、记录介质,例如,ROM/RAM、磁盘、光盘、闪存等存储器。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于超声RF信号的超声成像方法,其特征在于,所述方法包括:
接收采集的超声回波射频信号,沿所述超声回波射频信号对应的原始超声图像的深度方向,对所述原始超声图像进行分层,以得到多个分层图像;
对每个所述分层图像进行去噪处理,以从每个所述分层图像中提取对应的分层血流图像;
对所述每个所述分层血流图像进行插值,以得到插值后的高分辨插值图像;
将所述高分辨插值图像叠加重组成一幅图像,以得到所述超声回波射频信号对应的高分辨率超声图像。
2.如权利要求1所述的方法,其特征在于,对每个所述分层图像进行去噪处理,以从每个所述分层图像中提取对应的分层血流图像的步骤,包括:
通过差分处理去除每个所述分层图像中的部分静态组织和血管壁信号;
通过预设的壁滤波器对所述差分处理后的每个所述分层图像进行滤波,以去除所述差分处理后的每个所述分层图像中的杂波信号,以从每个所述分层图像中提取对应的分层血流信号,获得分层血流图像。
3.如权利要求1所述的方法,其特征在于,
对所述每个所述分层血流图像进行插值的步骤之后,将所述高分辨插值图像叠加重组成一幅图像的步骤之前,所述方法还包括:
根据每个所述高分辨插值图像中的每个像素对应的射频信号强度值、以及预设的强度阈值,从每个所述高分辨插值图像中提取对应的微泡定位图像;
将所述高分辨插值图像叠加重组成一幅图像,以得到所述超声回波射频信号对应的高分辨率超声图像的步骤,包括:
将提取到的多层所述微泡定位图像叠加重组成一幅图像,以得到所述超声回波射频信号对应的超声分辨率图像。
4.如权利要求3所述的方法,其特征在于,从每个所述高分辨插值图像中提取对应的微泡定位图像的步骤,包括:
将所述每个像素对应的射频信号强度值与所述强度阈值进行相减,将小于所述强度阈值的所述射频信号强度值对应的像素从所述每个分层血流图像中删除,以得到每个分层血流图像对应的微泡信号分布图;
对每个所述微泡信号分布图中的微泡进行定位,用单个像素对所述定位得到的每个微泡进行表示,以得到每个分层血流图像对应的微泡定位图像。
5.一种基于超声RF信号的超声成像装置,其特征在于,所述装置包括:
图像分层单元,用于接收采集的超声回波射频信号,沿所述超声回波射频信号对应的原始超声图像的深度方向,对所述原始超声图像进行分层,以得到多个分层图像;
图像去噪单元,用于对每个所述分层图像进行去噪处理,以从每个所述分层图像中提取对应的分层血流图像;
图像插值单元,用于对所述每个所述分层血流图像进行插值,以得到插值后的高分辨插值图像;以及
图像重建单元,用于将所述高分辨插值图像叠加重组成一幅图像,以得到所述超声回波射频信号对应的高分辨率超声图像。
6.如权利要求5所述的超声成像装置,其特征在于,所述图像去噪单元包括:
第一去噪子单元,用于通过差分处理去除每个所述分层图像中的部分静态组织和血管壁信号;以及
第二去噪子单元,用于通过预设的壁滤波器对所述差分处理后的每个所述分层图像进行滤波,以去除所述差分处理后的每个所述分层图像中的杂波信号,以从每个所述分层图像中提取对应的分层血流信号,获得分层血流图像。
7.如权利要求5所述的超声成像装置,其特征在于,
所述装置还包括:
图像提取单元,用于根据每个所述高分辨插值图像中的每个像素对应的射频信号强度值、以及预设的强度阈值,从每个所述高分辨插值图像中提取对应的微泡定位图像;
所述图像重建单元包括:
超声分辨率图像重建单元,用于将提取到的多层所述微泡定位图像叠加重组成一幅图像,以得到所述超声回波射频信号对应的超声分辨率图像。
8.如权利要求7所述的超声成像装置,其特征在于,所述图像提取单元包括:
微泡获取单元,用于将所述每个像素对应的射频信号强度值与所述强度阈值进行相减,将小于所述强度阈值的所述射频信号强度值对应的像素从所述每个分层血流图像中删除,以得到每个分层血流图像对应的微泡信号分布图;以及
微泡定位单元,用于对每个所述微泡信号分布图中的微泡进行定位,用单个像素对所述定位得到的每个微泡进行表示,以得到每个分层血流图像对应的微泡定位图像。
9.一种医疗设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现如权利要求1至4任一项所述方法的步骤。
10.一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1至4任一项所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810286743.8A CN108836392B (zh) | 2018-03-30 | 2018-03-30 | 基于超声rf信号的超声成像方法、装置、设备及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810286743.8A CN108836392B (zh) | 2018-03-30 | 2018-03-30 | 基于超声rf信号的超声成像方法、装置、设备及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108836392A true CN108836392A (zh) | 2018-11-20 |
CN108836392B CN108836392B (zh) | 2021-06-22 |
Family
ID=64211759
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810286743.8A Active CN108836392B (zh) | 2018-03-30 | 2018-03-30 | 基于超声rf信号的超声成像方法、装置、设备及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108836392B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110327077A (zh) * | 2019-07-09 | 2019-10-15 | 深圳开立生物医疗科技股份有限公司 | 一种血流显示方法、装置及超声设备和存储介质 |
CN110772285A (zh) * | 2019-10-31 | 2020-02-11 | 南京景瑞康分子医药科技有限公司 | 一种超声超分辨成像方法 |
CN112043308A (zh) * | 2020-08-31 | 2020-12-08 | 深圳市德力凯医疗设备股份有限公司 | 一种颅内三维脑血流重建方法、存储介质及超声设备 |
CN112370073A (zh) * | 2020-10-12 | 2021-02-19 | 深圳华声医疗技术股份有限公司 | 超声声速校正方法、超声信号处理设备及存储介质 |
CN112869767A (zh) * | 2021-01-11 | 2021-06-01 | 青岛海信医疗设备股份有限公司 | 一种超声图像存储方法、装置及其超声设备 |
CN112932539A (zh) * | 2021-01-23 | 2021-06-11 | 西安交通大学 | 一种经颅超声估计血管血流速度场的装置及方法 |
CN113180734A (zh) * | 2018-12-27 | 2021-07-30 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声血流成像方法及系统 |
WO2021233188A1 (zh) * | 2020-05-18 | 2021-11-25 | 南京超维景生物科技有限公司 | 三维超声造影图像的超分辨重建方法及装置 |
CN114073547A (zh) * | 2020-08-14 | 2022-02-22 | 深圳迈瑞生物医疗电子股份有限公司 | 用于超声成像的方法、装置和存储介质 |
CN114129190A (zh) * | 2021-11-30 | 2022-03-04 | 深圳先进技术研究院 | 全眼微血流的成像方法、装置及电子设备和存储介质 |
Citations (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1260730A (zh) * | 1997-06-19 | 2000-07-19 | 埃莱克特公司 | 自动化剂量规划法 |
JP2000287972A (ja) * | 1999-04-05 | 2000-10-17 | Toshiba Corp | 超音波診断装置および超音波診断装置の画像生成方法 |
US20030163040A1 (en) * | 2002-02-28 | 2003-08-28 | Philip Gildenberg | Audible feedback from positional guidance systems |
US6645147B1 (en) * | 1998-11-25 | 2003-11-11 | Acuson Corporation | Diagnostic medical ultrasound image and system for contrast agent imaging |
US20050090742A1 (en) * | 2003-08-19 | 2005-04-28 | Yoshitaka Mine | Ultrasonic diagnostic apparatus |
CN1792334A (zh) * | 2005-10-31 | 2006-06-28 | 西安交通大学 | 基于包膜微泡的灌注成像与超声控制释放的系统和方法 |
US20060241462A1 (en) * | 2005-03-14 | 2006-10-26 | Yi-Hong Chou | Method of intracranial ultrasound imaging and related system |
CN1903121A (zh) * | 2006-08-25 | 2007-01-31 | 北京源德生物医学工程有限公司 | 永磁磁共振图像导引体外高能聚焦超声系统和方法 |
CN1947148A (zh) * | 2004-04-30 | 2007-04-11 | 埃莱克塔公共有限公司 | 用于自动提高医学图片的可用性的方法和系统 |
CN101011264A (zh) * | 2006-01-25 | 2007-08-08 | 株式会社东芝 | 超声波探头和超声波图像取得装置 |
CN101061963A (zh) * | 2006-04-27 | 2007-10-31 | 通用电气公司 | 用于测量通过心脏瓣膜血流的方法和系统 |
CN101184428A (zh) * | 2005-04-01 | 2008-05-21 | 视声公司 | 利用超声波使血管结构三维可视化的系统和方法 |
CN101291635A (zh) * | 2005-10-20 | 2008-10-22 | 直观外科手术公司 | 医用机器人系统中的计算机显示屏上的辅助图像显示和操纵 |
CN101574267A (zh) * | 2008-05-08 | 2009-11-11 | 重庆医科大学超声影像学研究所 | 对超声微泡计数的方法 |
CN101577001A (zh) * | 2009-05-20 | 2009-11-11 | 电子科技大学 | 基于光线投射体绘制的三维医学图像剖分方法 |
CN101912273A (zh) * | 2010-07-23 | 2010-12-15 | 深圳市理邦精密仪器股份有限公司 | 超声动态图像处理方法及系统 |
CN202015163U (zh) * | 2010-12-30 | 2011-10-26 | 广州宝胆医疗器械科技有限公司 | 一种表层血管显示仪 |
CN102317873A (zh) * | 2008-12-10 | 2012-01-11 | 赫勒拉德有限责任公司 | 用于彩色运动全息术的系统和方法 |
CN102697523A (zh) * | 2011-03-25 | 2012-10-03 | 通用电气公司 | 用于显示体积超声图像上的交叉信息的方法和系统 |
CN103006538A (zh) * | 2012-11-09 | 2013-04-03 | 叶琳 | 一种用于关节炎治疗的超声微泡及其用途 |
CN103530885A (zh) * | 2013-10-23 | 2014-01-22 | 北京倍肯恒业科技发展有限责任公司 | 一种一维图像自适应分层边缘检测提取算法 |
CN104398271A (zh) * | 2014-11-14 | 2015-03-11 | 西安交通大学 | 血管与斑块的三维力学及组织特性成像检测方法 |
CN104517116A (zh) * | 2013-09-30 | 2015-04-15 | 北京三星通信技术研究有限公司 | 在图像中确定目标区域的装置和方法 |
US20150208039A1 (en) * | 2014-01-21 | 2015-07-23 | Kabushiki Kaisha Toshiba | Medical image diagnostic apparatus, image processing apparatus, and image processing method |
CN105160957A (zh) * | 2015-07-09 | 2015-12-16 | 宋雪燕 | 一种基于ct影像模拟与定位的超声培训控制系统 |
CN106290579A (zh) * | 2016-08-03 | 2017-01-04 | 华东交通大学 | 基于双目标遗传算法和非劣分层的超声无损检测探头分布优化 |
CN106473772A (zh) * | 2015-08-31 | 2017-03-08 | 精工爱普生株式会社 | 超声波器件、超声波模块以及超声波测量仪 |
CN106580369A (zh) * | 2016-11-30 | 2017-04-26 | 珠海威泓医疗科技有限公司 | 一种多角度平面波相干彩色多普勒成像方法 |
CN107361791A (zh) * | 2017-07-21 | 2017-11-21 | 北京大学 | 一种快速超分辨血流成像方法 |
CN107537044A (zh) * | 2017-08-29 | 2018-01-05 | 重庆医科大学 | 一种壳聚糖纳米微泡及其制备方法与用途 |
CN107664641A (zh) * | 2017-10-10 | 2018-02-06 | 湖南科技大学 | 一种基于光照强度的工程陶瓷磨削损伤检测方法 |
CN107753062A (zh) * | 2017-11-27 | 2018-03-06 | 西安交通大学 | 基于马尔科夫链蒙特卡洛多目标追踪的经颅超声脑血管造影超分辨率成像方法 |
-
2018
- 2018-03-30 CN CN201810286743.8A patent/CN108836392B/zh active Active
Patent Citations (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1260730A (zh) * | 1997-06-19 | 2000-07-19 | 埃莱克特公司 | 自动化剂量规划法 |
US6645147B1 (en) * | 1998-11-25 | 2003-11-11 | Acuson Corporation | Diagnostic medical ultrasound image and system for contrast agent imaging |
JP2000287972A (ja) * | 1999-04-05 | 2000-10-17 | Toshiba Corp | 超音波診断装置および超音波診断装置の画像生成方法 |
US20030163040A1 (en) * | 2002-02-28 | 2003-08-28 | Philip Gildenberg | Audible feedback from positional guidance systems |
US20050090742A1 (en) * | 2003-08-19 | 2005-04-28 | Yoshitaka Mine | Ultrasonic diagnostic apparatus |
CN1947148A (zh) * | 2004-04-30 | 2007-04-11 | 埃莱克塔公共有限公司 | 用于自动提高医学图片的可用性的方法和系统 |
US20060241462A1 (en) * | 2005-03-14 | 2006-10-26 | Yi-Hong Chou | Method of intracranial ultrasound imaging and related system |
CN101184428A (zh) * | 2005-04-01 | 2008-05-21 | 视声公司 | 利用超声波使血管结构三维可视化的系统和方法 |
CN101291635A (zh) * | 2005-10-20 | 2008-10-22 | 直观外科手术公司 | 医用机器人系统中的计算机显示屏上的辅助图像显示和操纵 |
CN1792334A (zh) * | 2005-10-31 | 2006-06-28 | 西安交通大学 | 基于包膜微泡的灌注成像与超声控制释放的系统和方法 |
CN101011264A (zh) * | 2006-01-25 | 2007-08-08 | 株式会社东芝 | 超声波探头和超声波图像取得装置 |
CN101061963A (zh) * | 2006-04-27 | 2007-10-31 | 通用电气公司 | 用于测量通过心脏瓣膜血流的方法和系统 |
CN1903121A (zh) * | 2006-08-25 | 2007-01-31 | 北京源德生物医学工程有限公司 | 永磁磁共振图像导引体外高能聚焦超声系统和方法 |
CN101574267A (zh) * | 2008-05-08 | 2009-11-11 | 重庆医科大学超声影像学研究所 | 对超声微泡计数的方法 |
CN102317873A (zh) * | 2008-12-10 | 2012-01-11 | 赫勒拉德有限责任公司 | 用于彩色运动全息术的系统和方法 |
CN101577001A (zh) * | 2009-05-20 | 2009-11-11 | 电子科技大学 | 基于光线投射体绘制的三维医学图像剖分方法 |
CN101912273A (zh) * | 2010-07-23 | 2010-12-15 | 深圳市理邦精密仪器股份有限公司 | 超声动态图像处理方法及系统 |
CN202015163U (zh) * | 2010-12-30 | 2011-10-26 | 广州宝胆医疗器械科技有限公司 | 一种表层血管显示仪 |
CN102697523A (zh) * | 2011-03-25 | 2012-10-03 | 通用电气公司 | 用于显示体积超声图像上的交叉信息的方法和系统 |
CN103006538A (zh) * | 2012-11-09 | 2013-04-03 | 叶琳 | 一种用于关节炎治疗的超声微泡及其用途 |
CN104517116A (zh) * | 2013-09-30 | 2015-04-15 | 北京三星通信技术研究有限公司 | 在图像中确定目标区域的装置和方法 |
CN103530885A (zh) * | 2013-10-23 | 2014-01-22 | 北京倍肯恒业科技发展有限责任公司 | 一种一维图像自适应分层边缘检测提取算法 |
US20150208039A1 (en) * | 2014-01-21 | 2015-07-23 | Kabushiki Kaisha Toshiba | Medical image diagnostic apparatus, image processing apparatus, and image processing method |
CN104398271A (zh) * | 2014-11-14 | 2015-03-11 | 西安交通大学 | 血管与斑块的三维力学及组织特性成像检测方法 |
CN105160957A (zh) * | 2015-07-09 | 2015-12-16 | 宋雪燕 | 一种基于ct影像模拟与定位的超声培训控制系统 |
CN106473772A (zh) * | 2015-08-31 | 2017-03-08 | 精工爱普生株式会社 | 超声波器件、超声波模块以及超声波测量仪 |
CN106290579A (zh) * | 2016-08-03 | 2017-01-04 | 华东交通大学 | 基于双目标遗传算法和非劣分层的超声无损检测探头分布优化 |
CN106580369A (zh) * | 2016-11-30 | 2017-04-26 | 珠海威泓医疗科技有限公司 | 一种多角度平面波相干彩色多普勒成像方法 |
CN107361791A (zh) * | 2017-07-21 | 2017-11-21 | 北京大学 | 一种快速超分辨血流成像方法 |
CN107537044A (zh) * | 2017-08-29 | 2018-01-05 | 重庆医科大学 | 一种壳聚糖纳米微泡及其制备方法与用途 |
CN107664641A (zh) * | 2017-10-10 | 2018-02-06 | 湖南科技大学 | 一种基于光照强度的工程陶瓷磨削损伤检测方法 |
CN107753062A (zh) * | 2017-11-27 | 2018-03-06 | 西安交通大学 | 基于马尔科夫链蒙特卡洛多目标追踪的经颅超声脑血管造影超分辨率成像方法 |
Non-Patent Citations (2)
Title |
---|
YOSHIKI YAMAKOSHI: "《Measurement of Secondary Ultrasonic Waves from Microbubbles by Holographic Image Reconstruction》", 《JAPANESE JOURNAL OF APPLIED PHYSICS》 * |
吴心云: "《超声微泡及其在分子影像诊断和靶向治疗中的应用》", 《现代诊断和治疗》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113180734A (zh) * | 2018-12-27 | 2021-07-30 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声血流成像方法及系统 |
CN110327077B (zh) * | 2019-07-09 | 2022-04-15 | 深圳开立生物医疗科技股份有限公司 | 一种血流显示方法、装置及超声设备和存储介质 |
CN110327077A (zh) * | 2019-07-09 | 2019-10-15 | 深圳开立生物医疗科技股份有限公司 | 一种血流显示方法、装置及超声设备和存储介质 |
CN110772285A (zh) * | 2019-10-31 | 2020-02-11 | 南京景瑞康分子医药科技有限公司 | 一种超声超分辨成像方法 |
CN110772285B (zh) * | 2019-10-31 | 2022-05-17 | 南京景瑞康分子医药科技有限公司 | 一种超声超分辨成像方法 |
US20220309613A1 (en) * | 2020-05-18 | 2022-09-29 | Nanjing Transcend Vivoscope Bio-Technology Co., Ltd | Super-resolution reconstruction method and apparatus for three-dimensional contrast-enhanced ultrasound images |
WO2021233188A1 (zh) * | 2020-05-18 | 2021-11-25 | 南京超维景生物科技有限公司 | 三维超声造影图像的超分辨重建方法及装置 |
CN114073547B (zh) * | 2020-08-14 | 2024-06-07 | 深圳迈瑞生物医疗电子股份有限公司 | 用于超声成像的方法、装置和存储介质 |
CN114073547A (zh) * | 2020-08-14 | 2022-02-22 | 深圳迈瑞生物医疗电子股份有限公司 | 用于超声成像的方法、装置和存储介质 |
CN112043308A (zh) * | 2020-08-31 | 2020-12-08 | 深圳市德力凯医疗设备股份有限公司 | 一种颅内三维脑血流重建方法、存储介质及超声设备 |
CN112370073A (zh) * | 2020-10-12 | 2021-02-19 | 深圳华声医疗技术股份有限公司 | 超声声速校正方法、超声信号处理设备及存储介质 |
CN112869767A (zh) * | 2021-01-11 | 2021-06-01 | 青岛海信医疗设备股份有限公司 | 一种超声图像存储方法、装置及其超声设备 |
CN112932539A (zh) * | 2021-01-23 | 2021-06-11 | 西安交通大学 | 一种经颅超声估计血管血流速度场的装置及方法 |
CN114129190A (zh) * | 2021-11-30 | 2022-03-04 | 深圳先进技术研究院 | 全眼微血流的成像方法、装置及电子设备和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN108836392B (zh) | 2021-06-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108836392A (zh) | 基于超声rf信号的超声成像方法、装置、设备及存储介质 | |
CN107361791B (zh) | 一种快速超分辨血流成像方法 | |
US7744533B2 (en) | Ultrasonic diagnostic apparatus, image processing apparatus and image processing method | |
US8202220B2 (en) | Ultrasound diagnostic apparatus | |
JP7139357B2 (ja) | 超音波を用いた組織微小脈管構造の可視化のためのシステム及び方法 | |
CN109803588A (zh) | 流体路径的增强分辨率超声成像 | |
US20100324420A1 (en) | Method and System for Imaging | |
US20090024033A1 (en) | Ultrasound diagnostic apparatus | |
Jensen et al. | Fast super resolution ultrasound imaging using the erythrocytes | |
Yan et al. | Transthoracic ultrasound localization microscopy of myocardial vasculature in patients | |
Taghavi et al. | In vivo 3D super-resolution ultrasound imaging of a rat kidney using a row-column array | |
Jensen et al. | In vivo super resolution ultrasound imaging using the erythrocytes-SURE | |
CN113573645B (zh) | 用于调整超声探头的视场的方法和系统 | |
Rohling et al. | Issues in 3-D free-hand medical ultrasound imaging | |
Wang et al. | Synthetic aperture focusing technique for 3-D ultrasound computed tomography | |
Yociss et al. | Null subtraction beamforming for improved vessel resolution in volumetric contrast-enhanced ultrasound | |
JP3702054B2 (ja) | 磁気共鳴映像装置 | |
US20240268793A1 (en) | Super resolution ultrasound imaging | |
US20220386999A1 (en) | Super Resolution Ultrasound Imaging | |
Kanoulas et al. | Realistic Super-Resolution Image Analysis for State of the Art 2D Contrast Enhanced Ultrasound Imaging | |
Koivukangas et al. | Completely computer‐focused ultrasound imaging. First clinical imaging results. | |
JP2024519328A (ja) | 超解像度超音波イメージング | |
Chabouh et al. | RF Channel-Based Adaptive Beamforming for 3D Ultrasound Localization Microscopy | |
Bureau et al. | Ultrasound matrix imaging for transcranial in-vivo localization microscopy | |
Lin et al. | Power Doppler: how it works, its clinical benefits, and recent technologic advances |
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 |