CN106063710A - 一种超声成像方法、系统及装置 - Google Patents

一种超声成像方法、系统及装置 Download PDF

Info

Publication number
CN106063710A
CN106063710A CN201610353267.8A CN201610353267A CN106063710A CN 106063710 A CN106063710 A CN 106063710A CN 201610353267 A CN201610353267 A CN 201610353267A CN 106063710 A CN106063710 A CN 106063710A
Authority
CN
China
Prior art keywords
signal
point
subimage
focus
focus point
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
Application number
CN201610353267.8A
Other languages
English (en)
Other versions
CN106063710B (zh
Inventor
潘晓畅
唐果
石学工
冯乃章
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sonoscape Medical Corp
Original Assignee
Sonoscape Medical Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Sonoscape Medical Corp filed Critical Sonoscape Medical Corp
Priority to CN201610353267.8A priority Critical patent/CN106063710B/zh
Publication of CN106063710A publication Critical patent/CN106063710A/zh
Priority to PCT/CN2017/085888 priority patent/WO2017202362A1/zh
Application granted granted Critical
Publication of CN106063710B publication Critical patent/CN106063710B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207Devices 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5238Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
    • A61B8/5246Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from the same or different imaging techniques, e.g. color Doppler and B-mode
    • A61B8/5253Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from the same or different imaging techniques, e.g. color Doppler and B-mode combining overlapping images, e.g. spatial compounding

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

本发明提供了一种超声成像方法、系统及装置,该方法包括:预设成像的区域,且预设多个接收超声波的聚焦点以及对应聚焦点的发射阵元的延时,聚焦点设置在成像区域下方;根据每个聚焦点对应的延时激活超声探头的发射阵元,以聚焦的方式采用为探头宽度的30%‑100%的孔径向聚焦点发射超声波;确定成像区域内的采样点与有效发射声场的重叠部分;计算成像区域内包含有效发射声场的采样点的散射强度,得到与聚焦点对应的重建后的子图像;将几次发射后重建的子图像按照像素的位置加权叠加生成一帧图像。本发明通过该方法提高了图像帧频,最少两次发射超声波就可以得到一幅超声图像,不同的聚焦点重建的子图像的叠加能显著减小散斑尺寸,提高了图像分辨率。

Description

一种超声成像方法、系统及装置
技术领域
本发明涉及医疗诊断超声成像技术领域,更具体的,本发明涉及一种超声成像方法、系统及装置。
背景技术
现代超声医学诊断对图像的实时性,即图像帧频的要求越来越高。早期的单波束成像技术发射一次仅能合成一条波束,需要发射许多次才能生成一帧超声图像;而现有的多波束技术通过一次发射即可同时合成m条波束(m通常为2,3或4),如图1所示。对于包含相同声束线数的图像,多波束技术的发射次数仅为早期单波束技术的1/m,从而其成像的速度能提高m倍。目前,多波束技术是现有的提高超声图像帧频的主要手段。
现代超声成像系统的多波束成像的流程图如图2所示。m波束技术的波束合成在每次超声发射后,能同时合成m条波束。对于每一次超声发射,采集所有阵元上的通道数据,按要合成的m条波束上的每个点进行延时和变迹,然后通过求和运算得到m条波束的射频信号;接着对射频信号进行包络检测,提取出基带信号,即包络信号;最后按照B超图像上所需的轴向分辨率所确定的成像区域所需的轴向采样点的个数,对包络信号进行下采样,即可得B超图像上的m条线;改变发射焦点的位置,再经过发射和波束合成后又能得到B超图像上的m条线,经过k次发射后一共能得到k×m条线,将所有的线拼接在一起,就得到了一帧B超图像。对于波束合成中的延时计算,多波束上的每个点的延时计算如图3所示,波束上的任意一点P相对于通道A的延时为其中为点P到探头的距离,为点P到通道A的距离,sound_speed为超声波在组织中的传播速度。
尽管多波束技术能显著提高图像的帧频,但其图像的分辨率却比单波束的图像要差,目前为了提高图像的分辨率,需要进行多次的聚焦发射,从而在一定程度上提高图像的分辨率,但是又会降低图像的帧频。
发明内容
本发明提供了一种超声成像方法、系统及装置,以解决现有技术中为了提高图像的分辨率,需要进行多次的聚焦发射,从而在一定程度上提高图像的分辨率,但是又会降低图像的帧频的问题。
为达到上述目的,本发明提供了以下技术方案:
本发明公开了一种超声成像方法,该方法包括:
预设成像的区域,在成像区域内进行离散采样,且预设多个接收超声波的聚焦点以及对应所述聚焦点的发射阵元的延时,所述聚焦点设置在成像区域下方;
根据每个所述聚焦点对应的所述延时激活超声探头的所述发射阵元,以聚焦的方式采用为探头宽度的30%-100%的孔径向所述聚焦点发射超声波;
确定所述成像区域内的采样点与有效发射声场采样点的重叠部分,获取所述超声波的所述有效发射声场内的采样点;
计算所述成像区域内包含所述有效发射声场的采样点的散射强度,得到与所述聚焦点对应的重建后的子图像;
将几次发射后重建的所述子图像按照像素的位置加权叠加生成一帧图像。
优选的,所述计算所述成像区域内包含所述有效发射声场的采样点的散射强度,得到与所述聚焦点对应的重建后的子图像具体为:
采集各个所述发射阵元对应的各个通道数据;
分别对所述通道数据进行包络检测,得到包络信号,并对所述包络信号进行降采样、聚焦延时、变迹以及求和,得到子图像内每个采样点的散射强度;
将所述不同聚焦点发射后重建的子图像按采样点的位置进行叠加,得到所述聚焦点对应的散射强度。
优选的,所述计算所述成像区域内包含所述有效发射声场的采样点的散射强度,得到与所述聚焦点对应的重建后的子图像具体为:
采集各个所述发射阵元对应的各个通道数据;
分别对所述通道数据进行IQ解调,得到IQ解析信号,所述解析信号包括I信号和Q信号,分别对所述I信号和所述Q信号进行降采样、聚焦延时、变迹、求和以及,得到子图像内每个采样点的复数信号;
对所述子图像内每个采样点的复数信号进行IQ取模得到所述聚焦点对应的散射强度。
其中,所述IQ解调为希尔伯特变换或数字正交解调。
在上述公开的方法基础上,本发明还公开了一种超声成像系统,该系统包括:
预设单元,用于预设成像的区域,在成像区域内进行离散采样,且预设多个接收超声波的聚焦点以及对应所述聚焦点的发射阵元的延时,所述聚焦点设置在成像区域下方;
发射单元,用于根据每个所述聚焦点对应的所述延时激活超声探头的所述发射阵元,以聚焦的方式采用为探头宽度的30%-100%的孔径向所述聚焦点发射超声波;
确定单元,用于确定所述成像区域内的采样点与有效发射声场的重叠部分,并获取所述超声波的所述有效发射声场内的采样点;
计算单元,用于计算所述成像区域内包含所述有效发射声场的采样点的散射强度,得到与所述聚焦点对应的重建后的子图像;
叠加单元,用于将几次发射后重建的所述子图像按照像素的位置加权叠加生成一帧图像。
优选的,所述计算单元具体为:
第一采集单元,用于采集各个所述发射阵元对应的各个通道数据;
第一计算单元,用于分别对所述通道数据进行包络检测,得到包络信号,并对所述包络信号进行降采样、聚焦延时、变迹和求和,得到子图像内每个采样点的散射强度;
第一叠加单元,用于将所述不同聚焦点发射后重建的子图像按采样点的位置进行叠加,得到所述聚焦点对应的散射强度。
优选的,所述计算单元具体为:
第二采集单元,用于采集各个所述发射阵元对应的各个通道数据;
第二计算单元,用于分别对所述通道数据进行IQ解调,得到IQ解析信号,所述解析信号包括I信号和Q信号,分别对所述I信号和所述Q信号进行降采样、聚焦延时、变迹以及求和,得到子图像内每个采样点的复数信号;
IQ取模单元,用于对所述子图像内每个采样点的复数信号进行IQ取模得到所述聚焦点对应的的散射强度。
其中,所述IQ解调为希尔伯特变换或数字正交解调。
本发明在上述公开的方法及系统的基础上,本发明还公开了一种超声成像装置,包括:超声换能器,前端发射和接收模块,前端处理模块以及超声主机系统,其中,
所述超声换能器包括多个阵元,用于发射超声波;
所述的前端发射和接收模块用于控制超声波的发射及接收超声回波信号;
所述前端处理模块用于对接收到的所述超声回波信号的预处理;
所述超声主机系统用于超声图像的重建。
经由上述的技术方案可知,与现有技术相比,本发明公开了一种超声成像方法、系统及装置,该方法包括:预设成像的区域,在成像区域内进行离散采样,且预设多个接收超声波的聚焦点以及对应聚焦点的发射阵元的延时,聚焦点设置在成像区域下方;根据每个聚焦点对应的延时激活超声探头的发射阵元,以聚焦的方式采用为探头宽度的30%-100%的孔径向聚焦点发射超声波;确定成像区域内的采样点与有效发射声场的重叠部分;计算成像区域内包含有效发射声场的采样点的散射强度,得到与聚焦点对应的重建后的子图像;将几次发射后重建的子图像按照像素的位置加权叠加生成一帧图像。本发明通过该方法提高了图像帧频,最少两次发射超声波就可以得到一幅超声图像,不同的聚焦点重建的子图像的叠加能显著减小散斑尺寸,提高了图像分辨率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1为现有技术中三波束技术原理示意图;
图2A-2B为现有技术中多波束合成方法的流程图;
图3为现有技术中波束合成的延时计算方式的示意图;
图4为本发明中单次发射生成一幅子图像的示意图;
图5为本发明实施例公开的一种超声成像方法的流程图;
图6为本发明实施例中超分辨率B超图像重建的流程图;
图7为图5中的步骤S504的一种具体实施方法流程图;
图8为本发明实施例中高分辨率B超图像重建的流程图;
图9为图5中的步骤S504的另一种具体实施方法流程图;
图10为本发明实施例中高分辨率B超图像重建的流程图;
图11为多幅子图像重建出B超图像的原理示意图;
图12A、12B展示了本发明中子图像叠加改善图像分辨率的示意图;
图13为本发明实施例公开的一种超声成像系统的结构示意图;
图14为图13中的U4的一种具体实施结构的示意图;
图15为图13中的U4的另一种具体实施结构的示意图;
图16为本发明实施例公开的一种超声成像装置的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
影响B超图像分辨率的主要的因素是成像系统的PSF(Point spread function,点扩散函数)的主瓣宽度,超声成像系统的点扩散函数对B超图像分辨率的影响与光学系统的点扩散函数对光学图像分辨率的影响相同,点扩散函数的主瓣越窄,相邻的两点在图像中就越容易分开(无论是对于光学图像还是超声散斑图像),从而图像的分辨率会更高。因此,提高成像系统所成的图像的分辨率的本质是减小成像系统点扩散函数的主瓣宽度。
请参阅附图5,为本发明实施例公开的一种超声成像方法的流程图。本发明公开了一种超声成像方法,如图5所示,该方法具体步骤包括:
S501、预设成像的区域,在成像区域内进行离散采样,且预设多个接收超声波的聚焦点以及对应聚焦点的发射阵元的延时,聚焦点设置在成像区域下方。
如图4所示,本发明的成像系统的超声发射方式,用不同的延时激励探头的发射阵元E1,各个阵元所发射超声波聚焦在焦点C处。
S502、根据每个聚焦点对应的延时激活超声探头的发射阵元,以聚焦的方式采用为探头宽度的30%-100%的孔径向聚焦点发射超声波。
本发明的发射方式采用大孔径发射,超声发射时的发射孔径占探头宽度的30%-100%。现有超声技术发射时的孔径通常较小,对于线阵探头,发射孔径通常小于探头宽度的30%。还有,本发明的成像系统在超声发射时将发射焦点C设置在成像区域R1的下方,如图4所示。本发明的大孔径发射和发射焦点设置在成像区域下方的目的是为了扩大有效声场的覆盖区域。
S503、确定成像区域内的采样点与有效发射声场的重叠部分,并获取超声波的有效发射声场内的采样点。
S504、计算成像区域内包含有效发射声场的采样点的散射强度,得到与聚焦点对应的重建后的子图像。
S505、将几次发射后重建的子图像按照像素的位置加权叠加生成一帧图像。
具体的,请参阅附图6,如图6所示,每次发射后进行子图像重建111a-111n,可得到不同的子图像113a-113n。最终图像重建过程115是将将所有的子图像113a-113n按像素的位置加权叠加在一起,最后就能重建出一帧高分辨率的B超图像。
本发明提供了一种超声成像方法,该方法包括:预设成像的区域,在成像区域内进行离散采样,且预设多个接收超声波的聚焦点以及对应聚焦点的发射阵元的延时,聚焦点设置在成像区域下方;根据每个聚焦点对应的延时激活超声探头的发射阵元,以聚焦的方式采用为探头宽度的30%-100%的孔径向聚焦点发射超声波;确定成像区域内的采样点与有效发射声场采样点的重叠部分,获取超声波的有效发射声场采样点;计算成像区域内包含有效发射声场的采样点的散射强度,得到与聚焦点对应的重建后的子图像;将几次发射后重建的子图像按照像素的位置加权叠加生成一帧图像。本发明通过该方法提高了图像帧频,最少两次发射超声波就可以得到一幅超声图像,不同的聚焦点重建的子图像的叠加能显著减小散斑尺寸,提高了图像分辨率。
下面以一个实例说明本发明的超声发射方式和现有超声技术的超声发射方式之间的区别。对于一个宽度为38mm的线阵探头,假设成像区域的尺寸为38×38mm(高度×宽度)。成像区域内的横向采样点数设为320,即波束合成后的超声图像的线数为320。对于本发明的超声成像技术,超声发射采用发射孔径为38mm,即所激活发射的阵元宽度为38mm。发射两次,发射焦点分别设置在(-19,76)和(19,76)处。则每次发射有效声场在所成像的区域内占的面积为3/4,仅通过两次发射就全部覆盖了所成像的区域,且还有1/2的成像区域的重合。对于现有的多波束技术,假设为四波束,发射的孔径为5mm,发射焦点的深度设置在图像中间,即19mm处。对于现有的多波束成像技术,一共需要发射80次,才能覆盖线数为320的成像区域。这样对比本发明的技术和现有多波束成像技术,本发明的超声发射次数仅为现有的多波束成像技术的1/40,理论上能将成像的速度提高40倍。
具体的,请参阅附图7,图7为图5中步骤S504的一种具体实施方法流程图。具体的步骤S504包括:
步骤S701、采集各个发射阵元对应的各个通道数据。
步骤S702、分别对通道数据进行包络检测,得到通道数据的包络信号。
步骤S703、对包络信号进行降采样、聚焦延时、变迹以及求和,得到子图像区域内每个采样点的散射强度。
步骤S704、将不同聚焦点发射后重建的子图像按采样点的位置进行求和叠加,得到聚焦点对应的散射强度。
上述步骤S504的具体步骤方法可以参阅附图8,如图8所示,本发明的图像重建技术发射超声波后,采集各个阵元所接收到的通道数据101a-101n。接着同时对这些通道信号进行包络检测103a-103n操作,得到每个通道信号的幅度信息。对于每次发射所要重建的子图像区域,是由有效的发射声场在成像区域内所覆盖的位置(图3)所决定的。对成像区域内包含有效发射声场的区域进行离散/采样,得到所要重建子图像内每个像素点/采样点的坐标位置。本发明流程图中的降采样,延时和变迹操作105a-105n的实现过程为:计算子图像区域内的所有采样点在通道数据中的延时,并根据这些采样点相对每个探头阵元的相对位置进行变迹。接着对每个采样点对应的延时和变迹后的所有通道数据求和操作107,这样就得到了每个采样点的散射系数。计算出所有采样点的散射系数后即重建出了一幅子图像109。
具体的,请参阅附图9,图9为图5中步骤S504的一种具体实施方法流程图。具体的步骤S504包括:
S901、采集各个发射阵元对应的各个通道数据。
S902、分别对通道数据进行IQ解调,得到IQ解析信号,解析信号包括I信号和Q信号;
S903、分别对I信号和Q信号进行降采样、聚焦延时、变迹以及求和,得到子图像内每个采样点的复数信号;
S904、对子图像内各个采样点的复数信号进行IQ取模得到聚焦点对应的散射强度。
上述步骤S504的具体步骤方法可以参阅附图10,如图10所示,图10是将包络检测操作103a-103n直接放在通道数据101a-101n后,且将降采样,延时和变迹等步骤105a-105n合并在一起操作,这样做的好处之一是成像过程中需要生成射频信号。对于早期的单波束成像和现代的多波束成像技术,合成射频信号是成像过程中必不可少的环节。如图2A所示,现有多波束成像技术是先对所合成的每条线上的每个点在进行延时和变迹22a-22n计算,然后经过求和24步骤,即可得到这条声束线所对应的射频信号26。现有多波束技术生成的射频信号的采样率必须满足乃奎斯特采样定律,否则会丢失信息,包络检测后得到的信号就会出错。然而最后得到的B超图像的轴向方向并不需要和射频信号一样高的采样率,所以现有技术中对射频信号26进行包络检测28后,得到的包络信号30需要再经过一次降采样32。本发明的成像技术将降采样环节放在前面,免掉了生成射频信号的环节,从而减少了许多不需要的采样点的散射强度的计算。另外,本发明的成像技术在重建一帧图像时,仅仅需要保存一系列采样率较低的子图像,而不需要保存采样率较高的原始数据(如通道数据)。
具体的,所述IQ解调为希尔伯特变换或数字正交解调;所述聚焦延时采用焦点延时,其计算公式为其中P代表子图像内的任意一个采样点,为聚焦点C到探头的距离,为点P到聚焦点C之间的距离,为点P到通道A的距离,sound_speed为超声波在组织中的传播速度。
本发明的成像技术所采用的延时为焦点延时,这种延时与传统延时的不同是还需要考虑发射焦点的位置。焦点延时的计算示意图如图11所示,图像内任意一点P相对于通道A的延时为其中为发射焦点C到探头的距离,为点P到发射焦点C之间的距离,为点P到通道A的距离,sound_speed为超声波在组织中的传播速度。不同组织中超声波传播速度是有一定的差异的,本发明中延时计算所采用的超声波在组织中的传播速度不是某一固定的值,而是可变的,该速度值会根据所成像组织中的实际超声传播速度而进行调整。
本发明在数字超声中包络解调的实现方式有两大类,一类是绝对值低通滤波法,另一类是IQ解调法。绝对值低通滤波法计算信号的包络的过程为,先对信号取绝对值,然后再通过低通滤波去掉高频成分。现代数字超声系统中有两种方式实现IQ解调法,一种是数字正交解调,另一种是希尔伯特变换。数字正交解调的原理是先对信号进行移频,然后将移频后的信号进行低通滤波,最后将滤波后的信号按反方向移频回去。IQ解调后的信号的实数部分为同相位(In-phase)信号,通常称为I信号;信号的虚数部分为正交相位(Quadrature-phase)信号,通常称为Q信号。除了数字正交解调法之外,希尔伯特变换也能计算出IQ信号。原信号是I信号,对原信号进行希尔伯特变换即可直接得到正交相位的信号,即Q信号。通过对I信号和Q信号进行取模运算即可得到原信号的幅度信息,即包络信号。
对于图10所展示的是实现本发明成像的另一种方案,与图8不同的是,先对通道数据进行IQ解调104a-104n,然后对IQ信号都进行降采样,延时,变迹和求和,最后进行IQ取模108得到子图像。图10的方案与图8的方案相比计算量更大,在降采样,延时,变迹105a-105n和求和107处理这几个操作的计算量增加了一倍。图10的方案的优势在于求和107后得到的IQ数据还能用于血流计算和弹性成像等其他高级成像功能。
图12示出的是多幅子图像重建出高分辨率图像的原理示意图。下面用一个实例说明本发明所提出的超分辨率成像系统是如何提高图像的分辨率的。图12A所示的是一次超声发射后成像区域内某一处的点扩散函数1110的形状示意图。点扩散函数的主瓣形状是各向异性的,其中对应发射焦点1112方向的点扩散函数的主瓣的轴向宽度窄,而点扩散函数的主瓣的横向宽度宽。点扩散函数的主瓣的轴向宽度要小于主瓣的横向宽度,这是B超图像中散斑为横向拉伸的形态,也是B超图像的轴向分辨率要高于横向分辨率的主要原因。
图12A的示意图表示,对于本发明的成像方法,每一次超声发射后所重建得到的子图像与发射方向一致的方向的分辨率高,而与发射方向垂直的方向的分辨率低。图12B的示意图反映的是,在所有不同的发射焦点117a-117e上进行超声发射后,重建的子图像上同一区域上不同的点扩散函数叠加1115的结果。不同子图像的点扩散函数在中心的叠加是相干(Coherent)叠加,所以点扩散函数的主瓣的中心会增强。点扩散函数中心之外的叠加是非相干(Incoherent)的叠加,点扩散函数的主瓣宽度会因为这种非相干的叠加而减弱变小。最后的结果是叠加后成像系统的点扩散函数的中心更突出,离中心越远的主瓣强度越低,效果就是主瓣宽度变窄。本发明的超分辨率成像系统通过子图像重建的方式,有效地减小了点扩散函数的主瓣宽度,从而提高了B超图像的横向分辨率。
不仅如此,子图像的叠加还能降低图像中出现的随机噪声的幅度,从而提高图像的信噪比。子图像的叠加相干地提高了散斑信号的质量,同时减少了噪声的幅度,最后叠加的图像的信噪比能得到提高,图像中更深区域的信噪比的提高意味着成像深度的提高。此外,现有超声技术成像生成的图像中,焦点区域的声场强度要显著高于其他区域,从而图像焦点区域的亮度要比其他区域更高。本发明得到的图像在声场强度高的区域叠加的少,在声场强度低的区域叠加得多,从而图像内的最终声场强度分布更加均匀,所呈现的图像的一致性要更好。
请参阅附图13,图13为本发明实施例公开的一种超声成像系统的结构示意图。如图13所示,本发明公开了一种超声成像系统,该系统基于快速成像技术和图像重建技术,该系统具体包括:预设单元U1、发射单元U2、确定单元U3、计算单元U4和叠加单元U5,其中:
预设单元U1,用于预设成像的区域,在成像区域内进行离散采样,且预设多个接收超声波的聚焦点以及对应聚焦点的发射阵元的延时,聚焦点设置在成像区域下方;发射单元U2,用于根据每个聚焦点对应的延时激活超声探头的发射阵元,以聚焦的方式采用为探头宽度的30%-100%的孔径向聚焦点发射超声波;确定单元U3,用于确定成像区域内的采样点与有效发射声场的重叠部分,并获取超声波的有效发射声场内的采样点;计算单元U4,用于计算成像区域内包含有效发射声场的采样点的散射强度,得到与聚焦点对应的重建后的子图像;叠加单元U5,用于将几次发射后重建的子图像按照像素的位置加权叠加生成一帧图像。
由于本实施例中的各单元能够执行图5所示的方法,本实施例未详细描述的部分,可参考对图5的相关说明。
本发明提供了一种超声成像系统,本发明通过该方法提高了图像帧频,最少两次发射超声波就可以得到一幅超声图像,不同的聚焦点重建的子图像的叠加能显著减小散斑尺寸,提高了图像分辨率。
具体的,请参阅附图14,图14为图13中计算单元的一种具体结构示意图。如图14所示,所述计算单元U4具体为:
第一采集单元U6,用于采集各个发射阵元对应的各个通道数据;
第一检测单元U7,用于分别对通道数据进行包络检测,得到通道数据的包络信号;
第一计算单元U8,用于对包络信号进行降采样、聚焦延时、变迹以及求和,得到子图像内每个采样点的散射强度;
第一叠加单元U9,用于用于将不同聚焦点发射后重建的子图像按采样点的位置进行叠加,得到聚焦点对应的散射强度。
由于本实施例中的各单元能够执行图7所示的方法,本实施例未详细描述的部分,可参考对图7的相关说明。
具体的,请参阅附图15,图15为图12中计算单元的另一种具体结构示意图。如图15所示,所述计算单元U4具体为:
第二采集单元U10,用于采集各个发射阵元对应的各个通道数据;
解调解析单元U11,用于分别对通道数据进行IQ解调,得到IQ解析信号,解析信号包括I信号和Q信号;
第二计算单元U12,用于分别对I信号和所述Q信号进行降采样、聚焦延时、变迹和求和,得到子图像内的每个采样点的复数信号;
IQ取模单元U13,用于对各个子图像内的每个采样点的复数信号进行IQ取模得到聚焦点对应子图像内每个采样点的散射强度。
由于本实施例中的各单元能够执行图9所示的方法,本实施例未详细描述的部分,可参考对图9的相关说明。
本发明在上述公开的方法及系统的基础上,本发明还公开了一种超声成像装置,具体的,请参阅附图16,图16为本发明实施例公开的一种超声成像装置的结构示意图。如图16所示,该超声成像装置具体包括:超声换能器P1,前端发射和接收模块P2,前端处理模块P3以及超声主机系统P4,其中,超声换能器P1包括多个阵元,用于发射超声波;前端发射和接收模块P2用于控制超声波的发射及接收超声回波信号;前端处理模块P3用于对接收到的超声回波信号的预处理;超声主机系统P4用于超声图像的重建。
本发明提供了一种超声成像装置,通过该装置提高了图像帧频,最少两次发射超声波就可以得到一幅超声图像,不同的聚焦点重建的子图像的叠加能显著减小散斑尺寸,提高了图像分辨率。
综上所述,本发明提供了一种超声成像方法、系统及装置,该方法包括:预设多个接收超声波的聚焦点以及对应聚焦点的发射阵元的延时发射信号,聚焦点设置在成像区域下方;根据每个聚焦点对应的延时发射信号激活超声探头的发射阵元,以聚焦的方式采用为探头宽度的30%-100%的孔径向聚焦点发射超声波;获取超声波的有效发射声场采样点,确定成像区域与有效发射声场采样点的重叠部分;计算成像区域内包含有效发射声场采样点的散射强度,得到聚焦点对应的重建后的子图像;将子图像按照像素的位置加权叠加生成一帧图像。本发明通过该方法提高了图像帧频,最少两次发射超声波就可以得到一幅超声图像,不同的聚焦点重建的子图像的叠加能显著减小散斑尺寸,提高了图像分辨率。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的系统,装置和单元的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
在本发明所提供的几个实施例中,应该理解到,所揭露的系统,装置和方法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如,多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
需要说明的是,本说明书中的各个实施例均采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似的部分互相参见即可。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (10)

1.一种超声成像方法,其特征在于,该方法包括:
预设成像的区域,在成像区域内进行离散采样,且预设多个接收超声波的聚焦点以及对应所述聚焦点的发射阵元的延时,所述聚焦点设置在成像区域下方;
根据每个所述聚焦点对应的所述延时激活超声探头的所述发射阵元,以聚焦的方式采用为探头宽度的30%-100%的孔径向所述聚焦点发射超声波;
确定所述成像区域内的采样点与有效发射声场的重叠部分,并获取所述超声波的所述有效发射声场内的采样点;
计算所述成像区域内包含所述有效发射声场的采样点的散射强度,得到与所述聚焦点对应的重建后的子图像;
将几次发射后重建的所述子图像按照像素的位置加权叠加生成一帧图像。
2.根据权利要求1所述的方法,其特征在于,所述计算所述成像区域内包含所述有效发射声场的采样点的散射强度,得到与所述聚焦点对应的重建后的子图像具体为:
采集各个所述发射阵元对应的各个通道数据;
分别对所述通道数据进行包络检测,得到通道数据的包络信号;
对所述包络信号进行降采样、聚焦延时、变迹以及求和,得到子图像区域内每个采样点的散射强度;
将所述不同聚焦点发射后重建的子图像按采样点的位置进行叠加,得到所述聚焦点对应的散射强度。
3.根据权利要求1所述的方法,其特征在于,所述计算所述成像区域内包含所述有效发射声场采样点的散射强度,得到所述聚焦点对应的重建后的子图像具体为:
采集各个所述发射阵元对应的各个通道数据;
分别对所述通道数据进行IQ解调,得到IQ解析信号,所述解析信号包括I信号和Q信号;
分别对所述I信号和所述Q信号进行降采样、聚焦延时、变迹以及求和,得到所述子图像内每个采样点的复数信号;
对所述子图像内每个采样点的复数信号进行IQ取模得到所述聚焦点对应的散射强度。
4.根据权利要求3所述的方法,其特征在于,所述IQ解调为希尔伯特变换或数字正交解调。
5.根据权利要求2或3所述的方法,其特征在于,所述聚焦延时采用焦点延时,其计算公式为其中P代表子图像内的任意一个采样点,为聚焦点C到探头的距离,为点P到聚焦点C之间的距离,为点P到通道A的距离,sound_speed为超声波在组织中的传播速度。
6.一种超声成像系统,其特征在于,该系统包括:
预设单元,用于预设成像区域,对成像区域进行离散采样,且预设多个接收超声波的聚焦点以及对应所述聚焦点的发射阵元的延时,所述聚焦点设置在成像区域下方;
发射单元,用于根据每个所述聚焦点对应的所述延时激活超声探头的所述发射阵元,以聚焦的方式采用为探头宽度的30%-100%的孔径向所述聚焦点发射超声波;
确定单元,用于确定所述成像区域内的采样点与有效发射声场的重叠部分,并获取所述超声波的所述有效发射声场内的采样点;
计算单元,用于计算所述成像区域内包含所述有效发射声场的采样点的散射强度,得到与所述聚焦点对应的重建后的子图像;
叠加单元,用于将几次发射后重建的所述子图像按照像素的位置加权叠加生成一帧图像。
7.根据权利要求6所述的系统,其特征在于,所述计算单元具体为:
第一采集单元,用于采集各个所述发射阵元对应的各个通道数据;
第一检测单元,用于分别对所述通道数据进行包络检测,得到所述通道数据的包络信号;
第一计算单元,用于对所述包络信号进行降采样、聚焦延时、变迹以及求和,得到子图像区域内每个采样点的散射强度;
第一叠加单元,用于将所述不同聚焦点发射后重建的子图像按采样点的位置进行叠加,得到所述聚焦点对应的散射强度。
8.根据权利要求6所述的系统,其特征在于,所述计算单元具体为:
第二采集单元,用于采集各个所述发射阵元对应的各个通道数据;
解调解析单元,用于分别对所述通道数据进行IQ解调,得到IQ解析信号,所述解析信号包括I信号和Q信号;
第二计算单元,用于分别对所述I信号和所述Q信号进行降采样、聚焦延时、变迹以及求和,得到子图像内每个采样点的复数信号;
IQ取模单元,用于对所述子图像内每个采样点的复数信号进行IQ取模得到所述聚焦点对应的散射强度。
9.根据权利要求8所述的方法,其特征在于,所述IQ解调为希尔伯特变换或数字正交解调。
10.一种超声成像装置,其特征在于,包括:超声换能器,前端发射和接收模块,前端处理模块以及超声主机系统,其中,
所述超声换能器包括多个阵元,用于发射超声波;
所述的前端发射和接收模块用于控制超声波的发射及接收超声回波信号;
所述前端处理模块用于对接收到的所述超声回波信号的预处理;
所述超声主机系统用于超声图像的重建。
CN201610353267.8A 2016-05-25 2016-05-25 一种超声成像方法、系统及装置 Active CN106063710B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201610353267.8A CN106063710B (zh) 2016-05-25 2016-05-25 一种超声成像方法、系统及装置
PCT/CN2017/085888 WO2017202362A1 (zh) 2016-05-25 2017-05-25 一种超声成像方法、系统及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610353267.8A CN106063710B (zh) 2016-05-25 2016-05-25 一种超声成像方法、系统及装置

Publications (2)

Publication Number Publication Date
CN106063710A true CN106063710A (zh) 2016-11-02
CN106063710B CN106063710B (zh) 2018-11-13

Family

ID=57420036

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610353267.8A Active CN106063710B (zh) 2016-05-25 2016-05-25 一种超声成像方法、系统及装置

Country Status (2)

Country Link
CN (1) CN106063710B (zh)
WO (1) WO2017202362A1 (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017202362A1 (zh) * 2016-05-25 2017-11-30 深圳开立生物医疗科技股份有限公司 一种超声成像方法、系统及装置
CN108732070A (zh) * 2018-03-21 2018-11-02 哈尔滨工程大学 一种基于脉冲管的悬浮泥沙散射强度测量装置
CN110974293A (zh) * 2019-12-11 2020-04-10 华中科技大学 一种基于c型探头的合成孔径成像方法
CN112263274A (zh) * 2020-11-18 2021-01-26 飞依诺科技(苏州)有限公司 基于多角度超声发射自适应成像方法、设备及存储介质
CN113693624A (zh) * 2021-09-03 2021-11-26 合肥工业大学 一种高质量的超声复合成像方法
CN114052786A (zh) * 2021-10-22 2022-02-18 武汉联影医疗科技有限公司 一种超声波发射方法和系统
WO2022041616A1 (en) * 2020-11-18 2022-03-03 Wuhan United Imaging Healthcare Co., Ltd. System and method for ultrasound imaging
CN114271851A (zh) * 2021-12-22 2022-04-05 武汉中旗生物医疗电子有限公司 基于凹阵探头的成像方法、装置、设备及存储介质
WO2022105843A1 (zh) * 2020-11-18 2022-05-27 武汉联影医疗科技有限公司 一种超声波成像方法、系统和存储介质

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110613476B (zh) * 2018-06-19 2022-09-23 青岛海信医疗设备股份有限公司 超声信号的处理方法和装置
WO2020084612A1 (en) * 2018-10-23 2020-04-30 Lensfree Ltd. System and method for use in x-ray imaging
CN112315500A (zh) * 2020-11-20 2021-02-05 深圳开立生物医疗科技股份有限公司 一种超声图像生成方法、装置、超声诊断设备及存储介质
CN112505668B (zh) * 2020-12-01 2024-02-13 中国人民解放军海军工程大学 激光致声发射相控阵聚焦系统
CN113332619B (zh) * 2021-05-28 2023-06-20 西安交通大学 一种用于相变纳米液滴药物载体的超声适形激活及监控成像系统
CN113842158B (zh) * 2021-08-09 2023-11-21 中南大学 基于定焦声场的光声/超声内窥探头及动态聚焦重建算法
CN117405779B (zh) * 2023-12-12 2024-02-13 广州多浦乐电子科技股份有限公司 一种提升fpga超声全聚焦帧率的方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104655728A (zh) * 2015-01-29 2015-05-27 中国科学院声学研究所 一种声学相控阵成像方法
CN104688273A (zh) * 2015-03-16 2015-06-10 哈尔滨工业大学 基于cpu+gpu异构架构的超高速超声成像装置及方法
WO2015091519A1 (en) * 2013-12-16 2015-06-25 Universität Bern Computed ultrasound tomography in echo mode (cute) for imaging speed of sound using pulse-echo sonography
US9198636B2 (en) * 2005-07-22 2015-12-01 Shenzhen Mindray Bio-Medical Electronics Co., Ltd. Continuous transmit focusing method and apparatus for ultrasound imaging system

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6773399B2 (en) * 2001-10-20 2004-08-10 Zonare Medical Systems, Inc. Block-switching in ultrasound imaging
CN101874743B (zh) * 2009-04-30 2013-10-02 深圳迈瑞生物医疗电子股份有限公司 一种多波束处理的方法和装置
JP2012105959A (ja) * 2010-10-20 2012-06-07 Toshiba Corp 超音波診断装置及び超音波送受信方法
CN103784166B (zh) * 2014-03-03 2015-08-19 哈尔滨工业大学 多功能一体化数字超声诊断系统
CN104055537B (zh) * 2014-07-16 2016-04-20 深圳开立生物医疗科技股份有限公司 用于超声设备的模式化设置的装置和相关方法
CN106063710B (zh) * 2016-05-25 2018-11-13 深圳开立生物医疗科技股份有限公司 一种超声成像方法、系统及装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9198636B2 (en) * 2005-07-22 2015-12-01 Shenzhen Mindray Bio-Medical Electronics Co., Ltd. Continuous transmit focusing method and apparatus for ultrasound imaging system
WO2015091519A1 (en) * 2013-12-16 2015-06-25 Universität Bern Computed ultrasound tomography in echo mode (cute) for imaging speed of sound using pulse-echo sonography
CN104655728A (zh) * 2015-01-29 2015-05-27 中国科学院声学研究所 一种声学相控阵成像方法
CN104688273A (zh) * 2015-03-16 2015-06-10 哈尔滨工业大学 基于cpu+gpu异构架构的超高速超声成像装置及方法

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017202362A1 (zh) * 2016-05-25 2017-11-30 深圳开立生物医疗科技股份有限公司 一种超声成像方法、系统及装置
CN108732070A (zh) * 2018-03-21 2018-11-02 哈尔滨工程大学 一种基于脉冲管的悬浮泥沙散射强度测量装置
CN110974293B (zh) * 2019-12-11 2021-05-18 华中科技大学 一种基于c型探头的合成孔径成像方法
CN110974293A (zh) * 2019-12-11 2020-04-10 华中科技大学 一种基于c型探头的合成孔径成像方法
WO2022105843A1 (zh) * 2020-11-18 2022-05-27 武汉联影医疗科技有限公司 一种超声波成像方法、系统和存储介质
WO2022041616A1 (en) * 2020-11-18 2022-03-03 Wuhan United Imaging Healthcare Co., Ltd. System and method for ultrasound imaging
CN112263274A (zh) * 2020-11-18 2021-01-26 飞依诺科技(苏州)有限公司 基于多角度超声发射自适应成像方法、设备及存储介质
CN112263274B (zh) * 2020-11-18 2023-01-31 飞依诺科技股份有限公司 基于多角度超声发射自适应成像方法、设备及存储介质
EP4228515A4 (en) * 2020-11-18 2024-04-03 Wuhan United Imaging Healthcare Co Ltd ULTRASONIC IMAGING SYSTEM AND METHOD
CN113693624A (zh) * 2021-09-03 2021-11-26 合肥工业大学 一种高质量的超声复合成像方法
CN113693624B (zh) * 2021-09-03 2023-07-18 合肥工业大学 一种高质量的超声复合成像方法
CN114052786A (zh) * 2021-10-22 2022-02-18 武汉联影医疗科技有限公司 一种超声波发射方法和系统
CN114271851A (zh) * 2021-12-22 2022-04-05 武汉中旗生物医疗电子有限公司 基于凹阵探头的成像方法、装置、设备及存储介质
CN114271851B (zh) * 2021-12-22 2023-08-29 武汉中旗生物医疗电子有限公司 基于凹阵探头的成像方法、装置、设备及存储介质

Also Published As

Publication number Publication date
WO2017202362A1 (zh) 2017-11-30
CN106063710B (zh) 2018-11-13

Similar Documents

Publication Publication Date Title
CN106063710A (zh) 一种超声成像方法、系统及装置
US9949718B2 (en) Method and system for controlling communication of data via digital demodulation in a diagnostic ultrasound system
US6527720B1 (en) Medical ultrasonic imaging method and system for spatial compounding
CN1593349B (zh) 用于实现斑点减少滤波器的系统和方法
CN100518656C (zh) 超声波摄像装置及其摄像数据处理方法以及处理程序
KR101511083B1 (ko) 합성 영상을 이용한 조직 도플러 영상 생성 장치 및 방법
JP6342212B2 (ja) 超音波診断装置
US6669641B2 (en) Method of and system for ultrasound imaging
CN110101411B (zh) 超声成像空间复合方法及系统
US9824442B2 (en) View direction adaptive volume ultrasound imaging
CN104902825B (zh) 超声波拍摄装置
JPH08336533A (ja) 3次元超音波画像作成方法および画像処理装置
CN105411626B (zh) 一种基于超声ct的合成孔径成像方法及系统
EP1927015A1 (en) Ultrasound imaging system and method for flow imaging using real-time spatial compounding
US20170238908A1 (en) Ultrasound diagnostic device
CN105074498A (zh) 具有对梯形的扇区的空间复合的超声诊断成像系统
EP3017767B1 (en) Ultrasonic diagnostic apparatus and method of controlling the same
US20160262729A1 (en) Systems and methods of reducing ultrasonic speckle using harmonic compounding
US6585648B1 (en) System, method and machine readable program for performing ultrasonic fat beam transmission and multiline receive imaging
US11408987B2 (en) Ultrasonic imaging with multi-scale processing for grating lobe suppression
CN107997783B (zh) 一种基于超声波方向性的自适应超声波束合成方法和系统
JP2022502113A (ja) 高い表示フレームレートで幅広集束送信ビームを用いる3d超音波イメージング
US6740034B2 (en) Three-dimensional ultrasound imaging system for performing receive-focusing at voxels corresponding to display pixels
KR101911470B1 (ko) 평면파를 이용하여 3차원 초음파 영상을 생성하는 장치 및 방법
CN111012379B (zh) 用于执行超声成像的方法和系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant